title: The Effective Planck Mass and the Scale of Inflation

abstract: Observable quantities in cosmology are dimensionless, and therefore independent of the units in which they are measured. This is true of all physical quantities associated with the primordial perturbations that source cosmic microwave background anisotropies such as their amplitude and spectral properties. However, if one were to try and infer an absolute energy scale for inflation-- a priori, one of the more immediate corollaries of detecting primordial tensor modes-- one necessarily makes reference to a particular choice of units, the natural choice for which is Planck units. In this note, we discuss various aspects of how inferring the energy scale of inflation is complicated by the fact that the effective strength of gravity as seen by inflationary quanta necessarily differs from that seen by gravitational experiments at presently accessible scales. The uncertainty in the former relative to the latter has to do with the unknown spectrum of universally coupled particles between laboratory scales and the putative scale of inflation. These intermediate particles could be in hidden as well as visible sectors or could also be associated with Kaluza-Klein resonances associated with a compactification scale below the scale of inflation. We discuss various implications for cosmological observables.


title: A decoupled unconditionally stable numerical scheme for the Cahn-Hilliard-Hele-Shaw system

abstract: We propose a novel decoupled unconditionally stable numerical scheme for the simulation of two-phase flow in a Hele-Shaw cell which is governed by the Cahn-Hilliard-Hele-Shaw system (CHHS) with variable viscosity. The temporal discretization of the Cahn-Hilliard equation is based on a convex-splitting of the associated energy functional. Moreover, the capillary forcing term in the Darcy equation is separated from the pressure gradient at the time discrete level by using an operator-splitting strategy. Thus the computation of the nonlinear Cahn-Hilliard equation is completely decoupled from the update of pressure. Finally, a pressure-stabilization technique is used in the update of pressure so that at each time step one only needs to solve a Poisson equation with constant coefficient. We show that the scheme is unconditionally stable. Numerical results are presented to demonstrate the accuracy and efficiency of our scheme.


title: Autocorrelation analysis for the unbiased determination of power-law exponents in single-quantum-dot blinking

abstract: We present an unbiased and robust analysis method for power-law blinking statistics in the photoluminescence of single nano-emitters, allowing us to extract both the bright- and dark-state power-law exponents from the emitters' intensity autocorrelation functions. As opposed to the widely-used threshold method, our technique therefore does not require discriminating the emission levels of bright and dark states in the experimental intensity timetraces. We rely on the simultaneous recording of 450 emission timetraces of single CdSe/CdS core/shell quantum dots at a frame rate of 250 Hz with single photon sensitivity. Under these conditions, our approach can determine ON and OFF power-law exponents with a precision of 3% from a comparison to numerical simulations, even for shot-noise-dominated emission signals with an average intensity below 1 photon per frame and per quantum dot. These capabilities pave the way for the unbiased, threshold-free determination of blinking power-law exponents at the micro-second timescale.


title: Characterization of 2D rational local conformal nets and its boundary conditions: the maximal case

abstract: Let $\mathcal{A}$ be a completely rational local M\"obius covariant net on $S^1$, which describes a set of chiral observables. We show that local M\"obius covariant nets $\mathcal{B}_2$ on 2D Minkowski space which contains $\mathcal{A}$ as chiral left-right symmetry are in one-to-one correspondence with Morita equivalence classes of Q-systems in the unitary modular tensor category $\mathrm{DHR}(\mathcal{A})$. The M\"obius covariant boundary conditions with symmetry $\mathcal{A}$ of such a net $\mathcal{B}_2$ are given by the Q-systems in the Morita equivalence class or by simple objects in the module category modulo automorphisms of the dual category. We generalize to reducible boundary conditions. To establish this result we define the notion of Morita equivalence for Q-systems (special symmetric $\ast$-Frobenius algebra objects) and non-degenerately braided subfactors. We prove a conjecture by Kong and Runkel, namely that Rehren's construction (generalized Longo-Rehren construction, $\alpha$-induction construction) coincides with the categorical full center. This gives a new view and new results for the study of braided subfactors.


title: Parton distributions for the LHC Run II

abstract: We present NNPDF3.0, the first set of parton distribution functions (PDFs) determined with a methodology validated by a closure test. NNPDF3.0 uses a global dataset including HERA-II deep-inelastic inclusive cross-sections, the combined HERA charm data, jet production from ATLAS and CMS, vector boson rapidity and transverse momentum distributions from ATLAS, CMS and LHCb, W+c data from CMS and top quark pair production total cross sections from ATLAS and CMS. Results are based on LO, NLO and NNLO QCD theory and also include electroweak corrections. To validate our methodology, we show that PDFs determined from pseudo-data generated from a known underlying law correctly reproduce the statistical distributions expected on the basis of the assumed experimental uncertainties. This closure test ensures that our methodological uncertainties are negligible in comparison to the generic theoretical and experimental uncertainties of PDF determination. This enables us to determine with confidence PDFs at different perturbative orders and using a variety of experimental datasets ranging from HERA-only up to a global set including the latest LHC results, all using precisely the same validated methodology. We explore some of the phenomenological implications of our results for the upcoming 13 TeV Run of the LHC, in particular for Higgs production cross-sections.


title: Local character of the highest antiferromagnetic Ce-system CeTi{1-x}Sc{x} Ge

abstract: The highest antiferromagnetic (AFM) temperature in Ce based compounds has been reported for CeScGe with Tn=47K, but its local or itinerant nature was not deeply investigated yet. In order to shed more light into this unusually high ordering temperature we have investigated structural, magnetic, transport and thermal properties of CeTi{1-x}Sc{x}Ge alloys within the range of stability of the CeScSi-type structure: 0.25<x<1. Along this concentration range, this strongly anisotropic system presents a complex magnetic phase diagram with a continuous modification of its magnetic behavior, from ferromagnetism (FM)for 0.25<x<0.50 (with 7K<Tc<16K) to AFM for 0.60<x<1 (with 19K<Tn<47K). The onset of the AFM phase is associated to a metamagnetic transition with a critical field increasing from H{cr}=0 at x~0.55 to ~6Tesla at x=1, coincident with an increasing contribution of the first excited crystal electric field doublet. At a critical point x{cr}=0.65 a second transition appears at Tl<Tn. In contrast to observations in itinerant systems like CeRh{2}Si{2} or CeRh{3}B{2}, no evidences for significant hybridization of the 4f electrons at large Sc contents were found. Therefore, the exceptionally large Tn of CeScGe can be attributed to an increasing RKKY inter-layer interaction as Sc content grows.


title: $C^{\infty}$ Stability, Canonical Maps, and Discrete Dynamics

abstract: We study a discrete dynamical system designed to find a 'most holomorphic' connection on a smooth complex vector bundle $E$. We examine the relation between the distance of the chern classes of $E$ from the $(p,p)$ axis of the Hodge diamond and singularity formation. Canonical connections and canonical metrics pulled back from Grassmannians play a major role, and we review their differential geometry. As an exercise in the geometry of canonical connections, we include an expression for the curvature in terms of heat kernels.


title: Decidable Models of Recursive Asynchronous Concurrency

abstract: Asynchronously communicating pushdown systems (ACPS) that satisfy the empty-stack constraint (a pushdown process may receive only when its stack is empty) are a popular decidable model for recursive programs with asynchronous atomic procedure calls. We study a relaxation of the empty-stack constraint for ACPS that permits concurrency and communication actions at any stack height, called the shaped stack constraint, thus enabling a larger class of concurrent programs to be modelled. We establish a close connection between ACPS with shaped stacks and a novel extension of Petri nets: Nets with Nested Coloured Tokens (NNCTs). Tokens in NNCTs are of two types: simple and complex. Complex tokens carry an arbitrary number of coloured tokens. The rules of NNCT can synchronise complex and simple tokens, inject coloured tokens into a complex token, and eject all tokens of a specified set of colours to predefined places. We show that the coverability problem for NNCTs is Tower-complete. To our knowledge, NNCT is the first extension of Petri nets, in the class of nets with an infinite set of token types, that has primitive recursive coverability. This result implies Tower-completeness of coverability for ACPS with shaped stacks.


title: Demonstrating Continuous Variable EPR Steering in spite of Finite Experimental Capabilities using Fano Steering Bounds

abstract: We show how one can demonstrate continuous-variable Einstein-Podolsky-Rosen (EPR) steering without needing to characterize entire measurement probability distributions. To do this, we develop a modified Fano inequality useful for discrete measurements of continuous variables, and use it to bound the conditional uncertainties in continuous-variable entropic EPR-steering inequalities. With these bounds, we show how one can hedge against experimental limitations including a finite detector size, dead space between pixels, and any such factors that impose an incomplete sampling of the true measurement probability distribution. Furthermore, we use experimental data from the position and momentum statistics of entangled photon pairs in parametric downconversion to show that this method is sufficiently sensitive for practical use.


title: Poisson structures on twistor spaces of hyperkaehler and HKT manifolds

abstract: We characterize HKT structure in terms of nondegenrate complex Poisson bivector on hypercomplex manifold. We extend the characterization to the twistor space. After considering the flat case in detail, we show that the twistor space of hyperkaehler manifold admits a holomorphic Poisson structure. We briefly mention the relation to quaternionic and hypercomplex deformations on tori and K3 surfaces


title: A Consistency Result for Bayes Classifiers with Censored Response Data

abstract: Naive Bayes classifiers have proven to be useful in many prediction problems with complete training data. Here we consider the situation where a naive Bayes classifier is trained with data where the response is right censored. Such prediction problems are for instance encountered in profiling systems used at National Employment Agencies. In this paper we propose the maximum collective conditional likelihood estimator for the prediction and show that it is strongly consistent under the usual identifiability condition.


title: Earth--Mars Transfers with Ballistic Capture

abstract: We construct a new type of transfer from the Earth to Mars, which ends in ballistic capture. This results in a substantial savings in capture $\Delta v$ from that of a classical Hohmann transfer under certain conditions. This is accomplished by first becoming captured at Mars, very distant from the planet, and then from there, following a ballistic capture transfer to a desired altitude within a ballistic capture set. This is achieved by manipulating the stable sets, or sets of initial conditions whose orbits satisfy a simple definition of stability. This transfer type may be of interest for Mars missions because of lower capture $\Delta v$, moderate flight time, and flexibility of launch period from the Earth.


