title: Ntuples for NLO Events at Hadron Colliders

abstract: We present an event-file format for the dissemination of next-to-leading-order (NLO) predictions for QCD processes at hadron colliders. The files contain all information required to compute generic jet-based infrared-safe observables at fixed order (without showering or hadronization), and to recompute observables with different factorization and renormalization scales. The files also make it possible to evaluate cross sections and distributions with different parton distribution functions. This in turn makes it possible to estimate uncertainties in NLO predictions of a wide variety of observables without recomputing the short-distance matrix elements. The event files allow a user to choose among a wide range of commonly-used jet algorithms and jet-size parameters. We provide event files for a $W$ or $Z$ boson accompanied by up to four jets, and for pure-jet events with up to four jets. The files are for the Large Hadron Collider with a center of mass energy of 7 or 8 TeV. A C++ library along with a Python interface for handling these files are also provided and described in this article. The library allows a user to read the event files and recompute observables transparently for different pdf sets and factorization and renormalization scales.


title: Neural perceptual model to global-local vision for recognition of the logical structure of administrative documents

abstract: This paper gives the definition of Transparent Neural Network "TNN" for the simulation of the globallocal vision and its application to the segmentation of administrative document image. We have developed and have adapted a recognition method which models the contextual effects reported from studies in experimental psychology. Then, we evaluated and tested the TNN and the multi-layer perceptron "MLP", which showed its effectiveness in the field of the recognition, in order to show that the TNN is clearer for the user and more powerful on the level of the recognition. Indeed, the TNN is the only system which makes it possible to recognize the document and its structure.


title: Hierarchical Clustering of Hyperspectral Images using Rank-Two Nonnegative Matrix Factorization

abstract: In this paper, we design a hierarchical clustering algorithm for high-resolution hyperspectral images. At the core of the algorithm, a new rank-two nonnegative matrix factorizations (NMF) algorithm is used to split the clusters, which is motivated by convex geometry concepts. The method starts with a single cluster containing all pixels, and, at each step, (i) selects a cluster in such a way that the error at the next step is minimized, and (ii) splits the selected cluster into two disjoint clusters using rank-two NMF in such a way that the clusters are well balanced and stable. The proposed method can also be used as an endmember extraction algorithm in the presence of pure pixels. The effectiveness of this approach is illustrated on several synthetic and real-world hyperspectral images, and shown to outperform standard clustering techniques such as k-means, spherical k-means and standard NMF.


title: Ranking basic belief assignments in decision making under uncertain environment

abstract: Dempster-Shafer theory (D-S theory) is widely used in decision making under the uncertain environment. Ranking basic belief assignments (BBAs) now is an open issue. Existing evidence distance measures cannot rank the BBAs in the situations when the propositions have their own ranking order or their inherent measure of closeness. To address this issue, a new ranking evidence distance (RED) measure is proposed. Compared with the existing evidence distance measures including the Jousselme's distance and the distance between betting commitments, the proposed RED measure is much more general due to the fact that the order of the propositions in the systems is taken into consideration. If there is no order or no inherent measure of closeness in the propositions, our proposed RED measure is reduced to the existing evidence distance. Numerical examples show that the proposed RED measure is an efficient alternative to rank BBAs in decision making under uncertain environment.


title: On Convergent Finite Difference Schemes for Variational - PDE Based Image Processing

abstract: We study an adaptive anisotropic Huber functional based image restoration scheme. By using a combination of L2-L1 regularization functions, an adaptive Huber functional based energy minimization model provides denoising with edge preservation in noisy digital images. We study a convergent finite difference scheme based on continuous piecewise linear functions and use a variable splitting scheme, namely the Split Bregman, to obtain the discrete minimizer. Experimental results are given in image denoising and comparison with additive operator splitting, dual fixed point, and projected gradient schemes illustrate that the best convergence rates are obtained for our algorithm.


title: Source Delay in Mobile Ad Hoc Networks

abstract: Source delay, the time a packet experiences in its source node, serves as a fundamental quantity for delay performance analysis in networks. However, the source delay performance in highly dynamic mobile ad hoc networks (MANETs) is still largely unknown by now. This paper studies the source delay in MANETs based on a general packet dispatching scheme with dispatch limit $f$ (PD-$f$ for short), where a same packet will be dispatched out up to $f$ times by its source node such that packet dispatching process can be flexibly controlled through a proper setting of $f$. We first apply the Quasi-Birth-and-Death (QBD) theory to develop a theoretical framework to capture the complex packet dispatching process in PD-$f$ MANETs. With the help of the theoretical framework, we then derive the cumulative distribution function as well as mean and variance of the source delay in such networks. Finally, extensive simulation and theoretical results are provided to validate our source delay analysis and illustrate how source delay in MANETs are related to network parameters.


title: Generalized Fitting subgroups of finite groups

abstract: In this paper we consider the Fitting subgroup $F(G)$ of a finite group $G$ and its generalizations: the quasinilpotent radical $F^*(G)$ and the generalized Fitting subgroup $\tilde{F}(G)$ defined by $\tilde{F}(G)\supseteq \Phi(G)$ and $\tilde{F}(G)/\Phi(G)=Soc(G/\Phi(G))$. We sum up known properties of $\tilde{F}(G)$ and suggest some new ones. Let $R$ be a subgroup of a group $G$. We shall call a subgroup $H$ of $G$ the $R$-subnormal subgroup if $H$ is subnormal in $ \langle H,R\rangle$. In this work the influence of $R$-subnormal subgroups (maximal, Sylow, cyclic primary) on the structure of finite groups are studied in the case when $R\in\{F(G), F^*(G),\tilde{F}(G)\}$.


title: Radial Transport of Large-Scale Magnetic Fields in Accretion Disks. I. Steady Solutions and an Upper Limit on the Vertical Field Strength

abstract: Large-scale magnetic fields are key ingredients of magnetically driven disk accretion. We study how large-scale poloidal fields evolve in accretion disks, with the primary aim of quantifying the viability of magnetic accretion mechanisms in protoplanetary disks. We employ a kinematic mean-field model for poloidal field transport and focus on steady states where inward advection of a field balances with outward diffusion due to effective resistivities. We analytically derive the steady-state radial distribution of poloidal fields in highly conducting accretion disks. The analytic solution reveals an upper limit on the strength of large-scale vertical fields attainable in steady states. Any excess poloidal field will be diffused away within a finite time, and we demonstrate this with time-dependent numerical calculations of the mean-field equations. We apply this upper limit to large-scale vertical fields threading protoplanetary disks. We find that the maximum attainable strength is about 0.1 G at 1 AU, and about 1 mG at 10 AU from the central star. When combined with recent magnetic accretion models, the maximum field strength translates into the maximum steady-state accretion rate of $\sim 10^{-7} M_\odot {\rm yr}^{-1}$, in agreement with observations. We also find that the maximum field strength is ~ 1 kG at the surface of the central star provided that the disk extends down to the stellar surface. This implies that any excess stellar poloidal field of strength >~ kG can be transported to the surrounding disk. This might in part resolve the magnetic flux problem in star formation.


title: Impulse Noise Removal In Speech Using Wavelets

abstract: A new method for removing impulse noise from speech in the wavelet transform domain is proposed. The method utilizes the multiresolution property of the wavelet transform, which provides finer time resolution at the higher frequencies than the short-time Fourier transform (STFT), to effectively identify and remove impulse noise. It uses two features of speech to discriminate speech from impulse noise: one is the slow time-varying nature of speech and the other is the Lipschitz regularity of the speech components. On the basis of these features, an algorithm has been developed to identify and suppress wavelet coefficients that correspond to impulse noise. Experiment results show that the new method is able to significantly reduce impulse noise without degrading the quality of the speech signal or introducing any audible artifacts.


title: An iterative algorithm for computed tomography image reconstruction from limited-angle projections

abstract: In application of tomography imaging, limited-angle problem is a quite practical and important issue. In this paper, an iterative reprojection-reconstruction (IRR) algorithm using a modified Papoulis-Gerchberg (PG) iterative scheme is developed for reconstruction from limited-angle projections which contain noise. The proposed algorithm has two iterative update processes, one is the extrapolation of unknown data, and the other is the modification of the known noisy observation data. And the algorithm introduces scaling factors to control the two processes, respectively. The convergence of the algorithm is guaranteed, and the method of choosing the scaling factors is given with energy constraints. The simulation result demonstrates our conclusions and indicates that the algorithm proposed in this paper can obviously improve the reconstruction quality.


title: Deadlock detection in linear recursive programs

abstract: Deadlock detection in recursive programs that admit dynamic resource creation is extremely complex and solutions either give imprecise answers or do not scale. We define an algorithm for detecting deadlocks of "linear recursive programs" of a basic model. The theory that underpins the algorithm is a generalization of the theory of permutations of names to so-called "mutations", which transform tuples by introducing duplicates and fresh names. Our algorithm realizes the back-end of deadlock analyzers for object-oriented programming languages, once the association programs/basic-model-programs has been defined as front-end.


title: Channel Flow of a Tensorial Shear-Thinning Maxwell Model: Lattice Boltzmann Simulations

abstract: We introduce a nonlinear generalized tensorial Maxwell-type constitutive equation to describe shear-thinning glass-forming fluids, motivated by a recent microscopic approach to the nonlinear rheology of colloidal suspensions. The model captures a nonvanishing dynamical yield stress at the glass transition and incorporates normal-stress differences. A modified lattice-Boltzmann (LB) simulation scheme is presented that includes non-Newtonian contributions to the stress tensor and deals with flow-induced pressure differences. We test this scheme in pressure-driven 2D Poiseuille flow of the nonlinear generalized Maxwell fluid. In the steady state, comparison with an analytical solution shows good agreement. The transient dynamics after startup and cessation of the pressure gradient are studied; the simulation reproduces a finite stopping time for the cessation flow of the yield-stress fluid in agreement with previous analytical estimates.


title: Emergence of pressure-induced metamagnetic-like state in Mn-doped CdGeAs2 chalcopyrite

