title: Survival analysis of DNA mutation motifs with penalized proportional hazards

abstract: Antibodies, an essential part of our immune system, develop through an intricate process to bind a wide array of pathogens. This process involves randomly mutating DNA sequences encoding these antibodies to find variants with improved binding, though mutations are not distributed uniformly across sequence sites. Immunologists observe this nonuniformity to be consistent with "mutation motifs", which are short DNA subsequences that affect how likely a given site is to experience a mutation. Quantifying the effect of motifs on mutation rates is challenging: a large number of possible motifs makes this statistical problem high dimensional, while the unobserved history of the mutation process leads to a nontrivial missing data problem. We introduce an $\ell_1$-penalized proportional hazards model to infer mutation motifs and their effects. In order to estimate model parameters, our method uses a Monte Carlo EM algorithm to marginalize over the unknown ordering of mutations. We show that our method performs better on simulated data compared to current methods and leads to more parsimonious models. The application of proportional hazards to mutation processes is, to our knowledge, novel and formalizes the current methods in a statistical framework that can be easily extended to analyze the effect of other biological features on mutation rates.


title: Not so many non-disjoint translations

abstract: We show that, consistently, there is a Borel set which has uncountably many pairwise very non-disjoint translations, but does not allow a perfect set of such translations.


title: On The Time Constant for Last Passage Percolation on Complete Graph

abstract: This paper focuses on the time constant for last passage percolation on complete graph. Let $G_n=([n],E_n)$ be the complete graph on vertex set $[n]=\{1,2,\ldots,n\}$, and i.i.d. sequence $\{X_e:e\in E_n\}$ be the passage times of edges. Denote by $W_n$ the largest passage time among all self-avoiding paths from 1 to $n$. First, it is proved that $W_n/n$ converges to constant $\mu$, where $\mu$ is called the time constant and coincides with the essential supremum of $X_e$. Second, when $\mu<\infty$, it is proved that the deviation probability $P(W_n/n\leq \mu-x)$ decays as fast as $e^{-\Theta(n^2)}$, and as a corollary, an upper bound for the variance of $W_n$ is obtained. Finally, when $\mu=\infty$, lower and upper bounds for $W_n/n$ are given.


title: A GPU Poisson-Fermi Solver for Ion Channel Simulations

abstract: The Poisson-Fermi model is an extension of the classical Poisson-Boltzmann model to include the steric and correlation effects of ions and water treated as nonuniform spheres in aqueous solutions. Poisson-Boltzmann electrostatic calculations are essential but computationally very demanding for molecular dynamics or continuum simulations of complex systems in molecular biophysics and electrochemistry. The graphic processing unit (GPU) with enormous arithmetic capability and streaming memory bandwidth is now a powerful engine for scientific as well as industrial computing. We propose two parallel GPU algorithms, one for linear solver and the other for nonlinear solver, for solving the Poisson-Fermi equation approximated by the standard finite difference method in 3D to study biological ion channels with crystallized structures from the Protein Data Bank, for example. Numerical methods for both linear and nonlinear solvers in the parallel algorithms are given in detail to illustrate the salient features of the CUDA (compute unified device architecture) software platform of GPU in implementation. It is shown that the parallel algorithms on GPU over the sequential algorithms on CPU (central processing unit) can achieve 22.8x and 16.9x speedups for the linear solver time and total runtime, respectively.


title: Going Further with Point Pair Features

abstract: Point Pair Features is a widely used method to detect 3D objects in point clouds, however they are prone to fail in presence of sensor noise and background clutter. We introduce novel sampling and voting schemes that significantly reduces the influence of clutter and sensor noise. Our experiments show that with our improvements, PPFs become competitive against state-of-the-art methods as it outperforms them on several objects from challenging benchmarks, at a low computational cost.


title: Mathematics of Isogeny Based Cryptography

