title: BLADYG: A Graph Processing Framework for Large Dynamic Graphs

abstract: Recently, distributed processing of large dynamic graphs has become very popular, especially in certain domains such as social network analysis, Web graph analysis and spatial network analysis. In this context, many distributed/parallel graph processing systems have been proposed, such as Pregel, GraphLab, and Trinity. These systems can be divided into two categories: (1) vertex-centric and (2) block-centric approaches. In vertex-centric approaches, each vertex corresponds to a process, and message are exchanged among vertices. In block-centric approaches, the unit of computation is a block, a connected subgraph of the graph, and message exchanges occur among blocks. In this paper, we are considering the issues of scale and dynamism in the case of block-centric approaches. We present bladyg, a block-centric framework that addresses the issue of dynamism in large-scale graphs. We present an implementation of BLADYG on top of akka framework. We experimentally evaluate the performance of the proposed framework.


title: Comprehensive analysis of large $\alpha$ yields observed in $^{6}$Li induced reactions

abstract: Background: Large $\alpha$ yields have been reported over the years in reactions with $^{6}$Li and $^{7}$Li projectiles. Previous theoretical analyses have shown that the elastic breakup (EBU) mechanism (i.e., projectile breakup leaving the target in its ground state) is able to account only for a small fraction of the total $\alpha$ inclusive breakup cross sections, pointing toward the dominance of non-elastic breakup (NEB) mechanisms. Purpose: We aim to provide a systematic study of the $\alpha$ inclusive cross sections observed in nuclear reactions induced by $^{6}$Li projectiles. In addition to estimating the total $\alpha$ singles cross sections, it is our goal to evaluate angular and energy distributions of these $\alpha$ particles and compare with experimental data, when available. Method: We compute separately the EBU and NEB components of the inclusive breakup cross sections. For the former, we use the continuum-discretized coupled-channels (CDCC) method, which treats this mechanism to all orders. For the NEB part, we employ the the model proposed in the eighties by Ichimura, Austern and Vincent [Phys. Rev. C32, 432 (1982)], within the DWBA approximation. Results: Overall, the sum of the computed EBU and NEB cross sections is found to reproduce very well the measured singles cross sections. In all cases analyzed, we find that the inclusive breakup cross section is largely dominated by the NEB component. Conclusions: The presented method provides a global and systematic description of inclusive breakup reactions induced by $^{6}$Li projectiles. It provides also a natural explanation of the previously observed underestimation of the measured $\alpha$ yields by CDCC calculations. The method used here can be extended to other weakly-bound projectiles, including halo nuclei.


title: Single shot, double differential spectral measurements of inverse Compton scattering in linear and nonlinear regimes

abstract: Inverse Compton scattering (ICS) is a unique mechanism for producing fast pulses - picosecond and below - of bright X- to gamma-rays. These nominally narrow spectral bandwidth electromagnetic radiation pulses are efficiently produced in the interaction between intense, well-focused electron and laser beams. The spectral characteristics of such sources are affected by many experimental parameters, such as the bandwidth of the laser, and the angles of both the electrons and laser photons at collision. The laser field amplitude induces harmonic generation and importantly, for the present work, nonlinear red shifting, both of which dilute the spectral brightness of the radiation. As the applications enabled by this source often depend sensitively on its spectra, it is critical to resolve the details of the wavelength and angular distribution obtained from ICS collisions. With this motivation, we present here an experimental study that greatly improves on previous spectral measurement methods based on X-ray K-edge filters, by implementing a multi-layer bent-crystal X-ray spectrometer. In tandem with a collimating slit, this method reveals a projection of the double-differential angular-wavelength spectrum of the ICS radiation in a single shot. The measurements enabled by this diagnostic illustrate the combined off-axis and nonlinear-field-induced red shifting in the ICS emission process. They reveal in detail the strength of the normalized laser vector potential, and provide a non-destructive measure of the temporal and spatial electron-laser beam overlap.


title: The size of the last merger and time reversal in $\Lambda$-coalescents

abstract: We consider the number of blocks involved in the last merger of a $\Lambda$-coalescent started with $n$ blocks. We give conditions under which, as $n \to \infty$, the sequence of these random variables a) is tight, b) converges in distribution to a finite random variable or c) converges to infinity in probability. Our conditions are optimal for $\Lambda$-coalescents that have a dust component. For general $\Lambda$, we relate the three cases to the existence, uniqueness and non-existence of quasi-invariant measures for the dynamics of the block-counting process, and in case b) investigate the time-reversal of the block-counting process back from the time of the last merger.


title: A Surrogate Model of Gravitational Waveforms from Numerical Relativity Simulations of Precessing Binary Black Hole Mergers

abstract: We present the first surrogate model for gravitational waveforms from the coalescence of precessing binary black holes. We call this surrogate model NRSur4d2s. Our methodology significantly extends recently introduced reduced-order and surrogate modeling techniques, and is capable of directly modeling numerical relativity waveforms without introducing phenomenological assumptions or approximations to general relativity. Motivated by GW150914, LIGO's first detection of gravitational waves from merging black holes, the model is built from a set of $276$ numerical relativity (NR) simulations with mass ratios $q \leq 2$, dimensionless spin magnitudes up to $0.8$, and the restriction that the initial spin of the smaller black hole lies along the axis of orbital angular momentum. It produces waveforms which begin $\sim 30$ gravitational wave cycles before merger and continue through ringdown, and which contain the effects of precession as well as all $\ell \in \{2, 3\}$ spin-weighted spherical-harmonic modes. We perform cross-validation studies to compare the model to NR waveforms \emph{not} used to build the model, and find a better agreement within the parameter range of the model than other, state-of-the-art precessing waveform models, with typical mismatches of $10^{-3}$. We also construct a frequency domain surrogate model (called NRSur4d2s_FDROM) which can be evaluated in $50\, \mathrm{ms}$ and is suitable for performing parameter estimation studies on gravitational wave detections similar to GW150914.


title: The weak rate of convergence for the Euler-Maruyama approximation of one-dimensional stochastic differential equations involving the local times of the unknown process

abstract: In this paper, we consider the weak convergence of the Euler-Maruyama approximation for one dimensional stochastic differential equations involving the local times of the unknown process. We use a transformation in order to remove the local time from the stochastic differential equations and we provide the approximation of Euler-maruyama for the stochastic differential equations without local time. After that, we conclude the approximation of Euler-maruyama for one dimensional stochastic differential equations involving the local times of the unknown process , and we provide the rate of weak convergence for any function G in a certain class.


abstract: In this paper, we introduce and study the concept of L-Dunford-Pettis sets and L-Dunford-Pettis property in Banach spaces. Next, we give a characterization of the L-Dunford-Pettis property with respect to some well-known geometric properties of Banach spaces. Finally, some complementability of operators on Banach spaces with the L-Dunford-Pettis property are also investigated


title: Existence, Stability & Dynamics of Nonlinear Modes in a 2d Partially $\mathcal{PT}$ Symmetric Potential

abstract: It is known that multidimensional complex potentials obeying $\mathcal{PT}$-symmetry may possess all real spectra and continuous families of solitons. Recently it was shown that for multi-dimensional systems these features can persist when the parity symmetry condition is relaxed so that the potential is invariant under reflection in only a single spatial direction. We examine the existence, stability and dynamical properties of localized modes within the cubic nonlinear Schr\"odinger equation in such a scenario of partially $\mathcal{PT}$-symmetric potential.


title: Exponential Convergence of a Distributed Algorithm for Solving Linear Algebraic Equations

abstract: In a recent paper, a distributed algorithm was proposed for solving linear algebraic equations of the form $Ax = b$ assuming that the equation has at least one solution. The equation is presumed to be solved by $m$ agents assuming that each agent knows a subset of the rows of the matrix $[A \; b]$, the current estimates of the equation's solution generated by each of its neighbors, and nothing more. Neighbor relationships are represented by a time-dependent directed graph $N(t)$ whose vertices correspond to agents and whose arcs characterize neighbor relationships. Sufficient conditions on $N(t)$ were derived under which the algorithm can cause all agents' estimates to converge exponentially fast to the same solution to $Ax = b$. These conditions were also shown to be necessary for exponential convergence, provided the data about $[A \; b]$ available to the agents is "non-redundant". The aim of this paper is to relax this "non-redundant" assumption. This is accomplished by establishing exponential convergence under conditions which are the weakest possible for the problem at hand; the conditions are based on a new notion of graph connectivity. An improved bound on the convergence rate is also derived.


title: The null-geodesic flow near horizons

abstract: This note describes the behavior of null-geodesics near nondegenerate Killing horizons in language amenable to the application of a general framework, due to Vasy and Hintz, for the analysis of both linear and nonlinear wave equations. Throughout, the viewpoint of Melrose's b-geometry on a suitable compactification of spacetime at future infinity is adopted.


title: Resynchronization of the Asynchronous Polar CD Ind

abstract: CD Ind is one of only four confirmed asynchronous polars (APs). APs are strongly magnetic cataclysmic variables of the AM Herculis subclass with the characteristic that their white dwarfs rotate a few per cent out of synchronism with their binary orbit. Theory suggests that nova eruptions disrupt previously synchronized states. Following the eruption, the system is expected to rapidly resynchronize over a timescale of centuries. The other three asynchronous polars - V1432 Aql, BY Cam and V1500 Cyg - have resynchronization time estimates ranging from 100 to more than 3500 years, with all but one being less than 1200 years. We report on the analysis of over 46000 observations of CD Ind taken between 2007 and 2016, combined with previous observations from 1996, and estimate a CD Ind resynchronization time of 6400 +/- 800 years. We also estimate an orbital period of 110.820(1) minutes and a current (2016.4) white dwarf spin period of 109.6564(1) minutes.


title: Discrete Optimal Global Convergence of an Evolutionary Algorithm for Clusters under the Potential of Lennard Jones

