title: Counting curves on Hirzebruch surfaces: tropical geometry and the Fock space

abstract: We study the stationary descendant Gromov-Witten theory of toric surfaces by combining and extending a range of techniques - tropical curves, floor diagrams, and Fock spaces. A correspondence theorem is established between tropical curves and descendant invariants on toric surfaces using maximal toric degenerations. An intermediate degeneration is then shown to give rise to floor diagrams, giving a geometric interpretation of this well-known bookkeeping tool in tropical geometry. In the process, we extend floor diagram techniques to include descendants in arbitrary genus. These floor diagrams are then used to connect tropical curve counting to the algebra of operators on the bosonic Fock space, and are shown to coincide with the Feynman diagrams of appropriate operators. This extends work of a number of researchers, including Block-G\"ottche, Cooper -Pandharipande, and Block-Gathmann-Markwig.


title: A Conceptual Model for Holistic Classification of Insider

abstract: The process through which an insider to an organization can be described or classified is lined within the orthodox paradigm of classification in which an organization considers only subject with requisite employee criterion as an insider to that organization. This is further clouded with the relative rigidity in operational security policies being implemented in organizations. Establishing investigation process in instances of misuse occurrence and or ascertaining the efficiency of staff member using such archaic paradigm is maligned with endless possibilities of uncertainties. This study, therefore, proposes a holistic model for which insider classification can be crystallized using the combination of qualitative research process and analysis of moment structure evaluation process. A full comprehension of this proposition could serve as a hinge through which insider misuse investigation can be thoroughly carried out. In addition, integrating this paradigm into existing operational security policies could serve as a metric upon which an organization can understand insider dynamics, in order to prevent misuses, and enhance staff management.


title: Optimizing Quantum Walk Search on a Reduced Uniform Complete Multi-Partite Graph

abstract: In a recent work by Novo et al. (Sci. Rep. 5, 13304, 2015), the invariant subspace method was applied to the study of continuous-time quantum walk (CTQW). The method helps to reduce a graph into a simpler version that allows more transparent analyses of the quantum walk model. In this work, we adopt the aforementioned method to investigate the optimality of a quantum walk search of a marked element on a complete multi-partite graph. We formulate the eigenbasis that would facilitate the transport between the two lowest energy eigenstates and demonstrate how to set the appropriate coupling factor to preserve the optimality.


title: Piecewise Constant Martingales and Lazy Clocks

abstract: This paper discusses the possibility to find and construct \textit{piecewise constant martingales}, that is, martingales with piecewise constant sample paths evolving in a connected subset of $\mathbb{R}$. After a brief review of standard possible techniques, we propose a construction based on the sampling of latent martingales $\tilde{Z}$ with \textit{lazy clocks} $\theta$. These $\theta$ are time-change processes staying in arrears of the true time but that can synchronize at random times to the real clock. This specific choice makes the resulting time-changed process $Z_t=\tilde{Z}_{\theta_t}$ a martingale (called a \textit{lazy martingale}) without any assumptions on $\tilde{Z}$, and in most cases, the lazy clock $\theta$ is adapted to the filtration of the lazy martingale $Z$. This would not be the case if the stochastic clock $\theta$ could be ahead of the real clock, as typically the case using standard time-change processes. The proposed approach yields an easy way to construct analytically tractable lazy martingales evolving on (intervals of) $\mathbb{R}$.


title: Recent advances in accelerated discovery through machine learning and statistical inference

abstract: Recent applications of machine learning and statistical inference provide case studies demonstrating how such approaches can accelerate the discovery process in physical chemistry and related fields. Examples discussed in this review include the introduction of automated approaches to systematically improve experimental design, increase the efficiency of computationally expensive molecular simulations, facilitate construction of predictive models for complex biological processes, and discover interparticle potentials that lead to materials which meet specified design goals. A common theme is the synergy between experiment and computation enabled by such approaches.


title: Mining Social Media to Inform Peatland Fire and Haze Disaster Management

abstract: Peatland fires and haze events are disasters with national, regional and international implications. The phenomena lead to direct damage to local assets, as well as broader economic and environmental losses. Satellite imagery is still the main and often the only available source of information for disaster management. In this article, we test the potential of social media to assist disaster management. To this end, we compare insights from two datasets: fire hotspots detected via NASA satellite imagery and almost all GPS-stamped tweets from Sumatra Island, Indonesia, posted during 2014. Sumatra Island is chosen as it regularly experiences a significant number of haze events, which affect citizens in Indonesia as well as in nearby countries including Malaysia and Singapore. We analyse temporal correlations between the datasets and their geo-spatial interdependence. Furthermore, we show how Twitter data reveals changes in users' behavior during severe haze events. Overall, we demonstrate that social media is a valuable source of complementary and supplementary information for haze disaster management. Based on our methodology and findings, an analytics tool to improve peatland fire and haze disaster management by the Indonesian authorities is under development.


title: On the absence of pulses from pulsars

abstract: The aim of this thesis is to investigate the absence of pulses in pulsar, a phenomena also know as pulsar nulling. The motivation for the work done in thesis is listed in Chapter 2. A detailed discussion on the previous work carried out to investigate pulsar nulling phenomena is also included in Chapter 2. Chapter 3 discusses in details all the observations carried out for the thesis. In Chapter 4, details regarding the nulling behaviour of various pulsars are listed. Chapter 4 highlights the fact that nulling is unpredictable. Chapter 5 presents a comparison between two pulsars with high NFs. It is shown that both the pulsars exhibit similar bursting behaviour, but strikingly different quasi-periodicities and emission properties. Chapter 6 discusses the details regarding the simultaneous observations of two nulling pulsars at four different frequencies. These observations suggest that nulling is a frequency independent phenomena simultaneous occuring over a decade of frequency. Chapter 7 summaries the results obtained in three different studies along with their implications. The future work to extend the work carried out in the thesis is also listed in Chapter 7. The main conclusions of the thesis are as follow, (1) NF is not an ideal parameter to quantify nulling behaviour, confirmed firmly by comparing low NF pulsars as well as high NF pulsars, (2) Nulling occurs randomly with unpredictable length durations. Quasi-periodicities seen in the high NF pulsars can also be explained by the Markov models with a forcing function, (3) Nulling is an extreme form of mode-changing phenomena which occurs on a global magnetospheric scale, (4) Geometric reasons are less favoured as a likely cause of nulling phenomena due to the randomness and broadband behaviour reported in this thesis.


title: Statistical Gravitational Waveform Models: What to Simulate Next?

abstract: Models of gravitational waveforms play a critical role in detecting and characterizing the gravitational waves (GWs) from compact binary coalescences. Waveforms from numerical relativity (NR), while highly accurate, are too computationally expensive to produce to be directly used with Bayesian parameter estimation tools like Markov-chain-Monte-Carlo and nested sampling. We propose a Gaussian process regression (GPR) method to generate accurate reduced-order-model waveforms based only on existing accurate (e.g. NR) simulations. Using a training set of simulated waveforms, our GPR approach produces interpolated waveforms along with uncertainties across the parameter space. As a proof of concept, we use a training set of IMRPhenomD waveforms to build a GPR model in the 2-d parameter space of mass ratio $q$ and equal-and-aligned spin $\chi_1=\chi_2$. Using a regular, equally-spaced grid of 120 IMRPhenomD training waveforms in $q\in[1,3]$ and $\chi_1 \in [-0.5,0.5]$, the GPR mean approximates IMRPhenomD in this space to mismatches below $4.3\times 10^{-5}$. Our approach can alternatively use training waveforms directly from numerical relativity. Beyond interpolation of waveforms, we also present a greedy algorithm that utilizes the errors provided by our GPR model to optimize the placement of future simulations. In a fiducial test case we find that using the greedy algorithm to iteratively add simulations achieves GPR errors that are $\sim 1$ order of magnitude lower than the errors from using Latin-hypercube or square training grids.


title: Electronic and magnetic properties of epitaxial SrRhO3 films

abstract: Strong interplay of fundamental order parameters in complex oxides are known to give rise to exotic physical phenomena. The 4d transition metal oxide SrRhO3 has generated much interest, but advances have been hindered by difficulties in preparing single crystalline phases. Here, we have epitaxially stabilized high quality single crystalline SrRhO3 films and investigated their structural, electronic, and magnetic properties. We determine that their properties significantly differ from the paramagnetic metallic ground state that governs bulk samples and are strongly related to rotations of the RhO6 octahedra.


title: On approximate Gauss-Lucas theorems

abstract: The Gauss--Lucas theorem states that any convex set $K\subset\mathbb{C}$ which contains all $n$ zeros of a degree $n$ polynomial $p\in\mathbb{C}[z]$ must also contain all $n-1$ critical points of $p$. In this paper we explore the following question: for which choices of positive integers $n$ and $k$, and positive real number $\epsilon$, will it follow that for every degree $n$ polynomial $p$ with at least $k$ zeros lying in $K$, $p$ will have at least $k-1$ critical points lying in the $\epsilon$-neighborhood of $K$. We supply an inequality relating $n$, $k$, and $\epsilon$ which, when satisfied, guarantees a positive answer to the above question.


title: A quantitative generalization of Prodanov-Stoyanov Theorem on minimal Abelian topological groups

abstract: A topological group $X$ is defined to have $compact$ $exponent$ if for some number $n\in\mathbb N$ the set $\{x^n:x\in X\}$ has compact closure in $X$. Any such number $n$ will be called a compact exponent of $X$. Our principal result states that a complete Abelian topological group $X$ has compact exponent (equal to $n\in\mathbb N$) if and only if for any injective continuous homomorphism $f:X\to Y$ to a topological group $Y$ and every $y\in \bar{f(X)}$ there exists a positive number $k$ (equal to $n$) such that $y^k\in f(X)$. This result has many interesting implications: (1) an Abelian topological group is compact if and only if it is complete in each weaker Hausdorff group topology; (2) each minimal Abelian topological group is precompact (this is the famous Prodanov-Stoyanov Theorem); (3) a topological group $X$ is complete and has compact exponent if and only if it is closed in each Hausdorff paratopological group containing $X$ as a topoloical subgroup (this confirms an old conjecture of Banakh and Ravsky).


