title: Gravitational entropy and the cosmological no-hair conjecture

abstract: The gravitational entropy and no-hair conjectures seems to predict contradictory future states of our Universe. The growth of the gravitational entropy is associated with the growth of inhomogeneity, while the no-hair conjecture argues that a universe dominated by dark energy should asymptotically approach a homogeneous and isotropic de Sitter state. The aim of this paper is to study these two conjectures. The investigation is based on the Simsilun simulation, which simulates the universe using the approximation of the Silent Universe. The Silent Universe is a solution to the Einstein equations that assumes irrotational, non-viscous, and insulated dust, with vanishing magnetic part of the Weyl curvature. The initial conditions for the Simsilun simulation are sourced from the Millennium simulation, which results with a realistically appearing but relativistic at origin simulation of a universe. The Simsilun simulation is evolved from the early universe (t = 25 Myr) till far future (t = 1000 Gyr). The results of this investigation show that both conjectures are correct. On global scales, a universe with a positive cosmological constant and non-positive spatial curvature does indeed approach the de Sitter state. At the same time it keeps generating the gravitational entropy.


title: The California-Kepler Survey. IV. Metal-rich Stars Host a Greater Diversity of Planets

abstract: Probing the connection between a star's metallicity and the presence and properties of any associated planets offers an observational link between conditions during the epoch of planet formation and mature planetary systems. We explore this connection by analyzing the metallicities of Kepler target stars and the subset of stars found to host transiting planets. After correcting for survey incompleteness, we measure planet occurrence: the number of planets per 100 stars with a given metallicity $M$. Planet occurrence correlates with metallicity for some, but not all, planet sizes and orbital periods. For warm super-Earths having $P = 10-100$ days and $R_P = 1.0-1.7~R_E$, planet occurrence is nearly constant over metallicities spanning $-$0.4 dex to +0.4 dex. We find 20 warm super-Earths per 100 stars, regardless of metallicity. In contrast, the occurrence of warm sub-Neptunes ($R_P = 1.7-4.0~R_E$) doubles over that same metallicity interval, from 20 to 40 planets per 100 stars. We model the distribution of planets as $d f \propto 10^{\beta M} d M$, where $\beta$ characterizes the strength of any metallicity correlation. This correlation steepens with decreasing orbital period and increasing planet size. For warm super-Earths $\beta = -0.3^{+0.2}_{-0.2}$, while for hot Jupiters $\beta = +3.4^{+0.9}_{-0.8}$. High metallicities in protoplanetary disks may increase the mass of the largest rocky cores or the speed at which they are assembled, enhancing the production of planets larger than 1.7 $R_E$. The association between high metallicity and short-period planets may reflect disk density profiles that facilitate the inward migration of solids or higher rates of planet-planet scattering.


title: How to navigate through obstacles?

abstract: Given a set of obstacles and two points, is there a path between the two points that does not cross more than $k$ different obstacles? This is a fundamental problem that has undergone a tremendous amount of work. It is known to be NP-hard, even when the obstacles are very simple geometric shapes (e.g., unit-length line segments). The problem can be generalized into the following graph problem: Given a planar graph $G$ whose vertices are colored by color sets, two designated vertices $s, t \in V(G)$, and $k \in \mathbb{N}$, is there an $s$-$t$ path in $G$ that uses at most $k$ colors? If each obstacle is connected, the resulting graph satisfies the color-connectivity property, namely that each color induces a connected subgraph. We study the complexity and design algorithms for the above graph problem with an eye on its geometric applications. We prove that without the color-connectivity property, the problem is W[SAT]-hard parameterized by $k$. A corollary of this result is that, unless W[2] $=$ FPT, the problem cannot be approximated in FPT time to within a factor that is a function of $k$. By describing a generic plane embedding of the graph instances, we show that our hardness results translate to the geometric instances of the problem. We then focus on graphs satisfying the color-connectivity property. By exploiting the planarity of the graph and the connectivity of the colors, we develop topological results to "represent" the valid $s$-$t$ paths containing subsets of colors from any vertex $v$. We employ these results to design an FPT algorithm for the problem parameterized by both $k$ and the treewidth of the graph, and extend this result to obtain an FPT algorithm for the parameterization by both $k$ and the length of the path. The latter result directly implies previous FPT results for various obstacle shapes, such as unit disks and fat regions.


title: Recursive computation of the invariant distributions of Feller processes: Revisited examples and new applications

abstract: In this paper, we show that the abstract framework developed in Pages & Rey (2017) and inspired by Lamberton & Pages (2002) can be used to build invariant distributions for Brownian diffusion processes using the Milstein scheme and for diffusion processes with censored jump using the Euler scheme. Both studies rely on a weakly mean reverting setting for both cases. For the Milstein scheme we prove the convergence for test functions with polynomial (Wasserstein convergence) and exponential growth. For the Euler scheme of diffusion processes with censored jump we prove the convergence for test functions with polynomial growth.


title: Choose your path wisely: gradient descent in a Bregman distance framework

abstract: We propose an extension of a special form of gradient descent -- in the literature known as linearised Bregman iteration -- to a larger class of non-convex functions. We replace the classical (squared) two norm metric in the gradient descent setting with a generalised Bregman distance, based on a proper, convex and lower semi-continuous function. The algorithm's global convergence is proven for functions that satisfy the Kurdyka-\L ojasiewicz property. Examples illustrate that features of different scale are being introduced throughout the iteration, transitioning from coarse to fine. This coarse-to-fine approach with respect to scale allows to recover solutions of non-convex optimisation problems that are superior to those obtained with conventional gradient descent, or even projected and proximal gradient descent. The effectiveness of the linearised Bregman iteration in combination with early stopping is illustrated for the applications of parallel magnetic resonance imaging, blind deconvolution as well as image classification with neural networks.


title: Character-Based Handwritten Text Transcription with Attention Networks

abstract: The paper approaches the task of handwritten text recognition (HTR) with attentional encoder-decoder networks trained on sequences of characters, rather than words. We experiment on lines of text from popular handwriting datasets and compare different activation functions for the attention mechanism used for aligning image pixels and target characters. We find that softmax attention focuses heavily on individual characters, while sigmoid attention focuses on multiple characters at each step of the decoding. When the sequence alignment is one-to-one, softmax attention is able to learn a more precise alignment at each step of the decoding, whereas the alignment generated by sigmoid attention is much less precise. When a linear function is used to obtain attention weights, the model predicts a character by looking at the entire sequence of characters and performs poorly because it lacks a precise alignment between the source and target. Future research may explore HTR in natural scene images, since the model is capable of transcribing handwritten text without the need for producing segmentations or bounding boxes of text in images.


title: Krylov integrators for Hamiltonian systems

abstract: We consider Arnoldi like processes to obtain symplectic subspaces for Hamiltonian systems. Large systems are locally approximated by ones living in low dimensional subspaces; we especially consider Krylov subspaces and some extensions. This will be utilized in two ways: solve numerically local small dimensional systems or in a given numerical, e.g. exponential, integrator, use the subspace for approximations of necessary functions. In the former case one can expect an excellent energy preservation. For the latter this is so for linear systems. For some second order exponential integrators we consider these two approaches are shown to be equivalent. In numerical experiments with nonlinear Hamiltonian problems their behaviour seems promising.


title: Cavs: A Vertex-centric Programming Interface for Dynamic Neural Networks

abstract: Recent deep learning (DL) models have moved beyond static network architectures to dynamic ones, handling data where the network structure changes every example, such as sequences of variable lengths, trees, and graphs. Existing dataflow-based programming models for DL---both static and dynamic declaration---either cannot readily express these dynamic models, or are inefficient due to repeated dataflow graph construction and processing, and difficulties in batched execution. We present Cavs, a vertex-centric programming interface and optimized system implementation for dynamic DL models. Cavs represents dynamic network structure as a static vertex function $\mathcal{F}$ and a dynamic instance-specific graph $\mathcal{G}$, and performs backpropagation by scheduling the execution of $\mathcal{F}$ following the dependencies in $\mathcal{G}$. Cavs bypasses expensive graph construction and preprocessing overhead, allows for the use of static graph optimization techniques on pre-defined operations in $\mathcal{F}$, and naturally exposes batched execution opportunities over different graphs. Experiments comparing Cavs to two state-of-the-art frameworks for dynamic NNs (TensorFlow Fold and DyNet) demonstrate the efficacy of this approach: Cavs achieves a near one order of magnitude speedup on training of various dynamic NN architectures, and ablations demonstrate the contribution of our proposed batching and memory management strategies.


title: Diffuse X-ray-emitting Gas in Major Mergers

abstract: Using archived data from the Chandra X-ray telescope, we have extracted the diffuse X-ray emission from 49 equal-mass interacting/merging galaxy pairs in a merger sequence, from widely separated pairs to merger remnants. After removal of contributions from unresolved point sources, we compared the diffuse thermal X-ray luminosity from hot gas (L(X)(gas)) with the global star formation rate (SFR). After correction for absorption within the target galaxy, we do not see strong trend of L(X)(gas)/SFR with SFR or merger stage for galaxies with SFR > 1 M(sun) yr^-1. For these galaxies, the median L(X)(gas)/SFR is 5.5 X 10^39 ((erg s^-1)/M(sun) yr^-1)), similar to that of normal spiral galaxies. These results suggest that stellar feedback in star forming galaxies reaches an approximately steady state condition, in which a relatively constant fraction of about 2% of the total energy output from supernovae and stellar winds is converted into X-ray flux. Three late-stage merger remnants with low SFRs and high K band luminosities (L(K)) have enhanced L(X)(gas)/SFR; their UV/IR/optical colors suggest that they are post-starburst galaxies, perhaps in the process of becoming ellipticals. Systems with L(K) < 10^10 L(sun) have lower L(X)(gas)/SFR ratios than the other galaxies in our sample, perhaps due to lower gravitational fields or lower metallicities. We see no relation between L(X)(gas)/SFR and Seyfert activity in this sample, suggesting that feedback from active galactic nuclei is not a major contributor to the hot gas in our sample galaxies.


title: Jet-medium interaction and the Gubser flow

abstract: We study the effects of expansion and viscous corrections on the hydrodynamical medium response to a high energy jet parton. More specifically, using a semi-analytical Gubser solution to relativistic fluid dynamics, modifications to the formed Mach cone, diffusive wake, and the momentum flow of the medium response along and perpendicular to the jet particle, are analyzed mode-by-mode. This should provide intuition and guidance for analyses of the experimentally measured jet sub-structure in heavy-ion collisions.