title: Measurement of the inclusive jet cross-section in proton-proton collisions at $\sqrt{s}=7$ TeV using 4.5 fb$^{-1}$ of data with the ATLAS detector

abstract: The inclusive jet cross-section is measured in proton-proton collisions at a centre-of-mass energy of 7 TeV using a data set corresponding to an integrated luminosity of 4.5 fb$^{-1}$ collected with the ATLAS detector at the Large Hadron Collider in 2011. Jets are identified using the anti-$k_t$ algorithm with radius parameter values of 0.4 and 0.6. The double-differential cross-sections are presented as a function of the jet transverse momentum and the jet rapidity, covering jet transverse momenta from 100 GeV to 2 TeV. Next-to-leading-order QCD calculations corrected for non-perturbative effects and electroweak effects, as well as Monte Carlo simulations with next-to-leading-order matrix elements interfaced to parton showering, are compared to the measured cross-sections. A quantitative comparison of the measured cross-sections to the QCD calculations using several sets of parton distribution functions is performed.


title: The Witten-Sakai-Sugimoto model: A brief review and some recent results

abstract: A brief review of the Witten-Sakai-Sugimoto model is given, which is a top-down holographic model of low-energy QCD with chiral quarks derived from type-IIA superstring theory. The main predictions of the model, in particular concerning meson spectra, the gluon condensate, the QCD string tension, the mass of the $\eta'$ and of baryons are discussed and compared quantitatively with available experimental and/or lattice results. Then some recent results of potential interest to the physics program at the future FAIR facility are presented: The spectrum of glueballs and their decay rates into pions, and the phase diagram of QCD at finite temperature, density, and magnetic field strength.


title: Parallel optimization for large eddy simulations

abstract: We developed a parallel Bayesian optimization algorithm for large eddy simulations. These simulations challenge optimization methods because they take hours or days to compute, and their objective function contains noise as turbulent statistics that are averaged over a finite time. Surrogate based optimization methods, including Bayesian optimization, have shown promise for noisy and expensive objective functions. Here we adapt Bayesian optimization to minimize drag in a turbulent channel flow and to design the trailing edge of a turbine blade to reduce turbulent heat transfer and pressure loss. Our optimization simultaneously runs several simulations, each parallelized to thousands of cores, in order to utilize additional concurrency offered by today's supercomputers.


title: String Bits and the Spin Vertex

abstract: We initiate a novel formalism for computing correlation functions of trace operators in the planar N=4 SYM theory. The central object in our formalism is the spin vertex, which is the weak coupling analogy of the string vertex in string field theory. We construct the spin vertex explicitly for all sectors at the leading order using a set of bosonic and fermionic oscillators. We prove that the vertex has trivial monodromy, or put in other words, it is a Yangian invariant. Since the monodromy of the vertex is the product of the monodromies of the three states, the Yangian invariance of the vertex implies an infinite exact symmetry for the three-point function. We conjecture that this infinite symmetry can be lifted to any loop order.


title: Causal Effect Estimation Methods

abstract: Relationship between two popular modeling frameworks of causal inference from observational data, namely, causal graphical model and potential outcome causal model is discussed. How some popular causal effect estimators found in applications of the potential outcome causal model, such as inverse probability of treatment weighted estimator and doubly robust estimator can be obtained by using the causal graphical model is shown. We confine to the simple case of binary outcome and treatment variables with discrete confounders and it is shown how to generalize results to cases of continuous variables.


title: The Corona Problem for Kernel Multiplier Algebras

abstract: We prove an alternate Toeplitz corona theorem for the algebras of pointwise kernel multipliers of Besov-Sobolev spaces on the unit ball in $\mathbb{C}^{n}$, and for the algebra of bounded analytic functions on certain strictly pseudoconvex domains and polydiscs in higher dimensions as well. This alternate Toeplitz corona theorem extends to more general Hilbert function spaces where it does not require the complete Pick property. Instead, the kernel functions $k_{x}\left(y\right)$ of certain Hilbert function spaces $\mathcal{H}$ are assumed to be invertible multipliers on $\mathcal{H}$, and then we continue a research thread begun by Agler and McCarthy in 1999, and continued by Amar in 2003, and most recently by Trent and Wick in 2009. In dimension $n=1$ we prove the corona theorem for the kernel multiplier algebras of Besov-Sobolev Banach spaces in the unit disk, extending the result for Hilbert spaces $H^\infty\cap Q_p$ by A. Nicolau and J. Xiao.


title: Parameter-dependent Gaussian $(z,N)$-generalized Yang-Baxter operators

abstract: We find unitary solutions $\tilde{R}(a)$ to the (multipicative parameter-dependent) $(z,N)$-generalized Yang-Baxter equation that carry the standard measurement basis to $m$-level $N$-partite states that generalize the Bell states corresponding to $\tilde{R}(0)$ in the case $m=N=2$. This is achieved by a careful study of solutions to the Yang-Baxter equation discovered by Fateev and Zamolodchikov in 1982.


title: Greedy Subspace Clustering

abstract: We consider the problem of subspace clustering: given points that lie on or near the union of many low-dimensional linear subspaces, recover the subspaces. To this end, one first identifies sets of points close to the same subspace and uses the sets to estimate the subspaces. As the geometric structure of the clusters (linear subspaces) forbids proper performance of general distance based approaches such as K-means, many model-specific methods have been proposed. In this paper, we provide new simple and efficient algorithms for this problem. Our statistical analysis shows that the algorithms are guaranteed exact (perfect) clustering performance under certain conditions on the number of points and the affinity between subspaces. These conditions are weaker than those considered in the standard statistical literature. Experimental results on synthetic data generated from the standard unions of subspaces model demonstrate our theory. We also show that our algorithm performs competitively against state-of-the-art algorithms on real-world applications such as motion segmentation and face clustering, with much simpler implementation and lower computational cost.


title: Digital Ecology: Coexistence and Domination among Interacting Networks

abstract: The overwhelming success of the web 2.0, with online social networks as key actors, has induced a paradigm shift in the nature of human interactions. The user-driven character of these services for the first time has allowed researchers to quantify large-scale social patterns. However, the mechanisms that determine the fate of networks at a system level are still poorly understood. For instance, the simultaneous existence of numerous digital services naturally raises the question under which conditions these services can coexist. In analogy to population dynamics, the digital world is forming a complex ecosystem of interacting networks whose fitnesses depend on their ability to attract and maintain users' attention, which constitutes a limited resource. In this paper, we introduce an ecological theory of the digital world which exhibits a stable coexistence of several networks as well as the domination of a single one, in contrast to the principle of competitive exclusion. Interestingly, our model also predicts that the most probable outcome is the coexistence of a moderate number of services, in agreement with empirical observations.

categories: physics.soc-ph cs.SI nlin.AO physics.data-an


title: Reconstructing the neutron-star equation of state with gravitational-wave detectors from a realistic population of inspiralling binary neutron stars

abstract: Gravitational-wave observations of inspiralling binary neutron star systems can be used to measure the neutron-star equation of state (EOS) through the tidally induced shift in the waveform phase that depends on the tidal deformability parameter $\lambda$. Previous work has shown that $\lambda$, a function of the neutron-star EOS and mass, is measurable by Advanced LIGO for a single event when including tidal information up to the merger frequency. In this work, we describe a method for stacking measurements of $\lambda$ from multiple inspiral events to measure the EOS. We use Markov chain Monte Carlo simulations to estimate the parameters of a 4-parameter piecewise polytrope EOS that matches theoretical EOS models to a few percent. We find that, for "realistic" event rates ($\sim 40$ binary neutron star inspiral events per year with signal-to-noise ratio $> 8$ in a single Advanced LIGO detector), combining a year of gravitational-wave data from a three-detector network with the constraints from causality and recent high mass neutron-star measurements, the EOS above nuclear density can be measured to better than a factor of two in pressure in most cases. We also find that in the mass range $1M_\odot$--$2M_\odot$, the neutron-star radius can be measured to better than $\pm 1$ km and the tidal deformability can be measured to better than $\pm 1 \times 10^{36}$ g cm$^2$ s$^2$ (10%--50% depending on the EOS and mass). The overwhelming majority of this information comes from the loudest $\sim 5$ events. Current uncertainties in the post-Newtonian waveform model, however, lead to systematic errors in the EOS measurement that are as large as the statistical errors, and more accurate waveform models are needed to minimize this error.


title: On Polar Magnetic Field Reversal and Surface Flux Transport During Solar Cycle 24

abstract: As each solar cycle progresses, remnant magnetic flux from active regions (ARs) migrates poleward to cancel the old-cycle polar field. We describe this polarity reversal process during Cycle 24 using four years (2010.33--2014.33) of line-of-sight magnetic field measurements from the Helioseismic and Magnetic Imager. The total flux associated with ARs reached maximum in the north in 2011, more than two years earlier than the south; the maximum is significantly weaker than Cycle 23. The process of polar field reversal is relatively slow, north-south asymmetric, and episodic. We estimate that the global axial dipole changed sign in October 2013; the northern and southern polar fields (mean above 60$^\circ$ latitude) reversed in November 2012 and March 2014, respectively, about 16 months apart. Notably, the poleward surges of flux in each hemisphere alternated in polarity, giving rise to multiple reversals in the north. We show that the surges of the trailing sunspot polarity tend to correspond to normal mean AR tilt, higher total AR flux, or slower mid-latitude near-surface meridional flow, while exceptions occur during low magnetic activity. In particular, the AR flux and the mid-latitude poleward flow speed exhibit a clear anti-correlation. We discuss how these features can be explained in a surface flux transport process that includes a field-dependent converging flow toward the ARs, a characteristic that may contribute to solar cycle variability.


title: Precinct Size Matters - The Large Precinct Bias in US Presidential Elections

abstract: Examination of precinct level data in US presidential elections reveals a correlation of large precincts and increased fraction of Republican votes. The large precinct bias is analyzed with respect to voter heterogeneity and voter inconvenience as precinct size increases. The analysis shows that voter inconvenience is a significant factor in election outcomes in certain states, and may significantly disadvantage Democratic candidates.


title: Resilience of Social Networks Under Different Attack Strategies