title: A Practical Algorithm for Enumerating Collinear Points

abstract: This paper studies the problem of enumerating all maximal collinear subsets of size at least three in a given set of $n$ points. An algorithm for this problem, besides solving degeneracy testing and the exact fitting problem, can also help with other problems, such as point line cover and general position subset selection. The classic \emph{topological sweeping} algorithm of Edelsbrunner and Guibas can find these subsets in $O(n^2)$ time in the dual plane. We present an alternative algorithm that, although asymptotically slower than their algorithm in the worst case, is simpler to implement and more amenable to parallelization. If the input points are decomposed into $m$ convex polygons, our algorithm has time complexity $O(n^2 \log m)$ and space complexity $O(n)$. Our algorithm can be parallelized on the CREW PRAM with time complexity $O(n \log m)$ using $n$ processors.


abstract: The NSF Workshop in Quantum Information and Computation for Chemistry assembled experts from directly quantum-oriented fields such as algorithms, chemistry, machine learning, optics, simulation, and metrology, as well as experts in related fields such as condensed matter physics, biochemistry, physical chemistry, inorganic and organic chemistry, and spectroscopy. The goal of the workshop was to summarize recent progress in research at the interface of quantum information science and chemistry as well as to discuss the promising research challenges and opportunities in the field. Furthermore, the workshop hoped to identify target areas where cross fertilization among these fields would result in the largest payoff for developments in theory, algorithms, and experimental techniques. The ideas can be broadly categorized in two distinct areas of research that obviously have interactions and are not separated cleanly. The first area is quantum information for chemistry, or how quantum information tools, both experimental and theoretical can aid in our understanding of a wide range of problems pertaining to chemistry. The second area is chemistry for quantum information, which aims to discuss the several aspects where research in the chemical sciences can aid progress in quantum information science and technology. The results of the workshop are summarized in this report.


title: Continuous easy-plane deconfined phase transition on the kagome lattice

abstract: We use large scale quantum Monte-Carlo simulations to study an extended Hubbard model of hardcore bosons on the kagome lattice. In the limit of strong nearest-neighbor interactions at 1/3 filling, the interplay between frustration and quantum fluctuations leads to a valence bond solid ground state. The system undergoes a quantum phase transition to a superfluid phase as the interaction strength is decreased. It is still under debate whether the transition is weakly first order or represents an unconventional continuous phase transition. We present a theory in terms of an easy-plane NCCP$^1$ gauge theory describing the phase transition at 1/3 filling. Utilizing large scale quantum Monte-Carlo simulations with parallel tempering in the canonical ensemble up to 15552 spins, we provide evidence that the phase transition is continuous at exactly 1/3 filling. A careful finite size scaling analysis reveals an unconventional scaling behavior hinting at deconfined quantum criticality.


title: Block-Matching Optical Flow for Dynamic Vision Sensor- Algorithm and FPGA Implementation

abstract: Rapid and low power computation of optical flow (OF) is potentially useful in robotics. The dynamic vision sensor (DVS) event camera produces quick and sparse output, and has high dynamic range, but conventional OF algorithms are frame-based and cannot be directly used with event-based cameras. Previous DVS OF methods do not work well with dense textured input and are designed for implementation in logic circuits. This paper proposes a new block-matching based DVS OF algorithm which is inspired by motion estimation methods used for MPEG video compression. The algorithm was implemented both in software and on FPGA. For each event, it computes the motion direction as one of 9 directions. The speed of the motion is set by the sample interval. Results show that the Average Angular Error can be improved by 30\% compared with previous methods. The OF can be calculated on FPGA with 50\,MHz clock in 0.2\,us per event (11 clock cycles), 20 times faster than a Java software implementation running on a desktop PC. Sample data is shown that the method works on scenes dominated by edges, sparse features, and dense texture.


title: Epidemiology of Objectively Measured Bedtime and Chronotype in the US adolescents and adults: NHANES 2003-2006

abstract: Background: We propose a method for estimating the timing of in-bed intervals using objective data in a large representative U.S. sample, and quantify the association between these intervals and age, sex, and day of the week. Methods: The study included 11,951 participants six years and older from the National Health and Nutrition Examination Survey (NHANES) 2003-2006, who wore accelerometers to measure physical activity for seven consecutive days. Participants were instructed to remove the device just before the nighttime sleep period and put it back on immediately after. This nighttime period of non-wear was defined in this paper as the objective bedtime (OBT), an objectively estimated record of the in-bed-interval. For each night of the week, we estimated two measures: the duration of the OBT (OBT-D) and, as a measure of the chronotype, the midpoint of the OBT (OBT-M). We estimated day-of-the-week-specific OBT-D and OBT-M using gender-specific population percentile curves. Differences in OBT-M (chronotype) and OBT-D (the amount of time spent in bed) by age and sex were estimated using regression models. Results: The estimates of OBT-M and their differences among age groups were consistent with the estimates of chronotype obtained via self-report in European populations. The average OBT-M varied significantly by age, while OBT-D was less variable with age. The most pronounced differences were observed between OBT-M of weekday and weekend nights. Conclusions: The proposed measures, OBT-D and OBT-M, provide useful information of time in bed and chronotype in NHANES 2003-2006. They identify within-week patterns of bedtime and can be used to study associations between the bedtime and the large number of health outcomes collected in NHANES 2003-2006.


title: Random walks on uniform and non-uniform combs and brushes

abstract: We consider random walks on comb- and brush-like graphs consisting of a base (of fractal dimension $D$) decorated with attached side-groups. The graphs are also characterized by the fractal dimension $D_a$ of a set of anchor points where side-groups are attached to the base. Two types of graphs are considered. Graphs of the first type are uniform in the sense that anchor points are distributed periodically over the base, and thus form a subset of the base with dimension $D_a=D$. Graphs of the second type are decorated with side-groups in a regular yet non-uniform way: the set of anchor points has fractal dimension smaller than that of the base, $D_a<D$. For uniform graphs, a qualitative method for evaluating the sub-diffusion exponent suggested by Forte et al. for combs ($D=1$) is extended for brushes ($D>1$) and numerically tested for the Sierpinski brush (with the base and anchor set built on the same Sierpinski gasket). As an example of nonuniform graphs we consider the Cantor comb composed of a one-dimensional base and side-groups, the latter attached to the former at anchor points forming the Cantor set. A peculiar feature of this and other nonuniform systems is a long-lived regime of super-diffusive transport when side-groups are of a finite size.


title: Model-based Iterative Reconstruction for Flat-Panel Cone-Beam CT with Focal Spot Blur, Detector Blur, and Correlated Noise

abstract: While model-based reconstruction methods have been successfully applied to flat-panel cone-beam CT (FP-CBCT) systems, typical implementations ignore both spatial correlations in the projection data as well as system blurs due to the detector and focal spot in the x-ray source. In this work, we develop a forward model for flat-panel-based systems that includes blur and noise correlation associated with finite focal spot size and an indirect detector (e.g., scintillator). This forward model is used to develop a staged reconstruction framework where projection data are deconvolved and log-transformed, followed by a generalized least-squares reconstruction that utilizes a non-diagonal statistical weighting to account for the correlation that arises from the acquisition and data processing chain. We investigate the performance of this novel reconstruction approach in both simulated data and in CBCT test-bench data. In comparison to traditional filtered backprojection and model-based methods that ignore noise correlation, the proposed approach yields a superior noise-resolution tradeoff. For example, for a system with 0.34 mm FWHM scintillator blur and 0.70 FWHM focal spot blur, using the a correlated noise model instead of an uncorrelated noise model increased resolution by 42% (with variance matched at 6.9 x 10^-8 mm^-2 ). While this advantage holds across a wide range of systems with differing blur characteristics, the improvements are greatest for systems where source blur is larger than detector blur.


title: A Distributed Cooperative Control Framework for Synchronized Reconnection of a Multi-Bus Microgrid

abstract: One critical value microgrids bring to power systems is resilience, the capability of being able to island from the main grid under certain conditions and connect back when necessary. Once islanded, a microgrid must be synchronized to the main grid before reconnection to prevent severe consequences. In general, synchronization of a single machine with the grid can be easily achieved using a synchronizer. The problem becomes more challenging when it comes to a multi-bus microgrid with multiple distributed generators (DGs) and dispersed loads. All distributed generators need to be properly controlled in a coordinated way to achieve synchronization. This paper presents a novel bi-level distributed cooperative control framework for a multi-bus microgrid. In this framework, DGs work collaboratively in a distributed manner using the minimum and sparse communication. The topology of the communication network can be flexible which supports the plug-and-play feature of microgrids. Fast and deterministic synchronization can be achieved with tolerance to communication latency. Experimental results obtained from Hardware-in-the-Loop (HIL) simulation demonstrate the effectiveness of the proposed approach.


title: Thermometry of Silicon Nanoparticles

abstract: Current thermometry techniques lack the spatial resolution required to see the temperature gradients in typical, highly-scaled modern transistors. As a step toward addressing this problem, we have measured the temperature dependence of the volume plasmon energy in silicon nanoparticles from room temperature to 1250$^\circ$C, using a chip-style heating sample holder in a scanning transmission electron microscope (STEM) equipped with electron energy loss spectroscopy (EELS). The plasmon energy changes as expected for an electron gas subject to the thermal expansion of silicon. Reversing this reasoning, we find that measurements of the plasmon energy provide an independent measure of the nanoparticle temperature consistent with that of the heater chip's macroscopic heater/thermometer to within the 5\% accuracy of the chip thermometer's calibration. Thus silicon has the potential to provide its own, high-spatial-resolution thermometric readout signal via measurements of its volume plasmon energy. Furthermore, nanoparticles in general can serve as convenient nanothermometers for \emph{in situ} electron microscopy experiments.


