title: Red Supergiants as Cosmic Abundance Probes: The Magellanic Clouds

abstract: Red Supergiants (RSGs) are cool (~4000K), highly luminous stars (L - 10^5 Lsun), and are among the brightest near-infrared (NIR) sources in star-forming galaxies. This makes them powerful probes of the properties of their host galaxies, such as kinematics and chemical abundances. We have developed a technique whereby metallicities of RSGs may be extracted from a narrow spectral window around 1{\mu}m from only moderate resolution data. The method is therefore extremely efficient, allowing stars at large distances to be studied, and so has tremendous potential for extragalactic abundance work. Here, we present an abundance study of the Large and Small Magellanic Clouds (LMC and SMC respectively) using samples of 9-10 RSGs in each. We find average abundances for the two galaxies of [Z]LMC = -0.37 +/- 0.14 and [Z]SMC = -0.53 +/- 0.16 (with respect to a Solar metallicity of Zsun=0.012). These values are consistent with other studies of young stars in these galaxies, and though our result for the SMC may appear high it is consistent with recent studies of hot stars which find 0.5-0.8dex below Solar. Our best-fit temperatures are on the whole consistent with those from fits to the optical-infrared spectral energy distributions, which is remarkable considering the narrow spectral range being studied. Combined with our recent study of RSGs in the Galactic cluster Per OB1, these results indicate that this technique performs well over a range of metallicities, paving the way for forthcoming studies of more distant galaxies beyond the Local Group.


title: On the chirality of the SM and the fermion content of GUTs

abstract: The Standard Model (SM) is a chiral theory, where right- and left-handed fermion fields transform differently under the gauge group. Extra fermions, if they do exist, need to be heavy otherwise they would have already been observed. With no complex mechanisms at work, such as confining interactions or extra-dimensions, this can only be achieved if every extra right-handed fermion comes paired with a left-handed one transforming in the same way under the Standard Model gauge group, otherwise the new states would only get a mass after electroweak symmetry breaking, which would necessarily be small ($\sim100\textrm{ GeV}$). Such a simple requirement severely constrains the fermion content of Grand Unified Theories (GUTs). It is known for example that three copies of the representations $\mathbf{\overline{5}}+\mathbf{10}$ of $SU(5)$ or three copies of the $\mathbf{16}$ of $SO(10)$ can reproduce the Standard Model's chirality, but how unique are these arrangements? In a systematic way, this paper looks at the possibility of having non-standard mixtures of fermion GUT representations yielding the correct Standard Model chirality. Family unification is possible with large special unitary groups --- for example, the $\mathbf{171}$ representation of $SU(19)$ may decompose as $3\left(\mathbf{16}\right)+\mathbf{120}+3\left(\mathbf{1}\right)$ under $SO(10)$.


title: Single photoelectron spin detection and angular momentum transfer in a gate defined quantum dot

abstract: Recent innovations in fabricating nanoscale confined spin systems have enabled investigation of fundamental quantum correlations between single quanta of photons and matter states. Realization of quantum state transfer from photon polarization to electron spin using gate defined quantum dots (QDs) may give evidence of preserved coherence of angular momentum basis states at the photon-spin interface. The interface would enlarge the concept of quantum information technology, in which single photogenerated electron spins are manipulated with the dots, but this remains a serious challenge. Here, we report the detection of single electron spins generated by polarized single photons via a double QD (DQD) to verify the angular momentum transfer from single photons to single electrons. Pauli spin blockade (PSB) is used to project the photoelectron spin state onto the up or down spin state. Our result promises the realization of coherent quantum state transfer and development of hybrid photon and spin quantum technology.


title: A dwarf galaxy's transformation and a massive galaxy's edge: autopsy of kill and killer in NGC 1097