abstract: The effect of hydrostatic pressure on resistivity and magnetic ac susceptibility has been studied in Mn-doped CdGeAs2 room-temperature (RT) ferromagnetic chalcopyrite with two types of MnAs micro-clusters. The slight increase of temperature by about 30 K in the region between RT and Curie temperature TC causes a significant change in the positions of pressure-induced semiconductor-metal transition and magnetic phase transitions in low pressure area. By conducting measurements of the anomalous Hall resistance in the field H \leq 5 kOe, we present experimental evidence for pressure-induced metamagnetic-like state during the paramagnetic phase at pressure P = 5 GPa.


title: How Much Should I Pay for Privacy Concerns in Truthful Online Crowd Sensing?

abstract: Crowd sensing is a new paradigm which leverages the pervasive smartphones to efficiently collect sensing data, enabling numerous novel applications. To achieve good service quality for a crowd sensing application, incentive mechanisms are indispensable to attract more user participation. Most of existing mechanisms only apply for the offline scenario, where the system has full information about the users' sensing profiles, i.e., a set of locations or mobility as well as the type of smartphones used, and their true costs. On the contrary, we focus on a more real scenario where users with their own privacy concerns arrive one by one online in a random order. We model the problem as a privacy-respecting online auction in which users are willing to negotiate access to certain private information and submit their sensing profiles satisfying privacy concerns to the platform (the provider of crowd sensing applications) over time, and the platform aims to the total total value of the services provided by selected users under a budget constraint. We then design two online mechanisms for a budgeted crowd sensing application, satisfying the computational efficiency, individual rationality, budget feasibility, truthfulness, consumer sovereignty, constant competitiveness and privacy concerns. Through extensive simulations, we evaluate the performance and validate the theoretical properties of our online mechanisms.


title: OutFlank Routing: Increasing Throughput in Toroidal Interconnection Networks

abstract: We present a new, deadlock-free, routing scheme for toroidal interconnection networks, called OutFlank Routing (OFR). OFR is an adaptive strategy which exploits non-minimal links, both in the source and in the destination nodes. When minimal links are congested, OFR deroutes packets to carefully chosen intermediate destinations, in order to obtain travel paths which are only an additive constant longer than the shortest ones. Since routing performance is very sensitive to changes in the traffic model or in the router parameters, an accurate discrete-event simulator of the toroidal network has been developed to empirically validate OFR, by comparing it against other relevant routing strategies, over a range of typical real-world traffic patterns. On the 16x16x16 (4096 nodes) simulated network OFR exhibits improvements of the maximum sustained throughput between 14% and 114%, with respect to Adaptive Bubble Routing.


title: Comparison of Gravitational Wave Detector Network Sky Localization Approximations

abstract: Gravitational waves emitted during compact binary coalescences are a promising source for gravitational-wave detector networks. The accuracy with which the location of the source on the sky can be inferred from gravitational wave data is a limiting factor for several potential scientific goals of gravitational-wave astronomy, including multi-messenger observations. Various methods have been used to estimate the ability of a proposed network to localize sources. Here we compare two techniques for predicting the uncertainty of sky localization -- timing triangulation and the Fisher information matrix approximations -- with Bayesian inference on the full, coherent data set. We find that timing triangulation alone tends to over-estimate the uncertainty in sky localization by a median factor of $4$ for a set of signals from non-spinning compact object binaries ranging up to a total mass of $20 M_\odot$, and the over-estimation increases with the mass of the system. We find that average predictions can be brought to better agreement by the inclusion of phase consistency information in timing-triangulation techniques. However, even after corrections, these techniques can yield significantly different results to the full analysis on specific mock signals. Thus, while the approximate techniques may be useful in providing rapid, large scale estimates of network localization capability, the fully coherent Bayesian analysis gives more robust results for individual signals, particularly in the presence of detector noise.


title: The Rigorous Stochastic Matrix Multiplication Scheme for the Calculations of Reduced Equilibrium Density Matrices of Open Multilevel Quantum Systems

abstract: Understanding the roles of the temporary and spatial structures of quantum functional noise in open multilevel quantum molecular systems attracts a lot of theoretical interests. I want to establish a rigorous and general framework for functional quantum noises from the constructive and computational perspectives, $\it{i.e.}$ how to generate the random trajectories to reproduce the kernel and path ordering of the influence functional with effective Monte Carlo methods for arbitrary spectral densities. This construction approach aims to unify the existing stochastic models to rigorously describe the temporary and spatial structure of Gaussian quantum noises. In this paper, I review the Euclidean imaginary time influence functional and propose the stochastic matrix multiplication scheme to calculate reduced equilibrium density matrices (REDM). In addition, I review and discuss the Feynman-Vernom influence functional according to the Gaussian quadratic integral, particularly its imaginary part which is critical to the rigorous description of the quantum detailed balance. As a result, I establish the conditions under which the influence functional can be interpreted as the average of exponential functional operator over real-valued Gaussian processes for open multilevel quantum systems. I also show the difference between the local and nonlocal phonons within this framework. With the stochastic matrix multiplication scheme, I compare the normalized REDM with the Boltzmann equilibrium distribution for open multilevel quantum systems.


title: Error Analysis of Nuclear Matrix Elements

abstract: We estimate the expected errors of nuclear matrix elements coming from the uncertainty on the NN interaction. We use a coarse grained (GR) interaction fitted to NN scattering data, with several prescriptions for the long-part of the interaction, including one pion exchange and chiral two-pion exchange interactions.


title: Measurements without Probabilities in the Final State Proposal

abstract: The black hole final state proposal reconciles the infalling vacuum with the unitarity of the Hawking radiation, but only for some experiments. We study experiments that first verify the exterior, then the interior purification of the same Hawking particle. (This is the same protocol that renders the firewall paradox operationally meaningful in standard quantum mechanics.) We show that the decoherence functional fails to be diagonal, even upon inclusion of external "pointer" systems. Hence, probabilities for outcomes of these measurements are not defined. We conclude that the final state proposal does not offer a consistent alternative to the firewall hypothesis.


title: Random walks generated by equilibrium contact processes

abstract: We consider dynamic random walks where the nearest neighbour jump rates are determined by an underlying supercritical contact process in equilibrium. This has previously been studied by den Hollander and dos Santos and den Hollander, dos Santos, Sidoravicius. We show the CLT for such a random walk, valid for all supercritical infection rates for the contact process environment.


title: Second Harmonic Generation in h-BN and MoS$_2$ monolayers: the role of electron-hole interaction

abstract: In this letter we show by means of first principle numerical simulations that electron-hole interaction significantly contributes to the second-harmonic generation spectrum of h-BN or MoS$_2$ monolayers. Specifically, it doubles the signal intensity at the excitonic resonances with respect to the contribution from independent electronic transitions. This result hints that the intensity of second-harmonic signal of those materials can be tuned by changing the dielectric screening, that controls the strength of the electron-hole interaction.


title: Direct optical sensing of single unlabeled small proteins and super-resolution microscopy of their binding sites

abstract: More than twenty years ago, scientists succeeded in pushing the limits of optical detection to single molecules using fluorescence. This breakthrough has revolutionized biophysical measurements, but restrictions in photophysics and labeling protocols have motivated many efforts to achieve fluorescence-free single-molecule sensitivity in biological studies. Although several interesting mechanisms using vibrational spectroscopy, photothermal detection, plasmonics or microcavities have been proposed for biosensing at the single-protein level, no method has succeeded in direct label-free detection of single proteins. Here, we present the first results using interferometric detection of scattering (iSCAT) from single proteins without the need for any label, optical nanostructure or microcavity. Furthermore, we demonstrate super-resolution imaging of protein binding with nanometer localization precision. The ease of iSCAT instrumentation promises a breakthrough for industrial usage as well as fundamental laboratory experiments.


title: Tuneable drug-loading capability of chitosan hydrogels with varied network architectures

abstract: Advanced bioactive systems with defined macroscopic properties and spatio-temporal sequestration of extracellular biomacromolecules are highly desirable for next generation therapeutics. Here, chitosan hydrogels were prepared with neutral or negatively-charged crosslinkers in order to promote selective electrostatic complexation with charged drugs. Chitosan (CT) was functionalised with varied dicarboxylic acids, such as tartaric acid (TA), poly(ethylene glycol) bis(carboxymethyl) ether (PEG), 1.4-Phenylenediacetic acid (4Ph) and 5-Sulfoisophthalic acid monosodium salt (PhS), whereby PhS was hypothesised to act as a simple mimetic of heparin. ATR FT-IR showed the presence of C=O amide I, N-H amide II and C=O ester bands, providing evidence of covalent network formation. The crosslinker content was reversely quantified by 1H-NMR on partially-degraded network oligomers, so that 18 mol% PhS was exemplarily determined. Swellability, compressability, material morphology, and drug-loading capability were successfully adjusted based on the selected network architecture. Here, hydrogel incubation with model drugs of varied electrostatic charge, i.e. allura red (AR, --), methyl orange (MO, -) or methylene blue (MB, +), resulted in direct hydrogel-dye electrostatic complexation. Importantly, the cationic compound, MB, showed different incorporation behaviours, depending on the electrostatic character of the selected crosslinker. In light of this tuneable drug-loading capability, these CT hydrogels would be highly attractive as drug reservoirs towards e.g. the fabrication of tissue models in vitro.


title: Asymptotic Properties of Bayes Risk of a General Class of Shrinkage Priors in Multiple Hypothesis Testing Under Sparsity

