title: Extracting the orbital axis from gravitational waves of precessing binary systems

abstract: We present a new method for extracting the instantaneous orbital axis only from gravitational wave strains of precessing binary systems observed from a particular observer direction. This method enables us to reconstruct the co-precessing frame waveforms only from observed quantities for the ideal case that the signal-to-noise ratio is high enough to analyze the waveforms directly. Specifically, we do not assume knowledge of the time evolution of the instantaneous orbital axis and the co-precessing waveforms before analyzing the data in our method. We test and measure the accuracy of our method using the numerical relativity simulation data of precessing binary black holes taken from the SXS Catalog. We show that the direction of the orbital axis is extracted within $\approx0.02~{\rm rad}$ error from gravitational waves emitted during the inspiral phase. The co-precessing waveforms are also reconstructed with high accuracy; the mismatch (assuming white noise) between them and the original co-precessing waveforms is typically a few times $10^{-3}$ including the merger-ringdown phase, and can be improved by an order of magnitude focusing only on the inspiral waveform. In this method, the co-precessing frame waveforms are not only the purely technical tools for understanding the complex nature of precessing waveforms but also direct observables.


title: Experience enrichment based task independent reward model

abstract: For most reinforcement learning approaches, the learning is performed by maximizing an accumulative reward that is expectedly and manually defined for specific tasks. However, in real world, rewards are emergent phenomena from the complex interactions between agents and environments. In this paper, we propose an implicit generic reward model for reinforcement learning. Unlike those rewards that are manually defined for specific tasks, such implicit reward is task independent. It only comes from the deviation from the agents' previous experiences.


title: Shallow Updates for Deep Reinforcement Learning

abstract: Deep reinforcement learning (DRL) methods such as the Deep Q-Network (DQN) have achieved state-of-the-art results in a variety of challenging, high-dimensional domains. This success is mainly attributed to the power of deep neural networks to learn rich domain representations for approximating the value function or policy. Batch reinforcement learning methods with linear representations, on the other hand, are more stable and require less hyper parameter tuning. Yet, substantial feature engineering is necessary to achieve good results. In this work we propose a hybrid approach -- the Least Squares Deep Q-Network (LS-DQN), which combines rich feature representations learned by a DRL algorithm with the stability of a linear least squares method. We do this by periodically re-training the last hidden layer of a DRL network with a batch least squares update. Key to our approach is a Bayesian regularization term for the least squares update, which prevents over-fitting to the more recent data. We tested LS-DQN on five Atari games and demonstrate significant improvement over vanilla DQN and Double-DQN. We also investigated the reasons for the superior performance of our method. Interestingly, we found that the performance improvement can be attributed to the large batch size used by the LS method when optimizing the last layer.


title: The Gelfand--Shilov type estimate for Green's function of the bounded solutions problem

abstract: An analog of the Gelfand--Shilov estimate of the matrix exponential is proved for Green's function of the problem of bounded solutions of the ordinary differential equation $x'(t)-Ax(t)=f(t)$.


title: Spatially Controlled Relay Beamforming: $2$-Stage Optimal Policies

abstract: The problem of enhancing Quality-of-Service (QoS) in power constrained, mobile relay beamforming networks, by optimally and dynamically controlling the motion of the relaying nodes, is considered, in a dynamic channel environment. We assume a time slotted system, where the relays update their positions before the beginning of each time slot. Modeling the wireless channel as a Gaussian spatiotemporal stochastic field, we propose a novel $2$-stage stochastic programming problem formulation for optimally specifying the positions of the relays at each time slot, such that the expected QoS of the network is maximized, based on causal Channel State Information (CSI) and under a total relay transmit power budget. This results in a schema where, at each time slot, the relays, apart from optimally beamforming to the destination, also optimally, predictively decide their positions at the next time slot, based on causally accumulated experience. Exploiting either the Method of Statistical Differentials, or the multidimensional Gauss-Hermite Quadrature Rule, the stochastic program considered is shown to be approximately equivalent to a set of simple subproblems, which are solved in a distributed fashion, one at each relay. Optimality and performance of the proposed spatially controlled system are also effectively assessed, under a rigorous technical framework; strict optimality is rigorously demonstrated via the development of a version of the Fundamental Lemma of Stochastic Control, and, performance-wise, it is shown that, quite interestingly, the optimal average network QoS exhibits an increasing trend across time slots, despite our myopic problem formulation. Numerical simulations are presented, experimentally corroborating the success of the proposed approach and the validity of our theoretical predictions.

categories: math.OC cs.IT math.IT math.ST stat.AP stat.ME stat.TH


title: Future Experimental Improvement for the Search of LNV Process in $e\mu$ Sector

abstract: Exploring the leptonic sector in frontier experiments is more of importance nowadays, since the conservation of lepton flavor and total lepton number are not guaranteed anymore in the Standard Model after the discovery of neutrino oscillations. $\mu^- + N(A,Z) \rightarrow e^+ + N(A,Z-2)$ conversion in a muonic atom is one of the most promising channels to investigate the lepton number violation process, and the measurement of this process is planned in future $\mu^--e^-$ conversion experiments with a muonic atom in a muon-stopping target. This paper discusses how to maximize the experimental sensitivity of the $\mu^--e^+$ conversion by introducing the new requirement of the mass relation of $M(A,Z-2)<M(A,Z-1)$, where $M(A,Z)$ is the mass of the muon-stopping target nucleus, to get rid of the background from radiative muon capture. The sensitivity of the $\mu^--e^+$ conversion is anticipated to have four orders of magnitude of improvement in forthcoming experiments using a proper target nucleus, which satisfies the mass relation. The most promising isotopes found are $^{40}$Ca and $^{32}$S.


title: Some Schemes for Implementation of Arithmetic Operations with Complex Numbers Using Squaring Units

abstract: In this paper, new schemes for a squarer, multiplier and divider of complex numbers are proposed. Traditional structural solutions for each of these operations require the presence some number of general-purpose binary multipliers. The advantage of our solutions is a removing of multiplications through replacing them by less costly squarers. We use Logan's trick and quarter square technique, which propose to replace the calculation of the product of two real numbers by summing the squares. Replacing usual multipliers on digital squares implies reducing power consumption as well as decreases hardware circuit complexity. The squarer requiring less area and power as compared to general-purpose multiplier, it is interesting to assess the use of squarers to implementation of complex arithmetic.


title: Reconstruction formulas for Photoacoustic Imaging in Attenuating Media

abstract: In this paper we study the problem of photoacoustic inversion in a weakly attenuating medium. We present explicit reconstruction formulas in such media and show that the inversion based on such formulas is moderately ill--posed. Moreover, we present a numerical algorithm for imaging and demonstrate in numerical experiments the feasibility of this approach.


title: Subcritical growth of electron phase-space holes in planetary radiation belts

abstract: The discovery of long-lived electrostatic coherent structures with large-amplitude electric fields ($1 \leq E \leq 500 $ mV/m) by the Van Allen Probes has revealed alternative routes through which planetary radiation belts' acceleration can take place. Following previous reports showing that small phase-space holes, with $q\phi /T^c_e\simeq 10^{-2}-10^{-3}$, could result from electron interaction with large-amplitude whistlers, we demonstrate one possible mechanism through which holes can grow nonlinearly (i.e. $\gamma \propto \sqrt{\phi}$) and subcritically as a result of momentum exchange between hot and cold electron populations. Our results provide an explanation for the common occurrence and fast growth of large-amplitude electron phase-space holes in the Earth's radiation belts.


title: Including surface ligand effects in continuum elastic models of nanocrystal vibrations

abstract: The measured low frequency vibrational energies of some quantum dots (QDs) deviate from the predictions of traditional elastic continuum models. Recent experiments have revealed that these deviations can be tuned by changing the ligands that passivate the QD surface. This observation has led to speculation that these deviations are due to a mass-loading effect of the surface ligands. In this article, we address this speculation by formulating a continuum elastic theory that includes the mass-loading effects of the surface ligands. We demonstrate that this model is capable of accurately reproducing the $l = 0$ phonon energy across variety of different QD samples, including cores with different ligand identities and epitaxially grown CdSe/CdS core/shell heterostructures. We highlight that our model performs well even in the small QD regime, where traditional elastic continuum models are especially prone to failure. Furthermore, we show that our model combined with Raman measurements can be used to infer the elastic properties of surface bound ligands, such as sound velocities and elastic moduli, that are otherwise challenging.


title: Improved Algorithms for Matrix Recovery from Rank-One Projections

abstract: We consider the problem of estimation of a low-rank matrix from a limited number of noisy rank-one projections. In particular, we propose two fast, non-convex \emph{proper} algorithms for matrix recovery and support them with rigorous theoretical analysis. We show that the proposed algorithms enjoy linear convergence and that their sample complexity is independent of the condition number of the unknown true low-rank matrix. By leveraging recent advances in low-rank matrix approximation techniques, we show that our algorithms achieve computational speed-ups over existing methods. Finally, we complement our theory with some numerical experiments.


title: Commensurability invariance for abelian splittings of right-angled Artin groups, braid groups and loop braid groups

abstract: We prove that if a right-angled Artin group $A_\Gamma$ is abstractly commensurable to a group splitting non-trivially as an amalgam or HNN-extension over $\mathbb{Z}^n$, then $A_\Gamma$ must itself split non-trivially over $\mathbb{Z}^k$ for some $k\le n$. Consequently, if two right-angled Artin groups $A_\Gamma$ and $A_\Delta$ are commensurable and $\Gamma$ has no separating $k$-cliques for any $k\le n$ then neither does $\Delta$, so "smallest size of separating clique" is a commensurability invariant. We also discuss some implications for issues of quasi-isometry. Using similar methods we also prove that for $n\ge 4$ the braid group $B_n$ is not abstractly commensurable to any group that splits non-trivially over a "free group-free" subgroup, and the same holds for $n\ge 3$ for the loop braid group $LB_n$. Our approach makes heavy use of the Bieri--Neumann--Strebel invariant.


title: Spatial Distribution of Element Abundances and Ionization States in Solar Energetic-Particle Events