abstract: (abridged) We present a dynamical analysis of the extended stellar stream encircling NGC 1097. Within a statistical framework, we model its surface brightness using mock streams as in Amorisco (2015) and deep imaging data from the CHART32 telescope (Stellar Tidal Stream Survey). We reconstruct the post-infall evolution of the progenitor, which has experienced 3 pericentric passages and lost more than 2 orders of magnitude in mass. At infall, $5.4\pm0.6$ Gyr ago, the progenitor was a disky dwarf with mass of $\log_{10}[m(<3.4\pm1 {\rm kpc})/ M_\odot]=10.35\pm0.25$. We illustrate how the 90$^\circ$ turn in the stream, identifying the `dog leg', is the signature of the progenitor's prograde rotation. Today, the remnant is a nucleated dwarf, with a LOS velocity of $v_{\rm p, los}^{\rm obs}=-30\pm 30$ kms$^{-1}$, and a luminosity of $3.3\times 10^7 L_{V,\odot}$ (Galianni et al. 2010). Our independent analysis predicts $v_{\rm p, los}=-51^{-17}_{+14}$ kms$^{-1}$, and measures $\log_{10}(m/ M_\odot)=7.4^{+0.6}_{-0.8}$, so that the compact nucleus is soon becoming a low-luminosity UCD. We find that NGC 1097 has a mass of $M_{200}=1.8^{+0.5}_{-0.4} \times 10^{12}\; M_{\odot}$, and its concentration $c_{200}=6.7^{+2.4}_{-1.3}$ is in agreement with LCDM. The stream is described almost down to the noise in a spherical host potential, we find this would not be possible if the halo was substantially triaxial at large radii. Its morphology shows that the slope of the total density profile bends from an inner $\gamma(r_{\rm peri})=1.5\pm0.15$. The progenitor's orbit reaches $r_{\rm apo}=150\pm 15$ kpc, more than a half of the virial radius of the host, so that, for the first time on an individual extragalactic halo, we measure the outer density slope, $\gamma(0.6r_{200,c})=3.9\pm0.5$. This demonstrates the promise of the newborn field of detailed, statistical modelling of extragalactic tidal streams.


title: A topologically twisted index for three-dimensional supersymmetric theories

abstract: We provide a general formula for the partition function of three-dimensional $\mathcal{N}=2$ gauge theories placed on $S^2 \times S^1$ with a topological twist along $S^2$, which can be interpreted as an index for chiral states of the theories immersed in background magnetic fields. The result is expressed as a sum over magnetic fluxes of the residues of a meromorphic form which is a function of the scalar zero-modes. The partition function depends on a collection of background magnetic fluxes and fugacities for the global symmetries. We illustrate our formula in many examples of 3d Yang-Mills-Chern-Simons theories with matter, including Aharony and Giveon-Kutasov dualities. Finally, our formula generalizes to $\Omega$-backgrounds, as well as two-dimensional theories on $S^2$ and four-dimensional theories on $S^2 \times T^2$. In particular this provides an alternative way to compute genus-zero A-model topological amplitudes and Gromov-Witten invariants.


title: Dust within the old nuclear star cluster in the Milky Way

abstract: The mean absolute extinction towards the central parsec of the Milky Way is A_K~3 mag, including both foreground and Galactic center dust. Here we present a measurement of dust extinction within the Galactic old nuclear star cluster (NSC), based on combining differential extinctions of NSC stars with their u_l proper motions along Galactic longitude. Extinction within the NSC preferentially affects stars at its far side, and because the NSC rotates, this causes higher extinctions for NSC stars with negative u_l, as well as an asymmetry in the u_l-histograms. We model these effects using an axisymmetric dynamical model of the NSC in combination with simple models for the dust distribution. Comparing the predicted asymmetry to data for ~7100 stars in several NSC fields, we find that dust associated with the Galactic center mini-spiral with extinction A_K~=0.15-0.8 mag explains most of the data. The largest extinction A_K~=0.8 mag is found in the region of the Western arm of the mini-spiral. Comparing with total A_K determined from stellar colors, we determine the extinction in front of the NSC. Finally, we estimate that for a typical extinction of A_K~=0.4 the statistical parallax of the NSC changes by ~0.4%.


title: Magnetic Untwisting in Solar Jets that Go into the Outer Corona in Polar Coronal Holes

abstract: We study 14 large solar jets observed in polar coronal holes. In EUV movies from SDO/AIA, each jet appears similar to most X-ray jets and EUV jets that erupt in coronal holes, but each is exceptional in that it goes higher than most, so high that it is observed in the outer corona beyond 2.2 RSun in images from the SOHO/LASCO/C2 coronagraph. From AIA He II 304 {\AA} movies and LASCO/C2 running-difference images of these high-reaching jets, we find: (1) the front of the jet transits the corona below 2.2 RSun at a speed typically several times the sound speed; (2) each jet displays an exceptionally large amount of spin as it erupts; (3) in the outer corona, most of the jets display measureable swaying and bending of a few degrees in amplitude; in three jets the swaying is discernibly oscillatory with a period of order 1 hour. These characteristics suggest that the driver in these jets is a magnetic-untwisting wave that is basically a large-amplitude (i.e., non-linear) torsional Alfven wave that is put into the reconnected open field in the jet by interchange reconnection as the jet erupts. From the measured spinning and swaying we estimate that the magnetic-untwisting wave loses most of its energy in the inner corona below 2.2 RSun. We point out that the torsional waves observed in Type-II spicules might dissipate in the corona in the same way as the magnetic-untwisting waves in our big jets and thereby power much of the coronal heating in coronal holes.


title: Probabilistic Clustering of Time-Evolving Distance Data

abstract: We present a novel probabilistic clustering model for objects that are represented via pairwise distances and observed at different time points. The proposed method utilizes the information given by adjacent time points to find the underlying cluster structure and obtain a smooth cluster evolution. This approach allows the number of objects and clusters to differ at every time point, and no identification on the identities of the objects is needed. Further, the model does not require the number of clusters being specified in advance -- they are instead determined automatically using a Dirichlet process prior. We validate our model on synthetic data showing that the proposed method is more accurate than state-of-the-art clustering methods. Finally, we use our dynamic clustering model to analyze and illustrate the evolution of brain cancer patients over time.


title: Energy Injection in Gamma-ray Burst Afterglows

abstract: We present multi-wavelength observations and modeling of Gamma-ray Bursts (GRBs) that exhibit a simultaneous re-brightening in their X-ray and optical light curves, and are also detected at radio wavelengths. We show that the re-brightening episodes can be modeled by injection of energy into the blastwave and that in all cases the energy injection rate falls within the theoretical bounds expected for a distribution of energy with ejecta Lorentz factor. Our measured values of the circumburst density, jet opening angle, and beaming corrected kinetic energy are consistent with the distribution of these parameters for long-duration GRBs at both z~1 and z>6, suggesting that the jet launching mechanism and environment of these events are similar to that of GRBs that do not have bumps in their light curves. However, events exhibiting re-brightening episodes have lower radiative efficiencies than average, suggesting that a majority of the kinetic energy of the outflow is carried by slow-moving ejecta, which is further supported by steep measured distributions of the ejecta energy as a function of Lorentz factor. We do not find evidence for reverse shocks over the energy injection period, implying that the onset of energy injection is a gentle process. We further show that GRBs exhibiting simultaneous X-ray and optical re-brightenings are likely the tail of a distribution of events with varying rates of energy injection, forming the most extreme events in their class. Future X-ray observations of GRB afterglows with Swift and its successors will thus likely discover several more such events, while radio follow-up and multi-wavelength modeling of similar events will unveil the role of energy injection in GRB afterglows.


title: Long-distance entanglement distribution using individual atoms in optical cavities

abstract: Individual atoms in optical cavities can provide an efficient interface between stationary qubits and flying qubits (photons), which is an essentiel building block for quantum communication. Furthermore, cavity assisted controlled-not (CNOT) gates can be used for swapping entanglement to long distances in a quantum repeater setup. Nonetheless, dissipation introduced by the cavity during the CNOT may increase the experimental difficulty in obtaining long-distance entanglement distribution using these systems. We analyse and compare a number of cavity-based repeater schemes combining various entanglement generation schemes and cavity assisted CNOT gates. We find that a scheme, where high-fidelity entanglement is first generated in a two-photon detection scheme and then swapped to long distances using a recently proposed heralded CZ-gate exhibits superior performance compared to the other schemes. The heralded gate moves the effect of dissipation from the fidelity to the success probability of the gate thereby enabling high-fidelity entanglement swapping. As a result, high-rate entanglement distribution can be achieved over long distances even for low cooperativities of the atom-cavity systems. This high-fidelity repeater is shown to outperform the other cavity-based schemes by up to two orders of magnitude in the rate for realistic parameters and large distances (1000 km).


title: Density-matrix based determination of low-energy model Hamiltonians from ab initio wavefunctions

abstract: We propose a way of obtaining effective low energy Hubbard-like model Hamiltonians from ab initio Quantum Monte Carlo calculations for molecular and extended systems. The Hamiltonian parameters are fit to best match the ab initio two-body density matrices and energies of the ground and excited states, and thus we refer to the method as ab initio density matrix based downfolding (AIDMD). For benzene (a finite system), we find good agreement with experimentally available energy gaps without using any experimental inputs. For graphene, a two dimensional solid (extended system) with periodic boundary conditions, we find the effective on-site Hubbard $U^{*}/t$ to be $1.3 \pm 0.2$, comparable to a recent estimate based on the constrained random phase approximation. For molecules, such parameterizations enable calculation of excited states that are usually not accessible within ground state approaches. For solids, the effective Hamiltonian enables large scale calculations using techniques designed for lattice models.


title: Racah Polynomials and Recoupling Schemes of $\mathfrak{su}(1,1)$

abstract: The connection between the recoupling scheme of four copies of $\mathfrak{su}(1,1)$, the generic superintegrable system on the 3 sphere, and bivariate Racah polynomials is identified. The Racah polynomials are presented as connection coefficients between eigenfunctions separated in different spherical coordinate systems and equivalently as different irreducible decompositions of the tensor product representations. As a consequence of the model, an extension of the quadratic algebra ${\rm QR}(3)$ is given. It is shown that this algebra closes only with the inclusion of an additional shift operator, beyond the eigenvalue operators for the bivariate Racah polynomials, whose polynomial eigenfunctions are determined. The duality between the variables and the degrees, and hence the bispectrality of the polynomials, is interpreted in terms of expansion coefficients of the separated solutions.


title: Higgs Dark Matter from a Warped Extra-Dimension -- the truncated-inert-doublet model

abstract: We construct a 5D $\mathbb{Z}_2$-symmetric model with three D3-branes: two IR ones with negative tension located at the ends of an extra-dimensional interval and a UV-brane with positive tension placed in the middle of the interval -- IR-UV-IR model. The background solutions for this geometric setup are found without and with taking into account the backreaction of the matter fields. A 5D $SU(2)$ Higgs doublet is employed as the Goldberger-Wise stabilizing field in this geometry and solutions of the 5D coupled scalar-gravity equations are found by using the superpotential method. Within this setup we investigate the low-energy (zero-mode) effective theory for the bulk Standard Model (SM) bosonic sector. The $\mathbb{Z}_2$-even zero-modes correspond to known standard degrees of freedom, whereas the $\mathbb{Z}_2$-odd zero modes might serve as a dark sector. The effective low-energy scalar sector contains a scalar which mimics the SM Higgs boson and a second stable scalar particle (dark-Higgs) is a dark matter candidate; the latter is a component of the zero-mode of the $\mathbb{Z}_2$-odd Higgs doublet. The model that results from the $\mathbb{Z}_2$-symmetric background geometry resembles the Inert Two Higgs Doublet Model. The effective theory turns out to have an extra residual $SU(2)\times U(1)$ global symmetry that is reminiscent of an underlying 5D gauge transformation for the odd degrees of freedom. At tree level the SM Higgs and the dark-Higgs have the same mass; however, when leading radiative corrections are taken into account the dark-Higgs turns out to be heavier than the SM Higgs. Implications for dark matter are discussed; it is found that the dark-Higgs can provide only a small fraction of the observed dark matter abundance.


title: Background Subtraction via Generalized Fused Lasso Foreground Modeling

abstract: Background Subtraction (BS) is one of the key steps in video analysis. Many background models have been proposed and achieved promising performance on public data sets. However, due to challenges such as illumination change, dynamic background etc. the resulted foreground segmentation often consists of holes as well as background noise. In this regard, we consider generalized fused lasso regularization to quest for intact structured foregrounds. Together with certain assumptions about the background, such as the low-rank assumption or the sparse-composition assumption (depending on whether pure background frames are provided), we formulate BS as a matrix decomposition problem using regularization terms for both the foreground and background matrices. Moreover, under the proposed formulation, the two generally distinctive background assumptions can be solved in a unified manner. The optimization was carried out via applying the augmented Lagrange multiplier (ALM) method in such a way that a fast parametric-flow algorithm is used for updating the foreground matrix. Experimental results on several popular BS data sets demonstrate the advantage of the proposed model compared to state-of-the-arts.


title: First-principles data set of 45,892 isolated and cation-coordinated conformers of 20 proteinogenic amino acids

abstract: We present a structural data set of the 20 proteinogenic amino acids and their amino-methylated and acetylated (capped) dipeptides. Different protonation states of the backbone (uncharged and zwitterionic) were considered for the amino acids as well as varied side chain protonation states. Furthermore, we studied amino acids and dipeptides in complex with divalent cations (Ca2+, Ba2+, Sr2+, Cd2+, Pb2+, and Hg2+). The database covers the conformational hierarchies of 280 systems in a wide relative energy range of up to 4 eV (390 kJ/mol), summing up to an overall of 45,892 stationary points on the respective potential-energy surfaces. All systems were calculated on equal first-principles footing, applying density-functional theory in the generalized gradient approximation corrected for long-range van der Waals interactions. We show good agreement to available experimental data for gas-phase ion affinities. Our curated data can be utilized, for example, for a wide comparison across chemical space of the building blocks of life, for the parametrization of protein force fields, and for the calculation of reference spectra for biophysical applications.


title: Linear metric and temperature fluctuations of a charged plasma in a primordial magnetic field

abstract: We discuss tensor metric perturbations in a magnetic field around the homogeneous Juttner equilibrium of massless particles in an expanding universe. We solve the Liouville equation and derive the energy-momentum tensor up to linear terms in the metric and in the magnetic field.The term linear in the magnetic field is different from zero if the total charge of the primordial plasma is non-zero. We obtain an analytic formula for temperature fluctuations treating the tensor metric perturbations and the magnetic field as independent random variables. Assuming a cutoff on large momenta of the magnetic spectral function we show that the presence of the magnetic field can discriminate only low multipoles in the multipole expansion of temperature fluctuations. In such a case the term linear in the magnetic field can be more important than the quadratic one (corresponding to the fluctuations of the pure magnetic field).


title: Sub-Riemannian mean curvature flow for image processing

abstract: In this paper we reconsider the sub-Riemannian cortical model of image completion introduced by Sarti and Citti in 2003. This model combines two mechanisms, the sub-Riemannian diffusion and the concentration, giving rise to a diffusion driven motion by curvature. In this paper we give a formal proof of the existence of viscosity solutions of the sub-Riemannian motion by curvature. Results of completion and enhancement on a number of natural images are shown and compared with other models.


title: Checking Interaction-Based Declassification Policies for Android Using Symbolic Execution

abstract: Mobile apps can access a wide variety of secure information, such as contacts and location. However, current mobile platforms include only coarse access control mechanisms to protect such data. In this paper, we introduce interaction-based declassification policies, in which the user's interactions with the app constrain the release of sensitive information. Our policies are defined extensionally, so as to be independent of the app's implementation, based on sequences of security-relevant events that occur in app runs. Policies use LTL formulae to precisely specify which secret inputs, read at which times, may be released. We formalize a semantic security condition, interaction-based noninterference, to define our policies precisely. Finally, we describe a prototype tool that uses symbolic execution to check interaction-based declassification policies for Android, and we show that it enforces policies correctly on a set of apps.


title: Measuring the Graph Concordance of Locally Dependent Observations

abstract: This paper introduces a simple measure of a concordance pattern among observed outcomes along a network, i.e., the pattern in which adjacent outcomes tend to be more strongly correlated than non-adjacent outcomes. The graph concordance measure can be generally used to quantify the empirical relevance of a network in explaining cross-sectional dependence of the outcomes, and as shown in the paper, can also be used to quantify the extent of homophily under certain conditions. When one observes a single large network, it is nontrivial to make inference about the concordance pattern. Assuming dependency graph, this paper develops a permutation-based confidence interval for the graph concordance measure. The confidence interval is valid in finite samples when the outcomes are exchangeable, and under the dependency graph assumption together with other regularity conditions, is shown to exhibit asymptotic validity. Monte Carlo simulation results show that the validity of the permutation method is more robust to various graph configurations than the asymptotic method.


title: Detecting Sponsored Recommendations

abstract: With a vast number of items, web-pages, and news to choose from, online services and the customers both benefit tremendously from personalized recommender systems. Such systems however provide great opportunities for targeted advertisements, by displaying ads alongside genuine recommendations. We consider a biased recommendation system where such ads are displayed without any tags (disguised as genuine recommendations), rendering them indistinguishable to a single user. We ask whether it is possible for a small subset of collaborating users to detect such a bias. We propose an algorithm that can detect such a bias through statistical analysis on the collaborating users' feedback. The algorithm requires only binary information indicating whether a user was satisfied with each of the recommended item or not. This makes the algorithm widely appealing to real world issues such as identification of search engine bias and pharmaceutical lobbying. We prove that the proposed algorithm detects the bias with high probability for a broad class of recommendation systems when sufficient number of users provide feedback on sufficient number of recommendations. We provide extensive simulations with real data sets and practical recommender systems, which confirm the trade offs in the theoretical guarantees.


title: The Next Generation Virgo Cluster Survey. XII. Stellar Populations and Kinematics of Compact, Low-Mass Early-Type Galaxies from Gemini GMOS-IFU Spectroscopy

abstract: We present Gemini GMOS-IFU data of eight compact low-mass early-type galaxies (ETGs) in the Virgo cluster. We analyse their stellar kinematics, stellar population, and present two-dimensional maps of these properties covering the central 5"x 7" region. We find a large variety of kinematics: from non- to highly-rotating objects, often associated with underlying disky isophotes revealed by deep images from the Next Generation Virgo Cluster Survey. In half of our objects, we find a centrally-concentrated younger and more metal-rich stellar population. We analyze the specific stellar angular momentum through the lambdaR parameter and find six fast-rotators and two slow-rotators, one having a thin counter-rotating disk. We compare the local galaxy density and stellar populations of our objects with those of 39 more extended low-mass Virgo ETGs from the SMAKCED survey and 260 massive ($M>10^{10}$\Msun) ETGs from the A3D sample. The compact low-mass ETGs in our sample are located in high density regions, often close to a massive galaxy and have, on average, older and more metal-rich stellar populations than less compact low-mass galaxies. We find that the stellar population parameters follow lines of constant velocity dispersion in the mass-size plane, smoothly extending the comparable trends found for massive ETGs. Our study supports a scenario where low-mass compact ETGs have experienced long-lived interactions with their environment, including ram-pressure stripping and gravitational tidal forces, that may be responsible for their compact nature.


title: On the Requirements of New Software Development

abstract: Changes, they use to say, are the only constant in life. Everything changes rapidly around us, and more and more key to survival is the ability to rapidly adapt to changes. This consideration applies to many aspects of our lives. Strangely enough, this nearly self-evident truth is not always considered by software engineers with the seriousness that it calls for: The assumptions we draw for our systems often do not take into due account that e.g., the run-time environments, the operational conditions, or the available resources will vary. Software is especially vulnerable to this threat, and with today's software-dominated systems controlling crucial services in nuclear plants, airborne equipments, health care systems and so forth, it becomes clear how this situation may potentially lead to catastrophes. This paper discusses this problem and defines some of the requirements towards its effective solution, which we call "New Software Development" as a software equivalent of the well-known concept of New Product Development. The paper also introduces and discusses a practical example of a software tool designed taking those requirements into account --- an adaptive data integrity provision in which the degree of redundancy is not fixed once and for all at design time, but rather it changes dynamically with respect to the disturbances experienced during the run time.


title: On hyperbolic equations and systems with non-regular time dependent coefficients

abstract: In this paper we study higher order weakly hyperbolic equations with time dependent non-regular coefficients. The non-regularity here means less than H\"older, namely bounded coefficients. As for second order equations in \cite{GR:14} we prove that such equations admit a `very weak solution' adapted to the type of solutions that exist for regular coefficients. The main idea in the construction of a very weak solution is the regularisation of the coefficients via convolution with a mollifier and a qualitative analysis of the corresponding family of classical solutions depending on the regularising parameter. Classical solutions are recovered as limit of very weak solutions. Finally, by using a reduction to block Sylvester form we conclude that any first order hyperbolic system with non-regular coefficients is solvable in the very weak sense.