abstract: These lectures notes were written for a summer school on Mathematics for post-quantum cryptography in Thi\`es, Senegal. They try to provide a guide for Masters' students to get through the vast literature on elliptic curves, without getting lost on their way to learning isogeny based cryptography. They are by no means a reference text on the theory of elliptic curves, nor on cryptography; students are encouraged to complement these notes with some of the books recommended in the bibliography. The presentation is divided in three parts, roughly corresponding to the three lectures given. In an effort to keep the reader interested, each part alternates between the fundamental theory of elliptic curves, and applications in cryptography. We often prefer to have the main ideas flow smoothly, rather than having a rigorous presentation as one would have in a more classical book. The reader will excuse us for the inaccuracies and the omissions.


title: Towards Planning and Control of Hybrid Systems with Limit Cycle using LQR Trees

abstract: We present a multi-query recovery policy for a hybrid system with goal limit cycle. The sample trajectories and the hybrid limit cycle of the dynamical system are stabilized using locally valid Time Varying LQR controller policies which probabilistically cover a bounded region of state space. The original LQR Tree algorithm builds such trees for non-linear static and non-hybrid systems like a pendulum or a cart-pole. We leverage the idea of LQR trees to plan with a continuous control set, unlike methods that rely on discretization like dynamic programming to plan for hybrid dynamical systems where it is hard to capture the exact event of discrete transition. We test the algorithm on a compass gait model by stabilizing a dynamic walking hybrid limit cycle with point foot contact from random initial conditions. We show results from the simulation where the system comes back to a stable behavior with initial position or velocity perturbation and noise.


title: On static and spherically symmetric solutions of Starobinsky model

abstract: We investigate the problem of static and spherically symmetric solutions in the Starobinsky gravity model. By extending the Lichnerowicz and Israel theorems, William Nelson have demonstrated that the Schwarzschild solution is the unique static, spherically symmetric, and asymptotically flat black hole solution in the Starobinsky model. However, Hong L{\"u} et al find that there are sign errors in the proof of Nelson. This raises the problem whether Nelson's proof is correct or not. In order to answer this question, we explore the corresponding solutions by using the Taylor series expansion method. We find that Nelson's conclusion is indeed correct despite the flaw in the proof.


title: Entanglement, non-Markovianity, and causal non-separability

abstract: Quantum mechanics, in principle, allows for processes with indefinite causal order. However, most of these causal anomalies have not yet been detected experimentally. We show that every such process can be simulated experimentally by means of non-Markovian dynamics with a measurement on additional degrees of freedom. Explicitly, we provide a constructive scheme to implement arbitrary acausal processes. Furthermore, we give necessary and sufficient conditions for open system dynamics with measurement to yield processes that respect causality locally, and find that tripartite entanglement and nonlocal unitary transformations are crucial requirements for the simulation of causally indefinite processes. These results show a direct connection between three counter-intuitive concepts: non-Markovianity, entanglement, and causal indefiniteness.


title: Communication Complexity of Discrete Fair Division

abstract: We initiate the study of the communication complexity of fair division with indivisible goods. We focus on some of the most well-studied fairness notions (envy-freeness, proportionality, and approximations thereof) and valuation classes (submodular, subadditive and unrestricted). Within these parameters, our results completely resolve whether the communication complexity of computing a fair allocation (or determining that none exist) is polynomial or exponential (in the number of goods), for every combination of fairness notion, valuation class, and number of players, for both deterministic and randomized protocols.


title: Downscaling data assimilation algorithm with applications to statistical solutions of the Navier-Stokes equations

abstract: Based on a previously introduced downscaling data assimilation algorithm, which employs a nudging term to synchronize the coarse mesh spatial scales, we construct a determining map for recovering the full trajectories from their corresponding coarse mesh spatial trajectories, and investigate its properties. This map is then used to develop a downscaling data assimilation scheme for statistical solutions of the two-dimensional Navier-Stokes equations, where the coarse mesh spatial statistics of the system is obtained from discrete spatial measurements. As a corollary, we deduce that statistical solutions for the Navier-Stokes equations are determined by their coarse mesh spatial distributions. Notably, we present our results in the context of the Navier-Stokes equations; however, the tools are general enough to be implemented for other dissipative evolution equations.


title: Reuters Tracer: Toward Automated News Production Using Large Scale Social Media Data

abstract: To deal with the sheer volume of information and gain competitive advantage, the news industry has started to explore and invest in news automation. In this paper, we present Reuters Tracer, a system that automates end-to-end news production using Twitter data. It is capable of detecting, classifying, annotating, and disseminating news in real time for Reuters journalists without manual intervention. In contrast to other similar systems, Tracer is topic and domain agnostic. It has a bottom-up approach to news detection, and does not rely on a predefined set of sources or subjects. Instead, it identifies emerging conversations from 12+ million tweets per day and selects those that are news-like. Then, it contextualizes each story by adding a summary and a topic to it, estimating its newsworthiness, veracity, novelty, and scope, and geotags it. Designing algorithms to generate news that meets the standards of Reuters journalists in accuracy and timeliness is quite challenging. But Tracer is able to achieve competitive precision, recall, timeliness, and veracity on news detection and delivery. In this paper, we reveal our key algorithm designs and evaluations that helped us achieve this goal, and lessons learned along the way.


title: Towards ECDSA key derivation from deep embeddings for novel Blockchain applications

abstract: In this work, we propose a straightforward method to derive Elliptic Curve Digital Signature Algorithm (ECDSA) key pairs from embeddings created using Deep Learning and Metric Learning approaches. We also show that these keys allows the derivation of cryptocurrencies (such as Bitcoin) addresses that can be used to transfer and receive funds, allowing novel Blockchain-based applications that can be used to transfer funds or data directly to domains such as image, text, sound or any other domain where Deep Learning can extract high-quality embeddings; providing thus a novel integration between the properties of the Blockchain-based technologies such as trust minimization and decentralization together with the high-quality learned representations from Deep Learning techniques.


title: Peer-to-Peer Control of Microgrids

abstract: In this paper, the motivation to develop microgrids as an effective solution for the control of distribution networks with high level penetration of Distributed Energy Resources(DERs) is discussed. As many different control methods for microgrids can be found in literature, this paper proposes a classification from highly centralized to distributed peer-to-peer control architectures. A peer-to-peer control paradigm is proposed as a way to control the distribution network with a high penetration of distributed energy resources. Different control algorithms suited for the proposed peer-to-peer control strategy are discussed.


title: KBGAN: Adversarial Learning for Knowledge Graph Embeddings

abstract: We introduce KBGAN, an adversarial learning framework to improve the performances of a wide range of existing knowledge graph embedding models. Because knowledge graphs typically only contain positive facts, sampling useful negative training examples is a non-trivial task. Replacing the head or tail entity of a fact with a uniformly randomly selected entity is a conventional method for generating negative facts, but the majority of the generated negative facts can be easily discriminated from positive facts, and will contribute little towards the training. Inspired by generative adversarial networks (GANs), we use one knowledge graph embedding model as a negative sample generator to assist the training of our desired model, which acts as the discriminator in GANs. This framework is independent of the concrete form of generator and discriminator, and therefore can utilize a wide variety of knowledge graph embedding models as its building blocks. In experiments, we adversarially train two translation-based models, TransE and TransD, each with assistance from one of the two probability-based models, DistMult and ComplEx. We evaluate the performances of KBGAN on the link prediction task, using three knowledge base completion datasets: FB15k-237, WN18 and WN18RR. Experimental results show that adversarial training substantially improves the performances of target embedding models under various settings.


title: On the solution of Stokes equation on regions with corners

abstract: In Stokes flow, the stream function associated with the velocity of the fluid satisfies the biharmonic equation. The detailed behavior of solutions to the biharmonic equation on regions with corners has been historically difficult to characterize. The problem was first examined by Lord Rayleigh in 1920; in 1973, the existence of infinite oscillations in the domain Green's function was proven in the case of the right angle by S.~Osher. In this paper, we observe that, when the biharmonic equation is formulated as a boundary integral equation, the solutions are representable by rapidly convergent series of the form $\sum_{j} ( c_{j} t^{\mu_{j}} \sin{(\beta_{j} \log{(t)})} + d_{j} t^{\mu_{j}} \cos{(\beta_{j} \log{(t)})} )$, where $t$ is the distance from the corner and the parameters $\mu_{j},\beta_{j}$ are real, and are determined via an explicit formula depending on the angle at the corner. In addition to being analytically perspicuous, these representations lend themselves to the construction of highly accurate and efficient numerical discretizations, significantly reducing the number of degrees of freedom required for the solution of the corresponding integral equations. The results are illustrated by several numerical examples.


title: Completion of the integrable coupling systems

abstract: In this paper, we proposed an procedure to construct the completion of the integrable system by adding a perturbation to the generalized matrix problem, which can be used to continuous integrable couplings, discrete integrable couplings and super integrable couplings. As example, we construct the completion of the Kaup-Newell (KN) integrable coupling, the Wadati-Konno-Ichikawa (WKI) integrable couplingsis, vector Ablowitz-Kaup-Newell-Segur (vAKNS) integrable couplings, the Volterra integrable couplings, Dirac type integrable couplings and NLS-mKdV type integrable couplings.


title: Fast forward of adiabatic spin dynamics of entangled states

abstract: We develop a scheme of fast forward of adiabatic spin dynamics of quantum entangled states. We settle the quasi-adiabatic dynamics by adding the regularization terms to the original Hamiltonian and then accelerate it with use of a large time-scaling factor. Assuming the experimentally-realizable candidate Hamiltonian consisting of the exchange interactions and magnetic field, we solved the regularization terms. These terms multiplied by the velocity function give rise to the state-dependent counter-diabatic terms. The scheme needs neither knowledge of full spectral properties of the system nor solving the initial and boundary value problem. Our fast forward Hamiltonian generates a variety of state-dependent counter-diabatic terms for each of adiabatic states, which can include the state-independent one. We highlight this fact by using minimum (two-spin) models for a simple transverse Ising model, quantum annealing and generation of entanglement.


title: Towards Automated ICD Coding Using Deep Learning

abstract: International Classification of Diseases(ICD) is an authoritative health care classification system of different diseases and conditions for clinical and management purposes. Considering the complicated and dedicated process to assign correct codes to each patient admission based on overall diagnosis, we propose a hierarchical deep learning model with attention mechanism which can automatically assign ICD diagnostic codes given written diagnosis. We utilize character-aware neural language models to generate hidden representations of written diagnosis descriptions and ICD codes, and design an attention mechanism to address the mismatch between the numbers of descriptions and corresponding codes. Our experimental results show the strong potential of automated ICD coding from diagnosis descriptions. Our best model achieves 0.53 and 0.90 of F1 score and area under curve of receiver operating characteristic respectively. The result outperforms those achieved using character-unaware encoding method or without attention mechanism. It indicates that our proposed deep learning model can code automatically in a reasonable way and provide a framework for computer-auxiliary ICD coding.


title: Differential Performance Debugging with Discriminant Regression Trees

abstract: Differential performance debugging is a technique to find performance problems. It applies in situations where the performance of a program is (unexpectedly) different for different classes of inputs. The task is to explain the differences in asymptotic performance among various input classes in terms of program internals. We propose a data-driven technique based on discriminant regression tree (DRT) learning problem where the goal is to discriminate among different classes of inputs. We propose a new algorithm for DRT learning that first clusters the data into functional clusters, capturing different asymptotic performance classes, and then invokes off-the-shelf decision tree learning algorithms to explain these clusters. We focus on linear functional clusters and adapt classical clustering algorithms (K-means and spectral) to produce them. For the K-means algorithm, we generalize the notion of the cluster centroid from a point to a linear function. We adapt spectral clustering by defining a novel kernel function to capture the notion of linear similarity between two data points. We evaluate our approach on benchmarks consisting of Java programs where we are interested in debugging performance. We show that our algorithm significantly outperforms other well-known regression tree learning algorithms in terms of running time and accuracy of classification.


title: Quantum Thermodynamics for Driven Dissipative Bosonic Systems

abstract: We investigate two prototypical dissipative bosonic systems under slow driving and arbitrary system-bath coupling strength, recovering their dynamic evolution as well as the heat and work rates, and we verify that thermodynamic laws are respected. Specifically, we look at the damped harmonic oscillator and the damped two-level system. For the former, we study independently the slow time- dependent perturbation in the oscillator frequency and in the coupling strength. For the latter, we concentrate on the slow modulation of the energy gap between the two levels. Importantly, we are able to find the entropy production rates for each case without explicitly defining nonequilibrium extensions for the entropy functional. This analysis also permits the definition of phenomenological friction coefficients in terms of structural properties of the system-bath composite.


title: Parkinson's Disease Digital Biomarker Discovery with Optimized Transitions and Inferred Markov Emissions

abstract: We search for digital biomarkers from Parkinson's Disease by observing approximate repetitive patterns matching hypothesized step and stride periodic cycles. These observations were modeled as a cycle of hidden states with randomness allowing deviation from a canonical pattern of transitions and emissions, under the hypothesis that the averaged features of hidden states would serve to informatively characterize classes of patients/controls. We propose a Hidden Semi-Markov Model (HSMM), a latent-state model, emitting 3D-acceleration vectors. Transitions and emissions are inferred from data. We fit separate models per unique device and training label. Hidden Markov Models (HMM) force geometric distributions of the duration spent at each state before transition to a new state. Instead, our HSMM allows us to specify the distribution of state duration. This modified version is more effective because we are interested more in each state's duration than the sequence of distinct states, allowing inclusion of these durations the feature vector.


title: Fine Grained Knowledge Transfer for Personalized Task-oriented Dialogue Systems

abstract: Training a personalized dialogue system requires a lot of data, and the data collected for a single user is usually insufficient. One common practice for this problem is to share training dialogues between different users and train multiple sequence-to-sequence dialogue models together with transfer learning. However, current sequence-to-sequence transfer learning models operate on the entire sentence, which might cause negative transfer if different personal information from different users is mixed up. We propose a personalized decoder model to transfer finer granularity phrase-level knowledge between different users while keeping personal preferences of each user intact. A novel personal control gate is introduced, enabling the personalized decoder to switch between generating personalized phrases and shared phrases. The proposed personalized decoder model can be easily combined with various deep models and can be trained with reinforcement learning. Real-world experimental results demonstrate that the phrase-level personalized decoder improves the BLEU over multiple sentence-level transfer baseline models by as much as 7.5%.


title: Compatible associative bialgebras

abstract: We introduce a non-symmetric operad $\mathcal{N}$, whose dimension in degree $n$ is given by the Catalan number $c_{n-1}$. It arises naturally in the study of coalgebra structures defined on compatible associative algebras. We prove that any free compatible associative algebra admits a compatible infinitesimal bialgebra structure, whose subspace of primitive elements is a $\mathcal{N}$-algebra. The data $({\rm As},{\rm As}^2, \mathcal{N})$ is a good triple of operads, in J.-L. Loday's sense. Our construction induces another triple of operads $({\rm As},{\rm As}_2,{\rm As})$, where ${\rm As}_2$ is the operad of matching dialgebras. Motivated by A. Goncharov's Hopf algebra of paths $P(S)$, we introduce the notion of bi-matching dialgebras and show that the Hopf algebra $P(S)$ is a bi-matching dialgebras.


title: On the second order derivative estimates for degenerate parabolic equations

abstract: We study the parabolic equation \begin{align} \notag &u_t(t,x)=a^{ij}(t)u_{x^ix^j}(t,x)+f(t,x), \quad (t,x) \in [0,T] \times \mathbf{R}^d \\ &u(0,x)=u_0(x) \label{main eqn} \end{align} with the full degeneracy of the leading coefficients, that is, \begin{align} (a^{ij}(t)) \geq \delta(t)I_{d\times d} \geq 0. \end{align} It is well known that if $f$ and $u_0$ are not smooth enough, say $f\in \mathbb{L}_p(T):=L_p([0,T] ; L_p(\mathbf{R}^d))$ and $u_0\in L_p(\mathbf{R}^d)$, then in general the solution is only in $C([0,T];L_p(\mathbf{R}^d))$, and thus derivative estimates are not possible. In this article we prove that $u_{xx}(t,\cdot)\in L_p(\mathbf{R}^d)$ on the set $\{t: \delta(t)>0 \}$ and \begin{align*} \int^T_0 \|u_{xx}(t)\|^p_{L_p} \delta(t)dt\leq N(d,p) \left(\int^T_0 \|f(t)\|^p_{L_p}\delta^{1-p}(t)dt + \|u_0\|^p_{B^{2-2/ p}_p} \right), \end{align*} where $B^{2-2/ p}_p$ is the Besov space of order $2-2/p$. We also prove that $u_{xx}(t,\cdot)\in L_p(\mathbf{R}^d)$ for all $t>0$ and \begin{equation} \label{10.13.3} \int^T_0 \|u_{xx}\|^p_{L_p(\mathbf{R}^d)}\,dt \leq N \|u_0\|^p_{B^{2-2/(\beta p)}_p}, \end{equation} if $f=0$, $\int^t_0 \delta(s)ds>0$ for each $t>0$, and a certain asymptotic behavior of $\delta(t)$ holds near $t=0$ (see (1.3)). Here $\beta>0$ is the constant related to the asymptotic behavior in (1.3). For instance, if $d=1$ and $a^{11}(t)=\delta(t)=1+\sin(1/t)$, then the estimate holds with $\beta=1$, which actually equals the maximal regularity of the heat equation $u_t=\Delta u$.


title: On epimorphisms of ordered algebras

abstract: We prove that epimorphisms are surjective in certain categories of ordered F-algebras. It then turns out that epimorphisms are also surjective in the category of all (unordered) algebras of type F.


title: From Near to Eternity: Spin-glass planting, tiling puzzles, and constraint satisfaction problems

abstract: We present a methodology for generating Ising Hamiltonians of tunable complexity and with a priori known ground states based on a decomposition of the model graph into edge-disjoint subgraphs. The idea is illustrated with a spin-glass model defined on a cubic lattice, where subproblems, whose couplers are restricted to the two values {-1,+1}, are specified on unit cubes and are parametrized by their local degeneracy. The construction is shown to be equivalent to a type of three-dimensional constraint satisfaction problem known as the tiling puzzle. By varying the proportions of subproblem types, the Hamiltonian can span a dramatic range of typical computational complexity, from fairly easy to many orders of magnitude more difficult than prototypical bimodal and Gaussian spin glasses in three space dimensions. We corroborate this behavior via experiments with different algorithms and discuss generalizations and extensions to different types of graphs.


title: Accelerating flux balance calculations in genome-scale metabolic models by localizing the application of loopless constraints

abstract: Background: Genome-scale metabolic network models and constraint-based modeling techniques have become important tools for analyzing cellular metabolism. Thermodynamically infeasible cy-cles (TICs) causing unbounded metabolic flux ranges are often encountered. TICs satisfy the mass balance and directionality constraints but violate the second law of thermodynamics. Current prac-tices involve implementing additional constraints to ensure not only optimal but also loopless flux distributions. However, the mixed integer linear programming problems required to solve become computationally intractable for genome-scale metabolic models. Results: We aimed to identify the fewest needed constraints sufficient for optimality under the loop-less requirement. We found that loopless constraints are required only for the reactions that share elementary flux modes representing TICs with reactions that are part of the objective function. We put forth the concept of localized loopless constraints (LLCs) to enforce this minimal required set of loopless constraints. By combining with a novel procedure for minimal null-space calculation, the computational time for loopless flux variability analysis is reduced by a factor of 10-150 compared to the original loopless constraints and by 4-20 times compared to the currently fastest method Fast-SNP with the percent improvement increasing with model size. Importantly, LLCs offer a scalable strategy for loopless flux calculations for multi-compartment/multi-organism models of very large sizes (e.g. >104 reactions) not feasible before. Matlab functions are available at https://github.com/maranasgroup/lll-FVA.


title: Symmetric weighted odd-power variations of fractional Brownian motion and applications

abstract: We prove a non-central limit theorem for the symmetric weighted odd-power variations of the fractional Brownian motion with Hurst parameter H< 1/2. As applications, we study the asymptotic behavior of the trapezoidal weighted odd-power variations of the fractional Brownian motion and the fractional Brownian motion in Brownian time Z_t:= X_{Y_t}, t >= 0, where X is a fractional Brownian motion and Y is an independent Brownian motion.


title: Radio and infrared study of southern H II regions G346.056-0.021 and G346.077-0.056

abstract: Aims. We present a multiwavelength study of two southern Galactic H II regions G346.056-0.021 and G346.077-0.056 which are located at a distance of 10.9 kpc. The distribution of ionized gas, cold and warm dust and the stellar population associated with the two H II regions are studied in detail using measurements at near-infrared, mid-infrared, far-infrared, submillimeter and radio wavelengths. Methods. The radio continuum maps at 1280 and 610 MHz were obtained using the Giant Metrewave Radio Telescope to probe the ionized gas. The dust temperature, column density and dust emissivity maps were generated by using modified blackbody fits in the far-infrared wavelength range 160 - 500 {\mu}m. Various near- and mid-infrared colour and magnitude criteria were adopted to identify candidate ionizing star(s) and the population of young stellar objects in the associated field. Results. The radio maps reveal the presence diffuse ionized emission displaying distinct cometary morphologies. The 1280 MHz flux densities translate to ZAMS spectral types in the range O7.5V - O7V and O8.5V - O8V for the ionizing stars of G346.056-0.021 and G346.077-0.056, respectively. A few promising candidate ionizing star(s) are identified using near-infrared photometric data. The column density map shows the presence of a large, dense dust clump enveloping G346.077-0.056. The dust temperature map shows peaks towards the two H II regions. The submillimetre image shows the presence of two additional clumps one being associated with G346.056-0.021. The masses of the clumps are estimated to range between {\sim} 1400 to 15250 M{\sun}. Based on simple analytic calculations and the correlation seen between the ionized gas distribution and the local density structure, the observed cometary morphology in the radio maps is better explained invoking the champagne-flow model.


title: Notes on Integrable Boundary Interactions of Open $SU(4)$ Alternating Spin Chains

abstract: In arXiv:1704.05807, it was shown that the planar flavored ABJM theory is integrable in the scalar sector at two-loop order using coordinate Bethe ansatz. A salient feature of this case is that the boundary reflection matrices are anti-diagonal with respect to the chosen basis. In this paper, we relax the coefficients of the boundary terms to be general constants in order to search for integrable systems among this class. We found that at each end of the spin chain, the only integrable boundary interaction, besides the one in arXiv:1704.05807, is the one with vanishing boundary interactions leading to diagonal reflection matrices. We also construct non-supersymmetric plannar flavored ABJM theory which leads to trivial boundary interaction at both ends of the open chain from two-loop anomalous dimension matrix in the scalar sector.


title: Kondo Signatures of a Quantum Magnetic Impurity in Topological Superconductors

abstract: We study the Kondo physics of a quantum magnetic impurity in two-dimensional topological superconductors (TSCs), either intrinsic or induced on the surface of a bulk topological insulator, using a numerical renormalization group technique. We show that, despite sharing the p + ip pairing symmetry, intrinsic and extrinsic TSCs host different physical processes that produce distinct Kondo signatures. Extrinsic TSCs harbor an unusual screening mechanism involving both electron and orbital degrees of freedom that produces rich and prominent Kondo phenomena, especially an intriguing pseudospin Kondo singlet state in the superconducting gap and a spatially anisotropic spin correlation. In sharp contrast, intrinsic TSCs support a robust impurity spin doublet ground state and an isotropic spin correlation. These findings advance fundamental knowledge of novel Kondo phenomena in TSCs and suggest experimental avenues for their detection and distinction.


title: Scattering theory for multistate Schr\"odinger operators

abstract: We study multistate Schr\"odinger operators related to molecular dynamics. We consider potentials which do not necessarily decay and prove absence of the singular continuous spectrum and propagation estimates which mean the scattering at speed larger than a positive constant and decay of the state with potentials higher than considered energy at infinity. We also consider the multistate Schr\"odinger operators with many-body structures. We obtain the Mourre estimate and the minimal velocity estimate for the many-body operators. The lower bound of the velocity is determined by the distance between the energy and thresholds below the energy.


title: MojiTalk: Generating Emotional Responses at Scale

abstract: Generating emotional language is a key step towards building empathetic natural language processing agents. However, a major challenge for this line of research is the lack of large-scale labeled training data, and previous studies are limited to only small sets of human annotated sentiment labels. Additionally, explicitly controlling the emotion and sentiment of generated text is also difficult. In this paper, we take a more radical approach: we exploit the idea of leveraging Twitter data that are naturally labeled with emojis. More specifically, we collect a large corpus of Twitter conversations that include emojis in the response, and assume the emojis convey the underlying emotions of the sentence. We then introduce a reinforced conditional variational encoder approach to train a deep generative model on these conversations, which allows us to use emojis to control the emotion of the generated text. Experimentally, we show in our quantitative and qualitative analyses that the proposed models can successfully generate high-quality abstractive conversation responses in accordance with designated emotions.


title: Maximizing Algebraic Connectivity of Constrained Graphs in Adversarial Environments

abstract: This paper aims to maximize algebraic connectivity of networks via topology design under the presence of constraints and an adversary. We are concerned with three problems. First, we formulate the concave maximization topology design problem of adding edges to an initial graph, which introduces a nonconvex binary decision variable, in addition to subjugation to general convex constraints on the feasible edge set. Unlike previous methods, our method is justifiably not greedy and capable of accommodating these additional constraints. We also study a scenario in which a coordinator must selectively protect edges of the network from a chance of failure due to a physical disturbance or adversarial attack. The coordinator needs to strategically respond to the adversary's action without presupposed knowledge of the adversary's feasible attack actions. We propose three heuristic algorithms for the coordinator to accomplish the objective and identify worst-case preventive solutions. Each algorithm is shown to be effective in simulation and we provide some discussion on their compared performance.


title: A Test for Isotropy on a Sphere using Spherical Harmonic Functions

abstract: Analysis of geostatistical data is often based on the assumption that the spatial random field is isotropic. This assumption, if erroneous, can adversely affect model predictions and statistical inference. Nowadays many applications consider data over the entire globe and hence it is necessary to check the assumption of isotropy on a sphere. In this paper, a test for spatial isotropy on a sphere is proposed. The data are first projected onto the set of spherical harmonic functions. Under isotropy, the spherical harmonic coefficients are uncorrelated whereas they are correlated if the underlying fields are not isotropic. This motivates a test based on the sample correlation matrix of the spherical harmonic coefficients. In particular, we use the largest eigenvalue of the sample correlation matrix as the test statistic. Extensive simulations are conducted to assess the Type I errors of the test under different scenarios. We show how temporal correlation affects the test and provide a method for handling temporal correlation. We also gauge the power of the test as we move away from isotropy. The method is applied to the near-surface air temperature data which is part of the HadCM3 model output. Although we do not expect global temperature fields to be isotropic, we propose several anisotropic models with increasing complexity, each of which has an isotropic process as model component and we apply the test to the isotropic component in a sequence of such models as a method of determining how well the models capture the anisotropy in the fields.


title: On the maximal saddle order of p:-q resonant saddle

abstract: In this paper, we obtain some estimations of the saddle order which is the sole topological invariant of the non-integrable resonant saddles of planar polynomial vector fields of arbitrary degree $n$. Firstly, we prove that, for any given resonance $p:-q$, $(p, q)=1$, and sufficiently big integer $n$, the maximal saddle order can grow at least as rapidly as $n^2$. Secondly, we show that there exists an integer $k_0$, which grows at least as rapidly as $3n^2/2$, such that $L_{k_0}$ does not belong to the ideal generated by the first $k_0-1$ saddle values $L_1, L_2, \cdots, L_{k_0-1}$, where $L_{k}$ means the $k$-th saddle value of the given system. In particular, if $p=1$ (or $q=1$), we obtain a sharper result that $k_0$ can grow at least as rapidly as $2 n^2$.


title: Enhancing Network Embedding with Auxiliary Information: An Explicit Matrix Factorization Perspective

abstract: Recent advances in the field of network embedding have shown the low-dimensional network representation is playing a critical role in network analysis. However, most of the existing principles of network embedding do not incorporate auxiliary information such as content and labels of nodes flexibly. In this paper, we take a matrix factorization perspective of network embedding, and incorporate structure, content and label information of the network simultaneously. For structure, we validate that the matrix we construct preserves high-order proximities of the network. Label information can be further integrated into the matrix via the process of random walk sampling to enhance the quality of embedding in an unsupervised manner, i.e., without leveraging downstream classifiers. In addition, we generalize the Skip-Gram Negative Sampling model to integrate the content of the network in a matrix factorization framework. As a consequence, network embedding can be learned in a unified framework integrating network structure and node content as well as label information simultaneously. We demonstrate the efficacy of the proposed model with the tasks of semi-supervised node classification and link prediction on a variety of real-world benchmark network datasets.


title: The energy change of the complete multipartite graph

abstract: The energy of a graph is defined as the sum of the absolute values of all eigenvalues of the graph. Akbari et al. \cite{S. Akbari} proved that for a complete multipartite graph $K_{t_1 ,\ldots,t_k}$, if $t_i\geq 2 \ (i=1,\ldots,k)$, then deleting any edge will increase the energy. A natural question is how the energy changes when $\min\{t_1 ,\ldots,t_k\}=1$. In this paper, we will answer this question and completely determine how the energy of a complete multipartite graph changes when one edge is removed.


title: Artificial Interference Aided Physical Layer Security in Cache-enabled Heterogeneous Networks

abstract: Caching popular contents is a promising way to offload the mobile data traffic in wireless networks, but so far the potential advantage of caching in improving physical layer security (PLS) is rarely considered. In this paper, we contribute to the design and theoretical understanding of exploiting the caching ability of users to improve the PLS in a wireless heterogeneous network (HetNet). In such network, the base station (BS) ensures the secrecy of communication by utilizing some of the available power to transmit a pre-cached file, such that only the eavesdropper's channel is degraded. Accordingly, the node locations of BSs, users and eavesdroppers are first modeled as mutually independent poisson point processes (PPPs) and the corresponding file access protocol is developed. We then derive analytical expressions of two metrics, average secrecy rate and secrecy coverage probability, for the proposed system. Numerical results are provided to show the significant security advantages of the proposed network and to characterize the impact of network resource on the secrecy metrics.


title: The CP-Violating 2HDM in Light of a Strong First Order Electroweak Phase Transition and Implications for Higgs Pair Production

abstract: We investigate the strength of the electroweak phase transition (EWPT) within the CP-violating 2-Higgs-Doublet Model (C2HDM). By applying a renormalisation scheme which allows efficient scans of the C2HDM parameter space, we analyse the possibility of a strong first order EWPT required for baryogenesis and study its phenomenological implications for the LHC. Like in the CP-conserving (real) 2HDM (R2HDM) we find that a strong EWPT favours mass gaps between the non-SM-like Higgs bosons. These lead to prominent final states comprised of gauge+Higgs bosons or pairs of Higgs bosons. In contrast to the R2HDM, the CP-mixing of the C2HDM also favours approximately mass degenerate spectra with dominant decays into SM particles. The requirement of a strong EWPT further allows us to distinguish the C2HDM from the R2HDM using the signal strengths of the SM-like Higgs boson. We additionally find that a strong EWPT requires an enhancement of the SM-like trilinear Higgs coupling at next-to-leading order (NLO) by up to a factor of 2.4 compared to the NLO SM coupling, establishing another link between cosmology and collider phenomenology. We provide several C2HDM benchmark scenarios compatible with a strong EWPT and all experimental and theoretical constraints. We include the dominant branching ratios of the non-SM-like Higgs bosons as well as the Higgs pair production cross section of the SM-like Higgs boson for every benchmark point. The pair production cross sections can be substantially enhanced compared to the SM and could be observable at the high-luminosity LHC, allowing access to the trilinear Higgs couplings.


title: Are the chiral based $\bar{K}N$ potentials really energy dependent?

abstract: It is shown, that the energy dependence of the chiral based $\bar{K}N$ potentials, responsible for the occurrence of two poles in the $I=0$ sector, is the consequence of applying the on-shell factorization introduced in [E. Oset, A. Ramos, Nucl. Phys. A 635(1998)99]. When the dynamical equation is solved without this approximation, the $T$-matrix has only one pole in the energy region of the $\Lambda(1405)$ resonance.


title: On aggregation of multitype Galton-Watson branching processes with immigration

abstract: Limit behaviour of temporal and contemporaneous aggregations of independent copies of a stationary multitype Galton-Watson branching process with immigration is studied in the so-called iterated and simultaneous cases, respectively. In both cases, the limit process is a zero mean Brownian motion with the same covariance function under third order moment conditions on the branching and immigration distributions. We specialize our results for generalized integer-valued autoregressive processes and single-type Galton-Watson processes with immigration as well.


title: Single-particle properties of the Hubbard model in a novel three-pole approximation

abstract: We study the 2D Hubbard model using the Composite Operator Method within a novel three-pole approximation. Motivated by the long-standing experimental puzzle of the single-particle properties of the underdoped cuprates, we include in the operatorial basis, together with the usual Hubbard operators, a field describing the electronic transitions dressed by the nearest-neighbor spin fluctuations, which play a crucial role in the unconventional behavior of the Fermi surface and of the electronic dispersion. Then, we adopt this approximation to study the single-particle properties in the strong coupling regime and find an unexpected behavior of the van Hove singularity that can be seen as a precursor of a pseudogap regime.


title: Recommender Systems with Random Walks: A Survey

abstract: Recommender engines have become an integral component in today's e-commerce systems. From recommending books in Amazon to finding friends in social networks such as Facebook, they have become omnipresent. Generally, recommender systems can be classified into two main categories: content based and collaborative filtering based models. Both these models build relationships between users and items to provide recommendations. Content based systems achieve this task by utilizing features extracted from the context available, whereas collaborative systems use shared interests between user-item subsets. There is another relatively unexplored approach for providing recommendations that utilizes a stochastic process named random walks. This study is a survey exploring use cases of random walks in recommender systems and an attempt at classifying them.


title: Constraints on scalar-tensor theory of gravity by the recent observational results on gravitational waves

abstract: The speed of gravitational waves provides us a new tool to test alternative theories of gravity. The constraint on the speed of gravitational waves from GW170817 and GRB170817A is used to test some classes of Horndeski theory. In particular, we consider the coupling of a scalar field to Einstein tensor and the coupling of the Gauss-Bonnet term to a scalar field. The coupling strength of the Gauss-Bonnet coupling is constrained to be in the order of $10^{-15}$. In the Horndeski theory we show that in order for this theory to satisfy the stringent constraint on the speed of GWs the mass scale $M$ introduced in the non-minimally derivative coupling is constrained to be in the range $10^{15}\text{GeV}\gg M \gtrsim 2\times 10^{-35}$GeV taking also under consideration the early times upper bound for the mass scale $M$. The large mass ranges require no fine-tuning because the effect of non-minimally derivative coupling is negligible at late times.


title: Optimal Two-Stage Programming for Integration of PHEV Parking Lots in Industrial Microgrids

abstract: With the progressive exhaustion of fossil energy and the enhanced awareness of environmental protection, more attention is being paid to plug in hybrid electric vehicles. Inappropriate siting and sizing of plug in hybrid electric vehicles parking lots could have negative effects on the development of plug in hybrid electric vehicles, the layout of the city traffic network, and the convenience of plug in hybrid electric vehicles drivers as well as lead to an increasing in the network losses and a degradation in voltage profiles at some nodes. Given this background, this paper aims to allocate parking lots in Industrial Micro-grids with the objective of minimizing system costs including investment cost, power loss and scheduling cost as possible objectives. A two-stage model has been designed for this purpose. The optimal siting and sizing of parking lots in order to minimize the investment cost of parking lots is performed in the first stage. At the second stage, the optimal plug in hybrid electric vehicles scheduling problem is solved considering market interactions to provide profit to the parking lots owner with taken into account various network constraints. Conclusions are duly drawn with a realistic example.


title: On the parameters of intertwining codes

abstract: Let $F$ be a field and let $F^{r\times s}$ denote the space of $r\times s$ matrices over $F$. Given equinumerous subsets $\mathcal{A}=\{A_i\mid i \in I\}\subseteq F^{r\times r}$ and $\mathcal{B}=\{B_i\mid i\in I\}\subseteq F^{s\times s}$ we call the subspace $C(\mathcal{A},\mathcal{B}):=\{X\in F^{r\times s}\mid A_iX=XB_i\ {\rm for }\ i\in I\}$ an \emph{intertwining code}. We show that if $C(\mathcal{A},\mathcal{B})\ne\{0\}$, then for each $i\in I$, the characteristic polynomials of $A_i$ and $B_i$ and share a nontrivial factor. We give an exact formula for $k=\dim(C(\mathcal{A},\mathcal{B}))$ and give upper and lower bounds. This generalizes previous work in this area. Finally we construct intertwining codes with large minimum distance when the field is not `too small'. We give examples of codes where $d=rs/k=1/R$ is large where the minimum distance, dimension, and rate of the linear code $C(\mathcal{A},\mathcal{B})$ by $d$, $k$, and $R=k/rs$, respectively.


