title: A Coinductive Approach to Computing with Compact Sets

abstract: Exact representations of real numbers such as the signed digit representation or more generally linear fractional representations or the infinite Gray code represent real numbers as infinite streams of digits. In earlier work by the first author it was shown how to extract certified algorithms working with the signed digit representations from constructive proofs. In this paper we lay the foundation for doing a similar thing with nonempty compact sets. It turns out that a representation by streams of finitely many digits is impossible and instead trees are needed.


title: The Homogeneous Spectrum of Milnor-Witt $K$-Theory

abstract: For any field $F$ (of characteristic not equal to 2), we determine the Zariski spectrum of homogeneous prime ideals in $K^{MW}_*(F)$, the Milnor-Witt $K$-theory ring of $F$. As a corollary, we recover Lorenz and Leicht's classical result on prime ideals in the Witt ring of $F$. Our computation can be seen as a first step in Balmer's program for studying the tensor triangular geometry of the stable motivic homotopy category.


title: Topologies of nodal sets of random band limited functions

abstract: It is shown that the topologies and nestings of the zero and nodal sets of random (Gaussian) band limited functions have universal laws of distribution. Qualitative features of the supports of these distributions are determined. In particular the results apply to random monochromatic waves and to random real algebraic hyper-surfaces in projective space.


title: Drinfel'd twist of multiplier Hopf algebras

abstract: This paper generalizes the Drinfel'd twist to the multiplier Hopf algebra case. For a multiplier Hopf algebra $A$ with a twist $J$, we construct a new multiplier Hopf algebra $A^{J}$. Furthermore, if $A$ is quasitriangular, then $A^{J}$ is also. Finally, for a counimodular algebraic quantum group $A$, $A^{J}$ is an algebraic quantum group.


title: Cyclic sieving and rational Catalan theory