abstract: A review of the properties that bond the particles under Lennard Jones Potential allow to states properties and conditions for building evolutive algorithms using the CB lattice with other different lattices. The new lattice is called CB lattice and it is based on small cubes. A set of propositions states convergence and optimal conditions over the CB lattice for an evolutionary algorithm. The evolutionary algorithm is a reload version of previous genetic algorithms based in phenotypes. The novelty using CB lattice, together with the other lattices, and ad-hoc cluster segmentation and enumeration, is to allow the combination of genotype (DNA coding for cluster using their particle's number) and phenotype (geometrical shapes using particle's coordinates in 3D). A parallel version of an evolutionary algorithm for determining the global optimality is depicted. The results presented are from a standalone program for a personal computer of the evolutionary algorithm, which can estimate all putative optimal Lennard Jones Clusters from 13 to 1612 particles. The novelty are the theoretical results for the evolutionary algorithm's efficiency, the strategies with phenotype or genotype, and the classification of the clusters based in an ad-hoc geometric algorithm for segmenting a cluster into its nucleus and layers. Also, the standalone program is not only capable to replicate the optimal Lennard Jones clusters in The Cambridge Cluster Database (CCD), but to find new ones.


title: Successive Convexification of Non-Convex Optimal Control Problems with State Constraints

abstract: This paper presents a Successive Convexification ($ \texttt{SCvx} $) algorithm to solve a class of non-convex optimal control problems with certain types of state constraints. Sources of non-convexity may include nonlinear dynamics and non-convex state/control constraints. To tackle the challenge posed by non-convexity, first we utilize exact penalty function to handle the nonlinear dynamics. Then the proposed algorithm successively convexifies the problem via a $ \textit{project-and-linearize} $ procedure. Thus a finite dimensional convex programming subproblem is solved at each succession, which can be done efficiently with fast Interior Point Method (IPM) solvers. Global convergence to a local optimum is demonstrated with certain convexity assumptions, which are satisfied in a broad range of optimal control problems. The proposed algorithm is particularly suitable for solving trajectory planning problems with collision avoidance constraints. Through numerical simulations, we demonstrate that the algorithm converges reliably after only a few successions. Thus with powerful IPM based custom solvers, the algorithm can be implemented onboard for real-time autonomous control applications.


title: Correction of exciton binding energy in monolayer transition metal dichalcogenides

abstract: We theoretically investigate the corrections of exciton binding energy in monolayer transition metal dichalcogenides (TMDs) due to the exciton-optical phonon coupling in the Fr$\ddot{o}$hlich interaction model by using the linear operator combined Lee-Low-Pines variational method. We not only consider the excitons couple with the intrinsic longitudinal optical (LO) phonon modes, but also the surface optical phonon modes that induced by the polar substrates underneath the TMDs. We find that exciton binding energies are corrected in a large scale due to these exciton-optical phonon couplings. We discuss the dependences of exciton binding energy on the cut-off wave vector of optical phonon modes, the polarization parameters of materials and the interlayer distance between the polar substrates and TMDs. These results provide potential explanations for the divergence of the exciton binding energy between experiment and theory in TMDs.


title: Modular representation theory in type A via Soergel bimodules

abstract: In this paper we express certain multiplicities in modular representation-theoretic categories of type A in terms of affine p-Kazhdan-Lusztig polynomials. The representation-theoretic categories we deal with include the categories of rational representations of GL(n), representations of the quantum group for gl(n), and representations of (degenerate) cyclotomic Hecke and Schur algebras, where the base field is an algebraically closed field of arbitrary prime characteristic. In order to approach this problem we define Soergel-theoretic versions of parabolic categories O in characteristic p. We show that these categories have many common features with the classical parabolic categories O; for example, they are highest weight. We produce a homomorphism from a (finite or affine) type A 2-Kac-Moody category to the diagrammatic version of the category of singular Soergel bimodules (again, of finite or affine type A). This leads to a categorical Kac-Moody action on the Soergel-theoretic categories O. Then we relate the representation-theoretic categories to Soergel-theoretic ones by proving a uniqueness result for highest weight categorical actions on Fock spaces.


title: Robust and Real-time Deep Tracking Via Multi-Scale Domain Adaptation

abstract: Visual tracking is a fundamental problem in computer vision. Recently, some deep-learning-based tracking algorithms have been achieving record-breaking performances. However, due to the high complexity of deep learning, most deep trackers suffer from low tracking speed, and thus are impractical in many real-world applications. Some new deep trackers with smaller network structure achieve high efficiency while at the cost of significant decrease on precision. In this paper, we propose to transfer the feature for image classification to the visual tracking domain via convolutional channel reductions. The channel reduction could be simply viewed as an additional convolutional layer with the specific task. It not only extracts useful information for object tracking but also significantly increases the tracking speed. To better accommodate the useful feature of the target in different scales, the adaptation filters are designed with different sizes. The yielded visual tracker is real-time and also illustrates the state-of-the-art accuracies in the experiment involving two well-adopted benchmarks with more than 100 test videos.


title: End-to-End Attention based Text-Dependent Speaker Verification

abstract: A new type of End-to-End system for text-dependent speaker verification is presented in this paper. Previously, using the phonetically discriminative/speaker discriminative DNNs as feature extractors for speaker verification has shown promising results. The extracted frame-level (DNN bottleneck, posterior or d-vector) features are equally weighted and aggregated to compute an utterance-level speaker representation (d-vector or i-vector). In this work we use speaker discriminative CNNs to extract the noise-robust frame-level features. These features are smartly combined to form an utterance-level speaker vector through an attention mechanism. The proposed attention model takes the speaker discriminative information and the phonetic information to learn the weights. The whole system, including the CNN and attention model, is joint optimized using an end-to-end criterion. The training algorithm imitates exactly the evaluation process --- directly mapping a test utterance and a few target speaker utterances into a single verification score. The algorithm can automatically select the most similar impostor for each target speaker to train the network. We demonstrated the effectiveness of the proposed end-to-end system on Windows $10$ "Hey Cortana" speaker verification task.


title: Entropy Formula for Random $\mathbb{Z}^k$-actions

abstract: In this paper, entropies, including measure-theoretic entropy and topological entropy, are considered for random $\mathbb{Z}^k$-actions which are generated by random compositions of the generators of $\mathbb{Z}^k$-actions. Applying Pesin's theory for commutative diffeomorphisms we obtain a measure-theoretic entropy formula of $C^{2}$ random $\mathbb{Z}^k$-actions via the Lyapunov spectra of the generators. Some formulas and bounds of topological entropy for certain random $\mathbb{Z}^k$(or $\mathbb{Z}_+^k$ )-actions generated by more general maps, such as Lipschitz maps, continuous maps on finite graphs and $C^{1}$ expanding maps, are also obtained. Moreover, as an application, we give a formula of Friedland's entropy for certain $C^{2}$ $\mathbb{Z}^k$-actions.


title: Area bound for a surface in a strong gravity region

abstract: For asymptotically flat spacetimes, using the inverse mean curvature flow, we show that any compact $2$-surface, $S_0$, whose mean curvature and its derivative for outward direction are positive in spacelike hypersurface with non-negative Ricci scalar satisfies the inequality $A_0 \leq 4 \pi (3Gm)^2$, where $A_0$ is the area of $S_0$ and $m$ is the total mass. The upper bound is realized when $S_0$ is the photon sphere in a hypersurface isometric to $t=$const. slice of the Schwarzschild spacetime.


title: Persistent Path Homology of Directed Networks

abstract: While standard persistent homology has been successful in extracting information from metric datasets, its applicability to more general data, e.g. directed networks, is hindered by its natural insensitivity to asymmetry. We study a construction of homology of digraphs due to Grigor'yan, Lin, Muranov and Yau, and extend this construction to the persistent framework. The result, which we call persistent path homology, can provide information about the digraph structure of a directed network at varying resolutions. Moreover, this method encodes a rich level of detail about the asymmetric structure of the input directed network. We test our method on both simulated and real-world directed networks and conjecture some of its characteristics.


title: Quantitative stability estimates for Fokker-Planck equations

abstract: We consider the Fokker--Planck equations with irregular coefficients. Two different cases are treated: in the degenerate case, the coefficients are assumed to be weakly differentiable, while in the non-degenerate case the drift satisfies only the Ladyzhenskaya--Prodi--Serrin condition. Using Trevisan's superposition principle which represents the solution as the marginal of the solution to the martingale problem of the diffusion operator, we establish quantitative stability estimates for the solutions of Fokker--Planck equations.


title: Monoidal categories enriched in braided monoidal categories

abstract: We introduce the notion of a monoidal category enriched in a braided monoidal category $\mathcal V$. We set up the basic theory, and prove a classification result in terms of braided oplax monoidal functors to the Drinfeld center of some monoidal category $\mathcal T$. Even the basic theory is interesting; it shares many characteristics with the theory of monoidal categories enriched in a symmetric monoidal category, but lacks some features. Of particular note, there is no cartesian product of braided-enriched categories, and the natural transformations do not form a 2-category, but rather satisfy a braided interchange relation. Strikingly, our classification is slightly more general than what one might have anticipated in terms of strong monoidal functors $\mathcal V \to Z(\mathcal T)$. We would like to understand this further; in a future paper we show that the functor is strong if and only if the enriched category is `complete' in a certain sense. Nevertheless it remains to understand what non-complete enriched categories may look like. One should think of our construction as a generalization of de-equivariantization, which takes a strong monoidal functor $\mathsf{Rep}(G) \to Z(\mathcal T)$ for some finite group $G$ and a monoidal category $\mathcal T$, and produces a new monoidal category $\mathcal T // G$. In our setting, given any braided oplax monoidal functor $\mathcal V \to Z(\mathcal T)$, for any braided $\mathcal V$, we produce $\mathcal T // \mathcal V$: this is not usually an `honest' monoidal category, but is instead $\mathcal V$-enriched. If $\mathcal V$ has a braided lax monoidal functor to $\mathsf{Vec}$, we can use this to reduce the enrichment to $\mathsf{Vec}$, and this recovers de-equivariantization as a special case.


title: Versatility of nodal affiliation to communities

abstract: Graph theoretical analysis of the community structure of networks attempts to identify the communities (or modules) to which each node affiliates. However, this is in most cases an ill-posed problem, as the affiliation of a node to a single community is often ambiguous. Previous solutions have attempted to identify all of the communities to which each node affiliates. Instead of taking this approach, we introduce versatility, $V$, as a novel metric of nodal affiliation: $V \sim 0$ means that a node is consistently assigned to a specific community; $V \gg 0$ means it is inconsistently assigned to different communities. Versatility works in conjunction with existing community detection algorithms, and it satisfies many theoretically desirable properties in idealised networks designed to maximise ambiguity of modular decomposition. The local minima of global mean versatility identified the resolution parameters of a hierarchical community detection algorithm that least ambiguously decomposed the community structure of a social (karate club) network and the mouse brain connectome. Our results suggest that nodal versatility is useful in quantifying the inherent ambiguity of modular decomposition.