abstract: We have studied the spatial and temporal distribution of abundances of chemical elements in large "gradual" solar energetic-particle (SEP) events, and especially the source plasma temperatures, derived from those abundances, using measurements from the Wind and Solar TErestrial RElations Observatory (STEREO) spacecraft, widely separated in solar longitude. A power-law relationship between abundance enhancements and mass-to-charge ratios [A/Q] of the ions can be used to determine Q-values and source plasma temperatures at remote spacecraft with instruments that were not designed for charge-state measurements. We search for possible source variations along the accelerating shock wave, finding one clear case where the accelerating shock wave appears to dispatch ions from 3.2+-0.8 MK plasma toward one spacecraft and those from 1.6+-0.2 MK plasma toward another, 116 deg away. The difference persists three days and then fades away. Three other SEP events show less-extreme variation in source temperatures at different spacecraft, in one case observed over 222 deg in longitude. This initial study shows how the power-law relation between abundance enhancements and ion A/Q-values pro-vides a new technique to determine Q and plasma temperatures in the seed population of SEP ions over a broad region of space using remote spacecraft with instruments that were not originally designed for measurements of ionization states.


title: On the Black's equation for the risk tolerance function

abstract: We analyze a nonlinear equation proposed by F. Black (1968) for the optimal portfolio function in a log-normal model. We cast it in terms of the risk tolerance function and provide, for general utility functions, existence, uniqueness and regularity results, and we also examine various monotonicity, concavity/convexity and S-shape properties. Stronger results are derived for utilities whose inverse marginal belongs to a class of completely monotonic functions.


title: Nonautonomous Young differential equations revisited

abstract: In this paper we prove that under weak conditions a nonautonomous Young differential equation possesses a unique solution which depends continuously on initial conditions. The proofs use estimates in p-variation norms, greedy time techniques, and Gronwall-type lemma with the help of Shauder theorem of fixed points.


title: Why are Big Data Matrices Approximately Low Rank?

abstract: Matrices of (approximate) low rank are pervasive in data science, appearing in recommender systems, movie preferences, topic models, medical records, and genomics. While there is a vast literature on how to exploit low rank structure in these datasets, there is less attention on explaining why the low rank structure appears in the first place. Here, we explain the effectiveness of low rank models in data science by considering a simple generative model for these matrices: we suppose that each row or column is associated to a (possibly high dimensional) bounded latent variable, and entries of the matrix are generated by applying a piecewise analytic function to these latent variables. These matrices are in general full rank. However, we show that we can approximate every entry of an $m \times n$ matrix drawn from this model to within a fixed absolute error by a low rank matrix whose rank grows as $\mathcal O(\log(m + n))$. Hence any sufficiently large matrix from such a latent variable model can be approximated, up to a small entrywise error, by a low rank matrix.


title: A Strongly Interacting Polaritonic Quantum Dot

abstract: Polaritons are an emerging platform for exploration of synthetic materials [1] and quantum information processing [2] that draw properties from two disparate particles: a photon and an atom. Cavity polaritons are particularly promising, as they are long-lived and their dispersion and mass are controllable through cavity geometry [3]. To date, studies of cavity polaritons have operated in the mean-field regime, using short-range interactions between their matter components [4]. Rydberg excitations have recently been demonstrated as a promising matter-component of polaritons [5], due to their strong interactions over distances large compared to an optical wavelength. In this work we explore, for the first time, the cavity quantum electrodynamics of Rydberg polaritons, combining the non-linearity of polaritonic quantum wires with the zero-dimensional strong coupling of an optical resonator. We assemble a quantum dot composed of $\sim 150$ strongly interacting, Rydberg-dressed $^{87}$Rb atoms in a cavity, and observe blockaded polariton transport as well as coherent quantum dynamics of a single polaritonic super-atom. This work establishes a new generation of photonic quantum information processors and quantum materials, along with a clear path to topological quantum matter [6].


title: On Stackelberg Mixed Strategies

abstract: It is sometimes the case that one solution concept in game theory is equivalent to applying another solution concept to a modified version of the game. In such cases, does it make sense to study the former separately (as it applies to the original representation of the game), or should we entirely subordinate it to the latter? The answer probably depends on the particular circumstances, and indeed the literature takes different approaches in different cases. In this article, I consider the specific example of Stackelberg mixed strategies. I argue that, even though a Stackelberg mixed strategy can also be seen as a subgame perfect Nash equilibrium of a corresponding extensive-form game, there remains significant value in studying it separately. The analysis of this special case may have implications for other solution concepts.


title: Statistical inference using SGD

abstract: We present a novel method for frequentist statistical inference in $M$-estimation problems, based on stochastic gradient descent (SGD) with a fixed step size: we demonstrate that the average of such SGD sequences can be used for statistical inference, after proper scaling. An intuitive analysis using the Ornstein-Uhlenbeck process suggests that such averages are asymptotically normal. From a practical perspective, our SGD-based inference procedure is a first order method, and is well-suited for large scale problems. To show its merits, we apply it to both synthetic and real datasets, and demonstrate that its accuracy is comparable to classical statistical methods, while requiring potentially far less computation.

categories: cs.LG cs.AI math.OC math.ST stat.ML stat.TH


title: Report of the HPC Correctness Summit, Jan 25--26, 2017, Washington, DC

abstract: Maintaining leadership in HPC requires the ability to support simulations at large scales and fidelity. In this study, we detail one of the most significant productivity challenges in achieving this goal, namely the increasing proclivity to bugs, especially in the face of growing hardware and software heterogeneity and sheer system scale. We identify key areas where timely new research must be proactively begun to address these challenges, and create new correctness tools that must ideally play a significant role even while ramping up toward exacale. We close with the proposal for a two-day workshop in which the problems identified in this report can be more broadly discussed, and specific plans to launch these new research thrusts identified.


title: Nonequilibrium dynamics of mixtures of active and passive colloidal particles

abstract: We develop a mesoscopic field theory for the collective nonequilibrium dynamics of multicomponent mixtures of interacting active (i.e., motile) and passive (i.e., nonmotile) colloidal particles with isometric shape in two spatial dimensions. By a stability analysis of the field theory, we obtain equations for the spinodal that describes the onset of a motility-induced instability leading to cluster formation in such mixtures. The prediction for the spinodal is found to be in good agreement with particle-resolved computer simulations. Furthermore, we show that in active-passive mixtures the spinodal instability can be of two different types. One type is associated with a stationary bifurcation and occurs also in one-component active systems, whereas the other type is associated with a Hopf bifurcation and can occur only in active-passive mixtures. Remarkably, the Hopf bifurcation leads to moving clusters. This explains recent results from simulations of active-passive particle mixtures, where moving clusters and interfaces that are not seen in the corresponding one-component systems have been observed.


title: A Closer Look at $R_D$ and $R_{D^\ast}$

abstract: The measurement of $R_D$ ($R_{D^*}$), the ratio of the branching fraction of $\overline{B} \to D \tau \bar{\nu}_\tau (\overline{B} \to D^* \tau \bar{\nu}_\tau)$ to that of $\overline{B} \to D l \bar{\nu}_l (\overline{B} \to D^* l \bar{\nu}_l)$, shows $1.9 \sigma$ $(3.3 \sigma)$ deviation from its Standard Model (SM) prediction. The combined deviation is at the level of $4 \sigma$ according to the Heavy Flavour Averaging Group (HFAG). We perform an effective field theory analysis (at the dimension 6 level) of these potential New Physics (NP) signals assuming $ SU(3)_{C} \times SU(2)_{L} \times U(1)_{Y}$ gauge invariance. We first show that, in general, $R_D$ and $R_{D^*}$ are theoretically independent observables and hence, their theoretical predictions are not correlated. We identify the operators that can explain the experimental measurements of $R_D$ and $R_{D^*}$ individually and also together. Motivated by the recent measurement of the $\tau$ polarisation in $\overline{B} \to D^* \tau \bar{\nu}_\tau$ decay, $P_\tau^{D^*}$ by the {\sc Belle} collaboration, we study the impact of a more precise measurement of $P_\tau^{D^*}$ (and a measurement of $P_\tau^D$) on the various possible NP explanations. Furthermore, we show that the measurement of $R_{D^*}$ in bins of $q^2$, the square of the invariant mass of the lepton neutrino system, along with the information on $\tau$ polarisation, can completely distinguish the various operator structures.


title: Spin light of neutrino in astrophysical environments

abstract: The ${\it {spin \ \ light \ \ of \ \ neutrino}}$ ($SL\nu$) is a new possible mechanism of electromagnetic radiation by a massive neutrino (with a nonzero magnetic moment) moving in media. Since the prediction of this mechanism, the question has been debated in a number of publications as whether the effect can be of any significance for realistic astrophysical conditions. Although this effect is strongly suppressed due to smallness of neutrino magnetic moment, for ultra-high energy neutrinos (PeV neutrinos recently observed by the IceCube collaboration, for instance) the $SL\nu$ might be of interest in the case of neutrinos propagating in dense matter. An advanced view on the $SL\nu$ in matter is given, and several astrophysical settings (a neutron star, supernova, Gamma-Ray Burst (GRB), and relic neutrino background) for which the effect can be realized are considered. Taking into account the threshold condition and also several competing processes, we determine conditions for which the $SL\nu$ mechanism is possible. We conclude that the most favorable case of the effect manifestation is provided by ultra dense matter of neutron stars and ultrahigh energy of the radiating neutrino, and note that these conditions can be met within galaxy clusters. It is also shown that due to the $SL\nu$ specific polarization properties this electromagnetic mechanism is of interest in the connection with the observed polarization of GRB emission.


title: Sharp geometric inequalities for the general $p$-affine capacity

abstract: In this article, we propose the notion of the general $p$-affine capacity and prove some basic properties for the general $p$-affine capacity, such as affine invariance and monotonicity. The newly proposed general $p$-affine capacity is compared with several classical geometric quantities, e.g., the volume, the $p$-variational capacity and the $p$-integral affine surface area. Consequently, several sharp geometric inequalities for the general $p$-affine capacity are obtained. These inequalities extend and strengthen many well-known (affine) isoperimetric and (affine) isocapacitary inequalities.