title: On the composition of finite rotations in $\mathbb{E}^4$

abstract: We achieve compositions rules for the geometric parameters of the composed rotations, which is in a certain sense analogous to the well known Rodrigues formula. We also obtain a necessary and sufficient condition for a composition of two simple rotations in $\mathbb{E}^4$ to be a simple rotation.


title: A simple and robust confidence interval for causal effects with possibly invalid instruments

abstract: Instrumental variables have been widely used to estimate the causal effect of a treatment on an outcome. Existing confidence intervals for causal effects based on instrumental variables assume that all of the putative instrumental variables are valid; a valid instrumental variable is a variable that affects the outcome only by affecting the treatment and is not related to unmeasured confounders. However, in practice, some of the putative instrumental variables are likely to be invalid. This paper presents a simple and general approach to construct a confidence interval that is robust to possibly invalid instruments. The robust confidence interval has theoretical guarantees on having the correct coverage and can also be used to assess the sensitivity of inference when instrumental variables assumptions are violated. The paper also shows that the robust confidence interval outperforms traditional confidence intervals popular in instrumental variables literature when invalid instruments are present. The new approach is applied to a developmental economics study of the causal effect of income on food expenditures.


title: Some New Directions for ACP Research

abstract: This paper lists some new directions for research related to the Algebra of Communicating Processes (ACP). Most of these directions have been inspired by work on SubScript, an ACP based extension to the programming language Scala. SubScript applies several new ideas that build on ACP, but currently these lack formal treatment. Some of these new ideas are rather fundamental. E.g. it appears that the theory of ACP may well apply to structures of any kind of items, rather than to just processes. The aim of this list is to raise awareness of the research community about these new ideas; this could help both the research area and the programming language SubScript.


title: Two fast X-ray transients in archival Chandra data

abstract: We present the discovery of two new X-ray transients in archival Chandra data. The first transient, XRT 110103, occurred in January 2011 and shows a sharp rise of at least three orders of magnitude in count rate in less than 10 s, a flat peak for about 20 s and decays by two orders of magnitude in the next 60 s. We find no optical or infrared counterpart to this event in preexisting survey data or in an observation taken by the SIRIUS instrument at the Infrared Survey Facility 2.1 yr after the transient, providing limiting magnitudes of J>18.1, H>17.6 and Ks>16.3. This event shows similarities to the transient previously reported in Jonker et al. which was interpreted as the possible tidal disruption of a white dwarf by an intermediate mass black hole. We discuss the possibility that these transients originate from the same type of event. If we assume these events are related a rough estimate of the rates gives 1.4*10^5 per year over the whole sky with a peak 0.3-7 keV X-ray flux greater than 2*10^-10 erg cm^-2 s^-1 . The second transient, XRT 120830, occurred in August 2012 and shows a rise of at least three orders of magnitude in count rate and a subsequent decay of around one order of magnitude all within 10 s, followed by a slower quasi-exponential decay over the remaining 30 ks of the observation. We detect a likely infrared counterpart with magnitudes J=16.70+/-0.06, H=15.92+/-0.04 and Ks=15.37+/-0.06 which shows an average proper motion of 74+/-19 milliarcsec per year compared to archival 2MASS observations. The JHKs magnitudes, proper motion and X-ray flux of XRT 120830 are consistent with a bright flare from a nearby late M or early L dwarf.


title: The Douglas-Rachford algorithm for two (not necessarily intersecting) affine subspaces

abstract: The Douglas--Rachford algorithm is a classical and very successful splitting method for finding the zeros of the sums of monotone operators. When the underlying operators are normal cone operators, the algorithm solves a convex feasibility problem. In this paper, we provide a detailed study of the Douglas--Rachford iterates and the corresponding {shadow sequence} when the sets are affine subspaces that do not necessarily intersect. We prove strong convergence of the shadows to the nearest generalized solution. Our results extend recent work from the consistent to the inconsistent case. Various examples are provided to illustrates the results.


title: The Distributions of Hilbert Space Frame Vectors and Frame Coefficients

abstract: The most fundamental notion for Hilbert space frames is the sequence of frame coefficients for a vector x in the space. Yet, we know little about the distribution of these coefficient sequences. In this paper, we make the first detailed study of the distribution of the frame coefficients for vectors in a Hilbert space, as well as the related notion of the square sums of the distances of vectors in the space from the frame vectors. We will give some surprisingly exact calculations for special frames such as unit norm tight frames and equiangular frames. This is a study which should have been done 20 years ago.


title: Bulk Duals for Generic Static, Scale-Invariant Holographic CFT States

abstract: Near horizon geometries have been widely studied, and have found many applications. Certain static, near horizon geometries are now understood to be bulk duals to CFTs with static scale-invariant sources under the AdS/CFT correspondence. However, static near-horizon geometries aren't just scale-invariant, they have extra `enhanced' symmetry. This means that they can only be the bulk duals for a special class of static, scale-invariant sources that share this enhanced symmetry. The purpose of this paper is to consider bulk duals for more generic static, scale-invariant sources, without this extra symmetry. These solutions are quite different to near-horizon geometries. In place of the extremal horizon they have a null singularity. We find specific examples of such bulk geometries numerically for the cases of pure gravity, and for an abelian gauge field.


title: Comparison of simple mass estimators for slowly rotating elliptical galaxies