title: The evolution of Giant Molecular Filaments

abstract: In recent years there has been a growing interest in studying giant molecular filaments (GMFs), which are extremely elongated (> 100pc in length) giant molecular clouds (GMCs). They are often seen as inter-arm features in external spiral galaxies, but have been tentatively associated with spiral arms when viewed in the Milky Way. In this paper, we study the time evolution of GMFs in a high-resolution section of a spiral galaxy simulation, and their link with spiral arm GMCs and star formation, over a period of 11Myrs. The GMFs generally survive the inter-arm passage, although they are subject to a number of processes (e.g. star formation, stellar feedback and differential rotation) which can break the giant filamentary structure into smaller sections. The GMFs are not gravitationally bound clouds as a whole, but are, to some extent, confined by external pressure. Once they reach the spiral arms, the GMFs tend to evolve into more substructured spiral arm GMCs, suggesting that GMFs may be precursors to arm GMCs. Here, they become incorporated into the more complex and almost continuum molecular medium that makes up the gaseous spiral arm. Instead of retaining a clear filamentary shape, their shapes are distorted both by their climb up the spiral potential and their interaction with the gas within the spiral arm. The GMFs do tend to become aligned with the spiral arms just before they enter them (when they reach the minimum of the spiral potential), which could account for the observations of GMFs in the Milky Way.


title: Angular power spectrum of galaxies in the 2MASS Redshift Survey

abstract: We present the measurement and interpretation of the angular power spectrum of nearby galaxies in the 2MASS Redshift Survey catalog with spectroscopic redshifts up to $z\approx 0.1$. We detect the angular power spectrum up to a multipole of $\ell\approx 1000$. We find that the measured power spectrum is dominated by galaxies living inside nearby galaxy clusters and groups. We use the halo occupation distribution (HOD) formalism to model the power spectrum, obtaining a fit with reasonable parameters. These HOD parameters are in agreement with the 2MASS galaxy distribution we measure toward the known nearby galaxy clusters, confirming validity of our analysis.


title: Weighted counting of solutions to sparse systems of equations

abstract: Given complex numbers $w_1, \ldots, w_n$, we define the weight $w(X)$ of a set $X$ of 0-1 vectors as the sum of $w_1^{x_1} \cdots w_n^{x_n}$ over all vectors $(x_1, \ldots, x_n)$ in $X$. We present an algorithm, which for a set $X$ defined by a system of homogeneous linear equations with at most $r$ variables per equation and at most $c$ equations per variable, computes $w(X)$ within relative error $\epsilon >0$ in $(rc)^{O(\ln n-\ln \epsilon)}$ time provided $|w_j| \leq \beta (r \sqrt{c})^{-1}$ for an absolute constant $\beta >0$ and all $j=1, \ldots, n$. A similar algorithm is constructed for computing the weight of a linear code over ${\Bbb F}_p$. Applications include counting weighted perfect matchings in hypergraphs, counting weighted graph homomorphisms, computing weight enumerators of linear codes with sparse code generating matrices, and computing the partition functions of the ferromagnetic Potts model at low temperatures and of the hard-core model at high fugacity on biregular bipartite graphs.


title: Proceedings of the Third International Conference on Cloud and Robotics (ICCR2016)

abstract: The 3rd International conference on Cloud and Robotics (ICCR2016) was held November 23-23, 2016 in Saint Quentin (France). ICCR is the premier gathering of practitioners and researchers interested in how to bring the power of Cloud computing and Robotics to a new cutting-edge domain: Cloud robotics. The objective of ICCR is a working conference, where together Cloud computing and robotics researchers/practitioners to exchange their ideas. Different with traditional academic conference, ICCR is also a forum for researchers and practitioners in the two disciplines, fostering the collaboration of Cloud computing with robotics, and for practitioners to show their working projects and to discuss and exchange their problems with researchers to find a solution.


title: Projective indecomposable modules and quivers for monoid algebras

abstract: We give a construction of the projective indecomposable modules and a description of the quiver for a large class of monoid algebras including the algebra of any finite monoid whose principal right ideals have at most one idempotent generator. Our results include essentially all families of finite monoids for which this has been done previously, for example, left regular bands, $\mathscr J$-trivial and $\mathscr R$-trivial monoids and left regular bands of groups.


title: Motion of a thin elliptic plate under symmetric and asymmetric orthotropic friction forces

abstract: Anisotropy of friction force is proved to be an important factor in various contact problems. We study dynamical behavior of thin plates with respect to symmetric and asymmetric orthotropic friction. Terminal motion of plates with circular and elliptic contact areas is mainly analyzed. Evaluation of friction forces for both symmetric and asymmetric orthotropic cases are shown. Regular pressure distribution is considered. Differential equations are formulated and solved numerically for a number of initial conditions. Examples show significant influence of friction force asymmetry on the motion.


title: Slow-Mode Oscillations of Hot Loops Excited at Flaring Footpoints

abstract: The analysis of a hot loop oscillation event using SOHO/SUMER, GOES/SXI, and RHESSI observations is presented. Damped Doppler shift oscillations were detected in the Fe XIX line by SUMER, and interpreted as a fundamental standing slow mode. The evolution of soft X-ray emission from GOES/SXI and hard X-ray sources from RHESSI suggests that the oscillations of a large loop are triggered by a small flare, which may be produced by interaction (local reconnection) of this large loop with a small loop at its footpoint. This study provides clear evidence supporting our early conjecture that the slow-mode standing waves in hot coronal loops are excited by impulsive heating (small or microflares) at the loop's footpoint.


title: Vortex Flows in the Solar Atmoshpere: Automated Identification and Statistical Analysis

abstract: Aims. Due to the fundamental importance of vortices on the photosphere, in this work we aim to fully automate the process of intensity vortex identification to facilitate a more robust statistical analysis of their properties. Methods. Using on-disk observational data of the Fe I continuum, the process of vortex identification is fully automated, for the first time in solar physics, with the help of an established method from hydrodynamics initially employed for the study of eddies in turbulent flows (Graftieaux et al. 2001). Results. We find that the expected lifetime of intensity vortices is much shorter (~ 17s) compared with previously observed magnetic bright point swirls. Our findings suggest that at any time there are 1.4e6 such small-scale intensity vortices covering about 2.8% of the total surface of the solar photosphere.


title: Modeling Biological Problems in Computer Science: A Case Study in Genome Assembly

abstract: As computer scientists working in bioinformatics/computational biology, we often face the challenge of coming up with an algorithm to answer a biological question. This occurs in many areas, such as variant calling, alignment, and assembly. In this tutorial, we use the example of the genome assembly problem to demonstrate how to go from a question in the biological realm to a solution in the computer science realm. We show the modeling process step-by-step, including all the intermediate failed attempts.


title: A new model of the chemistry of ionizing radiation in solids: CIRIS

abstract: The collisions between high-energy ions and solids can result in significant physical and chemical changes to the material. These effects are potentially important for better understanding the chemistry of interstellar and planetary bodies, which are exposed to cosmic radiation and the solar wind, respectively, however, modeling such collisions on a detailed microscopic basis has thus far been largely unsuccessful. To that end, a new model, entitled CIRIS: the Chemistry of Ionizing Radiation in Solids, was created to calculate the physical and chemical effects of the irradiation of solid materials. With the new code, we simulate O2 ice irradiated with 100 keV protons. Our models are able to reproduce the measured ozone abundances of a previous experimental study, as well as independently predict the approximate thickness of the ice used in that work.


title: Centralized Multi-Node Repair Regenerating Codes

abstract: In a distributed storage system, recovering from multiple failures is a critical and frequent task that is crucial for maintaining the system's reliability and fault-tolerance. In this work, we focus on the problem of repairing multiple failures in a centralized way, which can be desirable in many data storage configurations, and we show that a significant repair traffic reduction is possible. First, the fundamental tradeoff between the repair bandwidth and the storage size for functional repair is established. Using a graph-theoretic formulation, the optimal tradeoff is identified as the solution to an integer optimization problem, for which a closed-form expression is derived. Expressions of the extreme points, namely the minimum storage multi-node repair (MSMR) and minimum bandwidth multi-node repair (MBMR) points, are obtained. Second, we describe a general framework for converting single erasure minimum storage regenerating codes to MSMR codes. The repair strategy for $e$ failures is similar to that for single failure, however certain extra requirements need to be satisfied by the repairing functions for single failure. For illustration, the framework is applied to product-matrix codes and interference alignment codes. Furthermore, we prove that the functional MBMR point is not achievable for linear exact repair codes. We also show that exact-repair minimum bandwidth cooperative repair (MBCR) codes achieve an interior point, that lies near the MBMR point, when $k \equiv 1 \mod e$, $k$ being the minimum number of nodes needed to reconstruct the entire data. Finally, for $k> 2e, e\mid k$ and $e \mid d$, where $d$ is the number of helper nodes during repair, we show that the functional repair tradeoff is not achievable under exact repair, except for maybe a small portion near the MSMR point, which parallels the results for single erasure repair by Shah et al.


title: Obfuscation in Bitcoin: Techniques and Politics

abstract: In the cryptographic currency Bitcoin, all transactions are recorded in the blockchain - a public, global, and immutable ledger. Because transactions are public, Bitcoin and its users employ obfuscation to maintain a degree of financial privacy. Critically, and in contrast to typical uses of obfuscation, in Bitcoin obfuscation is not aimed against the system designer but is instead enabled by design. We map sixteen proposed privacy-preserving techniques for Bitcoin on an obfuscation-vs.-cryptography axis, and find that those that are used in practice tend toward obfuscation. We argue that this has led to a balance between privacy and regulatory acceptance.


title: Combining Stream Mining and Neural Networks for Short Term Delay Prediction