title: TuRF: Fast Data Collection for Fingerprint-based Indoor Localization

abstract: Many infrastructure-free indoor positioning systems rely on fine-grained location-dependent fingerprints to train models for localization. The site survey process to collect fingerprints is laborious and is considered one of the major obstacles to deploying such systems. In this paper, we propose TuRF, a fast path-based fingerprint collection mechanism for site survey. We demonstrate the feasibility to collect fingerprints for indoor localization during walking along predefined paths. A step counter is utilized to accommodate the variations in walking speed. Approximate location labels inferred from the steps are then used to train a Gaussian Process regression model. Extensive experiments show that TuRF can significantly reduce the required time for site survey, without compromising the localization performance.


title: On an equation involving fractional powers with prime numbers of a special type

abstract: We consider the equation $[p_{1}^{c}] + [p_{2}^{c}] + [p_{3}^{c}] = N$, where $N$ is a sufficiently large integer, and prove that if $1 < c < \frac{17}{16}$, then it has a solution in prime numbers $p_{1}$, $p_{2}$, $p_{3}$ such that each of the numbers $p_{1} + 2$, $p_{2} + 2$, $p_{3} + 2$ has at most $\left [ \frac{95}{17 - 16c} \right ]$ prime factors, counted with the multiplicity.


title: Shake-Shake regularization

abstract: The method introduced in this paper aims at helping deep learning practitioners faced with an overfit problem. The idea is to replace, in a multi-branch network, the standard summation of parallel branches with a stochastic affine combination. Applied to 3-branch residual networks, shake-shake regularization improves on the best single shot published results on CIFAR-10 and CIFAR-100 by reaching test errors of 2.86% and 15.85%. Experiments on architectures without skip connections or Batch Normalization show encouraging results and open the door to a large set of applications. Code is available at https://github.com/xgastaldi/shake-shake


title: A system for online beam emittance measurements and proton beam characterization

abstract: A system for online measurement of the transverse beam emittance was developed. It is named $^{4}$PrOB$\varepsilon$aM (4-Profiler Online Beam Emittance Measurement) and was conceived to measure the emittance in a fast and efficient way using the multiple beam profiler method. The core of the system is constituted by four consecutive UniBEaM profilers, which are based on silica fibers passing across the beam. The $^{4}$PrOB$\varepsilon$aM system was deployed for characterization studies of the 18~MeV proton beam produced by the IBA Cyclone 18 MeV cyclotron at Bern University Hospital (Inselspital). The machine serves daily radioisotope production and multi-disciplinary research, which is carried out with a specifically conceived Beam Transport Line (BTL). The transverse RMS beam emittance of the cyclotron was measured as a function of several machine parameters, such as the magnetic field, RF peak voltage, and azimuthal angle of the stripper. The beam emittance was also measured using the method based on the quadrupole strength variation. The results obtained with both techniques were compared and a good agreement was found. In order to characterize the longitudinal dynamics, the proton energy distribution was measured. For this purpose, a method was developed based on aluminum absorbers of different thicknesses, a UniBEaM detector, and a Faraday cup. The results were an input for a simulation of the BTL developed in the MAD-X software. This tool allows machine parameters to be tuned online and the beam characteristics to be optimized for specific applications.


title: A modified Newton iteration for finding nonnegative Z-eigenpairs of a nonnegative tensor

abstract: We propose a modified Newton iteration for finding some nonnegative Z-eigenpairs of a nonnegative tensor. When the tensor is irreducible, all nonnegative eigenpairs are known to be positive. We prove local quadratic convergence of the new iteration to any positive eigenpair of a nonnegative tensor, under the usual assumption guaranteeing the local quadratic convergence of the original Newton iteration. A big advantage of the modified Newton iteration is that it seems capable of finding a nonnegative eigenpair starting with any positive unit vector. Special attention is paid to transition probability tensors.


title: On cohomological Hall algebras of quivers : generators

abstract: We study the cohomological Hall algebra Y of a lagrangian substack of the moduli stack of representations of the preprojective algebra of an arbitrary quiver Q, and their actions on the cohomology of Nakajima quiver varieties. We prove that Y is pure and we compute its Poincare polynomials in terms of (nilpotent) Kac polynomials. We also provide a family of algebra generators. We conjecture that Y is equal, after a suitable extension of scalars, to the Yangian introduced by Maulik and Okounkov. As a corollary, we prove a variant of Okounkov's conjecture, which is a generalization of the Kac conjecture relating the constant term of Kac polynomials to root multiplicities of Kac-Moody algebras.


title: Dynamical depinning of chiral domain walls

abstract: The domain wall depinning field represents the minimum magnetic field needed to move a domain wall, typically pinned by samples' disorder or patterned constrictions. Conventionally, such field is considered independent on the Gilbert damping since it is assumed to be the field at which the Zeeman energy equals the pinning energy barrier (both damping independent). Here, we analyse numerically the domain wall depinning field as function of the Gilbert damping in a system with perpendicular magnetic anisotropy and Dzyaloshinskii-Moriya interaction. Contrary to expectations, we find that the depinning field depends on the Gilbert damping and that it strongly decreases for small damping parameters. We explain this dependence with a simple one-dimensional model and we show that the reduction of the depinning field is related to the internal domain wall dynamics, proportional to the Dzyaloshinskii-Moriya interaction, and the finite size of the pinning barriers.


title: MindDesktop: a general purpose brain computer interface

abstract: Recent advances in electroencephalography (EEG) and electromyography (EMG) enable communication for people with severe disabilities. In this paper we present a system that enables the use of regular computers using an off-the-shelf EEG/EMG headset, providing a pointing device and virtual keyboard that can be used to operate any Windows based system, minimizing the user effort required for interacting with a personal computer. Effectiveness of the proposed system is evaluated by a usability study, indicating decreasing learning curve for completing various tasks. The proposed system is available in the link provided.


title: On cohomological Hall algebras of quivers : Yangians

abstract: We consider the cohomological Hall algebra Y of a Lagrangian substack of the moduli stack of representations of the preprojective algebra of an arbitrary quiver Q, and its actions on the cohomology of quiver varieties. We conjecture that Y is equal, after a suitable extension of scalars, to the Yangian introduced by Maulik and Okounkov, and we construct an embedding of Y in the Yangian, intertwining the respective actions of both algebras on the cohomology of quiver varieties.


title: Parallel and in-process compilation of individuals for genetic programming on GPU

abstract: Three approaches to implement genetic programming on GPU hardware are compilation, interpretation and direct generation of machine code. The compiled approach is known to have a prohibitive overhead compared to other two. This paper investigates methods to accelerate compilation of individuals for genetic programming on GPU hardware. We apply in-process compilation to minimize the compilation overhead at each generation; and we investigate ways to parallelize in-process compilation. In-process compilation doesn't lend itself to trivial parallelization with threads; we propose a multiprocess parallelization using memory sharing and operating systems interprocess communication primitives. With parallelized compilation we achieve further reductions on compilation overhead. Another contribution of this work is the code framework we built in C# for the experiments. The framework makes it possible to build arbitrary grammatical genetic programming experiments that run on GPU with minimal extra coding effort, and is available as open source.


title: Queue Length Simulation for Signalized Arterial Networks and Steady State Computation under Fixed Time Control

abstract: We consider traffic flow dynamics for a network of signalized intersections, where the outflow from every link is constrained to be equal to a given capacity function if the queue length is positive, and equal to the minimum of cumulative inflow and capacity function otherwise. In spite of the resulting dynamics being discontinuous, recent work has proved existence and uniqueness of the resulting queue length trajectory if the inter-link travel times are strictly bounded away from zero. The proof, which also suggests a constructive procedure, relies on showing desired properties on contiguous time intervals of length equal to the minimum among all link travel times. We provide an alternate framework to obtain queue length trajectories by direct simulation of delay differential equations, where link outflows are obtained from the provably unique solution to a linear program. Existence and uniqueness of the solution to the proposed model for traffic flow dynamics is established for piecewise constant external inflow and capacity functions, and the proposed method does not require travel times to be bounded away from zero. Additionally, if the external inflow and capacity functions are periodic and satisfy a stability condition, then there exists a globally attractive periodic orbit. We provide an iterative procedure to compute this periodic orbit. A periodic trajectory is iteratively updated for every link based on updates to a specific time instant when its queue length transitions from being zero to being positive. The update for a given link is based on the periodic trajectories computed in the previous iteration for its upstream links. The resulting iterates are shown to converge uniformly monotonically to the desired periodic orbit.


title: Naturally graded Lie algebras (Carnot algebras) of slow growth

abstract: A nilpotent Lie algebra ${\mathfrak g}$ is said to be naturally graded if it is isomorphic to its associated graded Lie algebra ${\rm gr} \mathfrak{g}$ with respect to filtration by ideals of the lower central series. This concept is equivalent to the concept of the Carnot algebra arising in sub-Riemannian geometry and the geometric control theory. We classify finite-dimensional and infinite-dimensional naturally graded Lie algebras (Carnot algebras) ${\mathfrak g}=\oplus_{i=1}^{{+}\infty}{\mathfrak g}_i$ with properties $$ [{\mathfrak g}_1, {\mathfrak g}_i]={\mathfrak g}_{i{+}1}, \; \dim{{\mathfrak g}_i}+\dim{{\mathfrak g}_{i{+}1}} \le 3, \; i \ge 1. $$ For growth functions of such Lie algebras, we have the estimate $F(n) \le \frac{3}{2}n{+}1$.


title: Density Functional Theory for the Electron Gas and for Jellium

abstract: Density Functional Theory relies on universal functionals characteristic of a given system. Those functionals in general are different for the electron gas and for jellium (electron gas with uniform background). However, jellium is frequently used to construct approximate functionals for the electron gas (e.g., local density approximation, gradient expansions). The precise relationship of the exact functionals for the two systems is addressed here. In particular, it is shown that the exchange - correlation functionals for the inhomogeneous electron gas and inhomogeneous jellium are the same. This justifies theoretical and quantum Monte Carlo simulation studies of jellium to guide the construction of functionals for the electron gas. Related issues of the thermodynamic limit are noted as well.