title: Unfolding multi-particle quantum correlations hidden in decoherence

abstract: Quantum coherence is a fundamental characteristic to distinguish quantum systems from their classical counterparts. Though quantum coherence persists in isolated non-interacting systems, interactions inevitably lead to decoherence, which is in general believed to cause the lost of quantum correlations. Here, we show that, accompanying to the single-particle decoherence, interactions build up quantum correlations on the two-, three-, and multi-particle levels. Using the quantitative solutions of the quantum dynamics of a condensate occupying two modes, such as two bands of an optical lattice, we find out that such dynamically emergent multi-particle correlations not only reveal how interactions control the quantum coherence of a many-body system in a highly intriguing means, but also evince the rise of exotic fragmented condensates, which are difficult to access at the ground state. We further develop a generic interferometry that can be used in experiments to measure high order correlation functions directly.


title: Relativistic and Electromagnetic Molecular Dynamics Simulations for a Carbon-Gold Nanotube Accelerator

abstract: Relativistic molecular dynamics are described in ultra-high temperature and MeV energy behaviors. In strongly-coupled systems, the Coulomb electrostatic field is collected in the infinite space, and the electromagnetic fields are added in the coordinate space. Separation of the electromagnetic and electrostatic electric fields is a good approximation for short time periods. For a numerical application, a nanotube accelerator under an $E \times B$ pulse is studied. Positive ions are accelerated in the parallel direction, whereas the electrons proceed in the perpendicular direction. Rapid expansion to infinite space and short-range electromagnetic radiation cooperate for large intensities. At $10^{22}\rm{W/cm}^{2}$, pulsation oscillations for gold and carbon ions flare up and electrons acquire the relativistic velocities. They are observed in relativistic molecular dynamics simulation.