abstract: Consider the problem of simultaneous testing for the means of independent normal observations. In this paper, we study some asymptotic optimality properties of certain multiple testing rules induced by a general class of one-group shrinkage priors in a Bayesian decision theoretic framework, where the overall loss is taken as the number of misclassified hypotheses. We assume a two-groups normal mixture model for the data and consider the asymptotic framework adopted in Bogdan et al. (2011) who introduced the notion of asymptotic Bayes optimality under sparsity in the context of multiple testing. The general class of one-group priors under study is rich enough to include, among others, the families of three parameter beta, generalized double Pareto priors, and in particular the horseshoe, the normal-exponential-gamma and the Strawderman-Berger priors. We establish that within our chosen asymptotic framework, the multiple testing rules under study asymptotically attain the risk of the Bayes Oracle up to a multiplicative factor, with the constant in the risk close to the constant in the Oracle risk. This is similar to a result obtained in Datta and Ghosh (2013) for the multiple testing rule based on the horseshoe estimator introduced in Carvalho et al. (2009, 2010). We further show that under very mild assumption on the underlying sparsity parameter, the induced decision rules based on an empirical Bayes estimate of the corresponding global shrinkage parameter proposed by van der Pas et al. (2014), attain the optimal Bayes risk up to the same multiplicative factor asymptotically. We provide a unifying argument applicable for the general class of priors under study. In the process, we settle a conjecture regarding optimality property of the generalized double Pareto priors made in Datta and Ghosh (2013). Our work also shows that the result in Datta and Ghosh (2013) can be improved further.


title: Simulating cosmic reionization: How large a volume is large enough?

abstract: We present the largest-volume (425 Mpc/h=607 Mpc on a side) full radiative transfer simulation of cosmic reionization to date. We show that there is significant additional power in density fluctuations at very large scales. We systematically investigate the effects this additional power has on the progress, duration and features of reionization, as well as on selected reionization observables. We find that comoving simulation volume of ~100 Mpc/h per side is sufficient for deriving a convergent mean reionization history, but that the reionization patchiness is significantly underestimated. We use jackknife splitting to quantify the convergence of reionization properties with simulation volume for both mean-density and variable-density sub-regions. We find that sub-volumes of ~100 Mpc/h per side or larger yield convergent reionization histories, except for the earliest times, but smaller volumes of ~50 Mpc/h or less are not well converged at any redshift. Reionization history milestones show significant scatter between the sub-volumes, of Delta z=0.6-1 for ~50 Mpc/h volumes, decreasing to Delta z=0.3-0.5 for ~100 Mpc/h volumes, and $\Delta z$~0.1 for ~200 Mpc/h volumes. If we only consider mean-density sub-regions the scatter decreases, but remains at Delta z~0.1-0.2 for the different size sub-volumes. Consequently, many potential reionization observables like 21-cm rms, 21-cm PDF skewness and kurtosis all show good convergence for volumes of ~200 Mpc/h, but retain considerable scatter for smaller volumes. In contrast, the three-dimensional 21-cm power spectra at large scales (k<0.25 h/Mpc) do not fully converge for any sub-volume size. These additional large-scale fluctuations significantly enhance the 21-cm fluctuations, which should improve the prospects of detection considerably, given the lower foregrounds and greater interferometer sensitivity at higher frequencies. (abridged)


title: Asteroseismology of DAV star HS0507+0434B, including the core composition profiles

abstract: The DAV star, HS0507+0434B, was observed by Fu et al. (2013) in 2007 and from 2009 December to 2010 January. There were a total of six triplets with nearly equal split, which were identified as $l$ = 1 modes caused by rotation. In order to fit the six $l$ = 1 modes, grids of white dwarf models are generated by WDEC. For the core composition profiles, we choose the linear fittings to carbon profile of white dwarf models from MESA, which can be considered as results of real nuclear burning process. Coupled with diffusion equilibrium H/He and He/C mixtures, we make grids of models in WDEC and do asteroseismology works on HS0507+0434B. There is a total of 9.50 seconds error for our best-fitting model, which is smaller than the result (a total of 22.1 seconds error) of Fu et al. (2013), when fitting the six $l$ = 1 modes. The two other previous identified pulsation mode frequencies (286.1 s and 743.40 s) may also be well fitted by our best-fitting model. The model parameters are $T_{eff}$ = 11450 K, log$g$ = 8.088, $M_{*}$ = 0.640 $M_{\odot}$, log($M_{H}$/$M_{*}$) = -6, and log($M_{He}$/$M_{*}$) = -3.


title: Pion electromagnetic form factor in the Covariant Spectator Theory

abstract: The pion electromagnetic form factor at spacelike momentum transfer is calculated in relativistic impulse approximation using the Covariant Spectator Theory. The same dressed quark mass function and the equation for the pion bound-state vertex function as discussed in the companion paper are used for the calculation, together with a dressed quark current that satisfies the Ward-Takahashi identity. The results obtained for the pion form factor are in agreement with experimental data, they exhibit the typical monopole behavior at high-momentum transfer, and they satisfy some remarkable scaling relations.


title: Transport coefficients of the QGP

abstract: The FASTSUM collaboration presents a study on the temperature dependence of the electrical conductivity $\sigma$ in the quark-gluon plasma, using the methods of lattice QCD. Correlators of the exactly conserved vector current are measured at different temperatures across the deconfinement transition, using ensembles of $2+1$ flavours of dynamical fermions on anisotropic lattices. We then employ bayesian methods (MEM) to extract the relevant spectral functions, which are found to be consistent with $\sigma/T$ rising as a function of $T$. The robustness of the results is verified by a detailed analysis of the systematics involved in the bayesian reconstruction of the spectral functions.


title: Hypothesis Testing for Topological Data Analysis

abstract: Persistent homology is a vital tool for topological data analysis. Previous work has developed some statistical estimators for characteristics of collections of persistence diagrams. However, tools that provide statistical inference for observations that are persistence diagrams are limited. Specifically, there is a need for tests that can assess the strength of evidence against a claim that two samples arise from the same population or process. We propose the use of randomization-style null hypothesis significance tests (NHST) for these situations. The test is based on a loss function that comprises pairwise distances between the elements of each sample and all the elements in the other sample. We use this method to analyze a range of simulated and experimental data. Through these examples we experimentally explore the power of the p-values. Our results show that the randomization-style NHST based on pairwise distances can distinguish between samples from different processes, which suggests that its use for hypothesis tests upon persistence diagrams is reasonable. We demonstrate its application on a real dataset of fMRI data of patients with ADHD.


title: The Real Dynamics of Bieberbach's Example

abstract: Bieberbach constructed in 1933 domains in $ \bf {C}^2$ which were biholomorphic to $ \bf {C}^2$ but not dense. The existence of such domains was unexpected. The special domains Bieberbach considered are basins of attraction of a cubic H\'enon map. This classical method of construction is one of the first applications of dynamical systems to complex analysis. In this paper, the boundaries of the real sections of Bieberbach's domains will be calculated explicitly as the stable manifolds of the saddle points. The real filled Julia sets and the real Julia sets of Bieberbach's map will also be calculated explicitly and illustrated with computer generated graphics. Basic differences between real and the complex dynamics will be shown.


title: Mining the Temporal Evolution of the Android Bug Reporting Community via Sliding Windows

abstract: The open source development community consists of both paid and volunteer developers as well as new and experienced users. Previous work has applied social network analysis (SNA) to open source communities and has demonstrated value in expertise discovery and triaging. One problem with applying SNA directly to the data of the entire project lifetime is that the impact of local activities will be drowned out. In this paper we provide a method for aggregating, analyzing, and visualizing local (small time periods) interactions of bug reporting participants by using the SNA to measure the betweeness centrality of these participants. In particular we mined the Android bug repository by producing social networks from overlapping 30-day windows of bug reports, each sliding over by day. In this paper we define three patterns of participant behaviour based on their local centrality. We propose a method of analyzing the centrality of bug report participants both locally and globally, then we conduct a thorough case study of the bug reporter's activity within the Android bug repository. Furthermore, we validate the conclusions of our method by mining the Android version control system and inspecting the Android release history. We found that windowed SNA analysis elicited local behaviour that were invisible during global analysis.


title: Conductance through disordered graphene nanoribbons: Standard and anomalous electron localization

abstract: Conductance fluctuations produced by the presence of disorder in zigzag and armchair graphene nanoribbons are studied. We show that quantum transport in zigzag nanoribbons takes place via edge states which are exponentially localized, as in the standard Anderson localization problem, whereas for armchair nanoribbons the symmetry of the graphene sublattices produces anomalous localization, or delocalization. We show that these two different electron localizations lead to significant differences of the conductance statistics between zigzag and armchair nanoribbons. In particular, armchair nanoribbons show nonconventional large conductance fluctuations relative to those of Anderson-localized electrons. We calculate analytically the complete distribution of conductances for both types of ribbons. Without free fitting parameters, we verify our theoretical results by performing numerical simulations of disordered zigzag and armchair nanoribbons of experimentally achievable lengths and widths.


title: One-Loop Matching for Parton Distributions: Non-Singlet Case

abstract: We derive one-loop matching condition for non-singlet quark distributions in transverse-momentum cut-off scheme, including unpolarized, helicity and transversity distributions. The matching is between the quasi-distribution defined by static correlation at finite nucleon momentum and the light-cone distribution measurable in experiments. The result is useful for extracting the latter from the former in a lattice QCD calculation.


title: Semistability of Principal Bundles on a K\"ahler Manifold with a Non-Connected Structure Group

abstract: We investigate principal $G$-bundles on a compact K\"ahler manifold, where $G$ is a complex algebraic group such that the connected component of it containing the identity element is reductive. Defining (semi)stability of such bundles, it is shown that a principal $G$-bundle $E_G$ admits an Einstein-Hermitian connection if and only if $E_G$ is polystable. We give an equivalent formulation of the (semi)stability condition. A question is to compare this definition with that of math.AG/0506511.


title: Connectivity of confined 3D Networks with Anisotropically Radiating Nodes

abstract: Nodes in ad hoc networks with randomly oriented directional antenna patterns typically have fewer short links and more long links which can bridge together otherwise isolated subnetworks. This network feature is known to improve overall connectivity in 2D random networks operating at low channel path loss. To this end, we advance recently established results to obtain analytic expressions for the mean degree of 3D networks for simple but practical anisotropic gain profiles, including those of patch, dipole and end-fire array antennas. Our analysis reveals that for homogeneous systems (i.e. neglecting boundary effects) directional radiation patterns are superior to the isotropic case only when the path loss exponent is less than the spatial dimension. Moreover, we establish that ad hoc networks utilizing directional transmit and isotropic receive antennas (or vice versa) are always sub-optimally connected regardless of the environment path loss. We extend our analysis to investigate boundary effects in inhomogeneous systems, and study the geometrical reasons why directional radiating nodes are at a disadvantage to isotropic ones. Finally, we discuss multi-directional gain patterns consisting of many equally spaced lobes which could be used to mitigate boundary effects and improve overall network connectivity.