title: Irreducible Polynomials with Varying Constraints on Coefficients

abstract: We study the number of prime polynomials of degree $n$ over $\mathbb{F}_q$ in which the $i^{th}$ coefficient is either preassigned to be $a_i \in \mathbb{F}_q$ or outside a small set $S_i \subset \mathbb{F}_q$. This serves as a function field analogue of a recent work of Maynard, which counts integer primes that do not have specific digits in their base-$q$ expansion. Our work relates to Pollack's and Ha's work, which count the amount of prime polynomials with $\ll \sqrt{n}$ and $\ll n$ preassigned coefficients, respectively. Our result demonstrates how one can prove asymptotics of the number of prime polynomials with different types of constraints to each coefficient.


title: On the Lipman-Zariski conjecture for logarithmic vector fields on log canonical pairs

abstract: We consider a version of the Lipman-Zariski conjecture for logarithmic vector fields and logarithmic $1$-forms on pairs. Let $(X,D)$ be a pair consisting of a normal complex variety $X$ and an effective Weil divisor $D$ such that the sheaf of logarithmic vector fields (or dually the sheaf of reflexive logarithmic $1$-forms) is locally free. We prove that in this case the following holds: If $(X,D)$ is dlt, then $X$ is necessarily smooth and $\lfloor D\rfloor $ is snc. If $(X,D)$ is lc or the logarithmic $1$-forms are locally generated by closed forms, then $(X,\lfloor D\rfloor)$ is toroidal.


title: Cascading Failures as Continuous Phase-Space Transitions

abstract: In network systems, a local perturbation can amplify as it propagates, potentially leading to a large-scale cascading failure. Here we derive a continuous model to advance our understanding of cascading failures in power-grid networks. The model accounts for both the failure of transmission lines and the desynchronization of power generators, and incorporates the transient dynamics between successive steps of the cascade. In this framework, we show that a cascade event is a phase-space transition from an equilibrium state with high energy to an equilibrium state with lower energy, which can be suitably described in closed form using a global Hamiltonian-like function. From this function we show that a perturbed system cannot always reach the equilibrium state predicted by quasi-steady-state cascade models, which would correspond to a reduced number of failures, and may instead undergo a larger cascade. We also show that in the presence of two or more perturbations, the outcome depends strongly on the order and timing of the individual perturbations. These results offer new insights into the current understanding of cascading dynamics, with potential implications for control interventions.


title: Information Dissemination Speed in Delay Tolerant Urban Vehicular Networks in a Hyperfractal Setting

abstract: This paper studies the fundamental communication properties of urban vehicle networks by exploiting the self-similarity and hierarchical organization of modern cities. We use an innovative model called "hyperfractal" that captures the self-similarities of both the traffic and vehicle locations but avoids the extremes of regularity and randomness. We use analytical tools to derive theoretical upper and lower bounds for the information propagation speed in an urban delay tolerant network (i.e., a network that is disconnected at all time, and thus uses a store-carry-and-forward routing model). We prove that the average broadcast time behaves as $n^{1-\delta}$ times a slowly varying function, where $\delta$ depends on the precise fractal dimension. Furthermore, we show that the broadcast speedup is due in part to an interesting self-similar phenomenon, that we denote as {\em information teleportation}. This phenomenon arises as a consequence of the topology of the vehicle traffic, and triggers an acceleration of the broadcast time. We show that our model fits real cities where open traffic data sets are available. We present simulations confirming the validity of the bounds in multiple realistic settings, including scenarios with variable speed, using both QualNet and a discrete-event simulator in Matlab.


title: Many-body synchronization of interacting qubits by engineered ac-driving

abstract: In this work we introduce the many-body synchronization of an interacting qubit ensemble which allows one to switch dynamically from many-body-localized (MBL) to an ergodic state. We show that applying of $\pi$-pulses with altering phases, one can effectively suppress the MBL phase and, hence, eliminate qubits disorder. The findings are based on the analysis of the Loschmidt echo dynamics which shows a transition from a power-law decay to more rapid one indicating the dynamical MBL-to-ergodic transition. The technique does not require to know the microscopic details of the disorder.

categories: cond-mat.dis-nn cond-mat.mes-hall cond-mat.str-el quant-ph


title: Isentropic thermal instability in atomic surface layers of photodissociation regions

abstract: We consider the evolution of an isentropic thermal instability in the atomic zone of a photodissociation region (PDR). In this zone, gas heating and cooling are associated mainly with photoelectric emission from dust grains and fine-structure lines ([\ion{C}{ii}] 158, [\ion{O}{i}] 63, and [\ion{O}{i}] 146 {\micron}), respectively. The instability criterion has a multi-parametric dependence on the conditions of the interstellar medium. We found that instability occurs when the intensity of the incident far-ultraviolet field $G_0$ and gas density $n$ are high. For example, we have $3\times10^3<G_0<10^6$ and $4.5\times10^4<n<10^6$ {\cmc} at temperatures $360 <T<10^4$ K for typical carbon and oxygen abundances $\xi_{\rm C}=1.4\times10^{-4}$ and $\xi_{\rm O}=3.2\times10^{-4}$. The instability criterion depends on the relation between $\xi_{\rm C}$ and $\xi_{\rm O}$ abundances and line opacities. We also give examples of observed PDRs where instability could occur. For these PDRs, the characteristic perturbation growth time is $t_{\rm inst}\sim10^3$ -- $10^4$ yr and the distance characterizing the formation of secondary waves is $L\sim10^{-3}$ -- $5\times10^{-2}$ pc. For objects that are older than $t_{\rm inst}$ and have sizes of the atomic zone larger than $L$, we expect that instability influences the PDR structure significantly. The presence of multiple shock waves, turbulent velocities of several kilometers per second and inhomogeneities with higher density and temperature than the surrounding medium can characterize isentropic thermal instability in PDRs.


title: World-line formulation of chiral kinetic theory in topological background gauge fields

abstract: In heavy-ion collisions, an interesting question of phenomenological relevance is how the chiral imbalance generated at early times persists through a fluctuating background of sphalerons in addition to other "non-anomalous" interactions with the QGP. To address this question, we construct a relativistic chiral kinetic theory using the world-line formulation of quantum field theory. We outline how Berry's phase arises in this framework, and how its effects can be clearly distinguished from those arising from the chiral anomaly. We further outline how this framework can be matched to classical statistical simulations at early times and to anomalous chiral hydrodynamics at late times.


title: Distinguishing differential susceptibility, diathesis-stress and vantage sensitivity: beyond the single gene and environment model

abstract: Currently, two main approaches exist to distinguish differential susceptibility from diathesis-stress and vantage sensitivity in genotype x environment interaction (GxE) research: Regions of significance (RoS) and competitive-confirmatory approaches. Each is limited by their single-gene/single-environment foci given that most phenotypes are the product of multiple interacting genetic and environmental factors. We thus addressed these two concerns in a recently developed R package (LEGIT) for constructing GxE interaction models with latent genetic and environmental scores using alternating optimization. Herein we test, by means of computer simulation, diverse GxE models in the context of both single and multiple genes and environments. Results indicate that the RoS and competitive-confirmatory approaches were highly accurate when the sample size was large, whereas the latter performed better in small samples and for small effect sizes. The confirmatory approach generally had good accuracy (a) when effect size was moderate and N >= 500 and (b) when effect size was large and N >= 250, whereas RoS performed poorly. Computational tools to determine the type of GxE of multiple genes and environments are provided as extensions in our LEGIT R package.


title: Optimal Joint Routing and Scheduling in Millimeter-Wave Cellular Networks

abstract: Millimeter-wave (mmWave) communication is a promising technology to cope with the expected exponential increase in data traffic in 5G networks. mmWave networks typically require a very dense deployment of mmWave base stations (mmBS). To reduce cost and increase flexibility, wireless backhauling is needed to connect the mmBSs. The characteristics of mmWave communication, and specifically its high directional- ity, imply new requirements for efficient routing and scheduling paradigms. We propose an efficient scheduling method, so-called schedule-oriented optimization, based on matching theory that optimizes QoS metrics jointly with routing. It is capable of solving any scheduling problem that can be formulated as a linear program whose variables are link times and QoS metrics. As an example of the schedule-oriented optimization, we show the optimal solution of the maximum throughput fair scheduling (MTFS). Practically, the optimal scheduling can be obtained even for networks with over 200 mmBSs. To further increase the runtime performance, we propose an efficient edge-coloring based approximation algorithm with provable performance bound. It achieves over 80% of the optimal max-min throughput and runs 5 to 100 times faster than the optimal algorithm in practice. Finally, we extend the optimal and approximation algorithms for the cases of multi-RF-chain mmBSs and integrated backhaul and access networks.


title: A multi-parameter variant of the Erd\H{o}s distance problem

abstract: We study the following variant of the Erd\H{o}s distance problem. Given $E$ and $F$ a point sets in $\mathbb{R}^d$ and $p = (p_1, \ldots, p_q)$ with $p_1+ \cdots + p_q = d$ is an increasing partition of $d$ define $$ B_p(E,F)=\{(|x_1-y_1|, \ldots, |x_q-y_q|): x \in E, y \in F \},$$ where $x=(x_1, \ldots, x_q)$ with $x_i$ in $\mathbb{R}^{p_i}$. For $p_1 \geq 2$ it is not difficult to construct $E$ and $F$ such that $|B_{p}(E,F)|=1$. On the other hand, it is easy to see that if $\gamma_q$ is the best know exponent for the distance problem in $\mathbb{R}^{p_i}$ that $|B_p(E,E)| \geq C{|E|}^{\frac{\gamma_q}{q}}$. The question we study is whether we can improve the exponent $\frac{\gamma_q}{q}$. We first study partitions of length two in detail and prove the optimal result (up to logarithms) that $$ |B_{2,2}(E)| \gtrapprox |E|.$$ In the generalised two dimensional case for $B_{k,l}$ we need the stronger condition that $E$ is $s$-adaptable for $s<\frac{k}{2}+\frac{1}{3}$, letting $\gamma_m$ be the best known exponent for the Erd\H{o}s-distance problem in $\mathbb{R}^m$ for $k \neq l$ we gain a further optimal result of, $$ |B_{k,l}(E)| \gtrapprox |E|^{\gamma_l}.$$ When $k=l$ we use the explicit $\gamma_m=\frac{m}{2}-\frac{2}{m(m+2)}$ result due to Solymosi and Vu to gain $$ |B_{k,k}(E)| \gtrapprox |E|^{\frac{13}{14}\gamma_k}.$$ For a general partition, let $\gamma_i = \frac{2}{p_i}-\frac{2}{p_i(p_i+2)}$ and $\eta_i = \frac{2}{2d-(p_i-1)}$. Then if $E$ is $s$-adaptable with $s>d-\frac{p_1}{2}+\frac{1}{3}$ we have $$ B_p(E) \gtrapprox |E|^\tau \hspace{0.5cm} \text{where} \hspace{0.5cm} \tau = \gamma_q\left(\frac{\gamma_1+\eta_1}{\gamma_q+(q-1)(\gamma_1+\eta_1)}\right).$$ Where $p_i \sim \frac{d}{q}$ implies $\tau \sim \gamma_{q}\left(\frac{1}{q}+\frac{1}{dq}\right)$ and $p_q \sim d$ (with $q<<d$) implies $\tau \sim \gamma_{q}\left(\frac{1}{q}+\frac{1}{q^2}\right)$.