abstract: We compare the performance of mass estimators for elliptical galaxies that rely on the directly observable surface brightness and velocity dispersion profiles, without invoking computationally expensive detailed modeling. These methods recover the mass at a specific radius where the mass estimate is expected to be least sensitive to the anisotropy of stellar orbits. One method (Wolf et al. 2010) uses the total luminosity-weighted velocity dispersion and evaluates the mass at a 3D half-light radius $r_{1/2}$, i.e., it depends on the GLOBAL galaxy properties. Another approach (Churazov et al. 2010) estimates the mass from the velocity dispersion at a radius $R_2$ where the surface brightness declines as $R^{-2}$, i.e., it depends on the LOCAL properties. We evaluate the accuracy of the two methods for analytical models, simulated galaxies and real elliptical galaxies that have already been modeled by the Schwarzschild's orbit-superposition technique. Both estimators recover an almost unbiased circular speed estimate with a modest RMS scatter ($\lesssim 10 \%$). Tests on analytical models and simulated galaxies indicate that the local estimator has a smaller RMS scatter than the global one. We show by examination of simulated galaxies that the projected velocity dispersion at $R_2$ could serve as a good proxy for the virial galaxy mass. For simulated galaxies the total halo mass scales with $\sigma_p(R_2)$ as $M_{vir} \left[M_{\odot}h^{-1}\right] \approx 6\cdot 10^{12} \left( \frac{\sigma_p(R_2)}{200\, \rm km\, s^{-1}} \right)^{4}$ with RMS scatter $\approx 40 \%$.


title: An Algorithm for Global Maximization of Secrecy Rates in Gaussian MIMO Wiretap Channels

abstract: Optimal signaling for secrecy rate maximization in Gaussian MIMO wiretap channels is considered. While this channel has attracted a significant attention recently and a number of results have been obtained, including the proof of the optimality of Gaussian signalling, an optimal transmit covariance matrix is known for some special cases only and the general case remains an open problem. An iterative custom-made algorithm to find a globally-optimal transmit covariance matrix in the general case is developed in this paper, with guaranteed convergence to a \textit{global} optimum. While the original optimization problem is not convex and hence difficult to solve, its minimax reformulation can be solved via the convex optimization tools, which is exploited here. The proposed algorithm is based on the barrier method extended to deal with a minimax problem at hand. Its convergence to a global optimum is proved for the general case (degraded or not) and a bound for the optimality gap is given for each step of the barrier method. The performance of the algorithm is demonstrated via numerical examples. In particular, 20 to 40 Newton steps are already sufficient to solve the sufficient optimality conditions with very high precision (up to the machine precision level), even for large systems. Even fewer steps are required if the secrecy capacity is the only quantity of interest. The algorithm can be significantly simplified for the degraded channel case and can also be adopted to include the per-antenna power constraints (instead or in addition to the total power constraint). It also solves the dual problem of minimizing the total power subject to the secrecy rate constraint.


title: The INTEGRAL/SPI view of A0535+26 during the giant outburst of February 2011

abstract: A0535+26 is a slowly rotating pulsar accreting from the wind of a massive Be star, and that exhibits two cyclotron absorption lines in its X-ray spectrum, at about 45 and 100 keV, respectively. Unlike similar sources, no significant variations of the energy of its cyclotron lines with flux were observed to date. The bright outburst of February 2011 thus offers a unique occasion to probe this peculiar behavior at flux levels not yet observed with present-day instruments. Here we report on the spectral and timing analysis of the data from the spectrometer SPI on-board INTEGRAL collected during the outburst. At the peak of the outburst the estimated luminosity is ~4.9x10^37 erg/s. The fundamental cyclotron feature is detected at all flux levels, and its centroid energy is positively-correlated with the flux of the source, confirming that A0535+26 is accreting at a sub-critical regime. The correlation seems to fall off at ~10^37 erg/s, suggesting the transition from a Coulomb-stopping regime to a gas-mediated shock regime. From the timing analysis we found that the pulsar was spinning up during most of the outburst, and that the spin-up rate correlates with the flux of the source, albeit the correlation is steeper than the one expected from the standard disk accretion theory. Finally, we show that the pulse profile of the source changes dramatically as the flux increases. At high luminosity the profile is highly asymmetric, implying an asymmetry in the geometry of the accretion flow.


title: Cavity optomechanical transduction sensing of single molecules

abstract: We report narrow linewidth optomechanical oscillation of a silica microsphere immersed in a buffer solution. Through a novel optomechanical transduction sensing approach, single 10-nm-radius silica beads and Bovine serum albumin (BSA) protein molecules with a molecular weight of 66 kDalton were detected. This approach predicts the detection of 3.9 kDalton single molecules at a signal-to-noise ration above unity.


title: Robustness of power systems under a democratic fiber bundle-like model

abstract: We consider a power system with $N$ transmission lines whose initial loads (i.e., power flows) $L_1, \ldots, L_N$ are independent and identically distributed with $P_L(x)$. The capacity $C_i$ defines the maximum flow allowed on line $i$, and is assumed to be given by $C_i=(1+\alpha)L_i$, with $\alpha>0$. We study the robustness of this power system against random attacks (or, failures) that target a $p$-{\em fraction} of the lines, under a democratic fiber bundle-like model. Namely, when a line fails, the load it was carrying is redistributed equally among the remaining lines. Our contributions are as follows: i) we show analytically that the final breakdown of the system always takes place through a first-order transition at the critical attack size $p^{\star}=1-\frac{E[L]}{\max\{P(L>x)(\alpha x + E[L ~|~ L>x])\}}~~~$; ii) we derive conditions on the distribution $P_L(x)$ for which the first order break down of the system occurs abruptly without any preceding diverging rate of failure; iii) we provide a detailed analysis of the robustness of the system under three specific load distributions: Uniform, Pareto, and Weibull, showing that with the minimum load $L_{\textrm{min}}$ and mean load $E[L]$ fixed, Pareto distribution is the worst (in terms of robustness) among the three, whereas Weibull distribution is the best with shape parameter selected relatively large; iv) we provide numerical results that confirm our mean-field analysis; and v) we show that $p^{\star}$ is maximized when the load distribution is a Dirac delta function centered at $E[L]$, i.e., when all lines carry the same load; we also show that optimal $p^{\star}$ equals $\frac{\alpha}{\alpha+1}$. This last finding is particularly surprising given that heterogeneity is known to lead to high robustness against random failures in many other systems.


title: Logarithmic Correlations in Quantum Hall Plateau Transitions

abstract: The critical behavior of quantum Hall transitions in two-dimensional disordered electronic systems can be described by a class of complicated, non-unitary conformal field theories with logarithmic correlations. The nature and the physical origin of these logarithmic correlation functions remain however mysterious. Using the replica trick and the underlying symmetries of these quantum critical points, we show here how to construct non-perturbatively disorder-averaged observables in terms of Green's functions that scale logarithmically at criticality. In the case of the spin quantum Hall transition, which may occur in disordered superconductors with spin-rotation symmetry and broken time reversal invariance, we argue that our results are compatible with an alternative approach based on supersymmetry. The generalization to the Integer quantum Hall plateau transition is also discussed.


title: Optimal Training for Non-Feedback Adaptive PSAM over Time-Varying Rayleigh Fading Channels

abstract: Time-varying fast fading channels present a major challenge in the design of wireless communication systems. Pilot Symbol Assisted Modulation (PSAM) has been introduced to mitigate the effects of fading and allow coherent demodulation. Our work studies the performance of \emph{non-feedback} adaptive PSAM scheme over time-varying Rayleigh fading channels. A modular method is introduced for computing the rates in an efficient manner. Moreover, four transmission policies are analyzed and we show how optimal training in terms of duration and power allocation varies with the channel conditions and from one transmission policy to another. The performance of these schemes is measured in terms of achievable rates using binary signaling. We formally show that, for a causal estimation, placing all the power on the last pilot symbol is expected to be optimal. Furthermore, the autocorrelation of the fading process is based either on a stationary first order Gauss-Markov modeling of the process or on Jakes' model when higher orders of correlation are analyzed.


title: Safety enhancement through situation-aware user interfaces

abstract: Due to their privileged position halfway between the physical and the cyber universes, user interfaces may play an important role in preventing, tolerating, and learning from scenarios potentially affecting mission safety and the user's quality of experience. This vision is embodied here in the main ideas and a proof-of-concepts implementation of user interfaces that combine dynamic profiling with context- and situation-awareness and autonomic software adaptation.


title: Abelian Tensors

abstract: We analyze tensors in the tensor product of three m-dimensional vector spaces satisfying Strassen's equations for border rank m. Results include: two purely geometric characterizations of the Coppersmith-Winograd tensor, a reduction to the study of symmetric tensors under a mild genericity hypothesis, and numerous additional equations and examples. This study is closely connected to the study of the variety of m-dimensional abelian subspaces of the space of endomorphisms of an m-dimensional vector space, and the subvariety consisting of the Zariski closure of the variety of maximal tori, called the variety of reductions.


title: Switching-GAS Copula Models With Application to Systemic Risk

abstract: Recent financial disasters have emphasised the need to accurately predict extreme financial losses and their consequences for the institutions belonging to a given financial market. The ability of econometric models to predict extreme events strongly relies on their flexibility to account for the highly nonlinear and asymmetric dependence observed in financial returns. We develop a new class of flexible Copula models where the evolution of the dependence parameters follow a Markov-Switching Generalised Autoregressive Score (SGASC) dynamics. Maximum Likelihood estimation is consistently performed using the Inference Functions for Margins (IFM) approach and a version of the Expectation-Maximisation (EM) algorithm specifically tailored to this class of models. The SGASC models are then used to estimate the Conditional Value-at-Risk (CoVaR), which is defined as the VaR of a given asset conditional on another asset (or portfolio) being in financial distress, and the Conditional Expected Shortfall (CoES). Our empirical investigation shows that the proposed SGASC models are able to explain and predict the systemic risk contribution of several European countries. Moreover, we also find that the SGASC models outperform competitors using several CoVaR backtesting procedures.


title: Non-modal linear stability analysis of miscible viscous fingering in porous media