categories: cs.IT cond-mat.dis-nn cond-mat.stat-mech cs.NI math.IT


title: Cyclotomic Carter-Payne homomorphisms

abstract: We construct a new family of homomorphisms between (graded) Specht modules of the quiver Hecke algebras of type A. These maps have many similarities with the homomorphisms constructed by Carter and Payne in the special case of the symmetric groups, although the maps that we obtain are both more and less general than these.


title: Long-Lived Vortex Structures in Collisional Pure and Gas-Discharge Nonneutral Electron Plasmas

abstract: The analysis of experimental investigations of equilibrium, interaction and dynamics of vortex structures in pure electron and gas-discharge electron nonneutral plasmas during the time much more than the electron-neutral collision time has been carried out. The problem of long confinement of the column of pure electron plasma in Penning-Malmberg trap is considered. The mechanism of stability of long-lived vortex structure in gas-discharge nonneutral electron plasma is investigated. The collapse of electron sheath in gas-discharge nonneutral electron plasma in Penning cell at high pressures of neutral gas is described. The interaction between the stable vortex structure and the annular electron sheath, and the action of vortex structures on the transport of electrons along and across the magnetic field are discussed.


title: Strongly Koszul edge rings

abstract: We classify the finite connected simple graphs whose edge rings are strongly Koszul. From the classification, it follows that if the edge ring is strongly Koszul, then its toric ideal possesses a quadratic Gr\"obner basis.


title: Non-commutative integration, zeta functions and the Haar state for $SU_{q}(2)$

abstract: We study a notion of non-commutative integration, in the spirit of modular spectral triples, for the quantum group $SU_{q}(2)$. In particular we define the non-commutative integral as the residue at the spectral dimension of a zeta function, which is constructed using a Dirac operator and a weight. We consider the Dirac operator introduced by Kaad and Senior and a family of weights depending on two parameters, which are related to the diagonal automorphisms of $SU_{q}(2)$. We show that, after fixing one of the parameters, the non-commutative integral coincides with the Haar state of $SU_{q}(2)$. Moreover we can impose an additional condition on the zeta function, which also fixes the second parameter. For this unique choice the spectral dimension coincides with the classical dimension.


title: Phenomenology of neutrino oscillations and mixing

abstract: We review the status of three-neutrino mixing and the results of global analyses of short-baseline neutrino oscillation data in 3+1 and 3+2 neutrino mixing schemes.


title: Weak Convergence Rates of Population versus Single-Chain Stochastic Approximation MCMC Algorithms

abstract: In this paper, we establish the theory of weak convergence (toward a normal distribution) for both single-chain and population stochastic approximation MCMC algorithms. Based on the theory, we give an explicit ratio of convergence rates for the population SAMCMC algorithm and the single-chain SAMCMC algorithm. Our results provide a theoretic guarantee that the population SAMCMC algorithms are asymptotically more efficient than the single-chain SAMCMC algorithms when the gain factor sequence decreases slower than O(1/t), where t indexes the number of iterations. This is of interest for practical applications.


title: New identities for the polarized partitions and partitions with $d$-distant parts

abstract: In this paper we present a new class of integer partition identities. The number of partitions with d-distant parts can be represented as a sum of the number of partitions with 1-distant parts whose even parts are greater than twice the number of odd parts. We also provide a direct bijection between these classes of partitions.


title: McMullen polynomials and Lipschitz flows for free-by-cyclic groups

abstract: Consider a group G and an epimorphism u_0:G\to\Z inducing a splitting of G as a semidirect product ker(u_0)\rtimes_\varphi\Z with ker(u_0) a finitely generated free group and \varphi\in Out(ker(u_0)) representable by an expanding irreducible train track map. Building on our earlier work [Dynamics on free-by-cyclic groups, arXiv:1301.7739], in which we we realized G as \pi_1(X) for an Eilenberg-Maclane 2-complex X equipped with a semiflow \psi, and inspired by McMullen's Teichm\"uller polynomial for fibered hyperbolic 3-manifolds, we construct a polynomial invariant \m for (X,\psi) and investigate its properties. Specifically, \m determines a convex polyhedral cone \C_X in H^1(G;\R), a convex, real-analytic function \H:\C_X\to\R, and specializes to give an integral Laurent polynomial \m_u(\zeta) for each integral u\in\C_X. We show that \C_X is equal to the "cone of sections" of (X,\psi) (the convex hull of all cohomology classes dual to sections of of \psi), and that for each (compatible) cross section \Theta_u with first return map f_u:\Theta_u\to\Theta_u, the specialization \m_u(\zeta) encodes the characteristic polynomial of the transition matrix of f_u. More generally, for every class u\in\C_X there exists a geodesic metric d_u and a codimension-1 foliation \Omega_u of X transverse to \psi so that after reparametrizing the flow \psi^u_s maps leaves of \Omega_u to leaves via a local e^{s\H(u)}-homothety. Among other things, we additionally prove that \C_X is equal to (the cone over) the component of the BNS-invariant containing u_0 and that each primitive integral u\in\C_X induces a splitting of G as an ascending HNN-extension over a finite-rank free group along an injective endomorphism \phi_u. For any such splitting, we show that the stretch factor of \phi_u is exactly given by e^{\H(u)}. In particular, we see that \C_X and \H depend only on the group G and epimorphism u_0.


title: Primordial non-Gaussianity in the bispectra of large-scale structure

abstract: The statistics of large-scale structure in the Universe can be used to probe non-Gaussianity of the primordial density field, complementary to existing constraints from the cosmic microwave background. In particular, the scale dependence of halo bias, which affects the halo distribution at large scales, represents a promising tool for analyzing primordial non-Gaussianity of local form. Future observations, for example, may be able to constrain the trispectrum parameter gNL that is difficult to study and constrain using the CMB alone. We investigate how galaxy and matter bispectra can distinguish between the two non-Gaussian parameters fNL and gNL, whose effects give nearly degenerate contributions to the power spectra. We use a generalization of the univariate bias approach, making the hypothesis that the number density of halos forming at a given position is a function of the local matter density contrast and of its local higher-order statistics. Using this approach, we calculate the halo-matter bispectra and analyze their properties. We determine a connection between the sign of the halo bispectrum on large scales and the parameter gNL. We also construct a combination of halo and matter bispectra that is sensitive to fNL, with little contamination from gNL. We study both the case of single and multiple sources to the primordial gravitational potential, discussing how to extend the concept of stochastic halo bias to the case of bispectra. We use a specific halo mass-function to calculate numerically the bispectra in appropriate squeezed limits, confirming our theoretical findings.


title: Confirmation of the planet around HD 95086 by direct imaging

abstract: VLT/NaCo angular differential imaging at L' (3.8 microns) revealed a probable giant planet comoving with the young and early-type HD 95086 also known to harbor an extended debris disk. The discovery was based on the proper motion analysis of two datasets spanning 15 months. However, the second dataset suffered from bad atmospheric conditions, which limited the significance of the redetection at the 3 sigma level. In this Letter, we report new VLT/NaCo observations of HD 95086 obtained on 2013 June 26-27 at L' to recover the planet candidate. We unambiguously redetect the companion HD 95086 b with multiple independent pipelines at a signal-to-noise ratio greater than or equal to 5. Combined with previously reported measurements, our astrometry decisively shows that the planet is comoving with HD 95086 and inconsistent with a background object. With a revised mass of 5 pm 2 Jupiter masses, estimated from its L' photometry and "hot-start" models at 17 pm 4 Myr, HD 95086 b becomes a new benchmark for further physical and orbital characterization of young giant planets.


title: Multigraded Betti numbers of simplicial forests

abstract: We prove that multigraded Betti numbers of a simplicial forest are always either 0 or 1. Moreover a nonzero multidegree appears exactly at one homological degree in the resolution. Our work generalizes work of Bouchat on edge ideals of graph forests.


title: Color Differences between Clockwise and Counterclockwise Spiral Galaxies

abstract: While spiral galaxies observed from Earth clearly seem to spin in different directions, little is yet known about other differences between galaxies that spin clockwise and galaxies that spin counterclockwise. Here we compared the color of 64,399 spiral galaxies that spin clockwise to 63,215 spiral galaxies that spin counterclockwise. The results show that clockwise galaxies tend to be bluer than galaxies that spin counterclockwise. The probability that the color differences can be attributed to chance is ~0.019.


title: Quantum walk with a general coin: Exact solution and asymptotic properties

abstract: In this paper we present closed-form expressions for the wave function that governs the evolution of the discrete-time quantum walk on a line when the coin operator is arbitrary. The formulas were derived assuming that the walker can either remain put in the place or proceed in a fixed direction but never move backward, although they can be easily modified to describe the case in which the particle can travel in both directions. We use these expressions to explore the properties of magnitudes associated to the process, as the probability mass function or the probability current, even though we also consider the asymptotic behavior of the exact solution. Within this approximation, we will estimate upper and lower bounds, consider the origins of an emerging approximate symmetry, and deduce the general form of the stationary probability density of the relative location of the walker.

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: On the quantum stability of Q-balls

abstract: We consider the evolution and decay of Q-balls under the influence of quantum fluctuations. We argue that the most important effect resulting from these fluctuations is the modification of the effective potential in which the Q-ball evolves. This is in addition to spontaneous decay into elementary particle excitations and fission into smaller Q-balls previously considered in the literature, which -- like most tunnelling processes -- are likely to be strongly suppressed. We illustrate the effect of quantum fluctuations in a particular model $\phi^6$ potential, for which we implement the inhomogeneous Hartree approximation to quantum dynamics and solve for the evolution of Q-balls in 3+1 dimensions. We find that the stability range as a function of (field space) angular velocity $\omega$ is modified significantly compared to the classical case, so that small-$\omega$ Q-balls are less stable than in the classical limit, and large-$\omega$ Q-balls are more stable. This can be understood qualitatively in a simple way.