title: Local boundedness of solutions to nonlocal equations modeled on the fractional p-Laplacian

abstract: We state and prove estimates for the local boundedness of subsolutions of non-local, possibly degenerate, parabolic integro-differential equations of the form \begin{equation*} \partial_tu(x,t)+\mbox{P.V.}\int\limits_{\mathbb R^n}K(x,y,t) |u(x,t)-u(y,t) |^{p-2}(u(x,t)-u(y,t))\, dy,\end{equation*} $(x,t)\in\mathbb R^n\times\mathbb R$, where $\mbox{P.V.} $ means in the principle value sense, $p\in (1,\infty)$ and the kernel obeys $K(x,y,t)\approx |x-y |^{n+ps}$ for some $s\in (0,1)$, uniformly in $(x,y,t)\in\mathbb R^n\times \mathbb R^n\times\mathbb R$.


title: Distributed Bayesian Filtering using Logarithmic Opinion Pool for Dynamic Sensor Networks

abstract: The discrete-time Distributed Bayesian Filtering (DBF) algorithm is presented for the problem of tracking a target dynamic model using a time-varying network of heterogeneous sensing agents. In the DBF algorithm, the sensing agents combine their normalized likelihood functions in a distributed manner using the logarithmic opinion pool and the dynamic average consensus algorithm. We show that each agent's estimated likelihood function globally exponentially converges to an error ball centered on the joint likelihood function of the centralized multi-sensor Bayesian filtering algorithm. We rigorously characterize the convergence, stability, and robustness properties of the DBF algorithm. Moreover, we provide an explicit bound on the time step size of the DBF algorithm that depends on the time-scale of the target dynamics, the desired convergence error bound, and the modeling and communication error bounds. Furthermore, the DBF algorithm for linear-Gaussian models is cast into a modified form of the Kalman information filter. The performance and robust properties of the DBF algorithm are validated using numerical simulations.


title: Characterizing lab instructors' self-reported learning goals to inform development of an experimental modeling skills assessment

abstract: The ability to develop, use, and refine models of experimental systems is a nationally recognized learning outcome for undergraduate physics lab courses. However, no assessments of students' model-based reasoning exist for upper-division labs. This study is the first step toward development of modeling assessments for optics and electronics labs. In order to identify test objectives that are likely relevant across many institutional contexts, we interviewed 35 lab instructors about the ways they incorporate modeling in their course learning goals and activities. The study design was informed by the Modeling Framework for Experimental Physics. This framework conceptualizes modeling as consisting of multiple subtasks: making measurements, constructing system models, comparing data to predictions, proposing causes for discrepancies, and enacting revisions to models or apparatus. We found that each modeling subtask was identified by multiple instructors as an important learning outcome for their course. Based on these results, we argue that test objectives should include probing students' competence with most modeling subtasks, and test items should be designed to elicit students' justifications for choosing particular modeling pathways. In addition to discussing these and other implications for assessment, we also identify future areas of research related to the role of modeling in optics and electronics labs.


title: Extracting Persistent Clusters in Dynamic Data via M\"obius inversion

abstract: Identifying and representing clusters in time-varying network data is of particular importance when studying collective behaviors emerging in nature, in mobile device networks or in social networks. Based on combinatorial, categorical, and persistence theoretic viewpoints, we establish a stable functorial pipeline for the summarization of the evolution of clusters in a time-varying network. We first construct a complete summary of the evolution of clusters in a given time-varying network over a set of entities $X$ of which takes the form of a formigram. This formigram can be understood as a certain Reeb graph $\mathcal{R}$ which is labeled by subsets of $X$. By applying M\"obius inversion to the formigram in two different manners, we obtain two dual notions of diagram: the maximal group diagram and the persistence clustergram, both of which are in the form of an `annotated' barcode. The maximal group diagram consists of time intervals annotated by their corresponding maximal groups -- a notion due to Buchin et al., implying that we recognize the notion of maximal groups as a special instance of generalized persistence diagram by Patel. On the other hand, the persistence clustergram is mostly obtained by annotating the intervals in the zigzag barcode of the Reeb graph $\mathcal{R}$ with certain merging/disbanding events in the given time-varying network. We show that both diagrams are complete invariants of formigrams (or equivalently of trajectory grouping structure by Buchin et al.) and thus contain more information than the Reeb graph $\mathcal{R}$.


title: The Eigenoption-Critic Framework

abstract: Eigenoptions (EOs) have been recently introduced as a promising idea for generating a diverse set of options through the graph Laplacian, having been shown to allow efficient exploration. Despite its initial promising results, a couple of issues in current algorithms limit its application, namely: (1) EO methods require two separate steps (eigenoption discovery and reward maximization) to learn a control policy, which can incur a significant amount of storage and computation; (2) EOs are only defined for problems with discrete state-spaces and; (3) it is not easy to take the environment's reward function into consideration when discovering EOs. To addresses these issues, we introduce an algorithm termed eigenoption-critic (EOC) based on the Option-critic (OC) framework [Bacon17], a general hierarchical reinforcement learning (RL) algorithm that allows learning the intra-option policies simultaneously with the policy over options. We also propose a generalization of EOC to problems with continuous state-spaces through the Nystr\"om approximation. EOC can also be seen as extending OC to nonstationary settings, where the discovered options are not tailored for a single task.


title: DIAPHANE: a Portable Radiation Transport Library for Astrophysical Applications

abstract: One of the most computationally demanding aspects of the hydrodynamical modelling of Astrophysical phenomena is the transport of energy by radiation or relativistic particles. Physical processes involving energy transport are ubiquitous and of capital importance in many scenarios ranging from planet formation to cosmic structure evolution, including explosive events like core collapse supernova or gamma-ray bursts. Moreover, the ability to model and hence understand these processes has often been limited by the approximations and incompleteness in the treatment of radiation and relativistic particles. The DIAPHANE project has focused in developing a portable and scalable library that handles the transport of radiation and particles (in particular neutrinos) independently of the underlying hydrodynamic code. In this work, we present the computational framework and the functionalities of the first version of the DIAPHANE library, which has been successfully ported to three different smoothed-particle hydrodynamic codes, GADGET2, GASOLINE and SPHYNX. We also present validation of different modules solving the equations of radiation and neutrino transport using different numerical schemes.


title: Higher transitive quantum groups: theory and models

abstract: We investigate the notion of $k$-transitivity for the quantum permutation groups $G\subset S_N^+$, with a brief review of the known $k=1,2$ results, and with a study of what happens at $k\geq3$. We discuss then matrix modelling questions for the algebras $C(G)$, notably by introducing the related notions of double and triple flat matrix model. At the level of the examples, our main results concern the quantum groups coming from the complex Hadamard matrices, and from the Weyl matrices.


title: On radial Schroedinger operators with a Coulomb potential

abstract: This paper presents a thorough analysis of 1-dimensional Schroedinger operators whose potential is a linear combination of the Coulomb term 1/r and the centrifugal term 1/r^2. We allow both coupling constants to be complex. Using natural boundary conditions at 0, a two parameter holomorphic family of closed operators is introduced. We call them the Whittaker operators, since in the mathematical literature their eigenvalue equation is called the Whittaker equation. Spectral and scattering theory for Whittaker operators is studied. Whittaker operators appear in quantum mechanics as the radial part of the Schroedinger operator with a Coulomb potential.


title: Atmospheric circulation, chemistry, and infrared spectra of Titan-like exoplanets around different stellar types

abstract: With the discovery of ever smaller and colder exoplanets, terrestrial worlds with hazy atmospheres must be increasingly considered. Our Solar System's Titan is a prototypical hazy planet, whose atmosphere may be representative of a large number of planets in our Galaxy. As a step towards characterizing such worlds, we present simulations of exoplanets that resemble Titan, but orbit three different stellar hosts: G-, K-, and M-dwarf stars. We use general circulation and photochemistry models to explore the circulation and chemistry of these Titan-like planets under varying stellar spectra, in all cases assuming a Titan-like insolation. Due to the strong absorption of visible light by atmospheric haze, the redder radiation accompanying later stellar types produces more isothermal stratospheres, stronger meridional temperature gradients at mbar pressures, and deeper and stronger zonal winds. In all cases, the planets' atmospheres are strongly superrotating, but meridional circulation cells are weaker aloft under redder starlight. The photochemistry of hydrocarbon and nitrile species varies with stellar spectra, with variations in the FUV/NUV flux ratio playing an important role. Our results tentatively suggest that column haze production rates could be similar under all three hosts, implying that planets around many different stars could have similar characteristics to Titan's atmosphere. Lastly, we present theoretical emission spectra. Overall, our study indicates that, despite important and subtle differences, the circulation and chemistry of Titan-like exoplanets are relatively insensitive to differences in host star. These findings may be further probed with future space-based facilities, like WFIRST, LUVOIR, HabEx, and OST.


title: Tail asymptotics of light-tailed Weibull-like sums

abstract: We consider sums of $n$ i.i.d. random variables with tails close to $\exp\{-x^\beta\}$ for some $\beta>1$. Asymptotics developed by Rootz\'en (1987) and Balkema, Kl\"uppelberg & Resnick (1993) are discussed from the point of view of tails rather of densities, using a somewhat different angle, and supplemented with bounds, results on a random number $N$ of terms, and simulation algorithms.


title: Halo Substructure in the SDSS-Gaia Catalogue : Streams and Clumps