abstract: Recent years have seen the world become a closely connected society with the emergence of different types of social networks. Online social networks have provided a way to bridge long distances and establish numerous communication channels which were not possible earlier. These networks exhibit interesting behavior under intentional attacks and random failures where different structural properties influence the resilience in different ways. In this paper, we perform two sets of experiments and draw conclusions from the results pertaining to the resilience of social networks. The first experiment performs a comparative analysis of four different classes of networks namely small world networks, scale free networks, small world-scale free networks and random networks with four semantically different social networks under different attack strategies. The second experiment compares the resilience of these semantically different social networks under different attack strategies. Empirical analysis reveals interesting behavior of different classes of networks with different attack strategies.


title: Ergodic decompositions for folding and unfolding paths in Outer space

abstract: We relate ergodic-theoretic properties of a very small tree or lamination to the behavior of folding and unfolding paths in Outer space that approximate it, and we obtain a criterion for unique ergodicity in both cases. Our main result is that non-unique ergodicity gives rise to a transverse decomposition of the folding/unfolding path. It follows that non-unique ergodicity leads to distortion when projecting to the complex of free factors, and we give two applications of this fact. First, we show that if a subgroup $H$ of $Out(\FN)$ quasi-isometrically embeds into the complex of free factors via the orbit map, then the limit set of $H$ in the boundary of Outer space consists of trees that are uniquely ergodic and have uniquely ergodic dual lamination. Second, we describe the Poisson boundary for random walks coming from distributions with finite first moment with respect to the word metric on $Out(\FN)$: almost every sample path converges to a tree that is uniquely ergodic and that has a uniquely ergodic dual lamination, and the corresponding hitting measure on the boundary of Outer space is the Poisson boundary. This improves a recent result of Horbez. We also obtain sublinear tracking of sample paths with Lipschitz geodesic rays.


title: Polynomial partitioning for a set of varieties

abstract: Given a set $\Gamma$ of low-degree k-dimensional varieties in $\mathbb{R}^n$, we prove that for any $D \ge 1$, there is a non-zero polynomial $P$ of degree at most $D$ so that each component of $\mathbb{R}^n \setminus Z(P)$ intersects $O(D^{k-n} |\Gamma|)$ varieties of $\Gamma$.


title: Prognosis of Anterior Cruciate Ligament (ACL) Reconstruction: A Data Driven Approach

abstract: Individuals who suffer anterior cruciate ligament (ACL) injury are at higher risk of developing knee osteoarthritis (OA) and almost 50% display symptoms 10 to 20 years post injury. Anterior cruciate ligament reconstruction (ACLR) often does not protect against knee OA development. Accordingly, a multiscale formulation for Data Driven Prognosis (DDP) of post ACLR is developed. Unlike traditional predictive strategies that require controlled off-line measurements or training for determination of constitutive parameters to derive the transitional statistics, the proposed DDP algorithm relies solely on in situ measurements. The proposed DDP scheme is capable of predicting onset of instabilities. Since the need for off line testing (or training) is obviated, it can be easily implemented for ACLR, where such controlled a priori testing is almost impossible to conduct. The DDP algorithm facilitates hierarchical handling of the large data set, and can assess the state of recovery in post ACLR conditions based on data collected from stair ascent and descent exercises of subjects. The DDP algorithm identifies inefficient knee varus motion and knee rotation as primary difficulties experienced by some of the post ACLR population. In such cases, levels of energy dissipation rate at the knee, and its fluctuation may be used as measures for assessing progress after ACL reconstruction.


title: A Quick Empirical Reproof of the Asymptotic Normality of the Hirsch Citation Index (First proved by Canfield, Corteel, and Savage)

abstract: Once upon a time there was an esoteric and specialized notion, called "size of the Durfee square", of interest to at most 100 specialists in the whole world. Then it was kissed by a prince called Jorge Hirsch, and became the famous (and to quite a few people, infamous) h-index, of interest to every scientist, and scholar, since it tells you how productive a scientist (or scholar) you are! When Rodney Canfield, Sylvie Corteel, and Carla Savage wrote their beautiful 1998 article proving, rigorously, by a very deep and intricate analysis, the asymptotic normality of the random variable "size of Durfee square" defined on integer-partitions of n (as n goes to infinity), with precise asymptotics for the mean and variance, they did not dream that one day their result should be of interest to everyone who has ever published a paper. However Canfield et. al. had to work really hard to prove their deep result. Here we take an "empirical" shortcut, that proves the same thing much faster (modulo routine number- and symbol- crunching). More importantly, the empirical methodology should be useful in many other cases where rigorous proofs are either too hard, or not worth the trouble!


title: Cloud structure of the nearest brown dwarfs: Spectroscopic variability of Luhman 16AB from the Hubble Space Telescope

abstract: The binary brown dwarf WISE J104915.57$-$531906.1 (also Luhman 16AB), composed of a late L and early T dwarf, is a prototypical L/T transition flux reversal binary located at only 2 pc distance. Luhman 16B is a known variable whose light curves evolve rapidly. We present spatially resolved spectroscopic time-series of Luhman 16A and B covering 6.5 h using HST/WFC3 at 1.1 to 1.66 $\mu$m. The small, count-dependent variability of Luhman 16A at the beginning of the observations likely stems from instrumental systematics; Luhman 16A appears non-variable above $\approx$0.4%. Its spectrum is well fit by a single cloud layer with intermediate cloud thickness (f_sed=2, Teff=1200 K). Luhman 16B varies at all wavelengths with peak-to-valley amplitudes of 7-11%. The amplitude and light curve shape changes over only one rotation period. The lowest relative amplitude is found in the deep water absorption band at 1.4 $\mu$m, otherwise it mostly decreases gradually from the blue to the red edge of the spectrum. This is very similar to the other two known highly variable early T dwarfs. A two-component cloud model accounts for most of the variability, although small deviations are seen in the water absorption band. We fit the mean spectrum and relative amplitudes with a linear combination of two models of a warm, thinner cloud (Teff=1300 K, fsed=3) and a cooler, thicker cloud (Teff=1000-1100 K, f_sed=1), assuming out-of-equilibrium atmospheric chemistry. A cloud as for Luhman 16A but with holes cannot reproduce the variability of Luhman 16B, indicating more complex cloud evolution through the L/T transition. The projected separation of the binary has decreased by $\approx$0.3'' in 8 months.


title: Measurement of planet masses with transit timing variations due to synodic "chopping" effects

abstract: Gravitational interactions between planets in transiting exoplanetary systems lead to variations in the times of transit that are diagnostic of the planetary masses and the dynamical state of the system. Here we show that synodic "chopping" contributions to these transit timing variations (TTVs) can be used to uniquely measure the masses of planets without full dynamical analyses involving direct integration of the equations of motion. We present simple analytic formulae for the chopping signal, which are valid (generally <10% error) for modest eccentricities e <~ 0.1. Importantly, these formulae primarily depend on the mass of the perturbing planet, and therefore the chopping signal can be used to break the mass/free-eccentricity degeneracy which can appear for systems near first order mean motion resonances. Using a harmonic analysis, we apply these TTV formulae to a number of Kepler systems which had been previously analyzed with full dynamical analyses. We show that when chopping is measured, the masses of both planets can be determined uniquely, in agreement with previous results, but without the need for numerical orbit integrations. This demonstrates how mass measurements from TTVs may primarily arise from an observable chopping signal. The formula for chopping can also be used to predict the number of transits and timing precision required for future observations, such as those made by TESS or PLATO, in order to infer planetary masses through analysis of TTVs.


title: 3D Dirac semimetals: current materials, design principles and predictions of new materials

abstract: Design principles and novel predictions of new 3D Dirac semimetals are presented, along with the context of currently known materials. Current materials include those based on a topological to trivial phase transition, such as in TlBiSe$_{2-x}$S$_x$ and Hg$_{1-x}$Cd$_x$Te, Bi$_{1-x}$Sb$_x$, Bi$_{2-x}$In$_x$Se$_3$, and Pb$_{1-x}$Sn$_x$Se. Some more recently revealed materials, Na$_3$Bi and Cd$_3$As$_2$, are 3D Dirac semimetals in their native composition. The different design principles presented each yield novel predictions for new candidates. For Case I, 3D Dirac semimetals based on charge balanced compounds, BaAgBi, SrAgBi, YbAuSb, PtBi$_2$ and SrSn$_2$As$_2$ are identified as candidates. For Case II, 3D Dirac semi-metals in analogy to graphene, BaGa$_2$ is identified as a candidate, and BaPt and Li$_2$Pt are discussed. For Case III, 3D Dirac semi-metals based on glide planes and screw axes, TlMo$_3$Te$_3$ and the AMo$_3$X$_3$ family in general (A=K, Na, In, Tl, X=Se,Te) as well as the Group IVb trihalides such as HfI$_3$ are identified as candidates. Finally we discuss conventional intermetallic compounds with Dirac cones, and identify Cr$_2$B as a potentially interesting material.


title: Polarized light from Sgr A* in the near-infrared $K_{s}$-band

abstract: We present a statistical analysis of polarized near-infrared (NIR) light from Sgr A*, the radio source associated with the supermassive black hole at the center of the Milky Way. The observations have been carried out using the adaptive optics instrument NACO at the VLT UT4 in the infrared $K_\mathrm{s}$-band from 2004 to 2012. Several polarized flux excursions were observed during these years. Linear polarization at 2.2 $\mu m$, its statistics and time variation, can be used constrain the physical conditions of the accretion process onto this supermassive black hole. With an exponent of about 4 for the number density histogram of fluxes above 5~mJy, the distribution of polarized flux density is closely linked to the single state power-law distribution of the total $K_\mathrm{s}$-band flux densities reported earlier. We find typical polarization degrees of the order of 20\%$\pm$10\% and a preferred polarization angle of 13$^o$$\pm$15$^o$. Simulations show the uncertainties under a total flux density of $\sim 2\,{\rm mJy}$ are probably dominated by observational effects. At higher flux densities there are intrinsic variations of polarization degree and angle within rather well constrained ranges. Since the emission is most likely due to optically thin synchrotron radiation, this preferred polarization angle we find is very likely coupled to the intrinsic orientation of the Sgr A* system i.e. a disk or jet/wind scenario associated with the super massive black hole. If they are indeed linked to structural features of the source the data imply a rather stable geometry and accretion process for the Sgr A* system.


title: Rapid Adaptation of POS Tagging for Domain Specific Uses