title: Turbulence patterns and neutrino flavor transitions in high-resolution supernova models

abstract: During the shock-wave propagation in a core-collapse supernova (SN), matter turbulence may affect neutrino flavor conversion probabilities. Such effects have been usually studied by adding parametrized small-scale random fluctuations (with arbitrary amplitude) on top of coarse, spherically symmetric matter density profiles. Recently, however, two-dimensional (2D) SN models have reached a space resolution high enough to directly trace anisotropic density profiles, down to scales smaller than the typical neutrino oscillation length. In this context, we analyze the statistical properties of a large set of SN matter density profiles obtained in a high-resolution 2D simulation, focusing on a post-bounce time (2 s) suited to study shock-wave effects on neutrino propagation on scales as small as O(100) km and possibly below. We clearly find the imprint of a broken (Kolmogorov-Kraichnan) power-law structure, as generically expected in 2D turbulence spectra. We then compute the flavor evolution of SN neutrinos along representative realizations of the turbulent matter density profiles, and observe no or modest damping of the neutrino crossing probabilities on their way through the shock wave. In order to check the effect of possibly unresolved fluctuations at scales below O(100) km, we also apply a randomization procedure anchored to the power spectrum calculated from the simulation, and find consistent results within \pm 1 sigma fluctuations. These results show the importance of anchoring turbulence effects on SN neutrinos to realistic, fine-grained SN models.


title: Parallel Magnetic Resonance Imaging as Approximation in a Reproducing Kernel Hilbert Space

abstract: In Magnetic Resonance Imaging (MRI) data samples are collected in the spatial frequency domain (k-space), typically by time-consuming line-by-line scanning on a Cartesian grid. Scans can be accelerated by simultaneous acquisition of data using multiple receivers (parallel imaging), and by using more efficient non-Cartesian sampling schemes. As shown here, reconstruction from samples at arbitrary locations can be understood as approximation of vector-valued functions from the acquired samples and formulated using a Reproducing Kernel Hilbert Space (RKHS) with a matrix-valued kernel defined by the spatial sensitivities of the receive coils. This establishes a formal connection between approximation theory and parallel imaging. Theoretical tools from approximation theory can then be used to understand reconstruction in k-space and to extend the analysis of the effects of samples selection beyond the traditional g-factor noise analysis to both noise amplification and approximation errors. This is demonstrated with numerical examples.


title: Crystallizing hard-sphere glasses by doping with active particles

abstract: Crystallization and vitrification are two different routes to form a solid. Normally these two processes suppress each other, with the glass transition preventing crystallization at high density (or low temperature). This is even true for systems of colloidal hard spheres, which are commonly used as building blocks for novel functional materials with potential applications, e.g. photonic crystals. By performing Brownian dynamics simulations of glassy systems consisting of mixtures of active and passive hard spheres, we show that the crystallization of such hard-sphere glasses can be dramatically promoted by doping the system with small amounts of active particles. Surprisingly, even hard-sphere glasses of packing fraction up to $\phi = 0.635$ crystallize, which is around $0.5\%$ below the random close packing at $\phi \simeq 0.64$. Our results suggest a novel way of fabricating crystalline materials from (colloidal) glasses. This is particularly important for materials that get easily kinetically trapped in glassy states, and crystal nucleation hardly occurs.


title: Large enhancement of Forster resonance energy transfer on graphene platforms

abstract: In the view of the applications of Forster resonant energy transfer (FRET) in biological systems which especially require FRET in the inrared region we investigate the great advantage of graphene plasmonics in such studies. Focusing on the fundamental aspects of FRET between a donor-acceptor pair on a graphene platform showing that FRET mediated by the plasmons in graphene is broadband and enhanced by six orders of magnitude. We briefly discuss the impact of phonon-polaritonic substrates.


title: A coronal wave as CME footprint

abstract: I report a coronal wave observed in Fexii, soft X-ray, and Halpha on November 6, 2006, jointly with a CME and a flare to show the spatial and temporal relation of the wave and the CME. I also take advantage of the spectral resolution of the wave to obtain an approximation of its temperature. Finally, I compare the magnetic field topology to the wave location. The observations in two band passes show bright fronts that are co-spatial when observed at the same time, therefore I believe that the same wave is observed in different band passes. The wave front is not observable in Halpha, but just in Fexii and in soft X-ray, indicating that the temperature of the wave front is slightly higher than the chromospheric temperature. The ratio of emission through two SXI filters gives a temperature of 7 10^6 K for the wave front. The northern-most edge of the CME footprint related to this wave corresponds to the last location of the wave on the limb. The image processing reveals stationary brightenings produced on the passage of the wave front. The potential spherical extrapolation of the sun shows that those stationary brightenings are lying in jumps of magnetic field lines connectivity.


title: The merging cluster of galaxies Abell 3376: an optical view

abstract: Abell 3376 is a merging cluster of galaxies at redshift z=0.046, famous mostly for its giant radio arcs, and shows an elongated and highly substructured X-ray emission, but has not been analysed in detail at optical wavelengths. We have obtained wide field images of Abell 3376 in the B band and derive the GLF applying a statistical subtraction of the background in three regions: a circle of 0.29 deg radius (1.5 Mpc) encompassing the whole cluster, and two circles centered on each of the two brightest galaxies (BCG2, northeast, coinciding with the peak of X-ray emission, and BCG1, southwest) of radii 0.15 deg (0.775 Mpc). We also compute the GLF in the zone around BCG1, which is covered by the WINGS survey in the B and V bands, by selecting cluster members in the red sequence in a (B-V) versus V diagram. Finally, we discuss the dynamical characteristics of the cluster implied by a Serna & Gerbal analysis. The GLFs are not well fit by a single Schechter function, but satisfactory fits are obtained by summing a Gaussian and a Schechter function. The GLF computed by selecting galaxies in the red sequence in the region surrounding BCG1 can also be fit by a Gaussian plus a Schechter function. An excess of galaxies in the brightest bins is detected in the BCG1 and BCG2 regions. The dynamical analysis based on the Serna & Gerbal method shows the existence of a main structure of 82 galaxies which can be subdivided into two main substructures of 25 and 6 galaxies. A smaller structure of 6 galaxies is also detected. The B band GLFs of Abell 3376 are clearly perturbed, as already found in other merging clusters. The dynamical properties are consistent with the existence of several substructures, in agreement with a previously published X-ray analysis.


title: Isaac Newton's sinister heraldry

abstract: After Isaac Newton was knighted by Queen Anne in 1705 he adopted an unusual coat of arms: a pair of human tibiae crossed on a black background, like a pirate flag without the skull. After some general reflections on Newton's monumental scientific achievements and on his enigmatic life, we investigate the story of his coat of arms. We also discuss how its simple design illustrates the concept of chirality, which would later play an important role in the philosophical arguments about Newton's conception of space, as well as in the development of modern chemistry and particle physics.


title: Hypercharged Dark Matter and Direct Detection as a Probe of Reheating

abstract: The lack of new physics at the LHC so far weakens the argument for TeV scale thermal dark matter. On the other hand, heavier, non-thermal dark matter is generally difficult to test experimentally. Here we consider the interesting and generic case of hypercharged dark matter, which can allow for heavy dark matter masses without spoiling testability. Planned direct detection experiments will be able to see a signal for masses up to an incredible $10^{10}$ GeV, and this can further serve to probe the reheating temperature up to about $10^9$ GeV, as determined by the non-thermal dark matter relic abundance. The $Z$-mediated nature of the dark matter scattering may be determined in principle by comparing scattering rates on different detector nuclei, which in turn can reveal the dark matter mass. We will discuss the extent to which future experiments may be able to make such a determination.


title: Recollision scenario without tunneling : Role of the ionic core potential

abstract: The standard model of strong laser physics, the recollision scenario, omits the ionic core potential after tunneling. Strikingly, although the Coulomb interaction drives all stages of recollision, the maximum energy the electrons bring back to the core is found by ignoring it. We resolve this long-standing paradox by showing that this good agreement stems from a fortuitous cancellation at high intensities. Instead of the three step model, we find that the Coulomb interaction can be fully integrated into a purely classical scenario that explains recollisions without invoking tunneling.


title: Shock-induced $\mathcal{PT}$-symmetric potentials in gas-filled photonic crystal fibers

abstract: We have investigated the interaction between a strong soliton and a weak probe with certain configurations that allow optical trapping in gas-filled hollow-core photonic crystal fibers in the presence of the shock effect. We have shown theoretically and numerically that the shock term can lead to an unbroken parity-time $ \left(\mathcal{PT}\right) $ symmetry potential in these kinds of fibers. Reciprocity breaking, a remarkable feature of the $ \mathcal{PT} $ symmetry, is also demonstrated numerically. Our results will open different configurations and avenues for observing $ \mathcal{PT} $-symmetry breaking in optical fibers, without the need to resort to cumbersome dissipative structures.


title: On the Partial-Wave Analysis of Mesonic Resonances Decaying to Multiparticle Final States Produced by Polarized Photons

abstract: Meson spectroscopy is going through a revival with the advent of high statistics experiments and new advances in the theoretical predictions. The Constituent Quark Model (CQM) is finally being expanded considering more basic principles of field theory and using discrete calculations of Quantum Chromodynamics (lattice QCD). These new calculations are approaching predictive power for the spectrum of hadronic resonances and decay modes. It will be the task of the new experiments to extract the meson spectrum from the data and compare with those predictions. The goal of this report is to describe one particular technique for extracting resonance information from multiparticle final states. The technique described here, partial wave analysis based on the helicity formalism, has been used at Brookhaven National Laboratory (BNL) using pion beams, and Jefferson Laboratory (JLab) using photon beams. In particular this report broadens this technique to include production experiments using linearly polarized real photons or quasi-real photons. This article is of a didactical nature. We describe the process of analysis, detailing assumptions and formalisms, and is directed towards people interested in starting partial wave analysis.