abstract: The non-modal linear stability of miscible viscous fingering in a two dimensional homogeneous porous medium has been investigated. The linearized perturbed equations for Darcy's law coupled with a convection-diffusion equation is discretized using finite difference method. The resultant initial value problem is solved by fourth order Runge-Kutta method, followed by a singular value decomposition of the propagator matrix. Particular attention is given to the transient behavior rather than the long-time behavior of eigenmodes predicted by the traditional modal analysis. The transient behaviors of the response to external excitations and the response to initial conditions are studied by examining the $\epsilon-$pseudospectra structures and the largest energy growth function. With the help of non-modal stability analysis we demonstrate that at early times the displacement flow is dominated by diffusion and the perturbations decay. At later times, when convection dominates diffusion, perturbations grow. Furthermore, we show that the dominant perturbation that experiences the maximum amplification within the linear regime lead to the transient growth. These two important features were previously unattainable in the existing linear stability methods for miscible viscous fingering. To explore the relevance of the optimal perturbation obtained from non-modal analysis, we performed direct numerical simulations using a highly accurate pseudo-spectral method. Further, a comparison of the present stability analysis with existing modal and initial value approach is also presented. It is shown that the non-modal stability results are in better agreement, than the other existing stability analyses, with those obtained from direct numerical simulations.


title: Non-linear Ion-wake Excitation by Ultra-relativistic Electron Wakefields

abstract: The excitation of a non-linear ion-wake by a train of ultra-relativistic electron bubble wake [1]-[8] is modeled. The ion-wake is shown to be a driven non-linear ion-acoustic wave in the form of a cylindrical ion-soliton [11]-[14]. The phases of the oscillating radial electric fields of the slowly-propagating [21] electron bubble is asymmetric in time and excites time-averaged inertial ion motion radially. The electron compression in the back of the bubble sucks-in the ions and the space-charge within the bubble cavity expels them, driving a cylindrical ion-soliton structure with on-axis and bubble-edge density-spikes [5][6]. Once formed, the channel-edge density-spike is driven radially outwards by the thermal pressure of the wake energy [12]. Its channel-like structure due to the flat-residue left behind by the propagating ion-soliton, is independent of the energy-source driving the bubble [3][4] electron wake. We explore the use of the partially- filled channel formed by the cylindrical ion-soliton for a novel regime of positron acceleration [18]. OSIRIS PIC [25] simulations are used to study the ion-wake soliton structure, its driven propagation and its use for positron acceleration.

categories: physics.plasm-ph astro-ph.HE nlin.PS physics.acc-ph


title: Margination regimes and drainage transition in confined multicomponent suspensions

abstract: A mechanistic theory is developed to describe segregation in confined multicomponent suspensions such as blood. It incorporates the two key phenomena arising in these systems at low Reynolds number: hydrodynamic pair collisions and wall-induced migration. In simple shear flow, several regimes of segregation arise, depending on the value of a "margination parameter" M. Most importantly, there is a critical value of M below which a sharp "drainage transition" occurs: one component is completely depleted from the bulk flow to the vicinity of the walls. Direct simulations also exhibit this transition as the size or flexibility ratio of the components changes.


title: Hyperbolic Metamaterial Feasible for Fabrication with Direct Laser Writing Processes

abstract: Stimulated emission depletion microscopy inspired direct laser writing (STED-DLW) processes can offer diffraction-unlimited fabrication of 3D-structures, not possible with traditional electron-beam or optical lithography. We propose a hyperbolic metamaterial for fabrication with STED-DLW. First, we design meandering wire structures with three different magnetic dipoles which can be excited under different incidences of light. Then, based on effective parameters corresponding to normal incidence and lateral incidence, we find that the hyperbolic dispersion relation for five-layer structure appears between 15THz to 20THz. Finally, we investigate the influence of imaginary parts of the effective parameters on the metamaterial dispersion. The proposed metamaterial structure has also potential for three-dimensionally isotropic permeability despite geometric anisotropy.


title: Amplify-and-Forward Relaying in Two-Hop Diffusion-Based Molecular Communication Networks

abstract: This paper studies a three-node network in which an intermediate nano-transceiver, acting as a relay, is placed between a nano-transmitter and a nano-receiver to improve the range of diffusion-based molecular communication. Motivated by the relaying protocols used in traditional wireless communication systems, we study amplify-and-forward (AF) relaying with fixed and variable amplification factor for use in molecular communication systems. To this end, we derive a closed-form expression for the expected end-to-end error probability. Furthermore, we derive a closed-form expression for the optimal amplification factor at the relay node for minimization of an approximation of the expected error probability of the network. Our analytical and simulation results show the potential of AF relaying to improve the overall performance of nano-networks.


title: Acoustic phonon dynamics in thin-films of the topological insulator Bi2Se3

abstract: Transient reflectivity traces measured for nanometer-sized films of the topological insulator Bi2Se3 revealed GHz-range oscillations driven within the relaxation of hot carriers photoexcited with ultrashort laser pulses of 1.51 eV photon energy. These oscillations have been suggested to result from acoustic phonon dynamics, including coherent longitudinal acoustic phonons in the form of standing acoustic waves. An increase of oscillation frequency from ~35 to ~70 GHz with decreasing film thickness from 40 to 15 nm was attributed to the interplay between two different regimes employing traveling-acoustic-waves for films thicker than 40 nm and the film bulk acoustic wave resonator (FBAWR) modes for films thinner than 40 nm. The amplitude of oscillations decays rapidly for films below 15 nm thick when the indirect intersurface coupling in Bi2Se3 films switches the FBAWR regime to that of the Lamb wave excitation. The frequency range of coherent longitudinal acoustic phonons is in good agreement with elastic properties of Bi2Se3.


title: Cost of postselection in decision theory

abstract: Postselection is the process of discarding outcomes from statistical trials that are not the event one desires. Postselection can be useful in many applications where the cost of getting the wrong event is implicitly high. However, unless this cost is specified exactly, one might conclude that discarding all data is optimal. Here we analyze the optimal decision rules and quantum measurements in a decision theoretic setting where a prespecified cost is assigned to discarding data. Our scheme interpolates between unambiguous state discrimination (when the cost of postselection is zero) and a minimum error measurement (when the cost of postselection is maximal). We also relate our formulation to previous approaches which focus on minimizing the probability of indecision.


title: Jet propagation within a Linearized Boltzmann Transport Model

abstract: A Linear Boltzmann Transport (LBT) model has been developed for the study of jet propagation inside a quark-gluon plasma. Both leading and thermal recoiled partons are transported according to the Boltzmann equations to account for jet-induced medium excitations. In this talk, we present our study within the LBT model in which we implement the complete set of elastic parton scattering processes. We investigate elastic parton energy loss and their energy and length dependence. We further investigate elastic energy loss and transverse shape of reconstructed jets. Contributions from the recoiled thermal partons are found to have significant influences on the jet energy loss and transverse profile.


title: Finiteness of Teichm\"uller curves in non-arithmetic rank 1 orbit closures

abstract: We show that in any non-arithmetic rank 1 orbit closure of translation surfaces, there are only finitely many Teichm\"uller curves. We also show that in any non-arithmetic rank 1 orbit closure, any completely parabolic surface is Veech.


title: Lambda = 3 mm line survey of nearby active galaxies

abstract: We used the IRAM 30m telescope to observe the frequency range [86-116]GHz towards the central regions of the starburst galaxies M83, M82, and NGC253, the AGNs M51, NGC1068, and NGC7469, and the ULIRGs Arp220 and Mrk231. Assuming LTE conditions, we calculated the column densities of 27 molecules and 10 isotopologues. Among others, we report the first tentative detections of CH3CHO, HNCO, and NS in M82 and, for the first time in the extragalactic medium, HC5N in NGC253. Halpha recombination lines were only found in M82 and NGC253. Vibrationally excited lines of HC3N were only detected in Arp220. CH3CCH emission is only seen in the starburst-dominated galaxies. By comparison of the fractional abundances among the galaxies, we looked for the molecules that are best suited to characterise the chemistry of starbursts, AGNs and ULIRGs, as well as the differences among galaxies within the same group.


title: KERMAN: A Key Establishment Algorithm based on Harvesting Randomness in MANETs

abstract: Establishing secret common randomness between two or multiple devices in a network resides at the root of communication security. The problem is traditionally decomposed into a randomness generation stage (randomness purity is subject to employing often costly true random number generators) and a key-agreement information exchange stage, which can rely on public-key infrastructure or on key wrapping. In this paper, we propose KERMAN, an alternative key establishment algorithm for ad-hoc networks which works by harvesting randomness directly from the network routing metadata, thus achieving both pure randomness generation and (implicitly) secret-key agreement. Our algorithm relies on the route discovery phase of an ad-hoc network employing the Dynamic Source Routing protocol, is lightweight, and requires minimal communication overhead.


title: High Spectral and Spatial Resolution Observations of the PDR Emission in the NGC2023 Reflection Nebula with SOFIA and APEX

abstract: We have mapped the NGC 2023 reflection nebula in [CII] and CO(11--10) with the heterodyne receiver GREAT on SOFIA and obtained slightly smaller maps in 13CO(3--2), CO(3--2), CO(4--3), CO(6--5), and CO(7--6) with APEX in Chile. We use these data to probe the morphology, kinematics, and physical conditions of the C II region, which is ionized by FUV radiation from the B2 star HD37903. The [CII] emission traces an ellipsoidal shell-like region at a position angle of ~ -50 deg, and is surrounded by a hot molecular shell. In the southeast, where the C II region expands into a dense, clumpy molecular cloud ridge, we see narrow and strong line emission from high-J CO lines, which comes from a thin, hot molecular shell surrounding the [CII] emission. The [CII] lines are broader and show photo evaporating gas flowing into the C II region. Based on the strength of the [13CII] F=2--1 line, the [CII] line appears to be somewhat optically thick over most of the nebula with an optical depth of a few. We model the physical conditions of the surrounding molecular cloud and the PDR emission using both RADEX and simple PDR models. The temperature of the CO emitting PDR shell is ~ 90 -- 120 K, with densities of 10^5 -- 10^6 cm^-3, as deduced from RADEX modeling. Our PDR modeling indicates that the PDR layer where [CII] emission dominates has somewhat lower densities, 10^4 to a few times 10^5 cm^-3


title: A sparse reformulation of the Green's function formalism allows efficient simulations of partial differential equations on tree graphs

abstract: We prove that when a class of partial differential equations, generalized from the cable equation, is defined on tree graphs, and when the inputs are restricted to a spatially discrete, well chosen set of points, the Green's function (GF) formalism can be rewritten to scale as $O(n)$ with the number $n$ of input locations, contrary to the previously reported $O(n^2)$ scaling. We show that the linear scaling can be combined with an expansion of the remaining kernels as sums of exponentials, to allow efficient simulations of equations from the aforementioned class. We furthermore validate this simulation paradigm on models of nerve cells and explore its relation with more traditional finite difference approaches. Situations in which a gain in computational performance is expected, are discussed.


