title: A Relevant Content Filtering Based Framework For Data Stream Summarization

abstract: Social media platforms are a rich source of information these days, however, of all the available information, only a small fraction is of users' interest. To help users catch up with the latest topics of their interests from the large amount of information available in social media, we present a relevant content filtering based framework for data stream summarization. More specifically, given the topic or event of interest, this framework can dynamically discover and filter out relevant information from irrelevant information in the stream of text provided by social media platforms. It then further captures the most representative and up-to-date information to generate a sequential summary or event story line along with the evolution of the topic or event. Our framework does not depend on any labeled data, it instead uses the weak supervision provided by the user, which matches the real scenarios of users searching for information about an ongoing event. We experimented on two real events traced by a Twitter dataset from TREC 2011. The results verified the effectiveness of relevant content filtering and sequential summary generation of the proposed framework. It also shows its robustness of using the most easy-to-obtain weak supervision, i.e., trending topic or hashtag. Thus, this framework can be easily integrated into social media platforms such as Twitter to generate sequential summaries for the events of interest. We also make the manually generated gold-standard sequential summaries of the two test events publicly available for future use in the community.


title: An inverse problem from condense matter physics

abstract: We consider the problem of reconstructing the features of a weak anisotropic background potential by the trajectories of vortex dipoles in a nonlinear Gross-Pitaevskii equation. At leading order, the dynamics of vortex dipoles are given by a Hamiltonian system. If the background potential is sufficiently smooth and flat, the background can be reconstructed using ideas from the boundary and the lens rigidity problems. We prove that reconstructions are unique, derive an approximate reconstruction formula, and present numerical examples.


title: Local law for random Gram matrices

abstract: We prove a local law in the bulk of the spectrum for random Gram matrices $XX^*$, a generalization of sample covariance matrices, where $X$ is a large matrix with independent, centered entries with arbitrary variances. The limiting eigenvalue density that generalizes the Marchenko-Pastur law is determined by solving a system of nonlinear equations. Our entrywise and averaged local laws are on the optimal scale with the optimal error bounds. They hold both in the square case (hard edge) and in the properly rectangular case (soft edge). In the latter case we also establish a macroscopic gap away from zero in the spectrum of $XX^*$.


title: Constraining heavy decaying dark matter with the high energy gamma-ray limits

abstract: We consider decaying dark matter with masses $10^{7} \lesssim M \lesssim 10^{16}$ GeV, as a source of ultra-high energy (UHE) gamma rays. Using recent limits on UHE gamma-ray flux for energies $E_\gamma > 2 \cdot 10^{14}$ eV, provided by extensive air shower observatories, we put limits on masses and lifetimes of the dark matter. We also discuss possible dark matter decay origin of tentative 100 PeV photon flux detected with EAS-MSU experiment.


title: The ionization conjecture in Thomas-Fermi-Dirac-von Weizs\"acker theory

abstract: We prove that in Thomas-Fermi-Dirac-von Weizs\"acker theory, a nucleus of charge $Z>0$ can bind at most $Z+C$ electrons, where $C$ is a universal constant. This result is obtained through a comparison with Thomas-Fermi theory which, as a by-product, gives bounds on the screened nuclear potential and the radius of the minimizer. A key ingredient of the proof is a novel technique to control the particles in the exterior region, which also applies to the liquid drop model with a nuclear background potential.


title: Analyzing the Behavior of Visual Question Answering Models

abstract: Recently, a number of deep-learning based models have been proposed for the task of Visual Question Answering (VQA). The performance of most models is clustered around 60-70%. In this paper we propose systematic methods to analyze the behavior of these models as a first step towards recognizing their strengths and weaknesses, and identifying the most fruitful directions for progress. We analyze two models, one each from two major classes of VQA models -- with-attention and without-attention and show the similarities and differences in the behavior of these models. We also analyze the winning entry of the VQA Challenge 2016. Our behavior analysis reveals that despite recent progress, today's VQA models are "myopic" (tend to fail on sufficiently novel instances), often "jump to conclusions" (converge on a predicted answer after 'listening' to just half the question), and are "stubborn" (do not change their answers across images).


title: Parameter Optimisation of a Virtual Synchronous Machine in a Microgrid

abstract: Parameters of a virtual synchronous machine in a small microgrid are optimised. The dynamical behaviour of the system is simulated after a perturbation, where the system needs to return to its steady state. The cost functional evaluates the system behaviour for different parameters. This functional is minimised by Parallel Tempering. Two perturbation scenarios are investigated and the resulting optimal parameters agree with analytical predictions. Dependent on the focus of the optimisation different optima are obtained for each perturbation scenario. During the transient the system leaves the allowed voltage and frequency bands only for a short time if the perturbation is within a certain range.


title: Selection by Partitioning the Solution Paths

abstract: The performance of penalized likelihood approaches depends profoundly on the selection of the tuning parameter; however, there is no commonly agreed-upon criterion for choosing the tuning parameter. Moreover, penalized likelihood estimation based on a single value of the tuning parameter suffers from several drawbacks. This article introduces a novel approach for feature selection based on the entire solution paths rather than the choice of a single tuning parameter, which significantly improves the accuracy of the selection. Moreover, the approach allows for feature selection using ridge or other strictly convex penalties. The key idea is to classify variables as relevant or irrelevant at each tuning parameter and then to select all of the variables which have been classified as relevant at least once. We establish the theoretical properties of the method, which requires significantly weaker conditions than existing methods in the literature. We also illustrate the advantages of the proposed approach with simulation studies and a data example.


title: Magnetic eigenmaps for community detection in directed networks

abstract: Communities in directed networks have often been characterized as regions with a high density of links, or as sets of nodes with certain patterns of connection. Our approach for community detection combines the optimization of a quality function and a spectral clustering of a deformation of the combinatorial Laplacian, the so-called magnetic Laplacian. The eigenfunctions of the magnetic Laplacian, that we call magnetic eigenmaps, incorporate structural information. Hence, using the magnetic eigenmaps, dense communities including directed cycles can be revealed as well as "role" communities in networks with a running flow, usually discovered thanks to mixture models. Furthermore, in the spirit of the Markov stability method, an approach for studying communities at different energy levels in the network is put forward, based on a quantum mechanical system at finite temperature.


title: Low Power Wide Area Networks: An Overview

abstract: Low Power Wide Area (LPWA) networks are attracting a lot of attention primarily because of their ability to offer affordable connectivity to the low-power devices distributed over very large geographical areas. In realizing the vision of the Internet of Things (IoT), LPWA technologies complement and sometimes supersede the conventional cellular and short range wireless technologies in performance for various emerging smart city and machine-to-machine (M2M) applications. This review paper presents the design goals and the techniques, which different LPWA technologies exploit to offer wide-area coverage to low-power devices at the expense of low data rates. We survey several emerging LPWA technologies and the standardization activities carried out by different standards development organizations (e.g., IEEE, IETF, 3GPP, ETSI) as well as the industrial consortia built around individual LPWA technologies (e.g., LORa Alliance,WEIGHTLESS-SIG, and DASH7 Alliance). We further note that LPWA technologies adopt similar approaches, thus sharing similar limitations and challenges. This paper expands on these research challenges and identifies potential directions to address them. While the proprietary LPWA technologies are already hitting the market with large nationwide roll-outs, this paper encourages an active engagement of the research community in solving problems that will shape the connectivity of tens of billions of devices in the next decade.


title: Constraint Analysis of Two-Dimensional Quadratic Gravity from BF Theory

abstract: Quadratic gravity in two dimensions can be formulated as a Background Field (BF) theory plus an interaction term which is polynomial in both, the gauge and Background fields. This formulation is similar to the one given by Freidel and Starodubtsev to obtain MacDowell-Mansouri gravity in four dimensions. In this article we use the Dirac's Hamiltonian formalism to analyze the constraint structure of the two-dimensional Polynomial BF action. After we obtain the constraints of the theory, we proceed with the Batalin-Fradkin-Vilkovisky procedure to obtain the transition amplitude. We also compare our results with the ones obtained from generalized dilaton gravity.


title: Restricted developments in partizan misere game theory

abstract: Much progress has been made in misere game theory using the technique of restricted misere play, where games can be considered equivalent inside a restricted set of games without being equal in general. This paper provides a survey of recent results in this area, including two particularly interesting properties of restricted misere games: (1) a game can have an additive inverse that is not its negative, and (2) a position can be reversible through an end (a game with Left but not Right options, or vice versa). These properties are not possible in normal play and general misere play, respectively. Related open problems are discussed.


title: Coincidence Reidemeister trace and its generalization

abstract: We give a homotopy invariant construction of the Reidemeister trace for the coincidence of two maps between closed manifolds of not necessarily the same dimensions. It is realized as a homology class of the homotopy equalizer, which coincides with the Hurewicz image of Koschorke's stabilized bordism invariant. To define it, we use a kind of shriek maps appearing string topology. As an application, we compute the coincidence Reidemeister trace for the self-coincidence of the projections of $S^1$-bundles on $\mathbb{C}P^n$. We also mention how to relate our construction to the string topology operation called the loop coproduct.


title: Derived invariants for surface cut algebras II: the punctured case

abstract: For each algebra of global dimension 2 arising from the quiver with potential associated to a triangulation of an unpunctured surface, Amiot-Grimeland have defined an integer-valued function on the first singular homology group of the surface, and have proved that two such algebras of global dimension 2 are derived equivalent precisely when there exists an automorphism of the surface that makes their associated functions coincide. In the present paper we generalize the constructions and results of Amiot-Grimeland to the setting of arbitrarily punctured surfaces. As an application, we show that there always is a derived equivalence between any two algebras of global dimension 2 arising from the quivers with potential of (valency at least 2) triangulations of arbitrarily punctured polygons. While in the unpunctured case the quiver with potential of any triangulation admits cuts yielding algebras of global dimension at most 2, in the case of punctured surfaces the QPs of some triangulations do not admit cuts, and even when they do, the global dimension of the corresponding degree-0 algebra may exceed 2. In this paper we give a combinatorial characterization of each of these two situations.