title: Large polaron evolution in anatase TiO2 due to carrier and temperature dependence of electron-phonon coupling

abstract: The electronic and magneto transport properties of reduced anatase TiO2 epitaxial thin films are analyzed considering various polaronic effects. Unexpectedly, with increasing carrier concentration, the mobility increases, which rarely happens in common metallic systems. We find that the screening of the electron-phonon (e-ph) coupling by excess carriers is necessary to explain this unusual dependence. We also find that the magnetoresistance (MR) could be decomposed into a linear and a quadratic component, separately characterizing the transport and trap behavior of carriers as a function of temperature. The various transport behaviors could be organized into a single phase diagram which clarifies the nature of large polaron in this material.


title: Holomorphic aspects of moduli of representations of quivers

abstract: This article describes some complex-analytic aspects of the moduli space of the finite-dimensional complex representations of a finite quiver, which are stable with respect to a fixed rational weight. We construct a natural structure of a complex manifold on this moduli space, and a K\"ahler metric on the complex manifold. We then define a Hermitian holomorphic line bundle on the moduli space, and show that its curvature is a rational multiple of the K\"ahler form.


title: Natural exact covering systems and the reversion of the M\"obius series

abstract: We prove that the number of natural exact covering systems of cardinality $k$ is equal to the coefficient of $x^k$ in the reversion of the power series $\sum_{k \ge 1} \mu (k) x^k$, where $\mu(k)$ is the usual number-theoretic M\"obius function. Using this result, we deduce an asymptotic expression for the number of such systems.