abstract: Part-of-speech (POS) tagging is a fundamental component for performing natural language tasks such as parsing, information extraction, and question answering. When POS taggers are trained in one domain and applied in significantly different domains, their performance can degrade dramatically. We present a methodology for rapid adaptation of POS taggers to new domains. Our technique is unsupervised in that a manually annotated corpus for the new domain is not necessary. We use suffix information gathered from large amounts of raw text as well as orthographic information to increase the lexical coverage. We present an experiment in the Biological domain where our POS tagger achieves results comparable to POS taggers specifically trained to this domain.


title: Theory of Transport Phenomena in Coherent Quantum Hall Bilayers

abstract: We argue that all anomalous transport properties of coherent quantum quantum Hall bilayers can be understood in terms of a mean-field transport theory in which the condensate phase is nearly uniform across the sample, and the strength of condensate coupling to interlayer tunneling processes is substantially reduced compared to the predictions of disorder-free microscopic mean-field theory. These ingredients provide a natural explanation for recently established I-V characteristics which feature a critical current above which the tunneling resistance abruptly increases and a non-local interaction between interlayer tunneling at the inner and outer edges of Corbino rings. We propose a microscopic picture in which disorder is the main agent responsible for the reduction of the effective interlayer tunneling strength.


title: Modelling the supernova-driven ISM in different environments

abstract: We use hydrodynamical simulations in a $(256\;{\rm pc})^3$ periodic box to model the impact of supernova (SN) explosions on the multi-phase interstellar medium (ISM) for initial densities $n = 0.5-30$ cm$^{-3}$ and SN rates $1-720$ Myr$^{-1}$. We include radiative cooling, diffuse heating, and the formation of molecular gas using a chemical network. The SNe explode either at random positions, at density peaks, or both. We further present a model combining thermal energy for resolved and momentum input for unresolved SNe. Random driving at high SN rates results in hot gas ($T\gtrsim 10^6$ K) filling $> 90$% of the volume. This gas reaches high pressures ($10^4 < P/k_\mathrm{B} < 10^7$ K cm$^{-3}$) due to the combination of SN explosions in the hot, low density medium and confinement in the periodic box. These pressures move the gas from a two-phase equilibrium to the single-phase, cold branch of the cooling curve. The molecular hydrogen dominates the mass ($>50$%), residing in small, dense clumps. Such a model might resemble the dense ISM in high-redshift galaxies. Peak driving results in huge radiative losses, producing a filamentary ISM with virtually no hot gas, and a small molecular hydrogen mass fraction ($\ll 1$%). Varying the ratio of peak to random SNe yields ISM properties in between the two extremes, with a sharp transition for equal contributions. The velocity dispersion in HI remains $\lesssim 10$ km s$^{-1}$ in all cases. For peak driving the velocity dispersion in H$_\alpha$ can be as high as $70$ km s$^{-1}$ due to the contribution from young, embedded SN remnants.


title: Superfluid and Metamagnetic Phase Transitions in $\omega$-deformed Gauged Supergravity

abstract: We study non-supersymmetric truncations of $\omega$-deformed ${\cal N}=8$ gauged supergravity that retain a $U(1)$ gauge field and three scalars, of which two are neutral and one charged. We construct dyonic domain-wall and black hole solutions with AdS$_4$ boundary conditions when only one (neutral) scalar is non-vanishing, and examine their behavior as the magnetic field and temperature of the system are varied. In the infrared the domain-wall solutions approach either dyonic AdS$_2 \times \mathbb{R}^2$ or else Lifshitz-like, hyperscaling violating geometries. The scaling exponents of the latter are $z=3/2$ and $\theta = -2$, and are independent of the $\omega$-deformation. New $\omega$-dependent AdS$_4$ vacua are also identified. We find a rich structure for the magnetization of the system, including a line of metamagnetic first-order phase transitions when the magnetic field lies in a particular range. Such transitions arise generically in the $\omega$-deformed theories. Finally, we study the onset of a superfluid phase by allowing a fluctuation of the charged scalar field to condense, spontaneously breaking the abelian gauge symmetry. The mechanism by which the superconducting instability ceases to exist for strong magnetic fields is different depending on whether the field is positive or negative. Finally, such instabilities are expected to compete with spatially modulated phases.


title: Real-space quadrature: a convenient, efficient representation for multipole expansions

abstract: Multipolar expansions are a foundational tool for describing basis functions in quantum mechanics, many-body polarization, and other distributions on the unit sphere. Progress on these topics is often held back by complicated and competing formulas for calculating and using spherical harmonics. We present a complete representation for supersymmetric 3D tensors that replaces spherical harmonic basis functions by a dramatically simpler set of weights associated to discrete points in 3D space. This representation is shown to be space optimal. It reduces tensor contraction and the spherical harmonic decomposition of Poisson's operator to pairwise summations over the point set. Moreover, multiplication of spherical harmonic basis functions translates to a direct product in this representation.


title: Metastability and phase separation in a simple model of a superconductor with extremely short coherence length

abstract: We present studies of the atomic limit of the extended Hubbard model with pair hopping for arbitrary electron density and arbitrary chemical potential. The Hamiltonian consists of (i) the effective on-site interaction $U$ and (ii) the intersite charge exchange term $I$, determining the hopping of electron pairs between nearest-neighbour sites. In the analysis of the phase diagrams and thermodynamic properties of this model we treat the intersite interactions within the mean-field approximation. In this report we focus on metastable phases and determine their ranges of occurrence. Our investigations in the absence of the external magnetic field show that the system analysed exhibits tricritical behaviour. Two metastable phases (superconducting and nonordered) can exist inside the regions of the phase separated state stability and a first-order transition occurs between these metastable phases.


title: Insecurity of a relativistic quantum commitment scheme

abstract: We propose a cheating strategy to a relativistic quantum commitment scheme [Sci Rep 2014;4:6774] which was claimed to be unconditionally secure. It is shown that the sender Alice can cheat successfully with probability 100%, thus disproving the security claim.


title: Relating Spatially Resolved Maps of the Schottky Barrier Height to Metal/Semiconductor Interface Composition

abstract: The Schottky barrier height (SBH) is mapped with nanoscale resolution at pure Au/Si(001) and mixed Au/Ag/Si(001) interfaces utilizing ballistic electron emission microscopy (BEEM) by acquiring and fitting spectra every 11.7 nm over a 1 micron by 1 micron area. The energetic distribution of the SBH for the mixed interfaces contain several local maximums indicative of a mixture of metal species at the interface. To estimate the composition at the interface, the distributions are fit to multiple Gaussians that account for the species, "pinch-off" effects, and defects. This electrostatic composition is compared to Rutherford backscattering spectrometry (RBS) and x-ray photo-emission spectroscopy (XPS) measurements to relate it to the physical composition at the interface.


title: The MaxEnt extension of a quantum Gibbs family, convex geometry and geodesics

abstract: We discuss methods to analyze a quantum Gibbs family in the ultra-cold regime where the norm closure of the Gibbs family fails due to discontinuities of the maximum-entropy inference. The current discussion of maximum-entropy inference and irreducible correlation in the area of quantum phase transitions is a major motivation for this research. We extend a representation of the irreducible correlation from finite temperatures to absolute zero.


title: Spin-orbit coupling in fluorinated graphene

abstract: We report on theoretical investigations of the spin-orbit coupling effects in fluorinated graphene. First-principles density functional calculations are performed for the dense and dilute adatom coverage limits. The dense limit is represented by the single-side semifluorinated graphene, which is a metal with spin-orbit splittings of about 10 meV. To simulate the effects of a single adatom, we also calculate the electronic structure of a $10 \times 10$ supercell, with one fluorine atom in the top position. Since this dilute limit is useful to study spin transport and spin relaxation, we also introduce a realistic effective hopping Hamiltonian, based on symmetry considerations, which describes the supercell bands around the Fermi level. We provide the Hamiltonian parameters which are best fits to the first-principles data. We demonstrate that, unlike for the case of hydrogen adatoms, fluorine's own spin-orbit coupling is the principal cause of the giant induced local spin-orbit coupling in graphene. The $sp^3$ hybridization induced transfer of spin-orbit coupling from graphene's $\sigma$ bonds, important for hydrogenated graphene, contributes much less. Furthermore, the magnitude of the induced spin-orbit coupling due to fluorine adatoms is about $1000$ times more than that of pristine graphene, and 10 times more than that of hydrogenated graphene. Also unlike hydrogen, the fluorine adatom is not a narrow resonant scatterer at the Dirac point. The resonant peak in the density of states of fluorinated graphene in the dilute limit lies 260 meV below the Dirac point. The peak is rather broad, about 300 meV, making the fluorine adatom only a weakly resonant scatterer.


title: Rare Event Extinction on Stochastic Networks

abstract: We consider the problem of extinction processes on random networks with a given structure. For sufficiently large well-mixed populations, the process of extinction of one or more state variable components occurs in the tail of the quasi-stationary probability distribution, thereby making it a rare event. Here we show how to extend the theory of large deviations to random networks to predict extinction times. In particular, we use the theory to find the most probable path leading to extinction. We apply the methodology to epidemic models and discover how mean extinction times scale with epidemiological and network parameters in Erdos-Renyi networks. The results are shown to compare quite well with Monte Carlo simulations of the network in predicting both the most probable paths to extinction and mean extinction times.


title: Evidence for weakly bound electrons in non-irradiated alkane crystals. The electrons as a probe of structural differences in crystals

abstract: It is generally assumed that weakly bound (trapped) electrons in organic solids come only from radiolytical (or photochemical) processes like ionization caused by an excited positron entering the sample. This paper presents an evidence for the presence of these electrons in non-irradiated samples of docosane. We argue that these electrons can be located (trapped) either in interlamellar gaps or in spaces made by non-planar conformers. The electrons from the former ones are bound more weakly than those from the latter ones. The origin of Vis absorption for the samples is explained. These spectra can be used as a probe indicating differences in the solid structures of hydrocarbons.


title: Possible restoration of superconductivity in the quasi-one-dimensional conductor Li$_{0.9}$Mo$_6$O$_{17}$ in feasibly high pulsed magnetic fields, $H \simeq 100 \ T$