title: Slip length of confined liquid with small roughness of solid-liquid interfaces

abstract: We have studied the slip length of confined liquid with small roughness of solid-liquid interfaces. Dyadic Green function and perturbation expansion have been applied to get the slip length quantitatively. In the slip length, both effects of the roughness of the interfaces and the chemical interaction between the liquid and the solid surface are involved. For the numerical calculation, Monte Carlo method has been used to simulate the rough interfaces and the physical quantities are obtained statistically over the interfaces. Results show that the total slip length of the system is linearly proportional to the slip length contributed from the chemical interaction. And the roughness of the interfaces plays its role as the proportionality factor. For the roughness, the variance of the roughness decreases the total slip length while the correlation length of the roughness can enhance the slip length dramatically to a saturation value.


title: Okounkov bodies and transfinite diameter

abstract: We present an explicit calculation of an Okounkov body associated to an algebraic variety. This is used to derive a formula for transfinite diameter on the variety. We relate this formula to a recent result of D. Witt Nystrom.


title: Sutured Manifolds and Polynomial Invariants from Higher Rank Bundles

abstract: For each integer $N\geq 2$, Mari\~no and Moore defined generalized Donaldson invariants by the methods of quantum field theory, and made predictions about the values of these invariants. Subsequently, Kronheimer gave a rigorous definition of generalized Donaldson invariants using the moduli spaces of anti-self-dual connections on hermitian vector bundles of rank $N$. In this paper, Mari\~no and Moore's predictions are confirmed for simply connected elliptic surfaces without multiple fibers and certain surfaces of general type in the case that $N=3$. The primary motivation is to study 3-manifold instanton Floer homologies which are defined by higher rank bundles. In particular, the computation of the generalized Donaldson invariants are exploited to define a Floer homology theory for sutured 3-manifolds.


title: Generalized sigma model with dynamical antisymplectic potential and non-Abelian de Rham's differential

abstract: For topological sigma models, we propose that their local Lagragian density is allowed to depend non-linearly on the de Rham's "velocities" $D Z^{A}$. Then, by differentiating the Lagrangian density with respect to the latter de Rham's "velocities", we define a "dynamical" anti-symplectic potential, in terms of which a "dynamical" anti-symplectic metric is defined, as well. We define the local and the functional antibracket via the dynamical anti-symplectic metric. Finally, we show that the generalized action of the sigma model satisfies the functional master equation, as required.


title: Robust method for finding sparse solutions to linear inverse problems using an L2 regularization

abstract: We analyzed the performance of a biologically inspired algorithm called the Corrected Projections Algorithm (CPA) when a sparseness constraint is required to unambiguously reconstruct an observed signal using atoms from an overcomplete dictionary. By changing the geometry of the estimation problem, CPA gives an analytical expression for a binary variable that indicates the presence or absence of a dictionary atom using an L2 regularizer. The regularized solution can be implemented using an efficient real-time Kalman-filter type of algorithm. The smoother L2 regularization of CPA makes it very robust to noise, and CPA outperforms other methods in identifying known atoms in the presence of strong novel atoms in the signal.


title: Effects of mean-field and the softening of equation of state on elliptic flow in Au+Au collision at $\sqrt{s_{NN}}$ = 5 GeV from JAM model

abstract: We perform a systematic study of elliptic flow ($v_2$) in Au+Au collision at $\sqrt{s_{NN}}=5$ GeV by using a microscopic transport model JAM. The centrality, pseudorapidity, transverse momentum and beam energy dependence of $v_2$ for charged as well as identified hadrons are studied. We investigate the effects of both the hadronic mean-field and the softening of equation of state (EoS) on elliptic flow. The softening of EoS is realized by imposing attractive orbits in two body scattering, which can reduce the pressure of the system. We found that the softening of EoS leads to the enhancement of $v_2$, while the hadronic mean-field suppresses $v_2$ relative to the cascade mode. It indicates that elliptic flow at high baryon density regions is highly sensitive to the EoS and the enhancement of $v_2$ may probe the signature of a first-order phase transition in heavy-ion collisions at beam energies of a strong baryon stopping region.


title: Feasible superadiabatic-based shortcuts for fast generating 3D entanglement between two atoms

abstract: We propose a scheme to realize fast generation of three-dimensional entanglement between two atoms via superadiabatic-based shortcuts in an atom-cavity-fiber system. The scheme is experimentally feasible because of the same form of the counterdiabatic Hamiltonian as that of the effective Hamiltonian. Besides, numerical simulations are given to prove that the scheme is strongly robust against variations in various parameters and decoherence.


title: Shortcut Sequence Tagging

abstract: Deep stacked RNNs are usually hard to train. Adding shortcut connections across different layers is a common way to ease the training of stacked networks. However, extra shortcuts make the recurrent step more complicated. To simply the stacked architecture, we propose a framework called shortcut block, which is a marriage of the gating mechanism and shortcuts, while discarding the self-connected part in LSTM cell. We present extensive empirical experiments showing that this design makes training easy and improves generalization. We propose various shortcut block topologies and compositions to explore its effectiveness. Based on this architecture, we obtain a 6% relatively improvement over the state-of-the-art on CCGbank supertagging dataset. We also get comparable results on POS tagging task.


title: Explaining low l anomalies in the CMB power spectrum with resonant superstring excitations during inflation

abstract: We explore the possibility that both the suppression of the $\ell = 2$ multipole moment of the power spectrum of cosmic microwave background temperature fluctuations and the possible dip for $\ell = 10-30$ can be explained as well as a possible new dip for $\ell \approx 60$ as the result of the resonant creation of sequential excitations of a fermionic (or bosonic) closed superstring that couples to the inflaton field. We consider a D=26 closed bosonic string with one toroidal compact dimension as an illustration of how string excitations might imprint themselves on the CMB. We analyze the existence of successive momentum states, winding states or oscillations on the string as the source of the three possible dips in the power spectrum. Although the evidence of these dips are of marginal statistical significance, this might constitute the first observational evidence of successive superstring excitations in Nature.


title: The Unfolding and Control of Network Cascades

abstract: A characteristic property of networks is their ability to propagate influences, such as infectious diseases, behavioral changes, and failures. An especially important class of such contagious dynamics is that of cascading processes. These processes include, for example, cascading failures in infrastructure systems, extinctions cascades in ecological networks, and information cascades in social systems. In this review, we discuss recent progress and challenges associated with the modeling, prediction, detection, and control of cascades in networks.


title: A multifractal surrogate data generation algorithm that preserves pointwise Holder regularity structure, with initial applications to turbulence

abstract: An algorithm is described that can generate random variants of a time series or image while preserving the probability distribution of original values and the pointwise Holder regularity. Thus, it preserves the multifractal properties of the data. Our algorithm is similar in principle to well-known algorithms based on the preservation of the Fourier amplitude spectrum and original values of a time series. However, it is underpinned by a dual-tree complex wavelet transform rather than a Fourier transform. Our method, which we term the Iterated Amplitude Adjusted Wavelet Transform (IAAWT) method can be used to generate bootstrapped versions of multifractal data and, because it preserves the pointwise Holder regularity but not the local Holder regularity, it can be used to test hypotheses concerning the presence of oscillating singularities in a time series, an important feature of turbulence and econophysics data. Because the locations of the data values are randomized with respect to the multifractal structure, hypotheses about their mutual coupling can be tested, which is important for the velocity-intermittency structure of turbulence and self-regulating processes.


title: On an Enriques surface associated with a quartic Hessian surface

abstract: Let Y be a complex Enriques surface whose universal cover X is birational to a general quartic Hessian surface. Using the result on the automorphism group of X due to Dolgachev and Keum, we obtain a finite presentation of the automorphism group of Y. A fundamental domain of the action of the automorphism group on the nef cone of Y is described explicitly. The list of elliptic fibrations on Y and the list of combinations of rational double points that can appear on a surface birational to Y are presented. As an application, a set of generators of the automorphism group of the generic Enriques surface is calculated explicitly.


title: High Performance Parallel Image Reconstruction for New Vacuum Solar Telescope

abstract: Many technologies have been developed to help improve spatial resolution of observational images for ground-based solar telescopes, such as adaptive optics (AO) systems and post-processing reconstruction. As any AO system correction is only partial, it is indispensable to use post-processing reconstruction techniques. In the New Vacuum Solar Telescope (NVST), speckle masking method is used to achieve the diffraction limited resolution of the telescope. Although the method is very promising, the computation is quite intensive, and the amount of data is tremendous, requiring several months to reconstruct observational data of one day on a high-end computer. To accelerate image reconstruction, we parallelize the program package on a high performance cluster. We describe parallel implementation details for several reconstruction procedures. The code is written in C language using Message Passing Interface (MPI) and optimized for parallel processing in a multi-processor environment. We show the excellent performance of parallel implementation, and the whole data processing speed is about 71 times faster than before. Finally, we analyze the scalability of the code to find possible bottlenecks, and propose several ways to further improve the parallel performance. We conclude that the presented program is capable of executing in real-time reconstruction applications at NVST.


title: Transmission spectra and valley processing of graphene and carbon nanotube superlattices with inter-valley coupling

abstract: We numerically investigate the electronic transport properties of graphene nanoribbons and carbon nanotubes with inter-valley coupling, e.g., in \sqrt{3}N \times \sqrt{3}N and 3N \times 3N superlattices. By taking the \sqrt{3} \times \sqrt{3} graphene superlattice as an example, we show that tailoring the bulk graphene superlattice results in rich structural configurations of nanoribbons and nanotubes. After studying the electronic characteristics of the corresponding armchair and zigzag nanoribbon geometries, we find that the linear bands of carbon nanotubes can lead to the Klein tunnelling-like phenomenon, i.e., electrons propagate along tubes without backscattering even in the presence of a barrier. Due to the coupling between K and K' valleys of pristine graphene by \sqrt{3} \times \sqrt{3} supercells,we propose a valley-field-effect transistor based on the armchair carbon nanotube, where the valley polarization of the current can be tuned by applying a gate voltage or varying the length of the armchair carbon nanotubes.