title: From randomness in two symbols to randomness in three symbols

abstract: In 1909 Borel defined normality as a notion of randomness of the digits of the representation of a real number over certain base (fractional expansion). If we think the representation of a number over a base as an infinite sequence of symbols from a finite alphabet $A$, we can define normality directly for words of symbols of $A$: A word $x$ is normal to the alphabet $A$ if every finite block of symbols from $A$ appears with the same asymptotic frequency in $x$ as every other block of the same length. Many examples of normal words have been found since its definition, being Champernowne in 1933 the first to show an explicit and simple instance. Moreover, it has been characterized how we can select subsequences of a normal word $x$ preserving its normality, always leaving the alphabet $A$ fixed. In this work we consider the dual problem which consists of inserting symbols in infinite positions of a given word, in such a way that normality is preserved. Specifically, given a symbol $s$ that is not present on the original alphabet $A$ and given a word $x$ that is normal to the alphabet $A$ we solve how to insert the symbol $s$ in infinite positions of the word $x$ such that the resulting word is normal to the expanded alphabet $A\cup \{s\}$.


title: The Origin of the Beauregard--Suryanarayan Product on Pythagorean Triples

abstract: We show how the multiplicative structure on Pythagorean triples defined by Beauregard--Suryanarayan arises from automorphisms of the hyperbolic plane, and how the structure theorem of this product follow naturally from this connection. We also show how the natural involution on Pythagorean triples is related to Pell's equation with the parameter 2, and prove a generalization of this phenomenon.


title: A generalization of Bohr's Equivalence Theorem

abstract: Based on a generalization of Bohr's equivalence relation for general Dirichlet series, in this paper we study the sets of values taken by certain classes of equivalent almost periodic functions in their strips of almost periodicity. In fact, the main result of this paper consists of a result like Bohr's equivalence theorem extended to the case of these functions.


title: Biconservative ideal hypersurfaces in Euclidean spaces

abstract: A biconservative submanifold of a Riemannian manifold is a sub- manifold with divergence free stress-energy tensor with respect to bienergy. These are generalizations of biharamonic submanifolds. In 2013, B. Y. Chen and M.I. Munteanu proved that $\delta(2)$-ideal and $\delta(3)$-ideal biharmonic hypersurfaces in Euclidean space are minimal. In this paper, we generalize this result for $\delta(2)$-ideal and $\delta(3)$-ideal bisonservative hypersurfaces in Euclidean space. Also, we study $\delta(4)$-ideal biconservative hypersurfaces in Euclidean space $\mathbb{E}^{6}$ having constant scalar curvature. We prove that such a hypersurface must be of constant mean curvature.


title: Mobile Sensing of Two-Dimensional Bandlimited Fields on Random Paths

abstract: Mobile sensing has been recently proposed for sampling spatial fields, where mobile sensors record the field along various paths for reconstruction. Classical and contemporary sampling typically assumes that the sampling locations are approximately known. This work explores multiple sampling strategies along random paths to sample and reconstruct a two dimensional bandlimited field. Extensive simulations are carried out, with insights from sensing matrices and their properties, to evaluate the sampling strategies. Their performance is measured by evaluating the stability of field reconstruction from field samples. The effect of location unawareness on some sampling strategies is also evaluated by simulations.


title: Discovering conversational topics and emotions associated with Demonetization tweets in India

abstract: Social media platforms contain great wealth of information which provides us opportunities explore hidden patterns or unknown correlations, and understand people's satisfaction with what they are discussing. As one showcase, in this paper, we summarize the data set of Twitter messages related to recent demonetization of all Rs. 500 and Rs. 1000 notes in India and explore insights from Twitter's data. Our proposed system automatically extracts the popular latent topics in conversations regarding demonetization discussed in Twitter via the Latent Dirichlet Allocation (LDA) based topic model and also identifies the correlated topics across different categories. Additionally, it also discovers people's opinions expressed through their tweets related to the event under consideration via the emotion analyzer. The system also employs an intuitive and informative visualization to show the uncovered insight. Furthermore, we use an evaluation measure, Normalized Mutual Information (NMI), to select the best LDA models. The obtained LDA results show that the tool can be effectively used to extract discussion topics and summarize them for further manual analysis.


title: A Mean Field Games approach for multi-lane traffic management

abstract: In this work we discuss an Mean Field Games approach to traffic management on multi-lane roads. Such approach is particularly indicated to model self driven vehicles with perfect information of the domain. The mathematical interest of the problem is related to the fact that the system of partial differential equations obtained in this case is not in the classic form, but it consists of two continuity equations (one for each lane) and a variational inequality, coming from the Hamilton-Jacobi theory of the hybrid control.


title: Will Solar Cycles 25 and 26 Be Weaker than Cycle 24 ?