abstract: We present a theoretical study of restoration of superconductivity in the form of the triplet reentrant superconducting phase in a quasi-one-dimensional (Q1D) conductor. Substitution of known band and superconducting parameters of the presumably triplet Q1D superconductor Li$_{0.9}$Mo$_6$O$_{17}$ into our theoretical equations shows that such restoration can happen in feasibly high non-destructive pulsed magnetic field of the order of $H \simeq 100 \ T$. We investigate in detail how small inclinations of a direction of a magnetic field from its best experimental geometry decrease superconducting transition temperature of the reentrant phase, which is important for its possible experimental discovery.


title: Matrix product states for Hamiltonian lattice gauge theories

abstract: Over the last decade tensor network states (TNS) have emerged as a powerful tool for the study of quantum many body systems. The matrix product states (MPS) are one particular case of TNS and are used for the simulation of 1+1 dimensional systems. In [1] we considered the MPS formalism for the simulation of the Hamiltonian lattice gauge formulation of 1+1 dimensional one flavor quantum electrodynamics, also known as the massive Schwinger model. We deduced the ground state and lowest lying excitations. Furthermore, we performed a full quantum real-time simulation for a quench with a uniform background electric field. In this proceeding we continue our work on the Schwinger model. We demonstrate the advantage of working with gauge invariant MPS by comparing with MPS simulations on the full Hilbert space, that includes numerous non-physical gauge variant states. Furthermore, we compute the chiral condensate and recover the predicted UV-divergent behavior.


title: Dispersion Estimates for One-Dimensional Schr\"odinger and Klein-Gordon Equations Revisited

abstract: We show that for a one-dimensional Schr\"odinger operator with a potential whose first moment is integrable the scattering matrix is in the unital Wiener algebra of functions with integrable Fourier transforms. Then we use this to derive dispersion estimates for solutions of the associated Schr\"odinger and Klein-Gordon equations. In particular, we remove the additional decay conditions in the case where a resonance is present at the edge of the continuous spectrum.


title: Generalized Adaptive Dictionary Learning via Domain Shift Minimization

abstract: Visual data driven dictionaries have been successfully employed for various object recognition and classification tasks. However, the task becomes more challenging if the training and test data are from contrasting domains. In this paper, we propose a novel and generalized approach towards learning an adaptive and common dictionary for multiple domains. Precisely, we project the data from different domains onto a low dimensional space while preserving the intrinsic structure of data from each domain. We also minimize the domain-shift among the data from each pair of domains. Simultaneously, we learn a common adaptive dictionary. Our algorithm can also be modified to learn class-specific dictionaries which can be used for classification. We additionally propose a discriminative manifold regularization which imposes the intrinsic structure of class specific features onto the sparse coefficients. Experiments on image classification show that our approach fares better compared to the existing state-of-the-art methods.


title: Validation of Matching

abstract: We introduce a technique to compute probably approximately correct (PAC) bounds on precision and recall for matching algorithms. The bounds require some verified matches, but those matches may be used to develop the algorithms. The bounds can be applied to network reconciliation or entity resolution algorithms, which identify nodes in different networks or values in a data set that correspond to the same entity. For network reconciliation, the bounds do not require knowledge of the network generation process.


title: Robust sketching for multiple square-root LASSO problems

abstract: Many learning tasks, such as cross-validation, parameter search, or leave-one-out analysis, involve multiple instances of similar problems, each instance sharing a large part of learning data with the others. We introduce a robust framework for solving multiple square-root LASSO problems, based on a sketch of the learning data that uses low-rank approximations. Our approach allows a dramatic reduction in computational effort, in effect reducing the number of observations from $m$ (the number of observations to start with) to $k$ (the number of singular values retained in the low-rank model), while not sacrificing---sometimes even improving---the statistical performance. Theoretical analysis, as well as numerical experiments on both synthetic and real data, illustrate the efficiency of the method in large scale applications.


title: Optically Active Quantum Dots in Monolayer WSe$_2$

abstract: Semiconductor quantum dots have emerged as promising candidates for implementation of quantum information processing since they allow for a quantum interface between stationary spin qubits and propagating single photons. In the meanwhile, transition metal dichalcogenide (TMD) monolayers have moved to the forefront of solid-state research due to their unique band structure featuring a large band gap with degenerate valleys and non-zero Berry curvature. Here we report the observation of quantum dots in monolayer tungsten-diselenide with an energy that is 20 to 100 meV lower than that of two dimensional excitons. Photon antibunching in second-order photon correlations unequivocally demonstrates the zero-dimensional anharmonic nature of these quantum emitters. The strong anisotropic magnetic response of the spatially localized emission peaks strongly indicates that radiative recombination stems from localized excitons that inherit their electronic properties from the host TMD. The large $\sim$ 1 meV zero-field splitting shows that the quantum dots have singlet ground states and an anisotropic confinement most likely induced by impurities or defects in the host TMD. Electrical control in van der Waals heterostructures and robust spin-valley degree of freedom render TMD quantum dots promising for quantum information processing.


title: Sensitivity Increases for the TITAN Decay Spectroscopy Program

abstract: The TITAN facility at TRIUMF has recently initiated a program of performing decay spectroscopy measurements in an electron-beam ion-trap (EBIT). The unique environment of the EBIT provides backing-free storage of the radioactive ions, while guiding charged decay particles from the trap centre via the strong magnetic field. This measurement technique is able to provide a significant increase in detection sensitivity for photons which result from radioactive decay. A brief overview of this device is presented, along with methods of improving the signal-to-background ratio for photon detection by reducing Compton scattered events, and eliminating vibrational noise.


title: Constraining the Dust Coma Properties of Comet C/Siding Spring (2013 A1) at Large Heliocentric Distances

abstract: The close encounter of Comet C/2013 A1 (Siding Spring) with Mars on October 19, 2014 presented an extremely rare opportunity to obtain the first flyby quality data of the nucleus and inner coma of a dynamically new comet. However, the comet's dust tail potentially posed an impact hazard to those spacecraft. To characterize the comet at large heliocentric distances, study its long-term evolution, and provide critical inputs to hazard modeling, we imaged C/Siding Spring with the Hubble Space Telescope when the comet was at 4.58, 3.77, and 3.28 AU from the Sun. The dust production rate, parameterized by the quantity Af$\rho$, was 2500, 2100, and 1700 cm (5000-km radius aperture) for the three epochs, respectively. The color of the dust coma is 5.0$\pm$0.3$\%$/100 nm for the first two epochs, and 9.0$\pm$0.3$\%$/100 nm for the last epoch, and reddens with increasing cometocentric distance out to ~3000 km from the nucleus. The spatial distribution and the temporal evolution of the dust color are most consistent with the existence of icy grains in the coma. Two jet-like dust features appear in the north-northwest and southeast directions projected in the sky plane. Within each epoch of 1-2 hour duration, no temporal variations were observed for either feature, but the PA of the southeastern feature varied between the three epochs by ~30$^\circ$. The dust feature morphology suggests two possible orientations for the rotational pole of the nucleus, (RA, Dec) = (295$^\circ\pm$5$^\circ$, +43$^\circ\pm$2$^\circ$) and (190$^\circ\pm$10$^\circ$, 50$^\circ\pm$5$^\circ$), or their diametrically opposite orientations.


title: Lattice polytopes in coding theory

abstract: In this paper we discuss combinatorial questions about lattice polytopes motivated by recent results on minimum distance estimation for toric codes. We also prove a new inductive bound for the minimum distance of generalized toric codes. As an application, we give new formulas for the minimum distance of generalized toric codes for special lattice point configurations.


title: Effective uniform bounding in partial differential fields

abstract: Motivated by the effective bounds of ordinary differential equations, we prove an effective version of uniform bounding for partial differential fields with commuting derivations. More precisely, we provide an upper bound for the size of finite solution sets of partial differential polynomial equations in terms of data explicitly given in the equations and independent of parameters. Our methods also produce an upper bound for the degree of the Zariski closure of solution sets, whether they are finite or not.


title: A* Sampling

abstract: The problem of drawing samples from a discrete distribution can be converted into a discrete optimization problem. In this work, we show how sampling from a continuous distribution can be converted into an optimization problem over continuous space. Central to the method is a stochastic process recently described in mathematical statistics that we call the Gumbel process. We present a new construction of the Gumbel process and A* sampling, a practical generic sampling algorithm that searches for the maximum of a Gumbel process using A* search. We analyze the correctness and convergence time of A* sampling and demonstrate empirically that it makes more efficient use of bound and likelihood evaluations than the most closely related adaptive rejection sampling-based algorithms.


title: Likelihood-Based Inference for Discretely Observed Birth-Death-Shift Processes, with Applications to Evolution of Mobile Genetic Elements

abstract: Continuous-time birth-death-shift (BDS) processes are frequently used in stochastic modeling, with many applications in ecology and epidemiology. In particular, such processes can model evolutionary dynamics of transposable elements - important genetic markers in molecular epidemiology. Estimation of the effects of individual covariates on the birth, death, and shift rates of the process can be accomplished by analyzing patient data, but inferring these rates in a discretely and unevenly observed setting presents computational challenges. We propose a mutli-type branching process approximation to BDS processes and develop a corresponding expectation maximization (EM) algorithm, where we use spectral techniques to reduce calculation of expected sufficient statistics to low dimensional integration. These techniques yield an efficient and robust optimization routine for inferring the rates of the BDS process, and apply more broadly to multi-type branching processes where rates can depend on many covariates. After rigorously testing our methodology in simulation studies, we apply our method to study intrapatient time evolution of IS6110 transposable element, a frequently used element during estimation of epidemiological clusters of Mycobacterium tuberculosis infections.


title: Modelling the Transfer Function for the Dark Energy Survey

abstract: We present a forward-modelling simulation framework designed to model the data products from the Dark Energy Survey (DES). This forward-model process can be thought of as a transfer function -- a mapping from cosmological and astronomical signals to the final data products used by the scientists. Using output from the cosmological simulations (the Blind Cosmology Challenge), we generate simulated images (the Ultra Fast Image Simulator, Berge et al. 2013) and catalogs representative of the DES data. In this work we simulate the 244 sq. deg coadd images and catalogs in 5 bands for the DES Science Verification (SV) data. The simulation output is compared with the corresponding data to show that major characteristics of the images and catalogs can be captured. We also point out several directions of future improvements. Two practical examples, star/galaxy classification and proximity effects on object detection, are then used to demonstrate how one can use the simulations to address systematics issues in data analysis. With clear understanding of the simplifications in our model, we show that one can use the simulations side-by-side with data products to interpret the measurements. This forward modelling approach is generally applicable for other upcoming and future surveys. It provides a powerful tool for systematics studies which is sufficiently realistic and highly controllable.