title: Almost Rigidity of the Positive Mass Theorem for Asymptotically Hyperbolic Manifolds with Spherical Symmetry

abstract: We use the notion of intrinsic flat distance to address the almost rigidity of the positive mass theorem for asymptotically hyperbolic manifolds. In particular, we prove that a sequence of spherically symmetric asymptotically hyperbolic manifolds satisfying the conditions of the positive mass theorem converges to hyperbolic space in the intrinsic flat sense, if the limit of the mass along the sequence is zero.


title: Accelerating optimization-based computed tomography via sparse matrix approximations

abstract: Variational formulations of reconstruction in computed tomography have the notable drawback of requiring repeated evaluations of both the forward Radon transform and either its adjoint or an approximate inverse transform which are relatively expensive. We look at two methods for reducing the effect of this resulting computational bottleneck via approximating the transform evaluation with sparse matrix multiplications. The first method is applicable for general iterative optimization algorithms. The second is applicable in error-forgetting algorithms such as split Bregman. We demonstrate these approximations significantly reduce the needed computational time needed for the iterative algorithms needed to solve the reconstruction problem while still providing good reconstructions.


title: Angles of Gaussian primes

abstract: Fermat showed that every prime p = 1 mod 4 is a sum of two squares: $p = a^2 + b^2$. To any of the 8 possible representations (a,b) we associate an angle whose tangent is the ratio b/a. In 1919 Hecke showed that these angles are uniformly distributed as p varies, and in the 1950's Kubilius proved uniform distribution in somewhat short arcs. We study fine scale statistics of these angles, in particular the variance of the number of such angles in a short arc. We present a conjecture for this variance, motivated both by a random matrix model, and by a function field analogue of this problem, for which we prove an asymptotic form for the corresponding variance.


title: On the homotopy type of the space of Sullivan diagrams

abstract: We study the homotopy type of the harmonic compactification of the moduli space of a 2-cobordism S with one outgoing boundary component, or equivalently of the space of Sullivan diagrams of type S on one circle. Our results are of two types: vanishing and non-vanishing. In our vanishing results we are able to show that the connectivity of the harmonic compactification increases with the number of incoming boundary components. Moreover, we extend the genus stabilization maps of moduli spaces to the harmonic compactification and show that the connectivity of these maps increases with the genus and number of incoming boundary components. In our non-vanishing results we compute the non-trivial fundamental group of the harmonic compactification of the cobordism S of any genus with two unenumerated punctures and empty incoming boundary. Moreover, we construct five infinite families of non-trivial homology classes of the harmonic compactification, two of which correspond to non-trivial higher string topology operations.


title: The real plane Cremona group is a non-trivial amalgam

abstract: We show that the real Cremona group of the plane is a non-trivial amalgam of two groups amalgamated along their intersection and give an alternative proof of its abelianisation.


title: Roughness scattering induced insulator-metal-insulator transition in a quantum wire

abstract: We have theoretically investigated the influence of interface roughness scattering on the low temperature mobility of electrons in quantum wires when electrons fill one or many subbands. We find the Drude conductance of the wire with length $\mathcal{L}$ first increases with increasing linear concentration of electrons $\eta$ and then decreases at larger concentrations. For small radius $R$ of the wire with length $\mathcal{L}$ the peak of the conductance $G_{max}$ is below $e^2/h$ so that electrons are localized. The height of this peak grows as a large power of $R$, so that at large $R$ the conductance $G_{max}$ exceeds $e^2/h$ and a window of concentrations with delocalized states (which we call the metallic window) opens around the peak. Thus, we predict an insulator-metal-insulator transition with increasing concentration for large enough $R$. Furthermore, we show that the metallic domain can be sub-divided into three smaller domains: 1) single-subband ballistic conductor, 2) many-subband ballistic conductor 3) diffusive metal, and use our results to estimate the conductance in these domains. Finally we estimate the critical value of $R_c(\mathcal{L})$ at which the metallic window opens for a given length $\mathcal{L}$ and find it to be in reasonable agreement with experiment.


title: Permutation groups containing a regular abelian subgroup: the tangled history of two mistakes of Burnside

abstract: A group K is said to be a B-group if every permutation group containing K as a regular subgroup is either imprimitive or 2-transitive. In the second edition of his influential textbook on finite groups, Burnside published a proof that cyclic groups of composite prime-power degree are B-groups. Ten years later in 1921 he published a proof that every abelian group of composite degree is a B-group. Both proofs are character-theoretic and both have serious flaws. Indeed, the second result is false. In this note we explain these flaws and prove that every cyclic group of composite order is a B-group, using only Burnside's character-theoretic methods. We also survey the related literature, prove some new results on B-groups of prime-power order, state two related open problems and present some new computational data.


title: A QCD Lagrangian including renormalizable NJL terms

abstract: A local and gauge invariant version of QCD Lagrangian is introduced. The model includes Nambu-Jona-Lasinio (NJL) terms within its action in a surprisingly renormalizable form. This occurs thanks to the presence of action terms which at first sight, look as breaking power counting renormalizability. However, those terms also modify the quark propagators, to become more decreasing that the Dirac one at large momenta, indicating power counting renormalizability. The approach, can also be interpreted as generalized renormalization procedure for massless QCD. The free propagator, given by the substraction between a massive and a massless Dirac ones, in the Lee-Wick form, suggests that the theory also retains unitarity. The appearance of finite quark masses already in the tree approximation in the scheme is determined by the fact that the new action terms explicitly break chiral invariance. The approach looks as being able to implement the Fritzsch Democratic Symmetry breaking ideas about the quark mass hierarchy. Also, a link of the theory with the SM seems that can follow after employing the Zimmermann's couplings reduction scheme. The renormalized Feynman diagram expansion of the model is written here and the formula for the degree of divergence of the diagrams is derived. The primitive divergent graphs are identified. We start their evaluation through the calculation of the divergence of the gluon polarization operator. The divergent part of the gluon selfenergy results to be transverse as required by gauge invariance and also independent of the new couplings appearing. The calculations of the rest of primitive divergences will be presented elsewhere.


title: On some polynomials and series of Bloch-Polya Type

abstract: We will show that $(1-q)(1-q^2)\dots (1-q^m)$ is a polynomial in $q$ with coefficients from $\{-1,0,1\}$ iff $m=1,\ 2,\ 3,$ or $5$ and explore some interesting consequences of this result. We find explicit formulas for the $q$-series coefficients of $(1-q^2)(1-q^3)(1-q^4)(1-q^5)\dots$ and $(1-q^3)(1-q^4)(1-q^5)(1-q^6)\dots$. In doing so, we extend certain observations made by Sudler in 1964. We also discuss the classification of the products $(1-q)(1-q^2)\dots (1-q^m)$ and some related series with respect to their absolute largest coefficients.


title: Annealed Generative Adversarial Networks

abstract: We introduce a novel framework for adversarial training where the target distribution is annealed between the uniform distribution and the data distribution. We posited a conjecture that learning under continuous annealing in the nonparametric regime is stable irrespective of the divergence measures in the objective function and proposed an algorithm, dubbed {\ss}-GAN, in corollary. In this framework, the fact that the initial support of the generative network is the whole ambient space combined with annealing are key to balancing the minimax game. In our experiments on synthetic data, MNIST, and CelebA, {\ss}-GAN with a fixed annealing schedule was stable and did not suffer from mode collapse.


title: A Blind Deconvolution Technique Based on Projection Onto Convex Sets for Magnetic Particle Imaging

abstract: Magnetic Particle Imaging (MPI) is an emerging imaging modality that maps the spatial distribution of magnetic nanoparticles. The x-space reconstruction in MPI results in highly blurry images, where the resolution depends on both system parameters and nanoparticle type. Previous techniques to counteract this blurring rely on the knowledge of the imaging point spread function (PSF), which may not be available or may require additional measurements. This work proposes a blind deconvolution algorithm for MPI to recover the precise spatial distribution of nanoparticles. The proposed algorithm exploits the observation that the imaging PSF in MPI has zero phase in Fourier domain. Thus, even though the reconstructed images are highly blurred, phase remains unaltered. We leverage this powerful property to iteratively enforce consistency of phase and bounded l1 energy information, using an orthogonal Projections Onto Convex Sets (POCS) algorithm. To demonstrate the method, comprehensive simulations were performed without and with nanoparticle relaxation effects, and at various noise levels. In addition, imaging experiments were performed on an in-house MPI scanner using a three-vial phantom that contained different nanoparticle types. Image quality was compared with conventional deconvolution methods, Wiener deconvolution and Lucy-Richardson method, which explicitly rely on the knowledge of PSF. Both the simulation results and experimental imaging results show that the proposed blind deconvolution algorithm outperforms the conventional deconvolution methods. Without utilizing the imaging PSF, the proposed algorithm improves image quality and resolution even in the case of different nanoparticle types, while displaying reliable performance against loss of the fundamental harmonic, nanoparticle relaxation effects, and noise.


title: Analysis of Krylov Subspace Approximation to Large Scale Differential Riccati Equations

abstract: We consider a Krylov subspace approximation method for the symmetric differential Riccati equation $\dot{X} = AX + XA^T + Q - XSX$, $X(0)=X_0$. The method we consider is based on projecting the large scale equation onto a Krylov subspace spanned by the matrix $A$ and the low rank factors of $X_0$ and $Q$. We prove that the method is structure preserving in the sense that it preserves two important properties of the exact flow, namely the positivity of the exact flow, and also the property of monotonicity. We also provide a theoretical a priori error analysis which shows a superlinear convergence of the method. This behavior is illustrated in the numerical experiments. Moreover, we derive an efficient a posteriori error estimate as well as discuss multiple time stepping combined with a cut of the rank of the numerical solution.


title: A Comparative study of catalytic activity and lifetime of novel micro-meso porous catalysts in MTO process