abstract: The study of variations in solar activity is important for understanding the underlying mechanism of solar activity and for predicting the level of activity in view of the activity impact on space weather and global climate. Here we have used the amplitudes (the peak values of the 13-month smoothed international sunspot number) of Solar Cycles 1-24 to predict the relative amplitudes of the solar cycles during the rising phase of the upcoming Gleissberg cycle. We fitted a cosine function to the amplitudes and times of the solar cycles after subtracting a linear fit of the amplitudes. The best cosine fit shows overall properties (periods, maxima, minima, etc.) of Gleissberg cycles, but with large uncertainties. We obtain a pattern of the rising phase of the upcoming Gleissberg cycle, but there is considerable ambiguity. Using the epochs of violations of the Gnevyshev-Ohl rule (G-O rule) and the `tentative inverse G-O rule' of solar cycles during the period 1610-2015, and also using the epochs where the orbital angular momentum of the Sun is steeply decreased during the period 1600-2099, we infer that Solar Cycle 25 will be weaker than Cycle 24. Cycles 25 and 26 will have almost same strength, and their epochs are at the minimum between the current and upcoming Gleissberg cycles. In addition, Cycle 27 is expected to be stronger than Cycle 26 and weaker than Cycle 28, and Cycle 29 is expected to be stronger than both Cycles 28 and 30. The maximum of Cycle 29 is expected to represent the next Gleissberg maximum. Our analysis also suggests a much lower value (30-40) for the maximum amplitude of the upcoming Cycle 25.


title: Software for full-color 3D reconstruction of the biological tissues internal structure

abstract: A software for processing sets of full-color images of biological tissue histological sections is developed. We used histological sections obtained by the method of high-precision layer-by-layer grinding of frozen biological tissues. The software allows restoring the image of the tissue for an arbitrary cross-section of the tissue sample. Thus, our method is designed to create a full-color 3D reconstruction of the biological tissue structure. The resolution of 3D reconstruction is determined by the quality of the initial histological sections. The newly developed technology available to us provides a resolution of up to 5 - 10 {\mu}m in three dimensions.


title: Revealing Tripartite Quantum Discord with Tripartite Information Diagram

abstract: A new measure based on the tripartite information diagram is proposed for identifying quantum discord in tripartite systems. The proposed measure generalizes the mutual information underlying discord from bipartite to tripartite systems, and utilizes both one-particle and two-particle projective measurements to reveal the characteristics of the tripartite quantum discord. The feasibility of the proposed measure is demonstrated by evaluating the tripartite quantum discord for systems with states close to Greenberger-Horne-Zeilinger, W, and biseparable states. In addition, the connections between tripartite quantum discord and two other quantum correlations---namely genuine tripartite entanglement and genuine tripartite Einstein-Podolsky-Rosen steering---are briefly discussed. The present study considers the case of quantum discord in tripartite systems. However, the proposed framework can be readily extended to general N-partite systems.


title: Invariant surface area functionals and singular Yamabe problem in 3-dimensional CR geometry

abstract: We express two CR invariant surface area elements in terms of quantities in pseudohermitian geometry. We deduce the Euler-Lagrange equations of the associated energy functionals. Many solutions are given and discussed. In relation to the singular CR Yamabe problem, we show that one of the energy functionals appears as the coefficient (up to a constant multiple) of the log term in the associated volume renormalization.


title: Weakly Supervised Audio Source Separation via Spectrum Energy Preserved Wasserstein Learning

abstract: Separating audio mixtures into individual instrument tracks has been a long standing challenging task. We introduce a novel weakly supervised audio source separation approach based on deep adversarial learning. Specifically, our loss function adopts the Wasserstein distance which directly measures the distribution distance between the separated sources and the real sources for each individual source. Moreover, a global regularization term is added to fulfill the spectrum energy preservation property regardless separation. Unlike state-of-the-art weakly supervised models which often involve deliberately devised constraints or careful model selection, our approach need little prior model specification on the data, and can be straightforwardly learned in an end-to-end fashion. We show that the proposed method performs competitively on public benchmark against state-of-the-art weakly supervised methods.


title: Bohr's equivalence relation in the space of Besicovitch almost periodic functions

abstract: Based on Bohr's equivalence relation which was established for general Dirichlet series, in this paper we introduce a new equivalence relation on the space of almost periodic functions in the sense of Besicovitch, $B(\mathbb{R},\mathbb{C})$, defined in terms of polynomial approximations. From this, we show that in an important subspace $B^2(\mathbb{R},\mathbb{C})\subset B(\mathbb{R},\mathbb{C})$, where Parseval's equality and Riesz-Fischer theorem holds, its equivalence classes are sequentially compact and the family of translates of a function belonging to this subspace is dense in its own class.


title: OPE for XXX

abstract: We explain a new method for finding the correlation functions for the XXX model which is based on the concepts of Operator Product Expansion of Quantum Field Theory on one hand and of fermionic bases for the XXX spin chain on the other. With this method we are able to perform computations for up to 11 lattice sites. We show that these "experimental" data allow to guess exact formulae for the OPE coefficients.


title: The hydrodynamic singular vortex on the sphere and the Dirac monopole

abstract: An exact correspondence is established between mathematical description of the single "elementary vortex" (EV)velocity on the sphere (Zermelo,1902;Bogomolov,1977) and the Dirac magnetic monopole (DMM) vector potential (Dirac, 1931). Similar analogy with DMM was noted only for the vortices in quantum fluid He-3A (Blaha,1976;Volovik,Mineev,1976). Singular EV on a sphere is usually considered using compensating vortex field, uniformly distributed over the sphere. It is necessary to meet the Gauss-Kelvin theorem with zero integral vorticity over the sphere. However since Bogomolov (1977)(see also Dritshel,Boatto,2015)uniform vorticity impact on dynamics is neglected when only the velocity field of singular EV is taken into account. It is not true because the hydrodynamic equations do not allow existence of a solution in the form of an isolated EV but allow solutions in the form of the antipodal EV pairs. It is also suggested the possibility of DMM existence only in the form of the point magnetic dipoles, consisting of two DMM with different signs of the magnetic charge.


title: Stability and Stabilization of Fractional-order Systems with Different Derivative Orders: An LMI Approach

abstract: Stability and stabilization analysis of fractional-order linear time-invariant (FO-LTI) systems with different derivative orders is studied in this paper. First, by using an appropriate linear matrix function, a single-order equivalent system for the given different-order system is introduced by which a new stability condition is obtained that is easier to check in practice than the conditions known up to now. Then the stabilization problem of fractional-order linear systems with different fractional orders via a dynamic output feedback controller with a predetermined order is investigated, utilizing the proposed stability criterion. The linear matrix inequality based procedure of developing stabilizing output feedback control is preserved in spite of the complexity of assuming the most complete linear controller model, with direct feedthrough parameter. The proposed stability and stabilization theorems are applicable to FO-LTI systems with different fractional orders in one or both of and intervals. Eventually, some numerical examples are presented to confirm the obtained analytical results.


title: Adversarial Training for Disease Prediction from Electronic Health Records with Missing Data

abstract: Electronic health records (EHRs) have contributed to the computerization of patient records and can thus be used not only for efficient and systematic medical services, but also for research on biomedical data science. However, there are many missing values in EHRs when provided in matrix form, which is an important issue in many biomedical EHR applications. In this paper, we propose a two-stage framework that includes missing data imputation and disease prediction to address the missing data problem in EHRs. We compared the disease prediction performance of generative adversarial networks (GANs) and conventional learning algorithms in combination with missing data prediction methods. As a result, we obtained a level of accuracy of 0.9777, sensitivity of 0.9521, specificity of 0.9925, area under the receiver operating characteristic curve (AUC-ROC) of 0.9889, and F-score of 0.9688 with a stacked autoencoder as the missing data prediction method and an auxiliary classifier GAN (AC-GAN) as the disease prediction method. The comparison results show that a combination of a stacked autoencoder and an AC-GAN significantly outperforms other existing approaches. Our results suggest that the proposed framework is more robust for disease prediction from EHRs with missing data.


title: Behaviour of Charged Spinning Massless Particles

abstract: We revisit the classical theory of a relativistic massless charged point particle with spin and interacting with an external electromagnetic field. In particular, we give a proper definition of its kinetic energy and its total energy, the latter being conserved when the external field is stationary. We also write the conservation laws for the linear and angular momenta. Finally, we find that the particle's velocity may differ from $c$ as a result of the spin---electromagnetic field interaction, without jeopardizing Lorentz invariance.


title: Photonic Spin Hall Effect in Waveguides Composed of Two Types of Single-Negative Metamaterials

abstract: The polarization controlled optical signal routing has many important applications in photonics such as polarization beam splitter. By using two-dimensional transmission lines with lumped elements, we experimentally demonstrate the selective excitation of guided modes in waveguides composed of two kinds of single-negative metamaterials. A localized, circularly polarized emitter placed near the interface of the two kinds of single-negative metamaterials only couples with one guided mode with a specific propagating direction determined by the polarization handedness of the source. Moreover, this optical spin-orbit locking phenomenon, also called the photonic spin Hall effect, is robust against interface fluctuations, which may be very useful in the manipulation of electromagnetic signals.


title: Abelian ideals and amazing roots

abstract: Let $\mathfrak g$ be a simple Lie algebra with a Borel subalgebra $\mathfrak b$. To any long positive root $\gamma$, one associates two ideals of $\mathfrak b$: the abelian ideal $I(\gamma)_{max}$ and not necessarily abelian ideal $I\langle{\succcurlyeq}\gamma\rangle$. It is known that $I(\gamma)_{max} \subset I\langle{\succcurlyeq}\gamma\rangle$, and $\gamma$ is said to be amazing if the equality holds. The set of amazing roots, $\mathcal A$, is closed under the operation `$\vee$' in $\Delta^+$, and $\gamma\in\mathcal A$ is said to be primitive, if it cannot be written as $\gamma_1\vee\gamma_2$ with incomparable amazing roots $\gamma_1,\gamma_2$. We classify the amazing roots and notice that the number of primitive roots equals $\mathsf{rk}(\mathfrak g)$. Moreover, if $\Pi$ (resp. $\mathcal A_{\sf pr}$) is the set of simple (resp. primitive) roots, then there is a natural bijection $\Pi\longleftrightarrow \mathcal A_{\sf pr}$. We also describe the set $\mathcal A\cap{\mathcal H}$, where ${\mathcal H}$ is the Heisenberg subset of $\Delta^+$.


title: Phase imaging by spatial wavefront sampling

abstract: Phase imaging techniques extract the optical path-length information of a scene, whereas wavefront sensors provide the shape of an optical wavefront. Since these two applications have different technical requirements, they have developed their own specific technology. Here we show how to perform phase imaging combining wavefront sampling using a reconfigurable spatial light modulator with a beam position detector. The result is a time-multiplexed detection scheme, capable of being shortened considerably by compressive sensing. This robust reference-less method does not require the phase unwrapping algorithms demanded by conventional interferometry, and its lenslet-free nature removes tradeoffs usually found in Shack-Hartmann sensors.


title: On the Annihilation of Thin Sets

abstract: One says that a pair of sets $(S,Q)$ in $\mathbb{R}$ is 'annihilating' if no function can be concentrated on $S$ while having its Fourier transform concentrated on $Q$. One uses to distinguish between weak and strong annihilation types. It is well known that if both sets $S$ and $Q$ are of finite measure then they are strongly annihilating. In this paper we prove that if $S$ is a set of finite measure with periodic gaps, and $Q$ is a set of density zero, then weak annihilation holds. On the other hand a counter-example is constructed, showing that strong annihilation, in general, does not.