abstract: The systems monitoring the location of public transport vehicles rely on wireless transmission. The location readings from GPS-based devices are received with some latency caused by periodical data transmission and temporal problems preventing data transmission. This negatively affects identification of delayed vehicles. The primary objective of the work is to propose short term hybrid delay prediction method. The method relies on adaptive selection of Hoeffding trees, being stream classification technique and multilayer perceptrons. In this way, the hybrid method proposed in this study provides anytime predictions and eliminates the need to collect extensive training data before any predictions can be made. Moreover, the use of neural networks increases the accuracy of the predictions compared with the use of Hoeffding trees only.


title: The Pseudosynchronization of Binary Stars Undergoing Strong Tidal Interactions

abstract: Eccentric binaries known as heartbeat stars experience strong dynamical tides as the stars pass through periastron, providing a laboratory to study tidal interactions. We measure the rotation periods of 24 heartbeat systems, using the Kepler light curves to identify rotation peaks in the Fourier transform. Where possible, we compare the rotation period to the pseudosynchronization period derived by Hut 1981. Few of our heartbeat stars are pseudosynchronized with the orbital period. For four systems, we were able to identify two sets of rotation peaks, which we interpret as the rotation from both stars in the binary. The majority of the systems have a rotation period that is approximately 3/2 times the pseudosynchronization period predicted by Hut 1981, suggesting that other physical mechanisms influence the stars' rotation, or that stars typically reach tidal spin equilibrium at a rotation period slightly longer than predicted.


title: Solar neutrino physics on the beginning of 2017

abstract: This writeup is a review of current hot topics on solar neutrinos. It is based on a talk at the conference "Neutrinos: the quest for a new physics scale", held at the CERN on March 2017, where the Organizers entrusted me with a discussion of the provocative question "whether solar neutrino physics is over". Rather than providing a straight (negative) answer, in view of an audience consisting mostly of colleagues working in theoretical particle physics, I deemed it more useful providing a description of what is the current activity of the physicists working in solar neutrinos, leaving the listener free of forming his/her own opinion apropos.


title: Improving Distributed Gradient Descent Using Reed-Solomon Codes

abstract: Today's massively-sized datasets have made it necessary to often perform computations on them in a distributed manner. In principle, a computational task is divided into subtasks which are distributed over a cluster operated by a taskmaster. One issue faced in practice is the delay incurred due to the presence of slow machines, known as \emph{stragglers}. Several schemes, including those based on replication, have been proposed in the literature to mitigate the effects of stragglers and more recently, those inspired by coding theory have begun to gain traction. In this work, we consider a distributed gradient descent setting suitable for a wide class of machine learning problems. We adapt the framework of Tandon et al. (arXiv:1612.03301) and present a deterministic scheme that, for a prescribed per-machine computational effort, recovers the gradient from the least number of machines $f$ theoretically permissible, via an $O(f^2)$ decoding algorithm. We also provide a theoretical delay model which can be used to minimize the expected waiting time per computation by optimally choosing the parameters of the scheme. Finally, we supplement our theoretical findings with numerical results that demonstrate the efficacy of the method and its advantages over competing schemes.


title: All-optical logic gates in Stub photonic lattices

abstract: We experimentally study a Stub photonic lattice and excite their localized linear states originated from an isolated Flat Band at the center of the linear spectrum. By exciting these modes in different regions of the lattice, we observe that they do not diffract across the system and remain well trapped after propagating along the crystal. By using their wave nature, we are able to combine -- in phase and out of phase -- two neighbor states into a coherent superposition. These observations allow us to propose a novel setup for performing three different all-optical logical operations such as OR, AND, and XOR, positioning Flat Band systems as key setups to perform concrete applications at any level of power.


title: Numerical Simulation of the Hydrodynamical Combustion to Strange Quark Matter in the Trapped Neutrino Regime

abstract: We simulate and study the microphysics of combustion (flame burning) of two flavored quark matter (u, d) to three flavoured quark matter (u,d,s) in a trapped neutrino regime applicable to conditions prevailing in a hot proto-neutron star. The reaction-diffusion-advection equations for (u,d) to (u,d,s) combustion are coupled with neutrino transport, which is modelled through a flux-limited diffusion scheme. The flame speed is proportional to initial lepton fraction because of the release of electron chemical potential as heat, and reaches a steady-state burning speed of (0.001-0.008)c. We find that the burning speed is ultimately driven by the neutrino pressure gradient, given that the pressure gradient induced by quarks is opposed by the pressure gradients induced by electrons. This suggests, somewhat counter-intuitively, that the pressure gradients that drive the interface are controlled primarily by leptonic weak decays rather than by the quark Equation of State (EOS). In other words, the effects of the leptonic weak interaction, including the corresponding weak decay rates and the EOS of electrons and neutrinos, are at least as important as the uncertainties related to the EOS of high density matter. We find that for baryon number densities nB <= 0.35 fm-3, strong pressure gradients induced by leptonic weak decays drastically slow down the burning speed, which is thereafter controlled by the much slower burning process driven by backflowing downstream matter. We discuss the implications of our findings to proto-neutron stars.


title: Control Variates for Stochastic Gradient MCMC

abstract: It is well known that Markov chain Monte Carlo (MCMC) methods scale poorly with dataset size. A popular class of methods for solving this issue is stochastic gradient MCMC. These methods use a noisy estimate of the gradient of the log posterior, which reduces the per iteration computational cost of the algorithm. Despite this, there are a number of results suggesting that stochastic gradient Langevin dynamics (SGLD), probably the most popular of these methods, still has computational cost proportional to the dataset size. We suggest an alternative log posterior gradient estimate for stochastic gradient MCMC, which uses control variates to reduce the variance. We analyse SGLD using this gradient estimate, and show that, under log-concavity assumptions on the target distribution, the computational cost required for a given level of accuracy is independent of the dataset size. Next we show that a different control variate technique, known as zero variance control variates can be applied to SGMCMC algorithms for free. This post-processing step improves the inference of the algorithm by reducing the variance of the MCMC output. Zero variance control variates rely on the gradient of the log posterior; we explore how the variance reduction is affected by replacing this with the noisy gradient estimate calculated by SGMCMC.


title: Characterizing Lifting and Lowering Activities with Insole FSR sensors in Industrial Exoskeletons

abstract: This paper presents an insole FSR (Force Sensing Resistor) to dynamically detect weight variations in an exoskeleton system. The proposed methodology is intended for tasks of lifting and lowering heavy weights with an industrial exoskeleton to substantially reduce spinal loads during these manual handling activities. Instead of extensively placing high-dense force sensors by covering the whole plantar area, as most of commercial applications do, we integrate only a few force sensors in specific plantar area, so that the sensory system is not restricted to an individual foot size and shape, and on the other hand has relatively lower material cost. Industrial exoskeletons are intended to assist workers when handling heavy goods. With this in mind, wearers are not able to use their hands to control the exoskeleton since they use them to handle the goods. Therefore, the exoskeleton controller is required to indirectly infer how much and when the wearer requires assistance for lifting or lowering a heavy weight. Our approach of dynamically detect and characterize the increment/decrement of weight, as well as the rising/falling edge, enables the exoskeleton's controller to trigger the request of assistive force to the actuators.


title: Improved Convergence Rates for Distributed Resource Allocation

abstract: In this paper, we develop a class of decentralized algorithms for solving a convex resource allocation problem in a network of $n$ agents, where the agent objectives are decoupled while the resource constraints are coupled. The agents communicate over a connected undirected graph, and they want to collaboratively determine a solution to the overall network problem, while each agent only communicates with its neighbors. We first study the connection between the decentralized resource allocation problem and the decentralized consensus optimization problem. Then, using a class of algorithms for solving consensus optimization problems, we propose a novel class of decentralized schemes for solving resource allocation problems in a distributed manner. Specifically, we first propose an algorithm for solving the resource allocation problem with an $o(1/k)$ convergence rate guarantee when the agents' objective functions are generally convex (could be nondifferentiable) and per agent local convex constraints are allowed; We then propose a gradient-based algorithm for solving the resource allocation problem when per agent local constraints are absent and show that such scheme can achieve geometric rate when the objective functions are strongly convex and have Lipschitz continuous gradients. We have also provided scalability/network dependency analysis. Based on these two algorithms, we have further proposed a gradient projection-based algorithm which can handle smooth objective and simple constraints more efficiently. Numerical experiments demonstrates the viability and performance of all the proposed algorithms.


title: Probabilistic Jamming/Eavesdropping Attacks to Confuse a Buffer-Aided Transmitter-Receiver Pair

abstract: We assume that a buffer-aided transmitter communicates with a receiving node in the presence of an attacker. We investigate the impact of a radio-frequency energy-harvesting attacker that probabilistically operates as a jammer or an eavesdropper. We show that even without the need for an external energy source, the attacker can still degrade the security of the legitimate system. We show that the random data arrival behavior at the transmitter and the channel randomness of the legitimate link can improve the system's security. Moreover, we design a jamming scheme for the attacker and investigate its impact on the secure throughput of the legitimate system. The attacker designs his power splitting parameter and jamming/eavesdropping probability based on the energy state of the attacker's battery to minimize the secure throughput of the legitimate system.


title: Propagation of Solar Energetic Particles in Three-dimensional Interplanetary Magnetic Fields: Radial Dependence of Peak Intensities