abstract: We use the SDSS-Gaia Catalogue to identify six new pieces of halo substructure. SDSS-Gaia is an astrometric catalogue that exploits SDSS data release 9 to provide first epoch photometry for objects in the Gaia source catalogue. We use a version of the catalogue containing $245\,316$ stars with all phase space coordinates within a heliocentric distance of $\sim 10$ kpc. We devise a method to assess the significance of halo substructures based on their clustering in velocity space. The two most substantial structures are multiple wraps of a stream which has undergone considerable phase mixing (S1, with 94 members) and a kinematically cold stream (S2, with 61 members). The member stars of S1 have a median position of ($X,Y,Z$) = ($8.12, -0.22, 2.75$) kpc and a median metallicity of [Fe/H] $= -1.78$. The stars of S2 have median coordinates ($X,Y,Z$) = ($8.66, 0.30, 0.77$) kpc and a median metallicity of [Fe/H] $= -1.91$. They lie in velocity space close to some of the stars in the stream reported by Helmi et al. (1999). By modelling, we estimate that both structures had progenitors with virial masses $\approx 10^{10} M_\odot$ and infall times $\gtrsim 9$ Gyr ago. Using abundance matching, these correspond to stellar masses between $10^6$ and $10^7 M_\odot$. These are somewhat larger than the masses inferred through the mass-metallicity relation by factors of 5 to 15. Additionally, we identify two further substructures (S3 and S4 with 55 and 40 members) and two clusters or moving groups (C1 and C2 with 24 and 12) members. In all 6 cases, clustering in kinematics is found to correspond to clustering in both configuration space and metallicity, adding credence to the reliability of our detections.


title: Phonon dispersion and the competition between pairing and charge order

abstract: The Holstein Model (HM) describes the interaction between fermions and a collection of local (dispersionless) phonon modes. In the dilute limit, the phonon degrees of freedom dress the fermions, giving rise to polaron and bipolaron formation. At higher densities, the phonons mediate collective superconducting (SC) and charge density wave (CDW) phases. Quantum Monte Carlo (QMC) simulations have considered both these limits, but have not yet focused on the physics of more general phonon spectra. Here we report QMC studies of the role of phonon dispersion on SC and CDW order in such models. We quantify the effect of finite phonon bandwidth and curvature on the critical temperature $T_{\rm cdw}$ for CDW order, and also uncover several novel features of diagonal long range order in the phase diagram, including a competition between charge patterns at momenta ${\bf q}=(\pi,\pi)$ and ${\bf q}=(0,\pi)$ which lends insight into the relationship between Fermi surface nesting and the wavevector at which charge order occurs. We also demonstrate SC order at half-filling in situations where nonzero bandwidth sufficiently suppresses $T_{\rm cdw}$.


title: On a renormalizable class of gauge fixings for the gauge invariant operator $A_{\min }^{2}$

abstract: The dimension two gauge invariant non-local operator $A_{\min }^{2}$, obtained through the minimization of $\int d^4x A^2$ along the gauge orbit, allows to introduce a non-local gauge invariant configuration $A^h_\mu$ which can be employed to built up a class of Euclidean massive Yang-Mills models useful to investigate non-perturbative infrared effects of confining theories. A fully local setup for both $A_{\min }^{2}$ and $A^{h}_\mu$ can be achieved, resulting in a local and BRST invariant action which shares similarities with the Stueckelberg formalism. Though, unlike the case of the Stueckelberg action, the use of $A_{\min }^{2}$ gives rise to an all orders renormalizable action, a feature which will be illustrated by means of a class of covariant gauge fixings which, as much as 't Hooft's $R_\zeta$-gauge of spontaneously broken gauge theories, provide a mass for the Stueckelberg field.


title: Predicting the Spin Seebeck Voltage in Spin-polarized Materials: A Quantum Mechanical Transport Model Approach

abstract: The spin Seebeck effect has recently been demonstrated as a viable method of direct energy conversion that has potential to outperform energy conversion from the conventional Seebeck effect. In this study, a computational transport model is developed and validated that predicts the spin Seebeck voltage in spin-polarized materials using material parameter obtain from first principle ground state density functional calculations. The transport model developed is based on a 1D effective mass description coupled with a microscopic inverse spin Hall relationship. The model can predict both the spin current and voltage generated in a non-magnetic material placed on top of a ferromagnetic material in a transverse spin Seebeck configuration. The model is validated and verified with available experimental data of La:YIG. Future applications of this model include the high-throughput exploration of new spin-based thermoelectric materials.


title: Diagnosing added value of convection-permitting regional models using precipitation event identification and tracking

abstract: Dynamical downscaling with high-resolution regional climate models may offer the possibility of realistically reproducing precipitation and weather events in climate simulations. As resolutions fall to order kilometers, the use of explicit rather than parametrized convection may offer even greater fidelity. However, these increased model resolutions both allow and require increasingly complex diagnostics for evaluating model fidelity. In this study we use a suite of dynamically downscaled simulations of the summertime U.S. (WRF driven by NCEP) with systematic variations in parameters and treatment of convection as a test case for evaluation of model precipitation. In particular, we use a novel rainstorm identification and tracking algorithm that allocates essentially all rainfall to individual precipitation events (Chang et al. 2016). This approach allows multiple insights, including that, at least in these runs, model wet bias is driven by excessive areal extent of precipitating events. Biases are time-dependent, producing excessive diurnal cycle amplitude. We show that this effect is produced not by new production of events but by excessive enlargement of long-lived precipitation events during daytime, and that in the domain average, precipitation biases appear best represented as additive offsets. Of all model configurations evaluated, convection-permitting simulations most consistently reduced biases in precipitation event characteristics.


title: In a Nutshell -- The Sequential Parameter Optimization Toolbox

abstract: The performance of optimization algorithms relies crucially on their parameterizations. Finding good parameter settings is called algorithm tuning. The sequential parameter optimization (SPOT) package for R is a toolbox for tuning and understanding simulation and optimization algorithms. Model-based investigations are common approaches in simulation and optimization. Sequential parameter optimization has been developed, because there is a strong need for sound statistical analysis of simulation and optimization algorithms. SPOT includes methods for tuning based on classical regression and analysis of variance techniques; tree-based models such as CART and random forest; Gaussian process models (Kriging), and combinations of different meta-modeling approaches. Using a simple simulated annealing algorithm, we will demonstrate how optimization algorithms can be tuned using SPOT. The underling concepts of the SPOT approach are explained. This includes key techniques such as exploratory fitness landscape analysis and sensititvity analysis. Many examples illustrate how SPOT can be used for understanding the performance of algorithms and gaining insight into algorithm's behavior. Furthermore, we demonstrate how SPOT can be used as an optimizer and how a sophisticated ensemble approach is able to combine several meta models via stacking. This article exemplifies how SPOT can be used for automatic and interactive tuning.


title: Observations and light curve solutions of ultrashort-period eclipsing binaries

abstract: Photometric observations in V and I bands and low-dispersion spectra of ten ultrashort-period binaries (NSVS 2175434, NSVS 2607629, NSVS 5038135, NSVS 8040227, NSVS 9747584, NSVS 4876238, ASAS 071829-0336.7, SWASP 074658.62+224448.5, NSVS 2729229, NSVS 10632802) are presented. One of them, NSVS 2729229, is newly discovered target. The results from modeling and analysis of our observations revealed that: (i) Eight targets have overcontact configurations with considerable fillout factor (up to 0.5) while NSVS 4876238 and ASAS 0718-03 have almost contact configurations; (ii) NSVS 4876238 is rare ultrashort-period binary of detached type; (iii) all stellar components are late dwarfs; (iv) the temperature difference of the components of each target does not exceed 400 K; (v) NSVS 2175434 and SWASP 074658.62+224448.5 exhibit total eclipses and their parameters could be assumed as well-determined; (v) NSVS 2729229 shows emission in the H_{alpha} line. Masses, radii and luminosities of the stellar components were estimated by the empirical relation "period, orbital axis" for short- and ultrashort-period binaries. We found linear relations mass-luminosity and mass-radius for the stellar components of our targets.


title: Guidelines for Producing Useful Synthetic Data

abstract: We report on our experiences of helping staff of the Scottish Longitudinal Study to create synthetic extracts that can be released to users. In particular, we focus on how the synthesis process can be tailored to produce synthetic extracts that will provide users with similar results to those that would be obtained from the original data. We make recommendations for synthesis methods and illustrate how the staff creating synthetic extracts can evaluate their utility at the time they are being produced. We discuss measures of utility for synthetic data and show that one tabular utility measure is exactly equivalent to a measure calculated from a propensity score. The methods are illustrated by using the R package $synthpop$ to create synthetic versions of data from the 1901 Census of Scotland.


title: Interference Effects in a Very Precise Measurement of the Muon Charge Asymmetry at FCC-ee

abstract: At the future high luminosity electron-positron collider FCC-ee proposed for CERN, the precise measurement of the charge asymmetry in $e^+ e^- \to \mu^+ \mu^-$ near the Z resonance is of special interest. In particular, such a measurement at $M_Z \pm 3.5$ GeV may provide a very precise measurement of the electromagnetic coupling at the scale $\sim M_Z$, a fundamental constant of the Standard Model. However, this charge asymmetry is plagued by a large trivial contribution from the interference of photon emission from initial state electrons and final state muons. We address the question whether this interference can be reliably calculated and subtracted with the help of a resummed QED calculation.


title: Convexity in ordered matroids and the generalized external order

abstract: In 1980, Las Vergnas defined a notion of discrete convexity for oriented matroids, which Edelman subsequently related to the theory of anti-exchange closure functions and convex geometries. In this paper, we use generalized matroid activity to construct a convex geometry associated with an ordered, unoriented matroid. The construction in particular yields a new type of representability for an ordered matroid defined by the affine representability of its corresponding convex geometry. The lattice of convex sets of this convex geometry induces an ordering on the matroid independent sets which extends the external active order on matroid bases. We show that this generalized external order forms a supersolvable meet-distributive lattice refining the geometric lattice of flats, and we uniquely characterize the lattices isomorphic to the external order of a matroid. Finally, we introduce a new trivariate generating function generalizing the matroid Tutte polynomial.


title: Hypergraphs not containing a tight tree with a bounded trunk