title: Reply to "Failure to replicate long-range tunable attractions in colloidal system"

abstract: An arxiv paper, ref. [1] by Cao et al., claimed that the tunable attraction reported in our ref. [2] could not be detected. Ref. [1] was submitted to Nature in Apr. 2016 as a Comment on our ref. [2]. Our reply in May 2016 responded to ref. [1] and was reviewed by the editor of Nature and an external referee. Ref. [1] was rejected by Nature in Aug. 2016. We provide our reply, which answers all of the criticisms in ref. [1], in Appendix I. The main text contains a brief reply to the main criticisms.


title: Translocation time of a polymer chain through an energy gradient nanopore

abstract: The translocation time of a polymer chain through an interaction energy gradient nanopore was studied by Monte Carlo simulations and the Fokker-Planck equation with double-absorbing boundary conditions. Both the simulation and calculation revealed three different behaviors for polymer translocation. These behaviors can be explained qualitatively from free-energy landscapes obtained for polymer translocation at different parameters. Results show that the translocation time of a polymer chain through a nanopore can be tuned by suitably designing the interaction energy gradient.


title: A strengthened inequality of Alon-Babai-Suzuki's conjecture on set systems with restricted intersections modulo p

abstract: Let $K=\{k_1,k_2,\ldots,k_r\}$ and $L=\{l_1,l_2,\ldots,l_s\}$ be disjoint subsets of $\{0,1,\ldots,p-1\}$, where $p$ is a prime and $A=\{A_1,A_2,\ldots,A_m\}$ be a family of subsets of $[n]$ such that $|A_i|\pmod{p}\in K$ for all $A_i\in A$ and $|A_i\cap A_j|\pmod{p}\in L$ for $i\ne j$. In 1991, Alon, Babai and Suzuki conjectured that if $n\geq s+\max_{1\leq i\leq r} k_i$, then $|A|\leq {n\choose s}+{n\choose s-1}+\cdots+{n\choose s-r+1}$. In 2000, Qian and Ray-Chaudhuri proved the conjecture under the condition $n\geq 2s-r$. In 2015, Hwang and Kim verified the conjecture of Alon, Babai and Suzuki. In this paper, we will prove that if $n\geq 2s-2r+1$ or $n\geq s+\max_{1\leq i\leq r}k_i$, then \[ |A|\leq{n-1\choose s}+{n-1\choose s-1}+\cdots+{n-1\choose s-2r+1}. \] This result strengthens the upper bound of Alon, Babai and Suzuki's conjecture when $n\geq 2s-2$.


title: Characterizing K2 Candidate Planetary Systems Orbiting Low-Mass Stars I: Classifying Low-mass Host Stars Observed During Campaigns 1-7

abstract: We present near-infrared spectra for 144 candidate planetary systems identified during Campaigns 1-7 of the NASA K2 Mission. The goal of the survey was to characterize planets orbiting low-mass stars, but our IRTF/SpeX and Palomar/TripleSpec spectroscopic observations revealed that 49% of our targets were actually giant stars or hotter dwarfs reddened by interstellar extinction. For the 72 stars with spectra consistent with classification as cool dwarfs (spectral types K3 - M4), we refined their stellar properties by applying empirical relations based on stars with interferometric radius measurements. Although our revised temperatures are generally consistent with those reported in the Ecliptic Plane Input Catalog (EPIC), our revised stellar radii are typically 0.13 solar radii (39%) larger than the EPIC values, which were based on model isochrones that have been shown to underestimate the radii of cool dwarfs. Our improved stellar characterizations will enable more efficient prioritization of K2 targets for follow-up studies.


title: Interrogating the Escherichia coli cell cycle by cell dimension perturbations

abstract: Bacteria tightly regulate and coordinate the various events in their cell cycles to duplicate themselves accurately and to control their cell sizes. Growth of Escherichia coli, in particular, follows a relation known as Schaechter 's growth law. This law says that the average cell volume scales exponentially with growth rate, with a scaling exponent equal to the time from initiation of a round of DNA replication to the cell division at which the corresponding sister chromosomes segregate. Here, we sought to test the robustness of the growth law to systematic perturbations in cell dimensions achieved by varying the expression levels of mreB and ftsZ. We found that decreasing the mreB level resulted in increased cell width, with little change in cell length, whereas decreasing the ftsZ level resulted in increased cell length. Furthermore, the time from replication termination to cell division increased with the perturbed dimension in both cases. Moreover, the growth law remained valid over a range of growth conditions and dimension perturbations. The growth law can be quantitatively interpreted as a consequence of a tight coupling of cell division to replication initiation. Thus, its robustness to perturbations in cell dimensions strongly supports models in which the timing of replication initiation governs that of cell division, and cell volume is the key phenomenological variable governing the timing of replication initiation. These conclusions are discussed in the context of our recently proposed adder-per-origin model, in which cells add a constant volume per origin between initiations and divide a constant time after initiation.


title: Stellar & Planetary Parameters for K2's Late Type Dwarf Systems from C1 to C5

abstract: The NASA K2 mission uses photometry to find planets transiting stars of various types. M dwarfs are of high interest since they host more short period planets than any other type of main sequence stars and transiting planets around M dwarfs have deeper transits compared to other main sequence stars. In this paper, we present stellar parameters from K and M dwarfs hosting transiting planet candidates discovered by our team. Using the SOFI spectrograph on the European Southern Observatory's New Technology Telescope, we obtained R ~ 1000 J-, H-, and K-band (0.95 - 2.52 microns) spectra of 34 late-type K2 planet and candidate planet host systems and 12 bright K4-M5 dwarfs with interferometrically measured radii and effective temperatures. Out of our 34 late-type K2 targets, we identify 27 of these stars as M dwarfs. We measure equivalent widths of spectral features, derive calibration relations using stars with interferometric measurements, and estimate stellar radii, effective temperatures, masses, and luminosities for the K2 planet hosts. Our calibrations provide radii and temperatures with median uncertainties of 0.059 R_sun (16.09%) and 160 K (4.33%), respectively. We then reassess the radii and equilibrium temperatures of known and candidate planets based on our spectroscopically derived stellar parameters. Since a planet's radius and equilibrium temperature depends on the parameters of its host star, our study provides more precise planetary parameters for planets and candidates orbiting late-type stars observed with K2. We find a median planet radius and an equilibrium temperature of approximately 3R_earth and 500 K, with several systems (K2-18b and K2-72e) receiving near-Earth like levels of incident irradiation.


title: Two-photon photoassociation spectroscopy of an ultracold heteronuclear molecule

abstract: We report on two-photon photoassociation (PA) spectroscopy of ultracold heteronuclear LiRb molecules. This is used to determine the binding energies of the loosely bound levels of the electronic ground singlet and the lowest triplet states of LiRb. We observe strong two-photon PA lines with power broadened line widths greater than 20 GHz at relatively low laser intensity of 30 W/cm$^{2}$. The implication of this observation on direct atom to molecule conversion using stimulated Raman adiabatic passage (STIRAP) is discussed and the prospect for electronic ground state molecule production is theoretically analyzed.

categories: cond-mat.quant-gas physics.atom-ph physics.chem-ph physics.optics quant-ph


title: Chemical Composition of Young Stars in the Leading Arm of the Magellanic System

abstract: Chemical abundances of eight O- and B-type stars are determined from high-resolution spectra obtained with the MIKE instrument on the Magellan 6.5m Clay telescope. The sample is selected from 42 candidates of membership in the Leading Arm of the Magellanic System. Stellar parameters are measured by two independent grids of model atmospheres and analysis procedures, confirming the consistency of the stellar parameter results. Abundances of seven elements (He, C, N, O, Mg, Si, and S) are determined for the stars, as are their radial velocities and estimates of distances and ages. Among the seven B-type stars analyzed, the five that have radial velocities compatible with membership to the LA have an average [Mg/H] of $-0.42\pm0.16$, significantly lower than the average of the remaining two [Mg/H] = $-0.07\pm0.06$ that are kinematical members of the Galactic disk. Among the five LA members, four have individual [Mg/H] abundance compatible with that in the LMC. Within errors, we can not exclude the possibility that one of these stars has a [Mg/H] consistent with the more metal-poor, SMC-like material. The remaining fifth star has a [Mg/H] close to MW values. Distances to the LA members indicate that they are at the edge of the Galactic disk, while ages are of the order of $\sim 50-70$ Myr, lower than the dynamical age of the LA, suggesting a single star-forming episode in the LA. V$_{\rm LSR}$ the LA members decreases with decreasing Magellanic longitude, confirming the results of previous LA gas studies.


title: Andreev reflection in 2D relativistic materials with realistic tunneling transparency in normal-metal-superconductor junctions

abstract: The Andreev conductance across 2d normal metal (N)/superconductor (SC) junctions with relativistic Dirac spectrum is investigated theoretically in the Blonder-Tinkham-Klapwijk formalism. It is shown that for relativistic materials, due to the Klein tunneling instead of impurity potentials, the local strain in the junction is the key factor that determines the transparency of the junction. The local strain is shown to generate an effective Dirac $\delta$-gauge field. A remarkable suppression of the conductance are observed as the strength of the gauge field increases. The behaviors of the conductance are in well agreement with the results obtained in the case of 1d N/SC junction. We also study the Andreev reflection in a topological material near the chiral-to-helical phase transition in the presence of a local strain. The N side of the N/SC junction is modeled by the doped Kane-Mele (KM) model. The SC region is a doped correlated KM t-J (KMtJ) model, which has been shown to feature d+id'-wave spin-singlet pairing. With increasing intrinsic spin-orbit (SO) coupling, the doped KMtJ system undergoes a topological phase transition from the chiral d-wave superconductivity to the spin-Chern superconducting phase with helical Majorana fermions at edges. We explore the Andreev conductance at the two inequivalent Dirac points, respectively and predict the distinctive behaviors for the Andreev conductance across the topological phase transition. Relevance of our results for the adatom-doped graphene is discussed.


title: Comment on Jackson's analysis of electric charge quantization due to interaction with Dirac's magnetic monopole