title: Extracting black-hole rotational energy: The generalized Penrose process

abstract: In the case involving particles the necessary and sufficient condition for the Penrose process to extract energy from a rotating black hole is absorption of particles with negative energies and angular momenta. No torque at the black-hole horizon occurs. In this article we consider the case of arbitrary fields or matter described by an unspecified, general energy-momentum tensor $T_{\mu \nu}$ and show that the necessary and sufficient condition for extraction of a black hole's rotational energy is analogous to that in the mechanical Penrose process: absorption of negative energy and negative angular momentum. We also show that a necessary condition for the Penrose process to occur is for the Noether current (the conserved energy-momentum density vector) to be spacelike or past directed (timelike or null) on some part of the horizon. In the particle case, our general criterion for the occurrence of a Penrose process reproduces the standard result. In the case of relativistic jet-producing "magnetically arrested disks" we show that the negative energy and angular-momentum absorption condition is obeyed when the Blandford-Znajek mechanism is at work, and hence the high energy extraction efficiency up to $\sim 300\%$ found in recent numerical simulations of such accretion flows results from tapping the black hole's rotational energy through the Penrose process. We show how black-hole rotational energy extraction works in this case by describing the Penrose process in terms of the Noether current.


title: On-demand dark soliton train manipulation in a spinor polariton condensate

abstract: We theoretically demonstrate the generation of dark solitons trains in a one-dimensional exciton- polariton condensate within an experimentally trivial scheme. In particular we show that the fre- quency of the train can be finely tuned either electrically of fully optically to provide a stable and efficient output signal modulation. Taking the polarization degree of freedom into account we elucidate the possibility to additionally form on-demand half-solitons trains.


title: The Inferred Cardiogenic Gene Regulatory Network in the Mammalian Heart

abstract: Cardiac development is complex, multiscale process encompassing cell fate adoption, differentiation and morphogenesis. To elucidate pathways underlying this process, a recently developed algorithm to reverse engineer gene regulatory networks was applied to time-course microarray data obtained from the developing mouse heart. The algorithm generates many different putative network topologies that are capable of explaining the experimental data via model simulation. To cull specious network interactions, thousands of topologies are merged and filtered to generate a scale-free, hierarchical network. The network is validated with known gene interactions and used to identify regulatory pathways critical to the developing mammalian heart. The predicted gene interactions are prioritized using semantic similarity and gene profile uniqueness metrics. Using these metrics, the network is expanded to include all known mouse genes to form the most likely cardiogenic gene regulatory network. The method outlined herein provides an informative approach to network inference and leads to clear testable hypotheses related to gene regulation.


title: Positions of Pluto extracted from digitized Pulkovo photographic plates taken in 1930 - 1960

abstract: We present the results of determination of Pluto's positions derived from photographic plates taken in 1930 - 1960. Observations were made with Normal Astrograph at Pulkovo Observatory. Digitization of these plates was performed with high precision scanner at Royal Observatory of Belgium (ROB Digitizer). Mean values of standard errors of plate positions (x,y) lie between 12 and 18 mas. The UCAC4 catalogue was used as an astrometric calibrator. Standard errors of equatorial coordinates obtained are within 85 to 100 mas. Final table contains 63 positions of Pluto referred to the HCRF/UCAC4 frame.


title: Modelling binary rotating stars by new population synthesis code BONNFIRES

abstract: BONNFIRES, a new generation of population synthesis code, can calculate nuclear reaction, various mixing processes and binary interaction in a timely fashion. We use this new population synthesis code to study the interplay between binary mass transfer and rotation. We aim to compare theoretical models with observations, in particular the surface nitrogen abundance and rotational velocity. Preliminary results show binary interactions may explain the formation of nitrogen-rich slow rotators and nitrogen-poor fast rotators, but more work needs to be done to estimate whether the observed frequencies of those stars can be matched.


title: Analysis of interior penalty discontinuous Galerkin methods for the Allen-Cahn equation and the mean curvature flow

abstract: This paper develops and analyzes two fully discrete interior penalty discontinuous Galerkin (IP-DG) methods for the Allen-Cahn equation, which is a nonlinear singular perturbation of the heat equation and originally arises from phase transition of binary alloys in materials science, and its sharp interface limit (the mean curvature flow) as the perturbation parameter tends to zero. Both fully implicit and energy-splitting time-stepping schemes are proposed. The primary goal of the paper is to derive sharp error bounds which depend on the reciprocal of the perturbation parameter $\epsilon$ (also called "interaction length") only in some lower polynomial order, instead of exponential order, for the proposed IP-DG methods. The derivation is based on a refinement of the nonstandard error analysis technique first introduced in [12]. The centerpiece of this new technique is to establish a spectrum estimate result in totally discontinuous DG finite element spaces with a help of a similar spectrum estimate result in the conforming finite element spaces which was established in [12]. As a nontrivial application of the sharp error estimates, they are used to establish convergence and the rates of convergence of the zero level sets of the fully discrete IP-DG solutions to the classical and generalized mean curvature flow. Numerical experiment results are also presented to gauge the theoretical results and the performance of the proposed fully discrete IP-DG methods.


title: Quantifying age- and gender-related diabetes comorbidity risks using nation-wide big claims data

abstract: Currently emerging "big data" techniques are reshaping medical science into a data science. Medical claims data allow assessing an entire nation's health state in a quantitative way, in particular with regard to the occurrences and consequences of chronic and pandemic diseases like diabetes. We develop a quantitative, statistical approach to test for associations between the incidence of type 1 or type 2 diabetes and any possible other disease as provided by the ICD10 diagnosis codes using a complete set of Austrian inpatient data. With a new co-occurrence analysis the relative risks for each possible comorbidity are studied as a function of patient age and gender, a temporal analysis investigates whether the onset of diabetes typically precedes or follows the onset of the other disease. The samples is always of maximal size, i.e. contains all patients with that comorbidity within the country. The present study is an equivalent of almost 40,000 studies, all with maximum patient number available. Out of more than thousand possible associations, 123 comorbid diseases for type 1 or type 2 diabetes are identified at high significance levels. Well known diabetic comorbidities are recovered, such as retinopathies, hypertension, chronic kidney diseases, etc. This validates the method. Additionally, a number of comorbidities are identified which have only been recognized to a lesser extent, for example epilepsy, sepsis, or mental disorders. The temporal evolution, age, and gender-dependence of these comorbidities are discussed. The new statistical-network methodology developed here can be readily applied to other chronic diseases.


title: Fluctuation of similarity (FLUS) to detect transitions between distinct dynamical regimes in short time series

abstract: Recently a method which employs computing of fluctuations in a measure of nonlinear similarity based on local recurrence properties in a univariate time series, was introduced to identify distinct dynamical regimes and transitions between them in a short time series [1]. Here we present the details of the analytical relationships between the newly introduced measure and the well known concepts of attractor dimensions and Lyapunov exponents. We show that the new measure has linear dependence on the effective dimension of the attractor and it measures the variations in the sum of the Lyapunov spectrum. To illustrate the practical usefulness of the method, we employ it to identify various types of dynamical transitions in different nonlinear models. Also, we present testbed examples for the new method's robustness against the presence of noise and missing values in the time series. Furthermore, we use this method to analyze time series from the field of social dynamics, where we present an analysis of the US crime record's time series from the year 1975 to 1993. Using this method, we have found that dynamical complexity in robberies was influenced by the unemployment rate till late 1980's. We have also observed a dynamical transition in homicide and robbery rates in the late 1980's and early 1990's, leading to increase in the dynamical complexity of these rates.


title: Radiative pi^pm gamma transitions of excited light-quark mesons in the covariant oscillator quark model

abstract: The COMPASS collaboration, as a part of their hadron spectroscopy program, is going to measure the radiative decay widths of light-quark mesons via the Primakoff reactions. In this letter we study the photon couplings of light-quark q qbar states in the covariant oscillator quark model and evaluate the transition rates for {rho(770), b_1(1235), a_1(1260), a_2(1320), pi_2(1670), rho_3(1690), rho(1700)}^pm to pi^pm gamma, which are expected to be measured by the COMPASS. Such photon couplings could be useful not only for understanding the internal structures of observed light-quark mesons and their quark-model classification but also for the ongoing experimental studies by COMPASS.


title: Giant Collective Spin-Orbit Field in a Quantum Well: Fine Structure of Spin Plasmons

abstract: We employ inelastic light scattering with magnetic fields to study intersubband spin plasmons in a quantum well. We demonstrate the existence of a giant collective spin-orbit (SO) field that splits the spin-plasmon spectrum into a triplet. The effect is remarkable as each individual electron would be expected to precess in its own momentum-dependent SO field, leading to D'yakonov-Perel' dephasing. Instead, many-body effects lead to a striking organization of the SO fields at the collective level. The macroscopic spin moment is quantized by a uniform collective SO field, five times higher than the individual SO field. We provide a momentum-space cartography of this field.


title: Notes on Enhancement of Flavor Symmetry and 5d Superconformal Index

abstract: The UV fixed point theory of SU(2) gauge theory with N_f = 0,1,...,7 flavors is believed to have the enlarged E_{N_f +1} flavor symmetry. Actually it is not easy to check this conjecture because the UV theory is strongly-coupled, however, computation of certain SUSY protected quantities provides strong evidence for the enhancement of flavor symmetry. We study the superconformal index for SU(2) gauge theory with N_f = 0, 1 flavors in details, and we give a support for the enhancement by studying combinatorial structure of the superconformal indexes of these theories. We also give a nontrivial evidence that the local F_2 geometry leads to the E_1 superconformal field theory.


title: The Impact of Organic Friction Modifiers on Engine Oil Tribofilms