abstract: An $r$-uniform hypergraph is a tight $r$-tree if its edges can be ordered so that every edge $e$ contains a vertex $v$ that does not belong to any preceding edge and the set $e-v$ lies in some preceding edge. A conjecture of Kalai [Kalai], generalizing the Erd\H{o}s-S\'os Conjecture for trees, asserts that if $T$ is a tight $r$-tree with $t$ edges and $G$ is an $n$-vertex $r$-uniform hypergraph containing no copy of $T$ then $G$ has at most $\frac{t-1}{r}\binom{n}{r-1}$ edges. A trunk $T'$ of a tight $r$-tree $T$ is a tight subtree such that every edge of $T-T'$ has $r-1$ vertices in some edge of $T'$ and a vertex outside $T'$. For $r\ge 3$, the only nontrivial family of tight $r$-trees for which this conjecture has been proved is the family of $r$-trees with trunk size one in [FF] from 1987. Our main result is an asymptotic version of Kalai's conjecture for all tight trees $T$ of bounded trunk size. This follows from our upper bound on the size of a $T$-free $r$-uniform hypergraph $G$ in terms of the size of its shadow. We also give a short proof of Kalai's conjecture for tight $r$-trees with at most four edges. In particular, for $3$-uniform hypergraphs, our result on the tight path of length $4$ implies the intersection shadow theorem of Katona [Katona].


title: Engineering the eigenstates of coupled spin-1/2 atoms on a surface

abstract: Quantum spin networks having engineered geometries and interactions are eagerly pursued for quantum simulation and access to emergent quantum phenomena such as spin liquids. Spin-1/2 centers are particularly desirable because they readily manifest coherent quantum fluctuations. Here we introduce a controllable spin-1/2 architecture consisting of titanium atoms on a magnesium oxide surface. We tailor the spin interactions by atomic-precision positioning using a scanning tunneling microscope (STM), and subsequently perform electron spin resonance (ESR) on individual atoms to drive transitions into and out of quantum eigenstates of the coupled-spin system. Interactions between the atoms are mapped over a range of distances extending from highly anisotropic dipole coupling, to strong exchange coupling. The local magnetic field of the magnetic STM tip serves to precisely tune the superposition states of a pair of spins. The precise control of the spin-spin interactions and ability to probe the states of the coupled-spin network by addressing individual spins will enable exploration of quantum many-body systems based on networks of spin-1/2 atoms on surfaces.


title: Learning Compressible 360{\deg} Video Isomers

abstract: Standard video encoders developed for conventional narrow field-of-view video are widely applied to 360{\deg} video as well, with reasonable results. However, while this approach commits arbitrarily to a projection of the spherical frames, we observe that some orientations of a 360{\deg} video, once projected, are more compressible than others. We introduce an approach to predict the sphere rotation that will yield the maximal compression rate. Given video clips in their original encoding, a convolutional neural network learns the association between a clip's visual content and its compressibility at different rotations of a cubemap projection. Given a novel video, our learning-based approach efficiently infers the most compressible direction in one shot, without repeated rendering and compression of the source video. We validate our idea on thousands of video clips and multiple popular video codecs. The results show that this untapped dimension of 360{\deg} compression has substantial potential--"good" rotations are typically 8-10% more compressible than bad ones, and our learning approach can predict them reliably 82% of the time.


title: Fourier spectral approximation for the convective Cahn-Hilliard equation in 2D cas

abstract: In this paper, we consider the Fourier spectral method for numerically solving the 2D convective Cahn-Hilliard equation. The semi-discrete and fully discrete schemes are established. Moreover, the existence, uniqueness and the optimal error bound are also considered.


title: Geometric properties of minimizers in the planar three-body problem with two equal masses

abstract: It it shown that each lobe of the figure-eight orbit is star-shaped, which implies the polar angle is monotone in each lobe. In general, it is not clear when a minimizer is star-shaped. In this paper, we study minimizers connecting two fixed-ends (i.e. the Bolza problem) in the planar three-body problem with two equal masses. We show that if the Jacobi coordinates of the two fixed-ends are in adjacent closed quadrants, then the corresponding minimizer must stay in two adjacent closed quadrants. If we further assume the two Jacobi coordinates are orthogonal on one of the fixed-ends, then the polar angles of the Jacobi coordinates in the minimizer have at most one critical point. If the two Jacobi coordinates are orthogonal on both ends, then the two polar angles must be monotone. These geometric properties can be applied to show the existence of two sets of periodic orbits.


title: PacGAN: The power of two samples in generative adversarial networks

abstract: Generative adversarial networks (GANs) are innovative techniques for learning generative models of complex data distributions from samples. Despite remarkable recent improvements in generating realistic images, one of their major shortcomings is the fact that in practice, they tend to produce samples with little diversity, even when trained on diverse datasets. This phenomenon, known as mode collapse, has been the main focus of several recent advances in GANs. Yet there is little understanding of why mode collapse happens and why existing approaches are able to mitigate mode collapse. We propose a principled approach to handling mode collapse, which we call packing. The main idea is to modify the discriminator to make decisions based on multiple samples from the same class, either real or artificially generated. We borrow analysis tools from binary hypothesis testing---in particular the seminal result of Blackwell [Bla53]---to prove a fundamental connection between packing and mode collapse. We show that packing naturally penalizes generators with mode collapse, thereby favoring generator distributions with less mode collapse during the training process. Numerical experiments on benchmark datasets suggests that packing provides significant improvements in practice as well.


title: Local sensitivity analysis for the Cucker-Smale model with random inputs

abstract: We present pathwise flocking dynamics and local sensitivity analysis for the Cucker-Smale(C-S) model with random communications and initial data. For the deterministic communications, it is well known that the C-S model can model emergent local and global flocking dynamics depending on initial data and integrability of communication function. However, the communication mechanism between agents are not a priori clear and needs to be figured out from observed phenomena and data. Thus, uncertainty in communication is an intrinsic component in the flocking modeling of the C-S model. In this paper, we provide a class of admissible random uncertainties which allows us to perform the local sensitivity analysis for flocking and establish stability to the random C-S model with uncertain communication.


title: Zero-Modified Poisson-Lindley distribution with applications in zero-inflated and zero-deflated count data

abstract: The main object of this article is to present an extension of the zero-inflated Poisson-Lindley distribution, called of zero-modified Poisson-Lindley. The additional parameter $\pi$ of the zero-modified Poisson-Lindley has a natural interpretation in terms of either zero-deflated/inflated proportion. Inference is dealt with by using the likelihood approach. In particular the maximum likelihood estimators of the distribution's parameter are compared in small and large samples. We also consider an alternative bias-correction mechanism based on Efron's bootstrap resampling. The model is applied to real data sets and found to perform better than other competing models.


title: Regularity of Kleinian limit sets and Patterson-Sullivan measures

abstract: We consider several (related) notions of geometric regularity in the context of limit sets of geometrically finite Kleinian groups and associated Patterson-Sullivan measures. We begin by computing the upper and lower regularity dimensions of the Patterson-Sullivan measure, which involves controlling the relative measure of concentric balls. We then compute the Assouad and lower dimensions of the limit set, which involves controlling local doubling properties. Unlike the Hausdorff, packing, and box-counting dimensions, we show that the Assouad and lower dimensions are not necessarily given by the Poincar\'e exponent.


title: Shape (In)dependent Inequalities for Triangleland's Jacobi and Democratic-Linear Ellipticity Quantitities