title: Parallel SGD: When does averaging help?

abstract: Consider a number of workers running SGD independently on the same pool of data and averaging the models every once in a while -- a common but not well understood practice. We study model averaging as a variance-reducing mechanism and describe two ways in which the frequency of averaging affects convergence. For convex objectives, we show the benefit of frequent averaging depends on the gradient variance envelope. For non-convex objectives, we illustrate that this benefit depends on the presence of multiple globally optimal points. We complement our findings with multicore experiments on both synthetic and real data.


title: Origin of reversible photo-induced phase separation in hybrid perovskites

abstract: Nonequilibrium processes occurring in functional materials can significantly impact device efficiencies and are often difficult to characterize due to the broad range of length and time scales involved. In particular, mixed halide hybrid perovskites are promising for optoelectronics, yet the halides reversibly phase separate when photo-excited, significantly altering device performance. By combining nanoscale imaging and multiscale modeling, we elucidate the mechanism underlying this phenomenon, demonstrating that local strain induced by photo-generated polarons promotes halide phase separation and leads to nucleation of light-stabilized iodide-rich clusters. This effect relies on the unique electromechanical properties of hybrid materials, characteristic of neither their organic nor inorganic constituents alone. Exploiting photo-induced phase separation and other nonequilibrium phenomena in hybrid materials, generally, could enable new opportunities for expanding the functional applications in sensing, photoswitching, optical memory, and energy storage.


title: Thermodynamics of a Simple Three-Dimensional DNA Hairpin Model

abstract: We characterize the equation of state for a simple three-dimensional DNA hairpin model using a Metropolis Monte Carlo algorithm. This algorithm was run at constant temperature and fixed separation between the terminal ends of the strand. From the equation of state, we compute the compressibility, thermal expansion coefficient, and specific heat along with adiabatic path.


title: Transient cerebral hypoperfusion and hypertensive events during atrial fibrillation: a plausible mechanism for cognitive impairment

abstract: Atrial fibrillation (AF) is associated with an increased risk of dementia and cognitive decline, independent of strokes. Several mechanisms have been proposed to explain this association, but altered cerebral blood flow dynamics during AF has been poorly investigated: in particular, it is unknown how AF influences hemodynamic parameters of the distal cerebral circulation, at the arteriolar and capillary level. Two coupled lumped-parameter models (systemic and cerebrovascular circulations, respectively) were here used to simulate sinus rhythm (SR) and AF. For each simulation 5000 cardiac cycles were analyzed and cerebral hemodynamic parameters were calculated. With respect to SR, AF triggered a higher variability of the cerebral hemodynamic variables which increases proceeding towards the distal circulation, reaching the maximum extent at the arteriolar and capillary levels. This variability led to critical cerebral hemodynamic events of excessive pressure or reduced blood flow: 303 hypoperfusions occurred at the arteriolar level, while 387 hypertensive events occurred at the capillary level during AF. By contrast, neither hypoperfusions nor hypertensive events occurred during SR. Thus, the impact of AF per se on cerebral hemodynamics candidates as a relevant mechanism into the genesis of AF-related cognitive impairment/dementia.


title: Personalized Prognostic Models for Oncology: A Machine Learning Approach

abstract: We have applied a little-known data transformation to subsets of the Surveillance, Epidemiology, and End Results (SEER) publically available data of the National Cancer Institute (NCI) to make it suitable input to standard machine learning classifiers. This transformation properly treats the right-censored data in the SEER data and the resulting Random Forest and Multi-Layer Perceptron models predict full survival curves. Treating the 6, 12, and 60 months points of the resulting survival curves as 3 binary classifiers, the 18 resulting classifiers have AUC values ranging from .765 to .885. Further evidence that the models have generalized well from the training data is provided by the extremely high levels of agreement between the random forest and neural network models predictions on the 6, 12, and 60 month binary classifiers.


title: Almost balanced biased graph representations of frame matroids

abstract: Given a 3-connected biased graph $\Omega$ with a balancing vertex, and with frame matroid $F(\Omega)$ nongraphic and 3-connected, we determine all biased graphs $\Omega'$ with $F(\Omega') = F(\Omega)$. As a consequence, we show that if $M$ is a 4-connected nongraphic frame matroid represented by a biased graph $\Omega$ having a balancing vertex, then $\Omega$ essentially uniquely represents $M$. More precisely, all biased graphs representing $M$ are obtained from $\Omega$ by replacing a subset of the edges incident to its unique balancing vertex with unbalanced loops.


title: Model $I$--$V$ curves and figures of merit of underdamped deterministic Josephson ratchets

abstract: We propose simple models for the current-voltage characteristics of typical Josephson ratchets. We consider the case of a ratchet working against a constant applied counter force and derive analytical expressions for the key characteristics of such a ratchet: rectification curve, stopping force, input and output powers and rectification efficiency. Optimization of the ratchet performance is discussed.


title: Automatic Neuron Detection in Calcium Imaging Data Using Convolutional Networks

abstract: Calcium imaging is an important technique for monitoring the activity of thousands of neurons simultaneously. As calcium imaging datasets grow in size, automated detection of individual neurons is becoming important. Here we apply a supervised learning approach to this problem and show that convolutional networks can achieve near-human accuracy and superhuman speed. Accuracy is superior to the popular PCA/ICA method based on precision and recall relative to ground truth annotation by a human expert. These results suggest that convolutional networks are an efficient and flexible tool for the analysis of large-scale calcium imaging data.


title: VideoMCC: a New Benchmark for Video Comprehension

abstract: While there is overall agreement that future technology for organizing, browsing and searching videos hinges on the development of methods for high-level semantic understanding of video, so far no consensus has been reached on the best way to train and assess models for this task. Casting video understanding as a form of action or event categorization is problematic as it is not fully clear what the semantic classes or abstractions in this domain should be. Language has been exploited to sidestep the problem of defining video categories, by formulating video understanding as the task of captioning or description. However, language is highly complex, redundant and sometimes ambiguous. Many different captions may express the same semantic concept. To account for this ambiguity, quantitative evaluation of video description requires sophisticated metrics, whose performance scores are typically hard to interpret by humans. This paper provides four contributions to this problem. First, we formulate Video Multiple Choice Caption (VideoMCC) as a new well-defined task with an easy-to-interpret performance measure. Second, we describe a general semi-automatic procedure to create benchmarks for this task. Third, we publicly release a large-scale video benchmark created with an implementation of this procedure and we include a human study that assesses human performance on our dataset. Finally, we propose and test a varied collection of approaches on this benchmark for the purpose of gaining a better understanding of the new challenges posed by video comprehension.


title: Multi-Stage Temporal Difference Learning for 2048-like Games

abstract: Szubert and Jaskowski successfully used temporal difference (TD) learning together with n-tuple networks for playing the game 2048. However, we observed a phenomenon that the programs based on TD learning still hardly reach large tiles. In this paper, we propose multi-stage TD (MS-TD) learning, a kind of hierarchical reinforcement learning method, to effectively improve the performance for the rates of reaching large tiles, which are good metrics to analyze the strength of 2048 programs. Our experiments showed significant improvements over the one without using MS-TD learning. Namely, using 3-ply expectimax search, the program with MS-TD learning reached 32768-tiles with a rate of 18.31%, while the one with TD learning did not reach any. After further tuned, our 2048 program reached 32768-tiles with a rate of 31.75% in 10,000 games, and one among these games even reached a 65536-tile, which is the first ever reaching a 65536-tile to our knowledge. In addition, MS-TD learning method can be easily applied to other 2048-like games, such as Threes. Based on MS-TD learning, our experiments for Threes also demonstrated similar performance improvement, where the program with MS-TD learning reached 6144-tiles with a rate of 7.83%, while the one with TD learning only reached 0.45%.


title: A coupled $2\times2$D Babcock-Leighton solar dynamo model. II. Reference dynamo solutions

abstract: In this paper we complete the presentation of a new hybrid $2\times2$D flux transport dynamo (FTD) model of the solar cycle based on the Babcock-Leighton mechanism of poloidal magnetic field regeneration via the surface decay of bipolar magnetic regions (BMRs). This hybrid model is constructed by allowing the surface flux transport (SFT) simulation described in Lemerle et al. 2015 to provide the poloidal source term to an axisymmetric FTD simulation defined in a meridional plane, which in turn generates the BMRs required by the SFT. A key aspect of this coupling is the definition of an emergence function describing the probability of BMR emergence as a function of the spatial distribution of the internal axisymmetric magnetic field. We use a genetic algorithm to calibrate this function, together with other model parameters, against observed cycle 21 emergence data. We present a reference dynamo solution reproducing many solar cycle characteristics, including good hemispheric coupling, phase relationship between the surface dipole and the BMR-generating internal field, and correlation between dipole strength at cycle maximum and peak amplitude of the next cycle. The saturation of the cycle amplitude takes place through the quenching of the BMR tilt as a function of the internal field. The observed statistical scatter about the mean BMR tilt, built into the model, acts as a source of stochasticity which dominates amplitude fluctuations. The model thus can produce Dalton-like epochs of strongly suppressed cycle amplitude lasting a few cycles and can even shut off entirely following an unfavorable sequence of emergence events.


title: Flat structures on the deformations of Gepner chiral rings

abstract: We propose a simple method for the computation of the flat coordinates and Saito primitive forms on Frobenius manifolds of the deformations of Jacobi rings associated with isolated singularities. The method is based on using a conjecture about integral representations for the flat coordinates and on the Saito cohomology theory. This reduces the computation to a simple linear problem. We consider the case of the deformed Gepner chiral rings. The knowledge of the flat structures of Frobenius manifolds can be used for exact solution of the models of the topological conformal field theories corresponding to these chiral rings.