abstract: In J.D. Jackson's Classical Electrodynamics textbook, the analysis of Dirac's charge quantization condition in the presence of a magnetic monopole has a mathematical omission and an all too brief physical argument that might mislead some students. This paper presents a detailed derivation of Jackson's main result, explains the significance of the missing term, and highlights the close connection between Jackson's findings and Dirac's original argument.


title: HLA class I binding prediction via convolutional neural networks

abstract: Many biological processes are governed by protein-ligand interactions. One such example is the recognition of self and nonself cells by the immune system. This immune response process is regulated by the major histocompatibility complex (MHC) protein which is encoded by the human leukocyte antigen (HLA) complex. Understanding the binding potential between MHC and peptides can lead to the design of more potent, peptide-based vaccines and immunotherapies for infectious autoimmune diseases. We apply machine learning techniques from the natural language processing (NLP) domain to address the task of MHC-peptide binding prediction. More specifically, we introduce a new distributed representation of amino acids, name HLA-Vec, that can be used for a variety of downstream proteomic machine learning tasks. We then propose a deep convolutional neural network architecture, name HLA-CNN, for the task of HLA class I-peptide binding prediction. Experimental results show combining the new distributed representation with our HLA-CNN architecture achieves state-of-the-art results in the majority of the latest two Immune Epitope Database (IEDB) weekly automated benchmark datasets. We further apply our model to predict binding on the human genome and identify 15 genes with potential for self binding.


title: Reliability of Critical Infrastructure Networks: Challenges

abstract: Critical infrastructures form a technological skeleton of our world by providing us with water, food, electricity, gas, transportation, communication, banking, and finance. Moreover, as urban population increases, the role of infrastructures become more vital. In this paper, we adopt a network perspective and discuss the ever growing need for fundamental interdisciplinary study of critical infrastructure networks, efficient methods for estimating their reliability, and cost-effective strategies for enhancing their resiliency. We also highlight some of the main challenges arising on this way, including cascading failures, feedback loops, and cross-sector interdependencies.


title: Leveraging Multi-aspect Time-related Influence in Location Recommendation

abstract: Point-Of-Interest (POI) recommendation aims to mine a user's visiting history and find her/his potentially preferred places. Although location recommendation methods have been studied and improved pervasively, the challenges w.r.t employing various influences including temporal aspect still remain. Inspired by the fact that time includes numerous granular slots (e.g. minute, hour, day, week and etc.), in this paper, we define a new problem to perform recommendation through exploiting all diversified temporal factors. In particular, we argue that most existing methods only focus on a limited number of time-related features and neglect others. Furthermore, considering a specific granularity (e.g. time of a day) in recommendation cannot always apply to each user or each dataset. To address the challenges, we propose a probabilistic generative model, named after Multi-aspect Time-related Influence (MATI) to promote POI recommendation. We also develop a novel optimization algorithm based on Expectation Maximization (EM). Our MATI model firstly detects a user's temporal multivariate orientation using her check-in log in Location-based Social Networks(LBSNs). It then performs recommendation using temporal correlations between the user and proposed locations. Our method is adaptable to various types of recommendation systems and can work efficiently in multiple time-scales. Extensive experimental results on two large-scale LBSN datasets verify the effectiveness of our method over other competitors.


title: Constraints on the total coupling strength to bosons in iron based superconductors

abstract: At present, there is still no consistent interpretation of the normal and superconducting properties of Fe-based superconductors (FeSCs). The strength of the el-el interaction and the role of correlation effects are under debate. Here, we examine several common materials and illustrate various problems and concepts that are generic for all FeSCs. Based on empirical observations and qualitative insight from density functional theory, we show that the superconducting and low-energy thermodynamic properties of the FeSCs can be described semi-quantitively within multiband Eliashberg theory. We account for an important high-energy mass renormalization phenomenologically,and in agreement with constraints provided by thermodynamic, optical, and angle-resolved photoemission data. When seen in this way, all FeSCs with $T_\mathrm{c} <$ 40~K studied so far are found to belong to an {\it intermediate} coupling regime. This finding is in contrast to the strong coupling scenarios proposed in the early period of the FeSC history.We also discuss several related issues, including the role of band shifts as measured by the positions of van Hove singularities, and the nature of a recently suggested quantum critical point in the strongly hole-doped systems AFe$_2$As$_2$ (A = K, Rb, Cs). Using high-precision full relativistic GGA-band structure calculations, we arrive at a somewhat milder mass renormalization in comparison with previous studies. From the calculated mass anisotropies of all Fermi surface sheets, only the $\varepsilon$-pocket near the corner of the BZ is compatible with the experimentally observed anisotropy of the upper critical field. pointing to its dominant role in the superconductivity of these three compounds.


title: Deep Convolutional Neural Networks for Pairwise Causality

abstract: Discovering causal models from observational and interventional data is an important first step preceding what-if analysis or counterfactual reasoning. As has been shown before, the direction of pairwise causal relations can, under certain conditions, be inferred from observational data via standard gradient-boosted classifiers (GBC) using carefully engineered statistical features. In this paper we apply deep convolutional neural networks (CNNs) to this problem by plotting attribute pairs as 2-D scatter plots that are fed to the CNN as images. We evaluate our approach on the 'Cause- Effect Pairs' NIPS 2013 Data Challenge. We observe that a weighted ensemble of CNN with the earlier GBC approach yields significant improvement. Further, we observe that when less training data is available, our approach performs better than the GBC based approach suggesting that CNN models pre-trained to determine the direction of pairwise causal direction could have wider applicability in causal discovery and enabling what-if or counterfactual analysis.


title: Scalegenesis via dynamically induced multiple seesaws

abstract: We propose a model which accounts for the dynamical origin of the electroweak symmetry breaking (EWSB), directly linking to the mass generation of dark matter (DM) candidate and active neutrinos. The standard model (SM) is weakly charged under the $U(1)_{\rm B-L}$ gauge symmetry, in conjunction with newly introduced three right-handed Majorana neutrinos and the $U(1)_{\rm B-L}$ Higgs. The model is built on the classical scale invariance, that is dynamically broken by a new strongly coupled sector, what is called the hypercolor (HC) sector, which is also weakly coupled to the $U(1)_{\rm B-L}$ gauge. At the HC strong scale, the simultaneous breaking of the EW and $U(1)_{\rm B-L}$ gauge symmetries is triggered by dynamically induced multiple seesaw mechanisms, namely bosonic seesaw mechanisms. Thus, all the origins of masses are provided singly by the HC dynamics: that is what we call the dynamical scalegenesis. We also find that an HC baryon, with mass on the order of a few TeV, can be stabilized by the HC baryon number and the $U(1)_{\rm B-L}$ charge, so identified as a DM candidate. The relic abundance of the HC-baryon DM can be produced dominantly via the bosonic-seesaw portal process, and the HC-baryon DM can be measured through the large magnetic moment coupling generated from the HC dynamics, or the $U(1)_{\rm B-L}$-gauge boson portal in direct detection experiments.


title: AENet: Learning Deep Audio Features for Video Analysis

abstract: We propose a new deep network for audio event recognition, called AENet. In contrast to speech, sounds coming from audio events may be produced by a wide variety of sources. Furthermore, distinguishing them often requires analyzing an extended time period due to the lack of clear sub-word units that are present in speech. In order to incorporate this long-time frequency structure of audio events, we introduce a convolutional neural network (CNN) operating on a large temporal input. In contrast to previous works this allows us to train an audio event detection system end-to-end. The combination of our network architecture and a novel data augmentation outperforms previous methods for audio event detection by 16%. Furthermore, we perform transfer learning and show that our model learnt generic audio features, similar to the way CNNs learn generic features on vision tasks. In video analysis, combining visual features and traditional audio features such as MFCC typically only leads to marginal improvements. Instead, combining visual features with our AENet features, which can be computed efficiently on a GPU, leads to significant performance improvements on action recognition and video highlight detection. In video highlight detection, our audio features improve the performance by more than 8% over visual features alone.


title: On xD-Generalizations of Stirling Numbers and Lah Numbers via Graphs and Rooks

abstract: This paper studies the generalizations of the Stirling numbers of both kinds and the Lah numbers in association with the normal order problem in the Weyl algebra $W=\langle x,D|Dx-xD=1\rangle$. Any word $\omega\in W$ with $m$ $x$'s and $n$ $D$'s can be expressed in the normally ordered form $\omega=x^{m-n}\sum_{k\ge 0} {{\omega}\brace {k}} x^{k}D^{k}$, where ${{\omega}\brace {k}}$ is known as the Stirling number of the second kind for the word $\omega$. This study considers the expansions of restricted words $\omega$ in $W$ over the sequences $\{(xD)^{k}\}_{k\ge 0}$ and $\{xD^{k}x^{k-1}\}_{k\ge 0}$. Interestingly, the coefficients in individual expansions turn out to be generalizations of the Stirling numbers of the first kind and the Lah numbers. The coefficients will be determined through enumerations of some combinatorial structures linked to the words $\omega$, involving decreasing forest decompositions of quasi-threshold graphs and non-attacking rook placements on Ferrers boards. Extended to $q$-analogues, weighted refinements of the combinatorial interpretations are also investigated for words in the $q$-deformed Weyl algebra.


title: The gauge fixing theorem with applications to the Yang-Mills flow over Riemannian manifolds

abstract: In 1982, Uhlenbeck \cite {U2} established the well-known gauge fixing theorem, which has played a fundamental role for Yang-Mills theory. In this paper, we apply the idea of Uhlenbeck to establish a parabolic type of gauge fixing theorems for the Yang-Mills flow and prove existence of a weak solution of the Yang-Mills flow on a compact $n$-dimensional manifold with initial value $A_0$ in $W^{1,n/2}(M)$. When $n=4$, we improve a key lemma of Uhlenbeck (Lemma 2.7 of \cite {U2}) to prove uniqueness of weak solutions of the Yang-Mills flow on a four dimensional manifold.


title: A Molecular Implementation of the Least Mean Squares Estimator

abstract: In order to function reliably, synthetic molecular circuits require mechanisms that allow them to adapt to environmental disturbances. Least mean squares (LMS) schemes, such as commonly encountered in signal processing and control, provide a powerful means to accomplish that goal. In this paper we show how the traditional LMS algorithm can be implemented at the molecular level using only a few elementary biomolecular reactions. We demonstrate our approach using several simulation studies and discuss its relevance to synthetic biology.