abstract: Sides and medians are both Jacobi coordinate magnitudes, moreover then equably entering the spherical coordinates on Kendall's shape sphere and the Hopf coordinates. This motivates treating medians on the same footing as sides in triangle geometry and the resulting Shape Theory. In this paper, we consequently reformulate inequalities for the medians in terms of shape quantities, and proceed to find inequalities on the mass-weighted Jacobi coordinates. This work moreover identifies the $4/3$ -- powers of which occur frequently in the theory of medians -- as the ratio of Jacobi masses. One of the Hopf coordinates is tetra-area. Another is anisoscelesness, which parametrizes whether triangles are left-or-right leaning as bounded by isoscelesness itself. The third is ellipticity, which parametrizes tallness-or-flatness of triangles as bounded by regular triangles. Whereas tetra-area is clearly cluster choice invariant, Jacobi coordinates, anisoscelesness and ellipticity are cluster choice dependent but can be `democratized' by averaging over all clusters. Democratized ellipticity moreover trivializes, due to ellipticity being the difference of base-side and median second moments, whose averages are equal to each other. Thus we introduce a distinct `linear ellipticity' quantifier of tallness-or-flatness of triangles whose democratization is nontrivial, and find inequalities bounding this. Some of this paper's inequalities are shape-independent bounds, whereas others' bounds depend on the isoperimetric ratio and arithmetic-to-geometric side mean ratio shape variables.


title: On ancient solutions of the heat equation

abstract: An explicit representation formula for all positive ancient solutions of the heat equation in the Euclidean case is found. In the Riemannian case with nonnegative Ricci curvature, a similar but less explicit formula is also found. Here it is proven that any positive ancient solution is the standard Laplace transform of positive solutions of the family of elliptic operator $\Delta - s$ with $s>0$. Further relaxation of the curvature assumption is also possible. It is also shown that the linear space of ancient solutions of polynomial growth has finite dimension and these solutions are polynomials in time.


title: Null Geodesics and Wave Front Singularities in the Godel Space-time

abstract: We explore wave fronts of null geodesics in the Godel metric emitted from point sources both at, and away from, the origin. For constant time wave fronts emitted by sources away from the origin, we find cusp ridges as well as blue sky metamorphoses where spatially disconnected portions of the wave front appear, connect to the main wave front, and then later break free and vanish. These blue sky metamorphoses in the constant time wave fronts highlight the non-causal features of the Godel metric. We introduce a concept of physical distance along the null geodesics, and show that for wave fronts of constant physical distance, the reorganization of the points making up the wave front leads to the removal of cusp ridges.


title: Plasma Medicine and Osmosis

abstract: In this paper, attention is drawn to the importance of accounting for osmotic pressure when analyzing physiological effects on cellular structures in plasma medicine. Interaction of a weakly ionized plasma jet with a saline solution leads to a detectable changes in the saline's ion-molecular composition and hence changes in the osmotic pressure. This, in turn, leads to a stretching or compression of the membrane, depending on the difference of total external and internal pressures. The selective effect of plasma on cells, observed in experiments, is associated with the change in the mechanical properties of membranes (and thereby, a weakening of their protective properties). Corresponding estimates are given in the article.


title: Search for excited $B_c^+$ states

abstract: A search is performed in the invariant mass spectrum of the $B_c^{+}\pi^{+}\pi^{-}$ system for the excited $B_c^{+}$ states $B_c(2^{1}S_{0})^+$ and $B_c(2^{3}S_{1})^+$ using a data sample of $pp$ collisions collected by the LHCb experiment at the centre-of-mass energy of $\sqrt{s} = 8 \,{\mathrm{TeV}}$, corresponding to an integrated luminosity of $2 \,{\mathrm{fb^{-1}}}$. No evidence is seen for either state. Upper limits on the ratios of the production cross-sections of the $B_c(2^{1}S_{0})^+$ and $B_c(2^{3}S_{1})^+$ states times the branching fractions of ${B_c(2^{1}S_{0})^+} \to {B_c^{+}\pi^{+}\pi^{-}}$ and ${B_c(2^{3}S_{1})^+} \to {B_c^{*+}\pi^{+}\pi^{-}}$ over the production cross-section of the $B_c^{+}$ state are given as a function of their masses. They are found to be between 0.02 and 0.14 at $95\%$ confidence level for $B_c(2^{1}S_{0})^+$ and $B_c(2^{3}S_{1})^+$ in the mass ranges $[6830, 6890] \,{\mathrm{MeV}}/c^{2}$ and $[6795,6890] \,{\mathrm{MeV}}/c^{2}$, respectively.


title: Systems of Points with Coulomb Interactions

abstract: Large ensembles of points with Coulomb interactions arise in various settings of condensed matter physics, classical and quantum mechanics, statistical mechanics, random matrices and even approximation theory, and give rise to a variety of questions pertaining to calculus of variations, Partial Differential Equations and probability. We will review these as well as "the mean-field limit" results that allow to derive effective models and equations describing the system at the macroscopic scale. We then explain how to analyze the next order beyond the mean-field limit, giving information on the system at the microscopic level. In the setting of statistical mechanics, this allows for instance to observe the effect of the temperature and to connect with crystallization questions.


title: Deuterium Chemistry in the Young Massive Protostellar Core NGC 2264 CMM3

abstract: In this work we present the first attempt of modelling the deuterium chemistry in the massive young protostellar core NGC 2264 CMM3. We investigated the sensitivity of this chemistry to the physical conditions in its surrounding environment. The results showed that deuteration, in the protostellar gas, is affected by variations in the core density, the amount of gas depletion onto grain surfaces, the CR ionisation rate, but it is insensitive to variations in the H$_2$ ortho-to-para ratio. Our results, also, showed that deuteration is often enhanced in less-dense, partially depleted ($<$ 85\%), or cores that are exerted to high CR ionisation rates ($\ge$ 6.5 $\times$ 10$^{-17}$ s$^{-1}$ ). However, in NGC 2264 CMM3, decreasing the amount of gas depleted onto grains and enhancing the CR ionisation rate are often overestimating the observed values in the core. The best fit time to observations occurs around (1 - 5) $\times$ 10$^4$ yrs for core densities in the range (1 - 5) $\times$ 10$^6$ cm$^{-3}$ with CR ionisation rate between (1.7 - 6.5)$\times$ 10$^{-17}$ s$^{-1}$. These values are in agreement with the results of the most recent theoretical chemical model of CMM3, and the time range of best fit is, also, in-line with the estimated age of young protostellar objects. We conclude that deuterium chemistry in protostellar cores is: (i) sensitive to variations in the physical conditions in its environment, %(ii) deuteration in the protostellar gas is (ii) insensitive to changes in the H$_2$ ortho-to-para ratio. %, but is sensitive to these variations in the dark cloud phase. We also conclude that the core NGC 2264 CMM3 is in its early stages of chemical evolution with an estimated age of (1 - 5) $\times$ 10$^4$ yrs.


title: All-optical wavelength-tunable narrow-linewidth fiber laser

abstract: Parameter regulations of narrow-linewidth fiber lasers in frequency domain has drawn considerable interests for widespread applications in the light quantum computing, precise coherent detection, and generation of micro-waves. All-optical methods provide compact, precise and fast accesses to achieving these lasers with wavelength-tunability. Here, the optical-thermal effects of graphene is utilized to precisely control operations of free-running lasers with a tuning speed of 140 MHz/ms. Assisted by the single-longitude-mode operation and linewidth suppression of stimulated Brillouin backscattering, we obtain an optical-controllable ~750 Hz fiber laser with a wavelength-tuning range of 3.7 nm.


title: Normal Convergence Using Malliavin Calculus With Applications and Examples

abstract: We prove the chain rule in the more general framework of the Wiener-Poisson space, allowing us to obtain the so-called Nourdin-Peccati bound. From this bound we obtain a second-order Poincare-type inequality that is useful in terms of computations. For completeness we survey these results on the Wiener space, the Poisson space, and the Wiener-Poisson space. We also give several applications to central limit theorems with relevant examples: linear functionals of Gaussian subordinated fields (where the subordinated field can be processes like fractional Brownian motion or the solution of the Ornstein-Uhlenbeck SDE driven by fractional Brownian motion), Poisson functionals in the first Poisson chaos restricted to infinitely many \small" jumps (particularly fractional Levy processes) and the product of two Ornstein-Uhlenbeck processes (one in the Wiener space and the other in the Poisson space). We also obtain bounds for their rate of convergence to normality.


title: Towards a proof of the 24-cell conjecture

abstract: This review paper is devoted to the problems of sphere packings in 4 dimensions. The main goal is to find reasonable approaches for solutions to problems related to densest sphere packings in 4-dimensional Euclidean space. We consider two long-standing open problems: the uniqueness of maximum kissing arrangements in 4 dimensions and the 24-cell conjecture. Note that a proof of the 24-cell conjecture also proves that the checkerboard lattice packing D4 is the densest sphere packing in 4 dimensions.


title: The Investigation of Social Media Data Thresholds for Opinion Formation

abstract: The pervasive use of social media has grown to over two billion users to date, and is commonly utilized as a means to share information and shape world events. Evidence suggests that passive social media usage (i.e., viewing without taking action) has an impact on the user's perspective. This empirical influence over perspective could have significant impact on social events. Therefore, it is important to understand how social media contributes to the formation of an individual's perspective. A set of experimental tasks were designed to investigate empirically derived thresholds for opinion formation as a result of passive interactions with different social media data types (i.e., videos, images, and messages). With a better understanding of how humans passively interact with social media information, a paradigm can be developed that allows the exploitation of this interaction and plays a significant role in future military plans and operations.


title: Deep Reinforcement Learning Boosted by External Knowledge

abstract: Recent improvements in deep reinforcement learning have allowed to solve problems in many 2D domains such as Atari games. However, in complex 3D environments, numerous learning episodes are required which may be too time consuming or even impossible especially in real-world scenarios. We present a new architecture to combine external knowledge and deep reinforcement learning using only visual input. A key concept of our system is augmenting image input by adding environment feature information and combining two sources of decision. We evaluate the performances of our method in a 3D partially-observable environment from the Microsoft Malmo platform. Experimental evaluation exhibits higher performance and faster learning compared to a single reinforcement learning model.


title: 3D Characterisation of the Fe-rich intermetallic phases in Al-5%Cu alloys by synchrotron X-ray microtomography and skeletonisation

abstract: Synchrotron X-ray microtomography and skeletonisation method were used to study the true 3D network structures and morphologies of the Fe-rich intermetallic phases in Al-5.0%Cu-0.6%Mn alloys with 0.5% and 1.0% Fe. It was found that, the Fe-phases in the 1.0%Fe alloy have node lengths of 5-25m; while those in the 0.5%Fe alloy are of 3-17 m. The Fe-phases in the 1.0%Fe alloy also developed sharper mean curvature with wider distribution than those in the 0.5%Fe alloy. Combining SEM studies of the deeply-etched samples, the true 3D structures of 4 different type Fe-phases in both alloys are also revealed and demonstrated.


title: Search for the Signatures of a New-Born Black Hole from the Collapse of a Supra-massive Millisecond Magnetar in Short GRB Light Curves

abstract: `Internal plateau' followed by a sharp decay is commonly seen in short gamma-ray burst (GRB) light curves. The plateau component is usually interpreted as the dipole emission from a supra-massive magnetar, and the sharp decay may imply the collapse of the magnetar to a black hole (BH). Fall-back accretion onto the new-born BH could produce long-lasting activities via the Blandford-Znajek (BZ) process. The magnetic flux accumulated near the BH would be confined by the accretion disks for a period of time. As the accretion rate decreases, the magnetic flux is strong enough to obstruct gas infall, leading to a magnetically-arrested disk (MAD). Within this scenario, we show that the BZ process could produce two types of typical X-ray light curves: type I exhibits a long-lasting plateau, followed by a power-law decay with slopes ranging from 5/3 to 40/9; type II shows roughly a single power-law decay with slope of 5/3. The former requires low magnetic filed strength, while the latter corresponds to relatively high values. We search for such signatures of the new-born BH from a sample of short GRBs with an internal plateau, and find two candidates: GRB 101219A and GRB 160821B, corresponding to type II and type I light curve, respectively. It is shown that our model can explain the data very well.


title: Convergence Rates for Deterministic and Stochastic Subgradient Methods Without Lipschitz Continuity

abstract: We extend the classic convergence rate theory for subgradient methods to apply to non-Lipschitz functions. For the deterministic projected subgradient method, we present a global $O(1/\sqrt{T})$ convergence rate for any convex function which is locally Lipschitz around its minimizers. This approach is based on Shor's classic subgradient analysis and implies generalizations of the standard convergence rates for gradient descent on functions with Lipschitz or H\"older continuous gradients. Further, we show a $O(1/\sqrt{T})$ convergence rate for the stochastic projected subgradient method on convex functions with at most quadratic growth, which improves to $O(1/T)$ under either strong convexity or a weaker quadratic lower bound condition.


title: Dynamically Reconfigurable Sources for Arbitrary Gaussian States in Integrated Photonics Circuits

abstract: We present a modular design for integrated programmable multimode sources of arbitrary Gaussian states of light. The technique is based on current technologies, in particular recent demonstrations of on-chip photon manipulation and generation of highly squeezed vacuum states in semiconductors. While the design is generic and independent of the choice of integrated platform, we adopt recent experimental results on compound semiconductors as a demonstrative example. Such a device would be valuable as a source for many quantum protocols that range from imaging to communication and information processing.


title: Sparse Phase Retrieval via Sparse PCA Despite Model Misspecification: A Simplified and Extended Analysis