title: Split-sideband spectroscopy in slowly modulated optomechanics

abstract: Optomechanical coupling between the motion of a mechanical oscillator and a cavity represents a new arena for experimental investigation of quantum effects on the mesoscopic and macroscopic scale.The motional sidebands of the output of a cavity offer ultra-sensitive probes of the dynamics. We introduce a scheme whereby these sidebands split asymmetrically and show how they may be used as experimental diagnostics and signatures of quantum noise limited dynamics. We show split-sidebands with controllable asymmetry occur by simultaneously modulating the light-mechanical coupling $g$ and $\omega_M$ - slowly and out of-phase. Such modulations are generic but already occur in optically trapped set-ups where the equilibrium point of the oscillator is varied cyclically. We analyse recently observed, but overlooked, experimental split-sideband asymmetries; although not yet in the quantum regime, the data suggests that split sideband structures are easily accessible to future experiments.


title: Weibull Generalized Exponential Distribution

abstract: This paper introduces a new three-parameters model called the Weibull-G exponential distribution (WGED) distribution which exhibits bathtub-shaped hazard rate. Some of it's statistical properties are obtained including quantile, moments, generating functions, reliability and order statistics. The method of maximum likelihood is used for estimating the model parameters and the observed Fisher's information matrix is derived. We illustrate the usefulness of the proposed model by applications to real data.


title: Separately Radial and Radial Toeplitz Operators on the Projective Space and Representation Theory

abstract: We consider separately radial (with corresponding group $\mathbb{T}^n$) and radial (with corresponding group $\mathrm{U}(n))$ symbols on the projective space $\mathbb{P}^n(\mathbb{C})$, as well as the associated Toeplitz operators on the weighted Bergman spaces. It is known that the $C^*$-algebras generated by each family of such Toeplitz operators are commutative. We present a new representation theoretic proof of such commutativity. Our method is easier and more enlightening as it shows that the commutativity of the $C^*$-algebras is a consequence of the existence of multiplicity-free representations. Furthermore, our method shows how to extend the current formulas for the spectra of the corresponding Toeplitz operators to any closed group lying between $\mathbb{T}^n$ and $\mathrm{U}(n)$.


title: Variable-Sized Uncertainty and Inverse Problems in Robust Optimization

abstract: In robust optimization, the general aim is to find a solution that performs well over a set of possible parameter outcomes, the so-called uncertainty set. In this paper, we assume that the uncertainty size is not fixed, and instead aim at finding a set of robust solutions that covers all possible uncertainty set outcomes. We refer to these problems as robust optimization with variable-sized uncertainty. We discuss how to construct smallest possible sets of min-max robust solutions and give bounds on their size. A special case of this perspective is to analyze for which uncertainty sets a nominal solution ceases to be a robust solution, which amounts to an inverse robust optimization problem. We consider this problem with a min-max regret objective and present mixed-integer linear programming formulations that can be applied to construct suitable uncertainty sets. Results on both variable-sized uncertainty and inverse problems are further supported with experimental data.


title: Spread, volatility, and volume relationship in financial markets and market making profit optimization

abstract: We study the relationship between price spread, volatility and trading volume. We find that spread forms as a result of interplay between order liquidity and order impact. When trading volume is small adding more liquidity helps improve price accuracy and reduce spread, but after some point additional liquidity begins to deteriorate price. The model allows to connect the bid-ask spread and high-low bars to measurable microstructural parameters and express their dependence on trading volume, volatility and time horizon. Using the established relations, we address the operating spread optimization problem to maximize market-making profit.


title: Finite-time Partitions for Lagrangian Structure Identification in Gulf Stream Eddy Transport

abstract: We develop a methodology to identify finite-time Lagrangian structures from data and models using an extension of the Koopman operator-theoretic methods developed for velocity fields with simple (periodic, quasi-periodic) time-dependence. To achieve this, the notion of the Finite Time Ergodic (FiTER) partition is developed and rigorously justified. In combination with a clustering-based approach, the methodology enables identification of the temporal evolution of Lagrangian structures in a classic, benchmark, oceanographic transport problem, namely the cross-stream flux induced by the interaction of a meso- scale Gulf Stream Ring eddy with the main jet. We focus on a single mixing event driven by the interaction between an energetic cold core ring (a cyclone), the strong jet, and a number of smaller scale cyclones and anticyclones. The new methodology enab les reconstruction of Lagrangian structures in three dimensions and analysis of their time-evolution.


title: Network Infusion to Infer Information Sources in Networks

abstract: Several significant models have been developed that enable the study of diffusion of signals across biological, social and engineered networks. Within these established frameworks, the inverse problem of identifying the source of the propagated signal is challenging, owing to the numerous alternative possibilities for signal progression through the network. In real world networks, the challenge of determining sources is compounded as the true propagation dynamics are typically unknown, and when they have been directly measured, they rarely conform to the assumptions of any of the well-studied models. In this paper we introduce a method called Network Infusion (NI) that has been designed to circumvent these issues, making source inference practical for large, complex real world networks. The key idea is that to infer the source node in the network, full characterization of diffusion dynamics, in many cases, may not be necessary. This objective is achieved by creating a diffusion kernel that well-approximates standard diffusion models, but lends itself to inversion, by design, via likelihood maximization or error minimization. We apply NI for both single-source and multi-source diffusion, for both single-snapshot and multi-snapshot observations, and for both homogeneous and heterogeneous diffusion setups. We prove the mean-field optimality of NI for different scenarios, and demonstrate its effectiveness over several synthetic networks. Moreover, we apply NI to a real-data application, identifying news sources in the Digg social network, and demonstrate the effectiveness of NI compared to existing methods. Finally, we propose an integrative source inference framework that combines NI with a distance centrality-based method, which leads to a robust performance in cases where the underlying dynamics are unknown.


title: Robust Learning of Fixed-Structure Bayesian Networks

abstract: We investigate the problem of learning Bayesian networks in a robust model where an $\epsilon$-fraction of the samples are adversarially corrupted. In this work, we study the fully observable discrete case where the structure of the network is given. Even in this basic setting, previous learning algorithms either run in exponential time or lose dimension-dependent factors in their error guarantees. We provide the first computationally efficient robust learning algorithm for this problem with dimension-independent error guarantees. Our algorithm has near-optimal sample complexity, runs in polynomial time, and achieves error that scales nearly-linearly with the fraction of adversarially corrupted samples. Finally, we show on both synthetic and semi-synthetic data that our algorithm performs well in practice.


title: Nodes Localization in 3D Wireless Sensor Networks Based on Multidimensional Scaling Algorithm

abstract: In the recent years, there has been a huge advancement in wireless sensor computing technology. Today, wireless sensor network (WSN) has become a key technology for different types of smart environment. Nodes localization in WSN has arisen as a very challenging problem in the research community. Most of the applications for WSN are not useful without a priory known nodes positions. Adding GPS receivers to each node is an expensive solution and inapplicable for indoor environments. In this paper, we implemented and evaluated an algorithm based on multidimensional scaling (MDS) technique for three-dimensional (3D) nodes localization in WSN using improved heuristic method for distance calculation. Using extensive simulations we investigated our approach regarding various network parameters. We compared the results from the simulations with other approaches for 3D-WSN localization and showed that our approach outperforms other techniques in terms of accuracy.


title: Application of Wireless Sensor Networks for Indoor Temperature Regulation

abstract: Wireless sensor networks take a major part in our everyday lives by enhancing systems for home automation, healthcare, temperature control, energy consumption monitoring, and so forth. In this paper we focus on a system used for temperature regulation for residential, educational, industrial, and commercial premises, and so forth. We propose a framework for indoor temperature regulation and optimization using wireless sensor networks based on ZigBee platform. This paper considers architectural design of the system, as well as implementation guidelines. The proposed system favors methods that provide energy savings by reducing the amount of data transmissions through the network. Furthermore, the framework explores techniques for localization, such that the location of the nodes can be used by algorithms that regulate temperature settings.


title: Non-parametric estimation of morphological lopsidedness

abstract: Asymmetries in the neutral hydrogen gas distribution and kinematics of galaxies are thought to be indicators for both gas accretion and gas removal processes. These are of fundamental importance for galaxy formation and evolution. Upcoming large blind HI surveys will provide tens of thousands of galaxies for a study of these asymmetries in a proper statistical way. Due to the large number of expected sources and the limited resolution of the majority of objects, detailed modelling is not feasible for most detections. We need fast, automatic and sensitive methods to classify these objects in an objective way. Existing non-parametric methods suffer from effects like the dependence on signal to noise, resolution and inclination. Here we show how to correctly take these effects into account and show ways to estimate the precision of the methods. We will use existing and modelled data to give an outlook on the performance expected for galaxies observed in the various sky surveys planned for e.g. WSRT/APERTIF and ASKAP.


title: A taxonomy of localization techniques based on multidimensional scaling

abstract: Localization in Wireless Sensor Networks (WSNs) has been a challenging problem in the last decade. The most explored approaches for this purpose are based on multidimensional scaling (MDS) technique. The first algorithm that introduced MDS for nodes localization in sensor networks is well known as MDS-MAP. Since its appearance in 2003, many variations of MDS-MAP have been proposed in the literature. This paper aims to provide a comprehensive survey of the localization techniques that are based on MDS. We classify MDS-based algorithms according to different taxonomy features and different evaluation metrics.


title: Improved MDS-based Algorithm for Nodes Localization in Wireless Sensor Networks