title: The Bergman-Shilov boundary for subfamilies of $q$-plurisubharmonic functions

abstract: We introduce a notion of the Bergman-Shilov (or Shilov) boundary for some subclasses of upper-semicontinuous functions on a compact Hausdorff space. It is by definition the smallest closed subset of the given space on which all functions of that subclass attain their maximum. For certain subclasses with simple structure one can show the existence and uniqueness of the Shilov boundary. Then we provide its relation to the set of peak points and establish Bishop-type theorems. As an application we obtain a generalization of Bychkov's theorem which gives a geometric characterization of the Shilov boundary for $q$-plurisubharmonic functions on convex bounded domains. In the case of bounded pesudoconvex domains with smooth boundary we also show that some parts of the Shilov boundary for $q$-plurisubharmonic functions are foliated by $q$-dimensional complex submanifolds.


title: Hypercharge Flux in Heterotic Compactifications

abstract: We study heterotic Calabi-Yau models with hypercharge flux breaking, where the visible E8 gauge group is directly broken to the standard model group by a non-flat gauge bundle, rather than by a two-step process involving an intermediate grand unified theory and a Wilson line. It is shown that the required alternative E8 embeddings of hypercharge, normalized as required for gauge unification, can be found and we classify these possibilities. However, for all but one of these embeddings we prove a general no-go theorem which asserts that no suitable geometry and vector bundle leading to a standard model spectrum can be found. Intuitively, this happens due to the large number of index conditions which have to be imposed in order to obtain a correct physical spectrum in the absence of an underlying grand unified theory.


title: Higgs production constraints on anomalous fermion couplings

abstract: Certain anomalous fermion-gauge boson couplings, such as the flavor diagonal anomalous color magnetic (CMDM) and color electric (CEDM) dipole moments of quarks are not fully gauge invariant under the SM. Restoring gauge invariance with an elementary Higgs doublet implies that they also contribute to Higgs boson production at the LHC and we study the corresponding constraints that can be placed on them. In a similar manner we study the constraints that can be placed on the $\tau$-lepton anomalous magnetic moment, electric dipole moment, weak dipole moments, and dimension eight gluonic couplings at the LHC.


title: General Covariance Constraints on Cosmological Correlators

abstract: We study the extent to which diffeomorphism invariance restricts the properties of the primordial perturbations in single scalar field models. We derive a set of identities that constrain the connected correlators of the cosmological perturbations, as well as the one-particle-irreducible vertices of the theory in any gauge. These identities are the analogues of Slavnov-Taylor identities in gauge theories, and follow essentially from diffeomorphism invariance alone. Yet because quantization requires diffeomorphism invariance to be broken, they not only reflect invariance under diffeomorphisms, but also how the latter has been broken by gauge fixing terms. In order to not lose the symmetry altogether, we cannot simply set some fields to zero, as is usually done in cosmological perturbation theory, but need to decouple them smoothly and make sure that they do not contribute to cosmological correlators in the decoupling limit. We use these identities to derive a set of consistency relations between bispectra and power spectra of cosmological perturbations in different gauges. Without additional assumptions, these consistency relations just seem to reflect the redundancy implied by diffeomorphisms. But when combined with analyticity, in a formulation of the theory in which auxiliary fields have been integrated out, we recover novel and previously derived relations that follow from invariance under both time and spatial diffeomorphisms.


title: Advancing Nucleosynthesis in Self-consistent, Multidimensional Models of Core-Collapse Supernovae

abstract: We investigate core-collapse supernova (CCSN) nucleosynthesis in polar axisymmetric simulations using the multidimensional radiation hydrodynamics code CHIMERA. Computational costs have traditionally constrained the evolution of the nuclear composition in CCSN models to, at best, a 14-species $\alpha$-network. Such a simplified network limits the ability to accurately evolve detailed composition, neutronization and the nuclear energy generation rate. Lagrangian tracer particles are commonly used to extend the nuclear network evolution by incorporating more realistic networks in post-processing nucleosynthesis calculations. Limitations such as poor spatial resolution of the tracer particles, estimation of the expansion timescales, and determination of the "mass-cut" at the end of the simulation impose uncertainties inherent to this approach. We present a detailed analysis of the impact of these uncertainties on post-processing nucleosynthesis calculations and implications for future models.


title: The Far-Ultraviolet Ups and Downs of Alpha Centauri

abstract: Four years (2010-2014) of semiannual pointings by Hubble Space Telescope Imaging Spectrograph (STIS) on nearby Alpha Centauri have yielded a detailed time history of far-ultraviolet emissions of the solar-like primary (A: G2V) and the cooler, but more active, secondary (B: K1V). This period saw A climbing out of a prolonged coronal X-ray minimum, as documented by Chandra, while B was rising to, then falling from, a peak of its long-term (8 yr) starspot cycle. The FUV fluxes of the primary were steady over most of the STIS period, although the [Fe XII] 124 nm coronal forbidden line (T= 1.5 MK) partly mirrored the slowly rising X-ray fluxes. The FUV emissions of the secondary more closely tracked the rise and fall of its coronal luminosities, especially the "hot lines" like Si IV, C IV, and N V (T= 80,000-200,000 K), and coronal [Fe XII] itself. The hot lines of both stars were systematically redshifted, relative to narrow chromospheric emissions, by several km/s, showing little change in amplitude over the 4-year period; especially for Alpha Cen B, despite the significant evolution of its coronal activity. Further, the hot line profiles of both stars, individually and epoch-averaged, could be decomposed into two nearly equal components, one narrow (FWHM~ 25-45 km/s), the other broad (60-80 km/s). Not much variation of the component properties was seen over the 4-year period, even over the major cycle changes of B. This suggests that there is a dominant "quantum" of FUV surface activity that is relatively unchanged during the cycle, aside from the fractional area covered.


title: Conditioning moments of singular measures for entropy maximization II: Numerical examples

abstract: If moments of singular measures are passed as inputs to the entropy maximization procedure, the optimization algorithm might not terminate. The framework developed in our previous paper demonstrated how input moments of measures, on a broad range of domains, can be conditioned to ensure convergence of the entropy maximization. Here we numerically illustrate the developed framework on simplest possible examples: measures with one-dimensional, bounded supports. Three examples of measures are used to numerically compare approximations obtained through entropy maximization with and without the conditioning step.


title: The Slepian zero set, and Brownian bridge embedded in Brownian motion by a spacetime shift

abstract: This paper is concerned with various aspects of the Slepian process $(B_{t+1} - B_t, t \ge 0)$ derived from a one-dimensional Brownian motion $(B_t, t \ge 0 )$. In particular, we offer an analysis of the local structure of the Slepian zero set $\{t : B_{t+1} = B_t \}$, including a path decomposition of the Slepian process for $0 \le t \le 1$. We also establish the existence of a random time $T$ such that $T$ falls in the the Slepian zero set almost surely and the process $(B_{T+u} - B_T, 0 \le u \le 1)$ is standard Brownian bridge.


title: Patterns in random walks and Brownian motion

abstract: We ask if it is possible to find some particular continuous paths of unit length in linear Brownian motion. Beginning with a discrete version of the problem, we derive the asymptotics of the expected waiting time for several interesting patterns. These suggest corresponding results on the existence/non-existence of continuous paths embedded in Brownian motion. With further effort we are able to prove some of these existence and non-existence results by various stochastic analysis arguments. A list of open problems is presented.


title: Internal Proper-Motions in the Eskimo Nebula

abstract: We present measurements of internal proper motions at more than five hundred positions of NGC 2392, the Eskimo Nebula, based on images acquired with WFPC2 on board the Hubble Space Telescope at two epochs separated by 7.695 years. Comparison of the two observations shows clearly the expansion of the nebula. We measured the amplitude and direction of the motion of local structures in the nebula by determining their relative shift during that interval. In order to assess the potential uncertainties in the determination of proper motions in this object, and in general, the measurements were performed using two different methods, used previously in the literature. We compare the results from the two methods, and to perform the scientific analysis of the results we choose one, the cross-correlation method, as the more reliable. We go on to perform a "criss-cross" mapping analysis on the proper motion vectors which helps in the interpretation of the velocity pattern. Combining our results on the proper motions with radial velocity measurements obtained from high resolution spectroscopic observations, and employing an existing 3D model, we estimate the distance to the nebula as 1300 pc.


title: Nematic quantum liquid crystals of bosons in frustrated lattices

abstract: The problem of interacting bosons in frustrated lattices is an intricate one due to the absence of a unique minimum in the single-particle dispersion where macroscopic number of bosons can condense. Here we consider a family of tight-binding models with macroscopically degenerate lowest energy band, separated from other bands by a gap. We predict the formation of exotic states that spontaneously break rotational symmetry at relatively low filling. These states belong to three nematic phases: Wigner crystal, supersolid, and superfluid. The Wigner crystal phase is established exactly at low filling. Supersolid and superfluid phases, at larger filling, are obtained by making use of a projection onto the flat band, construction of an appropriate Wannier basis, and subsequent mean-field treatment. The nematic superfluid that we predict is uniform in real space but has an anisotropic momentum distribution, providing a novel scenario for Bose condensation with an additional nematic order. Our findings open up a promising direction of studying microscopic quantum liquid crystalline phases of bosons.

categories: cond-mat.mes-hall cond-mat.quant-gas cond-mat.soft cond-mat.str-el quant-ph


title: Pulsar timing constraints on narrow-band stochastic signals

abstract: We consider the sensitivity of the pulsar timing array (PTA) technique to specific kind of narrow-band stochastic signals in nano-Hz frequency range. Specifically, we examine the narrow-band signal produced by oscillating gravitational scalar potentials in the Galaxy (Gravitational Potential Background), which arise if an ultralight massive scalar field is the galactic dark matter. We have performed a Bayesian analysis of publicly available data on 12 pulsars obtained by the NANOGrav project. In the monochromatic approximation, the upper limit on the variable gravitational potential amplitude is $\Psi_c<1.14 \times 10^{-15}$, corresponding to the dimensionless strain amplitude $h_c=2\sqrt{3}\Psi_c < 4\times 10^{-15}$ at frequency $f=1.75\times 10^{-8} Hz$. In the narrow-band approximation, the upper limit on the energy density of GPB is found to be $\Omega_{GPB}<1.27 \times 10^{-9}$ at $f=6.2\times 10^{-9} Hz$. These limits are an order of magnitude higher than the theoretically expected values, if the ultralight scalar field with a mass of $\sim 10^{-23} eV$ is assumed to be the galactic dark matter with a local density of $\sim 0.3 GeV cm^{-3}$.