title: Android Permissions Remystified: A Field Study on Contextual Integrity

abstract: Due to the amount of data that smartphone applications can potentially access, platforms enforce permission systems that allow users to regulate how applications access protected resources. If users are asked to make security decisions too frequently and in benign situations, they may become habituated and approve all future requests without regard for the consequences. If they are asked to make too few security decisions, they may become concerned that the platform is revealing too much sensitive information. To explore this tradeoff, we instrumented the Android platform to collect data regarding how often and under what circumstances smartphone applications are accessing protected resources regulated by permissions. We performed a 36-person field study to explore the notion of "contextual integrity," that is, how often are applications accessing protected resources when users are not expecting it? Based on our collection of 27 million data points and exit interviews with participants, we examine the situations in which users would like the ability to deny applications access to protected resources. We found out that at least 80% of our participants would have preferred to prevent at least one permission request, and overall, they thought that over a third of requests were invasive and desired a mechanism to block them.


title: Minimal helix submanifolds and Minimal Riemannian foliations

abstract: We investigate minimal helix submanifolds of any dimension and codimension immersed in Euclidean space. Our main result proves that a ruled minimal helix submanifold is a cylinder. As an application we classify complex helix submanifolds of $\mathbb{C}^n$: They are extrinsic products with a complex line as a factor. The key tool is Corollary 1.3 which allows us to classify Riemannian foliations of open subsets of the Euclidean space with minimal leaves. Finally, we consider the case of a helix hypersurface with constant mean curvature and prove that it is either a cylinder or an open part of a hyperplane.


title: Galerkin v. least-squares Petrov--Galerkin projection in nonlinear model reduction