abstract: With the recent development of technology, wireless sensor networks (WSN) are becoming an important part of many applications. Knowing the exact location of each sensor in the network is very important issue. Therefore, the localization problem is a growing field of interest. Adding GPS receivers to each sensor node is costly solution and inapplicable on nodes with limited resources. Additionally, it is not suitable for indoor environments. In this paper, we propose an algorithm for nodes localization in WNS based on multidimensional scaling (MDS) technique. Our approach improves MDS by distance matrix refinement. Using extensive simulations we investigated in details our approach regarding different network topologies, various network parameters and performance issues. The results from simulations show that our improved MDS (IMDS) algorithm outperforms well known MDS-MAP algorithm [1] in terms of accuracy.


title: Gravitational redshift from rotating body having intense magnetic field

abstract: It is well known fact that gravitational field can alter the space-time structure and gravitational redshift is its one example. Electromagnetic field can also alter the space-time similar to gravitational field. So electromagnetic field can give rise to an additional effect on gravitational redshift. There are many objects in nature, like neutron stars, magnetars etc which have high amount of rotation and magnetic field. In the present paper we will derive the expression of gravitational redshift from rotating body having intense magnetic field by using the action function of the electromagnetic fields.


title: Atmospheric neutrinos, $nu_e-nu_s$ oscillations, and a novel neutrino evolution equation

abstract: If a sterile neutrino nu_s with an eV-scale mass and a sizeable mixing to the electron neutrino exists, as indicated by the reactor and gallium neutrino anomalies, a strong resonance enhancement of nu_e-nu_s oscillations of atmospheric neutrinos should occur in the TeV energy range. At these energies neutrino flavour transitions in the 3+1 scheme depend on just one neutrino mass squared difference and are fully described within a 3-flavour oscillation framework. We demonstrate that the flavour transitions of atmospheric nu_e can actually be very accurately described in a 2-flavour framework, with neutrino flavour evolution governed by an inhomogeneous Schroedinger-like equation. Evolution equations of this type have not been previously considered in the theory of neutrino oscillations.


title: On the decidability of the $\Sigma_2$ theories of the arithmetic and hyperarithmetic degrees as uppersemilattices

abstract: We establish the decidability of the $\Sigma_2$ theory of both the arithmetic and hyperarithmetic degrees in the language of uppersemilattices i.e. the language with $\leq, 0$ and $\sqcup$. This is achieved by using Kumabe-Slaman forcing - along with other known results - to show that given finite uppersemilattices $\mathcal{M}$ and $\mathcal{N}$, where $\mathcal{M}$ is a subuppersemilattice of $\mathcal{N}$, then for both degree structures, every embedding of $\mathcal{M}$ into the structure extends to one of $\mathcal{N}$ iff $\mathcal{N}$ is an end-extension of $\mathcal{M}$.


title: Artificial intelligence applied to the automatic analysis of absorption spectra. Objective measurement of the fine structure constant

abstract: A new and automated method is presented for the analysis of high-resolution absorption spectra. Three established numerical methods are unified into one "artificial intelligence" process: a genetic algorithm (GVPFIT); non-linear least-squares with parameter constraints (VPFIT); and Bayesian Model Averaging (BMA). The method has broad application but here we apply it specifically to the problem of measuring the fine structure constant at high redshift. For this we need objectivity and reproducibility. GVPFIT is also motivated by the importance of obtaining a large statistical sample of measurements of $\Delta\alpha/\alpha$. Interactive analyses are both time consuming and complex and automation makes obtaining a large sample feasible. In contrast to previous methodologies, we use BMA to derive results using a large set of models and show that this procedure is more robust than a human picking a single preferred model since BMA avoids the systematic uncertainties associated with model choice. Numerical simulations provide stringent tests of the whole process and we show using both real and simulated spectra that the unified automated fitting procedure out-performs a human interactive analysis. The method should be invaluable in the context of future instrumentation like ESPRESSO on the VLT and indeed future ELTs. We apply the method to the $z_{abs} = 1.8389$ absorber towards the $z_{em} = 2.145$ quasar J110325-264515. The derived constraint of $\Delta\alpha/\alpha = 3.3 \pm 2.9 \times 10^{-6}$ is consistent with no variation and also consistent with the tentative spatial variation reported in Webb et al (2011) and King et al (2012).


title: Deliberation favours social efficiency by making people disregard their relative shares: Evidence from US and India