abstract: Recently, two kinds of mesoporous catalysts with high propylene selectivity in propane dehydrogenation (PDH) process has been successfully synthesized. The first proposed catalyst is SAPO-34 molecular sieve with hierarchical tuned nanostructure. The second catalyst is a novel bi-phase SAPO-34/ZSM-5 zeolite hierarchical composite utilized with TPA-SAPO-34 exchanged core cystals being wrapped by ZSM-5 zeolite particles. In this contribution, the physico-chemical properties of the catalysts were analyzed by XRD, FESEM and N2 adsorption-desorption techniques and their catalytic activity and life time were investigated in MTO process. The results show that hierarchical SAPO-34 has a significant lifetime and selectivity to the light olefins compared with SAPO-34/ZSM-5. It shows a full conversion of MeOH during the first 200 min of reaction while composite sample conversion starts under 90 % and decreases during the time. However the deactivation phenomenon is hindered by using core-shell catalyst; hence the robustness can be named as a privilege employing composite zeolitic catalysts.


title: A truncation model for estimating Species Richness

abstract: We propose a truncation model for abundance distribution in the species richness estimation. This model is inherently semiparametric and incorporates an unknown truncation threshold between rare and abundant counts observations. Using the conditional likelihood, we derive a class of estimators for the parameters in the model by a stepwise maximisation. The species richness estimator is given by the integer maximising the binomial likelihood when all other parameters in the model are know. Under regularity conditions, we show that the estimators of the model parameters are asymptotically efficient. We recover the Chao$^{'}$s lower bound estimator of species richeness when the model is a unicomponent Poisson$^{'}$s model. So, it is an element of our class of estimators. In a simulation study, we show the performances of the proposed method and compare it to some others.


title: Search for Heavy Neutrinos in $K^+ \rightarrow \mu^+ \nu_{\mu}$ Decays

abstract: The NA62 experiment recorded a large sample of $K^+ \rightarrow \mu^+ \nu_{\mu}$ decays in 2007. A peak search has been performed in the reconstructed missing mass spectrum. In the absence of a signal, limits in the range $2 \times 10^{-6}$ to $10^{-5}$ have been set on the squared mixing matrix element $|U_{\mu4} |^2$ between muon and heavy neutrino states, for heavy neutrino masses in the range 300-375 MeV/$c^2$. The result extends the range of masses for which upper limits have been set on the value of $|U_{\mu4} |^2$ in previous production search experiments.


title: ARABIS: an Asynchronous Acoustic Indoor Positioning System for Mobile Devices

abstract: Acoustic ranging based indoor positioning solutions have the advantage of higher ranging accuracy and better compatibility with commercial-off-the-self consumer devices. However, similar to other time-domain based approaches using Time-of-Arrival and Time-Difference-of-Arrival, they suffer from performance degradation in presence of multi-path propagation and low received signal-to-noise ratio (SNR) in indoor environments. In this paper, we improve upon our previous work on asynchronous acoustic indoor positioning and develop ARABIS, a robust and low-cost acoustic positioning system (IPS) for mobile devices. We develop a low-cost acoustic board custom-designed to support large operational ranges and extensibility. To mitigate the effects of low SNR and multi-path propagation, we devise a robust algorithm that iteratively removes possible outliers by taking advantage of redundant TDoA estimates. Experiments have been carried in two testbeds of sizes 10.67m*7.76m and 15m*15m, one in an academic building and one in a convention center. The proposed system achieves average and 95% quantile localization errors of 7.4cm and 16.0cm in the first testbed with 8 anchor nodes and average and 95% quantile localization errors of 20.4cm and 40.0cm in the second testbed with 4 anchor nodes only.


title: Pyramid: Enhancing Selectivity in Big Data Protection with Count Featurization

abstract: Protecting vast quantities of data poses a daunting challenge for the growing number of organizations that collect, stockpile, and monetize it. The ability to distinguish data that is actually needed from data collected "just in case" would help these organizations to limit the latter's exposure to attack. A natural approach might be to monitor data use and retain only the working-set of in-use data in accessible storage; unused data can be evicted to a highly protected store. However, many of today's big data applications rely on machine learning (ML) workloads that are periodically retrained by accessing, and thus exposing to attack, the entire data store. Training set minimization methods, such as count featurization, are often used to limit the data needed to train ML workloads to improve performance or scalability. We present Pyramid, a limited-exposure data management system that builds upon count featurization to enhance data protection. As such, Pyramid uniquely introduces both the idea and proof-of-concept for leveraging training set minimization methods to instill rigor and selectivity into big data management. We integrated Pyramid into Spark Velox, a framework for ML-based targeting and personalization. We evaluate it on three applications and show that Pyramid approaches state-of-the-art models while training on less than 1% of the raw data.


title: Nonequilibrium diffusive gas dynamics: Poiseuille microflow

abstract: We test the recently developed hierarchy of diffusive moment closures for gas dynamics together with the near-wall viscosity scaling on the Poiseuille flow of argon and nitrogen in a one micrometer wide channel, and compare it against the corresponding Direct Simulation Monte Carlo computations. We find that the diffusive regularized Grad equations with viscosity scaling provide the most accurate approximation to the benchmark DSMC results. At the same time, the conventional Navier-Stokes equations without the near-wall viscosity scaling are found to be the least accurate among the tested closures.


title: A method for construction of rational points over elliptic curves

abstract: I provide a systematic construction of points (defined over number fields) on Legendre elliptic curves over $\mathbb{Q}$: for any odd integer $n\geq 3$ my method constructs $n$ points on the Legendre curve and I show that rank of the subgroup of the Mordell-Weil group they generate is $n$ if $n\geq 7$. I also show that every elliptic curve over any number field admits similar type of points after a finite base extension.


title: Effective Field Theory of Post-Newtonian Gravity Including Spins

abstract: We present in detail an Effective Field Theory (EFT) formulation for the essential case of spinning objects as the components of inspiralling compact binaries. We review its implementation, carried out in a series of works in recent years, which leveled the high post-Newtonian (PN) accuracy in the spinning sector to that, recently attained in the non-spinning sector. We note a public package, "EFTofPNG", that we recently created for high precision computation in the EFT of PN Gravity, which covers all sectors, and includes an observables pipeline.


title: Matter-antimatter coexistence method for finite density QCD toward a solution of the sign problem

abstract: Toward the lattice QCD calculation at finite density, we propose "matter-antimatter coexistence method", where matter and anti-matter systems are prepared on two parallel ${\bf R}^4$-sheets in five-dimensional Euclidean space-time. We put a matter system $M$ with a chemical potential $\mu \in {\bf C}$ on a ${\bf R}^4$-sheet, and also put an anti-matter system $\bar M$ with $-\mu^*$ on the other ${\bf R}^4$-sheet shifted in the fifth direction. Between the gauge variables $U_\nu \equiv e^{iagA_\nu}$ in $M$ and $\tilde U_\nu \equiv e^{iag \tilde A_\nu}$ in $\bar M$, we introduce a correlation term with a real parameter $\lambda$. In one limit of $\lambda \rightarrow \infty$, a strong constraint $\tilde U_\nu(x)=U_\nu(x)$ is realized, and therefore the total fermionic determinant becomes real and non-negative, due to the cancellation of the phase factors in $M$ and $\bar M$, although this system resembles QCD with an isospin chemical potential. In another limit of $\lambda \rightarrow 0$, this system goes to two separated ordinary QCD systems with the chemical potential of $\mu$ and $-\mu^*$. For a given finite-volume lattice, if one takes an enough large value of $\lambda$, $\tilde U_\nu(x) \simeq U_\nu(x)$ is realized and phase cancellation approximately occurs between two fermionic determinants in $M$ and $\bar M$, which suppresses the sign problem and is expected to make the lattice calculation possible. For the obtained gauge configurations of the coexistence system, matter-side quantities are evaluated through their measurement only for the matter part $M$. The physical quantities in finite density QCD are expected to be estimated by the calculations with gradually decreasing $\lambda$ and the extrapolation to $\lambda=0$. We also consider more sophisticated improvement of this method using an irrelevant-type correlation.


title: Free energy for non-equilibrium quasi-stationary states

abstract: We study a class of non-equilibrium quasi-stationary states for a Markov system interacting with two different thermal baths. We show that the work done under a slow, external change of parameters admits a potential, i.e., the free energy. Three conditions are needed for the existence of free energy in this non-equilibrium system: time-scale separation between variables of the system, partial controllability (external fields couple only with the slow variable), and an effective detailed balance. These conditions are facilitated in the continuous limit for the slow variable. In contrast to its equilibrium counterpart, the non-equilibrium free energy can increase with temperature. One example of this is that entropy reduction by means of external fields (cooling) can be easier (in the sense of the work cost) if it starts from a higher temperature.


title: Using a non-minimal coupling between matter and curvature to sequester the Cosmological Constant

abstract: We present a novel mechanism for generating a Cosmological Constant and suitably sequestering the vacuum contribution to it, so that the eponymous Cosmological Constant problem is avoided. We do so by resorting to a model endowed with a non-minimal coupling between curvature and matter in an appropriately defined relaxed regime, and show that this shares features with both Unimodular gravity as well as a recent proposal to sequester the vacuum contribution through the use of an external term to the action functional.


title: Sandpile Groups of Random Bipartite Graphs

abstract: We determine the asymptotic distribution of the p-rank of the sandpile groups of random bipartite graphs. We see that this depends on the ratio between the number of vertices on each side, with a threshold when the ratio between the sides is equal to 1/p. We follow the approach of Melanie Wood and consider random graphs as a special case of random matrices, and rely on a variant the definition of min-entropy given by Maples, in order to obtain useful results about these random matrices. Our results show that unlike the sandpile groups of Erdos-Renyi random graphs, the distribution of the sandpile groups of random bipartite graphs depends on the properties of the graph, rather than coming from some more general random group model.


title: Rewriting Context-free Families of String Diagrams