abstract: Least-squares Petrov--Galerkin (LSPG) model-reduction techniques such as the Gauss--Newton with Approximated Tensors (GNAT) method have shown promise, as they have generated stable, accurate solutions for large-scale turbulent, compressible flow problems where standard Galerkin techniques have failed. However, there has been limited comparative analysis of the two approaches. This is due in part to difficulties arising from the fact that Galerkin techniques perform optimal projection associated with residual minimization at the time-continuous level, while LSPG techniques do so at the time-discrete level. This work provides a detailed theoretical and computational comparison of the two techniques for two common classes of time integrators: linear multistep schemes and Runge--Kutta schemes. We present a number of new findings, including conditions under which the LSPG ROM has a time-continuous representation, conditions under which the two techniques are equivalent, and time-discrete error bounds for the two approaches. Perhaps most surprisingly, we demonstrate both theoretically and computationally that decreasing the time step does not necessarily decrease the error for the LSPG ROM; instead, the time step should be `matched' to the spectral content of the reduced basis. In numerical experiments carried out on a turbulent compressible-flow problem with over one million unknowns, we show that increasing the time step to an intermediate value decreases both the error and the simulation time of the LSPG reduced-order model by an order of magnitude.


title: Responding to complexity in socio-economic systems: How to build a smart and resilient society?

abstract: The world is changing at an ever-increasing pace. And it has changed in a much more fundamental way than one would think, primarily because it has become more connected and interdependent than in our entire history. Every new product, every new invention can be combined with those that existed before, thereby creating an explosion of complexity: structural complexity, dynamic complexity, functional complexity, and algorithmic complexity. How to respond to this challenge? And what are the costs?


title: Societal, Economic, Ethical and Legal Challenges of the Digital Revolution: From Big Data to Deep Learning, Artificial Intelligence, and Manipulative Technologies

abstract: In the wake of the on-going digital revolution, we will see a dramatic transformation of our economy and most of our societal institutions. While the benefits of this transformation can be massive, there are also tremendous risks to our society. After the automation of many production processes and the creation of self-driving vehicles, the automation of society is next. This is moving us to a tipping point and to a crossroads: we must decide between a society in which the actions are determined in a top-down way and then implemented by coercion or manipulative technologies (such as personalized ads and nudging) or a society, in which decisions are taken in a free and participatory way and mutually coordinated. Modern information and communication systems (ICT) enable both, but the latter has economic and strategic benefits. The fundaments of human dignity, autonomous decision-making, and democracies are shaking, but I believe that they need to be vigorously defended, as they are not only core principles of livable societies, but also the basis of greater efficiency and success.


title: The Physical Nature and Orbital Behavior of the Eclipsing System DK Cygni

abstract: New CCD photometry is presented for the hot overcontact binary DK Cyg, together with reasonable explanations for the light and period variations. Historical light and velocity curves from 1962 to 2012 were simultaneously analyzed with the Wilson-Devinney (W-D) synthesis code. The brightness disturbances were satisfactorily modeled by applying a magnetic cool spot on the primary star. Based on 261 times of minimum light including 116 new timings and spanning more than 87 yrs, a period study reveals that the orbital period has varied due to a periodic oscillation superposed on an upward parabola. The period and semi-amplitude of the modulation are about 78.1 yrs and 0.0037 d, respectively. This detail is interpreted as a light-travel-time effect due to a circumbinary companion with a minimum mass of $M_3$=0.065 $M_\odot$, within the theoretical limit of $\sim$0.07 M$_\odot$ for a brown dwarf star. The observed period increase at a fractional rate of $+$2.74 $\times $10$^{-10}$ is in excellent agreement with that calculated from our W-D synthesis. Mass transfer from the secondary to the primary component is mainly responsible for the secular period change. We examined the evolutionary status of the DK Cyg system from the absolute dimensions.


title: Measurement of the Atmospheric $\nu_e$ Spectrum with IceCube

abstract: We present a measurement of the atmospheric $\nu_e$ spectrum at energies between 0.1 TeV and 100 TeV using data from the first year of the complete IceCube detector. Atmospheric $\nu_e$ originate mainly from the decays of kaons produced in cosmic-ray air showers. This analysis selects 1078 fully contained events in 332 days of livetime, then identifies those consistent with particle showers. A likelihood analysis with improved event selection extends our previous measurement of the conventional $\nu_e$ fluxes to higher energies. The data constrain the conventional $\nu_e$ flux to be $1.3^{+0.4}_{-0.3}$ times a baseline prediction from a Honda's calculation, including the knee of the cosmic-ray spectrum. A fit to the kaon contribution ($\xi$) to the neutrino flux finds a kaon component that is $\xi =1.3^{+0.5}_{-0.4}$ times the baseline value. The fitted/measured prompt neutrino flux from charmed hadron decays strongly depends on the assumed astrophysical flux and shape. If the astrophysical component follows a power law, the result for the prompt flux is $0.0^{+3.0}_{-0.0}$ times a calculated flux based on the work by Enberg, Reno and Sarcevic.


title: Throughput and Delay Scaling of Content-Centric Ad Hoc and Heterogeneous Wireless Networks

abstract: We study the throughput and delay characteristics of wireless caching networks, where users are mainly interested in retrieving content stored in the network, rather than in maintaining source-destination communication. Nodes are assumed to be uniformly distributed in the network area. Each node has a limited-capacity content store, which it uses to cache contents. We propose an achievable caching and transmission scheme whereby requesters retrieve content from the caching point which is closest in Euclidean distance. We establish the throughput and delay scaling of the achievable scheme, and show that the throughput and delay performance are order-optimal within a class of schemes. We then solve the caching optimization problem, and evaluate the network performance for a Zipf content popularity distribution, letting the number of content types and the network size both go to infinity. Finally, we extend our analysis to heterogeneous wireless networks where, in addition to wireless nodes, there are a number of base stations uniformly distributed at random in the network area. We show that in order to achieve a better performance in a heterogeneous network in the order sense, the number of base stations needs to be greater than the ratio of the number of nodes to the number of content types. Furthermore, we show that the heterogeneous network does not yield performance advantages in the order sense if the Zipf content popularity distribution exponent exceeds 3/2.


title: Differential Chow varieties exist

abstract: Chow varieties are a parameter space for cycles of a given variety of a given codimension and degree. We construct their analog for differential algebraic varieties with differential algebraic subvarieties, answering a question of Gao, Li and Wei. The proof uses the construction of classical algebro-geometric Chow varieties, the theory of characteristic sets of differential varieties, the theory of prolongation spaces, and the theory of differential Chow forms. In the course of the proof several definability results from the theory of algebraically closed fields are required. Elementary proofs of these results are given in an appendix by William Johnson.


title: Invariants of a general branched cover of ${\bf P}^{1}$

abstract: We investigate the resolution of a general branched cover $\alpha : C \to {\bf P}^1$ in its relative canonical embedding $C \subset {\bf P} E$. We conjecture that the syzygy bundles appearing in the resolution are balanced for a general cover, provided that the genus is sufficiently large compared to the degree. We prove this for the Casnati-Ekedahl bundle, or "bundle of quadrics" $F$ - the first bundle appearing in the resolution of the ideal of the relative canonical embedding. Furthermore, we prove the conjecture for all syzygy bundles in the resolution when the genus satisfies $g = 1 \mod d$.


title: Strange Duality of Verlinde spaces for G_2 and F_4

abstract: We prove that the pull back of the canonical theta divisor for E_8-bundles at level one induces a strange duality between Verlinde spaces for G_2 and F_4 at level one on smooth curves of genus g. We also prove a parabolic generalization in terms of conformal blocks and write down identities between conformal blocks divisors in the Picard group of \bar-M_{g,n}.


title: On the number of edges in a graph with no $(k+1)$-connected subgraphs

abstract: Mader proved that for $k\geq 2$ and $n\geq 2k$, every $n$-vertex graph with no $(k+1)$-connected subgraphs has at most $(1+\frac{1}{\sqrt{2}})k(n-k)$ edges. He also conjectured that for $n$ large with respect to $k$, every such graph has at most $\frac{3}{2}\left(k - \frac{1}{3}\right)(n-k)$ edges. Yuster improved Mader's upper bound to $\frac{193}{120}k(n-k)$ for $n\geq\frac{9k}{4}$. In this note, we make the next step towards Mader's Conjecture: we improve Yuster's bound to $\frac{19}{12}k(n-k)$ for $n\geq\frac{5k}{2}$.


title: Stain-Induced Band Inversion and Topological Nontriviality in Antimonene

abstract: Antimonene is a novel two-dimensional (2D) semiconducting material of group V elements proposed in a recent literature [Zhang et al., Angew. Chem. Int. Ed. 54, 1-5 (2015)]. Using first-principles calculations, we demonstrated that the buckled configuration of antimonene enables it sustain large tensile strain up to 20%. Band inversion takes place in the vicinity of the Gamma point as the tensile strain is larger than 14.5%, leading to six tilted Dirac cones in the Brillouin zone. Spin-orbital coupling (SOC) effect opens up a topologically nontrivial bulk band gap at the Dirac points, exhibiting the features of 2D topological insulators characterized by a nonzero Z2 topological invariant. The tunable bulk band gap, 101-560 meV, make the antimonene a promising candidate material for achieving quantum spin Hall effect (QSH) at high temperatures which meet the requirement of future electronic devices with low power consumption.


title: New variables in M5 (NGC 5904) and some identification corrections

abstract: We report twelve variables not previously detected in the globular cluster M5 (NGC 5904); one SX Phe and eleven semi-regular variables (SR). Their identifications, equatorial coordinates, ephemerides, and light curves are given. Furthermore, we have explored the light curves of a group of stars whose variability has not been confirmed and that are marked as probable non- variables in the CVSGC. Finally, we offer detailed identifications for some of the known variables in crowded regions that were misidentified in previous studies. We shall also address the cases of the cataclysmic variable or U Gem type V101 and of the variable blue straggler V159.


title: Lower Bounds on Complexity of Lyapunov Functions for Switched Linear Systems

abstract: We show that for any positive integer $d$, there are families of switched linear systems---in fixed dimension and defined by two matrices only---that are stable under arbitrary switching but do not admit (i) a polynomial Lyapunov function of degree $\leq d$, or (ii) a polytopic Lyapunov function with $\leq d$ facets, or (iii) a piecewise quadratic Lyapunov function with $\leq d$ pieces. This implies that there cannot be an upper bound on the size of the linear and semidefinite programs that search for such stability certificates. Several constructive and non-constructive arguments are presented which connect our problem to known (and rather classical) results in the literature regarding the finiteness conjecture, undecidability, and non-algebraicity of the joint spectral radius. In particular, we show that existence of an extremal piecewise algebraic Lyapunov function implies the finiteness property of the optimal product, generalizing a result of Lagarias and Wang. As a corollary, we prove that the finiteness property holds for sets of matrices with an extremal Lyapunov function belonging to some of the most popular function classes in controls.


title: Attractors of Local Semiflows on Topological Spaces

abstract: In this paper we introduce a notion of an attractor for local semiflows on topological spaces, which in some cases seems to be more suitable than the existing ones in the literature. Based on this notion we develop a basic attractor theory on topological spaces under appropriate separation axioms. First, we discuss fundamental properties of attractors such as maximality and stability and establish some existence results. Then, we give a converse Lyapunov theorem. Finally, the Morse decomposition of attractors is also addressed.


title: Mutiscale Mapper: A Framework for Topological Summarization of Data and Maps

abstract: Summarizing topological information from datasets and maps defined on them is a central theme in topological data analysis. \textsf{Mapper}, a tool for such summarization, takes as input both a possibly high dimensional dataset and a map defined on the data, and produces a summary of the data by using a cover of the codomain of the map. This cover, via a pullback operation to the domain, produces a simplicial complex connecting the data points. The resulting view of the data through a cover of the codomain offers flexibility in analyzing the data. However, it offers only a view at a fixed scale at which the cover is constructed. Inspired by the concept, we explore a notion of a tower of covers which induces a tower of simplicial complexes connected by simplicial maps, which we call {\em multiscale mapper}. We study the resulting structure, its stability, and design practical algorithms to compute its associated persistence diagrams efficiently. Specifically, when the domain is a simplicial complex and the map is a real-valued piecewise-linear function, the algorithm can compute the exact persistence diagram only from the 1-skeleton of the input complex. For general maps, we present a combinatorial version of the algorithm that acts only on \emph{vertex sets} connected by the 1-skeleton graph, and this algorithm approximates the exact persistence diagram thanks to a stability result that we show to hold. We also relate the multiscale mapper with the \v{C}ech complexes arising from a natural pullback pseudometric defined on the input domain.


title: Calculation of a fluctuating entropic force by phase space sampling

abstract: A polymer chain pinned in space exerts a fluctuating force on the pin point in thermal equilibrium. The average of such fluctuating force is well understood from statistical mechanics as an entropic force, but little is known about the underlying force distribution. Here, we introduce two phase space sampling methods that can produce the equilibrium distribution of instantaneous forces exerted by a terminally pinned polymer. In these methods, both the positions and momenta of mass points representing a freely jointed chain are perturbed in accordance with the spatial constraints and the Boltzmann distribution of total energy. The constraint force for each conformation and momentum is calculated using Lagrangian dynamics. Using terminally pinned chains in space and on a surface, we show that the force distribution is highly asymmetric with both tensile and compressive forces. Most importantly, the mean of the distribution, which is equal to the entropic force, is not the most probable force even for long chains. Our work provides insights into the mechanistic origin of entropic forces, and an efficient computational tool for unbiased sampling of the phase space of a constrained system.


title: Free Probability for Pairs of Faces III: 2-Variables Bi-free Partial S- and T-Transforms

abstract: We introduce two 2-variables transforms: the partial bi-free S-transform and the partial bi-free T-transform. These transforms are the analogues for the bi-multiplicative and respectively for the additive-multiplicative bi-free convolution of the 2-variables partial bi-free R-transform in our previous paper in this series.


title: Null Spaces of Radon Transforms

abstract: We obtain new descriptions of the null spaces of several projectively equivalent transforms in integral geometry. The paper deals with the hyperplane Radon transform, the totally geodesic transforms on the sphere and the hyperbolic space, the spherical slice transform, and the Cormack-Quinto spherical mean transform for spheres through the origin. The consideration extends to the corresponding dual transforms and the relevant exterior/interior modifications. The method relies on new results for the Gegenbauer-Chebyshev integrals, which generalize Abel type fractional integrals on the positive half-line.


title: $\mathcal{PT}$-breaking threshold in spatially asymmetric Aubry-Andre Harper models: hidden symmetry and topological states

abstract: Aubry-Andre Harper (AAH) lattice models, characterized by reflection-asymmetric, sinusoidally varying nearest-neighbor tunneling profile, are well-known for their topological properties. We consider the fate of such models in the presence of balanced gain and loss potentials $\pm i\gamma$ located at reflection-symmetric sites. We predict that these models have a finite $\mathcal{PT}$ breaking threshold only for {\it specific locations} of the gain-loss potential, and uncover a hidden symmetry that is instrumental to the finite threshold strength. We also show that the topological edge-states remain robust in the $\mathcal{PT}$-symmetry broken phase. Our predictions substantially broaden the possible realizations of a $\mathcal{PT}$-symmetric system.


title: Low reheating temperatures in monomial and binomial inflationary potentials

abstract: We investigate the allowed range of reheating temperature values in light of the Planck 2015 results and the recent joint analysis of Cosmic Microwave Background (CMB) data from the BICEP2/Keck Array and Planck experiments, using monomial and binomial inflationary potentials. While the well studied $\phi^2$ inflationary potential is no longer favored by current CMB data, as well as $\phi^p$ with $p>2$, a $\phi^1$ potential and canonical reheating ($w_{re}=0$) provide a good fit to the CMB measurements. In this last case, we find that the Planck 2015 $68\%$ confidence limit upper bound on the spectral index, $n_s$, implies an upper bound on the reheating temperature of $T_{re}\lesssim 6\times 10^{10}\,{\rm GeV}$, and excludes instantaneous reheating. The low reheating temperatures allowed by this model open the possiblity that dark matter could be produced during the reheating period instead of when the Universe is radiation dominated, which could lead to very different predictions for the relic density and momentum distribution of WIMPs, sterile neutrinos, and axions. We also study binomial inflationary potentials and show the effects of a small departure from a $\phi^1$ potential. We find that as a subdominant $\phi^2$ term in the potential increases, first instantaneous reheating becomes allowed, and then the lowest possible reheating temperature of $T_{re}=4\,{\rm MeV}$ is excluded by the Planck 2015 $68\%$ confidence limit.


title: Causation Entropy Identifies Indirect Influences, Dominance of Neighbors and Anticipatory Couplings

abstract: Inference of causality is central in nonlinear time series analysis and science in general. A popular approach to infer causality between two processes is to measure the information flow between them in terms of transfer entropy. Using dynamics of coupled oscillator networks, we show that although transfer entropy can successfully detect information flow in two processes, it often results in erroneous identification of network connections under the presence of indirect interactions, dominance of neighbors, or anticipatory couplings. Such effects are found to be profound for time-dependent networks. To overcome these limitations, we develop a measure called causation entropy and show that its application can lead to reliable identification of true couplings.


title: Design Issues of JPQ: a Pattern-based Query Language for Document Databases

abstract: Document databases are becoming popular, but how to present complex document query to obtain useful information from the document remains an important topic to study. In this paper, we describe the design issues of a pattern-based document database query language named JPQ. JPQ uses various expressive patterns to extract and construct document fragments following a JSON-like document data model. It adopts tree-like extraction patterns with a coherent pattern composition mechanism to extract data elements from hierarchically structured documents and maintain the logical relationships among the elements. Based on these relationships, JPQ deploys a deductive mechanism to declaratively specify the data transformation requests and considers also data filtering on hierarchical data structure. We use various examples to show the features of the language and to demonstrate its expressiveness and declarativeness in presenting complex document queries.


title: Active galactic nuclei cores in infrared-faint radio sources: Very long baseline interferometry observations using the Very Long Baseline Array

abstract: Infrared-faint radio sources (IFRS) form a new class of galaxies characterised by radio flux densities between tenths and tens of mJy and faint or absent infrared counterparts. It has been suggested that these objects are radio-loud active galactic nuclei (AGNs) at significant redshifts (z >~ 2). Whereas the high redshifts of IFRS have been recently confirmed based on spectroscopic data, the evidence for the presence of AGNs in IFRS is mainly indirect. So far, only two AGNs have been unquestionably confirmed in IFRS based on very long baseline interferometry (VLBI) observations. In this work, we test the hypothesis that IFRS contain AGNs in a large sample of sources using VLBI. We observed 57 IFRS with the Very Long Baseline Array (VLBA) down to a detection sensitivity in the sub-mJy regime and detected compact cores in 35 sources. Our VLBA detections increase the number of VLBI-detected IFRS from 2 to 37 and provide strong evidence that most - if not all - IFRS contain AGNs. We find that IFRS have a marginally higher VLBI detection fraction than randomly selected sources with mJy flux densities at arcsec-scales. Moreover, our data provide a positive correlation between compactness - defined as the ratio of milliarcsec- to arcsec-scale flux density - and redshift for IFRS, but suggest a decreasing mean compactness with increasing arcsec-scale radio flux density. Based on these findings, we suggest that IFRS tend to contain young AGNs whose jets have not formed yet or have not expanded, equivalent to very compact objects. We found two IFRS that are resolved into two components. The two components are spatially separated by a few hundred milliarcseconds in both cases. They might be components of one AGN, a binary black hole, or the result of gravitational lensing.


title: Continuous decomposition of quantum measurements via Hamiltonian feedback

abstract: We characterize the set of generalized quantum measurements that can be decomposed into a continuous measurement process using a stream of probe qubits and a tunable interaction Hamilto- nian. Each probe in the stream interacts weakly with the target quantum system, then is measured projectively in a standard basis. This measurement result is used in a closed feedback loop to tune the interaction Hamiltonian for the next probe. The resulting evolution is a stochastic process with the structure of a one-dimensional random walk. To maintain this structure, and require that at long times the measurement outcomes be independent of the path, the allowed interaction Hamil- tonians must lie in a restricted set, such that the Hamiltonian terms on the target system form a finite dimensional Jordan algebra. This algebraic structure of the interaction Hamiltonians yields a large class of generalized measurements that can be continuously performed by our scheme, and we fully describe this set.


title: Permutation Symmetry Determines the Discrete Wigner Function

abstract: The Wigner function provides a useful quasiprobability representation of quantum mechanics, with applications in various branches of physics. Many nice properties of the Wigner function are intimately connected with the high symmetry of the underlying operator basis composed of phase point operators: any pair of phase point operators can be transformed to any other pair by a unitary symmetry transformation. We prove that, in the discrete scenario, this permutation symmetry is equivalent to the symmetry group being a unitary 2-design. Such a highly symmetric representation can only appear in odd prime power dimensions besides dimensions 2 and 8. It suffices to single out a unique discrete Wigner function among all possible quasiprobability representations. In the course of our study, we show that this discrete Wigner function is uniquely determined by Clifford covariance, while no Wigner function is Clifford covariant in any even prime power dimension.


title: Frequency-Bin Entanglement with Tunable Phase

abstract: We describe a technique to produce narrow-band photon pairs with frequency-bin entanglement, whose relative phase can be tuned using linear polarization optics. We show that, making use of the polarization-frequency coupling effect, the phase of a complex polarizer can be transferred into the frequency entanglement.


title: Invariant subspaces for some surface groups acting on A2-Euclidean buildings

abstract: This paper deals with non-Archimedean representations of punctured surface groups in PGL(3), associated actions on Euclidean buildings (of type A2), and degenerations of real convex projective structures on surfaces. The main result is that, under good conditions on Fock-Goncharov generalized shear parameters, non-Archimedean representations acting on the Euclidean building preserve a cocompact weakly convex subspace, which is part flat surface and part tree. In particular the eigenvalue and length(s) spectra are given by an explicit finite A2-complex. We use this result to describe degenerations of real convex projective structures on surfaces for an open cone of parameters. The main tool is a geometric interpretation of Fock-Goncharov parametrization in A2-buildings.


title: On the effects of self- and cross-phase modulation on photon purity for four-wave mixing photon-pair sources

abstract: We consider the effect of self-phase modulation and cross-phase modulation on the joint spectral amplitude of photon pairs generated by spontaneous four-wave mixing. In particular, the purity of a heralded photon from a pair is considered, in the context of schemes that aim to maximise the purity and minimise correlation in the joint spectral amplitude using birefringent phase-matching and short pump pulses. We find that non-linear phase modulation effects will be detrimental, and will limit the quantum interference visibility that can be achieved at a given generation rate. An approximate expression for the joint spectral amplitude with phase modulation is found by considering the group velocity walk-off between each photon and the pump, but neglecting the group-velocity dispersion at each wavelength. The group-velocity dispersion can also be included with a numerical calculation, and it is shown that it only has a small effect on the purity for the realistic parameters considered.


title: Near-Optimal Hybrid Processing for Massive MIMO Systems via Matrix Decomposition

abstract: For the practical implementation of massive multiple-input multiple-output (MIMO) systems, the hybrid processing (precoding/combining) structure is promising to reduce the high cost rendered by large number of RF chains of the traditional processing structure. The hybrid processing is performed through low-dimensional digital baseband processing combined with analog RF processing enabled by phase shifters. We propose to design hybrid RF and baseband precoders/combiners for multi-stream transmission in point-to-point massive MIMO systems, by directly decomposing the pre-designed unconstrained digital precoder/combiner of a large dimension. The constant amplitude constraint of analog RF processing results in the matrix decomposition problem non-convex. Based on an alternate optimization technique, the non-convex matrix decomposition problem can be decoupled into a series of convex sub-problems and effectively solved by restricting the phase increment of each entry in the RF precoder/combiner within a small vicinity of its preceding iterate. A singular value decomposition based technique is proposed to secure an initial point sufficiently close to the global solution of the original non-convex problem. Through simulation, the convergence of the alternate optimization for such a matrix decomposition based hybrid processing (MD-HP) scheme is examined, and the performance of the MD-HP scheme is demonstrated to be near-optimal.


title: End-to-end verifiability

abstract: This pamphlet describes end-to-end election verifiability (E2E-V) for a nontechnical audience: election officials, public policymakers, and anyone else interested in secure, transparent, evidence-based electronic elections. This work is part of the Overseas Vote Foundation's End-to-End Verifiable Internet Voting: Specification and Feasibility Assessment Study (E2E VIV Project), funded by the Democracy Fund.


title: Proof of Heisenberg's error-disturbance principle

abstract: In Heisenberg's error-disturbance relation for electron position measurement, the measurement error must be the one that determines the uncertainty in the electron position just after the measurement. It is the resolution $\epsilon(x_t)$ i.e., the measurement error of the post-measurement observable $\hat{x}_t$, not the precision $\epsilon(x_0)$ i.e., that of the pre-measurement observable $\hat{x}_0$ that determines the uncertainty of the observable $\hat{x}_t$. Therefore, Heisenberg's relation must be interpreted as one between the resolution $\epsilon(x_t)$ and the disturbance $\eta(y_0)$. The magnitude of the disturbance $\eta(y_0)$ is independent of the definition of the measurement value of the observable $\hat{x}_t$. Heisenberg's error-disturbance relation is proven to hold true in general, when the measurement value of $\hat{x}_t$ is defined in order to minimize the resolution $\epsilon(x_t)$.


title: Some remarks on Willmore surfaces embedded in $\mathbb{R}^3$

abstract: Let $f:\mathbb{C}\rightarrow \mathbb{R}^3$ be complete Willmore immersion with $\int_{\Sigma}|A_f|^2<+\infty$. We will show that if $f$ is the limit of an embedded surface sequence, then $f$ is a plane. As an application, we prove that if $\Sigma_k$ is a sequence of closed Willmore surface embedded in $\mathbb{R}^3$ with $W(\Sigma_k)<C$, and if the conformal class of $\Sigma_k$ converges in the moduli space, then we can find a M\"obius transformation $\sigma_k$, such that a subsequence of $\sigma_k(\Sigma_k)$ converges smoothly.


title: Large deviations for two scale chemical kinetic processes

abstract: We formulate the large deviations for a class of two scale chemical kinetic processes motivated from biological applications. The result is successfully applied to treat a genetic switching model with positive feedbacks. The corresponding Hamiltonian is convex with respect to the momentum variable as a by-product of the large deviation theory. This property ensures its superiority in the rare event simulations compared with the result obtained by formal WKB asymptotics. The result is of general interest to understand the large deviations for multiscale problems.


title: A Ratio of Alternants Formula for Loop Schur Functions

abstract: Lam and Pylyavskyy introduced loop symmetric functions as a generalization of symmetric functions. They defined loop Schur functions as generating functions over semistandard tableaux with respect to a `colored weight,' and they proved a Jacobi--Trudi-style determinantal formula for these generating functions. We prove that loop Schur functions can be expressed as a ratio of `loop alternants,' extending the analogy with Schur functions. As an application, we give a new proof of the loop version of the Murnaghan--Nakayama rule.