title: ThunderKAT: The MeerKAT Large Survey Project for Image-Plane Radio Transients

abstract: ThunderKAT is the image-plane transients programme for MeerKAT. The goal as outlined in 2010, and still today, is to find, identify and understand high-energy astrophysical processes via their radio emission (often in concert with observations at other wavelengths). Through a comprehensive and complementary programme of surveying and monitoring Galactic synchrotron transients (across a range of compact accretors and a range of other explosive phenomena) and exploring distinct populations of extragalactic synchrotron transients (microquasars, supernovae and possibly yet unknown transient phenomena) - both from direct surveys and commensal observations - we will revolutionise our understanding of the dynamic and explosive transient radio sky. As well as performing targeted programmes of our own, we have made agreements with the other MeerKAT large survey projects (LSPs) that we will also search their data for transients. This commensal use of the other surveys, which remains one of our key programme goals in 2016, means that the combined MeerKAT LSPs will produce by far the largest GHz-frequency radio transient programme to date.


title: Planar Hall effect in the Weyl semimetal GdPtBi

abstract: Observation of Weyl and Dirac Fermions in condensed matter systems is one of the most important discoveries. Among the very few available tools to characterize Weyl semimetals through electrical transport, negative magnetoresistance is most commonly used. Considering shortcomings of this method, new tools to characterize chiral anomaly in Weyl semimetals are desirable. We employ planar Hall effect as an effective technique in half Heusler Weyl semimetal GdPtBi to study chiral anomaly. This compound exhibits a large value of 1.5 mohm cm planar Hall resistivity at 2 K and in 9 T. Our analysis reveals that the observed amplitude is dominated by Berry curvature and chiral anomaly contributions. Through the angle dependent transport studies we establish that GdPtBi with relatively small orbital magnetoresistance is an ideal candidate to observe large planar Hall effect .


title: On Pure Degree Sequence Manipulations Forcing Long Cycles in Graphs

abstract: The well-known Hamiltonian sufficient conditions, proposed by Dirac, Faudree et al., P\'osa, Bondy, Chv\'atal are based on pure degree manipulations without any additional conditions. In this paper, we present two new types of pure degree manipulations that produce relatively simpler (in view of verification) results. The reverse versions (long-cycle versions) of the obtained results are presented as well. The importance of pure degree manipulations is motivated by their amenability (as starting points) to generalizations and new ideas in a great variety of ways.


title: State space models for non-stationary intermittently coupled systems: an application to the North Atlantic Oscillation

abstract: We develop Bayesian state space methods for modelling changes to the mean level or temporal correlation structure of an observed time series due to intermittent coupling with an unobserved process. Novel intervention methods are proposed to model the effect of repeated coupling as a single dynamic process. Latent time-varying autoregressive components are developed to model changes in the temporal correlation structure. Efficient filtering and smoothing methods are derived for the resulting class of models. We propose methods for quantifying the component of variance attributable to an unobserved process, the effect during individual coupling events, and the potential for skilful forecasts. The proposed methodology is applied to the study of winter-time variability in the dominant pattern of climate variation in the northern hemisphere, the North Atlantic Oscillation. Around 70% of the inter-annual variance in the winter (Dec-Jan-Feb) mean level is attributable to an unobserved process. Skilful forecasts for winter (Dec-Jan-Feb) mean are possible from the beginning of December.


title: A simple inductive proof of Levy-Steinitz theorem

abstract: We present a relatively simple inductive proof of the classical Levy-Steinitz Theorem saying that for a sequence $(x_n)_{n=1}^\infty$ in a finite-dimensional Banach space $X$ the set of all sums of rearranged series $\sum_{n=1}^\infty x_{\sigma(n)}$ is an affine subspace of $X$. This affine subspace is not empty if and only if for any linear functional $f:X\to \mathbb R$ the series $\sum_{n=1}^\infty f(x_{\sigma(n)})$ is convergent for some permutation $\sigma$ of $\mathbb N$. This gives an answer to a problem of Vaja Tarieladze, posed in Lviv Scottish Book in September, 2017. Also we construct a sequence $(x_n)_{n=1}^\infty$ in the torus $\mathbb T\times\mathbb T$ such that the series $\sum_{n=1}^\infty x_{\sigma(n)}$ is divergent for all permutations $\sigma$ of $\mathbb N$ but for any continuous homomorphism $f:\mathbb T^2\to\mathbb T$ to the circle group $\mathbb T:=\mathbb R/\mathbb Z$ the series $\sum_{n=1}^\infty f(x_{\sigma_f(n)})$ is convergent for some permutation $\sigma_f$ of $\mathbb N$. This example shows that the second part of Levy-Steinitz Theorem (characterizing sequences with non-empty set of potential sums) does not extend to locally compact Abelian groups.


title: Constraints on Born-Infeld gravity from the speed of gravitational waves after GW170817 and GRB 170817A

abstract: The observations of gravitational waves from the binary neutron star merger event GW170817 and the subsequent observation of its electromagnetic counterparts from the gamma-ray burst GRB 170817A provide us a significant opportunity to study theories of gravity beyond general relativity. An important outcome of these observations is that they constrain the difference between the speed of gravity and the speed of light to less than $10^{-15}c$. Also, the time delay between the arrivals of gravitational waves at different detectors constrains the speed of gravity at the Earth to be in the range $0.55c < v_{gw} < 1.42c$. We use these results to constrain a widely studied modified theory of gravity: Eddington-inspired Born-Infeld (EiBI) gravity. We show that, in EiBI theory, the speed of gravitational waves in matter deviates from $c$. From the time delay in arrival of gravitational wave signals at Earth-based detectors, we obtain the bound on the theory parameter $\kappa$ as $\vert\kappa\vert \lesssim 10^{21}\, m^2$. Similarly, from the time delay between the signals of GW170817 and GRB 170817A, in a background Friedmann-Robertson-Walker universe, we obtain $\vert \kappa \vert \lesssim 10^{37}\, m^2$. Although the bounds on $\kappa$ are weak compared to other earlier bounds from the study of neutron stars, stellar evolution, primordial nucleosynthesis, etc., our bounds are from the direct observations and thus worth noting.


title: Metal abundances in hot white dwarfs with signatures of a superionized wind

abstract: About a dozen hot white dwarfs with effective temperatures Teff = 65,000-120,000 K exhibit unusual absorption features in their optical spectra. These objects were tentatively identified as Rydberg lines of ultra-high excited metals in ionization stages V-X, indicating line formation in a dense environment with temperatures near one million Kelvin. Since some features show blueward extensions, it was argued that they stem from a superionized wind. A unique assignment of the lines to particular elements is not possible, although they probably stem from C, N, O, and Ne. To further investigate this phenomenon, we analyzed the ultraviolet spectra available from only three stars of this group; that is, two helium-rich white dwarfs, HE 0504-2408 and HS 0713+3958 with spectral type DO, and a hydrogen-rich white dwarf, HS 2115+1148 with spectral type DAO. We identified light metals (C, N, O, Si, P, and S) with generally subsolar abundances and heavy elements from the iron group (Cr, Mn, Fe, Co, Ni) with solar or oversolar abundance. The abundance patterns are not unusual for hot WDs and can be interpreted as the result of gravitational settling and radiative levitation of elements. As to the origin of the ultra-high ionized metals lines, we discuss the possible presence of a multicomponent radiatively driven wind that is frictionally heated.


title: On constraining projections of future climate using observations and simulations from multiple climate models

abstract: Numerical climate models are used to project future climate change due to both anthropogenic and natural causes. Differences between projections from different climate models are a major source of uncertainty about future climate. Emergent relationships shared by multiple climate models have the potential to constrain our uncertainty when combined with historical observations. We combine projections from 13 climate models with observational data to quantify the impact of emergent relationships on projections of future warming in the Arctic at the end of the 21st century. We propose a hierarchical Bayesian framework based on a coexchangeable representation of the relationship between climate models and the Earth system. We show how emergent constraints fit into the coexchangeable representation, and extend it to account for internal variability simulated by the models and natural variability in the Earth system. Our analysis shows that projected warming in some regions of the Arctic may be more than 2C lower and our uncertainty reduced by up to 30% when constrained by historical observations. A detailed theoretical comparison with existing multi-model projection frameworks is also provided. In particular, we show that projections may be biased if we do not account for internal variability in climate model predictions.


title: Surjectivity of Euler operators on temperate distributions

abstract: Euler operators are partial differential operators of the form $P(\theta)$ where $P$ is a polynomial and $\theta_j = x_j \partial/\partial x_j$. We show that every non-trivial Euler operator is surjective on the space of temperate distributions on $R^d$. This is in sharp contrast to the behaviour of such operators when acting on spaces of differentiable or analytic functions.


title: Truncated Polynomial Expansion Downlink Precoders and Uplink Detectors for Massive MIMO

abstract: In TDD reciprocity-based massive MIMO it is essential to be able to compute the downlink precoding matrix over all OFDM resource blocks within a small fraction of the uplink-downlink slot duration. Early implementation of massive MIMO are limited to the simple Conjugate Beamforming (ConjBF) precoding method, because of such computation latency limitation. However, it has been widely demonstrated by theoretical analysis and system simulation that Regularized Zero-Forcing (RZF) precoding is generally much more effective than ConjBF for a large but practical number of transmit antennas. In order to recover a significant fraction of the gap between ConjBF and RZF and yet meeting the very strict computation latency constraints, truncated polynomial expansion (TPE) methods have been proposed. In this paper we present a novel TPE method that outperforms all previously proposed methods in the general non-symmetric case of users with arbitrary antenna correlation. In addition, the proposed method is significantly simpler and more flexible than previously proposed methods based on deterministic equivalents and free probability in large random matrix theory. We consider power allocation with our TPE approach, and show that classical system optimization problems such as min-sum power and max-min rate can be easily solved. Furthermore, we provide a detailed computation latency analysis specifically targeted to a highly parallel FPGA hardware architecture.