abstract: A functional form I_{max}(R)=kR^{-\alpha}, where R is the radial distance of spacecraft, was usually used to model the radial dependence of peak intensities I_{max}(R) of solar energetic particles (SEPs). In this work, the five-dimensional Fokker-Planck transport equation incorporating perpendicular diffusion is numerically solved to investigate the radial dependence of SEP peak intensities. We consider two different scenarios for the distribution of spacecraft fleet: (1) along the radial direction line, (2) along the Parker magnetic field line. We find that the index \alpha in the above expression varies in a wide range, primarily depending on the properties (e.g., location, coverage) of SEP sources and on the longitudinal/latitudinal separations between the sources and the magnetic footpoints of the observers. Particularly, the situation that whether the magnetic footpoint of the observer is located inside or outside of the SEP source is a crucial factor determining the values of index \alpha. A two-phase phenomenon is found in the radial dependence of peak intensities. The "position" of the breakpoint (transition point/critical point) is determined by the magnetic connection status of the observers. This finding suggests that a very careful examination of magnetic connection between SEP source and each spacecraft should be taken in the observational studies. We obtain a lower limit of R^{-1.7\pm0.1} for empirically modelling the radial dependence of SEP peak intensities. Our findings in this work can be used to explain the majority of the previous multispacecraft survey results, and especially to reconcile the different/conflicting empirical values of index \alpha in the literature.

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


title: Character Values of Stanley Sequences

abstract: Stanley and Odlyzko proposed a method for greedily constructing sets with no 3-term arithmetic progressions. It is conjectured that there is a dichotomy between such sequences: those that have a periodic structure as the sequence satisfies certain recurrence relations while others appear to be chaotic. One large class of sequences that have these periodic behaviors are known as independent sequences that have two parameters, a character and a growth factor. It was conjectured by Rolnick that all but a finite set of integers can be achieved as characters of a independent sequences. Previously the only large class of integers known to be characters where those with base 3 representations consisting solely of the digits 0 and 2. This paper dramatically improves this result by demonstrating that all even integers not congruent to 244 mod 486 can be achieved as characters, therefore demonstrating that the set of all characters has a positive lower density.


title: A Stochastic Model for Short-Term Probabilistic Forecast of Solar Photo-Voltaic Power

abstract: In this paper, a stochastic model with regime switching is developed for solar photo-voltaic (PV) power in order to provide short-term probabilistic forecasts. The proposed model for solar PV power is physics inspired and explicitly incorporates the stochasticity due to clouds using different parameters addressing the attenuation in power.Based on the statistical behavior of parameters, a simple regime-switching process between the three classes of sunny, overcast and partly cloudy is proposed. Then, probabilistic forecasts of solar PV power are obtained by identifying the present regime using PV power measurements and assuming persistence in this regime. To illustrate the technique developed, a set of solar PV power data from a single rooftop installation in California is analyzed and the effectiveness of the model in fitting the data and in providing short-term point and probabilistic forecasts is verified. The proposed forecast method outperforms a variety of reference models that produce point and probabilistic forecasts and therefore portrays the merits of employing the proposed approach.


title: Adversarial Variational Bayes Methods for Tweedie Compound Poisson Mixed Models

abstract: The Tweedie Compound Poisson-Gamma model is routinely used for modeling non-negative continuous data with a discrete probability mass at zero. Mixed models with random effects account for the covariance structure related to the grouping hierarchy in the data. An important application of Tweedie mixed models is pricing the insurance policies, e.g. car insurance. However, the intractable likelihood function, the unknown variance function, and the hierarchical structure of mixed effects have presented considerable challenges for drawing inferences on Tweedie. In this study, we tackle the Bayesian Tweedie mixed-effects models via variational inference approaches. In particular, we empower the posterior approximation by implicit models trained in an adversarial setting. To reduce the variance of gradients, we reparameterize random effects, and integrate out one local latent variable of Tweedie. We also employ a flexible hyper prior to ensure the richness of the approximation. Our method is evaluated on both simulated and real-world data. Results show that the proposed method has smaller estimation bias on the random effects compared to traditional inference methods including MCMC; it also achieves a state-of-the-art predictive performance, meanwhile offering a richer estimation of the variance function.


title: High-resolution TNG spectra of T Tauri stars: Near-IR GIANO observations of the young variables XZ Tau and DR Tau

abstract: We used the TNG/GIANO instrument to obtain near-IR high-resolution spectra (R~50000) of XZ Tau and DR Tau, two actively accreting T Tauri stars classified as EXors. The analysis of the observed features provides information on the properties of the inner disk, the accretion columns, and the winds. Both sources display composite HI line profiles, with contributions from both accreting gas and high-velocity winds. These lines become more symmetric and narrower with increasing upper energy, which may be interpreted in terms of two components with different decrements or imputed to self-absorption effects. XZ Tau shows a relatively high state of activity with respect to literature observations. The variation of the HeI 1.08um line blue-shifted absorption suggests that the inner wind has undergone a dramatic change in its velocity structure, connected with a recent accretion event. DR Tau has a more stable wind as its HeI 1.08um absorption is stable in time in spite of strong variability of the emission component. The IR veiling can be interpreted as due to blackbody emission at temperatures of 1600K for XZ Tau and 2300K DR Tau, with emitting areas ~30 times larger than the star. For XZ Tau these conditions are consistent with emission from the inner rim of the dusty disk, while for DR Tau the fairly high temperature may suggest an origin from a thick gaseous disk within the dust sublimation radius. Strong and broad metallic lines, mainly from CI and FeI, are detected in XZ Tau, similar to those observed in other EXor sources during burst phases. DR Tau shows weaker and narrower metallic lines, despite its larger accretion luminosity. This suggests that accretion is not the only driver of metallic line excitation. Our observations show the potential of wide-band, high-resolution near-IR spectroscopy to simultaneously probe the phenomena that occur in the star-disk interaction region.


title: Advancements in Continuum Approximation Models for Logistics and Transportation Systems: 1996 - 2016

abstract: Continuum Approximation (CA) is an efficient and parsimonious technique for modeling complex logistics problems. In this paper,we review recent studies that develop CA models for transportation, distribution and logistics problems with the aim of synthesizing recent advancements and identifying current research gaps. This survey focuses on important principles and key results from CA models. In particular, we consider how these studies fill the gaps identified by the most recent literature reviews in this field. We observe that CA models are used in a wider range of applications, especially in the areas of facility location and integrated supply chain management. Most studies use CA as an alternative to exact solution approaches; however, CA can also be used in combination with exact approaches. We also conclude with promising areas of future work.


title: Moving mesh finite element simulation for phase-field modeling of brittle fracture and convergence of newton's iteration

abstract: A moving mesh finite element method is studied for the numerical solution of a phase-field model for brittle fracture. The moving mesh partial differential equation approach is employed to dynamically track crack propagation. Meanwhile, the decomposition of the strain tensor into tensile and compressive components is essential for the success of the phase-field modeling of brittle fracture but results in a non-smooth elastic energy and stronger nonlinearity in the governing equation. This makes the governing equation much more difficult to solve and, in particular, Newton's iteration often to fail to converge. Three regularization methods are proposed to smooth out the decomposition of the strain tensor. Numerical examples of fracture propagation under quasi-static load demonstrate that all of the methods can effectively improve the convergence of Newton's iteration for relatively small values of the regularization parameter but without comprising the accuracy of the numerical solution. They also show that the moving mesh finite element method is able to adaptively concentrate the mesh elements around propagating cracks and handle multiple and complex crack systems.


title: First Moment of Hecke $L$-functions with quartic characters at the central point

abstract: In this paper, we study the first moment of central values of Hecke $L$-functions associated with quartic characters.


title: Molecular Dynamics-Based Strength Estimates of Beta-Solenoid Proteins

abstract: The use of beta-solenoid proteins as functionalizable, nanoscale, self-assembling molecular building blocks may have many applications, including templating the growth of wires or higher-dimensional structures. By understanding their mechanical strengths, we can efficiently design the proteins for specific functions. We present a study of the mechanical properties of seven beta-solenoid proteins using GROMACS molecular dynamics software to produce force/torque-displacement data, implement umbrella sampling of bending/twisting trajectories, produce Potentials of Mean Force (PMFs), extract effective spring constants, and calculate rigidities for two bending and two twisting directions for each protein. We examine the differences between computing the strength values from force/torque-displacement data alone and PMF data, and show how higher precision estimates can be obtained from the former. In addition to the analysis of the methods, we report estimates for the bend/twist persistence lengths for each protein, which range from 0.5-3.4 $\mu$m. We note that beta-solenoid proteins with internal disulfide bridges do not enjoy enhanced bending or twisting strength, and that the strongest correlate with bend/twist rigidity is the number of hydrogen bonds per turn. In addition, we compute estimates of the Young's modulus ($Y$) for each protein, which range from $Y$ = 3.5 to 7.2 GPa.


title: A General Optimal Control Model of Human Movement Patterns I: Walking Gait

abstract: The biomechanics of the human body gives subjects a high degree of freedom in how they can execute movement. Nevertheless, subjects exhibit regularity in their movement patterns. One way to account for this regularity is to suppose that subjects select movement trajectories that are optimal in some sense. We adopt the principle that human movements are optimal and develop a general model for human movement patters that uses variational methods in the form of optimal control theory to calculate trajectories of movement trajectories of the body. We find that in this approach a constant of the motion that arises from the model and which plays a role in the optimal control model that is analogous to the role that the mechanical energy plays in classical physics. We illustrate how this approach works in practice by using it to develop a model of walking gait, making all the derivations and calculations in detail. We finally show that this optimal control model of walking gait recovers in an appropriate limit an existing model of walking gait which has been shown to provide good estimates of many observed characteristics of walking gait.


title: A note on two types of Lyapunov exponents and entropies for $\mathbb{Z}^k$-actions

abstract: In this paper, two types of Lyapunov exponents: random Lyapunov exponents and directional Lyapunov exponents, and the corresponding entropies: random entropy and directional entropy, are considered for smooth $\mathbb{Z}^k$-actions. The close relation among these quantities are investigated and the formulas of them are given via the Lyapunov exponents of the generators.


title: Constraints on the distance moduli, helium and metal abundances, and ages of globular clusters from their RR Lyrae and non-variable horizontal-branch stars. II. Multiple stellar populations in 47Tuc, M3, and M13