abstract: String diagrams provide a convenient graphical framework which may be used for equational reasoning about morphisms of monoidal categories. However, unlike term rewriting, rewriting string diagrams results in shorter equational proofs, because the string diagrammatic representation allows us to formally establish equalities modulo any rewrite steps which follow from the monoidal structure. Manipulating string diagrams by hand is a time-consuming and error-prone process, especially for large string diagrams. This can be ameliorated by using software proof assistants, such as Quantomatic. However, reasoning about concrete string diagrams may be limiting and in some scenarios it is necessary to reason about entire (infinite) families of string diagrams. When doing so, we face the same problems as for manipulating concrete string diagrams, but in addition, we risk making further mistakes if we are not precise enough about the way we represent (infinite) families of string diagrams. The primary goal of this thesis is to design a mathematical framework for equational reasoning about infinite families of string diagrams which is amenable to computer automation. We will be working with context-free families of string diagrams and we will represent them using context-free graph grammars. We will model equations between infinite families of diagrams using rewrite rules between context-free grammars. Our framework represents equational reasoning about concrete string diagrams and context-free families of string diagrams using double-pushout rewriting on graphs and context-free graph grammars respectively. We will prove that our representation is sound by showing that it respects the concrete semantics of string diagrammatic reasoning and we will show that our framework is appropriate for software implementation by proving important decidability properties.


title: Neon in ultrashort and intense x rays from free electron lasers

abstract: We theoretically examine neon atoms in ultrashort and intense x rays from free electron lasers and compare our results with data from experiments conducted at the Linac Coherent Light Source (LCLS). For this purpose, we treat in detail the electronic structure in all possible nonrelativistic cationic configurations using a relativistic multiconfiguration approach. The interaction with the x rays is described in rate-equation approximation. To understand the mechanisms of the interaction, a path analysis is devised which allows us to investigate what sequences of photoionization and decay processes lead to a specific configuration and with what probability. Thereby, we uncover a connection to the mathematics of graph theory and formal languages. In detail, we study the ion yields and find that plain rate equations do not provide a satisfactory description. We need to extend the rate equations for neon to incorporate double Auger decay of a $K$-shell vacancy and photoionization shake off for neutral neon. Shake off is included for valence and core ionization; the former has hitherto been overlooked but has important consequences for the ion yields from an x-ray energy below the core ionization threshold. Furthermore, we predict the photon yields from XUV and x-ray fluorescence; these allow one insights into the configurations populated by the interaction with the x rays. Finally, we discover that inaccuracies in those Auger decay widths employed in previous studies have only a minor influence on ion and photon yields.


title: Classification and Retrieval of Digital Pathology Scans: A New Dataset

abstract: In this paper, we introduce a new dataset, \textbf{Kimia Path24}, for image classification and retrieval in digital pathology. We use the whole scan images of 24 different tissue textures to generate 1,325 test patches of size 1000$\times$1000 (0.5mm$\times$0.5mm). Training data can be generated according to preferences of algorithm designer and can range from approximately 27,000 to over 50,000 patches if the preset parameters are adopted. We propose a compound patch-and-scan accuracy measurement that makes achieving high accuracies quite challenging. In addition, we set the benchmarking line by applying LBP, dictionary approach and convolutional neural nets (CNNs) and report their results. The highest accuracy was 41.80\% for CNN.


title: Jordan properties of automorphism groups of certain open algebraic varieties

abstract: Let $W$ be a quasiprojective variety over an algebraically closed field of characteristic zero. Assume that $W$ is birational to a product of a smooth projective variety $A$ and the projective line. We prove that if $A$ contains no rational curves then the automorphism group $G:=Aut(W)$ of $W$ is Jordan. That means that there is a positive integer $J=J(W)$ such that every finite subgroup $\mathcal{B}$ of ${G}$ contains a commutative subgroup $\mathcal{A}$ such that $\mathcal{A}$ is normal in $\mathcal{B}$ and the index $[\mathcal{B}:\mathcal{A}] \le J$ .


title: EDOS: Edge Assisted Offloading System for Mobile Devices

abstract: Offloading resource-intensive jobs to the cloud and nearby users is a promising approach to enhance mobile devices. This paper investigates a hybrid offloading system that takes both infrastructure-based networks and Ad-hoc networks into the scope. Specifically, we propose EDOS, an edge assisted offloading system that consists of two major components, an Edge Assistant (EA) and Offload Agent (OA). EA runs on the routers/towers to manage registered remote cloud servers and local service providers and OA operates on the users' devices to discover the services in proximity. We present the system with a suite of protocols to collect the potential service providers and algorithms to allocate tasks according to user-specified constraints. To evaluate EDOS, we prototype it on commercial mobile devices and evaluate it with both experiments on a small-scale testbed and simulations. The results show that EDOS is effective and efficient for offloading jobs.


title: Lorentz Gauge Theory of Gravity in Electron Positron Colliders

abstract: Lorentz gauge theory (LGT) is a feasible candidate for theory of quantum gravity in which routine field theory calculations can be carried out perturbatively without encountering too many divergences. In LGT spin of matter also gravitates. The spin-generated gravity is expected to be extremely stronger than that generated by mass and could be explored in current colliders. In this article the observable signals of the theory in an electron-positron collider is investigated. We specifically study pair annihilation into two gravitons, and LGT corrections to processes like $e^-+e^+\rightarrow \mu^-+\mu^+$ and $e^-+e^+\rightarrow e^-+e^+$.


title: Network Structure of Two-Dimensional Decaying Isotropic Turbulence

abstract: The present paper reports on our effort to characterize vortical interactions in complex fluid flows through the use of network analysis. In particular, we examine the vortex interactions in two-dimensional decaying isotropic turbulence and find that the vortical interaction network can be characterized by a weighted scale-free network. It is found that the turbulent flow network retains its scale-free behavior until the characteristic value of circulation reaches a critical value. Furthermore, we show that the two-dimensional turbulence network is resilient against random perturbations but can be greatly influenced when forcing is focused towards the vortical structures that are categorized as network hubs. These findings can serve as a network-analytic foundation to examine complex geophysical and thin-film flows and take advantage of the rapidly growing field of network theory, which complements ongoing turbulence research based on vortex dynamics, hydrodynamic stability, and statistics. While additional work is essential to extend the mathematical tools from network analysis to extract deeper physical insights of turbulence, an understanding of turbulence based on the interaction-based network-theoretic framework presents a promising alternative in turbulence modeling and control efforts.


title: Testing Degree Corrections in Stochastic Block Models

abstract: We study sharp detection thresholds for degree corrections in Stochastic Block Models in the context of a goodness of fit problem, and explore the effect of the unknown community assignment (a high dimensional nuisance parameter) and the graph density on testing for degree corrections. When degree corrections are relatively dense, a simple test based on the total number of edges is asymptotically optimal. For sparse degree corrections, the results undergo several changes in behavior depending on density of the underlying Stochastic Block Model. For graphs which are not extremely sparse, optimal tests are based on Higher Criticism or Maximum Degree type tests based on a linear combination of within and across (estimated) community degrees. In the special case of balanced communities, a simple degree based Higher Criticism Test (Mukherjee, Mukherjee, Sen 2016) is optimal in case the graph is not completely dense, while the more complicated linear combination based procedure is required in the completely dense setting. The ``necessity" of the two step procedure is demonstrated for the case of balanced communities by the failure of the ordinary Maximum Degree Test in achieving sharp constants. Finally for extremely sparse graphs the optimal rates change, and a version of the maximum degree test with a different rejection region is shown to be optimal.


title: The role of gap edge instabilities in setting the depth of planet gaps in protoplanetary discs

abstract: It is known that an embedded massive planet will open a gap in a protoplanetary disc via angular momentum exchange with the disc material. The resulting surface density profile of the disc is investigated for one dimensional and two dimensional disc models and, in agreement with previous work, it is found that one dimensional gaps are significantly deeper than their two dimensional counterparts for the same initial conditions. We find, by applying one dimensional torque density distributions to two dimensional discs containing no planet, that the excitement of the Rossby wave instability and the formation of Rossby vortices play a critical role in setting the equilibrium depth of the gap. Being a two dimensional instability, this is absent from one dimensional simulations and does not limit the equilibrium gap depth there. We find similar gap depths between two dimensional gaps formed by torque density distributions, in which the Rossby wave instability is present, and two dimensional planet gaps, in which no Rossby wave instability is present. This can be understood if the planet gap is maintained at marginal stability, even when there is no obvious Rossby wave instability present. Further investigation shows the final equilibrium gap depth is very sensitive to the form of the applied torque density distribution, and using improved one dimensional approximations from three dimensional simulations can go even further to reducing the discrepancy between one and two dimensional models, especially for lower mass planets. This behaviour is found to be consistent across discs with varying parameters.


title: Testing hypotheses on a tree: new error rates and controlling strategies

abstract: We introduce a multiple testing procedure (TreeBH) which addresses the challenge of controlling error rates at multiple levels of resolution. Conceptually, we frame this problem as the selection of hypotheses which are organized hierarchically in a tree structure. We describe a fast algorithm for the proposed sequential procedure, and prove that it controls relevant error rates given certain assumptions on the dependence among the p-values. Through simulations, we demonstrate that TreeBH offers the desired guarantees under a range of dependency structures (including one similar to that encountered in genome-wide association studies) and that it has the potential of gaining power over alternative methods. We also introduce a modified version of TreeBH which we prove to control the relevant error rates under any dependency structure. We conclude with two case studies: we first analyze data collected as part of the Genotype-Tissue Expression (GTEx) project, which aims to characterize the genetic regulation of gene expression across multiple tissues in the human body, and secondly, data examining the relationship between the gut microbiome and colorectal cancer.


title: Classification of toric manifolds over an $n$-cube with one vertex cut

abstract: We say that a complete nonsingular toric variety (called a toric manifold in this paper) is over $P$ if its quotient by the compact torus is homeomorphic to $P$ as a manifold with corners. Bott manifolds (or Bott towers) are toric manifolds over an $n$-cube $I^n$ and blowing them up at a fixed point produces toric manifolds over $\mathrm{vc}(I^n)$ an $n$-cube with one vertex cut. They are all projective. On the other hand, Oda's $3$-fold, the simplest non-projective toric manifold, is over $\mathrm{vc}(I^n)$. In this paper, we classify toric manifolds over $\mathrm{vc}(I^n)$ $(n\ge 3)$ as varieties and also as smooth manifolds. As a consequence, it turns out that (1) there are many non-projective toric manifolds over $\mathrm{vc}(I^n)$ but they are all diffeomorphic, and (2) toric manifolds over $\mathrm{vc}(I^n)$ in some class are determined by their cohomology rings as varieties among toric manifolds.