title: Computation of Standardized Residuals for MARSS Models

abstract: This report shows how to compute the variance of the joint conditional model and state residuals for multivariate autoregressive Gaussian state-space (MARSS) models. The MARSS model can be written: x(t)=Bx(t-1)+u+w(t), y(t)=Zx(t)+a+v(t), where w(t) and v(t) are multivariate normal error-terms with variance-covariance matrices Q and R respectively. The joint conditional residuals are the w(t) and v(t) conditioned on a set of, possibly incomplete, data y. Harvey, Koopman and Penzer (1998) show a recursive algorithm for these residuals. I show the equation for the residuals using the conditional variances of the states and the conditional covariance between unobserved data and states. This allows one to compute the variance of un-observed residuals, which could be useful for leave-one-out cross-validation tests. I also show how to modify the Harvey et al algorithm in the case of missing values and how to modify it to return the non-normalized conditional residuals.


title: The Dispersion Relation for Matter Waves in a Two-Phase Vacuum

abstract: The cosmological constant (lambda) of general relativity is a natural consequence of embedding Einstein's theory in a five-dimensional theory of the type needed for unification. The exact 5D solution for lambda less than 0 shows waves in ordinary 3D space with properties similar to those of de Broglie or matter waves. Here the dispersion relation is derived for matter waves in a toy two-phase model, where regions with lambda less than 0 and lambda greater than 0 average on the large scale to lambda = 0, thus providing in principle a resolution of the cosmological-constant problem. A striking result of the analysis is that the dispersion relation is bimodal, with a well-defined window of high-frequency transmission which effectively defines the speed of light.


title: Finite volume study of electric polarizabilities from lattice QCD

abstract: Knowledge of the electric polarizability is crucial to understanding the interactions of hadrons with electromagnetic fields. The neutron polarizability is very sensitive to the quark mass and is expected to diverge in the chiral limit. Here we present results for the electric polarizability of the neutron, neutral pion, and neutral kaon on eight ensembles with nHYP-smeared clover dynamical fermions with two different pion masses (227 and 306 MeV). These are currently the lightest pion masses used in polarizability studies. For each pion mass we compute the polarizability at four different volumes and perform an infinite volume extrapolation for the three hadrons. Along with the infinite volume extrapolation we conduct a chiral extrapolation for the kaon polarizability to the physical point. We compare our results for the neutron polarizability to predictions from chiral perturbation theory.


title: Fusion Tree Sorting

abstract: The sorting problem is one of the most relevant problems in computer science. Within the scope of modern computer science it has been studied for more than 70 years. In spite of these facts, new sorting algorithms have been developed in recent years. Among several types of sorting algorithms, some are quicker; others are more economic in relation to space, whereas others insert a few restrictions in relation to data input. This paper is aimed at explaining the fusion tree data structure, which is responsible for the first sorting algorithm with complexity time smaller than nlgn. The nlgn time complexity has led to some confusion and generated the wrong belief in part of the community of being the minimum possible for this type of problem.


title: New interaction solutions from Lax pair related symmetry of the Generalized fifth order KdV equation

abstract: The nonlocal symmetry of the generalized fifth order KdV equation (FOKdV) is first obtained by using the related Lax pair and then localized in a new enlarged system by introducing some new variables. On this basis, new Backlund transformation is obtained through Lie's first theorem. Furthermore, the general form of Lie point symmetry for the enlarged FOKdV system is found and new interaction solutions for the FOKdV equation are explored by using classical symmetry reduction method.


title: A Suzaku Search for Dark Matter Emission Lines in the X-ray Brightest Galaxy Clusters

abstract: We present the results of a search for unidentified emission lines in deep Suzaku X-ray spectra for the central regions of the four X-ray brightest galaxy clusters: Perseus, Coma, Virgo and Ophiuchus. We employ an optimized energy range for our analysis (3.2-5.3 keV) that is relatively free of instrumental features, and a baseline plasma emission model that incorporates the abundances of elements with the strongest expected emission lines at these energies (S, Ar, Ca) as free parameters. For the Perseus Cluster core, employing this baseline model, we find evidence for an additional emission feature at an energy $3.51^{+0.02}_{-0.01}$ keV with a flux of ~$2.87\times10^{-7}$ ph/s/cm^2/arcmin^2. At slightly larger radii, we detect an emission line at 3.59+/-0.02 keV with a flux of ~$4.8\times10^{-8}$ ph/s/cm^2/arcmin^2. The energies and fluxes of these features are broadly consistent with previous claims, although the radial variation of the line strength appears in tension with standard dark matter (DM) model predictions. Assuming a decaying DM origin for the Perseus emission features allows us to predict the energies and line fluxes for the other clusters in our sample. Critically, we do not detect an emission feature at the predicted energy and line flux in the Coma, Virgo and Ophiuchus clusters. The formal 99.5 per cent upper limits on the strengths of an emission line in each cluster are below the decaying DM model predictions, scaling from the Perseus Cluster center, apparently ruling the model out. In the light of these results, we search for other explanations for the ~3.55 keV emission feature in Perseus. Our results suggest that systematic effects associated with modeling the complex spectra for the Perseus Cluster core, details of the assumed ionization balance, and errors in the predicted emissivities of the spectral lines may in part be responsible for the ~3.5 keV feature.


title: Thick-Target Neutron Yield from the 19F(alpha,n) Reaction

abstract: Thick-target neutron yields from the 19F(alpha,n) reaction are reported for E(alpha) = 3.5 - 10.0 MeV.


title: ContactTrees: A Technique for Studying Personal Network Data

abstract: Network visualization allows a quick glance at how nodes (or actors) are connected by edges (or ties). A conventional network diagram of "contact tree" maps out a root and branches that represent the structure of nodes and edges, often without further specifying leaves or fruits that would have grown from small branches. By furnishing such a network structure with leaves and fruits, we reveal details about "contacts" in our ContactTrees that underline ties and relationships. Our elegant design employs a bottom-up approach that resembles a recent attempt to understand subjective well-being by means of a series of emotions. Such a bottom-up approach to social-network studies decomposes each tie into a series of interactions or contacts, which help deepen our understanding of the complexity embedded in a network structure. Unlike previous network visualizations, ContactTrees can highlight how relationships form and change based upon interactions among actors, and how relationships and networks vary by contact attributes. Based on a botanical tree metaphor, the design is easy to construct and the resulting tree-like visualization can display many properties at both tie and contact levels, a key ingredient missing from conventional techniques of network visualization. We first demonstrate ContactTrees using a dataset consisting of three waves of 3-month contact diaries over the 2004-2012 period, then compare ContactTrees with alternative tools and discuss how this tool can be applied to other types of datasets.


title: Renormalization constants for $N_{\rm f}=2+1+1$ twisted mass QCD

abstract: We summarize recent non-perturbative results obtained for the renormalization constants computed in the RI'-MOM scheme for $N_{\rm f}=2+1+1$ twisted mass QCD. Our implementation employs the Iwasaki gauge action and four dynamical degenerate twisted mass fermions. Renormalization constants for scalar, pseudo-scalar, vector and axial operators, as well as the quark propagator renormalization, are computed at three different values of the lattice spacing, two different volumes and several values of the twisted mass. Our method allows for a precise cross-check of the running, because of the particular proper treatment of the hypercubic artifacts. Preliminary results for twist-2 operators are also presented.


title: Cognitive image processing: the time is right to recognize that the world does not rest more on turtles and elephants

abstract: Traditional image processing is a field of science and technology developed to facilitate human-centered image management. But today, when huge volumes of visual data inundate our surroundings (due to the explosive growth of image-capturing devices, proliferation of Internet communication means and video sharing services over the World Wide Web), human-centered handling of Big-data flows is impossible anymore. Therefore, it has to be replaced with a machine (computer) supported counterpart. Of course, such an artificial counterpart must be equipped with some cognitive abilities, usually characteristic for a human being. Indeed, in the past decade, a new computer design trend - Cognitive Computer development - is become visible. Cognitive image processing definitely will be one of its main duties. It must be specially mentioned that this trend is a particular case of a much more general movement - the transition from a "computational data-processing paradigm" to a "cognitive information-processing paradigm", which affects today many fields of science, technology, and engineering. This transition is a blessed novelty, but its success is hampered by the lack of a clear delimitation between the notion of data and the notion of information. Elaborating the case of cognitive image processing, the paper intends to clarify these important research issues.


title: Dipoles and Pixie dust

abstract: Every closed subset of the Riemann sphere can be approximated in the Hausdorff topology by the Julia set of a rational map.


title: Semiconducting Black Phosphorus: Synthesis, Transport Properties and Electronic Applications

abstract: Phosphorus is one of the most abundant elements preserved in earth, constructing with a fraction of ~0.1% of the earth crust. In general, phosphorus has several allotropes. The two most commonly seen allotropes, white and red phosphorus, are widely used in explosives and safety matches. In addition, black phosphorus, though rarely mentioned, is a layered semiconductor and have great potentials in optical and electronic applications. Remarkably, this layered material can be reduced to one single atomic layer in the vertical direction owing to the van der Waals structure, known as phosphorene, where the physical properties can be tremendously different from its bulk counterpart. In this review article, we trace back to the 100 years research history on black phosphorus from the synthesis to material properties, and extend the topic from black phosphorus to phosphorene. The physical and transport properties are highlighted, aiming at further applications in electronic and optoelectronics devices.


title: Complex Hadamard matrices contained in a Bose-Mesner algebra