abstract: We present a new set of horizontal-branch (HB) models computed with the MESA stellar evolution code. The models adopt $\alpha$-enhanced \cite{ags09} metals mixtures and include the gravitational settling of He. They are used in our HB population synthesis tool to generate theoretical distributions of HB stars in order to describe the multiple stellar populations in the globular clusters 47Tuc, M3, and M13. The observed HB in 47Tuc is reproduced very well by our simulations for [Fe/H] $= -0.70$ and [$\alpha$/Fe] $= +0.4$ if the initial helium mass fraction varies by $\Delta Y_0 \sim 0.03$ and approximately 21%, 37%, and 42% of the stars have $Y_0 = 0.257$, 0.270, and 0.287, respectively. These simulations yield $(m-M)_V = 13.27$, implying an age near 13.0 Gyr. In the case of M3 and M13, our synthetic HBs for [Fe/H] $= -1.55$ and [$\alpha$/Fe] $= 0.4$ match the observed ones quite well if M3 has $\Delta Y_0 \sim 0.01$ and $(m-M)_V = 15.02$, resulting in an age of 12.6 Gyr, whereas M13 has $\Delta Y_0 \sim 0.08$ and $(m-M)_V = 14.42$, implying an age of 12.9 Gyr. Mass loss during giant-branch evolution and $\Delta Y_0$ appear to be the primary second parameters for M3 and M13. New observations for 7 of the 9 known RR Lyrae in M13 are also reported. Surprisingly, periods predicted for the $c$-type variables tend to be too high (by up to $\sim 0.1$~d).


title: Unbound motion on a Schwarzschild background: Practical approaches to frequency domain computations

abstract: Gravitational perturbations due to a point particle moving on a static black hole background are naturally described in Regge-Wheeler gauge. The first-order field equations reduce to a single master wave equation for each radiative mode. The master function satisfying this wave equation is a linear combination of the metric perturbation amplitudes with a source term arising from the stress-energy tensor of the point particle. The original master functions were found by Regge and Wheeler (odd parity) and Zerilli (even parity). Subsequent work by Moncrief and then Cunningham, Price and Moncrief introduced new master variables which allow time domain reconstruction of the metric perturbation amplitudes. Here I explore the relationship between these different functions and develop a general procedure for deriving new higher-order master functions from ones already known. The benefit of higher-order functions is that their source terms always converge faster at large distance than their lower-order counterparts. This makes for a dramatic improvement in both the speed and accuracy of frequency domain codes when analyzing unbound motion.


title: On the normality of the null-fiber of the moment map for $\theta$- and tori representations

abstract: Let (G, V) be a representation with either G a torus or (G, V) a locally free stable $\theta$-representation. We study the fiber at 0 of the associated moment map, which is a commuting variety in the latter case. We characterize the cases where this fiber is normal. The quotient (i.e. the symplectic reduction) turns out to be a specific orbifold when the representation is polar. In the torus case, this confirms a conjecture stated by C. Lehn, M. Lehn, R. Terpereau and the author in a former article. In the $\theta$-case, the conjecture was already known but the present approach yield another proof.


title: Eigenvalue Asymptotics of Narrow Domains

abstract: In this paper, we considered the spectrum of the Dirichlet Laplacian $\Delta_\epsilon$ on $\Omega_\epsilon=\{(x,y): -l_1<x<l_2, 0<y<\epsilon h(x)]\}$ where $ l_1,l_2>0$ and $h(x)$ is a positive analytic function having $0$ the only point where it achieves its global maximum $M$. In particular we studied in details about the full asymptotics of the eigenvalues.


title: An Investigation into Exoplanet Transits and Uncertainties

abstract: A simple transit model is described along with tests of this model against published results for 4 exoplanet systems (Kepler-1, 2, 8, and 77). Data from the Kepler mission are used. The Markov Chain Monte Carlo (MCMC) method is applied to obtain realistic error estimates. Optimization of limb darkening coefficients is subject to data quality. It is more likely for MCMC to derive an empirical limb darkening coefficient for light curves with S/N (signal to noise) above 15. Finally, the model is applied to Kepler data for 4 Kepler candidate systems (KOI 760.01, 767.01, 802.01, and 824.01) with previously unpublished results. Error estimates for these systems are obtained via the MCMC method.


title: Inferring probabilistic stellar rotation periods using Gaussian processes

abstract: Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasi-periodic --- spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalising over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine-fitting periodogram and an autocorrelation function method. We also demonstrate that it works well on real data, by inferring rotation periods for 275 Kepler stars with previously measured periods. We provide a table of rotation periods for these 1132 Kepler objects of interest and their posterior probability density function samples. Because this method delivers posterior probability density functions, it will enable hierarchical studies involving stellar rotation, particularly those involving population modelling, such as inferring stellar ages, obliquities in exoplanet systems, or characterising star-planet interactions. The code used to implement this method is available online.


title: Strongly regular Cayley graphs from partitions of subdifference sets of the Singer difference sets

abstract: In this paper, we give a new lifting construction of "hyperbolic" type of strongly regular Cayley graphs. Also we give new constructions of strongly regular Cayley graphs over the additive groups of finite fields based on partitions of subdifference sets of the Singer difference sets. Our results unify some recent constructions of strongly regular Cayley graphs related to $m$-ovoids and $i$-tight sets in finite geometry. Furthermore, some of the strongly regular Cayley graphs obtained in this paper are new or nonisomorphic to known strongly regular graphs with the same parameters.


title: Truly Multi-modal YouTube-8M Video Classification with Video, Audio, and Text

abstract: The YouTube-8M video classification challenge requires teams to classify 0.7 million videos into one or more of 4,716 classes. In this Kaggle competition, we placed in the top 3% out of 650 participants using released video and audio features. Beyond that, we extend the original competition by including text information in the classification, making this a truly multi-modal approach with vision, audio and text. The newly introduced text data is termed as YouTube-8M-Text. We present a classification framework for the joint use of text, visual and audio features, and conduct an extensive set of experiments to quantify the benefit that this additional mode brings. The inclusion of text yields state-of-the-art results, e.g. 86.7% GAP on the YouTube-8M-Text validation dataset.


title: State observation and sensor selection for nonlinear networks

abstract: A large variety of dynamical systems, such as chemical and biomolecular systems, can be seen as networks of nonlinear entities. Prediction, control, and identification of such nonlinear networks require knowledge of the state of the system. However, network states are usually unknown, and only a fraction of the state variables are directly measurable. The observability problem concerns reconstructing the network state from this limited information. Here, we propose a general optimization-based approach for observing the states of nonlinear networks and for optimally selecting the observed variables. Our results reveal several fundamental limitations in network observability, such as the trade-off between the fraction of observed variables and the observation length on one side, and the estimation error on the other side. We also show that owing to the crucial role played by the dynamics, purely graph- theoretic observability approaches cannot provide conclusions about one's practical ability to estimate the states. We demonstrate the effectiveness of our methods by finding the key components in biological and combustion reaction networks from which we determine the full system state. Our results can lead to the design of novel sensing principles that can greatly advance prediction and control of the dynamics of such networks.


title: Automatized photometric monitoring of active galactic nuclei with the 46 cm telescope of the Wise observatory

abstract: We present the first results of an ongoing variability monitoring program of active galactic nuclei (AGNs) using the 46\,cm telescope of the Wise observatory in Israel. The telescope has a field of view of $1.25^{\circ} \times 0.84^{\circ}$ and is specially equipped with five narrow band filters at 4300, 5200, 5700, 6200 and 7000\,\AA\, in order to perform photometric reverberation mapping studies of the central engine of AGNs. The program aims to observe a sample of 27 AGNs (V $<17$ mag) selected according to tentative continuum and line time delay measurements obtained in previous works. We describe the autonomous operation of the telescope together with the fully automatic pipeline used to achieve high-performance unassisted observations, data reduction, and light curves extraction using different photometric methods. The science verification data presented here, demonstrates the performance of the monitoring program in particular for efficiently photometric reverberation mapping of AGNs with additional capabilities to carry out complementary studies of other transient and variable phenomena such as variable stars studies.


title: A Thermodynamic Perspective of Negative-capacitance Field-effect-transistors

abstract: Physical phenomena underlying operation of ferroelectric field-effect transistors (FeFETs) is treated within a unified simulation framework. The framework incorporates the Landau mean-field treatment of free energy of a ferroelectric and the polarization dynamics according to Landau-Khalatnikov (LK) equation. These equations are self-consistently solved with the one-dimensional metal-oxide-semiconductor (MOS) structure electrostatics and the drift-diffusion solution for the current in the semiconductor channel. Numerical simulations demonstrate, depending on the ferroelectric (FE) thickness, both regimes of hysteresis switching (relevant for a non-volatile memory) and of higher on-currents and steeper subthreshold slope (SS) with a negligible hysteresis (relevant for logic) via the negative capacitance effect.


title: Plasmonic modes of polygonal particles calculated using a quantum hydrodynamics method

abstract: Plasmonic resonances of nanoparticles have drawn lots of attentions due to their interesting and useful properties such as strong field enhancements. These systems are typically studied using either classical electrodynamics or fully quantum theory. Each approach can handle some aspects of plasmonic systems accurately and efficiently, while having its own limitation. The self-consistent hydrodynamics model has the advantage that it can incorporate the quantum effect of the electron gas into classical electrodynamics in a consistent way. We use the method to study the plasmonic response of polygonal particles under the influence of an external electromagnetic wave, and we pay particular attention to the size and shape of the particle and the effect of charging. We find that the particles support edge modes, face modes and hybrid modes. The charges induced by the external field in the edge (face) modes mainly localize at the edges (faces), while the induced charges in the hybrid modes are distributed nearly evenly in both the edges and faces. The edge modes are less sensitive to particle size than the face modes, but are sensitive to the corner angles of the edges. When the number of sides of regular polygons increases, the edge and face modes gradually change into the classical dipole plasmonic mode of a cylinder. The hybrid modes are found to be the precursor of the Bennett mode, which cannot be found in classical electrodynamics.


title: Geodesic Motion of Neutral Particles around a Kerr-Newman Black Hole