title: Corrupted Sensing with Sub-Gaussian Measurements

abstract: This paper studies the problem of accurately recovering a structured signal from a small number of corrupted sub-Gaussian measurements. We consider three different procedures to reconstruct signal and corruption when different kinds of prior knowledge are available. In each case, we provide conditions for stable signal recovery from structured corruption with added unstructured noise. The key ingredient in our analysis is an extended matrix deviation inequality for isotropic sub-Gaussian matrices.


title: Persistent Flows and Non-Reciprocal Interactions in Deterministic Networks

abstract: This paper studies deterministic consensus networks with discrete-time dynamics under persistent flows and non-reciprocal agent interactions. An arc describing the interaction strength between two agents is said to be persistent if its weight function has an infinite $l_1$ norm. We discuss two balance conditions on the interactions between agents which generalize the arc-balance and cut-balance conditions in the literature respectively. The proposed conditions require that such a balance should be satisfied over each time window of a fixed length instead of at each time instant. We prove that in both cases global consensus is reached if and only if the persistent graph, which consists of all the persistent arcs, contains a directed spanning tree. The convergence rates of the system to consensus are also provided in terms of the interactions between agents having taken place. The results are obtained under a weak condition without assuming the existence of a positive lower bound of all the nonzero weights of arcs and are compared with the existing results. Illustrative examples are provided to show the critical importance of the nontrivial lower boundedness of the self-confidence of the agents.


title: A Note on the Information-Theoretic-(in)Security of Fading Generated Secret Keys

abstract: In this work we explore the security of secret keys generated via the electromagnetic reciprocity of the wireless fading channel. Identifying a new sophisticated colluding attack, we explore the information-theoretic-security for such keys in the presence of an all-powerful adversary constrained only by the laws of quantum mechanics. Specifically, we calculate the reduction in the conditional mutual information between transmitter and receiver that can occur when an adversary with unlimited computational and communication resources places directional antenna interceptors at chosen locations. Such locations, in principal, can be arbitrarily far from the intended receiver yet still influence the secret key rate.


title: Random walks among time increasing conductances: heat kernel estimates

abstract: For any graph having a suitable uniform Poincare inequality and volume growth regularity, we establish two-sided Gaussian transition density estimates and parabolic Harnack inequality, for constant speed continuous time random walks evolving via time varying, uniformly elliptic conductances, provided the vertex conductances (i.e. reversing measures), increase in time. Such transition density upper bounds apply for discrete time uniformly lazy walks, with the matching lower bounds holding once the parabolic Harnack inequality is proved.


title: Evading Classifiers by Morphing in the Dark

abstract: Learning-based systems have been shown to be vulnerable to evasion through adversarial data manipulation. These attacks have been studied under assumptions that the adversary has certain knowledge of either the target model internals, its training dataset or at least classification scores it assigns to input samples. In this paper, we investigate a much more constrained and realistic attack scenario wherein the target classifier is minimally exposed to the adversary, revealing on its final classification decision (e.g., reject or accept an input sample). Moreover, the adversary can only manipulate malicious samples using a blackbox morpher. That is, the adversary has to evade the target classifier by morphing malicious samples "in the dark". We present a scoring mechanism that can assign a real-value score which reflects evasion progress to each sample based on the limited information available. Leveraging on such scoring mechanism, we propose an evasion method -- EvadeHC -- and evaluate it against two PDF malware detectors, namely PDFRate and Hidost. The experimental evaluation demonstrates that the proposed evasion attacks are effective, attaining $100\%$ evasion rate on the evaluation dataset. Interestingly, EvadeHC outperforms the known classifier evasion technique that operates based on classification scores output by the classifiers. Although our evaluations are conducted on PDF malware classifier, the proposed approaches are domain-agnostic and is of wider application to other learning-based systems.


title: Integrable structure of products of finite complex Ginibre random matrices

abstract: We consider the squared singular values of the product of $M$ standard complex Gaussian matrices. Since the squared singular values form a determinantal point process with a particular Meijer G-function kernel, the gap probabilities are given by a Fredholm determinant based on this kernel. It was shown by Strahov \cite{St14} that a hard edge scaling limit of the gap probabilities is described by Hamiltonian differential equations which can be formulated as an isomonodromic deformation system similar to the theory of the Kyoto school. We generalize this result to the case of finite matrices by first finding a representation of the finite kernel in integrable form. As a result we obtain the Hamiltonian structure for a finite size matrices and formulate it in terms of a $(M+1) \times (M+1)$ matrix Schlesinger system. The case $M=1$ reproduces the Tracy and Widom theory which results in the Painlev\'e V equation for the $(0,s)$ gap probability. Some integrals of motion for $M = 2$ are identified, and a coupled system of differential equations in two unknowns is presented which uniquely determines the corresponding $(0,s)$ gap probability.


title: A note on Li-Yau type gradient estimate

abstract: In this paper, we obtain Li-Yau type gradient estimates with time dependent parameter for positive solutions of the heat equation that are different with the estimates by Li-Xu \cite{LX} and Qian \cite{Qi}. As an application of the estimate, we also obtained improvements of Davies' Li-Yau type gradient estimate.


title: Infrastructure for Usable Machine Learning: The Stanford DAWN Project