title: Gravitational Instabilities of the Cosmic Neutrino Background with Non-zero Lepton Number

abstract: We argue that a cosmic neutrino background that carries non-zero lepton charge develops gravitational instabilities. Fundamentally, these instabilities are related to the mixed gravity-lepton number anomaly. We have explicitly computed the gravitational Chern-Simons term which is generated quantum-mechanically in the effective action in the presence of a lepton number asymmetric neutrino background. The induced Chern-Simons term has a twofold effect: (i) gravitational waves propagating in such a neutrino background exhibit birefringent behaviour leading to an enhancement/suppression of the gravitational wave amplitudes depending on the polarisation, where the magnitude of this effect is related to the size of the lepton asymmetry; (ii) Negative energy graviton modes are induced in the high frequency regime, which leads to very fast vacuum decay producing, e.g., positive energy photons and negative energy gravitons. From the constraint on the present radiation energy density, we obtain an interesting bound on the lepton asymmetry of the universe.


title: Electronic correlations of element cerium under high pressure

abstract: The electronic structures of element cerium under high pressure remain unclear all the time. We tried to calculate the electronic structures of $\alpha'$, $\alpha"$, and $\epsilon$-Ce which only exist in the presence of pressure, by using the combination of traditional density functional theory and single-site dynamical mean-field theory. The momentum-resolved band structures, total and partial density of states, 4$f$ electronic configurations of these phases were exhaustively studied. We found that the 4$f$ electrons tend to be itinerant, and the hybridizations between the 4$f$ and $spd$ conduction electrons are remarkable. In addition, the fluctuations among the 4$f$ atomic eigenstates are prominent, especially for the $\epsilon$ phase, which leads to a slight modification of its 4$f$ occupancy.


title: On monomial linearisation and supercharacters of pattern subgroups

abstract: Column closed pattern subgroups $U$ of the finite upper unitriangular groups $U_n(q)$ are defined as sets of matrices in $U_n(q)$ having zeros in a prescribed set of columns besides the diagonal ones. We explain Jedlitschky's construction of monomial linearisation and apply this to $C U$ yielding a generalisation of Yan's coadjoint cluster representations. Then we give a complete classification of the resulting supercharacters, by describing the resulting orbits and determining the Hom-spaces between orbit modules.


title: Witnessing nonclassical correlations via a single-shot experiment on an ensemble of spins using NMR

abstract: A bipartite quantum system in a mixed state can exhibit nonclassical correlations, which can go beyond quantum entanglement. While quantum discord is the standard measure of quantifying such general quantum correlations, the nonclassicality can be determined by simpler means via the measurement of witness operators. We experimentally construct a positive map to witness nonclassicality of two qubits in an NMR system. The map can be decomposed so that a single run of an experiment on an ensemble of spins suffices to detect the nonclassicality in the state, if present. We let the state evolve in time and use the map to detect nonclassicality as a function of time. To evaluate the efficacy of the witness operator as a means to detect nonclassicality, we measure quantum discord by performing full quantum state tomography at each time point and obtained a fairly good match between the two methods.


title: Self-consistent hyperfluid

abstract: Hyperfluid model is reconstructed on the basis of its action free from any external constraints, regarding the hyperfluid as a self-consistent classical field. Intrinsic hypermomentum is no more a given variable, but arises purely from the diffeomorphism covariance of the dynamical field, which allows the hypermomentum to be exchanged via fundamental interactions of various classes; fluid-fluid interactions, SU(n)-gauge coupling with non-minimal interaction, and coupling with the metric-affine gravity are all successfully formulated in the classical regime.


title: The vanishing cycles of curves in toric surfaces I

abstract: This article is the first in a series of two in which we study the vanishing cycles of curves in toric surfaces. We give a list of possible obstructions to contract vanishing cycles within a given complete linear system. Using tropical means, we show that any non-separating simple closed curve is a vanishing cycle whenever none of the listed obstructions appears.


title: Akid: A Library for Neural Network Research and Production from a Dataism Approach

abstract: Neural networks are a revolutionary but immature technique that is fast evolving and heavily relies on data. To benefit from the newest development and newly available data, we want the gap between research and production as small as possibly. On the other hand, differing from traditional machine learning models, neural network is not just yet another statistic model, but a model for the natural processing engine --- the brain. In this work, we describe a neural network library named {\texttt akid}. It provides higher level of abstraction for entities (abstracted as blocks) in nature upon the abstraction done on signals (abstracted as tensors) by Tensorflow, characterizing the dataism observation that all entities in nature processes input and emit out in some ways. It includes a full stack of software that provides abstraction to let researchers focus on research instead of implementation, while at the same time the developed program can also be put into production seamlessly in a distributed environment, and be production ready. At the top application stack, it provides out-of-box tools for neural network applications. Lower down, akid provides a programming paradigm that lets user easily build customized models. The distributed computing stack handles the concurrency and communication, thus letting models be trained or deployed to a single GPU, multiple GPUs, or a distributed environment without affecting how a model is specified in the programming paradigm stack. Lastly, the distributed deployment stack handles how the distributed computing is deployed, thus decoupling the research prototype environment with the actual production environment, and is able to dynamically allocate computing resources, so development (Devs) and operations (Ops) could be separated. Please refer to http://akid.readthedocs.io/en/latest/ for documentation.


title: Change detection in complex dynamical systems using intrinsic phase and amplitude synchronization

abstract: We present an approach for the detection of sharp change points (short-lived and persistent) in nonlinear and nonstationary dynamic systems under high levels of noise by tracking the local phase and amplitude synchronization among the components of a univariate time series signal. The signal components are derived via Intrinsic Time scale Decomposition (ITD)--a nonlinear, non-parametric analysis method. We show that the signatures of sharp change points are retained across multiple ITD components with a significantly higher probability as compared to random signal fluctuations. Theoretical results are presented to show that combining the change point information retained across a specific set of ITD components offers the possibility of detecting sharp transitions with high specificity and sensitivity. Subsequently, we introduce a concept of mutual agreement to identify the set of ITD components that are most likely to capture the information about dynamical changes of interest and define an InSync statistic to capture this local information. Extensive numerical, as well as real-world case studies involving benchmark neurophysiological processes and industrial machine sensor data, suggest that the present method can detect sharp change points, on an average 62% earlier (in terms of average run length) as compared to other contemporary methods tested.


title: Eichler-Shimura isomorphism and group cohomology on arithmetic groups

abstract: The Eichler-Shimura isomorphism realizes the automorphic representation generated by an automorphic newform in certain cohomology of an arithmetic group. In this short note, we give a cohomological interpretation of the Eichler-Shimura isomorphism as a connection morphism of certain exact sequence of G-modules, for the multiplicative group G of any quaternion algebra over a totally real field.


title: Citation indices and dimensional homogeneity

abstract: The importance of dimensional analysis and dimensional homogeneity in bibliometric studies is always overlooked. In this paper, we look at this issue systematically and show that most h-type indices have the dimensions of [P], where [P] is the basic dimensional unit in bibliometrics which is the unit publication or paper. The newly introduced Euclidean index, based on the Euclidean length of the citation vector has the dimensions [P3/2]. An empirical example is used to illustrate the concepts.


title: Phase diagrams of the extended Bose-Hubbard model in one dimension by Monte-Carlo simulation with the help of a stochastic-series expansion

abstract: In this paper, we study phase diagrams of the extended Bose-Hubbard model (EBHM) in one dimension by means of the quantum Monte-Carlo (QMC) simulation using the stochastic-series expansion (SSE).In the EBHM, there exists a nearest-neighbor repulsion as well as the on-site repulsion. In the SSE-QMC simulation, the highest particle number at each site, $n_c$, is also a controllable parameter, and we found that the phase diagrams depend on the value of $n_c$. It is shown that in addition to the Mott insulator, superfluid, density wave, the phase so-called Haldane insulator and supersolid appear in the phase diagrams, and their locations in the phase diagrams are clarified.


title: Coloring graphs of various maximum degree from random lists

abstract: Let $G=G(n)$ be a graph on $n$ vertices with maximum degree $\Delta=\Delta(n)$. Assign to each vertex $v$ of $G$ a list $L(v)$ of colors by choosing each list independently and uniformly at random from all $k$-subsets of a color set $\mathcal{C}$ of size $\sigma= \sigma(n)$. Such a list assignment is called a \emph{random $(k,\mathcal{C})$-list assignment}. In this paper, we are interested in determining the asymptotic probability (as $n \to \infty$) of the existence of a proper coloring $\varphi$ of $G$, such that $\varphi(v) \in L(v)$ for every vertex $v$ of $G$, a so-called $L$-coloring. We give various lower bounds on $\sigma$, in terms of $n$, $k$ and $\Delta$, which ensures that with probability tending to 1 as $n \to \infty$ there is an $L$-coloring of $G$. In particular, we show, for all fixed $k$ and growing $n$, that if $\sigma(n) = \omega(n^{1/k^2} \Delta^{1/k})$ and $\Delta=O\left(n^{\frac{k-1}{k(k^3+ 2k^2 - k +1)}}\right)$, then the probability that $G$ has an $L$-coloring tends to 1 as $n \rightarrow \infty$. If $k\geq 2$ and $\Delta= \Omega(n^{1/2})$, then the same conclusion holds provided that $\sigma=\omega(\Delta)$. We also give related results for other bounds on $\Delta$, when $k$ is constant or a strictly increasing function of $n$.


title: More satellites around the nearby spiral galaxy NGC 7331

abstract: We report the detection of two dwarf galaxies in a projected distance of ~50 kpc from NGC 7331 and suspect the physical nature of dwarfs of this spiral galaxy.


title: Multivariable Conformable Fractional Calculus

abstract: Conformable fractional derivative is introduced by the authors Khalil et al. In this study we develop their concept and introduce multivariable conformable derivative for a vector valued function with several variables.


title: Mean field limit for survival probability of the high-dimensional contact process