abstract: We consider the problem of high-dimensional misspecified phase retrieval. This is where we have an $s$-sparse signal vector $\mathbf{x}_*$ in $\mathbb{R}^n$, which we wish to recover using sampling vectors $\textbf{a}_1,\ldots,\textbf{a}_m$, and measurements $y_1,\ldots,y_m$, which are related by the equation $f(\left<\textbf{a}_i,\textbf{x}_*\right>) = y_i$. Here, $f$ is an unknown link function satisfying a positive correlation with the quadratic function. This problem was analyzed in a recent paper by Neykov, Wang and Liu, who provided recovery guarantees for a two-stage algorithm with sample complexity $m = O(s^2\log n)$. In this paper, we show that the first stage of their algorithm suffices for signal recovery with the same sample complexity, and extend the analysis to non-Gaussian measurements. Furthermore, we show how the algorithm can be generalized to recover a signal vector $\textbf{x}_*$ efficiently given geometric prior information other than sparsity.


title: Attack Vulnerability of Complex Networks in Center-Based Strategies

abstract: Central nodes are critical in establishing structural connectivity in a complex network. Attacking such nodes can create real havoc in a complex system. We propose attack strategies based on four types of centers, namely betweenness center, degree center, median, and center. We study the vulnerability of synthetic as well as real-world networks in these center-based attacks. These attacks are node-removal attacks which involve identifying the central node set and removing them from the network. We observed that the attacks based on recalculated network information are more efficient than ones based on initial network information. This work shows that the median-based attack, which is a novel strategy proposed in this work, is highly destructive in real as well as synthetic networks.


title: Incremental View Maintenance for Property Graph Queries

abstract: This paper discusses the challenges of incremental view maintenance for property graph queries. We select a subset of property graph queries and present an approach that uses nested relational algebra to allow incremental evaluation.


title: Im2Flow: Motion Hallucination from Static Images for Action Recognition

abstract: Existing methods to recognize actions in static images take the images at their face value, learning the appearances---objects, scenes, and body poses---that distinguish each action class. However, such models are deprived of the rich dynamic structure and motions that also define human activity. We propose an approach that hallucinates the unobserved future motion implied by a single snapshot to help static-image action recognition. The key idea is to learn a prior over short-term dynamics from thousands of unlabeled videos, infer the anticipated optical flow on novel static images, and then train discriminative models that exploit both streams of information. Our main contributions are twofold. First, we devise an encoder-decoder convolutional neural network and a novel optical flow encoding that can translate a static image into an accurate flow map. Second, we show the power of hallucinated flow for recognition, successfully transferring the learned motion into a standard two-stream network for activity recognition. On seven datasets, we demonstrate the power of the approach. It not only achieves state-of-the-art accuracy for dense optical flow prediction, but also consistently enhances recognition of actions and dynamic scenes.


title: Observing galaxy mergers at the epoch of reionization

abstract: The galaxies with photometric redshifts observed in a close angular proximity might be either projection coincidences, strongly lensed images of the same galaxy, or separate galaxies that are in a stage of merging. We search for the groups of galaxies in the Hubble Ultra Deep Field (HUDF09) in $z\sim7$ and $z\sim8$ drop-out samples. We find no close pairs among 50 galaxies in the $z\sim7$ sample, while in the $z\sim8$ sample we find that 6 out of 22 galaxies have a companion within $\sim1'$ (3 pairs). Adopting a numerical simulation and performing forward modeling we show that even though mergers are unlikely to have such a high fraction, the projection coincidences and the strong lensing are even less likely mechanisms to account for all of three pairs. Alternatively, there is a possibility of the contamination in the drop-out catalog from lower redshifts, which potentially can account for all of the groups. Finally, we make projection on the sensitivity to mergers of the James Webb Space Telescope, and discuss the possible applications of the high-redshift merging galaxies for decreasing cosmic variance effect on the luminosity function and for improving the accuracy of photometric redshifts in general.


title: Tunable Graphene Metasurface Reflectarray for Cloaking, Illusion and Focusing

abstract: We present a graphene-based metasurface that can be actively tuned between different regimes of operation, such as anomalous beam steering and focusing, cloaking and illusion optics, by applying electrostatic gating without modifying the geometry of the metasurface. The metasurface is designed by placing graphene nano-ribbons (GNRs) on a dielectric cavity resonator, where interplay between geometric plasmon resonances in the ribbons and Fabry-Perot resonances in the cavity is used to achieve 2$\pi$ phase shift. As a proof of the concept, we demonstrate that wavefront of the field reflected from a triangular bump covered by the metasurface can be tuned by applying electric bias so as to resemble that of bare plane and of a spherical object. Moreover, reflective focusing and change of the reflection direction for the above-mentioned cases are also shown.


title: Theoretical Model and Characteristics of Mitochondrial Thermogenesis

abstract: Mitochondria of brown adipocyte (BA) are the main intracellular sites for thermogenesis, which have been targeted for therapy to reduce obesity. However, there are long-standing critique and debates about the ability of raising cellular temperature by endogenous thermogenesis. Currently, the wrong theoretical model gave about the five orders of magnitude less than facts. Here, based on the first law of thermodynamics and thermal diffusion equation, the deduced theoretical model of mitochondrial thermogenesis satisfies Laplace equation, and is a special case for thermal diffusion equation. The model settles the long-standing questioning about the ability of raising cellular temperature by endogenous thermogenesis, and explains the thermogenic characteristics of brown adipocyte. The model and calculations also suggest that the number of free available proton is the major limiting factor for endogenous thermogenesis and its speed.


title: Direct observation of heterogeneous valence state in Yb-based quasicrystalline approximants

abstract: We study the electronic structure of Tsai-type cluster-based quasicrystalline approximants, Au$_{64}$Ge$_{22}$Yb$_{14}$ (AGY-I), Au$_{63.5}$Ge$_{20.5}$Yb$_{16}$ (AGY-II), and Zn$_{85.4}$Yb$_{14.6}$ (Zn-Yb), by means of photoemission spectroscopy. In the valence band hard x-ray photoemission spectra of AGY-II and Zn-Yb, we separately observe a fully occupied Yb 4$f$ state and a valence fluctuation derived Kondo resonance peak, reflecting two inequivalent Yb sites, a single Yb atom in the cluster center and its surrounding Yb icosahedron, respectively. The fully occupied 4$f$ signal is absent in AGY-I containing no Yb atom in the cluster center. The results provide direct evidence for a heterogeneous valence state in AGY-II and Zn-Yb.


title: ESASky v.2.0: all the skies in your browser