abstract: Let $a < b$ be coprime positive integers. Armstrong, Rhoades, and Williams defined a set $\mathsf{NC}(a,b)$ of `rational noncrossing partitions', which form a subset of the ordinary noncrossing partitions of $\{1, 2, \dots, b-1\}$. Confirming a conjecture of Armstrong et. al., we prove that $\mathsf{NC}(a,b)$ is closed under rotation and prove an instance of the cyclic sieving phenomenon for this rotational action. We also define a rational generalization of the $\mathfrak{S}_a$-noncrossing parking functions of Armstrong, Reiner, and Rhoades.


title: Ground-state energy, density profiles, and momentum distribution of attractively interacting 1D Fermi gases with hard-wall boundaries: a Monte Carlo study

abstract: Motivated by the realization of hard-wall boundary conditions in experiments with ultracold atoms, we investigate the ground-state properties of spin-1/2 fermions with attractive interactions in a one-dimensional box. We use lattice Monte Carlo methods to determine essential quantities like the energy, which we compute as a function of coupling strength and particle number in the regime from few to many particles. Many-fermion systems bound by hard walls display non-trivial density profiles characterized by so-called Friedel oscillations (which are similar to those observed in harmonic traps). In non-interacting systems, the characteristic length scale of the oscillations is set by (2 kF)^(-1), where kF is the Fermi momentum, while repulsive interactions tend to generate Wigner-crystal oscillations of period (4 kF)^(-1). Based on the non-interacting result, we find a remarkably simple parametrization of the density profiles of the attractively interacting case, which we generalize to the one-body density matrix. While the total momentum is not a conserved quantity in the presence of hard walls, the magnitude of the momentum does provide a good quantum number. We are therefore able to provide a detailed characterization of the (quasi-)momentum distribution, which displays rather robust discontinuity at the Fermi surface. In addition, we determine the spatially varying on-site density-density correlation, which in turn yields Tan's contact density and, upon integration, Tan's contact. As is well known, the latter fully determines the short-range correlations and plays a crucial role in a multitude of equilibrium and non-equilibrium sum rules.


title: Delaunay-type singular solutions for the fractional Yamabe Problem

abstract: We construct Delaunay-type solutions for the fractional Yamabe problem with an isolated singularity $(-\Delta)^\gamma w = c_{n, \gamma} w^{\frac{n+2\gamma}{n-2\gamma}}, w>0 \ \mbox{in} \ \mathbb{R}^n \backslash \{0\}$ We follow a variational approach, in which the key is the computation of the fractional Laplacian in polar coordinates.


title: Ghost Cities Analysis Based on Positioning Data in China

abstract: Real estate projects are developed excessively in China in this decade. Many new housing districts are built, but they far exceed the actual demand in some cities. These cities with a high housing vacancy rate are called ghost cities. The real situation of vacant housing areas in China has not been studied in previous research. This study, using Baidu positioning data, presents the spatial distribution of the vacant housing areas in China and classifies cities with a large vacant housing area as cities or tourism sites. To the best of our knowledge, it is the first time that we detected and analyzed the ghost cities in China at such fine scale. To understand the human dynamic in ghost cities, we select one city and one tourism sites as cases to analyze the features of human dynamics. This study illustrates the capability of big data in sensing our cities objectively and comprehensively.


title: Regular Symmetry Patterns (Technical Report)

abstract: Symmetry reduction is a well-known approach for alleviating the state explosion problem in model checking. Automatically identifying symmetries in concurrent systems, however, is computationally expensive. We propose a symbolic framework for capturing symmetry patterns in parameterised systems (i.e. an infinite family of finite-state systems): two regular word transducers to represent, respectively, parameterised systems and symmetry patterns. The framework subsumes various types of symmetry relations ranging from weaker notions (e.g. simulation preorders) to the strongest notion (i.e. isomorphisms). Our framework enjoys two algorithmic properties: (1) symmetry verification: given a transducer, we can automatically check whether it is a symmetry pattern of a given system, and (2) symmetry synthesis: we can automatically generate a symmetry pattern for a given system in the form of a transducer. Furthermore, our symbolic language allows additional constraints that the symmetry patterns need to satisfy to be easily incorporated in the verification/synthesis. We show how these properties can help identify symmetry patterns in examples like dining philosopher protocols, self-stabilising protocols, and prioritised resource-allocator protocol. In some cases (e.g. Gries's coffee can problem), our technique automatically synthesises a safety-preserving finite approximant, which can then be verified for safety solely using a finite-state model checker.


title: Low-Complexity Channel Reconstruction Methods Based on SVD-ZF Precoding in Massive 3D-MIMO Systems

abstract: In this paper, we study the low-complexity channel reconstruction methods for downlink precoding in massive multiple-Input multiple-Output (MIMO) systems. When the user is allocated less streams than the number of its antennas, the base station (BS) or user usually utilizes the singular value decomposition (SVD) to get the effective channels, whose dimension is equal to the number of streams. This process is called channel reconstruction and done in BS for time division duplex (TDD) mode. However, with the increasing of antennas in BS, the computation burden of SVD is getting incredible. Here, we propose a series of novel low-complexity channel reconstruction methods for downlink precoding in 3D spatial channel model. We consider different correlations between elevation and azimuth antennas, and divide the large dimensional matrix SVD into two kinds of small-size matrix SVD. The simulation results show that the proposed methods only produce less than 10% float computation than the traditional SVD zero-forcing (SVD-ZF) precoding method, while keeping nearly the same performance of 1Gbps.


title: An Evaluation of Sparse Inverse Covariance Models for Group Functional Connectivity in Molecular Imaging

abstract: Evaluating the functional relationships between brain regions measured with neuroimaging provides insight into how the brain is sharing information at a macro scale. Many functional connectivity methods have been developed for dynamic imaging modalities such as functional MRI (fMRI), while less work has focused on models for static molecular imaging techniques such as FDG-PET and Tc-99m HMPAO SPECT across groups of individuals. In this work we provide a quantitative assessment of how well three functional connec- tivity models based on sparse inverse covariance estimation can accurately recover gold standard connectivity patterns across multiple cohorts and data set sizes. We compare the accuracies of learning regularized inverse covariance matrices across cohorts independently with those learned using two different group-based regular- ization models. By using large cohorts of SPECT scans, we are able to provide a quantitative assessment of the accuracy of the models in recovering the gold standard functional conn


title: Nonlinear thermoelectric response due to energy-dependent transport properties of a quantum dot

abstract: Quantum dots are useful model systems for studying quantum thermoelectric behavior because of their highly energy-dependent electron transport properties, which are tunable by electrostatic gating. As a result of this strong energy dependence, the thermoelectric response of quantum dots is expected to be nonlinear with respect to an applied thermal bias. However, until now this effect has been challenging to observe because, first, it is experimentally difficult to apply a sufficiently large thermal bias at the nanoscale and, second, it is difficult to distinguish thermal bias effects from purely temperature-dependent effects due to overall heating of a device. Here we take advantage of a novel thermal biasing technique and demonstrate a nonlinear thermoelectric response in a quantum dot which is defined in a heterostructured semiconductor nanowire. We also show that a theoretical model based on the Master equations fully explains the observed nonlinear thermoelectric response given the energy-dependent transport properties of the quantum dot.


title: Architectural Consistency Checking in Plugin-Based Software Systems

abstract: Manually ensuring that the implementation of a software system is consistent with the software architecture is a laborious and error-prone task. Thus, a variety of approaches towards automated consistency checking have been developed to counteract architecture erosion. However, these approaches lack means to define and check architectural restrictions concerning plugin dependencies, which is required for plugin-based software systems. In this paper, we propose a domain-specific language called Dependency Constraint Language (DepCoL) to facilitate the definition of constraints concerning plugin dependencies. Using DepCoL, it is possible to define constraints affecting groups of plugins, reducing the required specification effort, to formulate constraints for specific plugins only and to refine constraints. Moreover, we provide an Eclipse plugin, which checks whether the software system under development is consistent with the modeled constraints. This enables a seamless integration into the development process to effortless check consistency during development of the software system. In this way, developers are informed about dependency violations immediately and this supports developers in counteracting architecture erosion.


title: Lower bound of the asymptotic complexity of self-similar fractal graphs

abstract: We study the asymptotic complexity constant of the sequence of approximating graphs to a fully symmetric self-similar structure on a finitely ramified fractal $K$. We show how full symmetry implies existence of the asymptotic complexity constant and obtain a sharp lower bound thereby answering two conjectures by Anema.


title: Robust Gaussian Graphical Modeling with the Trimmed Graphical Lasso

abstract: Gaussian Graphical Models (GGMs) are popular tools for studying network structures. However, many modern applications such as gene network discovery and social interactions analysis often involve high-dimensional noisy data with outliers or heavier tails than the Gaussian distribution. In this paper, we propose the Trimmed Graphical Lasso for robust estimation of sparse GGMs. Our method guards against outliers by an implicit trimming mechanism akin to the popular Least Trimmed Squares method used for linear regression. We provide a rigorous statistical analysis of our estimator in the high-dimensional setting. In contrast, existing approaches for robust sparse GGMs estimation lack statistical guarantees. Our theoretical results are complemented by experiments on simulated and real gene expression data which further demonstrate the value of our approach.


title: Time-delayed models of gene regulatory networks

abstract: In this review we discuss different mathematical models of gene regulatory networks as relevant to the onset and development of cancer. After discussion of alternative modelling approaches, we use a paradigmatic two-gene network to focus on the role played by time delays in the dynamics of gene regulatory networks. We contrast the dynamics of the reduced model arising in the limit of fast mRNA dynamics with that of the full model. The review concludes with the discussion of some open problems.


title: Affleck-Dine leptogenesis and its backreaction to inflaton dynamics

abstract: We investigate the backreaction of the Affleck-Dine leptogenesis to inflaton dynamics in the F-term hybrid and chaotic inflation models in supergravity. We determine the lightest neutrino mass in both models so that the predictions of spectral index, tensor-to-scalar ratio, and baryon abundance are consistent with observations.


title: Chaos in the Honeycomb Optical Lattice Unit Cell

abstract: Natural and artificial honeycomb lattices are of great interest because the band structure of these lattices, if properly constructed, contains a Dirac point. Such lattices occur naturally in the form of graphene and carbon nanotubes. They have been created in the lab in the form of semiconductor 2DEGs, optical lattices, and photonic crystals. We show that, over a wide energy range, gases (of electrons, atoms, or photons) that propagate through these lattices are Lorentz gases and the corresponding classical dynamics is chaotic. Thus, honeycomb lattices are also of interest for understanding eigenstate thermalization and the conductor-insulator transition due to dynamic Anderson localization.


title: Nonextensive Nambu Jona-Lasinio model of QCD matter revisited

abstract: We present a revisited version of the nonextensive QCD-based Nambu - Jona-Lasinio (NJL) model describing the behavior of strongly interacting matter proposed by us some time ago. As before, it is based on the nonextensive generalization of the Boltzmann-Gibbs (BG) statistical mechanics used in the NJL model to its nonextensive version based on Tsallis statistics, but this time it fulfils the basic requirements of thermodynamical consistency. Different ways in which this can be done, connected with different possible choices of the form of the corresponding nonextensive entropies, are presented and discussed in detail. The corresponding results are compared, discussed and confronted with previous findings.


title: Algorithmic Analysis of Qualitative and Quantitative Termination Problems for Affine Probabilistic Programs

abstract: In this paper, we consider termination of probabilistic programs with real-valued variables. The questions concerned are: 1. qualitative ones that ask (i) whether the program terminates with probability 1 (almost-sure termination) and (ii) whether the expected termination time is finite (finite termination); 2. quantitative ones that ask (i) to approximate the expected termination time (expectation problem) and (ii) to compute a bound B such that the probability to terminate after B steps decreases exponentially (concentration problem). To solve these questions, we utilize the notion of ranking supermartingales which is a powerful approach for proving termination of probabilistic programs. In detail, we focus on algorithmic synthesis of linear ranking-supermartingales over affine probabilistic programs (APP's) with both angelic and demonic non-determinism. An important subclass of APP's is LRAPP which is defined as the class of all APP's over which a linear ranking-supermartingale exists. Our main contributions are as follows. Firstly, we show that the membership problem of LRAPP (i) can be decided in polynomial time for APP's with at most demonic non-determinism, and (ii) is NP-hard and in PSPACE for APP's with angelic non-determinism; moreover, the NP-hardness result holds already for APP's without probability and demonic non-determinism. Secondly, we show that the concentration problem over LRAPP can be solved in the same complexity as for the membership problem of LRAPP. Finally, we show that the expectation problem over LRAPP can be solved in 2EXPTIME and is PSPACE-hard even for APP's without probability and non-determinism (i.e., deterministic programs). Our experimental results demonstrate the effectiveness of our approach to answer the qualitative and quantitative questions over APP's with at most demonic non-determinism.


title: Competeing orders in spin-1 and spin-3/2 XXZ Kagome antiferromagnets: A series expansion study

abstract: We study the competition between $\sqrt{3} \times \sqrt{3}$ (RT3) and $q=0$ (Q0) magnetic orders in spin-one and spin-$3/2$ Kagome-lattice XXZ antiferromagnets with varying XY anisotropy parameter $\Delta$, using series expansion methods. The Hamiltonian is split into two parts: an $H_0$ which favors the classical order in the desired pattern and an $H_1$, which is treated in perturbation theory by a series expansion. We find that the ground state energy series for the RT3 and Q0 phases are identical up to sixth order in the expansion, but ultimately a selection occurs, which depends on spin and the anisotropy $\Delta$. Results for ground state energy and the magnetization are presented. These results are compared with recent spin-wave theory and coupled-cluster calculations. The series results for the phase diagram are close to the predictions of spin-wave theory. For the spin-one model at the Heisenberg point ($\Delta=1$), our results are consistent with a vanishing order parameter, that is an absence of a magnetically ordered phase. We also develop series expansions for the ground state energy of the spin-one Heisenberg model in the trimerized phase. We find that the ground state energy in this phase is lower than those of magnetically ordered ones, supporting the existence of a spontaneously trimerized phase in this model.


title: Stochastic Optimally-Tuned Ranged-Separated Hybrid Density Functional Theory

abstract: We develop a stochastic formulation of the optimally-tuned range-separated hybrid density functional theory which enables significant reduction of the computational effort and scaling of the non-local exchange operator at the price of introducing a controllable statistical error. Our method is based on stochastic representations of the Coulomb convolution integral and of the generalized Kohn-Sham density matrix. The computational cost of the approach is similar to that of usual Kohn-Sham density functional theory, yet it provides much more accurate description of the quasiparticle energies for the frontier orbitals. This is illustrated for a series of silicon nanocrystals up to sizes exceeding 3000 electrons. Comparison with the stochastic GW many-body perturbation technique indicates excellent agreement for the fundamental band gap energies, good agreement for the band-edge quasiparticle excitations, and very low statistical errors in the total energy for large systems. The present approach has a major advantage over one-shot GW by providing a self-consistent Hamiltonian which is central for additional post-processing, for example in the stochastic Bethe-Salpeter approach.

categories: physics.chem-ph cond-mat.mes-hall cond-mat.mtrl-sci physics.comp-ph


title: Learning with $\ell^{0}$-Graph: $\ell^{0}$-Induced Sparse Subspace Clustering

abstract: Sparse subspace clustering methods, such as Sparse Subspace Clustering (SSC) \cite{ElhamifarV13} and $\ell^{1}$-graph \cite{YanW09,ChengYYFH10}, are effective in partitioning the data that lie in a union of subspaces. Most of those methods use $\ell^{1}$-norm or $\ell^{2}$-norm with thresholding to impose the sparsity of the constructed sparse similarity graph, and certain assumptions, e.g. independence or disjointness, on the subspaces are required to obtain the subspace-sparse representation, which is the key to their success. Such assumptions are not guaranteed to hold in practice and they limit the application of sparse subspace clustering on subspaces with general location. In this paper, we propose a new sparse subspace clustering method named $\ell^{0}$-graph. In contrast to the required assumptions on subspaces for most existing sparse subspace clustering methods, it is proved that subspace-sparse representation can be obtained by $\ell^{0}$-graph for arbitrary distinct underlying subspaces almost surely under the mild i.i.d. assumption on the data generation. We develop a proximal method to obtain the sub-optimal solution to the optimization problem of $\ell^{0}$-graph with proved guarantee of convergence. Moreover, we propose a regularized $\ell^{0}$-graph that encourages nearby data to have similar neighbors so that the similarity graph is more aligned within each cluster and the graph connectivity issue is alleviated. Extensive experimental results on various data sets demonstrate the superiority of $\ell^{0}$-graph compared to other competing clustering methods, as well as the effectiveness of regularized $\ell^{0}$-graph.


title: Caustic Structures and Detectability of Circumbinary Planets in Microlensing

abstract: Recent discoveries of circumbinary planets in Kepler data show that there is a viable channel of planet formation around binary main sequence stars. Motivated by these discoveries, we have investigated the caustic structures and detectability of circumbinary planets in microlensing events. We have produced a suite of animations of caustics as a function of the projected separation and angle of the binary host to efficiently explore caustic structures over the entire circumbinary parameter space. Aided by these animations, we have derived a semi-empirical analytic expression for the location of planetary caustics, which are displaced in circumbinary lenses relative to those of planets with a single host. We have used this expression to show that the dominant source of caustic motion will be due to the planet's orbital motion and not that of the binary star. Finally, we estimate the fraction of circumbinary microlensing events that are recognizable as such to be significant (5-50 percent) for binary projected separations in the range 0.1-0.5 in units of Einstein radii.


title: The VMC Survey - XX. Identification of new Cepheids in the Small Magellanic Cloud

abstract: We present Ks -band light curves for 299 Cepheids in the Small Magellanic Cloud (SMC) of which 288 are new discoveries that we have identified using multi-epoch near-infrared photometry obtained by the VISTA survey of the Magellanic Clouds system (VMC). The new Cepheids have periods in the range from 0.34 to 9.1 days and cover the magnitude interval 12.9 <= Ks <= 17.6 mag. Our method was developed using variable stars previously identified by the optical microlensing survey OGLE. We focus on searching new Cepheids in external regions of the SMC for which complete VMC Ks-band observations are available and no comprehensive identification of different types of variable stars from other surveys exists yet.


title: Energy dissipation of energetic electrons in the inhomogeneous intergalactic medium during the epoch of reionization

abstract: We explore a time-dependent energy dissipation of the energetic electrons in the inhomogeneous intergalactic medium (IGM) during the epoch of cosmic reionization. In addition to the atomic processes we take into account the Inverse Compton (IC) scattering of the electrons on the comic microwave background (CMB) photons, which is the dominant channel of energy loss for the electrons with energies above a few MeV. We show that: (1) the effect on the IGM has both local (atomic processes) and non-local (IC radiation) components; (2) the energy distribution between Hydrogen and Helium ionizations depends on the initial electron energy; (3) the local baryon overdensity significantly affects the fractions of energy distributed in each channel; and (4) the relativistic effect of atomic cross section become important during the epoch of cosmic reionization. We release our code as open source for further modification by the community.


title: Modeling Boyciana-fish-human Interaction with Partial Differential Algebraic Equations

abstract: With human social behaviors influence, some boyciana-fish reaction-diffusion system coupled with elliptic human distribution equation is considered. Firstly, under homogeneous Neumann boundary conditions and ratio-dependent functional response the system can be described as a nonlinear partial differential algebraic equations (PDAEs) and the corresponding linearized system is discussed with singular system theorem. In what follows we discuss the elliptic subsystem and show that the three kinds of nonnegative are corresponded to three different human interference conditions: human free, overdevelopment and regular human activity. Next we examine the system persistence properties: absorbtion region and the stability of positive steady states of three systems. And the diffusion-driven unstable property is also discussed. Moreover, we propose some energy estimation discussion to reveal the dynamic property among the boyciana-fish-human interaction systems.Finally, using the realistic data collected in the past fourteen years, by PDAEs model parameter optimization, we carry out some predicted results about wetland boyciana population. The applicability of the proposed approaches are confirmed analytically and are evaluated in numerical simulations.


title: Automatic Synthesis of Geometry Problems for an Intelligent Tutoring System

abstract: This paper presents an intelligent tutoring system, GeoTutor, for Euclidean Geometry that is automatically able to synthesize proof problems and their respective solutions given a geometric figure together with a set of properties true of it. GeoTutor can provide personalized practice problems that address student deficiencies in the subject matter.


title: Gravitational lensing by Reissner-Nordstr\"om black holes with topological defects

abstract: Using a new geometrical approach introduced by Gibbons and Werner we study the deflection angle by Reissner-Nordstr\"{o}m black holes in the background spacetimes with a global monopole and a cosmic string. By calculating the corresponding optical Gaussian curvature and applying the Gauss-Bonnet theorem to the optical metric we find the leading terms of the deflection angle in the weak limit approximation. We find that the deflection angle increases due to the presence of topological defects.


title: $Z_2$ SIMP Dark Matter

abstract: Dark matter with strong self-interactions provides a compelling solution to several small-scale structure puzzles. Under the assumption that the coupling between dark matter and the Standard Model particles is suppressed, such strongly interacting massive particles (SIMPs) allow for a successful thermal freeze-out through N-to-N' processes, where N dark matter particles annihilate to N' of them. In the most common scenarios, where dark matter stability is guaranteed by a $\mathbb{Z}_2$ symmetry, the seemingly leading annihilating channel, i.e. 3-to-2 process, is forbidden, so the 4-to-2 one dominate the production of the dark matter relic density. Moreover, cosmological observations require that the dark matter sector is colder than the thermal bath of Standard Model particles, a condition that can be dynamically generated via a small portal between dark matter and Standard Model particles, \`a la freeze-in. This scenario is exemplified in the context of the Singlet Scalar dark matter model.


title: On Equivariant Elliptic Genera of Toric Calabi-Yau 3-folds

abstract: We show that equivariant elliptic genera of toric Calabi-Yau 3-folds are generalized weak Jacobi forms. We also introduce a notion of averaged equivariant elliptic genera of toric Calabi-Yau 3-folds, and show that they are ordinary weak Jacobi forms given by an explicit formula predicted by Eguchi and Sugawara.


title: A Digital Binomial Theorem for Sheffer Sequences

abstract: We extend the digital binomial theorem to Sheffer polynomial sequences by demonstrating that their corresponding Sierpi\'nski matrices satisfy a multiplication property that is equivalent to the convolution identity for Sheffer sequences.


title: SAID: A Control Protocol for Scalable and Adaptive Information Dissemination in ICN

abstract: Information dissemination applications (video, news, social media, etc.) with large number of receivers need to be efficient but also have limited loss tolerance. The new Information-Centric Networks (ICN) paradigm offers an alternative approach for reliably delivering data by naming content and exploiting data available at any intermediate point (e.g., caches). However, receivers are often heterogeneous, with widely varying receive rates. When using existing ICN congestion control mechanisms with in-sequence delivery, a particularly thorny problem of receivers going out-of-sync results in inefficiency and unfairness with heterogeneous receivers. We argue that separating reliability from congestion control leads to more scalable, efficient and fair data dissemination, and propose SAID, a Control Protocol for Scalable and Adaptive Information Dissemination in ICN. To maximize the amount of data transmitted at the first attempt, receivers request any next packet (ANP) of a flow instead of next-in-sequence packet, independent of the provider's transmit rate. This allows providers to transmit at an application-efficient rate, without being limited by the slower receivers. SAID ensures reliable delivery to all receivers eventually, by cooperative repair, while preserving privacy without unduly trusting other receivers.


title: New Threats to SMS-Assisted Mobile Internet Services from 4G LTE: Lessons Learnt from Distributed Mobile-Initiated Attacks towards Facebook and Other Services

abstract: Mobile Internet is becoming the norm. With more personalized mobile devices in hand, many services choose to offer alternative, usually more convenient, approaches to authenticating and delivering the content between mobile users and service providers. One main option is to use SMS (i.e., short messaging service). Such carrier-grade text service has been widely used to assist versatile mobile services, including social networking, banking, to name a few. Though the text service can be spoofed via certain Internet text service providers which cooperated with carriers, such attacks haven well studied and defended by industry due to the efforts of research community. However, as cellular network technology advances to the latest IP-based 4G LTE, we find that these mobile services are somehow exposed to new threats raised by this change, particularly on 4G LTE Text service (via brand-new distributed Mobile-Initiated Spoofed SMS attack which is not available in legacy 2G/3G systems). The reason is that messaging service over LTE shifts from the circuit-switched (CS) design to the packet-switched (PS) paradigm as 4G LTE supports PS only. Due to this change, 4G LTE Text Service becomes open to access. However, its shields to messaging integrity and user authentication are not in place. As a consequence, such weaknesses can be exploited to launch attacks (e.g., hijack Facebook accounts) against a targeted individual, a large scale of mobile users and even service providers, from mobile devices. Current defenses for Internet-Initiated Spoofed SMS attacks cannot defend the unprecedented attack. Our study shows that 53 of 64 mobile services over 27 industries are vulnerable to at least one threat. We validate these proof-of-concept attacks in one major US carrier which supports more than 100 million users. We finally propose quick fixes and discuss security insights and lessons we have learnt.


title: The Singular Value Decomposition, Applications and Beyond

abstract: The singular value decomposition (SVD) is not only a classical theory in matrix computation and analysis, but also is a powerful tool in machine learning and modern data analysis. In this tutorial we first study the basic notion of SVD and then show the central role of SVD in matrices. Using majorization theory, we consider variational principles of singular values and eigenvalues. Built on SVD and a theory of symmetric gauge functions, we discuss unitarily invariant norms, which are then used to formulate general results for matrix low rank approximation. We study the subdifferentials of unitarily invariant norms. These results would be potentially useful in many machine learning problems such as matrix completion and matrix data classification. Finally, we discuss matrix low rank approximation and its recent developments such as randomized SVD, approximate matrix multiplication, CUR decomposition, and Nystrom approximation. Randomized algorithms are important approaches to large scale SVD as well as fast matrix computations.


title: A Hubble Space Telescope Study of the Enigmatic Milky Way Halo Globular Cluster Crater

abstract: We analyze the resolved stellar populations of the faint stellar system, Crater, based on deep optical imaging taken with the Hubble Space Telescope. The HST/ACS-based color-magnitude diagram (CMD) of Crater extends $\sim$4 magnitudes below the oldest main sequence turnoff, providing excellent leverage on Crater's physical properties. Structurally, Crater has a half-light radius of $\sim$20 pc and shows no evidence for tidal distortions. Crater is well-described by a simple stellar population with an age of $\sim$7.5 Gyr, [M/H]$\sim-1.65$, a M$_{\star}\sim10^4$ M$_{\odot}$, M$_{\rm V}\sim -5.3$, located at a distance of d$\sim$ 145 kpc, with modest uncertainties in these properties due to differences in the underlying stellar evolution models. The sparse sampling of stars above the turnoff and sub-giant branch are likely to be 1.0-1.4 M$_{\odot}$ binary star systems (blue stragglers) and their evolved descendants, as opposed to intermediate age main sequence stars. Confusion of these populations highlights a substantial challenge in accurately characterizing sparsely populated stellar systems. Our analysis shows that Crater is not a dwarf galaxy, but instead is an unusually young cluster given its location in the Milky Way's very outer stellar halo. Crater is similar to SMC cluster Lindsay 38, and its position and velocity are in good agreement with observations and models of the Magellanic stream debris, suggesting it may have accreted from the Magellanic Clouds. However, its age and metallicity are also in agreement with the age-metallicity relationships of lower mass dwarf galaxies such as Leo I or Carina. Despite uncertainty over its progenitor system, Crater appears to have been incorporated into the Galaxy more recently than $z\sim1$ (8 Gyr ago), providing an important new constraint on the accretion history of the Milky Way. [abridged]


title: The Differential Geometry of Homogeneity Spaces Across Effect Scales

abstract: If an effect measure is more homogeneous than others, then its value is more likely to be stable across different subgroups or subpopulations. Therefore, it is of great importance to find a more homogeneous effect measure that allows for transportability of research results. For a binary outcome, applied researchers often claim that the risk difference is more heterogeneous than the risk ratio or odds ratio, because they find, based on evidence from surveys of meta-analyses, that the null hypotheses of homogeneity are rejected more often for the risk difference than for the risk ratio and odds ratio. However, the evidence for these claims are far from satisfactory, because of different statistical powers of the homogeneity tests under different effect scales. For binary treatment, covariate and outcome, we theoretically quantify the homogeneity of different effect scales. Because when homogeneity holds the four outcome probabilities lie in a three dimensional sub-space of the four dimensional space, we can use results from differential geometry to compute the volumes of these three dimensional spaces to compare the relative homogeneity of the risk difference, risk ratio, and odds ratio. We demonstrate that the homogeneity space for the risk difference has the smallest volume, and the homogeneity space for the odds ratio has the largest volume, providing some further evidence for the previous claim that the risk difference is more heterogeneous than the risk ratio and odds ratio.


title: On the Covering Radius of the Second Order Reed-Muller Code of Length 128

abstract: In 1981, Schatz proved that the covering radius of the binary Reed-Muller code $RM(2,6)$ is 18. For $RM(2,7)$, we only know that its covering radius is between 40 and 44. In this paper, we prove that the covering radius of the binary Reed-Muller code $RM(2,7)$ is at most 42. Moreover, we give a sufficient and necessary condition for Boolean functions of 7-variable to achieve the second-order nonlinearity 42.


title: The Connection between Supernova Remnants and the Galactic Magnetic Field: A Global Radio Study of the Axisymmetric Sample

abstract: The study of supernova remnants (SNRs) is fundamental to understanding the chemical enrichment and magnetism in galaxies, including our own Milky Way. In an effort to understand the connection between the morphology of SNRs and the Galactic magnetic field (GMF), we have examined the radio images of all known SNRs in our Galaxy and compiled a large sample that have an "axisymmetric" morphology, which we define to mean SNRs with a "bilateral" or "barrel"-shaped morphology, in addition to one-sided shells. We selected the cleanest examples and model each of these at their appropriate Galactic position using two GMF models, those of Jansson & Farrar (2012a), which includes a vertical halo component, and Sun et al. (2008) that is oriented entirely parallel to the plane. Since the magnitude and relative orientation of the magnetic field changes with distance from the sun, we analyse a range of distances, from 0.5 to 10 kpc in each case. Using a physically motivated model of a SNR expanding into the ambient GMF, we find the models using Jansson & Farrar (2012a) are able to reproduce observed morphologies of many SNRs in our sample. These results strongly support the presence of an off-plane, vertical component to the GMF, and the importance of the Galactic field on SNR morphology. Our approach also provides a potential new method for determining distances to SNRs, or conversely, distances to features in the large-scale GMF if SNR distances are known.


title: The frequency-localization technique and minimal decay-regularity for Euler-Maxwell equations

abstract: Dissipative hyperbolic systems of \textit{regularity-loss} have been recently received increasing attention. Usually, extra higher regularity is assumed to obtain the optimal decay estimates, in comparison with that for the global-in-time existence of solutions. In this paper, we develop a new frequency-localization time-decay property, which enables us to overcome the technical difficulty and improve the minimal decay-regularity for dissipative systems. As an application, it is shown that the optimal decay rate of $L^1(\mathbb{R}^3)$-$L^2(\mathbb{R}^3)$ is available for Euler-Maxwell equations with the critical regularity $s_{c}=5/2$, that is, the extra higher regularity is not needed.


title: "Hard" crystalline lattice in the Weyl semimetal NbAs

abstract: We report the effect of hydrostatic pressure on the magnetotransport properties of the Weyl semimetal NbAs. Subtle changes can be seen in the $\rho_{xx}(T)$ profiles with pressure up to 2.31 GPa. The Fermi surfaces undergo an anisotropic evolution under pressure: the extremal areas slightly increase in the $\mathbf{k_x}$-$\mathbf{k_y}$ plane, but decrease in the $\mathbf{k_z}$-$\mathbf{k_y}$($\mathbf{k_x}$) plane. The topological features of the two pockets observed at atmospheric pressure, however, remain unchanged at 2.31 GPa. No superconductivity can be seen down to 0.3 K for all the pressures measured. By fitting the temperature dependence of specific heat to the Debye model, we obtain a small Sommerfeld coefficient $\gamma_0=$ 0.09(1) mJ/(mol$\cdot$K$^2$) and a large Debye temperature, $\Theta_D=$ 450(9) K, confirming a "hard" crystalline lattice that is stable under pressure. We also studied the Kadowaki-Woods ratio of this low-carrier-density massless system, $R_{KW}=$ 3.2$\times 10^4$ $\mu\Omega$ cm mol$^2$ K$^2$ J$^{-2}$. After accounting for the small carrier density in NbAs, this $R_{KW}$ indicates a suppressed transport scattering rate relative to other metals.


title: There is Individualized Treatment. Why Not Individualized Inference?

abstract: Doctors use statistics to advance medical knowledge; we use a medical analogy to introduce statistical inference "from scratch" and to highlight an improvement. Your doctor, perhaps implicitly, predicts the effectiveness of a treatment for you based on its performance in a clinical trial; the trial patients serve as controls for you. The same logic underpins statistical inference: to identify the best statistical procedure to use for a problem, we simulate a set of control problems and evaluate candidate procedures on the controls. Now for the improvement: recent interest in personalized/individualized medicine stems from the recognition that some clinical trial patients are better controls for you than others. Therefore, treatment decisions for you should depend only on a subset of relevant patients. Individualized statistical inference implements this idea for control problems (rather than patients). Its potential for improving data analysis matches personalized medicine's for improving healthcare. The central issue--for both individualized medicine and individualized inference--is how to make the right relevance robustness trade-off: if we exercise too much judgement in determining which controls are relevant, our inferences will not be robust. How much is too much? We argue that the unknown answer is the Holy Grail of statistical inference.


title: Hints of a PeV Scale Supersymmetric Neutrino Sector: Neutrino Masses, Dark Matter, keV and PeV Anomalies

abstract: Neutrino masses and light (keV-GeV) sterile neutrinos can arise naturally via a modified, low energy seesaw mechanism if the right-handed neutrinos are charged under a new symmetry broken by a PeV scale vacuum expectation value, presumably tied to supersymmetry breaking. The additional field content also allows for freeze-in production of sterile neutrino dark matter. This framework can accommodate the recently observed 3.5 keV X-ray line, while a straightforward extension of the framework, using the new symmetry and the PeV energy scale, can explain the PeV energy neutrino events at IceCube. Together, these can therefore be taken as hints of the existence of a PeV scale supersymmetric neutrino sector.


title: Measuring $\eta$ Carinae's High Mass Ejecta in the Infrared and Sub-millimeter

abstract: I address uncertainties on the spatial distribution and mass of the dust formed in $\eta$ Carinae's Homunculus nebula with data being combined from several space- and ground-based facilities spanning near-infrared to sub-mm wavelengths, in terms of observational constraints and modeling. Until these aspects are better understood, the mass loss history and mechanisms responsible for $\eta$ Car's enormous eruption(s) remain poorly constrained.


title: Multi-scale structure and topological anomaly detection via a new network statistic: The onion decomposition

abstract: We introduce a new network statistic that measures diverse structural properties at the micro-, meso-, and macroscopic scales, while still being easy to compute and easy to interpret at a glance. Our statistic, the onion spectrum, is based on the onion decomposition, which refines the k-core decomposition, a standard network fingerprinting method. The onion spectrum is exactly as easy to compute as the k-cores: It is based on the stages at which each vertex gets removed from a graph in the standard algorithm for computing the k-cores. But the onion spectrum reveals much more information about a network, and at multiple scales; for example, it can be used to quantify node heterogeneity, degree correlations, centrality, and tree- or lattice-likeness of the whole network as well as of each k-core. Furthermore, unlike the k-core decomposition, the combined degree-onion spectrum immediately gives a clear local picture of the network around each node which allows the detection of interesting subgraphs whose topological structure differs from the global network organization. This local description can also be leveraged to easily generate samples from the ensemble of networks with a given joint degree-onion distribution. We demonstrate the utility of the onion spectrum for understanding both static and dynamic properties on several standard graph models and on many real-world networks.

categories: physics.soc-ph cond-mat.dis-nn cs.DM cs.SI math.CO


title: Micromagnetic simulations with periodic boundary conditions: Hard-soft nanocomposites

abstract: We developed a micromagnetic method for modeling magnetic systems with periodic boundary conditions along an arbitrary number of dimensions. The main feature is an adaptation of the Ewald summation technique for evaluation of long-range dipolar interactions. The method was applied to investigate the hysteresis process in hard-soft magnetic nanocomposites with various geometries. The dependence of the results on different micromagnetic parameters was studied. We found that for layered structures with an out-of-plane hard phase easy axis the hysteretic properties are very sensitive to the strength of the interlayer exchange coupling, as long as the spontaneous magnetization for the hard phase is significantly smaller than for the soft phase. The origin of this behavior was discussed. Additionally, we investigated the soft phase size optimizing the energy product of hard-soft nanocomposites.


title: Empirical Output Distribution of Good Delay-Limited Codes for Quasi-Static Fading Channels

abstract: This paper considers delay-limited communication over quasi-static fading channels under a long-term power constraint. A sequence of length-$n$ delay-limited codes for a quasi-static fading channel is said to be capacity-achieving if the codes achieve the delay-limited capacity, which is defined to be the maximum rate achievable by delay-limited codes. The delay-limited capacity is sometimes referred to as the zero-outage capacity in wireless communications. The delay-limited capacity is the appropriate choice of performance measure for delay-sensitive applications such as voice and video over fading channels. It is shown that for any sequence of capacity-achieving delay-limited codes with vanishing error probabilities, the normalized relative entropy between the output distribution induced by the length-$n$ code and the $n$-fold product of the capacity-achieving output distribution, denoted by $\frac{1}{n}D(p_{Y^n}\|p_{Y^n}^*)$, converges to zero. Additionally, we extend our convergence result to capacity-achieving delay-limited codes with non-vanishing error probabilities.


title: High Energy Physics Forum for Computational Excellence: Working Group Reports (I. Applications Software II. Software Libraries and Tools III. Systems)

abstract: Computing plays an essential role in all aspects of high energy physics. As computational technology evolves rapidly in new directions, and data throughput and volume continue to follow a steep trend-line, it is important for the HEP community to develop an effective response to a series of expected challenges. In order to help shape the desired response, the HEP Forum for Computational Excellence (HEP-FCE) initiated a roadmap planning activity with two key overlapping drivers -- 1) software effectiveness, and 2) infrastructure and expertise advancement. The HEP-FCE formed three working groups, 1) Applications Software, 2) Software Libraries and Tools, and 3) Systems (including systems software), to provide an overview of the current status of HEP computing and to present findings and opportunities for the desired HEP computational roadmap. The final versions of the reports are combined in this document, and are presented along with introductory material.


title: On Differentially Private Online Collaborative Recommendation Systems

abstract: In collaborative recommendation systems, privacy may be compromised, as users' opinions are used to generate recommendations for others. In this paper, we consider an online collaborative recommendation system, and we measure users' privacy in terms of the standard differential privacy. We give the first quantitative analysis of the trade-offs between recommendation quality and users' privacy in such a system by showing a lower bound on the best achievable privacy for any non-trivial algorithm, and proposing a near-optimal algorithm. From our results, we find that there is actually little trade-off between recommendation quality and privacy for any non-trivial algorithm. Our results also identify the key parameters that determine the best achievable privacy.


title: A Smooth and Locally Sparse Estimator for Functional Linear Regression via Functional SCAD Penalty

abstract: In this paper, we propose a new regularization technique called "functional SCAD". We then combine this technique with the smoothing spline method to develop a smooth and locally sparse (i.e., zero on some sub-regions) estimator for the coefficient function in functional linear regression. The functional SCAD has a nice shrinkage property that enables our estimating procedure to identify the null subregions of the coefficient function without over shrinking the non-zero values of the coefficient function. Additionally, the smoothness of our estimated coefficient function is regularized by a roughness penalty rather than by controlling the number of knots. Our method is more theoretically sound and is computationally simpler than the other available methods. An asymptotic analysis shows that our estimator is consistent and can identify the null region with the probability tending to one. Furthermore, simulation studies show that our estimator has superior numerical performance. Finally, the practical merit of our method is demonstrated on two real applications.


title: Chiral Symmetry transition in the Linear Sigma Model with quarks: Counting effective QCD degrees of freedom from low to high temperature

abstract: We use the linear sigma model coupled to quarks, together with a plausible location of the critical end point (CEP), to study the chiral symmetry transition in the QCD phase diagram. We compute the effective potential at finite temperature and density up to the contribution of the ring diagrams, both in the low and high temperature limits, and use it to compute the pressure and the position of the CEP. In the high temperature regime, by comparing to results from extrapolated lattice data, we determine the model coupling constants. Demanding that the CEP remains in the same location when described in the high temperature limit, we determine again the couplings and the pressure for the low temperature regime. We show that this procedure gives an average description of the lattice QCD results for the pressure and that the change from the low to the high temperature domains in this quantity can be attributed to the change in the coupling constants which in turn we link to the change in the effective degrees of freedom.


title: Ballistic annihilation with superimposed diffusion in one dimension

abstract: We consider a one-dimensional system with particles having either positive or negative velocity, which annihilate on contact. To the ballistic motion of the particle, a diffusion is superimposed. The annihilation may represent a reaction in which the two particles yield an inert species. This model has been the object of previous work, in which it was shown that the particle concentration decays faster than either the purely ballistic or the purely diffusive case. We report on previously unnoticed behaviour for large times, when only one of the two species remains and also unravel the underlying fractal structure present in the system. We also consider in detail the case in which the initial concentration of right-going particles is $1/2+\varepsilon$, with $\varepsilon\neq0$. It is shown that a remarkably rich behaviour arises, in which two crossover times are observed as $\varepsilon \to 0$.


title: Allowances for evolving coastal flood risk under uncertain local sea-level rise

abstract: Sea-level rise (SLR) causes estimates of flood risk made under the assumption of stationary mean sea level to be biased low. However, adjustments to flood return levels made assuming fixed increases of sea level are also inaccurate when applied to sea level that is rising over time at an uncertain rate. To accommodate both the temporal dynamics of SLR and their uncertainty, we develop an Average Annual Design Life Level (AADLL) metric and associated Design Life SLR (DL-SLR) allowances. The AADLL is the flood level corresponding to a time-integrated annual expected probability of occurrence (AEP) under uncertainty over the design life of an asset; DL-SLR allowances are the adjustment from 2000 levels that maintain current average probability over the design life. Given non-stationary and uncertain sea-level rise, AADLL flood levels and DL-SLR allowances provide estimates of flood protection heights and offsets for different planning horizons and different levels of confidence in SLR projections in coastal areas. Here we employ probabilistic sea-level rise projections to illustrate the calculation of AADLL flood levels and DL-SLR allowances for a set of long-duration tide gauges along U.S. coastlines.


title: Scheduling under Linear Constraints

abstract: We introduce a parallel machine scheduling problem in which the processing times of jobs are not given in advance but are determined by a system of linear constraints. The objective is to minimize the makespan, i.e., the maximum job completion time among all feasible choices. This novel problem is motivated by various real-world application scenarios. We discuss the computational complexity and algorithms for various settings of this problem. In particular, we show that if there is only one machine with an arbitrary number of linear constraints, or there is an arbitrary number of machines with no more than two linear constraints, or both the number of machines and the number of linear constraints are fixed constants, then the problem is polynomial-time solvable via solving a series of linear programming problems. If both the number of machines and the number of constraints are inputs of the problem instance, then the problem is NP-Hard. We further propose several approximation algorithms for the latter case.


title: Spherically symmetric solution of $f(R,\mathcal{G})$ gravity at low energy

abstract: The weak-field and slow-motion limit of $f(R,\mathcal{G})$ gravity is developed up to $(v/c)^{4}$ order in a spherically symmetric background. Considering the Taylor expansion of a general function $f$ around vanishing values of $R$ and $\mathcal{G}$, we present general vacuum solutions up to $(v/c)^{4}$ order for the gravitational field generated by a ball-like source. The spatial behaviors at $(v/c)^{2}$ order are the same for $f(R,\mathcal{G})$ gravity and $f(R)$ gravity, and their corresponding real valued static behaviors are presented and compared with the one in general relativity. The static Yukawa-like behavior is proved to be compatible with the previous result of the most general fourth-order theory. At $(v/c)^{4}$ order, the static corrections to the Yukawa-like behavior for $f(R,\mathcal{G})$ gravity, $f(R)$ gravity, and the Starobinsky gravity are presented and compared with the one in general relativity.


title: Preheating after Multifield Inflation with Nonminimal Couplings, I: Covariant Formalism and Attractor Behavior

abstract: This is the first of a three-part series of papers, in which we study the preheating phase for multifield models of inflation involving nonminimal couplings. In this paper, we study the single-field attractor behavior that these models exhibit during inflation and quantify its strength and parameter dependence. We further demonstrate that the strong single-field attractor behavior persists after the end of inflation. Preheating in such models therefore generically avoids the "de-phasing" that typically affects multifield models with minimally coupled fields, allowing efficient transfer of energy from the oscillating inflaton condensate(s) to coupled perturbations across large portions of parameter space. We develop a doubly-covariant formalism for studying the preheating phase in such models and identify several features specific to multifield models with nonminimal couplings, including effects that arise from the nontrivial field-space manifold. In papers II and III, we apply this formalism to study how the amplification of adiabatic and isocurvature perturbations varies with parameters, highlighting several distinct regimes depending on the magnitude of the nonminimal couplings $\xi_I$.


title: "We're on the Same Page": A Usability Study of Secure Email Using Pairs of Novice Users

abstract: Secure email is increasingly being touted as usable by novice users, with a push for adoption based on recent concerns about government surveillance. To determine whether secure email is for grassroots adoption, we employ a laboratory user study that recruits pairs of novice to install and use several of the latest systems to exchange secure messages. We present quantitative and qualitative results from 25 pairs of novice users as they use Pwm, Tutanota, and Virtru. Participants report being more at ease with this type of study and better able to cope with mistakes since both participants are "on the same page". We find that users prefer integrated solutions over depot-based solutions, and that tutorials are important in helping first-time users. Hiding the details of how a secure email system provides security can lead to a lack of trust in the system. Participants expressed a desire to use secure email, but few wanted to use it regularly and most were unsure of when they might use it.


title: Why Johnny Still, Still Can't Encrypt: Evaluating the Usability of a Modern PGP Client

abstract: This paper presents the results of a laboratory study involving Mailvelope, a modern PGP client that integrates tightly with existing webmail providers. In our study, we brought in pairs of participants and had them attempt to use Mailvelope to communicate with each other. Our results shown that more than a decade and a half after \textit{Why Johnny Can't Encrypt}, modern PGP tools are still unusable for the masses. We finish with a discussion of pain points encountered using Mailvelope, and discuss what might be done to address them in future PGP systems.


title: Enumerative geometry of elliptic curves on toric surfaces

abstract: We establish the equality of classical and tropical curve counts for elliptic curves on toric surfaces with fixed $j$-invariant, refining results of Mikhalkin and Nishinou--Siebert. As an application, we determine a formula for such counts on $\mathbb P^2$ and all Hirzebruch surfaces. This formula relates the count of elliptic curves with the number of rational curves on the surface satisfying a small number of tangency conditions with the toric boundary. Furthermore, the combinatorial tropical multiplicities of Kerber and Markwig for counts in $\mathbb P^2$ are derived and explained algebro-geometrically, using Berkovich geometry and logarithmic Gromov--Witten theory. As a consequence, a new proof of Pandharipande's formula for counts of elliptic curves in $\mathbb P^2$ with fixed $j$-invariant is obtained.


title: Magnetic brane solutions of Lovelock gravity with nonlinear electrodynamics

abstract: In this paper, we consider logarithmic and exponential forms of nonlinear electrodynamics as a source and obtain magnetic brane solutions of the Lovelock gravity. Although these solutions have no curvature singularity and no horizon, they have a conic singularity with a deficit angle. We investigate the effects of nonlinear electrodynamics and the Lovelock gravity on the value of deficit angle and find that various terms of Lovelock gravity do not affect deficit angle. Next, we generalize our solutions to spinning cases with maximum rotating parameters in arbitrary dimensions and calculate the conserved quantities of the solutions. Finally, we consider nonlinear electrodynamics as a correction of the Maxwell theory and investigate the properties of the solutions.


title: Study of the Tracking Method and Expected Performance of the Silicon Pixel Inner Tracker Applied in BESIII

abstract: The inner drift chamber of the BESIII is encountering serious aging problem after five year's running. For the first layer, the decrease in gas gain is about 26% from 2009 to 2013. The upgrade of the inner tracking detector has become an urgent problem for the BESIII experiment. An inner tracker using CMOS pixel sensors is an important candidate because of its great advantages on spatial resolution and radiation hardness. In order to carry out a Monte Carlo study on the expected performance, a Geant4-based full simulation for the silicon pixel detector has been implemented. The tracking method combining the silicon pixel inner tracker and outer drift chamber has been studied and a preliminary reconstruction software was developed. The Monte Carlo study shows that the performances including momentum resolution, vertex resolution and the tracking efficiency are significantly improved due to the good spatial resolution and moderate material budget of the silicon pixel detector.


title: Techniques in Iterative Proton CT Image Reconstruction

abstract: This is a review paper on some of the physics, modeling, and iterative algorithms in proton computed tomography (pCT) image reconstruction. The primary challenge in pCT image reconstruction lies in the degraded spatial resolution resulting from multiple Coulomb scattering within the imaged object. Analytical models such as the most likely path (MLP) have been proposed to predict the scattered trajectory from measurements of individual proton location and direction before and after the object. Iterative algorithms provide a flexible tool with which to incorporate these models into image reconstruction. The modeling leads to a large and sparse linear system of equations that can efficiently be solved by projection methods-based iterative algorithms. Such algorithms perform projections of the iterates onto the hyperlanes that are represented by the linear equations of the system. They perform these projections in possibly various algorithmic structures, such as block-iterative projections (BIP), string-averaging projections (SAP). These algorithmic schemes allow flexibility of choosing blocks, strings, and other parameters. They also cater for parallel implementations which are apt to further save clock time in computations. Experimental results are presented which compare some of those algorithmic options.


title: Why Random Reshuffling Beats Stochastic Gradient Descent

abstract: We analyze the convergence rate of the random reshuffling (RR) method, which is a randomized first-order incremental algorithm for minimizing a finite sum of convex component functions. RR proceeds in cycles, picking a uniformly random order (permutation) and processing the component functions one at a time according to this order, i.e., at each cycle, each component function is sampled without replacement from the collection. Though RR has been numerically observed to outperform its with-replacement counterpart stochastic gradient descent (SGD), characterization of its convergence rate has been a long standing open question. In this paper, we answer this question by showing that when the component functions are quadratics or smooth and the sum function is strongly convex, RR with iterate averaging and a diminishing stepsize $\alpha_k=\Theta(1/k^s)$ for $s\in (1/2,1)$ converges at rate $\Theta(1/k^{2s})$ with probability one in the suboptimality of the objective value, thus improving upon the $\Omega(1/k)$ rate of SGD. Our analysis draws on the theory of Polyak-Ruppert averaging and relies on decoupling the dependent cycle gradient error into an independent term over cycles and another term dominated by $\alpha_k^2$. This allows us to apply law of large numbers to an appropriately weighted version of the cycle gradient errors, where the weights depend on the stepsize. We also provide high probability convergence rate estimates that shows decay rate of different terms and allows us to propose a modification of RR with convergence rate ${\cal O}(\frac{1}{k^2})$.


title: Magneto-transport behaviour of Bi2Se3-xTex: Role of disorder

abstract: Magnetoresistance and Hall resistance measurements have been carried out in fastcooled single crystals of Bi2Se3-xTex (x: 0 to 2) in 4 to 300 K temperature range, under magnetic fields up to 15 T. The variation of resistivity with temperature that points to a metallic behaviour in Bi2Se3, shows an upturn at low temperatures in the Te doped samples. Magnetoresistance measurements in Bi2Se3 show clear signatures of Shubnikov de Hass oscillations that gets suppressed in the Te doped samples. In the Bi2SeTe2 sample, the magneto-resistance shows a cusp like positive magneto-resistance at low magnetic fields and low temperatures, a feature associated with weak antilocalisation (WAL), that crosses over to negative magneto-resistance at higher fields. The qualitatively different magnetotransport behaviour seen in Bi2SeTe2 as compared to Bi2Se3 is rationalised in terms of the disorder, through an estimate of the carrier density, carrier mobility and an analysis in terms of the Ioffe Regel criterion with support from Hall Effect measurements.


title: Convergence Rate of Incremental Gradient and Newton Methods

abstract: The incremental gradient method is a prominent algorithm for minimizing a finite sum of smooth convex functions, used in many contexts including large-scale data processing applications and distributed optimization over networks. It is a first-order method that processes the functions one at a time based on their gradient information. The incremental Newton method, on the other hand, is a second-order variant which exploits additionally the curvature information of the underlying functions and can therefore be faster. In this paper, we focus on the case when the objective function is strongly convex and present fast convergence results for the incremental gradient and incremental Newton methods under the constant and diminishing stepsizes. For a decaying stepsize rule $\alpha_k = \Theta(1/k^s)$ with $s \in (0,1]$, we show that the distance of the IG iterates to the optimal solution converges at rate ${\cal O}(1/k^{s})$ (which translates into ${\cal O}(1/k^{2s})$ rate in the suboptimality of the objective value). For $s>1/2$, this improves the previous ${\cal O}(1/\sqrt{k})$ results in distances obtained for the case when functions are non-smooth. We show that to achieve the fastest ${\cal O}(1/k)$ rate, incremental gradient needs a stepsize that requires tuning to the strong convexity parameter whereas the incremental Newton method does not. The results are based on viewing the incremental gradient method as a gradient descent method with gradient errors, devising efficient upper bounds for the gradient error to derive inequalities that relate distances of the consecutive iterates to the optimal solution and finally applying Chung's lemmas from the stochastic approximation literature to these inequalities to determine their asymptotic behavior. In addition, we construct examples to show tightness of our rate results.


title: Surface Recombination Limited Lifetimes of Photoexcited Carriers in Few-Layer Transition Metal Dichalcogenide MoS2

abstract: We present results on photoexcited carrier lifetimes in few-layer transition metal dichalcogenide MoS2 using nondegenerate ultrafast optical pump-probe technique. Our results show a sharp increase of the carrier lifetimes with the number of layers in the sample. Carrier lifetimes increase from few tens of picoseconds in monolayer samples to more than a nanosecond in 10-layer samples. The inverse carrier lifetime was found to scale according to the probability of the carriers being present at the surface layers, as given by the carrier wavefunction in few layer samples, which can be treated as quantum wells. The carrier lifetimes were found to be largely independent of the temperature and the inverse carrier lifetimes scaled linearly with the photoexcited carrier density. These observations are consistent with defect-assisted carrier recombination, in which the capture of electrons and holes by defects occurs via Auger scatterings. Our results suggest that carrier lifetimes in few-layer samples are surface recombination limited due to the much larger defect densities at the surface layers compared to the inner layers.


title: Build your own clarithmetic I: Setup and completeness

abstract: Clarithmetics are number theories based on computability logic (see http://www.csc.villanova.edu/~japaridz/CL/ ). Formulas of these theories represent interactive computational problems, and their "truth" is understood as existence of an algorithmic solution. Various complexity constraints on such solutions induce various versions of clarithmetic. The present paper introduces a parameterized/schematic version CLA11(P1,P2,P3,P4). By tuning the three parameters P1,P2,P3 in an essentially mechanical manner, one automatically obtains sound and complete theories with respect to a wide range of target tricomplexity classes, i.e. combinations of time (set by P3), space (set by P2) and so called amplitude (set by P1) complexities. Sound in the sense that every theorem T of the system represents an interactive number-theoretic computational problem with a solution from the given tricomplexity class and, furthermore, such a solution can be automatically extracted from a proof of T. And complete in the sense that every interactive number-theoretic problem with a solution from the given tricomplexity class is represented by some theorem of the system. Furthermore, through tuning the 4th parameter P4, at the cost of sacrificing recursive axiomatizability but not simplicity or elegance, the above extensional completeness can be strengthened to intensional completeness, according to which every formula representing a problem with a solution from the given tricomplexity class is a theorem of the system. This article is published in two parts. The present Part I introduces the system and proves its completeness, while Part II is devoted to proving soundness.


title: Attention with Intention for a Neural Network Conversation Model

abstract: In a conversation or a dialogue process, attention and intention play intrinsic roles. This paper proposes a neural network based approach that models the attention and intention processes. It essentially consists of three recurrent networks. The encoder network is a word-level model representing source side sentences. The intention network is a recurrent network that models the dynamics of the intention process. The decoder network is a recurrent network produces responses to the input from the source side. It is a language model that is dependent on the intention and has an attention mechanism to attend to particular source side words, when predicting a symbol in the response. The model is trained end-to-end without labeling data. Experiments show that this model generates natural responses to user inputs.


title: Build your own clarithmetic II: Soundness

abstract: Clarithmetics are number theories based on computability logic (see http://www.csc.villanova.edu/~japaridz/CL/ ). Formulas of these theories represent interactive computational problems, and their "truth" is understood as existence of an algorithmic solution. Various complexity constraints on such solutions induce various versions of clarithmetic. The present paper introduces a parameterized/schematic version CLA11(P1,P2,P3,P4). By tuning the three parameters P1,P2,P3 in an essentially mechanical manner, one automatically obtains sound and complete theories with respect to a wide range of target tricomplexity classes, i.e. combinations of time (set by P3), space (set by P2) and so called amplitude (set by P1) complexities. Sound in the sense that every theorem T of the system represents an interactive number-theoretic computational problem with a solution from the given tricomplexity class and, furthermore, such a solution can be automatically extracted from a proof of T. And complete in the sense that every interactive number-theoretic problem with a solution from the given tricomplexity class is represented by some theorem of the system. Furthermore, through tuning the 4th parameter P4, at the cost of sacrificing recursive axiomatizability but not simplicity or elegance, the above extensional completeness can be strengthened to intensional completeness, according to which every formula representing a problem with a solution from the given tricomplexity class is a theorem of the system. This article is published in two parts. The previous Part I has introduced the system and proved its completeness, while the present Part II is devoted to proving soundness.


title: Location-Based Beamforming for Rician Wiretap Channels

abstract: We propose a location-based beamforming scheme for wiretap channels, where a source communicates with a legitimate receiver in the presence of an eavesdropper. We assume that the source and the eavesdropper are equipped with multiple antennas, while the legitimate receiver is equipped with a single antenna. We also assume that all channels are in a Rician fading environment, the channel state information from the legitimate receiver is perfectly known at the source, and that the only information on the eavesdropper available at the source is her location. We first describe how the beamforming vector that minimizes the secrecy outage probability of the system is obtained, illustrating its dependence on the eavesdropper's location. We then derive an easy-to-compute expression for the secrecy outage probability when our proposed location-based beamforming is adopted. Finally, we investigate the impact location uncertainty has on the secrecy outage probability, showing how our proposed solution can still allow for secrecy even when the source has limited information on the eavesdropper's location.


title: Feature-Based Diversity Optimization for Problem Instance Classification

abstract: Understanding the behaviour of heuristic search methods is a challenge. This even holds for simple local search methods such as 2-OPT for the Traveling Salesperson problem. In this paper, we present a general framework that is able to construct a diverse set of instances that are hard or easy for a given search heuristic. Such a diverse set is obtained by using an evolutionary algorithm for constructing hard or easy instances that are diverse with respect to different features of the underlying problem. Examining the constructed instance sets, we show that many combinations of two or three features give a good classification of the TSP instances in terms of whether they are hard to be solved by 2-OPT.


title: Exact Lyapunov exponents of the generalized Boole transformations

abstract: The generalized Boole transformations have rich behavior ranging from the \textit{mixing} phase with the Cauchy invariant measure to the \textit{dissipative} phase through the \textit{infinite ergodic} phase with the Lebesgue measure. In this Letter, by giving the proof of mixing property for $0<\alpha<1$ we show an \textit{analytic} formula of the Lyapunov exponents $\lambda$ which are explicitly parameterized in terms of the parameter $\alpha$ of the generalized Boole transformations for the whole region $\alpha>0$ and bridge those three phase \textit{continuously}. We found the different scale behavior of the Lyapunov exponent near $\alpha=1$ using analytic formula with the parameter $\alpha$. In particular, for $0<\alpha<1$, we then prove an existence of extremely sensitive dependency of Lyapunov exponents, where the absolute values of the derivative of Lyapunov exponents with respect to the parameter $\alpha$ diverge to infinity in the limit of $\alpha\to 0$, and $\alpha \to 1$. This result shows the computational complexity on the numerical simulations of the Lyapunov exponents near $\alpha \simeq$ 0, 1.


title: A Semismooth Newton Method for Tensor Eigenvalue Complementarity Problem

abstract: In this paper, we consider the tensor eigenvalue complementarity problem which is closely related to the optimality conditions for polynomial optimization, as well as a class of differential inclusions with nonconvex processes. By introducing an NCP-function, we reformulate the tensor eigenvalue complementarity problem as a system of nonlinear equations. We show that this function is strongly semismooth but not differentiable, in which case the classical smoothing methods cannot apply. Furthermore, we propose a damped semismooth Newton method for tensor eigenvalue complementarity problem. A new procedure to evaluate an element of the generalized Jocobian is given, which turns out to be an element of the B-subdifferential under mild assumptions. As a result, the convergence of the damped semismooth Newton method is guaranteed by existing results. The numerical experiments also show that our method is efficient and promising.


title: Self-Organization and Heating by Inward Diffusion in Magnetospheric Plasmas

abstract: Through the process of inward diffusion, a strongly localized clump of plasma is created in a magnetosphere. The creation of the density gradient, instead of the usual flattening by a diffusion process, can be explained by the topological constraints given by the adiabatic invariants of magnetized particles. After developing a canonical formalism for the standard guiding center dynamics in a dipole magnetic field, we complete our attempt to build a statistical mechanics on a constrained phase space by discussing the construction principles of the associated diffusion operator. We then investigate the heating mechanism associated with inward diffusion: as particles move toward regions of higher magnetic field, they experience preferential heating of the perpendicular (with respect to the magnetic field) temperature in order to preserve the magnetic moment. A relationship between conservation of bounce action and temperature isotropy emerged. We further show that this behavior is scaled by the diffusion parameter of the Fokker-Planck equation. These results are confirmed by numerical simulations.


title: Radiation Hard & High Light Yield Scintillator Search for CMS Phase II Upgrade

abstract: The CMS detector at the LHC requires a major upgrade to cope with the higher instantaneous luminosity and the elevated radiation levels. The active media of the forward backing hadron calorimeters is projected to be radiation-hard, high light yield scintillation materials or similar alternatives. In this context, we have studied various radiation-hard scintillating materials such as Polyethylene Terephthalate (PET), Polyethylene Naphthalate (PEN), High Efficiency Mirror (HEM) and quartz plates with various coatings. The quartz plates are pure Cerenkov radiators and their radiation hardness has been confirmed. In order to increase the light output, we considered organic and inorganic coating materials such as p-Terphenyl (pTp), Anthracene and Gallium-doped Zinc Oxide (ZnO:Ga) that are applied as thin layers on the surface of the quartz plates. Here, we present the results of the related test beam activities, laboratory measurements and recent developments.


title: On the convergence analysis of the optimized gradient method

abstract: This paper considers the problem of unconstrained minimization of smooth convex functions having Lipschitz continuous gradients with known Lipschitz constant. We recently proposed an optimized gradient method (OGM) for this problem and showed that it has a worst-case convergence bound for the cost function decrease that is twice as small as that of Nesterov's fast gradient method (FGM), yet has a similarly efficient practical implementation. Drori showed recently that OGM has optimal complexity over the general class of first-order methods. This optimality makes it important to study fully the convergence properties of OGM. The previous worst-case convergence bound for OGM was derived for only the last iterate of a secondary sequence. This paper provides an analytic convergence bound for the primary sequence generated by OGM. We then discuss additional convergence properties of OGM, including the interesting fact that OGM has two types of worst-case functions: a piecewise affine-quadratic function and a quadratic function. These results help complete the theory of optimal first-order methods for smooth convex minimization.


title: Charge-Transfer Dynamics and Nonlocal Dielectric Permittivities Tuned with Metamaterial Structures as Solvent Analogues

abstract: Charge transfer (CT) is a fundamental and ubiquitous mechanism in biology, physics and chemistry. Here, we evidence that CT dynamics can be altered by multi-layered hyperbolic metamaterial (HMM) substrates. Taking triphenylene:perylene diimide dyad supramolecular self-assemblies as a model system, we reveal longer-lived CT states in presence of HMM structures with both charge separation and recombination characteristic times increased by factors of 2.4 and 1.7, i.e. relative variations of 140 and 73 %, respectively. To rationalise these experimental results in terms of driving force, we successfully introduce image dipole interactions in Marcus theory. The non-local effect herein demonstrated is directly linked to the number of metal-dielectric pairs, can be formalised in the dielectric permittivity, and is presented as a solid analogue to local solvent polarizability effects. This model and extra PH3T:PC60BM results show the generality of this non-local phenomenon and that a wide range of kinetic tailoring opportunities can arise from substrate engineering. This work paves the way toward the design of artificial substrates to control CT dynamics of interest for applications in optoelectronics and chemistry.


title: A high obliquity orbit for the hot-Jupiter HATS-14b transiting a 5400K star

abstract: We report a spin-orbit misalignment for the hot-Jupiter HATS-14b, measuring a projected orbital obliquity of |lambda|= 76 -5/+4 deg. HATS-14b orbits a high metallicity, 5400 K G dwarf in a relatively short period orbit of 2.8 days. This obliquity was measured via the Rossiter-McLaughlin effect, obtained with observations from Keck-HIRES. The velocities were extracted using a novel technique, optimised for low signal-to-noise spectra, achieving a high precision of 4 m/s point-to-point scatter. However, we caution that our uncertainties may be underestimated. Due to the low rotational velocity of the star, the detection significance is dependent on the vsini prior that is imposed in our modelling. Based on trends observed in the sample of hot Jupiters with obliquity measurements, it has been suggested that these planets modify the spin axes of their host stars, with an efficiency that depends on the stellar type and orbital period of the system. In this framework, short-period planets around stars with surface convective envelopes, like HATS-14b, are expected to have orbits that are aligned with the spin axes of their host stars. HATS-14b, however, is a significant outlier from this trend, challenging the effectiveness of the tidal realignment mechanism.


title: 2-Normed D-Modules and Hahn-Banach Theorem for D-linear 2-functionals

abstract: In this paper we introduce the notion of D-valued 2-norm on hy- perbolic or D-valued modules. Further, we define D-linear 2-functional on these modules and consider some of their properties. We also establish the Hahn- Banach type extension theorem for D-linear 2-functionals.


title: The U-Lagrangian of a prox-regular function

abstract: When restricted to a subspace, a nonsmooth function can be differentiable. It is known that for a nonsmooth convex function f and a point x, the Euclidean space can be decomposed into two subspaces: U, over which a special Lagrangian can be defined and has nice smooth properties and V, the orthogonal complement subspace of U. In this paper we generalize the definition of UV-decomposition and U-Lagrangian to the context of nonconvex functions, specifically that of a prox-regular function.


title: My Reflections on the First Man vs. Machine No-Limit Texas Hold 'em Competition

abstract: The first ever human vs. computer no-limit Texas hold 'em competition took place from April 24-May 8, 2015 at River's Casino in Pittsburgh, PA. In this article I present my thoughts on the competition design, agent architecture, and lessons learned.


title: Initialization-free Distributed Algorithms for Optimal Resource Allocation with Feasibility Constraints and its Application to Economic Dispatch of Power Systems

abstract: In this paper, the distributed resource allocation optimization problem is investigated. The allocation decisions are made to minimize the sum of all the agents' local objective functions while satisfying both the global network resource constraint and the local allocation feasibility constraints. Here the data corresponding to each agent in this separable optimization problem, such as the network resources, the local allocation feasibility constraint, and the local objective function, is only accessible to individual agent and cannot be shared with others, which renders new challenges in this distributed optimization problem. Based on either projection or differentiated projection, two classes of continuous-time algorithms are proposed to solve this distributed optimization problem in an initialization-free and scalable manner. Thus, no re-initialization is required even if the operation environment or network configuration is changed, making it possible to achieve a "plug-and-play" optimal operation of networked heterogeneous agents. The algorithm convergence is guaranteed for strictly convex objective functions, and the exponential convergence is proved for strongly convex functions without local constraints. Then the proposed algorithm is applied to the distributed economic dispatch problem in power grids, to demonstrate how it can achieve the global optimum in a scalable way, even when the generation cost, or system load, or network configuration, is changing.


title: Primal-Dual Algorithm for Distributed Constrained Optimization

abstract: The paper studies a distributed constrained optimization problem, where multiple agents connected in a network collectively minimize the sum of individual objective functions subject to a global constraint being an intersection of the local constraint sets assigned to the agents. Based on the augmented Lagrange method, a distributed primal-dual algorithm with a projection operation included is proposed to solve the problem. It is shown that with appropriately chosen constant step size, the local estimates derived at all agents asymptotically reach a consensus at an optimal solution. In addition, the value of the cost function at the time-averaged estimate converges with rate $O(\frac{1}{k})$ to the optimal value for the unconstrained problem. By these properties the proposed primal-dual algorithm is distinguished from the existing algorithms for distributed constrained optimization. The theoretical analysis is justified by numerical simulations.


title: Composition of topological correspondences

abstract: In the previous article, we proved that a topological correspondence $(X,\lambda)$ from a locally compact groupoid with a Haar system $(G,\alpha)$ to another one, $(H,\beta)$, produces a $C^*$-correspondence $\mathcal{H}(X)$ from $C^*(G,\alpha)$ to $C^*(H,\beta)$. In the present article, we describe how to form a composite of two topological correspondences when the bispaces are Hausdorff and second countable in addition to being locally compact.


title: Kinetic Freeze-out Spectra of Identified Particles Produced in p-Pb Collisions at root(sNN) = 5.02 TeV

abstract: We study the transverse momentum spectra of identified pions, Kaons, protons and Lambdas produced at midrapidity (0 < y < 0.5)in most central p-Pb collisions at root(sNN) = 5.02 TeV in comparison with a Unified Statistical Thermal Freeze-out Model (USTFM). The measurements of pions are reported upto pT = 3 GeV, the kaons (K- + K+) are reported upto pT = 2.5 GeV, K0s is reported upto pT = 7 GeV and the baryons (protons and Lambdas) are reported upto pT = 3.5 GeV. A good agreement is seen between the calculated results and the experimental data points taken from the ALICE experiment. The transverse momentum spectra are found to be flatter for heavy particles than for light particles. Bulk freeze-out properties in terms of kinetic freezeout temperature and the transverse collective flow velocity are extracted from the fits of the transverse momentum spectra of these hadrons. The effect of resonance decay contributions has also been taken care of.


title: Privacy Prediction of Images Shared on Social Media Sites Using Deep Features

abstract: Online image sharing in social media sites such as Facebook, Flickr, and Instagram can lead to unwanted disclosure and privacy violations, when privacy settings are used inappropriately. With the exponential increase in the number of images that are shared online every day, the development of effective and efficient prediction methods for image privacy settings are highly needed. The performance of models critically depends on the choice of the feature representation. In this paper, we present an approach to image privacy prediction that uses deep features and deep image tags as feature representations. Specifically, we explore deep features at various neural network layers and use the top layer (probability) as an auto-annotation mechanism. The results of our experiments show that models trained on the proposed deep features and deep image tags substantially outperform baselines such as those based on SIFT and GIST as well as those that use "bag of tags" as features.


title: Electroweak interaction of particles with accelerated matter and astrophysical applications

abstract: The description of physical processes in accelerated frames opens a window to numerous new phenomena. One can encounter these effects both in the subatomic world and on a macroscale. In the present work we review our recent results on the study of the electroweak interaction of particles with an accelerated background matter. In our analysis we choose the noninertial comoving frame, where matter is at rest. Our study is based on the solution of the Dirac equation, which exactly takes into account both the interaction with matter and the nonintertial effects. First, we study the interaction of ultrarelativistic neutrinos, electrons and quarks with the rotating matter. We consider the influence of the matter rotation on the resonance in neutrino oscillations and the generation of anomalous electric current of charged particles along the rotation axis. Then, we study the creation of neutrino-antineutrino pairs in a linearly accelerated matter. The applications of the obtained results for elementary particle physics and astrophysics are discussed.


title: The impact of the orbital decay of the LAGEOS satellites on the frame-dragging tests

abstract: The laser-tracked geodetic satellites LAGEOS, LAGEOS II and LARES are currently employed, among other things, to measure the general relativistic Lense-Thirring effect in the gravitomagnetic field of the spinning Earth with the hope of providing a more accurate test of such a prediction of the Einstein's theory of gravitation than the existing ones. The secular decay $\dot a$ of the semimajor axes $a$ of such spacecrafts, recently measured in an independent way to a $\sigma_{\dot a}\approx 0.1-0.01$ m yr$^{-1}$ accuracy level, may indirectly impact the proposed relativistic experiment through its connection with the classical orbital precessions induced by the Earth's oblateness $J_2$. \textcolor{black}{Indeed,} the systematic bias due to the current measurement errors $\sigma_{\dot a}$ is of the same order of magnitude of, or even larger than, the expected relativistic signal itself; moreover, it grows linearly with the time span $T$ of the analysis. \textcolor{black}{Therefore, the parameter-fitting algorithms must be properly updated in order to suitably cope with such a new source of systematic uncertainty. Otherwise,} an improvement of one-two orders of magnitude in measuring the orbital decay of the satellites of the LAGEOS family would be required to reduce this source of systematic uncertainty to a percent fraction of the Lense-Thirring signature.

categories: gr-qc astro-ph.EP physics.geo-ph physics.space-ph


title: Anatomy of the Akhmediev breather: cascading instability, first formation time and Fermi-Pasta-Ulam recurrence

abstract: By invoking Bogoliubov's spectrum, we show that for the nonlinear Schrodinger equation, the modulation instability (MI) of its n = 1 Fourier mode on a finite background automatically triggers a further cascading instability, forcing all the higher modes to grow exponentially in locked-step with the n = 1 mode. This fundamental insight, the enslavement of all higher modes to the n = 1 mode, explains the formation of a triangular-shaped spectrum which generates the Akhmediev breather, predicts its formation time analytically from the initial modulation amplitude, and shows that the Fermi-Pasta-Ulam (FPU) recurrence is just a matter of energy conservation with a period twice the breather's formation time. For higher order MI with more than one initial unstable modes, while most evolutions are expected to be chaotic, we show that it is possible to have isolated cases of "super-recurrence", where the FPU period is much longer than that of a single unstable mode.


title: Existence, uniqueness and approximation for $L^p$ solutions of reflected BSDEs under weaker assumptions

abstract: We put forward and prove several existence and uniqueness results for $L^p\ (p>1)$ solutions of reflected BSDEs with continuous barriers and generators satisfying a one-sided Osgood condition together with a general growth condition in $y$ and a uniform continuity condition or a linear growth condition in $z$. A necessary and sufficient condition with respect to the growth of barrier is also explored to ensure the existence of a solution. And, we show that the solutions may be approximated by the penalization method and by some sequences of solutions of reflected BSDEs. Our results improve considerably some known works.


title: Spectrally enhancing near-field radiative transfer between gold gratings by exciting magnetic polariton in nanometric vacuum gaps

abstract: In the present work, we theoretically demonstrate that near field radiative transport between one dimensional periodic grating microstructures separated by nanometer vacuum gaps can be spectrally enhanced by exciting magnetic polariton. Fluctuational electrodynamics that incorporates scattering matrix theory with rigorous coupled wave analysis is employed to exactly calculate the near field radiative flux between two gold gratings. Besides the well known coupled surface plasmon polaritons, the radiative flux can be also spectrally enhanced due to magnetic polariton, which is excited in the gap between gold ridges. The mechanisms of magnetic polariton in the near field radiative transport are elucidated in detail, while the unusual enhancement cannot be predicted by either the Derjaguin or effective medium approximations. The effects of vacuum gap distance and grating geometry parameters between the two gratings are investigated. The findings will open up a new way to control near field radiative transfer by magnetic polariton with micro or nanostructured metamaterials.