abstract: We examine the dynamics of a neutral particle around a Kerr-Newman black hole, and in particular focus on the effects of the charge of the spinning black hole on the motion of the particle. We first consider the innermost stable circular orbits (ISCO) on the equatorial plane. It is found that the presence of the charge of the black hole leads to the effective potential of the particle with stronger repulsive effects as compared with the Kerr black hole. As a result, the radius of ISCO decreases as charge $Q$ of the black hole increases for a fixed value of black hole's angular momentum $a$. We then consider a kick on the particle from its initial orbit out of the equatorial motion. The perturbed motion of the particle will eventually be bounded, or unbounded so that it escapes to spatial infinity. Even more, the particle will likely be captured by the black hole. Thus we analytically and numerically determine the parameter regions of the corresponding motions, in terms of the initial radius of the orbital motion and the strength of the kick. The comparison will be made with the motion of a neutral particle in the Kerr black hole.


title: B-ducted Heating of Black Widow Companions

abstract: The companions of evaporating binary pulsars (black widows and related systems) show optical emission suggesting strong heating. In a number of cases large observed temperatures and asymmetries are inconsistent with direct radiative heating for the observed pulsar spindown power and expected distance. Here we describe a heating model in which the pulsar wind sets up an intrabinary shock (IBS) against the companion wind and magnetic field, and a portion of the shock particles duct along this field to the companion magnetic poles. We show that a variety of heating patterns, and improved fits to the observed light curves, can be obtained at expected pulsar distances and luminosities, at the expense of a handful of model parameters. We test this `IBS-B' model against three well observed binaries and comment on the implications for system masses.


title: Adversarial Network Coding

abstract: A combinatorial framework for adversarial network coding is presented. Channels are described by specifying the possible actions that one or more (possibly coordinated) adversaries may take. Upper bounds on three notions of capacity (the one-shot capacity, the zero-error capacity, and the compound zero-error capacity) are obtained for point-to-point channels, and generalized to corresponding capacity regions appropriate for multi-source networks. A key result of this paper is a general method by which bounds on these capacities in point-to-point channels may be ported to networks. This technique is illustrated in detail for Hamming-type channels with multiple adversaries operating on specific coordinates, which correspond, in the context of networks, to multiple adversaries acting on specific network edges. Capacity-achieving coding schemes are described for some of the considered adversarial models.


title: Majorana-Hubbard model on the square lattice

abstract: We study a tight-binding model of interacting Majorana (Hermitian) modes on a square lattice. The model may have an experimental realization in a superconducting-film--topological-insulator heterostructure in a magnetic field. We find a rich phase diagram, as a function of interaction strength, including a superfluid phase with spontaneous breaking of an emergent $U(1)$ symmetry, separated by a supersymmetric transition from a gapless normal phase.


title: Cosmology in the laboratory: an analogy between hyperbolic metamaterials and the Milne universe

abstract: This article shows that the compactified Milne universe geometry, a toy model for the big crunch/big bang transition, can be realized in hyperbolic metamaterials, a new class of nanoengineered systems which have recently found its way as an experimental playground for cosmological ideas. On one side, Klein-Gordon particles, as well as tachyons, are used as probes of the Milne geometry. On the other side, the propagation of light in two versions of a liquid crystal-based metamaterial provides the analogy. It is shown that ray and wave optics in the metamaterial mimic, respectively, the classical trajectories and wave function propagation, of the Milne probes, leading to the exciting perspective of realizing experimental tests of particle tunneling through the cosmic singularity, for instance.


title: Strong ordered Abelian groups and dp-rank

abstract: We provide an algebraic characterization of strong ordered Abelian groups: An ordered Abelian group is strong iff it has bounded regular rank and almost finite dimension. Moreover, we show that any strong ordered Abelian group has finite Dp-rank. We also provide a formula that computes the exact valued of the Dp-rank of any ordered Abelian group. In particular characterizing those ordered Abelian groups with Dp-rank equal to $n$. We also show the Dp-rank coincides with the Vapnik-Chervonenkis density.


title: Dynamic scaling in the 2D Ising spin glass with Gaussian couplings

abstract: We carry out simulated annealing and employ a generalized Kibble-Zurek scaling hypothesis to study the 2D Ising spin glass with normal-distributed couplings. The system has an equilibrium glass transition at temperature $T=0$. From a scaling analysis when $T\rightarrow 0$ at different annealing velocities, we extract the dynamic critical exponent $z$, i.e., the exponent relating the relaxation time $\tau$ to the system length $L$; $\tau\sim L^z$. We find $z=13.6 \pm 0.4$ for both the Edwards-Anderson spin-glass order parameter and the excess energy. This is different from a previous study of the system with bimodal couplings [S. J. Rubin, N. Xu, and A. W. Sandvik, Phys. Rev. E {\bf 95}, 052133 (2017)] where the dynamics is faster and the above two quantities relax with different exponents (and that of the energy is larger). We here argue that the different behaviors arise as a consequence of the different low-energy landscapes---for normal-distributed couplings the ground state is unique (up to a spin reflection) while the system with bimodal couplings is massively degenerate. Our results reinforce the conclusion of anomalous entropy-driven relaxation behavior in the bimodal Ising glass. In the case of a continuous coupling distribution, our results presented here indicate that, although Kibble-Zurek scaling holds, the perturbative behavior normally applying in the slow limit breaks down, likely due to quasi-degenerate states, and the scaling function takes a different form.


title: Large-type Artin groups are systolic

abstract: We prove that Artin groups from a class containing all large-type Artin groups are systolic. This provides a concise yet precise description of their geometry. Immediate consequences are new results concerning large-type Artin groups: biautomaticity; existence of $EZ$-boundaries; the Novikov conjecture; descriptions of finitely presented subgroups, of virtually solvable subgroups, and of centralizers for infinite order elements; the Burghelea conjecture; existence of low-dimensional models for classifying spaces for some families of subgroups.


title: Invariant vacuum

abstract: We apply the Lewis-Riesenfeld invariant method for the harmonic oscillator with time dependent mass and frequency to the modes of a charged scalar field that propagates in a curved, homogeneous and isotropic spacetime. We recover the Bunch-Davies vacuum in the case of a flat DeSitter spacetime, the equivalent one in the case of a closed DeSitter spacetime and the invariant vacuum in a curved spacetime that evolves adiabatically. In the three cases, it is computed the thermodynamical magnitudes of entanglement between the modes of the particles and antiparticles of the invariant vacuum, and the modification of the Friedmann equation caused by the existence of the energy density of entanglement. The amplitude of the vacuum fluctuations are also computed.


title: Fractional Bhatnagar-Gross-Krook kinetic equation

abstract: The linear Boltzmann equation approach is generalized to describe fractional superdiffusive transport of the Levy walk type in external force fields. The time distribution between scattering events is assumed to have a finite mean value and infinite variance. It is completely characterized by the two scattering rates, one fractional and a normal one, which defines also the mean scattering rate. We formulate a general fractional linear Boltzmann equation approach and exemplify it with a particularly simple case of the Bohm and Gross scattering integral leading to a fractional generalization of the Bhatnagar, Gross and Krook kinetic equation. Here, at each scattering event the particle velocity is completely randomized and takes a value from equilibrium Maxwell distribution at a given fixed temperature. We show that the retardation effects are indispensable even in the limit of infinite mean scattering rate and argue that this novel fractional kinetic equation provides a viable alternative to the fractional Kramers-Fokker-Planck (KFP) equation by Barkai and Silbey and its generalization by Friedrich et al. based on the picture of divergent mean time between scattering events. The case of divergent mean time is also discussed at length and compared with the earlier results obtained within the fractional KFP.


title: An Efficient Probabilistic Approach for Graph Similarity Search

abstract: Graph similarity search is a common and fundamental operation in graph databases. One of the most popular graph similarity measures is the Graph Edit Distance (GED) mainly because of its broad applicability and high interpretability. Despite its prevalence, exact GED computation is proved to be NP-hard, which could result in unsatisfactory computational efficiency on large graphs. However, exactly accurate search results are usually unnecessary for real-world applications especially when the responsiveness is far more important than the accuracy. Thus, in this paper, we propose a novel probabilistic approach to efficiently estimate GED, which is further leveraged for the graph similarity search. Specifically, we first take branches as elementary structures in graphs, and introduce a novel graph similarity measure by comparing branches between graphs, i.e., Graph Branch Distance (GBD), which can be efficiently calculated in polynomial time. Then, we formulate the relationship between GED and GBD by considering branch variations as the result ascribed to graph edit operations, and model this process by probabilistic approaches. By applying our model, the GED between any two graphs can be efficiently estimated by their GBD, and these estimations are finally utilized in the graph similarity search. Extensive experiments show that our approach has better accuracy, efficiency and scalability than other comparable methods in the graph similarity search over real and synthetic data sets.


title: Bayesian Conditional Generative Adverserial Networks

abstract: Traditional GANs use a deterministic generator function (typically a neural network) to transform a random noise input $z$ to a sample $\mathbf{x}$ that the discriminator seeks to distinguish. We propose a new GAN called Bayesian Conditional Generative Adversarial Networks (BC-GANs) that use a random generator function to transform a deterministic input $y'$ to a sample $\mathbf{x}$. Our BC-GANs extend traditional GANs to a Bayesian framework, and naturally handle unsupervised learning, supervised learning, and semi-supervised learning problems. Experiments show that the proposed BC-GANs outperforms the state-of-the-arts.


title: Magnon-phonon coupling effects on the indirect K-edge resonant inelastic X-ray scattering spectrum of a 2D Heisenberg antiferromagnet

abstract: We compute the effects of magnon-phonon coupling on the indirect K-edge bimagnon resonant inelastic x-ray scattering (RIXS) intensity spectrum of a square lattice Heisenberg antiferromagnet. We analyze the effects of competing nearest and next--nearest magnetic and magnon-phonon coupling interaction in the RIXS spectrum, for both the antiferromagnetic (AF) and the collinear antiferromagnetic (CAF) phases of the model. Utilizing the Dyson-Maleev representation of spin operators, the Bethe-Salpeter ladder approximation scheme for the bimagnon interacting channel, and considering the lowest order magnon-phonon-magnon scattering interaction we highlight distinct features in the X-ray spectrum. Considering damping effects, arising due to the presence of phonons, we find that in the AF phase the RIXS intensity spectrum attains a maximum value primarily localized around the K $\left(\pm\frac{\pi}{2}, \pm \frac{\pi}{2}\right)$ - point. For the CAF phase the intensity is broadly distributed with a significant scattering intensity located around the Y $\left(\pm\frac{\pi}{2}, 0\right)$ - point. Furthermore, in the CAF phase for suitable anisotropy, nearest-, and next-nearest neighbor interaction parameters the phonon effects can manifest itself as a distinct peak both below and above the bimagnon peak. Such a feature is in contrast to the antiferromagnetic spectrum where the effect due to the phonon peak was located consistently beyond the bimagnon peak in the high energy end of the spectrum. Additionally, in the CAF phase we find the RIXS bimagnon-phonon spectrum to be more sensitive to anisotropy compared to its antiferromagnetic counterpart. We conclude that the ultimate effect of magnon-phonon effects in the indirect K-edge RIXS spectrum, in both the antiferromagnetic and the collinear antiferromagnetic phase, is an observable effect.


title: A Data Envelopment Analysis (DEA)-Based Model for Power Interruption Cost Estimation for Industrial Companies

abstract: In this paper, a new model based on Data Envelopment Analysis (DEA) and Inverse Data Envelopment Analysis (IDEA) is presented for estimating the effect of electricity on the output of industrial companies. To this end, the effect of electricity deviation, which serves as one input that can influence a manufacturing company's final product, is evaluated. Intuitively, it is known that a direct relationship exists between electricity consumption and the output of the manufacturing company. However, finding a function that accurately represents this relationship is not easy. To check the applicability of the proposed method, it is tested on data from eight major vehicle manufacturing companies. In this model, labor hours, electrical energy consumption, and the value of raw materials are used as inputs, and the sales value is used as the proxy for the output of the company. These input and output data are used to find the efficiency of each company. Then, by changing the electricity consumption level, the output changes are derived. To calculate the outage cost, the deviation of the output is divided by the deviation of the electricity consumption, and the outage cost is estimated.


title: On the conjecture of Je\'smanowicz

abstract: We give a survey on some results covering the last 60 years concerning Je\'smanowicz' conjecture. Moreover, we conclude the survey with a new result by showing that the special Diophantine equation $$(20k)^x+(99k)^y=(101k)^z$$ has no solution other than $(x,y,z)=(2,2,2)$.


title: Foreground effect on the $J$-factor estimation of ultra-faint dwarf spheroidal galaxies

abstract: Dwarf spheroidal galaxies (dSphs) are promising targets for the gamma-ray dark matter (DM) search. In particular, DM annihilation signal is expected to be strong in some of the recently discovered nearby ultra-faint dSphs, which potentially give stringent constraints on the $\mathcal{O}(1)$ TeV WIMP DM. However, various non-negligible systematic uncertainties complicate the estimation of the astrophysical factors relevant for the DM search in these objects. Among them, the effects of foreground stars particularly attract attention because the contamination is unavoidable even for the future kinematical survey. In this article, we assess the effects of the foreground contamination on the astrophysical $J$-factor estimation by generating mock samples of stars in the four ultra-faint dSphs and using a model of future spectrographs. We investigate various data cuts to optimize the quality of the data and apply a likelihood analysis which takes member and foreground stellar distributions into account. We show that the foreground star contaminations in the signal region (the region of interest) and their statistical uncertainty can be estimated by interpolating the foreground star distribution in the control region where the foreground stars dominate the member stars. Such regions can be secured at future spectroscopic observations utilizing a multiple object spectrograph with a large field of view; e.g. the Prime Focus Spectrograph mounted on Subaru Telescope. The above estimation has several advantages: The data-driven estimation of the contamination makes the analysis of the astrophysical factor stable against the complicated foreground distribution. Besides, foreground contamination effect is considered in the likelihood analysis.


title: Nitrate removal from water using electrostatic regeneration of functionalized adsorbent

abstract: Nitrate is an important pollutant in drinking water worldwide, and a number of methods exist for the removal of nitrate from water including ion exchange and reverse osmosis. However, these approaches suffer from a variety of disadvantages including requirements for supply and disposal of brine used for regeneration in ion exchange and low water recovery ratio for reverse osmosis. Here, we demonstrate the use of high surface area activated carbon electrodes functionalized with moieties having high affinity for adsorption of nitrate from aqueous solution, such as those used in ion exchange. Adsorption of surfactant molecules having a quaternary amine ionic group to the activated carbon surfaces provides functionalization of the surfaces without complex chemistries. The functionalized electrodes have adsorption capacities of about 80 mg NaNO$_3$ per gram of activated carbon material. Unlike a traditional ion exchanger, the functionalized surfaces can be repeatedly regenerated by the application of an electrostatic potential which displaces the bound NO$_3^-$ while leaving an excess of electronic charge on the electrode. The cell is completed by an inert counter electrode where Faradaic reactions occur during regeneration. Up to approximately 40% of the initial capacity of the electrode is accessible following electrical regeneration.


title: A central limit theorem for the gossip process

abstract: The Aldous gossip process represents the dissemination of information in geographical space as a process of locally deterministic spread, augmented by random long range transmissions. Starting from a single initially informed individual, the proportion of individuals informed follows an almost deterministic path, but for a random time shift, caused by the stochastic behaviour in the very early stages of development. In this paper, it is shown that, even with the extra information available after a substantial development time, this broad description remains accurate to first order. However, the precision of the prediction is now much greater, and the random time shift is shown to have an approximately normal distribution, with mean and variance that can be computed from the current state of the process.


title: Electron-phonon coupling in the spin-split valence band of single layer WS$_2$

abstract: The absence of inversion symmetry leads to a strong spin-orbit splitting of the upper valence band of semiconducting single layer transition metal dichalchogenides such as MoS$_2$ or WS$_2$. This permits a direct comparison of the electron-phonon coupling strength in states that only differ by their spin. Here, the electron-phonon coupling in the valence band maximum of single-layer WS$_2$ is studied by first principles calculations and angle-resolved photoemission. The coupling strength is found to be drastically different for the two spin-split branches, with calculated values of $\lambda_K=$0.0021 and 0.40 for the upper and lower spin-split valence band of the free-standing layer, respectively. This difference is somewhat reduced when including scattering processes involving the Au(111) substrate present in the experiment and the experimental results confirm the strongly branch-dependent coupling strength.


title: Competing interactions in artificial spin chains

abstract: The low-energy magnetic configurations of artificial frustrated spin chains are investigated using magnetic force microscopy and micromagnetic simulations. Contrary to most studies on two-dimensional artificial spin systems where frustration arises from the lattice geometry, here magnetic frustration originates from competing interactions between neighboring spins. By tuning continuously the strength and sign of these interactions, we show that different magnetic phases can be stabilized. Comparison between our experimental findings and predictions from the one-dimensional Anisotropic Next-Nearest-Neighbor Ising (ANNNI) model reveals that artificial frustrated spin chains have a richer phase diagram than initially expected. Besides the observation of several magnetic orders and the potential extension of this work to highly-degenerated artificial spin chains, our results suggest that the micromagnetic nature of the individual magnetic elements allows observation of metastable spin configurations.


title: Parametric Inference for Discretely Observed Subordinate Diffusions

abstract: Subordinate diffusions are constructed by time changing diffusion processes with an independent L\'{e}vy subordinator. This is a rich family of Markovian jump processes which exhibit a variety of jump behavior and have found many applications. This paper studies parametric inference of discretely observed ergodic subordinate diffusions. We solve the identifiability problem for these processes using spectral theory and propose a two-step estimation procedure based on estimating functions. In the first step, we use an estimating function that only involves diffusion parameters. In the second step, a martingale estimating function based on eigenvalues and eigenfunctions of the subordinate diffusion is used to estimate the parameters of the L\'{e}vy subordinator and the problem of how to choose the weighting matrix is solved. When the eigenpairs do not have analytical expressions, we apply the constant perturbation method with high order corrections to calculate them numerically and the martingale estimating function can be computed efficiently. Consistency and asymptotic normality of our estimator are established considering the effect of numerical approximation. Through numerical examples, we show that our method is both computationally and statistically efficient. A subordinate diffusion model for VIX (CBOE volatility index) is developed which provides good fit to the data.


title: Random recursive trees and preferential attachment trees are random split trees

abstract: We consider linear preferential attachment trees, and show that they can be regarded as random split trees in the sense of Devroye (1999), although with infinite potential branching. In particular, this applies to the random recursive tree and the standard preferential attachment tree. An application is given to the sum over all pairs of nodes of the common number of ancestors.


title: $d_{xz/yz}$ Orbital Subband Structures and Chiral Orbital Angular Momentum in the (001) Surface States of SrTiO$_3$

abstract: We have performed angle resolved photoemission spectroscopy (ARPES) experiments on the surface states of SrTiO$_3$(001) using linearly and circularly polarized light to investigate the subband structures of out-of-plane $d_{xz/yz}$ orbitals and chiral orbital angular momentum (OAM). The data taken in the first Brillouin zone reveal new subbands for $d_{xz/yz}$ orbitals with Fermi wave vectors of 0.25 and 0.45 $\mathrm{\AA}^{-1}$ in addition to the previously reported ones. As a result, there are at least two subbands for all the Ti 3d t$_{2g}$ orbitals. Our circular dichroism ARPES data is suggestive of a chiral OAM structure in the surface states and may provide clues to the origin of the linear Rashba-like surface band splitting.