abstract: With the goal of simplifying the access to science data to scientists and citizens, ESA recently released ESASky (http://sky.esa.int), a new open-science easy-to-use portal with the science-ready Astronomy data from ESA and other major data providers. In this presentation, we announced version 2.0 of the application, which includes access to all science-ready images, catalogues and spectra, a feature to help planning of future JWST observations, the possibility to search for data of all (targeted and serendipitously observed) Solar System Objects in Astronomy images, a first support to mobile devices and several other smaller usability features. We also discussed the future evolution of the portal and the lessons learnt from the 1+ year of operations from the point of view of access, visualization and manipulation of big datasets (all sky maps, also called HiPS) and large catalogues (like e.g. the Gaia DR1 catalogues or the Hubble Source Catalogue) and the design and validation principles for the development of friendly GUIs for thin layer web clients aimed at scientists.


title: Electronic structure, charge transfer, and intrinsic luminescence of gadolinium oxide nanoparticles: Experiment and theory

abstract: The cubic (c) and monoclinic (m) polymorphs of Gd2O3 were studied using the combined analysis of several materials science techniques - X-ray diffraction (XRD), scanning electron microscopy (SEM), X-ray photoelectron spectroscopy (XPS), and photoluminescence (PL) spectroscopy. Density functional theory (DFT) based calculations for the samples under study were performed as well. The cubic phase of gadolinium oxide (c-Gd2O3) synthesized using a precipitation method exhibits spheroidal-like nanoclusters with well-defined edges assembled from primary nanoparticles with an average size of 50 nm, whereas the monoclinic phase of gadolinium oxide (m-Gd2O3) deposited using explosive pyrolysis has a denser structure compared with natural gadolinia. This phase also has a structure composed of three-dimensional complex agglomerates without clear-edged boundaries that are ~21 nm in size plus a cubic phase admixture of only 2 at. % composed of primary edge-boundary nanoparticles ~15 nm in size. These atomic features appear in the electronic structure as different defects ([Gd...O-OH] and [Gd...O-O]) and have dissimilar contributions to the charge-transfer processes among the appropriate electronic states with ambiguous contributions in the Gd 5p - O 2s core-like levels in the valence band structures. The origin of [Gd...O-OH] defects found by XPS was well-supported by PL analysis. The electronic and atomic structures of the synthesized gadolinias calculated using DFT were compared and discussed on the basis of the well-known joint OKT-van der Laan model, and good agreement was established.


title: A Novel Document Generation Process for Topic Detection based on Hierarchical Latent Tree Models

abstract: We propose a novel document generation process based on hierarchical latent tree models (HLTMs) learned from data. An HLTM has a layer of observed word variables at the bottom and multiple layers of latent variables on top. For each document, we first sample values for the latent variables layer by layer via logic sampling, then draw relative frequencies for the words conditioned on the values of the latent variables, and finally generate words for the document using the relative word frequencies. The motivation for the work is to take word counts into consideration with HLTMs. In comparison with LDA-based hierarchical document generation processes, the new process achieves drastically better model fit with much fewer parameters. It also yields more meaningful topics and topic hierarchies. It is the new state-of-the-art for the hierarchical topic detection.


title: The Calculus of Democratization and Development

abstract: In accordance with "Democracy's Effect on Development: More Questions than Answers", we seek to carry out a study in following the description in the 'Questions for Further Study.' To that end, we studied 33 countries in the Sub-Saharan Africa region, who all went through an election which should signal a "step-up" for their democracy, one in which previously homogenous regimes transfer power to an opposition party that fairly won the election. After doing so, liberal-democracy indicators and democracy indicators were evaluated in the five years prior to and after the election took place, and over that ten-year period, we examine the data for trends. If we see positive or negative trends over this time horizon, we are able to conclude that it was the recent increase in the quality of their democracy which led to it. Having investigated examples of this in depth, there seem to be three main archetypes which drive the results. Countries with positive results to their democracy from the election have generally positive effects on their development, countries with more "plateau" like results also did well, but countries for whom the descent to authoritarianism was continued by this election found more negative results.


title: Neural Component Analysis for Fault Detection

abstract: Principal component analysis (PCA) is largely adopted for chemical process monitoring and numerous PCA-based systems have been developed to solve various fault detection and diagnosis problems. Since PCA-based methods assume that the monitored process is linear, nonlinear PCA models, such as autoencoder models and kernel principal component analysis (KPCA), has been proposed and applied to nonlinear process monitoring. However, KPCA-based methods need to perform eigen-decomposition (ED) on the kernel Gram matrix whose dimensions depend on the number of training data. Moreover, prefixed kernel parameters cannot be most effective for different faults which may need different parameters to maximize their respective detection performances. Autoencoder models lack the consideration of orthogonal constraints which is crucial for PCA-based algorithms. To address these problems, this paper proposes a novel nonlinear method, called neural component analysis (NCA), which intends to train a feedforward neural work with orthogonal constraints such as those used in PCA. NCA can adaptively learn its parameters through backpropagation and the dimensionality of the nonlinear features has no relationship with the number of training samples. Extensive experimental results on the Tennessee Eastman (TE) benchmark process show the superiority of NCA in terms of missed detection rate (MDR) and false alarm rate (FAR). The source code of NCA can be found in https://github.com/haitaozhao/Neural-Component-Analysis.git.


title: 200x Low-dose PET Reconstruction using Deep Learning

abstract: Positron emission tomography (PET) is widely used in various clinical applications, including cancer diagnosis, heart disease and neuro disorders. The use of radioactive tracer in PET imaging raises concerns due to the risk of radiation exposure. To minimize this potential risk in PET imaging, efforts have been made to reduce the amount of radio-tracer usage. However, lowing dose results in low Signal-to-Noise-Ratio (SNR) and loss of information, both of which will heavily affect clinical diagnosis. Besides, the ill-conditioning of low-dose PET image reconstruction makes it a difficult problem for iterative reconstruction algorithms. Previous methods proposed are typically complicated and slow, yet still cannot yield satisfactory results at significantly low dose. Here, we propose a deep learning method to resolve this issue with an encoder-decoder residual deep network with concatenate skip connections. Experiments shows the proposed method can reconstruct low-dose PET image to a standard-dose quality with only two-hundredth dose. Different cost functions for training model are explored. Multi-slice input strategy is introduced to provide the network with more structural information and make it more robust to noise. Evaluation on ultra-low-dose clinical data shows that the proposed method can achieve better result than the state-of-the-art methods and reconstruct images with comparable quality using only 0.5% of the original regular dose.


title: GibbsNet: Iterative Adversarial Inference for Deep Graphical Models

abstract: Directed latent variable models that formulate the joint distribution as $p(x,z) = p(z) p(x \mid z)$ have the advantage of fast and exact sampling. However, these models have the weakness of needing to specify $p(z)$, often with a simple fixed prior that limits the expressiveness of the model. Undirected latent variable models discard the requirement that $p(z)$ be specified with a prior, yet sampling from them generally requires an iterative procedure such as blocked Gibbs-sampling that may require many steps to draw samples from the joint distribution $p(x, z)$. We propose a novel approach to learning the joint distribution between the data and a latent code which uses an adversarially learned iterative procedure to gradually refine the joint distribution, $p(x, z)$, to better match with the data distribution on each step. GibbsNet is the best of both worlds both in theory and in practice. Achieving the speed and simplicity of a directed latent variable model, it is guaranteed (assuming the adversarial game reaches the virtual training criteria global minimum) to produce samples from $p(x, z)$ with only a few sampling iterations. Achieving the expressiveness and flexibility of an undirected latent variable model, GibbsNet does away with the need for an explicit $p(z)$ and has the ability to do attribute prediction, class-conditional generation, and joint image-attribute modeling in a single model which is not trained for any of these specific tasks. We show empirically that GibbsNet is able to learn a more complex $p(z)$ and show that this leads to improved inpainting and iterative refinement of $p(x, z)$ for dozens of steps and stable generation without collapse for thousands of steps, despite being trained on only a few steps.


title: On relative Gromov--Witten invariants of projective completions of vector bundles

abstract: It was proved by Fan--Lee and Fan that the absolute Gromov--Witten invariants of two projective bundles $\mathbb P(V_i)\rightarrow X$ are identified canonically when the total Chern classes $c(V_1)=c(V_2)$ for two bundles $V_1$ and $V_2$ over a smooth projective variety $X$. In this note we show that for the two projective completions $\mathbb P(V_i\oplus\mathcal O)$ of $V_i$ and their infinity divisors $\mathbb P(V_i)$, the relative Gromov--Witten invariants of $(\mathbb P(V_i\oplus\mathcal O),\mathbb P(V_i))$ are identified canonically when $c(V_1)=c(V_2)$.


title: Performance guarantees for greedy maximization of non-submodular controllability metrics

abstract: A key problem in emerging complex cyber-physical networks is the design of information and control topologies, including sensor and actuator selection and communication network design. These problems can be posed as combinatorial set function optimization problems to maximize a dynamic performance metric for the network. Some systems and control metrics feature a property called submodularity, which allows simple greedy algorithms to obtain provably near-optimal topology designs. However, many important metrics lack submodularity and therefore lack provable guarantees for using a greedy optimization approach. Here we show that performance guarantees can be obtained for greedy maximization of certain non-submodular functions of the controllability and observability Gramians. Our results are based on two key quantities: the submodularity ratio, which quantifies how far a set function is from being submodular, and the curvature, which quantifies how far a set function is from being supermodular.


title: Search for $B^{-}\to\mu^{-}\bar\nu_\mu$ Decays at the Belle Experiment

abstract: We report the result of a search for the decay $B^{-}\to\mu^{-}\bar\nu_\mu$. The signal events are selected based on the presence of a high momentum muon and the topology of the rest of the event showing properties of a generic $B$-meson decay, as well as the missing energy and momentum being consistent with the hypothesis of a neutrino from the signal decay. We find a 2.4 standard deviation excess above background including systematic uncertainties, which corresponds to a branching fraction of ${\cal B}(B^{-}\to\mu^{-}\bar\nu_\mu) =(6.46 \pm 2.22 \pm 1.60 )\times10^{-7}$ or a frequentist 90% confidence level interval on the $B^{-}\to\mu^{-}\bar\nu_\mu$ branching fraction of $[2.9, 10.7]\times 10^{-7}$. This result is obtained from a $711\ \text{fb}^{-1}$ data sample that contains $772 \times 10^6$ $B\bar{B}$ pairs, collected near the $\Upsilon(4S)$ resonance with the Belle detector at the KEKB asymmetric-energy $e^+ e^-$ collider.


title: Ergodic Capacity of Composite Fading Channels in Cognitive Radios with the Product of $\kappa$-$\mu$ and $\alpha$-$\mu$ Variates

abstract: In this study, the product of two independent and non-identically distributed (i.n.i.d.) random variables (RVs) for \k{appa}-{\mu} fading distribution and {\alpha}-{\mu} fading distribution is considered. The method of the product model of RVs has been widely applied in numerous of communications fields, such as cascaded fading channels, multiple input multiple output (MIMO) systems, radar communications and cognitive radio networks (CRs). The exact series expressions of the product of two i.n.i.d. RVs X for \k{appa}-{\mu} variates and Y for {\alpha}-{\mu} variates are derived instead of Fox H-function to solve the problem that Fox H-function in the RVs product could not be implemented in popular mathematical software packages as Mathematica and Maple. Novel Exact close-form expressions of probability density function (PDF) and cumulative distribution function (CDF) of proposed models are deduced to present the series expressions of product and generalized composite multipath shadowing models. Furthermore, novel exact expressions of the ergodic channel capacity (ECC) are obtained under optimal rate adaptation with constant transmit power (ORA). At last, these analytical results are confirmed with monte-carlo simulations to evaluate spectrum efficiency over generalized composite shadowing fading scenarios in CRs.


title: On homologically locally connected spaces

abstract: We provide some properties and characterizations of homologically $UV^n$-maps and $lc^n_G$-spaces. We show that there is a parallel between recently introduced by Cauty algebraic $ANR$'s and homologically $lc^n_G$-metric spaces, and this parallel is similar to the parallel between ordinary $ANR$'s and $LC^n$-metric spaces. We also show that there is a similarity between the properties of $LC^n$-spaces and $lc^n_G$-spaces. Some open questions are raised.


title: A new scheme for heavy nuclei: proxy-SU(3)

abstract: The SU(3) symmetry realized by J. P. Elliott in the sd nuclear shell is destroyed in heavier shells by the strong spin-orbit interaction. However, the SU(3) symmetry has been used for the description of heavy nuclei in terms of bosons in the framework of the Interacting Boson Approximation, as well as in terms of fermions using the pseudo-SU(3) approximation. We introduce a new fermionic approximation, called the proxy-SU(3), and we discuss how some of its novel predictions come out as a consequence of the short range of the nucleon-nucleon interaction and the Pauli principle.


title: On the existence of a cherry-picking sequence

abstract: Recently, the minimum number of reticulation events that is required to simultaneously embed a collection P of rooted binary phylogenetic trees into a so-called temporal network has been characterized in terms of cherry-picking sequences. Such a sequence is a particular ordering on the leaves of the trees in P. However, it is well-known that not all collections of phylogenetic trees have a cherry-picking sequence. In this paper, we show that the problem of deciding whether or not P has a cherry-picking sequence is NP-complete for when P contains at least eight rooted binary phylogenetic trees. Moreover, we use automata theory to show that the problem can be solved in polynomial time if the number of trees in P and the number of cherries in each such tree are bounded by a constant.


title: Substitution- and Strain-induced Magnetic Phase Transition in Iron Carbide

abstract: Cementite-type carbides are of interest for magnetocaloric applications owing to their temperature- or pressure-induced magnetic phase transition. Here, using first-principles calculations, we investigate the magnetism and the magnetic phase transition in iron carbide (Fe3C) with the substitution of Cr atoms at Fe sites with the strain effect. The presence of Cr atoms is found to give rise to a second-order magnetic phase transition from a ferromagnetic phase for Fe3C to a nonmagnetic phase in chromium carbide (Cr3C).While the ternary Fe2CrC and Cr2FeC compounds prefer the ferrimagnetic ground state, the magnitudes of both the Fe and Cr spin moments, which are antiparallel in orientation, decrease as x increases in Fe3-xCrxC (x = 0, 1, 2, and 3). Furthermore, the fixed spin-moment calculations indicate that the magnetization of Fe3-xCrxC compounds can be delicately altered via the strain effect and that the magnetic-nonmagnetic phase transition occurs at an early stage of Cr substitution, x = 2.