abstract: In this paper we are concerned with the contact process on the squared lattice. The contact process intuitively describes the spread of the infectious disease on a graph, where an infectious vertex becomes healthy at a constant rate while a healthy vertex is infected at rate proportional to the number of infectious neighbors. As the dimension of the lattice grows to infinity, we give a mean field limit for the survival probability of the process conditioned on the event that only the origin of the lattice is infected at t=0. The binary contact path process is a main auxiliary tool for our proof.


title: First-principles modeling of the Invar effect in Fe65Ni35 by the spin-wave method

abstract: Thermal lattice expansion of the Invar Fe65Ni35 alloy is investigated in first-principles calculations using the spin-wave method, which is generalized here for the ferromagnetic state with short range order. It is shown that magnetic short-range order effects make substantial contribution to the equilibrium lattice constant and cannot be neglected in the accurate ab initio modeling of the thermal expansion in Fe-Ni alloys. We also demonstrate that at high temperatures, close and above the magnetic transition, magnetic entropy associated with transverse and longitudinal spin fluctuations yields a noticeable contribution to the equilibrium lattice constant. The obtained theoretical results for the temperature dependent lattice constant are in semiquantitative agreement with the experimental data apart from the region close the magnetic transition.


title: Diagnosing Coronal Heating Processes with Spectrally Resolved Soft X-ray Measurements

abstract: Decades of astrophysical observations have convincingly shown that soft X-ray (SXR; ~0.1--10 keV) emission provides unique diagnostics for the high temperature plasmas observed in solar flares and active regions. SXR observations critical for constraining models of energy release in these phenomena can be provided using instruments that have already been flown on sounding rockets and CubeSats, including miniaturized high-resolution photon-counting spectrometers and a novel diffractive spectral imager. These instruments have relatively low cost and high TRL, and would complement a wide range of mission concepts. In this white paper, we detail the scientific background and open questions motivating these instruments, the measurements required, and the instruments themselves that will make groundbreaking progress in answering these questions.

categories: astro-ph.IM astro-ph.HE astro-ph.SR physics.plasm-ph physics.space-ph


title: Vertical perimeter versus horizontal perimeter

abstract: The discrete Heisenberg group $\mathbb{H}_{\mathbb{Z}}^{2k+1}$ is the group generated by $a_1,b_1,\ldots,a_k,b_k,c$, subject to the relations $[a_1,b_1]=\ldots=[a_k,b_k]=c$ and $[a_i,a_j]=[b_i,b_j]=[a_i,b_j]=[a_i,c]=[b_i,c]=1$ for every distinct $i,j\in \{1,\ldots,k\}$. Denote $S=\{a_1^{\pm 1},b_1^{\pm 1},\ldots,a_k^{\pm 1},b_k^{\pm 1}\}$. The horizontal boundary of $\Omega\subset \mathbb{H}_{\mathbb{Z}}^{2k+1}$, denoted $\partial_{h}\Omega$, is the set of all $(x,y)\in \Omega\times (\mathbb{H}_{\mathbb{Z}}^{2k+1}\setminus \Omega)$ such that $x^{-1}y\in S$. The horizontal perimeter of $\Omega$ is $|\partial_{h}\Omega|$. For $t\in \mathbb{N}$, define $\partial^t_{v} \Omega$ to be the set of all $(x,y)\in \Omega\times (\mathbb{H}_{\mathsf{Z}}^{2k+1}\setminus \Omega)$ such that $x^{-1}y\in \{c^t,c^{-t}\}$. The vertical perimeter of $\Omega$ is defined by $|\partial_{v}\Omega|= \sqrt{\sum_{t=1}^\infty |\partial^t_{v}\Omega|^2/t^2}$. It is shown here that if $k\ge 2$, then $|\partial_{v}\Omega|\lesssim \frac{1}{k} |\partial_{h}\Omega|$. The proof of this "vertical versus horizontal isoperimetric inequality" uses a new structural result that decomposes sets of finite perimeter in the Heisenberg group into pieces that admit an "intrinsic corona decomposition." This allows one to deduce an endpoint $W^{1,1}\to L_2(L_1)$ boundedness of a certain singular integral operator from a corresponding lower-dimensional $W^{1,2}\to L_2(L_2)$ boundedness. The above inequality has several applications, including that any embedding into $L_1$ of a ball of radius $n$ in the word metric on $\mathbb{H}_{\mathbb{Z}}^{5}$ incurs bi-Lipschitz distortion that is at least a constant multiple of $\sqrt{\log n}$. It follows that the integrality gap of the Goemans--Linial semidefinite program for the Sparsest Cut Problem on inputs of size $n$ is at least a constant multiple of $\sqrt{\log n}$.


abstract: In this paper, we investigate *-DMP elements in $*$-semigroups and $*$-rings. The notion of *-DMP element was introduced by Patr\'{i}cio in 2004. An element $a$ is *-DMP if there exists a positive integer $m$ such that $a^{m}$ is EP. We first characterize *-DMP elements in terms of the \{1,3\}-inverse, Drazin inverse and pseudo core inverse, respectively. Then, we give the pseudo core decomposition utilizing the pseudo core inverse, which extends the core-EP decomposition introduced by Wang for matrices to an arbitrary $*$-ring; and this decomposition turns to be a useful tool to characterize *-DMP elements. Further, we extend Wang's core-EP order from matrices to $*$-rings and use it to investigate *-DMP elements. Finally, we give necessary and sufficient conditions for two elements $a,~b$ in $*$-rings to have $aa^{\scriptsize\textcircled{\tiny D}}=bb^{\scriptsize\textcircled{\tiny D}}$, which contribute to investigate *-DMP elements.


title: Knowledge Engineering for Hybrid Deductive Databases

abstract: Modern knowledge base systems frequently need to combine a collection of databases in different formats: e.g., relational databases, XML databases, rule bases, ontologies, etc. In the deductive database system DDBASE, we can manage these different formats of knowledge and reason about them. Even the file systems on different computers can be part of the knowledge base. Often, it is necessary to handle different versions of a knowledge base. E.g., we might want to find out common parts or differences of two versions of a relational database. We will examine the use of abstractions of rule bases by predicate dependency and rule predicate graphs. Also the proof trees of derived atoms can help to compare different versions of a rule base. Moreover, it might be possible to have derivations joining rules with other formalisms of knowledge representation. Ontologies have shown their benefits in many applications of intelligent systems, and there have been many proposals for rule languages compatible with the semantic web stack, e.g., SWRL, the semantic web rule language. Recently, ontologies are used in hybrid systems for specifying the provenance of the different components.


title: Bottom-Up Evaluation of Datalog: Preliminary Report

abstract: Bottom-up evaluation of Datalog has been studied for a long time, and is standard material in textbooks. However, if one actually wants to develop a deductive database system, it turns out that there are many implementation options. For instance, the sequence in which rule instances are applied is not given. In this paper, we study a method that immediately uses a derived tuple to derive more tuples (called the Push method). In this way, storage space for intermediate results can be reduced. The main contribution of our method is the way in which we minimize the copying of values at runtime, and do much work already at compile-time.


title: A Practical View on Renaming

abstract: We revisit variable renaming from a practitioner's point of view, presenting concepts we found useful in dealing with operational semantics of pure Prolog. A concept of relaxed core representation is introduced, upon which a concept of prenaming is built. Prenaming formalizes the intuitive practice of renaming terms by just considering the necessary bindings, where now some passive "bindings" x/x may be necessary as well. As an application, a constructive version of variant lemma for implemented Horn clause logic has been obtained. There, prenamings made it possible to incrementally handle new (local) variables.


title: Limit theorems for Hilbert space-valued linear processes under long range dependence

abstract: Let $(X_{k})_{k \in \mathbb Z }$ be a linear process with values in a separable Hilbert space $\mathbb{H}$ given by $X_{k} =\sum_{j=0}^{\infty} (j+1)^{-N}\varepsilon_{k-j}$ for each $k \in \mathbb Z$, where $N:\mathbb{H} \to \mathbb{H}$ is a bounded, linear normal operator and $(\varepsilon_{k})_{ k \in \mathbb Z }$ is a sequence of independent, identically distributed $\mathbb{H}$-valued random variables with $E\varepsilon_{0}=0$ and $E\| \varepsilon_{0} \|^2<\infty$. We investigate the central and the functional central limit theorem for $(X_{k})_{k \in \mathbb Z }$ when the series of operator norms $\sum_{j=0}^{\infty} \|(j+1)^{-N}\|_{op}$ diverges. Furthermore we show that the limit process in case of the functional central limit theorem generates an operator self-similar process.


title: Implementing GraphQL as a Query Language for Deductive Databases in SWI-Prolog Using DCGs, Quasi Quotations, and Dicts

abstract: The methods to access large relational databases in a distributed system are well established: the relational query language SQL often serves as a language for data access and manipulation, and in addition public interfaces are exposed using communication protocols like REST. Similarly to REST, GraphQL is the query protocol of an application layer developed by Facebook. It provides a unified interface between the client and the server for data fetching and manipulation. Using GraphQL's type system, it is possible to specify data handling of various sources and to combine, e.g., relational with NoSQL databases. In contrast to REST, GraphQL provides a single API endpoint and supports flexible queries over linked data. GraphQL can also be used as an interface for deductive databases. In this paper, we give an introduction of GraphQL and a comparison to REST. Using language features recently added to SWI-Prolog 7, we have developed the Prolog library GraphQL.pl, which implements the GraphQL type system and query syntax as a domain-specific language with the help of definite clause grammars (DCG), quasi quotations, and dicts. Using our library, the type system created for a deductive database can be validated, while the query system provides a unified interface for data access and introspection.


title: Experiences with Some Benchmarks for Deductive Databases and Implementations of Bottom-Up Evaluation

abstract: OpenRuleBench is a large benchmark suite for rule engines, which includes deductive databases. We previously proposed a translation of Datalog to C++ based on a method that "pushes" derived tuples immediately to places where they are used. In this paper, we report performance results of various implementation variants of this method compared to XSB, YAP and DLV. We study only a fraction of the OpenRuleBench problems, but we give a quite detailed analysis of each such task and the factors which influence performance. The results not only show the potential of our method and implementation approach, but could be valuable for anybody implementing systems which should be able to execute tasks of the discussed types.


title: Immortal homogeneous Ricci flows