title: A Sepak Takraw - Based Molecular Model For C60 - A Mathematical Study Of A 60-Atom Molecule

abstract: We mathematically investigate four molecular models of buckminsterfullerene (C60) discussing the strengths and weaknesses of each, and a new orthorhombic 20-dodecahedron model is proposed to replace the traditional truncated icosahedron model. This representation as a sepak takraw ball rather than a soccer ball allows for capturing the positional relations and electron orbits of the 60 atoms comprising the full molecule.


title: Single Transverse Spin Asymmetries in Semi-inclusive Deep Inelastic Scattering in a Spin-1 Diquark Model

abstract: The observed results for the azimuthal single spin asymmetries (SSAs) of the proton, measured in the semi-inclusive deep inelastic scattering (SIDIS), can be explained by the final-state interaction (FSI) from the gluon exchange between the outgoing quark and the target spectator system. SSAs require a phase difference between two amplitudes coupling the target with opposite spins to the same final state. We have used the model of light front wave functions (LFWFs) consisting of a spin-$\frac{1}{2}$ system as a composite of a spin-$\frac{1}{2}$ fermion and a spin-1 vector boson to estimate the SSAs. The implications of such a model have been investigated in detail by considering different coupling constants. The FSIs also produce a complex phase which can be included in the LFWFs to calculate the Sivers and Boer-Mulders distribution functions of the nucleon.