abstract: Organic friction modifiers (OFMs) are important additives in the lubrication of machines and especially of car engines where performance improvements are constantly sought-after. Together with zinc dialkyldithiophosphates (ZDDPs) antiwear additives, OFMs have a predominant impact on the tribological behaviour of the lubricant. In the current study, the influence of OFMs on the generation, tribological properties and chemistry of ZDDP tribofilms has been investigated by combining tribological experiments (MTM) with in-situ film thickness measurements through optical interference imaging (SLIM), Alicona profilometry and X-ray photoelectron spectroscopy. OFMs and antiwear additives have been found to competitively react/adsorb on the rubbing ferrous substrates in a tribological contact. The formation and removal (through wear) of tribofilms are dynamic processes which result from the simultaneous interaction of these two additives with the surface of the wear track. By carefully selecting the chemistry of OFMs, the formulator can achieve lubricants that generate ZDDP antiwear films of optimum thickness, morphology and friction according to the application-specific requirements.


title: A Criterion for the Bipartite Separability of Bell Diagonal States

abstract: A decomposition form is introduced in this report to establish a criterion for the bi-partite separability of Bell diagonal states. A such criterion takes a quadratic form of the coefficients of a given Bell diagonal states and can be derived via a simple algorithmic calculation of its invariants. In addition, the criterion can be extended to a quantum system of higher dimension.


title: A Decomposition Form of the Werner State

abstract: In this report, a scheme different from the PT and Wootters concurrence is developed to acquire a criterion to investigate the bipartite separability of the Werner state.


title: Existence of a persistent hub in the convex preferential attachment model

abstract: A vertex of a randomly growing graph is called a persistent hub if at all but finitely many moments of time it has the maximal degree in the graph. We establish the existence of a persistent hub in the Barab\'asi--Albert random graph model with probability one. We also extend this result to the class of convex preferential attachment graphs, where a vertex of degree $k$ gets a new edge with probability proportional to some convex function of $k$.


title: Algebraic Quantum Error-Correction Codes

abstract: Based on the group structure of a unitary Lie algebra, a scheme is provided to systematically and exhaustively generate quantum error correction codes, including the additive and nonadditive codes. The syndromes in the process of error-correction distinguished by different orthogonal vector subspaces, the coset subspaces. Moreover, the generated codes can be classified into four types with respect to the spinors in the unitary Lie algebra and a chosen initial quantum state.


title: Precise determination of the spectroscopic g-factor using broadband ferromagnetic resonance spectroscopy

abstract: We demonstrate that the spectroscopic g-factor can be determined with high precision and accuracy by broadband ferromagnetic resonance measurements and applying an asymptotic analysis to the data. Spectroscopic data used to determine the g-factor is always obtained over a finite range of frequencies, which can result in significant errors in the fitted values of the spectroscopic g-factor. We show that by applying an asymptotic analysis to broadband datasets, precise values of the intrinsic g-factor can be determined with errors well below 1 %, even when the exact form of the Kittel equation (which describes the relationship between the frequency and resonance field) is unknown. We demonstrate this methodology with measured data obtained for sputtered Ni80Fe20 ("Permalloy") thin films of varied thicknesses, where we determine the bulk g-factor value to be 2.109 +/- 0.003. Such an approach is further validated by application to simulated data that includes both noise and an anisotropy that is not included in the Kittel equation that was used in the analysis. Finally, we show a correlation of thickness and interface structure to the magnitude of the asymptotic behavior, which provide insight into additional mechanisms that may induce deviations from the Kittel equation.


title: Coupling between bulk- and surface chemistry in suspensions of charged colloids

abstract: The ionic composition and pair correlations in fluid phases of realistically salt-free charged colloidal sphere suspensions are calculated in the primitive model. We obtain the number densities of all ionic species in suspension, including low-molecular weight microions, and colloidal macroions with acidic surface groups, from a self-consistent solution of a coupled physicochemical set of nonlinear algebraic equations and non-mean-field liquid integral equations. Here, we study suspensions of colloidal spheres with sulfonate or silanol surface groups, suspended in demineralized water that is saturated with carbon dioxide under standard atmosphere. The only input required for our theoretical scheme are the acidic dissociation constants pKa, and effective sphere diameters of all involved ions. Our method allows for an ab initio calculation of colloidal bare and effective charges, at high numerical efficiency.


title: The Impact of Magnification and Size Bias on Weak Lensing Power Spectrum and Peak Statistics

abstract: The weak lensing power spectrum is a powerful tool to probe cosmological parameters. Additionally, lensing peak counts contain cosmological information beyond the power spectrum. Both of these statistics can be affected by the preferential selection of source galaxies in patches of the sky with high magnification, as well as by the dilution in the source galaxy surface density in such regions. If not accounted for, these biases introduce systematic errors for cosmological measurements. Here we quantify these systematic errors, using convergence maps from a suite of ray-tracing N-body simulations. At the cut-off magnitude m of on-going and planned major weak lensing surveys, the logarithmic slope of the cumulative number counts s = dlog[n(>m)]/dlog(m) is in the range 0.1 < s < 0.5. At s = 0.2, expected in the I band for LSST, the inferred values of Omega_m, w and sigma_8 are biased by many sigma (where sigma denotes the marginalized error) and therefore the biases will need to be carefully modeled. We also find that the parameters are biased differently in the (Omega_m, w, sigma_8) parameter space when the power spectrum and when the peak counts are used. In particular, w derived from the power spectrum is less affected than w derived from peak counts, while the opposite is true for the best-constrained combination of [sigma_8 Omega_m^gamma] (with gamma=0.62 from the power spectrum and gamma = 0.48 from peak counts). This suggests that the combination of the power spectrum and peak counts can help mitigate the impact of magnification and size biases.


title: Obtaining localization properties efficiently using the Kubo-Greenwood formalism

abstract: We establish, through numerical calculations and comparisons with a recursive Green's function based implementation of the Landauer-B\"uttiker formalism, an efficient method for studying Anderson localization in quasi-one-dimensional and two-dimensional systems using the Kubo-Greenwood formalism. Although the recursive Green's function method can be used to obtain the localization length of a mesoscopic conductor, it is numerically very expensive for systems that contain a large number of atoms transverse to the transport direction. On the other hand, linear-scaling has been achieved with the Kubo-Greenwood method, enabling the study of effectively two-dimensional systems. While the propagating length of the charge carriers will eventually saturate to a finite value in the localized regime, the conductances given by the Kubo-Greenwood method and the recursive Green's function method agree before the saturation. The converged value of the propagating length is found to be directly proportional to the localization length obtained from the exponential decay of the conductance.


title: Nuclear three-body problem in the complex energy plane: Complex-Scaling-Slater method

abstract: The physics of open quantum systems is an interdisciplinary area of research. The nuclear "openness" manifests itself through the presence of the many-body continuum representing various decay, scattering, and reaction channels. As the radioactive nuclear beam experimentation extends the known nuclear landscape towards the particle drip lines, the coupling to the continuum space becomes exceedingly more important. Of particular interest are weakly bound and unbound nuclear states appearing around particle thresholds. Theories of such nuclei must take into account their open quantum nature. To describe open quantum systems, we introduce a Complex Scaling (CS) approach in the Slater basis. We benchmark it with the complex-energy Gamow Shell Model (GSM) by studying energies and wave functions of the bound and unbound states of the two-neutron halo nucleus 6He viewed as an $\alpha$+ n + n cluster system. In the CS approach, we use the Slater basis, which exhibits the correct asymptotic behavior at large distances. To extract particle densities from the back-rotated CS solutions, we apply the Tikhonov regularization procedure, which minimizes the ultraviolet numerical noise. While standard applications of the inverse complex transformation to the complex-rotated solution provide unstable results, the stabilization method fully reproduces the GSM benchmark. We also propose a method to determine the smoothing parameter of the Tikhonov regularization. The combined suite of CS-Slater and GSM techniques has many attractive features when applied to nuclear problems involving weakly-bound and unbound states. While both methods can describe energies, total widths, and wave functions of nuclear states, the CS-Slater method, if it can be applied, can provide an additional information about partial energy widths associated with individual thresholds.


title: Principle of Maximum Entanglement Entropy and Local Physics of Correlated many-body Electron-Systems

abstract: We argue that, because of the quantum-entanglement, the local physics of the strongly-correlated materials at zero temperature is described in very good approximation by a simple generalized Gibbs distribution, which depends on a relatively small number local quantum thermodynamical potentials. We demonstrate that our statement is exact in certain limits, and we perform numerical calculations of the iron compounds FeSe and FeTe and of the elemental cerium by employing the Gutzwiller Approximation (GA) that strongly support our theory in general.


title: Electron Pulse Compression with a Practical Reflectron Design for Ultrafast Electron Diffraction

abstract: Ultrafast electron diffraction (UED) is a powerful method for studying time-resolved structural changes. Currently, space charge induced temporal broadening prevents obtaining high brightness electron pulses with sub-100 fs durations limiting the range of phenomena that can be studied with this technique. We review the state of the the art of UED in this respect and propose a practical design for reflectron based pulse compression which utilizes only electro-static optics and has a tunable temporal focal point. Our simulation shows that this scheme is capable of compressing an electron pulse containing 100,000 electrons with 60:1 temporal compression ratio.


title: Protrusion fluctuations direct cell motion

abstract: Many physiological phenomena involve directional cell migration. It is usually attributed to chemical gradients in vivo. Recently, other cues have been shown to guide cells in vitro, including stiffness/adhesion gradients or micro-patterned adhesive motifs. However, the cellular mechanism leading to these biased migrations remains unknown, and, often, even the direction of motion is unpredictable. In this study, we show the key role of fluctuating protrusions on ratchet-like structures in driving NIH3T3 cell migration. We identified the concept of efficient protrusion and an associated direction index. Our analysis of the protrusion statistics facilitated the quantitative prediction of cell trajectories in all investigated conditions. We varied the external cues by changing the adhesive patterns. We also modified the internal cues using drug treatments, which modified the protrusion activity. Stochasticity affects the short- and long-term steps. We developed a theoretical model showing that an asymmetry in the protrusion fluctuations is sufficient for predicting all measures associated with the long-term motion, which can be described as a biased persistent random walk.


title: Quantum Shuffles and Quantum Supergroups of Basic Type