abstract: We show that for an immortal homogeneous Ricci flow solution any sequence of parabolic blow-downs subconverges to a homogeneous expanding Ricci soliton. This is established by constructing a new Lyapunov function based on curvature estimates which come from real geometric invariant theory.


title: Constraint Logic Programming over Infinite Domains with an Application to Proof

abstract: We present a CLP(FD)-based constraint solver able to deal with unbounded domains. It is based on constraint propagation, resorting to enumeration if all other methods fail. An important aspect is detecting when enumeration was complete and if this has an impact on the soundness of the result. We present a technique which guarantees soundness in the following way: if the constraint solver finds a solution it is guaranteed to be correct; if the constraint solver fails to find a solution it can either return the result "definitely false" in case it knows enumeration was exhaustive, or "unknown" in case it was aborted. The technique can deal with nested universal and existential quantifiers. It can easily be extended to set comprehensions and other operators introducing new quantified variables. We show applications in data validation and proof.


title: Source Code Verification for Embedded Systems using Prolog

abstract: System relevant embedded software needs to be reliable and, therefore, well tested, especially for aerospace systems. A common technique to verify programs is the analysis of their abstract syntax tree (AST). Tree structures can be elegantly analyzed with the logic programming language Prolog. Moreover, Prolog offers further advantages for a thorough analysis: On the one hand, it natively provides versatile options to efficiently process tree or graph data structures. On the other hand, Prolog's non-determinism and backtracking eases tests of different variations of the program flow without big effort. A rule-based approach with Prolog allows to characterize the verification goals in a concise and declarative way. In this paper, we describe our approach to verify the source code of a flash file system with the help of Prolog. The flash file system is written in C++ and has been developed particularly for the use in satellites. We transform a given abstract syntax tree of C++ source code into Prolog facts and derive the call graph and the execution sequence (tree), which then are further tested against verification goals. The different program flow branching due to control structures is derived by backtracking as subtrees of the full execution sequence. Finally, these subtrees are verified in Prolog. We illustrate our approach with a case study, where we search for incorrect applications of semaphores in embedded software using the real-time operating system RODOS. We rely on computation tree logic (CTL) and have designed an embedded domain specific language (DSL) in Prolog to express the verification goals.


title: A Typeful Integration of SQL into Curry

abstract: We present an extension of the declarative programming language Curry to support the access to data stored in relational databases via SQL. Since Curry is statically typed, our emphasis on this SQL integration is on type safety. Our extension respects the type system of Curry so that run-time errors due to ill-typed data are avoided. This is obtained by preprocessing SQL statements at compile time and translating them into type-safe database access operations. As a consequence, the type checker of the Curry system can spot type errors in SQL statements at compile time. To generate appropriately typed access operations, the preprocessor uses an entity-relationship (ER) model describing the structure of the relational data. In addition to standard SQL, SQL statements embedded in Curry can include program expressions and also relationships specified in the ER model. The latter feature is useful to avoid the error-prone use of foreign keys. As a result, our SQL integration supports a high-level and type-safe access to databases in Curry programs.


title: A Simulation Tool for tccp Programs

abstract: The Timed Concurrent Constraint Language tccp is a declarative synchronous concurrent language, particularly suitable for modelling reactive systems. In tccp, agents communicate and synchronise through a global constraint store. It supports a notion of discrete time that allows all non-blocked agents to proceed with their execution simultaneously. In this paper, we present a modular architecture for the simulation of tccp programs. The tool comprises three main components. First, a set of basic abstract instructions able to model the tccp agent behaviour, the memory model needed to manage the active agents and the state of the store during the execution. Second, the agent interpreter that executes the instructions of the current agent iteratively and calculates the new agents to be executed at the next time instant. Finally, the constraint solver components which are the modules that deal with constraints. In this paper, we describe the implementation of these components and present an example of a real system modelled in tccp.


title: A Framework for Extending microKanren with Constraints

abstract: We present a framework for building CLP languages with symbolic constraints based on microKanren, a domain-specific logic language shallowly embedded in Racket. We rely on Racket's macro system to generate a constraint solver and other components of the microKanren embedding. The framework itself and the constraints' implementations amounts to just over 100 lines of code. Our framework is both a teachable implementation for CLP as well as a test-bed and prototyping tool for symbolic constraint systems.


title: A Practical Study of Control in Objected-Oriented--Functional--Logic Programming with Paisley

abstract: Paisley is an extensible lightweight embedded domain-specific language for nondeterministic pattern matching in Java. Using simple APIs and programming idioms, it brings the power of functional-logic processing of arbitrary data objects to the Java platform, without constraining the underlying object-oriented semantics. Here we present an extension to the Paisley framework that adds pattern-based control flow. It exploits recent additions to the Java language, namely functional interfaces and lambda expressions, for an explicit and transparent continuation-passing style approach to control. We evaluate the practical impact of the novel features on a real-world case study that reengineers a third-party open-source project to use Paisley in place of conventional object-oriented data query idioms. We find the approach viable for incremental refactoring of legacy code, with significant qualitative improvements regarding separation of concerns, clarity and intentionality, thus making for easier code understanding, testing and debugging.


title: An Agglomeration Law for Sorting Networks and its Application in Functional Programming

abstract: In this paper we will present a general agglomeration law for sorting networks. Agglomeration is a common technique when designing parallel programmes to control the granularity of the computation thereby finding a better fit between the algorithm and the machine on which the algorithm runs. Usually this is done by grouping smaller tasks and computing them en bloc within one parallel process. In the case of sorting networks this could be done by computing bigger parts of the network with one process. The agglomeration law in this paper pursues a different strategy: The input data is grouped and the algorithm is generalized to work on the agglomerated input while the original structure of the algorithm remains. This will result in a new access opportunity to sorting networks well-suited for efficient parallelization on modern multicore computers, computer networks or GPGPU programming. Additionally this enables us to use sorting networks as (parallel or distributed) merging stages for arbitrary sorting algorithms, thereby creating new hybrid sorting algorithms with ease. The expressiveness of functional programming languages helps us to apply this law to systematically constructed sorting networks, leading to efficient and easily adaptable sorting algorithms. An application example is given, using the Eden programming language to show the effectiveness of the law. The implementation is compared with different parallel sorting algorithms by runtime behaviour.


title: Proving Non-Deterministic Computations in Agda

abstract: We investigate proving properties of Curry programs using Agda. First, we address the functional correctness of Curry functions that, apart from some syntactic and semantic differences, are in the intersection of the two languages. Second, we use Agda to model non-deterministic functions with two distinct and competitive approaches incorporating the non-determinism. The first approach eliminates non-determinism by considering the set of all non-deterministic values produced by an application. The second approach encodes every non-deterministic choice that the application could perform. We consider our initial experiment a success. Although proving properties of programs is a notoriously difficult task, the functional logic paradigm does not seem to add any significant layer of difficulty or complexity to the task.


title: On Upper Bounds on the Church-Rosser Theorem

abstract: The Church-Rosser theorem in the type-free lambda-calculus is well investigated both for beta-equality and beta-reduction. We provide a new proof of the theorem for beta-equality with no use of parallel reductions, but simply with Takahashi's translation (Gross-Knuth strategy). Based on this, upper bounds for reduction sequences on the theorem are obtained as the fourth level of the Grzegorczyk hierarchy.


title: Confluence of Conditional Term Rewrite Systems via Transformations

abstract: Conditional term rewriting is an intuitive yet complex extension of term rewriting. In order to benefit from the simpler framework of unconditional rewriting, transformations have been defined to eliminate the conditions of conditional term rewrite systems. Recent results provide confluence criteria for conditional term rewrite systems via transformations, yet they are restricted to CTRSs with certain syntactic properties like weak left-linearity. These syntactic properties imply that the transformations are sound for the given CTRS. This paper shows how to use transformations to prove confluence of operationally terminating, right-stable deterministic conditional term rewrite systems without the necessity of soundness restrictions. For this purpose, it is shown that certain rewrite strategies, in particular almost U-eagerness and innermost rewriting, always imply soundness.


title: An Extension of Proof Graphs for Disjunctive Parameterised Boolean Equation Systems

abstract: A parameterised Boolean equation system (PBES) is a set of equations that defines sets as the least and/or greatest fixed-points that satisfy the equations. This system is regarded as a declarative program defining functions that take a datum and returns a Boolean value. The membership problem of PBESs is a problem to decide whether a given element is in the defined set or not, which corresponds to an execution of the program. This paper introduces reduced proof graphs, and studies a technique to solve the membership problem of PBESs, which is undecidable in general, by transforming it into a reduced proof graph. A vertex X(v) in a proof graph represents that the data v is in the set X, if the graph satisfies conditions induced from a given PBES. Proof graphs are, however, infinite in general. Thus we introduce vertices each of which stands for a set of vertices of the original ones, which possibly results in a finite graph. For a subclass of disjunctive PBESs, we clarify some conditions which reduced proof graphs should satisfy. We also show some examples having no finite proof graph except for reduced one. We further propose a reduced dependency space, which contains reduced proof graphs as sub-graphs if a proof graph exists. We provide a procedure to construct finite reduced dependency spaces, and show the soundness and completeness of the procedure.


title: An Environment for Analyzing Space Optimizations in Call-by-Need Functional Languages

abstract: We present an implementation of an interpreter LRPi for the call-by-need calculus LRP, based on a variant of Sestoft's abstract machine Mark 1, extended with an eager garbage collector. It is used as a tool for exact space usage analyses as a support for our investigations into space improvements of call-by-need calculi.


title: 3d supersymmetric gauge theories and Hilbert series

abstract: The Hilbert series is a generating function that enumerates gauge invariant chiral operators of a supersymmetric field theory with four supercharges and an R-symmetry. In this article I review how counting dressed 't Hooft monopole operators leads to a formula for the Hilbert series of a 3d $\mathcal{N}\geq 2$ gauge theory, which captures precious information about the chiral ring and the moduli space of supersymmetric vacua of the theory.


title: Finding Risk-Averse Shortest Path with Time-dependent Stochastic Costs

abstract: In this paper, we tackle the problem of risk-averse route planning in a transportation network with time-dependent and stochastic costs. To solve this problem, we propose an adaptation of the A* algorithm that accommodates any risk measure or decision criterion that is monotonic with first-order stochastic dominance. We also present a case study of our algorithm on the Manhattan, NYC, transportation network.