abstract: Groups make decisions on both the production and the distribution of resources. These decisions typically involve a tension between increasing the total level of group resources (i.e. social efficiency) and distributing these resources among group members (i.e. individuals' relative shares). This is the case because the redistribution process may destroy part of the resources, thus resulting in socially inefficient allocations. Here we apply a dual-process approach to understand the cognitive underpinnings of this fundamental tension. We conducted a set of experiments to examine the extent to which different allocation decisions respond to intuition or deliberation. In a newly developed approach, we assess intuition and deliberation at both the trait level (using the Cognitive Reflection Test, henceforth CRT) and the state level (through the experimental manipulation of response times). To test for robustness, experiments were conducted in two countries: the US and India. Despite absolute level differences across countries, in both locations we show that: (i) time pressure and low CRT scores are associated with individuals' concerns for their relative shares; (ii) time delay and high CRT scores are associated with individuals' concerns for social efficiency. These findings demonstrate that deliberation favours social efficiency by overriding individuals' intuitive tendency to focus on relative shares.


title: Syzygies over the Polytope Semiring

abstract: Tropical geometry and its applications indicate a "theory of syzygies" over polytope semirings. Taking cue from this indication, we study a notion of syzygies over the polytope semiring. We begin our exploration with the concept of Newton basis, an analogue of Gr\"obner basis that captures the image of an ideal under the Newton polytope map. The image ${\rm New}(I)$ of a graded ideal $I$ under the Newton polytope is a graded sub-semimodule of the polytope semiring. Analogous to the Hilbert series, we define the notion of Newton-Hilbert series that encodes the rank of each graded piece of ${\rm New}(I)$. We prove the rationality of the Newton-Hilbert series for sub-semimodules that satisfy a property analogous to Cohen-Macaulayness. We define notions of regular sequence of polytopes and syzygies of polytopes. We show an analogue of the Koszul property characterizing the syzygies of a regular sequence of polytopes.


title: Fast Multi-Layer Laplacian Enhancement

abstract: A novel, fast and practical way of enhancing images is introduced in this paper. Our approach builds on Laplacian operators of well-known edge-aware kernels, such as bilateral and nonlocal means, and extends these filter's capabilities to perform more effective and fast image smoothing, sharpening and tone manipulation. We propose an approximation of the Laplacian, which does not require normalization of the kernel weights. Multiple Laplacians of the affinity weights endow our method with progressive detail decomposition of the input image from fine to coarse scale. These image components are blended by a structure mask, which avoids noise/artifact magnification or detail loss in the output image. Contributions of the proposed method to existing image editing tools are: (1) Low computational and memory requirements, making it appropriate for mobile device implementations (e.g. as a finish step in a camera pipeline), (2) A range of filtering applications from detail enhancement to denoising with only a few control parameters, enabling the user to apply a combination of various (and even opposite) filtering effects.


title: Swift observations of the 2015 outburst of AG Peg -- from slow nova to classical symbiotic outburst

abstract: Symbiotic stars often contain white dwarfs with quasi-steady shell burning on their surfaces. However, in most symbiotics, the origin of this burning is unclear. In symbiotic slow novae, however, it is linked to a past thermonuclear runaway. In June 2015, the symbiotic slow nova AG Peg was seen in only its second optical outburst since 1850. This recent outburst was of much shorter duration and lower amplitude than the earlier eruption, and it contained multiple peaks -- like outbursts in classical symbiotic stars such as Z And. We report Swift X-ray and UV observations of AG Peg made between June 2015 and January 2016. The X-ray flux was markedly variable on a time scale of days, particularly during four days near optical maximum, when the X-rays became bright and soft. This strong X-ray variability continued for another month, after which the X-rays hardened as the optical flux declined. The UV flux was high throughout the outburst, consistent with quasi-steady shell burning on the white dwarf. Given that accretion disks around white dwarfs with shell burning do not generally produce detectable X-rays (due to Compton-cooling of the boundary layer), the X-rays probably originated via shocks in the ejecta. As the X-ray photo-electric absorption did not vary significantly, the X-ray variability may directly link to the properties of the shocked material. AG Peg's transition from a slow symbiotic nova (which drove the 1850 outburst) to a classical symbiotic star suggests that shell burning in at least some symbiotic stars is residual burning from prior novae.


title: Qualitative changes in spike-based neural coding and synchronization at the saddle-node loop bifurcation

abstract: Information processing in the brain crucially depends on encoding properties of single neurons, with particular relevance of the spike-generation mechanism. The latter hinges upon the bifurcation type at the transition point between resting state and limit cycle spiking. Prominent qualitative changes in encoding have previously been attributed to a specific switch of such a bifurcation at the Bogdanov-Takens (BT) point. This study unveils another, highly relevant and so far underestimated transition point: the saddle-node loop bifurcation. As we show, this bifurcation turns out to induce even more drastic changes in spike-based coding than the BT transition. This result arises from a direct effect of the saddle-node loop bifurcation on the limit cycle and hence spike dynamics, in contrast to the BT bifurcation, whose immediate influence is exerted upon the subthreshold dynamics and hence only indirectly relates to spiking. We specifically demonstrate that the saddle-node loop bifurcation (i) ubiquitously occurs in planar neuron models with a saddle-node on invariant cycle onset bifurcation, and (ii) results in a symmetry breaking of the system's phase-response curve. The latter entails close to optimal coding and synchronization properties in event-based information processing units, such as neurons. The saddle-node loop bifurcation leads to a peak in synchronization range and provides an attractive mechanism for the so far unresolved facilitation of high frequencies in neuronal processing. The derived bifurcation structure is of interest in any system for which a relaxation limit is admissible, such as Josephson junctions and chemical oscillators. On the experimental side, our theory applies to optical stimulation of nerve cells, and reveals that these techniques could manipulate a variety of information processing characteristics in nerve cells beyond pure activation.


title: jInv -- a flexible Julia package for PDE parameter estimation

abstract: Estimating parameters of Partial Differential Equations (PDEs) from noisy and indirect measurements often requires solving ill-posed inverse problems. These so called parameter estimation or inverse medium problems arise in a variety of applications such as geophysical, medical imaging, and nondestructive testing. Their solution is computationally intense since the underlying PDEs need to be solved numerous times until the reconstruction of the parameters is sufficiently accurate. Typically, the computational demand grows significantly when more measurements are available, which poses severe challenges to inversion algorithms as measurement devices become more powerful. In this paper we present jInv, a flexible framework and open source software that provides parallel algorithms for solving parameter estimation problems with many measurements. Being written in the expressive programming language Julia, jInv is portable, easy to understand and extend, cross-platform tested, and well-documented. It provides novel parallelization schemes that exploit the inherent structure of many parameter estimation problems and can be used to solve multiphysics inversion problems as is demonstrated using numerical experiments motivated by geophysical imaging.


title: Rashba-Dresselhaus spin-splitting in the bulk ferroelectric oxide BiAlO$_3$

abstract: It has been recently suggested that the coexistence of ferroelectricity and Rashba-like spin-splitting effects due to spin-orbit coupling in a single material may allow for a non-volatile electric control of spin degrees of freedom. In the present work, we compared the structural and ferroelectric properties of tetragonal and rhombohedral phases of ferroelectric BiAlO$_3$ by means of density-functional calculations. In both phases, we carefully investigated Rashba and Dresselhaus effects, giving rise to spin-splitting in their bulk electronic structure, particularly near the conduction band minimum, supplementing our first-principles results with an effective $ k\cdot p$ model analysis. The full reversal of the spin texture with ferroelectric polarization switching was also predicted. BiAlO$_3$ can therefore be considered as the first known oxide to exhibit a coexistence of ferroelectricity and Rashba-Dresselhaus effects.


title: On specification and measure expansiveness

abstract: We relate the local specification and periodic shadowing properties. We also clarify the relation between local weak specification and local specification if the system is measure expansive. The notion of strong measure expansiveness is introduced, and an example of a non-expansive systems with the strong measure expansive property is given. Moreover, we find a family of examples with the $N$-expansive property, which are not strong measure expansive. We finally show a spectral decomposition theorem for strong measure expansive dynamical systems with shadowing.


title: Dispersing artifacts in FT-STS: a comparison of set point effects across acquisition modes

abstract: Fourier transform scanning tunneling spectroscopy (FT-STS), or quasiparticle interference (QPI), has become an influential tool for the study of a wide range of important materials in condensed matter physics. However, FT-STS in complex materials is often challenging to interpret, requiring significant theoretical input in many cases, making it crucial to understand potential artifacts of the measurement. Here, we compare the most common modes of acquiring FT-STS data and show through both experiment and simulations that artifact features can arise that depend on how the tip height is stabilized throughout the course of the measurement. The most dramatic effect occurs when a series of dI/dV maps at different energies are acquired with simultaneous constant current feedback; here a feature that disperses in energy appears that is not observed in other measurement modes. Such artifact features are similar to those arising from real physical processes in the sample and are susceptible to misinterpretation.


title: Do alien particles exist, and can they be detected?

abstract: We may call "alien particles" those particles belonging to the matter/field content of a $d$-dimensional brane other than the $3$-brane (or stack of branes) sweeping the space-time in which we live. They can appear in our space-time at the regions of intersection between our and their brane. They can be identified (or not) as alien matter depending on their properties, on the physical laws governing their evolution in the "homeland" brane, and on the details of our detection techniques.


title: Lasing in dark and bright modes of a finite-sized plasmonic lattice

abstract: Lasing at the nanometre scale promises strong light-matter interactions and ultrafast operation. Plasmonic resonances supported by metallic nanoparticles have extremely small mode volumes and high field enhancements, making them an ideal platform for studying nanoscale lasing. At visible frequencies, however, the applicability of plasmon resonances is limited due to strong ohmic and radiative losses. Intriguingly, plasmonic nanoparticle arrays support non-radiative dark modes that offer longer life-times but are inaccessible to far-field radiation. Here, we show lasing both in dark and bright modes of an array of silver nanoparticles combined with optically pumped dye molecules. Linewidths of 0.2 nanometers at visible wavelengths and room temperature are observed. Access to the dark modes is provided by a coherent out-coupling mechanism based on the finite size of the array. The results open a route to utilize all modes of plasmonic lattices, also the high-Q ones, for studies of strong light-matter interactions, condensation and photon fluids.


title: The Formation of Bulges, Discs and Two Component Galaxies in the CANDELS Survey at z < 3

abstract: We examine a sample of 1495 galaxies in the CANDELS fields to determine the evolution of two component galaxies, including bulges and discs, within massive galaxies at the epoch 1 < z < 3 when the Hubble sequence forms. We fit all of our galaxies' light profiles with a single S\'ersic fit, as well as with a combination of exponential and S\'ersic profiles. The latter is done in order to describe a galaxy with an inner and an outer component, or bulge and disc component. We develop and use three classification methods (visual, F-test and the RFF) to separate our sample into 1-component galaxies (disc/spheroids-like galaxies) and 2-component galaxies (galaxies formed by an 'inner part' or bulge and an 'outer part' or disc). We then compare the results from using these three different ways to classify our galaxies. We find that the fraction of galaxies selected as 2-component galaxies increases on average 50 per cent from the lowest mass bin to the most massive galaxies, and decreases with redshift by a factor of four from z = 1 to z = 3. We find that single S\'ersic 'disc-like' galaxies have the highest relative number densities at all redshifts, and that 2-component galaxies have the greatest increase and become at par with S\'ersic discs by z = 1. We also find that the systems we classify as 2-component galaxies have an increase in the sizes of their outer components, or 'discs', by about a factor of three from z = 3 to z = 1.5, while the inner components or 'bulges' stay roughly the same size. This suggests that these systems are growing from the inside out, whilst the bulges or protobulges are in place early in the history of these galaxies. This is also seen to a lesser degree in the growth of single 'disc-like' galaxies vs. 'spheroid-like' galaxies over the same epoch.


title: Polymer-Fourier quantization of the scalar field revisited

abstract: The Polymer Quantization of the Fourier modes of the real scalar field is studied within algebraic scheme. We replace the positive linear functional of the standard Poincar\'e invariant quantization by a singular one. This singular positive linear functional is constructed as mimicking the singular limit of the complex structure of the Poincar\'e invariant Fock quantization. The resulting symmetry group of such Polymer Quantization is the subgroup $\mbox{SDiff}(\mathbb{R}^4)$ which is a subgroup of $\mbox{Diff}(\mathbb{R}^4)$ formed by spatial volume preserving diffeomorphisms. In consequence, this yields an entirely different irreducible representation of the Canonical Commutation Relations, non-unitary equivalent to the standard Fock representation. We also compared the Poincar\'e invariant Fock vacuum with the Polymer Fourier vacuum.


title: High-Dimensional Sparse Fourier Algorithms

abstract: In this paper, we discuss the development of a sublinear sparse Fourier algorithm for high-dimensional data. In ``Adaptive Sublinear Time Fourier Algorithm" by D. Lawlor, Y. Wang and A. Christlieb (2013), an efficient algorithm with $\Theta(k\log k)$ average-case runtime and $\Theta(k)$ average-case sampling complexity for the one-dimensional sparse FFT was developed for signals of bandwidth $N$, where $k$ is the number of significant modes such that $k\ll N$. In this work we develop an efficient algorithm for sparse FFT for higher dimensional signals, extending some of the ideas in the paper mentioned above. Note a higher dimensional signal can always be unwrapped into a one dimensional signal, but when the dimension gets large, unwrapping a higher dimensional signal into a one dimensional array is far too expensive to be realistic. Our approach here introduces two new concepts: ``partial unwrapping'' and ``tilting''. These two ideas allow us to efficiently compute the sparse FFT of higher dimensional signals.


title: Birch's theorem: if f(n) is multiplicative and has a non-decreasing normal order then f(n)=n^{alpha}

abstract: For pedagogical purposes (inclusion in lecture notes) we review the proof of the theorem stated in the title. At the end we state a problem.


title: Raman tensor elements of $\beta\text{-Ga}_2\text{O}_3$

abstract: The Raman spectrum and particularly the Raman scattering intensities of monoclinic $\beta\text{-Ga}_2\text{O}_3$ are investigated by experiment and theory. The low symmetry of $\beta\text{-Ga}_2\text{O}_3$ results in a complex dependence of the Raman intensity for the individual phonon modes on the scattering geometry which is additionally affected by birefringence. We measured the Raman spectra in dependence on the polarization direction for backscattering on three crystallographic planes of $\beta\text{-Ga}_2\text{O}_3$ and modeled these dependencies using a modified Raman tensor formalism which takes birefringence into account. The spectral position of all 15 Raman-active phonon modes and the Raman tensor elements of 13 modes were determined and are compared to results from ab-initio calculations.


title: Full Optical Control of Topological Transitions in Polariton Chern Insulator Analog

abstract: Exciton-polariton lattices allow to implement topologically protected photonic edge states at optical frequencies. Taking advantage of the interacting character of polaritons, we show that several topological phases, belonging to the Quantum Anomalous Hall family, can occur in polariton graphene with the TE-TM spin-orbit coupling and an effective Zeeman field. This Zeeman field can be entirely controlled by the non-linearity of the macro-occupied mode created by circularly polarized resonant pumping. The topological phases found are distinct from the ones of normal graphene in presence of a combination of Rashba and Zeeman field.


title: Fragmentation of a Jet with Small Radius

abstract: In this paper we consider the fragmentation of a parton into a jet with small jet radius $R$. Perturbatively, logarithms of $R$ can appear, which for narrow jets can lead to large corrections. Using soft-collinear effective theory, we introduce the fragmentation function to a jet (FFJ), which describes the fragmentation of a parton into a jet. We discuss how these objects are related to the standard jet functions. Calculating the FFJ to next-to-leading order, we show that these objects satisfy the standard Dokshitzer-Gribov-Lipatov-Altarelli-Parisi evolution equations, with a natural scale that depends upon $R$. By using the standard renormalization group evolution, we can therefore resum logarithms of $R$. We further use the soft-collinear effective theory to prove a factorization theorem where the FFJs naturally appear, for the fragmentation of a hadron within a jet with small $R$. Finally, we also show how this formalism can be used to resum the ratio of jet radii for a subjet to be emitted from within a fat jet.


title: Magnetic Skyrmions on a Two-Lane Racetrack

abstract: Magnetic skyrmions are particle-like textures in the magnetization, characterized by a topological winding number. Nanometer-scale skyrmions have been observed at room temperature in magnetic multilayer structures. The combination of small size, topological quantization, and their efficient electric manipulation makes them interesting candidates for information carriers in high-performance memory devices. A skyrmion racetrack memory has been suggested where information is encoded in the distance between skyrmions moving in a one-dimensional nanostrip. Here, we propose an alternative design where skyrmions move in two (or more) parallel lanes and the information is stored in the lane number of each skyrmion. Such a multilane track can be constructed by controlling the height profile of the nanostrip. Repulsive skyrmion-skyrmion interactions in narrow nanostrips guarantee that skyrmions on different lanes cannot pass each other. Current pulses can be used to induce a lane change. Combining these elements provides a robust, efficient design of skyrmion-based storage devices.


title: Parallelizing quantum circuit synthesis

abstract: Quantum circuit synthesis is the process in which an arbitrary unitary operation is decomposed into a sequence of gates from a universal set, typically one which a quantum computer can implement both efficiently and fault-tolerantly. As physical implementations of quantum computers improve, the need is growing for tools which can effectively synthesize components of the circuits and algorithms they will run. Existing algorithms for exact, multi-qubit circuit synthesis scale exponentially in the number of qubits and circuit depth, leaving synthesis intractable for circuits on more than a handful of qubits. Even modest improvements in circuit synthesis procedures may lead to significant advances, pushing forward the boundaries of not only the size of solvable circuit synthesis problems, but also in what can be realized physically as a result of having more efficient circuits. We present a method for quantum circuit synthesis using deterministic walks. Also termed pseudorandom walks, these are walks in which once a starting point is chosen, its path is completely determined. We apply our method to construct a parallel framework for circuit synthesis, and implement one such version performing optimal $T$-count synthesis over the Clifford+$T$ gate set. We use our software to present examples where parallelization offers a significant speedup on the runtime, as well as directly confirm that the 4-qubit 1-bit full adder has optimal $T$-count 7 and $T$-depth 3.


title: Multiplierless 16-point DCT Approximation for Low-complexity Image and Video Coding

abstract: An orthogonal 16-point approximate discrete cosine transform (DCT) is introduced. The proposed transform requires neither multiplications nor bit-shifting operations. A fast algorithm based on matrix factorization is introduced, requiring only 44 additions---the lowest arithmetic cost in literature. To assess the introduced transform, computational complexity, similarity with the exact DCT, and coding performance measures are computed. Classical and state-of-the-art 16-point low-complexity transforms were used in a comparative analysis. In the context of image compression, the proposed approximation was evaluated via PSNR and SSIM measurements, attaining the best cost-benefit ratio among the competitors. For video encoding, the proposed approximation was embedded into a HEVC reference software for direct comparison with the original HEVC standard. Physically realized and tested using FPGA hardware, the proposed transform showed 35% and 37% improvements of area-time and area-time-squared VLSI metrics when compared to the best competing transform in the literature.


title: Find your Way by Observing the Sun and Other Semantic Cues

abstract: In this paper we present a robust, efficient and affordable approach to self-localization which does not require neither GPS nor knowledge about the appearance of the world. Towards this goal, we utilize freely available cartographic maps and derive a probabilistic model that exploits semantic cues in the form of sun direction, presence of an intersection, road type, speed limit as well as the ego-car trajectory in order to produce very reliable localization results. Our experimental evaluation shows that our approach can localize much faster (in terms of driving time) with less computation and more robustly than competing approaches, which ignore semantic information.


title: Noncommutative Brownian motion

abstract: We investigate the classical Brownian motion of a particle in a two-dimensional noncommutative (NC) space. Using the standard NC algebra embodied by the sympletic Weyl-Moyal formalism we find that noncommutativity induces a non-vanishing correlation between both coordinates at different times. The effect stands out as a signature of spatial noncommutativity and thus could offer a way to experimentally detect the phenomena. We further discuss some limiting scenarios and the trade-off between the scale imposed by the NC structure and the parameters of the Brownian motion itself.


title: Axial acoustic radiation force on rigid oblate and prolate spheroids in Bessel vortex beams of progressive, standing and quasi-standing waves

abstract: The analysis using the partial-wave series expansion (PWSE) method in spherical coordinates is extended to evaluate the acoustic radiation force experienced by rigid oblate and prolate spheroids centered on the axis of wave propagation of high-order Bessel vortex beams composed of progressive, standing and quasi-standing waves, respectively. A coupled system of linear equations is derived after applying the Neumann boundary condition for an immovable surface in a non-viscous fluid, and solved numerically by matrix inversion after performing a single numerical integration procedure. The system of linear equations depends on the partial-wave index n and the order of the Bessel vortex beam m using truncated but converging PWSEs in the least-squares sense. Numerical results for the radiation force function, which is the radiation force per unit energy density and unit cross-sectional surface, are computed with particular emphasis on the amplitude ratio describing the transition from the progressive to the pure standing waves cases, the aspect ratio (i.e., the ratio of the major axis over the minor axis of the spheroid), the half-cone angle and order of the Bessel vortex beam, as well as the dimensionless size parameter. A generalized expression for the radiation force function is derived for cases encompassing the progressive, standing and quasi-standing waves of Bessel vortex beams. This expression can be reduced to other types of beams/waves such as the zeroth-order Bessel non-vortex beam or the infinite plane wave case by appropriate selection of the beam parameters. The results for progressive waves reveal a tractor beam behavior, characterized by the emergence of an attractive pulling force acting in opposite direction of wave propagation. Moreover, the transition to the quasi-standing and pure standing wave cases shows the acoustical tweezers behavior in dual-beam Bessel vortex beams.


title: Priority-based Riemann solver for traffic flow on networks

abstract: In this article we introduce a new Riemann solver for traffic flow on networks. The Priority Riemann solver (PRS) provides a solution at junctions by taking into consideration priorities for the incoming roads and maximization of through flux. We prove existence of solutions for the solver for junctions with up to two incoming and two outgoing roads and show numerically the comparison with previous Riemann solvers. Additionally, we introduce a second version of the solver that considers the priorities as softer constraints and illustrate numerically the differences between the two solvers.


title: Learning to Poke by Poking: Experiential Learning of Intuitive Physics

abstract: We investigate an experiential learning paradigm for acquiring an internal model of intuitive physics. Our model is evaluated on a real-world robotic manipulation task that requires displacing objects to target locations by poking. The robot gathered over 400 hours of experience by executing more than 100K pokes on different objects. We propose a novel approach based on deep neural networks for modeling the dynamics of robot's interactions directly from images, by jointly estimating forward and inverse models of dynamics. The inverse model objective provides supervision to construct informative visual features, which the forward model can then predict and in turn regularize the feature space for the inverse model. The interplay between these two objectives creates useful, accurate models that can then be used for multi-step decision making. This formulation has the additional benefit that it is possible to learn forward models in an abstract feature space and thus alleviate the need of predicting pixels. Our experiments show that this joint modeling approach outperforms alternative methods.


title: On Calabi-Yau generalized complete intersections from Hirzebruch varieties and novel K3-fibrations

abstract: We consider the construction of Calabi-Yau varieties recently generalized to where the defining equations may have negative degrees over some projective space factors in the embedding space. Within such "generalized complete intersection" Calabi-Yau ("gCICY") three-folds, we find several sequences of distinct manifolds. These include both novel elliptic and K3-fibrations and involve Hirzebruch surfaces and their higher dimensional analogues. En route, we generalize the standard techniques of cohomology computation to these generalized complete intersection Calabi-Yau varieties.


title: The Albedo Distribution of Near Earth Asteroids

abstract: The cryogenic WISE mission in 2010 was extremely sensitive to asteroids and not biased against detecting dark objects. The albedos of 428 Near Earth Asteroids (NEAs) observed by WISE during its fully cryogenic mission can be fit quite well by a 3 parameter function that is the sum of two Rayleigh distributions. The Rayleigh distribution is zero for negative values, and follows $f(x) = x \exp[-x^2/(2\sigma^2)]/\sigma^2$ for positive x. The peak value is at x=\sigma, so the position and width are tied together. The three parameters are the fraction of the objects in the dark population, the position of the dark peak, and the position of the brighter peak. We find that 25.3% of the NEAs observed by WISE are in a very dark population peaking at $p_V = 0.03$, while the other 74.7% of the NEAs seen by WISE are in a moderately dark population peaking at $p_V = 0.168$. A consequence of this bimodal distribution is that the Congressional mandate to find 90% of all NEAs larger than 140 m diameter cannot be satisfied by surveying to H=22 mag, since a 140 m diameter asteroid at the very dark peak has H=23.7 mag, and more than 10% of NEAs are darker than p_V = 0.03.


title: Joint product numerical range and geometry of reduced density matrices

abstract: The reduced density matrices of a many-body quantum system form a convex set, whose three-dimensional projection $\Theta$ is convex in $\mathbb{R}^3$. The boundary $\partial\Theta$ of $\Theta$ may exhibit nontrivial geometry, in particular ruled surfaces. Two physical mechanisms are known for the origins of ruled surfaces: symmetry breaking and gapless. In this work, we study the emergence of ruled surfaces for systems with local Hamiltonians in infinite spatial dimension, where the reduced density matrices are known to be separable as a consequence of the quantum de Finetti's theorem. This allows us to identify the reduced density matrix geometry with joint product numerical range $\Pi$ of the Hamiltonian interaction terms. We focus on the case where the interaction terms have certain structures, such that ruled surface emerge naturally when taking a convex hull of $\Pi$. We show that, a ruled surface on $\partial\Theta$ sitting in $\Pi$ has a gapless origin, otherwise it has a symmetry breaking origin. As an example, we demonstrate that a famous ruled surface, known as the oloid, is a possible shape of $\Theta$, with two boundary pieces of symmetry breaking origin separated by two gapless lines.


title: Aharonov-Bohm effect in a helical ring with long-range hopping: Effects of Rashba spin-orbit interaction and disorder

abstract: We study Aharonov-Bohm effect in a two-terminal helical ring with long-range hopping in presence of Rashba spin-orbit interaction. We explore how the spin polarization behavior changes depending on the applied magnetic flux and the incoming electron energy. The most interesting feature that we articulate in this system is that zero-energy crossings appear in the energy spectra at $\Phi=0$ and also at integer multiples of half-flux quantum values ($n \Phi_0/2$, $n$ being an integer) of the applied magnetic flux. We investigate the transport properties of the ring using Green's function formalism and find that the zero energy transmission peaks corresponding to those zero energy crossings vanish in presence of Rashba spin-orbit interaction. We also incorporate static random disorder in our system and show that the zero energy crossings and transmission peaks are not immune to disorder even in absence of Rashba spin-orbit interaction. The latter prevents the possibility of behaving these helical states in the ring like topological insulator edge states.


title: Enhanced production of multi-strange hadrons in high-multiplicity proton-proton collisions

abstract: At sufficiently high temperature and energy density, nuclear matter undergoes a transition to a phase in which quarks and gluons are not confined: the Quark-Gluon Plasma (QGP) [1]. Such an extreme state of strongly-interacting QCD (Quantum Chromo-Dynamics) matter is produced in the laboratory with high-energy collisions of heavy nuclei, where an enhanced production of strange hadrons is observed [2-6]. Strangeness enhancement, originally proposed as a signature of QGP formation in nuclear collisions [7], is more pronounced for multi-strange baryons. Several effects typical of heavy-ion phenomenology have been observed in high-multiplicity proton-proton (pp) collisions [8,9]. Yet, enhanced production of multi-strange particles has not been reported so far. Here we present the first observation of strangeness enhancement in high-multiplicity pp collisions. We find that the integrated yields of strange and multi-strange particles relative to pions increases significantly with the event charged-particle multiplicity. The measurements are in remarkable agreement with p-Pb collision results [10,11] indicating that the phenomenon is related to the final system created in the collision. In high-multiplicity events strangeness production reaches values similar to those observed in Pb-Pb collisions, where a QGP is formed.


title: Generalized Preconditioning and Network Flow Problems

abstract: We consider approximation algorithms for the problem of finding $x$ of minimal norm $\|x\|$ satisfying a linear system $\mathbf{A} x = \mathbf{b}$, where the norm $\|\cdot \|$ is arbitrary and generally non-Euclidean. We show a simple general technique for composing solvers, converting iterative solvers with residual error $\|\mathbf{A} x - \mathbf{b}\| \leq t^{-\Omega(1)}$ into solvers with residual error $\exp(-\Omega(t))$, at the cost of an increase in $\|x\|$, by recursively invoking the solver on the residual problem $\tilde{\mathbf{b}} = \mathbf{b} - \mathbf{A} x$. Convergence of the composed solvers depends strongly on a generalization of the classical condition number to general norms, reducing the task of designing algorithms for many such problems to that of designing a \emph{generalized preconditioner} for $\mathbf{A}$. The new ideas significantly generalize those introduced by the author's earlier work on maximum flow, making them more widely applicable. As an application of the new technique, we present a nearly-linear time approximation algorithm for uncapacitated minimum-cost flow on undirected graphs. Given an undirected graph with $m$ edges labelled with costs, and $n$ vertices labelled with demands, the algorithm takes $\epsilon^{-2}m^{1+o(1)}$-time and outputs a flow routing the demands with total cost at most $(1+\epsilon)$ times larger than minimal, along with a dual solution proving near-optimality. The generalized preconditioner is obtained by embedding the cost metric into $\ell_1$, and then considering a simple hierarchical routing scheme in $\ell_1$ where demands initially supported on a dense lattice are pulled from a sparser lattice by randomly rounding unaligned coordinates to their aligned neighbors. Analysis of the generalized condition number for the preconditioner follows that of the classical multigrid algorithm for lattice Laplacian systems.


title: On the Poisson relation for compact Lie groups

abstract: Intuition drawn from quantum mechanics and geometric optics raises the following long-standing question: can the length spectrum of a closed Riemannian manifold be recovered from its Laplace spectrum? The Poisson relation states that for any closed Riemannian manifold $(M,g)$ the singular support of the trace of its wave group---a spectrally determined tempered distribution---is contained in the set consisting of $\pm \tau$, where $\tau$ is the length of a smoothly closed geodesic in $(M,g)$. Therefore, in cases where the Poisson relation is an equality, we obtain a method for retrieving the length spectrum of a manifold from its Laplace spectrum. The Poisson relation is known to be an equality for sufficiently "bumpy" Riemannian manifolds and there are no known counterexamples. We demonstrate that the Poisson relation is an equality for a compact Lie group equipped with a generic bi-invariant metric. Consequently, the length spectrum of a generic bi-invariant metric (and the rank of its underlying Lie group) can be recovered from its Laplace spectrum. Furthermore, we exhibit a substantial collection $\mathscr{G}$ of compact Lie groups---including those that are either tori, simple, simply-connected, or products thereof---with the property that for each group $U \in \mathscr{G}$ the Laplace spectrum of any bi-invariant metric $g$ carried by $U$ encodes the length spectrum of $g$ and the rank of $U$. The preceding statements are special cases of results concerning compact globally symmetric spaces for which the semi-simple part of the universal cover is split-rank. The manifolds considered herein join a short list of families of non-"bumpy" Riemannian manifolds for which the Poisson relation is known to be an equality.


title: Special Values of Motivic $L$-Functions and Zeta-Polynomials for Symmetric Powers of Elliptic Curves

abstract: Let $\mathcal{M}$ be a pure motive over $\mathbb{Q}$ of odd weight $w\geq 3$, even rank $d\geq 2$, and global conductor $N$ whose $L$-function $L(s,\mathcal{M})$ coincides with the $L$-function of a self-dual algebraic tempered cuspidal symplectic representation of $\mathrm{GL}_{d}(\mathbb{A}_{\mathbb{Q}})$. We show that a certain polynomial which generates special values of $L(s,\mathcal{M})$ (including all of the critical values) has all of its zeros equidistributed on the unit circle, provided that $N$ or $w$ are sufficiently large with respect to $d$. These special values have arithmetic significance in the context of the Bloch-Kato conjecture. We focus on applications to symmetric powers of semistable elliptic curves over $\mathbb{Q}$. Using the Rodriguez-Villegas transform, we use these results to construct large classes of "zeta-polynomials" (in the sense of Manin) arising from symmetric powers of semistable elliptic curves; these polynomials have a functional equation relating $s\mapsto 1-s$, and all of their zeros on the line $\operatorname{Re}(s)=1/2$.


title: An integral equation formulation for rigid bodies in Stokes flow in three dimensions

abstract: We present a new derivation of a boundary integral equation (BIE) for simulating the three-dimensional dynamics of arbitrarily-shaped rigid particles of genus zero immersed in a Stokes fluid, on which are prescribed forces and torques. Our method is based on a single-layer representation and leads to a simple second-kind integral equation. It avoids the use of auxiliary sources within each particle that play a role in some classical formulations. We use a spectrally accurate quadrature scheme to evaluate the corresponding layer potentials, so that only a small number of spatial discretization points per particle are required. The resulting discrete sums are computed in $\mathcal{O}(n)$ time, where $n$ denotes the number of particles, using the fast multipole method (FMM). The particle positions and orientations are updated by a high-order time-stepping scheme. We illustrate the accuracy, conditioning and scaling of our solvers with several numerical examples.


title: Constructing and proving the ground state of a generalized Ising model by the cluster tree optimization algorithm

abstract: Generalized Ising models, also known as cluster expansions, are an important tool in many areas of condensed-matter physics and materials science, as they are often used in the study of lattice thermodynamics, solid-solid phase transitions, magnetic and thermal properties of solids, and fluid mechanics. However, the problem of finding the global ground state of generalized Ising model has remained unresolved, with only a limited number of results for simple systems known. We propose a method to efficiently find the periodic ground state of a generalized Ising model of arbitrary complexity by a new algorithm which we term cluster tree optimization. Importantly, we are able to show that even in the case of an aperiodic ground state, our algorithm produces a sequence of states with energy converging to the true ground state energy, with a provable bound on error. Compared to the current state-of-the-art polytope method, this algorithm eliminates the necessity of introducing an exponential number of variables to counter frustration, and thus significantly improves tractability. We believe that the cluster tree algorithm offers an intuitive and efficient approach to finding and proving ground states of generalized Ising Hamiltonians of arbitrary complexity, which will help validate assumptions regarding local vs. global optimality in lattice models, as well as offer insights into the low-energy behavior of highly frustrated systems.


title: Coupled-Beam and Coupled-Bunch Instabilities

abstract: A problem of coupled-beam instability is solved for two multibunch beams with slightly different revolution frequencies, as in the Fermilab Recycler Ring (RR). Sharing of the inter-bunch growth rates between the intra-bunch modes is described. The general analysis is applied to the RR; possibilities to stabilize the beams by means of chromaticity and feedback are considered.


title: Cluster dynamics largely shapes protoplanetary disc sizes

abstract: It is still on open question to what degree the cluster environment influences the sizes of protoplanetary discs surrounding young stars. Particularly so for the short-lived clusters typical for the solar neighbourhood in which the stellar density and therefore the influence of the cluster environment changes considerably over the first 10 Myr. In previous studies often the effect of the gas on the cluster dynamics has been neglected, this is remedied here. Using the code NBody6++ we study the stellar dynamics in different developmental phases - embedded, expulsion, expansion - including the gas and quantify the effect of fly-bys on the disc size. We concentrate on massive clusters ($M_{\text{cl}} \geq 10^3 - 6 \cdot 10^4 M_{\text{Sun}}$), which are representative for clusters like the Orion Nebula Cluster (ONC) or NGC 6611. We find that not only the stellar density but also the duration of the embedded phase matters. The densest clusters react fastest to the gas expulsion and drop quickly in density, here 98% of relevant encounters happen before gas expulsion. By contrast, discs in sparser clusters are initially less affected but as they expand slower 13% of discs are truncated after gas expulsion. For ONC-like clusters we find that usually discs larger than 500 AU are affected by the environment, which corresponds to the observation that 200 AU-sized discs are common. For NGC 6611-like clusters disc sizes are cut-down on average to roughly 100 AU. A testable hypothesis would be that the discs in the centre of NGC 6611 should be on average ~20 AU and therefore considerably smaller than in the ONC.


title: Ultracold neutron detection with 6Li-doped glass scintillators, NANOSC: a fast ultracold neutron detector for the nEDM experiment at the Paul Scherrer Institute

abstract: This paper summarizes the results from measurements aiming to characterize ultracold neutron detection with 6Li-doped glass scintillators. Single GS10 or GS20 scintillators, with a thickness of 100-200 micrometer, fulfill the ultracold neutron detection requirements with an acceptable neutron-gamma discrimination. This discrimination is clearly improved with a stack of two scintillators: a 6Li-depleted glass bonded to a 6Li-enriched glass. The optical contact bonding is used between the scintillators in order to obtain a perfect optical contact. The scintillator's detection efficiency is similar to that of a 3He Strelkov gas detector. Coupled to a digital data acquisition system, counting rates up to a few 10^5 counts/s can be handled. A detector based on such a scintillator stack arrangement was built and has been used in the neutron electric dipole moment experiment at the Paul Scherrer Institute since 2010. Its response for the regular runs of the neutron electric dipole moment experiment is presented.


title: Deep XMM-Newton Observations of the NW Radio Relic Region of Abell 3667

abstract: The results of long XMM-Newton X-ray observations of the NW radio relic of Abell 3667 are presented. A shock is detected at the sharp outer edge of the radio relic, both in the X-ray surface brightness and the temperature profiles. The Mach number is M = 2.54^+0.80_-0.43. The temperature jump at the shock is larger than expected from the density jump, which may indicate that a dynamically important magnetic field aligned primarily parallel to the shock front is present. The gas temperature rises gradually over several arc minutes within the shock region. This could indicate that the shock energy is initially dissipated into some mix of thermal and nonthermal (e.g., turbulence) components, and that the nonthermal energy decays into heat in the post-shock region. The observed radio relic can be powered if ~0.2% of the energy dissipated in the shock goes into the (re)acceleration of relativistic electrons. We show that the observed steepening of the radio spectrum with distance behind the shock is consistent with radiative losses by the radio-emitting electrons. However, the radio spectrum immediately behind the shock is flatter than expected for linear diffusive shock acceleration of thermal electrons. This suggests that the shock re-accelerates a pre-existing population of relativistic electrons. We also detect a bright, cool region (the "Mushroom") to the south of the radio relic, which we propose is the remnant cool core of a merging subcluster, and that this subcluster was the driver for the observed NW shock. In this model, the properties of Abell 3667 are mainly the result of an offset binary merger, and the cluster is being observed about 1 Gyr after first core passage. We predict that deeper X-ray or SZ observations of the SE radio relic will reveal a second merger shock at the outer edge.


title: Cosmological Constraints with Clustering-Based Redshifts

abstract: We demonstrate that observations lacking reliable redshift information, such as photometric and radio continuum surveys, can produce robust measurements of cosmological parameters when empowered by clustering-based redshift estimation. This method infers the redshift distribution based on the spatial clustering of sources, using cross-correlation with a reference dataset with known redshifts. Applying this method to the existing SDSS photometric galaxies, and projecting to future radio continuum surveys, we show that sources can be efficiently divided into several redshift bins, increasing their ability to constrain cosmological parameters. We forecast constraints on the dark-energy equation-of-state and on local non-gaussianity parameters. We explore several pertinent issues, including the tradeoff between including more sources versus minimizing the overlap between bins, the shot-noise limitations on binning, and the predicted performance of the method at high redshifts. Remarkably, we find that, once this technique is implemented, constraints on dynamical dark energy from the SDSS imaging catalog can be competitive with, or better than, those from the spectroscopic BOSS survey and even future planned experiments. Further, constraints on primordial non-Gaussianity from future large-sky radio-continuum surveys can outperform those from the Planck CMB experiment, and rival those from future spectroscopic galaxy surveys. The application of this method thus holds tremendous promise for cosmology.


title: The Lyman-alpha luminosity function at z=5.7-6.6 and the steep drop of the faint end: implications for reionization

abstract: We present new results from the widest narrow band survey search for Lyman-alpha (Lya) emitters at z=5.7, just after reionization. We survey a total of 7 deg$^2$ spread over the COSMOS, UDS and SA22 fields. We find over 11,000 line emitters, out of which 514 are robust Lya candidates at z=5.7 within a volume of 6.3x10$^6$ Mpc$^3$. Our Lya emitters span a wide range in Lya luminosities, from faint to bright (L$_{\rm Ly\alpha}\sim10^{42.5-44}$ erg s$^{-1}$) and rest-frame equivalent widths (EW$_0$~25-1000 \AA) in a single, homogeneous data-set. By combining all our fields we find that the faint end slope of the z=5.7 Lya luminosity function is very steep, with $\alpha=-2.3^{+0.4}_{-0.3}$. We also present an updated z=6.6 Lya luminosity function, based on comparable volumes and obtained with the same methods, which we directly compare with that at z=5.7. We find a significant decline of the number density of faint Lya emitters from z=5.7 to z=6.6 (by $0.5\pm0.1$ dex), but no evolution at the bright end/no evolution in L*. Faint Lya emitters at z=6.6 show much more extended haloes than those at z=5.7, suggesting that neutral Hydrogen plays an important role, increasing the scattering and leading to observations missing faint Lya emission within the epoch of reionization. All together, our results suggest that we are observing patchy reionization which happens first around the brightest Lya emitters, allowing the number densities of those sources to remain unaffected by the increase of neutral Hydrogen fraction from z~5 to z~7.


title: Equilibrium model prediction for the scatter in the star-forming main sequence

abstract: The analytic "equilibrium model" for galaxy evolution using a mass balance equation is able to reproduce mean observed galaxy scaling relations between stellar mass, halo mass, star formation rate (SFR) and metallicity across the majority of cosmic time with a small number of parameters related to feedback. Here we aim to test this data-constrained model to quantify deviations from the mean relation between stellar mass and SFR, i.e. the star-forming galaxy main sequence (MS). We implement fluctuation in halo accretion rates parameterised from merger-based simulations, and quantify the intrinsic scatter introduced into the MS under the assumption that fluctuations in star formation follow baryonic inflow fluctuations. We predict the 1-sigma MS scatter to be ~ 0.2 - 0.25 dex over the stellar mass range 10^8 Mo to 10^11 Mo and a redshift range 0.5 < z < 3 for SFRs averaged over 100 Myr. The scatter increases modestly at z > 3, as well as by averaging over shorter timescales. The contribution from merger-induced star formation is generally small, around 5% today and 10 - 15% during the peak epoch of cosmic star formation. These results are generally consistent with available observations, suggesting that deviations from the MS primarily reflect stochasticity in the inflow rate owing to halo mergers.


title: Orbital eccentricities in primordial black holes binaries

abstract: It was recently suggested that the merger of $\sim30\,M_\odot$ primordial black holes (PBHs) may provide a significant number of events in gravitational-wave observatories over the next decade, if they make up an appreciable fraction of the dark matter. Here we show that measurement of the eccentricities of the inspiralling binary black holes can be used to distinguish these binaries from those produced by more traditional astrophysical mechanisms. These PBH binaries are formed on highly eccentric orbits and can then merge on timescales that in some cases are years or less, retaining some eccentricity in the last seconds before the merger. This is to be contrasted with massive-stellar-binary, globular-cluster, or other astrophysical origins for binary black holes (BBHs) in which the orbits have very effectively circularized by the time the BBH enters the observable LIGO window. Here we discuss the features of the gravitational-wave signals that indicate this eccentricity and forecast the sensitivity of LIGO and the Einstein Telescope to such effects. We show that if PBHs make up the dark matter, then roughly one event should have a detectable eccentricity given LIGO's expected sensitivity and observing time of six years. The Einstein Telescope should see $O(10)$ such events after ten years.


title: Secular dynamics of multiplanet systems: implications for the formation of hot and warm Jupiters via high-eccentricity migration

abstract: Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.