abstract: We initiate the study of several distinguished bases for the positive half of a quantum supergroup $U_q$ associated to a general super Cartan datum $(\mathrm{I}, (\cdot,\cdot))$ of basic type inside a quantum shuffle superalgebra. The combinatorics of words for an arbitrary total ordering on $\mathrm{I}$ is developed in connection with the root system associated to $\mathrm{I}$. The monomial, Lyndon, and PBW bases of $U_q$ are constructed, and moreover, a direct proof of the orthogonality of the PBW basis is provided within the framework of quantum shuffles. Consequently, the canonical basis is constructed for $U_q$ associated to the standard super Cartan datum of type $\mathfrak{gl}(n|1)$, $\mathfrak{osp}(1|2n)$, or $\mathfrak{osp}(2|2n)$ or an arbitrary non-super Cartan datum. In the non-super case, this refines Leclerc's work and provides a new self-contained construction of canonical bases. The canonical bases of $U_q$, of its polynomial modules, as well as of Kac modules in the case of quantum $\mathfrak{gl}(2|1)$ are explicitly worked out.


title: Phase noise of oscillators with unsaturated amplifiers

abstract: We study the role of amplifier saturation in eliminating feedback noise in self-sustained oscillators. We extend previous works that use a saturated amplifier to quench fluctuations in the feedback magnitude, while simultaneously tuning the oscillator to an operational point at which the resonator nonlinearity cancels fluctuations in the feedback phase. We consider a generalized model which features an amplitude-dependent amplifier gain function. This allows us to determine the total oscillator phase noise in realistic configurations due to noise in both quadratures of the feedback, and to show that it is not necessary to drive the resonator to large oscillation amplitudes in order to eliminate noise in the phase of the feedback.


title: Coding theorems for compound problems via quantum R\'enyi divergences

abstract: Recently, a new notion of quantum R\'enyi divergences has been introduced by M\"uller-Lennert, Dupuis, Szehr, Fehr and Tomamichel, J.Math.Phys. 54:122203, (2013), and Wilde, Winter, Yang, Commun.Math.Phys. 331:593--622, (2014), that has found a number of applications in strong converse theorems. Here we show that these new R\'enyi divergences are also useful tools to obtain coding theorems in the direct domain of various problems. We demonstrate this by giving new and considerably simplified proofs for the achievability parts of Stein's lemma with composite null hypothesis, universal state compression, and the classical capacity of compound classical-quantum channels, based on single-shot error bounds already available in the literature, and simple properties of the quantum R\'enyi divergences. The novelty of our proofs is that the composite/compound coding theorems can be almost directly obtained from the single-shot error bounds, with essentially the same effort as for the case of simple null-hypothesis/single source/single channel.


title: KCH representations, augmentations, and $A$-polynomials

abstract: We describe a correspondence between augmentations and certain representations of the knot group. The correspondence makes the 2-variable augmentation polynomial into a generalization of the classical $A$-polynomial. It also associates to an augmentation a rank, which is bounded by the bridge number and shares its behavior under connect sums. We also study augmentations with rank equal to the braid index.


title: Knowledge-fused differential dependency network models for detecting significant rewiring in biological networks

abstract: Modeling biological networks serves as both a major goal and an effective tool of systems biology in studying mechanisms that orchestrate the activities of gene products in cells. Biological networks are context specific and dynamic in nature. To systematically characterize the selectively activated regulatory components and mechanisms, the modeling tools must be able to effectively distinguish significant rewiring from random background fluctuations. We formulated the inference of differential dependency networks that incorporates both conditional data and prior knowledge as a convex optimization problem, and developed an efficient learning algorithm to jointly infer the conserved biological network and the significant rewiring across different conditions. We used a novel sampling scheme to estimate the expected error rate due to random knowledge and based on which, developed a strategy that fully exploits the benefit of this data-knowledge integrated approach. We demonstrated and validated the principle and performance of our method using synthetic datasets. We then applied our method to yeast cell line and breast cancer microarray data and obtained biologically plausible results.


title: A diabatic representation of the two lowest electronic states of Li3

abstract: Using the Multi-Reference Configuration Interaction method, the adiabatic potential energy surfaces of Li3 are computed. The two lowest electronic states are bound and exhibit a conical intersection. By fitting the calculated potential energy surfaces to the cubic Exe Jahn-Teller model we extract the effective Jahn-Teller parameters corresponding to Li3. These are used to set up the transformation matrix which transforms from the adiabatic to a diabatic representation. This diabatization method gives a Hamiltonian for Li3 which is free from singular adiabatic couplings and should be accurate for large internuclear distances, and it thereby allows for bound dynamics in the vicinity of the conical intersection to be explored.


title: Successive Nonnegative Projection Algorithm for Robust Nonnegative Blind Source Separation

abstract: In this paper, we propose a new fast and robust recursive algorithm for near-separable nonnegative matrix factorization, a particular nonnegative blind source separation problem. This algorithm, which we refer to as the successive nonnegative projection algorithm (SNPA), is closely related to the popular successive projection algorithm (SPA), but takes advantage of the nonnegativity constraint in the decomposition. We prove that SNPA is more robust than SPA and can be applied to a broader class of nonnegative matrices. This is illustrated on some synthetic data sets, and on a real-world hyperspectral image.


title: Exact model structures and recollements

abstract: We show how to obtain recollements of triangulated categories using the theory of exact model structures. After noting how the theory relates to well-known notions in the simplest case of Frobenius categories, we apply these ideas to recollements of homotopy categories of chain complexes. In short, we give model categorical explanations for the classical Verdier localization recollement as well as several recollement situations due to Neeman and Murfet.


title: Derivatives of the tree function

abstract: We study some sequences of polynomials that appear when we consider the successive derivatives of the tree function (or Lambert's W function). We show in particular that they are related with a generalization of Cayley trees, called Greg trees. Besides the combinatorial result in itself, it is interesting to see how this is related with previous work: similar problems were considered first by Ramanujan, and more recently in the theory of completely monotonic functions and its link with probability. Also of great interest is the fact that these Greg trees were introduced in a problem of textual criticism, as a kind a genealogical trees, where they had a priori no mathematical meaning.


title: Matchmaker, Matchmaker, Make Me a Match: Migration of Populations via Marriages in the Past

abstract: The study of human mobility is both of fundamental importance and of great potential value. For example, it can be leveraged to facilitate efficient city planning and improve prevention strategies when faced with epidemics. The newfound wealth of rich sources of data---including banknote flows, mobile phone records, and transportation data---has led to an explosion of attempts to characterize modern human mobility. Unfortunately, the dearth of comparable historical data makes it much more difficult to study human mobility patterns from the past. In this paper, we present an analysis of long-term human migration, which is important for processes such as urbanization and the spread of ideas. We demonstrate that the data record from Korean family books (called "jokbo") can be used to estimate migration patterns via marriages from the past 750 years. We apply two generative models of long-term human mobility to quantify the relevance of geographical information to human marriage records in the data, and we find that the wide variety in the geographical distributions of the clans poses interesting challenges for the direct application of these models. Using the different geographical distributions of clans, we quantify the "ergodicity" of clans in terms of how widely and uniformly they have spread across Korea, and we compare these results to those obtained using surname data from the Czech Republic. To examine population flow in more detail, we also construct and examine a population-flow network between regions. Based on the correlation between ergodicity and migration in Korea, we identify two different types of migration patterns: diffusive and convective. We expect the analysis of diffusive versus convective effects in population flows to be widely applicable to the study of mobility and migration patterns across different cultures.

categories: physics.soc-ph cond-mat.dis-nn cs.CE nlin.AO q-bio.PE


title: Digraphs and cycle polynomials for free-by-cyclic groups

abstract: Let $\phi \in \mbox{Out}(F_n)$ be a free group outer automorphism that can be represented by an expanding, irreducible train-track map. The automorphism $\phi$ determines a free-by-cyclic group $\Gamma=F_n \rtimes_\phi \mathbb Z,$ and a homomorphism $\alpha \in H^1(\Gamma; \mathbb Z)$. By work of Neumann, Bieri-Neumann-Strebel and Dowdall-Kapovich-Leininger, $\alpha$ has an open cone neighborhood $\mathcal A$ in $H^1(\Gamma;\mathbb R)$ whose integral points correspond to other fibrations of $\Gamma$ whose associated outer automorphisms are themselves representable by expanding irreducible train-track maps. In this paper, we define an analog of McMullen's Teichm\"uller polynomial that computes the dilatations of all outer automorphism in $\mathcal A$.


title: Order of Magnitude Smaller Limit on the Electric Dipole Moment of the Electron

abstract: The Standard Model (SM) of particle physics fails to explain dark matter and why matter survived annihilation with antimatter following the Big Bang. Extensions to the SM, such as weak-scale Supersymmetry, may explain one or both of these phenomena by positing the existence of new particles and interactions that are asymmetric under time-reversal (T). These theories nearly always predict a small, yet potentially measurable ($10^{-27}$-$10^{-30}$ $e$ cm) electron electric dipole moment (EDM, $d_e$), which is an asymmetric charge distribution along the spin ($\vec{S}$). The EDM is also asymmetric under T. Using the polar molecule thorium monoxide (ThO), we measure $d_e = (-2.1 \pm 3.7_\mathrm{stat} \pm 2.5_\mathrm{syst})\times 10^{-29}$ $e$ cm. This corresponds to an upper limit of $|d_e| < 8.7\times 10^{-29}$ $e$ cm with 90 percent confidence, an order of magnitude improvement in sensitivity compared to the previous best limits. Our result constrains T-violating physics at the TeV energy scale.


title: Improving Photoacoustic-guided Focusing in Scattering Media by Spectrally Filtered Detection

abstract: We experimentally and numerically study the potential of photoacoustic-guiding for light focusing through scattering samples via wavefront-shaping and iterative optimization. We experimentally demonstrate that the focusing efficiency on an extended absorber can be improved by iterative optimization of the high frequency components of the broadband photoacoustic signal detected with a spherically focused transducer. We demonstrate more than 8-fold increase in the photoacoustic signal generated by a 30 microns wire using a narrow frequency band around 60MHz. We numerically confirm that such optimization leads to a smaller optical focus than using the low frequency content of the photoacoustic feedback.