abstract: Despite incredible recent advances in machine learning, building machine learning applications remains prohibitively time-consuming and expensive for all but the best-trained, best-funded engineering organizations. This expense comes not from a need for new and improved statistical models but instead from a lack of systems and tools for supporting end-to-end machine learning application development, from data preparation and labeling to productionization and monitoring. In this document, we outline opportunities for infrastructure supporting usable, end-to-end machine learning applications in the context of the nascent DAWN (Data Analytics for What's Next) project at Stanford.


title: On the Phase Transition of Corrupted Sensing

abstract: In \cite{FOY2014}, a sharp phase transition has been numerically observed when a constrained convex procedure is used to solve the corrupted sensing problem. In this paper, we present a theoretical analysis for this phenomenon. Specifically, we establish the threshold below which this convex procedure fails to recover signal and corruption with high probability. Together with the work in \cite{FOY2014}, we prove that a sharp phase transition occurs around the sum of the squares of spherical Gaussian widths of two tangent cones. Numerical experiments are provided to demonstrate the correctness and sharpness of our results.


title: An Overview of Massive MIMO Research at the University of Bristol

abstract: Massive MIMO has rapidly gained popularity as a technology crucial to the capacity advances required for 5G wireless systems. Since its theoretical conception six years ago, research activity has grown exponentially, and there is now a developing industrial interest to commercialise the technology. For this to happen effectively, we believe it is crucial that further pragmatic research is conducted with a view to establish how reality differs from theoretical ideals. This paper presents an overview of the massive MIMO research activities occurring within the Communication Systems & Networks Group at the University of Bristol centred around our 128-antenna real-time testbed, which has been developed through the BIO programmable city initiative in collaboration with NI and Lund University. Through recent preliminary trials, we achieved a world first spectral efficiency of 79.4 bits/s/Hz, and subsequently demonstrated that this could be increased to 145.6 bits/s/Hz. We provide a summary of this work here along with some of our ongoing research directions such as large-scale array wave-front analysis, optimised power control and localisation techniques.


title: Learning from Complementary Labels

abstract: Collecting labeled data is costly and thus a critical bottleneck in real-world classification tasks. To mitigate this problem, we propose a novel setting, namely learning from complementary labels for multi-class classification. A complementary label specifies a class that a pattern does not belong to. Collecting complementary labels would be less laborious than collecting ordinary labels, since users do not have to carefully choose the correct class from a long list of candidate classes. However, complementary labels are less informative than ordinary labels and thus a suitable approach is needed to better learn from them. In this paper, we show that an unbiased estimator to the classification risk can be obtained only from complementarily labeled data, if a loss function satisfies a particular symmetric condition. We derive estimation error bounds for the proposed method and prove that the optimal parametric convergence rate is achieved. We further show that learning from complementary labels can be easily combined with learning from ordinary labels (i.e., ordinary supervised learning), providing a highly practical implementation of the proposed method. Finally, we experimentally demonstrate the usefulness of the proposed methods.


title: Categorical relations between Langlands dual quantum affine algebras: Doubly laced types

abstract: We prove that the Grothendieck rings of category $\mathcal{C}^{(t)}_Q$ over quantum affine algebras $U_q'(\g^{(t)})$ $(t=1,2)$ associated to each Dynkin quiver $Q$ of finite type $A_{2n-1}$ (resp. $D_{n+1}$) is isomorphic to one of category $\mathcal{C}_{\mQ}$ over the Langlands dual $U_q'({^L}\g^{(2)})$ of $U_q'(\g^{(2)})$ associated to any twisted adapted class $[\mQ]$ of $A_{2n-1}$ (resp. $D_{n+1}$). This results provide partial answers of conjectures of Frenkel-Hernandez on Langlands duality for finite-dimensional representation of quantum affine algebras.


title: Building Emotional Machines: Recognizing Image Emotions through Deep Neural Networks

abstract: An image is a very effective tool for conveying emotions. Many researchers have investigated in computing the image emotions by using various features extracted from images. In this paper, we focus on two high level features, the object and the background, and assume that the semantic information of images is a good cue for predicting emotion. An object is one of the most important elements that define an image, and we find out through experiments that there is a high correlation between the object and the emotion in images. Even with the same object, there may be slight difference in emotion due to different backgrounds, and we use the semantic information of the background to improve the prediction performance. By combining the different levels of features, we build an emotion based feed forward deep neural network which produces the emotion values of a given image. The output emotion values in our framework are continuous values in the 2-dimensional space (Valence and Arousal), which are more effective than using a few number of emotion categories in describing emotions. Experiments confirm the effectiveness of our network in predicting the emotion of images.


title: Pre-inflationary universe in loop quantum cosmology

abstract: The evolutions of the flat FLRW universe and its linear perturbations are studied systematically in the dressed metric approach of LQC. When it is dominated by the kinetic energy of the inflaton at the quantum bounce, the evolution of the background can be divided into three different phases prior to the preheating, {\em bouncing, transition and slow-roll inflation}. During the bouncing phase, the evolution is independent of not only the initial conditions, but also the inflationary potentials. In particular, the expansion factor can be well described by the same exact solution in all the cases considered. In contrast, in the potential dominated case such a universality is lost. It is because of this universality that the linear perturbations are also independent of the inflationary models and obtained exactly. During the transition phase, the evolutions of the background and its linear perturbations are found explicitly, and then matched to the ones given in the other two phases. Hence, once the initial conditions are imposed, the linear scalar and tensor perturbations will be uniquely determined. Considering two different sets of initial conditions, one imposed during the contracting phase and the other at the bounce, we calculate the Bogoliubov coefficients and find that the two sets yield the same results and all lead to particle creations at the onset of the inflation. Due to the pre-inflationary dynamics, the scalar and tensor power spectra become scale-dependent. Comparing with the Planck 2015 data, we find constraints on the total e-folds that the universe must have expanded since the bounce, in order to be consistent with current observations.


title: A note on the number of edges in a Hamiltonian graph with no repeated cycle length

abstract: Let $G$ be an $n$-vertex graph obtained by adding chords to a cycle of length $n$. Markstr\"{o}m asked for the maximum number of edges in $G$ if there are no two cycles in $G$ with the same length. A simple counting argument shows that such a graph can have at most $n + \sqrt{2n} +1 $ edges. Using difference sets in $\mathbb{Z}_n$, we show that for infinitely many $n$, there is an $n$-vertex Hamiltonian graph with $n + \sqrt{n - 3/4} - 3/2$ edges and no repeated cycle length.


title: Half-Integral Weight Modular Forms and Modular forms for Weil representations

abstract: We establish an isomorphism between certain complex-valued and vector-valued modular form spaces of half-integral weight, generalizing the well-known isomorphism between modular forms for $\Gamma_0(4)$ with Kohnen's plus condition and modular forms for the Weil representation associated to the discriminant form for the lattice with Gram matrix $(2)$. With such an isomorphism, we prove the Zagier duality and write down the Borcherds lifts explicitly.


title: Hydrogen-modified inverse Spin Hall Effect in palladium-cobalt bi-layer films

abstract: The influence of hydrogen gas absorption by the Pd layer of bi-layered films containing Pd and Co layers on the inverse Spin Hall Effect (iSHE) in the material is measured. iSHE is driven by ferromagnetic resonance in the cobalt layer. In these conditions, the iSHE is seen as a d.c. voltage across the Pd layer. In the presence of hydrogen gas the iSHE peak shifts downwards in the applied field together with the ferromagnetic resonance absorption peak for the material. In parallel, an increase in the iSHE peak height and in its width is observed. Our analysis suggests that these observations can potentially be explained as a reduction in the spin Hall angle for the Palladium layer in the presence of hydrogen gas.


title: $K^+$-nucleus elastic scattering revisited from perspective of partial restoration of chiral symmetry

abstract: The $K^{+}$ meson properties in the nuclear medium are investigated by considering the wavefunction renormalization as a first step to reveal the in-medium properties of the $K^{+}$ meson in the context of partial restoration of chiral symmetry. The $K^{+}N$ elastic scattering amplitude is constructed using chiral perturbation theory up to the next-to-leading order. Using the constructed amplitude, we calculate the wavefunction renormalization in the Thomas-Fermi approximation. We obtained a good description of the $K^{+}N$ elastic scattering amplitude. The obtained wavefunction renormalization factor suggests the 2 to 6% enhancement of the $K^{+}N$ interaction at the saturation density. We conclude that the wavefunction renormalization could be one of the important medium effects for the $K^{+}$ meson.


title: The GIT moduli of semistable pairs consisting of a cubic curve and a line on ${\mathbb P}^{2}$

abstract: We discuss the GIT moduli of semistable pairs consisting of a cubic curve and a line on the projective plane. We study in some detail this moduli and compare it with another moduli suggested by Alexeev. It is the moduli of pairs (with no specified semi-abelian action) consisting of a cubic curve with at worst nodal singularities and a line which does not pass through singular points of the cubic curve. Meanwhile, we make a comparison between Nakamura's compactification of the moduli of level three elliptic curves and these two moduli spaces.


title: Local bifurcations in differential equations with state-dependent delay

abstract: A common task when analysing dynamical systems is the determination of normal forms near local bifurcations of equilibria. As most of these normal forms have been classified and analysed, finding which particular class of normal form one encounters in a numerical bifurcation study guides follow-up computations. This paper builds on normal form algorithms for equilibria of delay differential equations with constant delay that were developed and implemented in DDE-Biftool recently. We show how one can extend these methods to delay-differential equations with state-dependent delay (sd-DDEs). Since higher degrees of regularity of local center manifolds are still open for sd-DDEs, we give an independent (still only partial) argument which phenomena from the truncated normal must persist in the full sd-DDE. In particular, we show that all invariant manifolds with a sufficient degree of normal hyperbolicity predicted by the normal form exist also in the full sd-DDE.


title: Parameterized Complexity of the List Coloring Reconfiguration Problem with Graph Parameters

abstract: Let $G$ be a graph such that each vertex has its list of available colors, and assume that each list is a subset of the common set consisting of $k$ colors. For two given list colorings of $G$, we study the problem of transforming one into the other by changing only one vertex color assignment at a time, while at all times maintaining a list coloring. This problem is known to be PSPACE-complete even for bounded bandwidth graphs and a fixed constant $k$. In this paper, we study the fixed-parameter tractability of the problem when parameterized by several graph parameters. We first give a fixed-parameter algorithm for the problem when parameterized by $k$ and the modular-width of an input graph. We next give a fixed-parameter algorithm for the shortest variant when parameterized by $k$ and the size of a minimum vertex cover of an input graph. As corollaries, we show that the problem for cographs and the shortest variant for split graphs are fixed-parameter tractable even when only $k$ is taken as a parameter. On the other hand, we prove that the problem is W[1]-hard when parameterized only by the size of a minimum vertex cover of an input graph.


title: Predicted MAX phase Sc2InC: Dynamical stability, vibrational and optical properties

abstract: First principles pseudopotential calculations have been performed for the first time to investigate the phonon dispersion, thermodynamic and optical properties including charge density, Fermi surface, Mulliken population analysis, theoretical Vickers hardness of predicted MAX phase Sc2InC. We revisited the structural, elastic and electronic properties of the compound which assessed the reliability of our calculations. The analysis of the elastic constants and the phonon dispersion along with phonon density of states indicates the mechanical stability and dynamical stability of the MAX phase. The Helmholtz free energy, internal energy, entropy specific heat capacity and Debye temperature have also been calculated from the phonon density of states. Mulliken population analysis indicates the existence of prominent covalency in chemical bonding of Sc2InC. The electronic charge density mapping shows a combination of ionic, covalent and metallic bonding in the compound. The Fermi surface is comprised due to the low dispersive Sc 3d and C 2p states from the [ScC] blocks. The phase is expected to be a soft material and easily mechinable due to its low Vicker hardness value. Furthermore, the analysis of various optical properties (such as dielectric function, refractive index, photoconductivity, absorption coefficients, loss function and reflectivity) suggests that the nanolaminate Sc2InC is a promising candidate for optoelectronic devices in the visible and ultraviolet energy regions and as a coating material to avoid solar heating.


title: A Multi-telescope Campaign on FRB 121102: Implications for the FRB Population

abstract: We present results of the coordinated observing campaign that made the first subarcsecond localization of a Fast Radio Burst, FRB 121102. During this campaign, we made the first simultaneous detection of an FRB burst by multiple telescopes: the VLA at 3 GHz and the Arecibo Observatory at 1.4 GHz. Of the nine bursts detected by the Very Large Array at 3 GHz, four had simultaneous observing coverage at other observatories. We use multi-observatory constraints and modeling of bursts seen only at 3 GHz to confirm earlier results showing that burst spectra are not well modeled by a power law. We find that burst spectra are characterized by a ~500 MHz envelope and apparent radio energy as high as $10^{40}$ erg. We measure significant changes in the apparent dispersion between bursts that can be attributed to frequency-dependent profiles or some other intrinsic burst structure that adds a systematic error to the estimate of DM by up to 1%. We use FRB 121102 as a prototype of the FRB class to estimate a volumetric birth rate of FRB sources $R_{FRB} \approx 5x10^{-5}/N_r$ Mpc$^{-3}$ yr$^{-1}$, where $N_r$ is the number of bursts per source over its lifetime. This rate is broadly consistent with models of FRBs from young pulsars or magnetars born in superluminous supernovae or long gamma-ray bursts, if the typical FRB repeats on the order of thousands of times during its lifetime.


title: Dynamical Evolution of the Debris Disk after a Satellite Catastrophic Disruption around Saturn

abstract: The hypothesis of a recent origin of Saturn's rings and its mid-sized moons is actively debated. It was suggested that a proto-Rhea and a proto-Dione might have collided recently, giving birth to the modern system of mid-sized moons. It is also suggested that the rapid viscous spreading of the debris may have implanted mass inside Saturn's Roche limit, giving birth to the modern Saturn's ring system. However, this scenario has been only investigated in very simplified way for the moment. This paper investigates it in detail to assess its plausibility by using $N$-body simulations and analytical arguments. When the debris disk is dominated by its largest remnant, $N$-body simulations show that the system quickly re-accrete into a single satellite without significant spreading. On the other hand, if the disk is composed of small particles, analytical arguments suggest that the disk experiences dynamical evolutions in three steps. The disk starts significantly excited after the impact and collisional damping dominates over the viscous spreading. After the system flattens, the system can become gravitationally unstable when particles are smaller than $\sim$ 100 m. However, the particles grow faster than spreading. Then, the system becomes gravitationally stable again and accretion continues at a slower pace, but spreading is inhibited. Therefore, the debris is expected to re-accrete into several large bodies. In conclusion, our results show that such a scenario may not form the today's ring system. In contrast, our results suggest that today's mid-sized moons are likely re-accreted from such a catastrophic event.