title: Beurling's Theorem for the Two-sided Quaternion Fourier Transform

abstract: The two-sided quaternion Fourier transform satisfies some uncertainty principles similar to the Euclidean Fourier transform. A generalization of Beurling's theorem, Hardy, Cowling-Price and Gelfand-Shilov theorems, is obtained for the two-sided quaternion Fourier transform.


title: Entropic patchiness drives multi-phase coexistence in discotic colloid-depletant mixtures

abstract: Entropy-driven equilibrium phase behaviour of hard particle dispersions can be understood from excluded volume arguments only. While monodisperse hard spheres only exhibit a fluid-solid phase transition, anisotropic hard particles such as rods, discs, cuboids or boards exhibit various multi-phase equilibria. Ordering of such anisotropic particles increases the free volume entropy by reducing the excluded volume between them. The addition of depletants gives rise to an entropic patchiness represented by orientation-dependent attraction resulting in non-trivial phase behaviour. We show that free volume theory is a simple, generic and tractable framework that enables to incorporate these effects and rationalise various experimental findings. Plate-shaped particles constitute the main building blocks of clays, asphaltenes and chromonic liquid crystals that find widespread use in the food, cosmetics and oil industry. We demonstrate that mixtures of platelets and ideal depletants exhibit a strikingly rich phase behaviour containing several types of three-phase coexistence areas and even a quadruple region with four coexisting phases.


title: Boundedness of the Stationary Solution to the Boltzmann Equation with Spatial Smearing, Diffusive Boundary Conditions, and Lions' Collision Kernel

abstract: We investigate the Boltzmann equation with spatial smearing, diffusive boundary conditions, and Lions' collision kernel. Both, the physical as well as the velocity space, are assumed to be bounded. Existence and uniqueness of a stationary solution, which is a probability density, has been demonstrated in [3] under a certain smallness assumption on the collision term. We prove that whenever there is a stationary solution then it is a.e. positively bounded from below and above.


title: Minimax estimation in linear models with unknown design over finite alphabets

abstract: We provide a minimax optimal estimation procedure for F and W in matrix valued linear models Y = F W + Z where the parameter matrix W and the design matrix F are unknown but the latter takes values in a known finite set. The proposed finite alphabet linear model is justified in a variety of applications, ranging from signal processing to cancer genetics. We show that this allows to separate F and W uniquely under weak identifiability conditions, a task which is not doable, in general. To this end we quantify in the noiseless case, that is, Z = 0, the perturbation range of Y in order to obtain stable recovery of F and W. Based on this, we derive an iterative Lloyd's type estimation procedure that attains minimax estimation rates for W and F for Gaussian error matrix Z. In contrast to the least squares solution the estimation procedure can be computed efficiently and scales linearly with the total number of observations. We confirm our theoretical results in a simulation study and illustrate it with a genetic sequencing data example.


title: Quasinormal modes of black holes in Einstein-power-Maxwell theory

abstract: In the present work we compute the spectrum of quasinormal frequencies of four-dimensional charged black holes in Einstein-power-Maxwell theory. In particular we study scalar perturbations adopting the 6th order WKB approximation. We analyze in detail the behaviour of the spectrum depending on the charge of the black hole, the quantum number of angular momentum and the overtone number. In addition, a comparison is made between the results obtained here and the results valid for charged black holes in other theories as well as for the Reissner-Nordstr${\"o}$m black hole. Finally, we have provided an analytical expression for the quasinormal spectrum in the eikonal limit.


title: Deep Residual Text Detection Network for Scene Text

abstract: Scene text detection is a challenging problem in computer vision. In this paper, we propose a novel text detection network based on prevalent object detection frameworks. In order to obtain stronger semantic feature, we adopt ResNet as feature extraction layers and exploit multi-level feature by combining hierarchical convolutional networks. A vertical proposal mechanism is utilized to avoid proposal classification, while regression layer remains working to improve localization accuracy. Our approach evaluated on ICDAR2013 dataset achieves F-measure of 0.91, which outperforms previous state-of-the-art results in scene text detection.


title: A Riemann-Hilbert correspondence in positive characteristic

abstract: We explain a version of the Riemann-Hilbert correspondence for $p$-torsion \'etale sheaves on an arbitrary $\mathbf{F}_p$-scheme.


title: Broadcast in radio networks: time vs. energy tradeoffs

abstract: In wireless networks, consisting of battery-powered devices, energy is a costly resource and most of it is spent on transmitting and receiving messages. Broadcast is a problem where a message needs to be transmitted from one node to all other nodes of the network. We study algorithms that can work under limited energy measured as the maximum number of transmissions by a single station. The goal of the paper is to study tradeoffs between time and energy complexity of broadcast problem in multi-hop radio networks. We consider a model where the topology of the network is unknown and if two neighbors of a station are transmitting in the same discrete time slot, then the signals collide and the receiver cannot distinguish the collided signals from silence. We observe that existing, time efficient, algorithms are not optimized with respect to energy expenditure. We then propose and analyse two new randomized energy-efficient algorithms. Our first algorithm works in time $O((D+\varphi)\cdot n^{1/\varphi}\cdot \varphi)$ with high probability and uses $O(\varphi)$ energy per station for any $\varphi \leq \log n/(2\log\log n)$ for any graph with $n$ nodes and diameter $D$. Our second algorithm works in time $O((D+\log n)\log n)$ with high probability and uses $O(\log n/\log\log n)$ energy. We prove that our algorithms are almost time-optimal for given energy limits for graphs with constant diameters by constructing lower bound on time of $\Omega(n^{1/\varphi} \cdot \varphi)$. The lower bound shows also that any algorithm working in polylogaritmic time in $n$ for all graphs needs energy $\Omega(\log n/\log\log n)$.


title: STWalk: Learning Trajectory Representations in Temporal Graphs

abstract: Analyzing the temporal behavior of nodes in time-varying graphs is useful for many applications such as targeted advertising, community evolution and outlier detection. In this paper, we present a novel approach, STWalk, for learning trajectory representations of nodes in temporal graphs. The proposed framework makes use of structural properties of graphs at current and previous time-steps to learn effective node trajectory representations. STWalk performs random walks on a graph at a given time step (called space-walk) as well as on graphs from past time-steps (called time-walk) to capture the spatio-temporal behavior of nodes. We propose two variants of STWalk to learn trajectory representations. In one algorithm, we perform space-walk and time-walk as part of a single step. In the other variant, we perform space-walk and time-walk separately and combine the learned representations to get the final trajectory embedding. Extensive experiments on three real-world temporal graph datasets validate the effectiveness of the learned representations when compared to three baseline methods. We also show the goodness of the learned trajectory embeddings for change point detection, as well as demonstrate that arithmetic operations on these trajectory representations yield interesting and interpretable results.


title: Second order splitting for a class of fourth order equations

abstract: We formulate a well-posedness and approximation theory for a class of generalised saddle point problems. In this way we develop an approach to a class of fourth order elliptic partial differential equations using the idea of splitting into coupled second order equations. Our main motivation is to treat certain fourth order surface equations arising in the modelling of biomembranes but the approach may be applied more generally. In particular, we are interested in equations with non-smooth right hand sides and operators which have non-trivial kernels.The theory for well posedness and approximation is presented in an abstract setting. Several examples are described together with some numerical experiments.


title: No static sphericaly symmetric wormholes in Horndeski theory

abstract: We consider the Horndeski theory in four-dimensional space-time. We show that this theory does not admit stable, static, spherically symmetric, asymptotically flat, Lorentzian wormholes.


title: Inverse Design of Simple Pair Potentials for the Self-Assembly of Complex Structures

abstract: The synthesis of complex materials through the self-assembly of particles at the nanoscale provides opportunities for the realization of novel material properties. However, the inverse design process to create experimentally feasible interparticle interaction strategies is uniquely challenging. Standard methods for the optimization of isotropic pair potentials tend toward overfitting, resulting in solutions with too many features and length scales that are challenging to map to mechanistic models. Here we introduce a method for the optimization of simple pair potentials that minimizes the relative entropy of the complex target structure while directly considering only those length scales most relevant for self-assembly. Our approach maximizes the relative information of a target pair distribution function with respect to an ansatz distribution function via an iterative update process. During this process, we filter high frequencies from the Fourier spectrum of the pair potential, resulting in interaction potentials that are smoother and simpler in real space, and therefore likely easier to make. We show that pair potentials obtained by this method assemble their target structure more robustly with respect to optimization method parameters than potentials optimized without filtering.


title: Interpretable probabilistic embeddings: bridging the gap between topic models and neural networks

abstract: We consider probabilistic topic models and more recent word embedding techniques from a perspective of learning hidden semantic representations. Inspired by a striking similarity of the two approaches, we merge them and learn probabilistic embeddings with online EM-algorithm on word co-occurrence data. The resulting embeddings perform on par with Skip-Gram Negative Sampling (SGNS) on word similarity tasks and benefit in the interpretability of the components. Next, we learn probabilistic document embeddings that outperform paragraph2vec on a document similarity task and require less memory and time for training. Finally, we employ multimodal Additive Regularization of Topic Models (ARTM) to obtain a high sparsity and learn embeddings for other modalities, such as timestamps and categories. We observe further improvement of word similarity performance and meaningful inter-modality similarities.


title: Deterministic parallel analysis: An improved method for selecting factors and principal components

abstract: Factor analysis and principal component analysis (PCA) are used in many application areas. The first step, choosing the number of components, remains a serious challenge. Our work proposes improved methods for this important problem. One of the most popular state-of-the-art methods is Parallel Analysis (PA), which compares the observed factor strengths to simulated ones under a noise-only model. This paper proposes improvements to PA. We first de-randomize it, proposing Deterministic Parallel Analysis (DPA), which is faster and more reproducible than PA. Both PA and DPA are prone to a shadowing phenomenon in which a strong factor makes it hard to detect smaller but more interesting factors. We propose deflation to counter shadowing. We also propose to raise the decision threshold to improve estimation accuracy. We prove several consistency results for our methods, and test them in simulations. We also illustrate our methods on data from the Human Genome Diversity Project, where they significantly improve the accuracy.