abstract: A complex Hadamard matrix is a square matrix H with complex entries of absolute value 1 satisfying $HH^*= nI$, where $*$ stands for the Hermitian transpose and I is the identity matrix of order $n$. In this paper, we first determine the image of a certain rational map from the $d$-dimensional complex projective space to $\mathbb{C}^{d(d+1)/2}$. Applying this result with $d=3$, we give constructions of complex Hadamard matrices, and more generally, type-II matrices, in the Bose-Mesner algebra of a certain 3-class symmetric association scheme. In particular, we recover the complex Hadamard matrices of order 15 found by Ada Chan. We compute the Haagerup sets to show inequivalence of resulting type-II matrices, and determine the Nomura algebras to show that the resulting matrices are not decomposable into generalized tensor products.


title: Ferromagnetic interactions and martensitic transformation in Fe doped Ni-Mn-In shape memory alloys

abstract: The structure, magnetic and martensitic properties of Fe doped Ni-Mn-In magnetic shape memory alloys have been studied by differential scanning calorimetry, magnetization, resistivity, X-ray diffraction (XRD) and EXAFS. While Ni$_{2}$MnIn$_{1-x}$Fe$_{x}$ ($0 \le x \le 0.6$) alloys are ferromagnetic and non martensitic, the martensitic transformation temperature in Ni$_{2}$Mn$_{1.5} $In$_{1-y}$Fe$_{y}$ and Ni$_{2}$Mn$ _{1.6} $In$_{1-y}$Fe$_{y}$ increases for lower Fe concentrations ($y \le 0.05$) before decreasing sharply for higher Fe concentrations. XRD analysis reveals presence of cubic and tetragonal structural phases in Ni$_{2}$MnIn$_{1-x}$Fe$_{x}$ at room temperature with tetragonal phase content increasing with Fe doping. Even though the local structure around Mn and Ni in these Fe doped alloys is similar to martensitic Mn rich Ni-Mn-In alloys, presence of ferromagnetic interactions and structural disorder induced by Fe affect Mn-Ni-Mn antiferromagnetic interactions resulting in suppression of martensitic transformation in these Fe doped alloys.


title: The Burden of Risk Aversion in Mean-Risk Selfish Routing

abstract: Considering congestion games with uncertain delays, we compute the inefficiency introduced in network routing by risk-averse agents. At equilibrium, agents may select paths that do not minimize the expected latency so as to obtain lower variability. A social planner, who is likely to be more risk neutral than agents because it operates at a longer time-scale, quantifies social cost with the total expected delay along routes. From that perspective, agents may make suboptimal decisions that degrade long-term quality. We define the {\em price of risk aversion} (PRA) as the worst-case ratio of the social cost at a risk-averse Wardrop equilibrium to that where agents are risk-neutral. For networks with general delay functions and a single source-sink pair, we show that the PRA depends linearly on the agents' risk tolerance and on the degree of variability present in the network. In contrast to the {\em price of anarchy}, in general the PRA increases when the network gets larger but it does not depend on the shape of the delay functions. To get this result we rely on a combinatorial proof that employs alternating paths that are reminiscent of those used in max-flow algorithms. For {\em series-parallel} (SP) graphs, the PRA becomes independent of the network topology and its size. As a result of independent interest, we prove that for SP networks with deterministic delays, Wardrop equilibria {\em maximize} the shortest-path objective among all feasible flows.


title: Secrecy in Cascade Networks

abstract: We consider a cascade network where a sequence of nodes each send a message to their downstream neighbor to enable coordination, the first node having access to an information signal. An adversary also receives all of the communication as well as additional side-information. The performance of the system is measured by a payoff function evaluated on actions produced at each of the nodes, including the adversary. The challenge is to effectively use a secret key to infuse some level of privacy into the encoding, in order thwart the adversary's attempt to reduce the payoff. We obtain information-theoretic inner and outer bounds on performance, and give examples where they are tight. From these bounds, we also derive the optimal equivocation for this setting as a special case.


title: Component periods of non-principal-axis rotation and their manifestations in the lightcurves of asteroids and bare cometary nuclei

abstract: Lightcurve observations of asteroids and bare cometary nuclei are the most widely used observational tool to derive the rotational parameters. Therefore, an in-depth understanding of how component periods of dynamically excited non-principal axis (NPA) rotators manifest in lightcurves is a crucial step in this process. We investigated this with the help of numerically generated lightcurves of NPA rotators with component periods known a priori. The component periods of NPA rotation were defined in terms of two widely used yet complementary conventions. We derive the relationships correlating the component rotation periods in the two conventions. These relationships were then used to interpret the periodicity signatures present in the simulated lightcurves and rationalize them in either convention.


title: Algorithms for Maximum Agreement Forest of Multiple General Trees

abstract: The Maximum Agreement Forest (Maf) problem is a well-studied problem in evolutionary biology, which asks for a largest common subforest of a given collection of phylogenetic trees with identical leaf label-set. However, the previous work about the Maf problem are mainly on two binary phylogenetic trees or two general (i.e., binary and non-binary) phylogenetic trees. In this paper, we study the more general version of the problem: the Maf problem on multiple general phylogenetic trees. We present a parameterized algorithm of running time $O(3^k n^2m)$ and a 3-approximation algorithm for the Maf problem on multiple rooted general phylogenetic trees, and a parameterized algorithm of running time $O(4^k n^2m)$ and a 4-approximation algorithm for the Maf problem on multiple unrooted general phylogenetic trees. We also implement the parameterized algorithm and approximation algorithm for the Maf problem on multiple rooted general phylogenetic trees, and test them on simulated data and biological data.


title: Top quark pair production and top quark properties at CDF

abstract: We present the most recent measurements of top quark pairs production and top quark properties in proton-antiproton collisions with center-of-mass energy of 1.96 TeV using CDF II detector at the Tevatron. The combination of top pair production cross section measurements and the direct measurement of top quark width are reported. The test of Standard Model predictions for top quark decaying into $b$-quarks, performed by measuring the ratio $R$ between the top quark branching fraction to $b$-quark and the branching fraction to any type of down quark is shown. The extraction of the CKM matrix element $|V_{tb}|$ from the ratio $R$ is discussed. We also present the latest measurements on the forward-backward asymmetry ($A_{FB}$) in top anti-top quark production. With the full CDF Run II data set, the measurements are performed in top anti-top decaying to final states that contain one or two charged leptons (electrons or muons). In addition, we combine the results of the leptonic forward-backward asymmetry in $t\bar t$ system between the two final states. All the results show deviations from the next-to-leading order (NLO) standard model (SM) calculation.


title: ALID: Scalable Dominant Cluster Detection

abstract: Detecting dominant clusters is important in many analytic applications. The state-of-the-art methods find dense subgraphs on the affinity graph as the dominant clusters. However, the time and space complexity of those methods are dominated by the construction of the affinity graph, which is quadratic with respect to the number of data points, and thus impractical on large data sets. To tackle the challenge, in this paper, we apply Evolutionary Game Theory (EGT) and develop a scalable algorithm, Approximate Localized Infection Immunization Dynamics (ALID). The major idea is to perform Localized Infection Immunization Dynamics (LID) to find dense subgraph within local range of the affinity graph. LID is further scaled up with guaranteed high efficiency and detection quality by an estimated Region of Interest (ROI) and a carefully designed Candidate Infective Vertex Search method (CIVS). ALID only constructs small local affinity graphs and has a time complexity of O(C(a^*+ {\delta})n) and a space complexity of O(a^*(a^*+ {\delta})), where a^* is the size of the largest dominant cluster and C << n and {\delta} << n are small constants. We demonstrate by extensive experiments on both synthetic data and real world data that ALID achieves state-of-the-art detection quality with much lower time and space cost on single machine. We also demonstrate the encouraging parallelization performance of ALID by implementing the Parallel ALID (PALID) on Apache Spark. PALID processes 50 million SIFT data points in 2.29 hours, achieving a speedup ratio of 7.51 with 8 executors.


title: Hlawka-Popoviciu inequalities on positive definite tensors

abstract: We prove inequalities on symmetric tensor sums of positive definite operators. In particular, we prove multivariable operator inequalities inspired by generalizations to the well-known Hlawka and Popoviciu inequalities. As corollaries, we obtain generalized Hlawka and Popoviciu inequalities for determinants, permanents, and generalized matrix functions. The new operator inequalities and their corollaries contain a few recently published inequalities on positive definite matrices as special cases.


title: Analysis of the divide-and-conquer method for electronic structure calculations

abstract: We study the accuracy of the divide-and-conquer method for electronic structure calculations. The analysis is conducted for a prototypical subdomain problem in the method. We prove that the pointwise difference between electron densities of the global system and the subsystem decays exponentially as a function of the distance away from the boundary of the subsystem, under the gap assumption of both the global system and the subsystem. We show that gap assumption is crucial for the accuracy of the divide-and-conquer method by numerical examples. In particular, we show examples with the loss of accuracy when the gap assumption of the subsystem is invalid.


title: A transformed framework for dynamic correlation in multireference problems

abstract: We describe how multirefence dynamic correlation theories can be naturally obtained as single-reference correlation theories in a canonically transformed frame. Such canonically transformed correlation theories are very simple and involve identical expressions to their single-reference counterparts. The corresponding excitations involve quasiparticles rather than the bare particles of the system. High-order density matrices (or their approximations) and the numerical metric instabilities common to multireference correlation theories do not appear. As an example, we formulate the Bogoliubov canonically transformed version of second-order M{\o}ller-Plesset perturbation theory and demonstrate its performance in hydrogen, water, nitrogen, and BeH$_2$ bond dissociation.


title: Gradient Flow Analysis on MILC HISQ Ensembles

abstract: We report on a preliminary scale determination with gradient-flow techniques on the $N_f = 2 + 1 + 1$ HISQ ensembles generated by the MILC collaboration. The ensembles include four lattice spacings, ranging from 0.15 to 0.06 fm, and both physical and unphysical values of the quark masses. The scales $\sqrt{t_0}/a$ and $w_0/a$ are computed using Symanzik flow and the cloverleaf definition of $\langle E \rangle$ on each ensemble. Then both scales and the meson masses $aM_\pi$ and $aM_K$ are adjusted for mistunings in the charm mass. Using a combination of continuum chiral perturbation theory and a Taylor series ansatz in the lattice spacing, the results are simultaneously extrapolated to the continuum and interpolated to physical quark masses. Our preliminary results are $\sqrt{t_0} = 0.1422(7)$fm and $w_0 = 0.1732(10)$fm. We also find the continuum mass-dependence of $w_0$.