title: PPFNet: Global Context Aware Local Features for Robust 3D Point Matching

abstract: We present PPFNet - Point Pair Feature NETwork for deeply learning a globally informed 3D local feature descriptor to find correspondences in unorganized point clouds. PPFNet learns local descriptors on pure geometry and is highly aware of the global context, an important cue in deep learning. Our 3D representation is computed as a collection of point-pair-features combined with the points and normals within a local vicinity. Our permutation invariant network design is inspired by PointNet and sets PPFNet to be ordering-free. As opposed to voxelization, our method is able to consume raw point clouds to exploit the full sparsity. PPFNet uses a novel $\textit{N-tuple}$ loss and architecture injecting the global information naturally into the local descriptor. It shows that context awareness also boosts the local feature representation. Qualitative and quantitative evaluations of our network suggest increased recall, improved robustness and invariance as well as a vital step in the 3D descriptor extraction performance.


title: Stochastic control in microscopic nonequilibrium systems

abstract: Quantifying energy flows at nanometer scales promises to guide future research in a variety of disciplines, from microscopic control and manipulation, to autonomously operating molecular machines. A general understanding of the thermodynamic costs of nonequilibrium processes would illuminate the design principles for efficient microscopic machines. Considerable effort has gone into finding and classifying the deterministic control protocols that drive a system rapidly between states at minimum energetic cost. But for autonomous microscopic systems, driving processes are themselves stochastic. Here we generalize a linear-response framework to incorporate such protocol variability, deriving a lower bound on the work that is realized at finite protocol duration, far from the quasistatic limit. Our findings are confirmed in model systems. This theory provides a thermodynamic rationale for rapid operation, independent of functional incentives.


title: Upper bound for the minimal quantifier depth of the first part of a monadic second-order sentence without asymptotic probability

abstract: In this paper we found an upper bound for the minimal quantifier depth of the first part of a monadic second-order sentence without asymptotic probability described by Jerzy Tyszkiewicz, which express the extension grid axiom in the Erd\H{o}s-R\'enyi model of random graphs $G(n,n^{-\alpha})$ for some irrational $\alpha$.


title: Theory and Ab Initio Computation of the Anisotropic Light Emission in Monolayer Transition Metal Dichalcogenides

abstract: Monolayer transition metal dichalcogenides (TMDCs) are direct gap semiconductors with unique potential for ultrathin light emitters. Yet, their photoluminescence (PL) is not completely understood. We compute the radiative recombination rate in monolayer TMDCs as a function of photon emission direction and polarization, and obtain polar plots of the PL for different excitation scenarios using the ab initio Bethe-Salpeter equation. We show that excitons in a quantum superposition state of the K and K' inequivalent valleys emit light anisotropically upon recombination. Our results can explain the PL anisotropy and polarization dependence measured in recent experiments, and predict new light emission regimes. When averaged over emission angle and exciton momentum, our new treatment recovers the temperature dependent radiative lifetimes we previously derived. Our work provides a first-principles approach to study light emission in two-dimensional materials.


title: Position-Based Multi-Agent Dynamics for Real-Time Crowd Simulation (MiG paper)

abstract: Exploiting the efficiency and stability of Position-Based Dynamics (PBD), we introduce a novel crowd simulation method that runs at interactive rates for hundreds of thousands of agents. Our method enables the detailed modeling of per-agent behavior in a Lagrangian formulation. We model short-range and long-range collision avoidance to simulate both sparse and dense crowds. On the particles representing agents, we formulate a set of positional constraints that can be readily integrated into a standard PBD solver. We augment the tentative particle motions with planning velocities to determine the preferred velocities of agents, and project the positions onto the constraint manifold to eliminate colliding configurations. The local short-range interaction is represented with collision and frictional contact between agents, as in the discrete simulation of granular materials. We incorporate a cohesion model for modeling collective behaviors and propose a new constraint for dealing with potential future collisions. Our new method is suitable for use in interactive games.


title: Efficient collective swimming by harnessing vortices through deep reinforcement learning

abstract: Fish in schooling formations navigate complex flow-fields replete with mechanical energy in the vortex wakes of their companions. Their schooling behaviour has been associated with evolutionary advantages including collective energy savings. How fish harvest energy from their complex fluid environment and the underlying physical mechanisms governing energy-extraction during collective swimming, is still unknown. Here we show that fish can improve their sustained propulsive efficiency by actively following, and judiciously intercepting, vortices in the wake of other swimmers. This swimming strategy leads to collective energy-savings and is revealed through the first ever combination of deep reinforcement learning with high-fidelity flow simulations. We find that a `smart-swimmer' can adapt its position and body deformation to synchronise with the momentum of the oncoming vortices, improving its average swimming-efficiency at no cost to the leader. The results show that fish may harvest energy deposited in vortices produced by their peers, and support the conjecture that swimming in formation is energetically advantageous. Moreover, this study demonstrates that deep reinforcement learning can produce navigation algorithms for complex flow-fields, with promising implications for energy savings in autonomous robotic swarms.


title: Does the black hole shadow probe the event horizon geometry?

abstract: There is an exciting prospect of obtaining the shadow of astrophysical black holes (BHs) in the near future with the Event Horizon Telescope. As a matter of principle, this justifies asking how much one can learn about the BH horizon itself from such a measurement. Since the shadow is determined by a set of special photon orbits, rather than horizon properties, it is possible that different horizon geometries yield similar shadows. One may then ask how sensitive is the shadow to details of the horizon geometry? As a case study, we consider the double Schwarzschild BH and analyse the impact on the lensing and shadows of the conical singularity that holds the two BHs in equilibrium -- herein taken to be a strut along the symmetry axis in between the two BHs. Whereas the conical singularity induces a discontinuity of the scattering angle of photons, clearly visible in the lensing patterns along the direction of the strut's location, it produces no observable effect on the shadows, whose edges remain everywhere smooth. The latter feature is illustrated by examples including both equal and unequal mass BHs. This smoothness contrasts with the intrinsic geometry of the (spatial sections of the) horizon of these BHs, which is not smooth, and provides a sharp example on how BH shadows are insensitive to some horizon geometry details. This observation, moreover, suggests that for the study of their shadows, this static double BH system may be an informative proxy for a dynamical binary.


title: Normal elements of completed group algebras over ${\rm SL}_3(\mathbb{Z}_p) $

abstract: Let $p$ be a prime integer and $\mathbb{Z}_p$ be the ring of $p$-adic integers. By a purely computational approach we prove that each nonzero normal element of a completed group algebra over the special linear group ${\rm SL}_3(\mathbb{Z}_p)$ is a unit. This give a positive answer to an open question in \cite{WeiBian2} and make up for an earlier mistake in \cite{WeiBian1} simultaneously.


title: Clustering Gene Expression Time Series with Coregionalization: Speed propagation of ALS

abstract: Clustering of gene expression time series gives insight into which genes may be coregulated, allowing us to discern the activity of pathways in a given microarray experiment. Of particular interest is how a given group of genes varies with different model conditions or genetic background. Amyotrophic lateral sclerosis (ALS), an irreversible diverse neurodegenerative disorder showed consistent phenotypic differences and the disease progression is heterogeneous with significant variability. This paper demonstrated about finding some significant gene expression profiles and its associated or co-regulated cluster of gene expressions from four groups of data with different genetic background or models conditions. Gene enrichment score analysis and pathway analysis of judicially selected clusters lead toward identifying features underlying the differential speed of disease progression. Gene ontology overrepresentation analysis showed clusters from the proposed method are less likely to be clustered just by chance. In this paper, we develop a new clustering method that allows each cluster to be parameterised according to whether the behaviour of the genes across conditions is correlated or anti-correlated. Our proposed method unveil the potency of latent information shared between multiple model conditions and their replicates during modelling gene expression data.


title: Biological Mechanisms for Learning: A Computational Model of Olfactory Learning in the Manduca sexta Moth, with Applications to Neural Nets

abstract: The insect olfactory system, which includes the antennal lobe (AL), mushroom body (MB), and ancillary structures, is a relatively simple neural system capable of learning. Its structural features, which are widespread in biological neural systems, process olfactory stimuli through a cascade of networks where large dimension shifts occur from stage to stage and where sparsity and randomness play a critical role in coding. Learning is partly enabled by a neuromodulatory reward mechanism of octopamine stimulation of the AL, whose increased activity induces rewiring of the MB through Hebbian plasticity. Enforced sparsity in the MB focuses Hebbian growth on neurons that are the most important for the representation of the learned odor. Based upon current biophysical knowledge, we have constructed an end-to-end computational model of the Manduca sexta moth olfactory system which includes the interaction of the AL and MB under octopamine stimulation. Our model is able to robustly learn new odors, and our simulations of integrate-and-fire neurons match the statistical features of in-vivo firing rate data. From a biological perspective, the model provides a valuable tool for examining the role of neuromodulators, like octopamine, in learning, and gives insight into critical interactions between sparsity, Hebbian growth, and stimulation during learning. Our simulations also inform predictions about structural details of the olfactory system that are not currently well-characterized. From a machine learning perspective, the model yields bio-inspired mechanisms that are potentially useful in constructing neural nets for rapid learning from very few samples. These mechanisms include high-noise layers, sparse layers as noise filters, and a biologically-plausible optimization method to train the network based on octopamine stimulation, sparse layers, and Hebbian growth.


title: A Semi-Supervised Two-Stage Approach to Learning from Noisy Labels

abstract: The recent success of deep neural networks is powered in part by large-scale well-labeled training data. However, it is a daunting task to laboriously annotate an ImageNet-like dateset. On the contrary, it is fairly convenient, fast, and cheap to collect training images from the Web along with their noisy labels. This signifies the need of alternative approaches to training deep neural networks using such noisy labels. Existing methods tackling this problem either try to identify and correct the wrong labels or reweigh the data terms in the loss function according to the inferred noisy rates. Both strategies inevitably incur errors for some of the data points. In this paper, we contend that it is actually better to ignore the labels of some of the data points than to keep them if the labels are incorrect, especially when the noisy rate is high. After all, the wrong labels could mislead a neural network to a bad local optimum. We suggest a two-stage framework for the learning from noisy labels. In the first stage, we identify a small portion of images from the noisy training set of which the labels are correct with a high probability. The noisy labels of the other images are ignored. In the second stage, we train a deep neural network in a semi-supervised manner. This framework effectively takes advantage of the whole training set and yet only a portion of its labels that are most likely correct. Experiments on three datasets verify the effectiveness of our approach especially when the noisy rate is high.


title: On the inverse problem of source reconstruction from coherence measurements

abstract: We consider an inverse source problem for partially coherent light propagating in the Fresnel regime. The data is the coherence of the field measured away from the source. The reconstruction is based on a minimum residue formulation, which uses the authors' recent closed-form approximation formula for the coherence of the propagated field. The developed algorithms require a small data sample for convergence and yield stable inversion by exploiting information in the coherence as opposed to intensity-only measurements. Examples with both simulated and experimental data demonstrate the ability of the proposed approach to simultaneously recover complex sources in different planes transverse to the direction of propagation.


title: Towards A Systems Approach To Distributed Programming

abstract: It is undeniable that most developers today are building distributed applications. However, most of these applications are developed by composing existing systems together through unspecified APIs exposed to the application developer. Systems are not going away: they solve a particular problem and most applications today need to rely on several of these systems working in concert. Given this, we propose a research direction where higher-level languages with well defined semantics target underlying systems infrastructure as a middle-ground.


title: A diffusion generated method for computing Dirichlet partitions

abstract: A Dirichlet $k$-partition of a closed $d$-dimensional surface is a collection of $k$ pairwise disjoint open subsets such that the sum of their first Laplace-Beltrami-Dirichlet eigenvalues is minimal. In this paper, we develop a simple and efficient diffusion generated method to compute Dirichlet $k$-partitions for $d$-dimensional flat tori and spheres. For the $2d$ flat torus, for most values of $k=3$-9,11,12,15,16, and 20, we obtain hexagonal honeycombs. For the $3d$ flat torus and $k=2,4,8,16$, we obtain the rhombic dodecahedral honeycomb, the Weaire-Phelan honeycomb, and Kelvin's tessellation by truncated octahedra. For the $4d$ flat torus, for $k=4$, we obtain a constant extension of the rhombic dodecahedral honeycomb along the fourth direction and for $k=8$, we obtain a 24-cell honeycomb. For the $2d$ sphere, we also compute Dirichlet partitions for $k=3$-7,9,10,12,14,20. Our computational results agree with previous studies when a comparison is available. As far as we are aware, these are the first published results for Dirichlet partitions of the $4d$ flat torus.


title: Prediction of Shared Bicycle Demand with Wavelet Thresholding

abstract: Consumers are creatures of habit, often periodic, tied to work, shopping and other schedules. We analyzed one month of data from the world's largest bike-sharing company to elicit demand behavioral cycles, initially using models from animal tracking that showed large customers fit an Ornstein-Uhlenbeck model with demand peaks at periodicities of 7, 12, 24 hour and 7-days. Lorenz curves of bicycle demand showed that the majority of customer usage was infrequent, and demand cycles from time-series models would strongly overfit the data yielding unreliable models. Analysis of thresholded wavelets for the space-time tensor of bike-sharing contracts was able to compress the data into a 56-coefficient model with little loss of information, suggesting that bike-sharing demand behavior is exceptionally strong and regular. Improvements to predicted demand could be made by adjusting for 'noise' filtered by our model from air quality and weather information and demand from infrequent riders.


title: Gaussian binomial coefficients with negative arguments

abstract: Loeb showed that a natural extension of the usual binomial coefficient to negative (integer) entries continues to satisfy many of the fundamental properties. In particular, he gave a uniform binomial theorem as well as a combinatorial interpretation in terms of choosing subsets of sets with a negative number of elements. We show that all of this can be extended to the case of Gaussian binomial coefficients. Moreover, we demonstrate that several of the well-known arithmetic properties of binomial coefficients also hold in the case of negative entries. In particular, we show that Lucas' Theorem on binomial coefficients modulo $p$ not only extends naturally to the case of negative entries, but even to the Gaussian case.


title: Stubborn Transaction Reduction (with Proofs)

abstract: The exponential explosion of parallel interleavings remains a fundamental challenge to model checking of concurrent programs. Both partial-order reduction (POR) and transaction reduction (TR) decrease the number of interleavings in a concurrent system. Unlike POR, transactions also reduce the number of intermediate states. Modern POR techniques, on the other hand, offer more dynamic ways of identifying commutative behavior, a crucial task for obtaining good reductions. We show that transaction reduction can use the same dynamic commutativity as found in stubborn set POR. We also compare reductions obtained by POR and TR, demonstrating with several examples that these techniques complement each other. With an implementation of the dynamic transactions in the model checker LTSmin, we compare its effectiveness with the original static TR and two POR approaches. Several inputs, including realistic case studies, demonstrate that the new dynamic TR can surpass POR in practice.


title: Thermodynamically Favorable Computation via Tile Self-assembly

abstract: The recently introduced Thermodynamic Binding Networks (TBN) model was developed with the purpose of studying self-assembling systems by focusing on their thermodynamically favorable final states, and ignoring the kinetic pathways through which they evolve. The model was intentionally developed to abstract away not only the notion of time, but also the constraints of geometry. Collections of monomers with binding domains which allow them to form polymers via complementary bonds are analyzed to determine their final, stable configurations, which are those which maximize the number of bonds formed (i.e. enthalpy) and the number of independent components (i.e. entropy). In this paper, we first develop a definition of what it means for a TBN to perform a computation, and then present a set of constructions which are capable of performing computations by simulating the behaviors of space-bounded Turing machines and boolean circuits. In contrast to previous TBN results, these constructions are robust to great variability in the counts of monomers existing in the systems and the numbers of polymers that form in parallel. Although the Turing machine simulating TBNs are inefficient in terms of the numbers of unique monomer types required, as compared to algorithmic self-assembling systems in the abstract Tile Assembly Model (aTAM), we then show that a general strategy of porting those aTAM system designs to TBNs produces TBNs which incorrectly simulate computations. Finally, we present a refinement of the TBN model which we call the Geometric Thermodynamic Binding Networks (GTBN) model in which monomers are defined with rigid geometries and form rigid bonds. Utilizing the constraints imposed by geometry, we then provide a GTBN construction capable of simulating Turing machines as efficiently as in the aTAM.


title: Stability of Kramers Majorana doublets: the effects of interactions and disorders

abstract: In this work we study the effects of interactions and disorder on $1D$ DIII topological superconductors and the Majorana Kramers doublets (MKDs). In contract to the case without the time-reversal symmetry, the Umklapp interaction plays important roles in this system. The underlying phases due to the Umklapp interaction and disorder are found by using a perturbative renormaliztion analysis based on the Abelian Bosonization. Importantly, the stable topological regime can be found within a rather wide parameter space. Furthermore, the degeneracy splitting of the MKDs is shown to be still exponentially dependent on the length of the wire in the presence of both the Umklapp interaction and disorder, when the Luttinger parameter $K_0>\sqrt{2}/2$. The differences caused by the Umklapp interaction are highlighted in contrast to the time-reversal breaking cases.


title: Exact Semidefinite Formulations for a Class of (Random and Non-Random) Nonconvex Quadratic Programs

abstract: We study a class of quadratically constrained quadratic programs (QCQPs), called {\em diagonal QCQPs\/}, which contain no off-diagonal terms $x_j x_k$ for $j \ne k$, and we provide a sufficient condition on the problem data guaranteeing that the basic Shor semidefinite relaxation is exact. Our condition complements and refines those already present in the literature and can be checked in polynomial time. We then extend our analysis from diagonal QCQPs to general QCQPs, i.e., ones with no particular structure. By reformulating a general QCQP into diagonal form, we establish new, polynomial-time-checkable sufficient conditions for the semidefinite relaxations of general QCQPs to be exact. Finally, these ideas are extended to show that a class of random general QCQPs has exact semidefinite relaxations with high probability as long as the number of constraints grows no faster than a fixed polynomial in the number of variables. To the best of our knowledge, this is the first result establishing the exactness of the semidefinite relaxation for random general QCQPs.


title: Digital Data Archives as Knowledge Infrastructures: Mediating Data Sharing and Reuse

abstract: Digital archives are the preferred means for open access to research data. They play essential roles in knowledge infrastructures - robust networks of people, artifacts, and institutions - but little is known about how they mediate information exchange between stakeholders. We open the "black box" of data archives by studying DANS, the Data Archiving and Networked Services institute of The Netherlands, which manages 50+ years of data from the social sciences, humanities, and other domains. Our interviews, weblogs, ethnography, and document analyses reveal that a few large contributors provide a steady flow of content, but most are academic researchers who submit datasets infrequently and often restrict access to their files. Consumers are a diverse group that overlaps minimally with contributors. Archivists devote about half their time to aiding contributors with curation processes and half to assisting consumers. Given the diversity and infrequency of usage, human assistance in curation and search remains essential. DANS' knowledge infrastructure encompasses public and private stakeholders who contribute, consume, harvest, and serve their data - many of whom did not exist at the time the DANS collections originated - reinforcing the need for continuous investment in digital data archives as their communities, technologies, and services evolve.


title: Driver Gaze Zone Estimation using Convolutional Neural Networks: A General Framework and Ablative Analysis

abstract: Driver gaze has been shown to be an excellent surrogate for driver attention in intelligent vehicles. With the recent surge of highly autonomous vehicles, driver gaze can be useful for determining the handoff time to a human driver. While there has been significant improvement in personalized driver gaze zone estimation systems, a generalized system which is invariant to different subjects, perspectives and scales is still lacking. We take a step towards this generalized system using Convolutional Neural Networks (CNNs). We finetune 4 popular CNN architectures for this task, and provide extensive comparisons of their outputs. We additionally experiment with different input image patches, and also examine how image size affects performance. For training and testing the networks, we collect a large naturalistic driving dataset comprising of 11 long drives, driven by 10 subjects in two different cars. Our best performing model achieves an accuracy of 95.18% during cross-subject testing, outperforming current state of the art techniques for this task. Finally, we evaluate our best performing model on the publicly available Columbia Gaze Dataset comprising of images from 56 subjects with varying head pose and gaze directions. Without any training, our model successfully encodes the different gaze directions on this diverse dataset, demonstrating good generalization capabilities.


title: A Bayesian Approach to Multi-State Hidden Markov Models: Application to Dementia Progression

abstract: People are living longer than ever before, and with this arises new complications and challenges for humanity. Among the most pressing of these challenges is of understanding the role of aging in the development of dementia. This paper is motivated by the Mayo Clinic Study of Aging data for 4742 subjects since 2004, and how it can be used to draw inference on the role of aging in the development of dementia. We construct a hidden Markov model (HMM) to represent progression of dementia from states associated with the buildup of amyloid plaque in the brain, and the loss of cortical thickness. A hierarchical Bayesian approach is taken to estimate the parameters of the HMM with a truly time-inhomogeneous infinitesimal generator matrix, and response functions of the continuous-valued biomarker measurements are cut-point agnostic. A Bayesian approach with these features could be useful in many disease progression models. Additionally, an approach is illustrated for correcting a common bias in delayed enrollment studies, in which some or all subjects are not observed at baseline. Standard software is incapable of accounting for this critical feature, so code to perform the estimation of the model described below is made available online.


title: Longitudinal structure of the photospheric magnetic field in the Carrington system

abstract: The observations of the Sun have been performed over the years, even centuries- Whether are there active longitudes? If yes how stable are they? One of the first The Wilcox Solar Observatory data taken over three cycles N 21, N 22, N 23 have been used to reveal the longitudinal structure of the photospheric magnetic field. Mean over three cycles magnetic field distribution has been calculated in the North and in the South hemispheres as well as at 30 levels of latitude from -75 to 75 degrees. This study was performed using observations of the magnetic field taking into account its polarity or only intensity. The longitudinal structure of the magnetic field was calculated in the coordinate system rotating with Carrington rate. These structures were compared with a model of random longitudinal distribution of the magnetic field. Random character of the longitudinal structure was refused. The results agree with the presence of two active meridians seen in different phenomena of solar activity at longitudes separated by 150-170 degrees in the Carrington coordinate system.


title: Gamification: a Game Changer for Managing Technical Debt? A Design Study

abstract: Context: Technical debt management is challenging for software engineers due to poor tool support and a lack of knowledge on how to prioritize technical debt repayment and prevention activities. Furthermore, when there is a large backlog of debt, developers often lack the motivation to address it. Objective: In this paper, we describe a design study to investigate how gamification can support Technical Debt Management in a large legacy software system of an industrial company. Our study leads to a novel tool (named Themis) that combines technical debt support, version control, and gamification features. In addition to gamification features, Themis provides suggestions for developers on where to focus their effort, and visualizations for managers to track technical debt activities. Method: We describe how Themis was refined and validated in an iterative deployment with the company, finally conducting a qualitative study to investigate how the features of Themis affect technical debt management behavior. We consider the impact on both developers and managers. Results: Our results show that it achieves increased developer motivation, and supports managers in monitoring and influencing developer behaviors. We show how our findings may be transferable to other contexts by proposing guidelines on how to apply gamification. Conclusions: With this case, gamification appears as a promising solution to help technical debt management, although it needs to be carefully designed and implemented to avoid its possible negative effects.


title: Monotone Operator Theory in Convex Optimization

abstract: Several aspects of the interplay between monotone operator theory and convex optimization are presented. The crucial role played by monotone operators in the analysis and the numerical solution of convex minimization problems is emphasized. We review the properties of subdifferentials as maximally monotone operators and, in tandem, investigate those of proximity operators as resolvents. In particular, we study new transformations which map proximity operators to proximity operators, and establish connections with self-dual classes of firmly nonexpansive operators. In addition, new insights and developments are proposed on the algorithmic front.


title: The physics of climate change: simple models in climate science

abstract: There is a perception that climate science can only be approached with complex computer simulations. But working climate scientists often use simple models to understand their simulations and make order-of-magnitude estimates. This article presents some of these simple models with the goal of making climate science more accessible and comprehensible.


title: Improving the Universality and Learnability of Neural Programmer-Interpreters with Combinator Abstraction

abstract: To overcome the limitations of Neural Programmer-Interpreters (NPI) in its universality and learnability, we propose the incorporation of combinator abstraction into neural programing and a new NPI architecture to support this abstraction, which we call Combinatory Neural Programmer-Interpreter (CNPI). Combinator abstraction dramatically reduces the number and complexity of programs that need to be interpreted by the core controller of CNPI, while still allowing the CNPI to represent and interpret arbitrary complex programs by the collaboration of the core with the other components. We propose a small set of four combinators to capture the most pervasive programming patterns. Due to the finiteness and simplicity of this combinator set and the offloading of some burden of interpretation from the core, we are able construct a CNPI that is universal with respect to the set of all combinatorizable programs, which is adequate for solving most algorithmic tasks. Moreover, besides supervised training on execution traces, CNPI can be trained by policy gradient reinforcement learning with appropriately designed curricula.


title: 3D Meteoroid Trajectories

abstract: Meteoroid modelling of fireball data typically uses a one dimensional model along a straight line triangulated trajectory. The assumption of a straight line trajectory has been considered an acceptable simplification for fireballs, but it has not been rigorously tested. The unique capability of the Desert Fireball Network (DFN) to triangulate discrete observation times gives the opportunity to investigate the deviation of a meteoroid's position to different model fits. Here we assess the viability of a straight line assumption for fireball data in two meteorite-dropping test cases observed by the Desert Fireball Network (DFN) in Australia -- one over 21 seconds (\textit{DN151212\_03}), one under 5 seconds (\textit{DN160410\_03}). We show that a straight line is not valid for these two meteorite dropping events and propose a three dimensional particle filter to model meteoroid positions without any straight line constraints. The single body equations in three dimensions, along with the luminosity equation, are applied to the particle filter methodology described by \citet{Sansom2017}. Modelling fireball camera network data in three dimensions has not previously been attempted. This allows the raw astrometric, line-of-sight observations to be incorporated directly. In analysing these two DFN events, the triangulated positions based on a straight line assumption result in the modelled meteoroid positions diverging up to $3.09\, km$ from the calculated observed point (for \textit{DN151212\_03}). Even for the more typical fireball event, \textit{DN160410\_03}, we see a divergence of up to $360$\,m. As DFN observations are typically precise to $<100$\,m, it is apparent that the assumption of a straight line is an oversimplification that will affect orbit calculations and meteorite search regions for a significant fraction of events.


title: More Efficient Estimation for Logistic Regression with Optimal Subsample

abstract: In this paper, we propose improved estimation method for logistic regression based on subsamples taken according the optimal subsampling probabilities developed in Wang et al. 2018 Both asymptotic results and numerical results show that the new estimator has a higher estimation efficiency. We also develop a new algorithm based on Poisson subsampling, which does not require to approximate the optimal subsampling probabilities all at once. This is computationally advantageous when available random-access memory is not enough to hold the full data. Interestingly, asymptotic distributions also show that Poisson subsampling produces a more efficient estimator if the sampling rate, the ratio of the subsample size to the full data sample size, does not converge to zero. We also obtain the unconditional asymptotic distribution for the estimator based on Poisson subsampling. The proposed approach requires to use a pilot estimator to correct biases of un-weighted estimators. We further show that even if the pilot estimator is inconsistent, the resulting estimators are still consistent and asymptotically normal if the model is correctly specified.


title: Immediate Causality Network of Stock Markets

abstract: A financial system contains many elements networked by their relationships. Extensive works show that topological structure of the network stores rich information on evolutionary behaviors of the system such as early warning signals of collapses and/or crises. Existing works focus mainly on the network structure within a single stock market, while a collapse/crisis occurs in a macro-scale covering several or even all markets in the world. This mismatch of scale leads to unacceptable noise to the topological structure, and lack of information stored in relationships between different markets. In this work by using the transfer entropy we reconstruct the influential network between ten typical stock markets distributed in the world. Interesting findings include, before a financial crisis the connection strength reaches a maxima, which can act as an early warning signal of financial crises; The markets in America are mono-directionally and strongly influenced by that in Europe and act as the center; Some strongly linked pairs have also close correlations. The findings are helpful in understanding the evolution and modelling the dynamical process of the global financial system.


title: ODINI : Escaping Sensitive Data from Faraday-Caged, Air-Gapped Computers via Magnetic Fields

abstract: Air-gapped computers are computers which are kept isolated from the Internet, because they store and process sensitive information. When highly sensitive data is involved, an air-gapped computer might also be kept secluded in a Faraday cage. The Faraday cage prevents the leakage of electromagnetic signals emanating from various computer parts, which may be picked up by an eavesdropping adversary remotely. The air-gap separation, coupled with the Faraday shield, provides a high level of isolation, preventing the potential leakage of sensitive data from the system. In this paper, we show how attackers can bypass Faraday cages and air-gaps in order to leak data from highly secure computers. Our method is based on an exploitation of the magnetic field generated by the computer CPU. Unlike electromagnetic radiation (EMR), low frequency magnetic radiation propagates though the air, penetrating metal shielding such as Faraday cages (e.g., compass still works inside Faraday cages). We introduce a malware code-named ODINI that can control the low frequency magnetic fields emitted from the infected computer by regulating the load of the CPU cores. Arbitrary data can be modulated and transmitted on top of the magnetic emission and received by a magnetic receiver (bug) placed nearby. We provide technical background and examine the characteristics of the magnetic fields. We implement a malware prototype and discuss the design considerations along with the implementation details. We also show that the malicious code does not require special privileges (e.g., root) and can successfully operate from within isolated virtual machines (VMs) as well.


title: Restate the reference for EEG microstate analysis

abstract: Despite the decades of efforts, the choice of EEG reference is still a debated fundamental issue. Non-neutral reference can inevitably inject the uncontrolled temporal biases into all EEG recordings, which may influence the spatiotemporal analysis of brain activity. A method, termed microstates, identifying spatiotemporal EEG features as the quasi-stable topography states in milliseconds, suggests its potential as biomarkers of neurophysiological disease. As reference electrode standardization technique (REST) could reconstruct an infinity reference approximately, it is a question whether REST or the other references will be more reliable than average reference (AR) for the microstates analysis. In this study, we design the microstate-based EEG forward model, and apply different references for microstates analysis. The spatial similarity between the generated and assumed cluster maps is mainly investigated. Furthermore, the real EEG data by the parametric bootstrap method is used to validate the performance of the references. Finally, we find that REST is robust to recover more similar cluster maps to the assumption than AR in the simulation, and the cluster maps between REST and AR on the real EEG data are quite different. This study may indicate that REST contributes to identifying more objective microstates features than AR.


title: On Quantizer Design to Exploit Common Information in Layered Coding of Vector Sources

abstract: This paper studies a layered coding framework with a relaxed hierarchical structure. Advances in wired/wireless communication and consumer electronic devices have created a requirement for serving the same content at different quality levels. The key challenge is to optimally encode all the required quality levels with efficient usage of storage and networking resources. The approach to store and transmit independent copies for every required quality level is highly wasteful in resources. Alternatively, conventional scalable coding has inherent loss due to its structure. This paper studies a layered coding framework with a relaxed hierarchical structure to transmit information common to different quality levels along with individual bit streams for each quality level. The flexibility of sharing only a properly selected subset of information from a lower quality level with the higher quality level, enables achieving operating points between conventional scalable coding and independent coding, to control the layered coding penalty. Jointly designing common and individual layers' coders overcomes the limitations of conventional scalable coding and non-scalable coding, by providing the flexibility of transmitting common and individual bit-streams for different quality levels. It extracts the common information between different quality levels with negligible performance penalty. Simulation results for practically important sources, confirm the superiority of the work.


title: A strongly interacting Sarma superfluid near orbital Feshbach resonances

abstract: We investigate the nature of superfluid pairing in a strongly interacting Fermi gas near orbital Feshbach resonances with spin-population imbalance in three dimensions, which can be well described by a two-band or two-channel model. We show that a Sarma superfluid with gapless single-particle excitations is favored in the closed channel at large imbalance. It is thermodynamically stable against the formation of an inhomogeneous Fulde\textendash Ferrell\textendash Larkin\textendash Ovchinnikov superfluid and features a well-defined Goldstone-Anderson-Bogoliubov phonon mode and a massive Leggett mode as collective excitations at low momentum. At large momentum, the Leggett mode disappears and the phonon mode becomes damped at zero temperature, due to the coupling to the particle-hole excitations. We discuss possible experimental observation of a strongly interacting Sarma superfluid with ultracold alkaline-earth-metal Fermi gases.


title: Folding mechanisms at finite temperature

abstract: Folding mechanisms are zero elastic energy motions essential to the deployment of origami, linkages, reconfigurable metamaterials and robotic structures. In this paper, we determine the fate of folding mechanisms when such structures are miniaturized so that thermal fluctuations cannot be neglected. First, we identify geometric and topological design strategies aimed at minimizing undesired thermal energy barriers that generically obstruct kinematic mechanisms at the microscale. Our findings are illustrated in the context of a quasi one-dimensional linkage structure that harbors a topologically protected mechanism. However, thermal fluctuations can also be exploited to deliberately lock a reconfigurable metamaterial into a fully expanded configuration, a process reminiscent of order by disorder transitions in magnetic systems. We demonstrate that this effect leads certain topological mechanical structures to exhibit an abrupt change in the pressure -- a bulk signature of the underlying topological invariant at finite temperature. We conclude with a discussion of anharmonic corrections and potential applications of our work to the the engineering of DNA origami devices and molecular robots.


title: Frobenius powers

abstract: This article extends the notion of a Frobenius power of an ideal in prime characteristic to allow arbitrary nonnegative real exponents. These generalized Frobenius powers are closely related to test ideals in prime characteristic, and multiplier ideals over fields of characteristic zero. For instance, like these well-known families of ideals, Frobenius powers also give rise to jumping exponents that we call critical Frobenius exponents. In fact, the Frobenius powers of a principal ideal coincides with its test ideals, but appear to be a more refined measure of singularities in general. Herein, we develop the theory of Frobenius powers in regular domains, and apply it to study singularities, especially those of generic hypersurfaces. These applications illustrate one way in which multiplier ideals behave more like Frobenius powers than like test ideals.


title: Coded Caching with Heterogeneous Cache Sizes and Link Qualities: The Two-User Case

abstract: Centralized coded caching problem is studied for the two-user scenario, considering heterogeneous cache capacities at the users and private channels from the server to the users, in addition to a shared channel. Optimal caching and delivery strategies that minimize the worst-case delivery latency are presented for an arbitrary number of files. The converse proof follows from the sufficiency of file-index-symmetric caching and delivery codes, while the achievability is obtained through memory-sharing among a number of special memory capacity pairs. The optimal scheme is shown to exploit the private link capacities by transmitting part of the corresponding user`s request in an uncoded fashion. When there are no private links, the results presented here improve upon the two known results in the literature, namely, i) equal cache capacities and arbitrary number of files; and ii) unequal cache capacities and $N=2$ files. The results are then extended to the caching problem with heterogeneous distortion requirements.


title: An interview based study of pioneering experiences in teaching and learning Complex Systems in Higher Education

abstract: Due to the interdisciplinary nature of complex systems as a field, students studying complex systems at University level have diverse disciplinary backgrounds. This brings challenges (e.g. wide range of computer programming skills) but also opportunities (e.g. facilitating interdisciplinary interactions and projects) for the classroom. However, there is little published regarding how these challenges and opportunities are handled in teaching and learning Complex Systems as an explicit subject in higher education, and how this differs in comparison to other subject areas. We seek to explore these particular challenges and opportunities via an interview-based study of pioneering teachers and learners (conducted amongst the authors) regarding their experiences. We compare and contrast those experiences, and analyse them with respect to the educational literature. Our discussions explored: approaches to curriculum design, how theories/models/frameworks of teaching and learning informed decisions and experience, how diversity in student backgrounds was addressed, and assessment task design. We found a striking level of commonality in the issues expressed as well as the strategies to handle them, for example a significant focus on problem-based learning, and the use of major student-led creative projects for both achieving and assessing learning outcomes.


title: Some remarks on the non-real roots of polynomials

abstract: Let $f \in { \mathbb R} ( t) [x]$ be given by $ f(t, x) = x^n + t \cdot g(x) $ and $\beta_1 < \dots < \beta_m$ the distinct real roots of the discriminant $\Delta_{(f, x)} (t)$ of $f(t, x)$ with respect to $x$. Let $\gamma$ be the number of real roots of $g(x)=\sum_{k=0}^s t_{s-k} x^{s-k}$. For any $\xi > | \beta_m |$, if $n-s$ is odd then the number of real roots of $f(\xi, x)$ is $\gamma+1$, and if $n-s$ is even then the number of real roots of $f(\xi, x)$ is $\gamma$, $\gamma+2$ if $t_s>0$ or $t_s < 0$ respectively. A special case of the above result is constructing a family of totally complex polynomials which are reducible over $\mathbb Q$.


title: Layered Coding of Hidden Markov Sources

abstract: The paper studies optimal coding of hidden Markov sources (HMS), which represent a broad class of practical sources obtained through noisy acquisition processes, beside their explicit modeling use in speech processing and recognition, image understanding and sensor networks. A new fundamental source coding approach for HMS is proposed, based on tracking an estimate of the state probability distribution, and is shown to be optimal. Practical encoder and decoder schemes that leverage the main concepts are introduced. An iterative approach is developed for optimizing the system. It also focuses on a significant extension of the optimal HMS quantization paradigm. It proposes a new approach for scalable coding of HMS which accounts for all the available information while coding a given layer. Simulation results confirm that these approaches significantly reduce the reconstructed distortion and substantially outperform existing techniques.


title: A low phase noise microwave frequency synthesizer based on parameters optimized NLTL for Cs fountain clock

abstract: We report on the development and phase noise performance of a 9.1926 GHz microwave frequency synthesizer to be used as the local oscillator for a Cs fountain clock. It is based on frequency multiplication and synthesis from an ultralow phase noise 5 MHz Oven Controlled Crystal Oscillator (OCXO) and 100 MHz Voltage Controlled Crystal Oscillator (VCXO).The key component of the frequency multiplication is a non-linear transmission-line (NLTL) used as a frequency comb generator. The phase noise of the synthesizer is improved by carefully optimizing the input power, the input and output impedances of the NLTL. The absolute phase noises of the 9.1926 GHz output signal are measured to be -64 dBc/Hz, -83 dBc/Hz, -92 dBc/Hz, -117 dBc/Hz and -119 dBc/Hz at 1 Hz, 10Hz, 100Hz, 1 kHz and 10 kHz offset frequencies, respectively. The residual phase noise of the synthesizer is measured to be -82 dBc/Hz at 1 Hz offset frequency. The measurement result shows that the absolute phase noise at the frequency range of 1 - 100 Hz is mainly limited by the phase noise of the OCXO. The contribution of the absolute phase noise to the fountain clock short-term frequency stability is calculated to be 7.0x10^(-14). The residual frequency stability of the synthesizer is measured to be1.5x10^(-14), which is consistent with the calculated frequency stability due to the residual phase noise of the synthesizer. Meanwhile we designed and realized an interferometric microwave switch in the synthesizer to eliminate the frequency shifts induced by the microwave leakage. The extinction ratio of the switch is measured to be more than 50 dB. In the scheme, we use only commercially available components to build the microwave frequency synthesizer with excellent phase noise performance for high-performance Cs fountain clocks.


title: Visualizing Heavy Fermion Confinement and Pauli-Limited Superconductivity in Layered CeCoIn5

abstract: Layered material structures play a key role in enhancing electron-electron interactions to create correlated metallic phases that can transform into unconventional superconducting states. The quasi-two-dimensional electronic properties of such compounds are often inferred indirectly through examination of their bulk properties. Here we use scanning tunneling microscopy and spectroscopy to directly probe in cross section the quasi-two-dimensional correlated electronic states of the heavy fermion superconductor CeCoIn5. Our measurements reveal the strong confined nature of heavy quasi-particles, anisotropy of tunneling characteristics, and layer-by-layer modulated behavior of the precursor pseudogap gap phase in this compound. Examining the interlayer coupled superconducting state at low temperatures, we find that the orientation of line defects relative to the d-wave order parameter determines whether in-gap states form due to scattering. Spectroscopic imaging of the anisotropic magnetic vortex cores directly characterizes the short interlayer superconducting coherence length and shows an electronic phase separation near the upper critical in-plane magnetic field, consistent with a Pauli-limited first-order phase transition into a pseudogap phase.


title: New Equation For Describing Time Dependence of Moon Orbit Radius

abstract: I derived an equation to describe the dependence of Moon orbit radius around the Earth. I obtained the radius changes with time according to a scaling equation, Moon orbit radius is proportional to power 2/5 of the Moon age. Using the equation I have been able to predict accurately some quantities that have been well accepted


title: Direct observation of spin excitation anisotropy in the paramagnetic orthorhombic state of BaFe$_{2-x}$Ni$_x$As$_2$

abstract: We use transport and inelastic neutron scattering measurements to investigate single crystals of iron pnictide BaFe$_{2-x}$Ni$_{x}$As$_{2}$ ($x=0,0.03$), which exhibit a tetragonal-to-orthorhombic structural transition at $T_s$ and stripe antiferromagnetic order at $T_N$ ($T_s\geq T_N$). Using a tunable uniaxial pressure device, we detwin the crystals and study their transport and spin excitation properties at antiferromagnetic wave vector $S_1(1,0)$ and its 90$^\circ$ rotated wave vector $S_2(0,1)$ under different pressure conditions. We find that uniaxial pressure necessary to detwin and maintain single domain orthorhombic antiferromagnetic phase of BaFe$_{2-x}$Ni$_{x}$As$_{2}$ induces resistivity and spin excitation anisotropy at temperatures above zero pressure $T_s$. In uniaxial pressure-free detwinned sample, spin excitation anisotropy between $S_1(1,0)$ and $S_2(0,1)$ first appear in the paramagnetic orthorhombic phase below $T_s$. These results are consistent with predictions of spin nematic theory, suggesting the absence of structural or nematic phase transition above $T_s$ in iron pnictides.


title: Earth-like and Tardigrade survey of exoplanets

abstract: Finding life on other worlds is a fascinating area of astrobiology and planetary sciences. Presently, over 3500 exoplanets, representing a very wide range of physical and chemical environments, are known. Tardigrades (water bears) are microscopic invertebrates that inhabit almost all terrestrial, freshwater and marine habitats, from the highest mountains to the deepest oceans. Thanks to their ability to live in a state of cryptobiosis, which is known to be an adaptation to unpredictably fluctuating environmental conditions, these organisms are able to survive when conditions are not suitable for active life; consequently, tardigrades are known as the toughest animals on Earth. In their cryptobiotic state, they can survive extreme conditions, such as temperatures below -250{\deg}C and up to 150{\deg}C, high doses of ultraviolet and ionising radiation, up to 30 years without liquid water, low and high atmospheric pressure, and exposure to many toxic chemicals. Active tardigrades are also resistant to a wide range of unfavourable environmental conditions, which makes them an excellent model organism for astrobiological studies. In our study, we have established a metric tool for distinguishing the potential survivability of active and cryptobiotic tardigrades on rocky-water and water-gas planets in our solar system and exoplanets, taking into consideration the geometrical means of surface temperature and surface pressure of the considered planets. The Active Tardigrade Index (ATI) and Cryobiotic Tardigrade Index (CTI) are two metric indices with minimum value 0 (= tardigrades cannot survive) and maximum 1 (= tardigrades will survive in their respective state). Values between 0 and 1 indicate a percentage chance of the active or cryptobiotic tardigrades surviving on a given exoplanet.


title: Gromov-Hyperbolicity of the ray graph and quasimorphisms on a big mapping class group

abstract: These notes are the English version of the paper "Hyperbolicit\'e du graphe des rayons et quasi-morphismes sur un gros groupe modulaire". The mapping class group Gamma of the complement of a Cantor set in the plane arises naturally in dynamics. We show that the ray graph, which is the analog of the complex of curves for this surface of infinite type, has infinite diameter and is hyperbolic. We use the action of Gamma on this graph to find an explicit non trivial quasimorphism on Gamma and to show that this group has infinite dimensional second bounded cohomology. Finally we give an example of a hyperbolic element of Gamma with vanishing stable commutator length. This carries out a program proposed by Danny Calegari.


title: On the small-$x$ behavior of the orbital angular momentum distributions in QCD

abstract: We present the numerical solution of the leading order QCD evolution equation for the orbital angular momentum distributions of quarks and gluons and discuss its implications for the nucleon spin sum rule. We observe that at small-$x$, the gluon helicity and orbital angular momentum distributions are roughly of the same magnitude but with opposite signs, indicating a significant cancellation between them. A similar cancellation occurs also in the quark sector. We explain analytically the reason for this cancellation.


title: Paraffin coated rubidium cell with an internal atomic vapor source

abstract: We present the results of a study on relaxation and diffusion processes of rubidium atoms in a rubidium cell with an internal vapor source. The cell is an evacuated glass bulb, which is characterized in that the source of atomic vapors in the form of a metal film Rb is evenly distributed throughout the inner surface of the bulb, and the paraffin film is uniformly distributed over the entire area and over the metal surface. By using laser optical pumping, we performed measurements of the relaxation time and the average number of bounces of optical pumped rubidium atoms in the bulb. We have measured the adsorption time of rubidium atoms by paraffin coating and rubidium atoms diffusion coefficient in paraffin used. A simple model of the pumping and atomic diffusion processes in the cell is discussed as well.


title: General Strong Polarization

abstract: Arikan's exciting discovery of polar codes has provided an altogether new way to efficiently achieve Shannon capacity. Given a (constant-sized) invertible matrix $M$, a family of polar codes can be associated with this matrix and its ability to approach capacity follows from the {\em polarization} of an associated $[0,1]$-bounded martingale, namely its convergence in the limit to either $0$ or $1$. Arikan showed polarization of the martingale associated with the matrix $G_2 = \left(\begin{matrix} 1& 0 1& 1\end{matrix}\right)$ to get capacity achieving codes. His analysis was later extended to all matrices $M$ that satisfy an obvious necessary condition for polarization. While Arikan's theorem does not guarantee that the codes achieve capacity at small blocklengths, it turns out that a "strong" analysis of the polarization of the underlying martingale would lead to such constructions. Indeed for the martingale associated with $G_2$ such a strong polarization was shown in two independent works ([Guruswami and Xia, IEEE IT '15] and [Hassani et al., IEEE IT '14]), resolving a major theoretical challenge of the efficient attainment of Shannon capacity. In this work we extend the result above to cover martingales associated with all matrices that satisfy the necessary condition for (weak) polarization. In addition to being vastly more general, our proofs of strong polarization are also simpler and modular. Specifically, our result shows strong polarization over all prime fields and leads to efficient capacity-achieving codes for arbitrary symmetric memoryless channels. We show how to use our analyses to achieve exponentially small error probabilities at lengths inverse polynomial in the gap to capacity. Indeed we show that we can essentially match any error probability with lengths that are only inverse polynomial in the gap to capacity.


title: Fourier Analysis and Evaluation of DG, FD and Compact Difference Methods for Conservation Laws

abstract: Large eddy simulation (LES) has been increasingly used to tackle vortex-dominated turbulent flows. In LES, the quality of the simulation results hinges upon the quality of the numerical discretizations in both space and time. It is in this context we perform a Fourier analysis of several popular methods in LES including the discontinuous Galerkin (DG), finite difference (FD), and compact difference (CD) methods. We begin by reviewing the semi-discrete versions of all methods under-consideration, followed by a fully-discrete analysis with explicit Runge-Kutta (RK) time integration schemes. In this regard, we are able to unravel the true dispersion/dissipation behavior of DG and Runge-Kutta DG (RKDG) schemes for the entire wavenumber range. The physical-mode is verified to be a good approximation for the asymptotic behavior of these DG schemes in the low wavenumber range. After that, we proceed to compare the DG, FD, and CD methods in dispersion and dissipation properties. Numerical tests are conducted using the linear advection equation to verify the analysis. In comparing different methods, it is found that the overall numerical dissipation strongly depends on the time step. Compact difference (CD) and central finite difference (FD) schemes, in some particular settings, can have more numerical dissipation than the DG scheme with an upwind flux. This claim is then verified through a numerical test using the Burgers' equation.


title: A fast C++ implementation of thermal functions

abstract: We provide a small C++ library with Mathematica and Python interfaces for computing thermal functions, defined $$ J_\text{B/F}(y^2) \equiv \Re \int_0^\infty x^2 \log\left[1 \mp e^{-\sqrt{x^2 + y^2}} \right] \,\text{d}x, $$ which appear in finite-temperature quantum field theory and play a role in phase-transitions in the early Universe, including baryogenesis, electroweak symmetry breaking and the Higgs mechanism.


title: Deep Image Super Resolution via Natural Image Priors

abstract: Single image super-resolution (SR) via deep learning has recently gained significant attention in the literature. Convolutional neural networks (CNNs) are typically learned to represent the mapping between low-resolution (LR) and high-resolution (HR) images/patches with the help of training examples. Most existing deep networks for SR produce high quality results when training data is abundant. However, their performance degrades sharply when training is limited. We propose to regularize deep structures with prior knowledge about the images so that they can capture more structural information from the same limited data. In particular, we incorporate in a tractable fashion within the CNN framework, natural image priors which have shown to have much recent success in imaging and vision inverse problems. Experimental results show that the proposed deep network with natural image priors is particularly effective in training starved regimes.


title: B-metric spaces, fixed points and Lipschitz functions

abstract: The paper is concerned with b-metric and generalized b-metric spaces. One proves the existence of the completion of a generalized b-metric space and some fixed point results. The behavior of Lipschitz functions on b-metric spaces of homogeneous type, as well as of Lipschitz functions defined on, or with values in quasi-Banach spaces, is studied.


title: Nevanlinna theory and value distribution in the unicritical polynomials family

abstract: In the space $\mathbb{C}$ of the parameters $\lambda$ of the unicritical polynomials family $f(\lambda,z)=f_\lambda(z)=z^d+\lambda$ of degree $d>1$, we establish a quantitative equidistribution result towards the bifurcation current (indeed measure) $T_f$ of $f$ as $n\to\infty$ on the averaged distributions of all parameters $\lambda$ such that $f_\lambda$ has a superattracting periodic point of period $n$ in $\mathbb{C}$, with a concrete error estimate for $C^2$-test functions on $\mathbb{P}^1$. In the proof, not only complex dynamics but also a standard argument from the Nevanlinna theory play key roles.


title: Primal-dual stochastic gradient method for convex programs with many functional constraints

abstract: Stochastic gradient method (SGM) has been popularly applied to solve optimization problems with objective that is stochastic or an average of many functions. Most existing works on SGMs assume that the underlying problem is unconstrained or has an easy-to-project constraint set. In this paper, we consider problems that have a stochastic objective and also many functional constraints. For such problems, it could be extremely expensive to project a point to the feasible set, or even compute subgradient and/or function value of all constraint functions. To find solutions of these problems, we propose a novel (adaptive) SGM based on the classical augmented Lagrangian function. Within every iteration, it inquires a stochastic subgradient of the objective, and a subgradient and the function value of one randomly sampled constraint function. Hence, the per-iteration complexity is low. We establish its convergence rate for convex problems and also problems with strongly convex objective. It can achieve the optimal $O(1/\sqrt{k})$ convergence rate for convex case and nearly optimal $O\big((\log k)/k\big)$ rate for strongly convex case. Numerical experiments on a sample approximation problem of the robust portfolio selection and quadratically constrained quadratic programming are conducted to demonstrate its efficiency.


title: Composable security analysis of continuous-variable measurement-device-independent quantum key distribution with squeezed states for coherent attacks

abstract: Measurement-device-independent quantum key distribution protocol, whose security analysis does not rely on any assumption on the detection system, can immune the attacking against detectors. We give a first composable security analysis for continuous-variable measurement-device-independent quantum key distribution using squeezed states against general coherent attacks. The security analysis is derived based on the entanglement-based scheme considering finite size effect. A version of entropic uncertainty relation is exploited to give a lower bound on the conditional smooth min-entropy by trusting Alice's and Bob's devices. The simulation results indicate that, in the universal composable security framework, the protocol can tolerate 2.5 dB and 6.5 dB channel loss against coherent attacks with direct and reverse reconciliation, respectively.


title: A Simple proof for the algorithms of relaxed $(u, v)$-cocoercive mappings and $\alpha$-inverse strongly monotone mappings

abstract: In this paper, a simple proof is presented for the convergence of the algorithms for the class of relaxed $(u, v)$-cocoercive mappings and $\alpha$-inverse strongly monotone mappings. Based on $\alpha$-expansive maps, for example, a simple proof of the convergence of the recent iterative algorithms by relaxed $(u, v)$-cocoercive mappings due to Kumam-Jaiboon is provided. Also a simple proof for the convergence of the iterative algorithms by inverse-strongly monotone mappings due to Iiduka-Takahashi in a special case is provided. These results are an improvement as well as a refinement of previously known results.


title: On the Packet Decoding Delay of Linear Network Coded Wireless Broadcast

abstract: We apply linear network coding (LNC) to broadcast a block of data packets from one sender to a set of receivers via lossy wireless channels, assuming each receiver already possesses a subset of these packets and wants the rest. We aim to characterize the average packet decoding delay (APDD), which reflects how soon each individual data packet can be decoded by each receiver on average, and to minimize it while achieving optimal throughput. To this end, we first derive closed-form lower bounds on the expected APDD of all LNC techniques under random packet erasures. We then prove that these bounds are NP-hard to achieve and, thus, that APDD minimization is an NP-hard problem. We then study the performance of some existing LNC techniques, including random linear network coding (RLNC) and instantly decodable network coding (IDNC). We proved that all throughput-optimal LNC techniques can approximate the minimum expected APDD with a ratio between 4/3 and 2. In particular, the ratio of RLNC is exactly 2. We then prove that all IDNC techniques are only heuristics in terms of throughput optimization and {cannot guarantee an APDD approximation ratio for at least a subset of the receivers}. Finally, we propose hyper-graphic linear network coding (HLNC), a novel throughput-optimal and APDD-approximating LNC technique based on a hypergraph model of receivers' packet reception state. We implement it under different availability of receiver feedback, and numerically compare its performance with RLNC and a heuristic general IDNC technique. The results show that the APDD performance of HLNC is better under all tested system settings, even if receiver feedback is only collected intermittently.


title: Dark energy scenario consistent with GW170817 in theories beyond Horndeski gravity

abstract: The Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories up to quartic order are the general scheme of scalar-tensor theories allowing the possibility for realizing the tensor propagation speed $c_t$ equivalent to 1 on the isotropic cosmological background. We propose a dark energy model in which the late-time cosmic acceleration occurs by a simple k-essence Lagrangian analogous to the ghost condensate with cubic and quartic Galileons in the framework of GLPV theories. We show that a wide variety of the variation of the dark energy equation of state $w_{\rm DE}$ including the entry to the region $w_{\rm DE}<-1$ can be realized without violating conditions for the absence of ghosts and Laplacian instabilities. The approach to the tracker equation of state $w_{\rm DE}=-2$ during the matter era, which is disfavored by observational data, can be avoided by the existence of a quadratic k-essence Lagrangian $X^2$. We study the evolution of nonrelativistic matter perturbations for the model $c_t^2=1$ and show that the two quantities $\mu$ and $\Sigma$, which are related to the Newtonian and weak lensing gravitational potentials respectively, are practically equivalent to each other, such that $\mu \simeq \Sigma>1$. For the case in which the deviation of $w_{\rm DE}$ from $-1$ is significant at a later cosmological epoch, the values of $\mu$ and $\Sigma$ tend to be larger at low redshifts. We also find that our dark energy model can be consistent with the bounds on the deviation parameter $\alpha_{\rm H}$ from Horndeski theories arising from the modification of gravitational law inside massive objects.


title: Transverse momentum distributions of electron in simulated QED model

abstract: In the present work, we have studied the transverse momentum distributions (TMDs) for the electron in simulated QED model. We have used the overlap representation of light-front wave functions (LFWFs) where the spin-1/2 relativistic composite system consists of spin-1/2 fermion and spin-1 vector boson. The results have been obtained for T-even TMDs in transverse momentum plane for fixed value of longitudinal momentum fraction $x$.


title: Representation and Characterization of Non-Stationary Processes by Dilation Operators and Induced Shape Space Manifolds

abstract: We have introduce a new vision of stochastic processes through the geometry induced by the dilation. The dilation matrices of a given processes are obtained by a composition of rotations matrices, contain the measure information in a condensed way. Particularly interesting is the fact that the obtention of dilation matrices is regardless of the stationarity of the underlying process. When the process is stationary, it coincides with the Naimark Dilation and only one rotation matrix is computed, when the process is non-stationary, a set of rotation matrices are computed. In particular, the periodicity of the correlation function that may appear in some classes of signal is transmitted to the set of dilation matrices. These rotation matrices, which can be arbitrarily close to each other depending on the sampling or the rescaling of the signal are seen as a distinctive feature of the signal. In order to study this sequence of matrices, and guided by the possibility to rescale the signal, the correct geometrical framework to use with the dilation's theoretic results is the space of curves on manifolds, that is the set of all curve that lies on a base manifold. To give a complete sight about the space of curve, a metric and the derived geodesic equation are provided. The general results are adapted to the more specific case where the base manifold is the Lie group of rotation matrices. The notion of the shape of a curve can be formalized as the set of equivalence classes of curves given by the quotient space of the space of curves and the increasing diffeomorphisms. The metric in the space of curve naturally extent to the space of shapes and enable comparison between shapes.


title: Topologically Controlled Lossy Compression

abstract: This paper presents a new algorithm for the lossy compression of scalar data defined on 2D or 3D regular grids, with topological control. Certain techniques allow users to control the pointwise error induced by the compression. However, in many scenarios it is desirable to control in a similar way the preservation of higher-level notions, such as topological features , in order to provide guarantees on the outcome of post-hoc data analyses. This paper presents the first compression technique for scalar data which supports a strictly controlled loss of topological features. It provides users with specific guarantees both on the preservation of the important features and on the size of the smaller features destroyed during compression. In particular, we present a simple compression strategy based on a topologically adaptive quantization of the range. Our algorithm provides strong guarantees on the bottleneck distance between persistence diagrams of the input and decompressed data, specifically those associated with extrema. A simple extension of our strategy additionally enables a control on the pointwise error. We also show how to combine our approach with state-of-the-art compressors, to further improve the geometrical reconstruction. Extensive experiments, for comparable compression rates, demonstrate the superiority of our algorithm in terms of the preservation of topological features. We show the utility of our approach by illustrating the compatibility between the output of post-hoc topological data analysis pipelines, executed on the input and decompressed data, for simulated or acquired data sets. We also provide a lightweight VTK-based C++ implementation of our approach for reproduction purposes.


title: The Higher-Order Prover Leo-III (Extended Version)

abstract: The automated theorem prover Leo-III for classical higher-order logic with Henkin semantics and choice is presented. Leo-III is based on extensional higher-order paramodulation and accepts every common TPTP dialect (FOF, TFF, THF), including their recent extensions to rank-1 polymorphism (TF1, TH1). In addition, the prover natively supports almost every normal higher-order modal logic. Leo-III cooperates with first-order reasoning tools using translations to many-sorted first-order logic and produces verifiable proof certificates. The prover is evaluated on heterogeneous benchmark sets.


title: From Hashing to CNNs: Training BinaryWeight Networks via Hashing

abstract: Deep convolutional neural networks (CNNs) have shown appealing performance on various computer vision tasks in recent years. This motivates people to deploy CNNs to realworld applications. However, most of state-of-art CNNs require large memory and computational resources, which hinders the deployment on mobile devices. Recent studies show that low-bit weight representation can reduce much storage and memory demand, and also can achieve efficient network inference. To achieve this goal, we propose a novel approach named BWNH to train Binary Weight Networks via Hashing. In this paper, we first reveal the strong connection between inner-product preserving hashing and binary weight networks, and show that training binary weight networks can be intrinsically regarded as a hashing problem. Based on this perspective, we propose an alternating optimization method to learn the hash codes instead of directly learning binary weights. Extensive experiments on CIFAR10, CIFAR100 and ImageNet demonstrate that our proposed BWNH outperforms current state-of-art by a large margin.


title: Is the Spiral Galaxy a Cosmic Hurricane?

abstract: It is discussed that the formation of the spiral galaxies is driven by the cosmic background rotation, not a result of an isolated evolution proposed by the density wave theory. To analyze the motions of the galaxies, a simple double particle galaxy model is considered and the Coriolis force formed by the rotational background is introduced. The numerical analysis shows that not only the trajectory of the particle is the spiral shape, but also the relationship between the velocity and the radius reveals both the existence of spiral arm and the change of the arm number. In addition, the results of the three-dimensional simulation also give the warped structure of the spiral galaxies, and shows that the disc surface of the warped galaxy, like a spinning coins on the table, exists a whole overturning movement. Through the analysis, it can be concluded that the background environment of the spiral galaxies have a large-scale rotation, and both the formation and evolution of hurricane-like spiral galaxies are driven by this background rotation.


title: A new presentation of the plane Cremona group

abstract: We give a presentation of the plane Cremona group over an algebraically closed field with respect to the generators given by the Theorem of Noether and Castelnuovo. This presentation is particularly simple and can be used for explicit calculations.


title: Autonomous Power Allocation based on Distributed Deep Learning for Device-to-Device Communication Underlaying Cellular Network

abstract: For Device-to-device (D2D) communication of Internet-of-Things (IoT) enabled 5G system, there is a limit to allocating resources considering a complicated interference between different links in a centralized manner. If D2D link is controlled by an enhanced node base station (eNB), and thus, remains a burden on the eNB and it causes delayed latency. This paper proposes a fully autonomous power allocation method for IoT-D2D communication underlaying cellular networks using deep learning. In the proposed scheme, an IoT-D2D transmitter decides the transmit power independently from an eNB and other IoT-D2D devices. In addition, the power set can be nearly optimized by deep learning with distributed manner to achieve higher cell throughput. We present a distributed deep learning architecture in which the devices are trained as a group but operate independently. The deep learning can attain near optimal cell throughput while suppressing interference to eNB.


title: The dynamics of disappearing pulses in a singularly perturbed reaction-diffusion system with parameters that vary in time and space

abstract: We consider the evolution of multi-pulse patterns in an extended Klausmeier equation with parameters that change in time and/or space. We formally show that the full PDE dynamics of a $N$-pulse configuration can be reduced to a $N$-dimensional dynamical system describing the dynamics on a $N$-dimensional manifold $\mathcal{M}_N$. Next, we determine the local stability of $\mathcal{M}_N$ via the quasi-steady spectrum associated to evolving $N$-pulse patterns, which provides provides explicit information on the boundary $\partial\mathcal{M}_N$. Following the dynamics on $\mathcal{M}_N$, a $N$-pulse pattern may move through $\partial\mathcal{M}_N$ and `fall off' $\mathcal{M}_N$. A direct nonlinear extrapolation of our linear analysis predicts the subsequent fast PDE dynamics as the pattern `jumps' to another invariant manifold $\mathcal{M}_M$, and specifically predicts the number $N-M$ of pulses that disappear. Combining the asymptotic analysis with numerical simulations of the dynamics on the various invariant manifolds yields a hybrid asymptotic-numerical method describing the full process that starts with a $N$-pulse pattern and typically ends in the trivial homogeneous state without pulses. We extensively test this method against PDE simulations and deduce general conjectures on the nature of pulse interactions with disappearing pulses. We especially consider the differences between the evolution of irregular and regular patterns. In the former case, the disappearing process is gradual: irregular patterns loose their pulses one by one. In contrast, regular, spatially periodic, patterns undergo catastrophic transitions in which either half or all pulses disappear. However, making a precise distinction between these two drastically different processes is quite subtle, since irregular $N$-pulse patterns that do not cross $\partial\mathcal{M}_N$ typically evolve towards regularity.


title: The topology of the set of non-escaping endpoints

abstract: There are several classes of transcendental entire functions for which the Julia set consists of an uncountable union of disjoint curves each of which joins a finite endpoint to infinity. Many authors have studied the topological properties of this set of finite endpoints. It was recently shown that, for certain functions in the exponential family, there is a strong dichotomy between the topological properties of the set of endpoints which escape and those of the set of endpoints which do not escape. In this paper, we show that this result holds for large families of functions in the Eremenko-Lyubich class. We also show that this dichotomy holds for a family of functions, outside that class, which includes the much-studied Fatou function defined by $f(z) := z + 1+ e^{-z}.$ Finally, we show how our results can be used to demonstrate that various sets are spiders' webs, generalising results such as those in a recent paper of the first author.


title: Transverse Extension of Partons in the Proton probed by Deeply Virtual Compton Scattering

abstract: We report on the first measurement of exclusive single-photon muoproduction on the proton by COMPASS using 160 GeV/$c$ polarized $\mu^+$ and $\mu^-$ beams of the CERN SPS impinging on a liquid hydrogen target. We determine the dependence of the average of the measured $\mu^+$ and $\mu^-$ cross sections for deeply virtual Compton scattering on the squared four-momentum transfer $t$ from the initial to the final final proton. The slope $B$ of the $t$-dependence is fitted with a single exponential function, which yields $B=(4.3 \ \pm \ 0.6_{\text{stat}}\_{- \ 0.3}^{+ \ 0.1}\big\rvert_{\text{sys}}) (\text{GeV}/c)^{-2}$. This result can be converted into an average transverse extension of partons in the proton, $\sqrt{\langle r_{\perp}^2 \rangle} = (0.58 \ \pm \ 0.04_{\text{stat}}\_{- \ 0.02}^{+ \ 0.01}\big\rvert_{\text{sys}})\text{fm}$. For this measurement, the average virtuality of the photon mediating the interaction is $\langle Q^2 \rangle = 1.8\,(\text{GeV/}c)^2$ and the average value of the Bjorken variable is $\langle x_{\text{Bj}} \rangle = 0.056$.


title: Relativistic coupled-cluster theory analysis of energies, hyperfine structure constants, and dipole polarizabilities of Cd$^{+}$

abstract: Roles of electron correlation effects in the determination of attachment energies, magnetic dipole hyperfine structure constants and electric dipole (E1) matrix elements of the low-lying states in the singly charged cadmium ion (Cd$^+$) have been analyzed. We employ the singles and doubles approximated relativistic coupled-cluster (RCC) method to calculate these properties. Intermediate results from the Dirac-Hartree-Fock approximation, second-order many-body perturbation theory and considering only the linear terms of the RCC method are given to demonstrate propagation of electron correlation effects in this ion. Contributions from important RCC terms are also given to highlight importance of various correlation effects in the evaluation of these properties. At the end, we also determine E1 polarizabilities ($\alpha^{E1}$) of the ground and $5p \ ^2P_{1/2;3/2}$ states of Cd$^+$ in the {\it ab initio} approach. We estimate them again by replacing some of the E1 matrix elements and energies from the measurements to reduce their uncertainties so that they can be used in the high precision experiments of this ion.


title: Average number of zeros and mixed symplectic volume of Finsler sets

abstract: Let $X$ be an $n$-dimensional manifold and $V_1, \ldots, V_n \subset C^\infty(X, \mathbb R)$ finite-dimensional vector spaces with Euclidean metric. We assign to each $V_i$ a Finsler ellipsoid, i.e., a family of ellipsoids in the fibers of the cotangent bundle of $X$. We prove that the average number of isolated common zeros of $f_1 \in V_1, \ldots, f_n \in V_n$ is equal to the mixed symplectic volume of these Finsler ellipsoids. If $X$ is a homogeneous space of a compact Lie group and all vector spaces $V_i$ and their Euclidean metrics are invariant, then the average numbers of zeros satisfy the inequalities, similar to Hodge inequalities for intersection numbers of divisors on a projective variety. This is applied to the eigenspaces of Laplace operator of an invariant Riemannian metric. The proofs are based on a construction of the ring of normal densities on $X$, an analogue of the ring of differential forms. In particular, this construction is used for a generalization of Crofton formula to the product of spheres.


title: Ground-based lightcurve observation campaign of (25143) Itokawa between 2001 and 2004

abstract: The asteroid (25143) Itokawa is a target object of the Japanese sample return mission, HAYABUSA. We have observed Itokawa in optical wave- length (R-band) with the 1.05-m Schmidt telescope at the Kiso Observatory, the 2.24-m telescope of University of Hawaii, and the 1.05-m telescope at the Misato Observatory since 2001. From the analysis of the data, we present the relationship between brightness and the solar phase angle, 6.9 to 87.8 deg. We obtained the absolute magnitude H_R(0) = 19.09+-0.37, and the slope parameter G_R = 0.25 +- 0.29. The rotational period of Itokawa is 12.1324 +- 0.0001 hours.


title: Crossroads of the mesoscale circulation

abstract: Quantifying the mechanisms of tracer dispersion in the ocean remains a central question in oceanography, for problems ranging from nutrient delivery to phytoplankton, to the early detection of contaminants. Most analyses have been based on Lagrangian concepts of transport, focusing on the identification of features minimizing fluid exchange among regions, or more recently on network tools which focus on connectivity and transport pathways. Neither of these approaches allows ranking the geographical sites of major water passage and selecting them so that they monitor waters coming from separate parts of the ocean. These are instead key criteria when deploying an observing network. Here we address this issue by estimating at any point the extent of the ocean surface which transits through it in a given time window. With such information we are able to rank the sites with major fluxes that intercept waters originating from different regions. We show that this allows us to optimize an observing network, where a set of sampling sites can be chosen for monitoring the largest flux of water dispersing out of a given region. When the analysis is performed backward in time, this method allows us to identify the major sources which feed a target region. The method is first applied to a minimalistic model of a mesoscale eddy field, and then to realistic satellite-derived ocean currents in the Kerguelen area. In this region we identify the optimal location of fixed stations capable of intercepting the trajectories of 43 surface drifters, along with statistics on the temporal persistence of the stations determined in this way. We then identify possible hotspots of micro-nutrient enrichment for the recurrent spring phytoplanktonic bloom occuring here. Promising applications to other fields, such as larval connectivity, marine spatial planning or contaminant detection, are then discussed.


title: Minimizing Latency in Online Ride and Delivery Services

abstract: Motivated by the popularity of online ride and delivery services, we study natural variants of classical multi-vehicle minimum latency problems where the objective is to route a set of vehicles located at depots to serve request located on a metric space so as to minimize the total latency. In this paper, we consider point-to-point requests that come with source-destination pairs and release-time constraints that restrict when each request can be served. The point-to-point requests and release-time constraints model taxi rides and deliveries. For all the variants considered, we show constant-factor approximation algorithms based on a linear programming framework. To the best of our knowledge, these are the first set of results for the aforementioned variants of the minimum latency problems. Furthermore, we provide an empirical study of heuristics based on our theoretical algorithms on a real data set of taxi rides.


title: Learning Inductive Biases with Simple Neural Networks

abstract: People use rich prior knowledge about the world in order to efficiently learn new concepts. These priors - also known as "inductive biases" - pertain to the space of internal models considered by a learner, and they help the learner make inferences that go beyond the observed data. A recent study found that deep neural networks optimized for object recognition develop the shape bias (Ritter et al., 2017), an inductive bias possessed by children that plays an important role in early word learning. However, these networks use unrealistically large quantities of training data, and the conditions required for these biases to develop are not well understood. Moreover, it is unclear how the learning dynamics of these networks relate to developmental processes in childhood. We investigate the development and influence of the shape bias in neural networks using controlled datasets of abstract patterns and synthetic images, allowing us to systematically vary the quantity and form of the experience provided to the learning algorithms. We find that simple neural networks develop a shape bias after seeing as few as 3 examples of 4 object categories. The development of these biases predicts the onset of vocabulary acceleration in our networks, consistent with the developmental process in children.


title: Rank Revealing Gaussian Elimination by the Maximum Volume Concept

abstract: A Gaussian elimination algorithm is presented that reveals the numerical rank of a matrix by yielding small entries in the Schur complement. The algorithm uses the maximum volume concept to find a square nonsingular submatrix of maximum dimension. The bounds on the revealed singular values are similar to the best known bounds for rank revealing LU factorization, but in contrast to existing methods the algorithm does not make use of the normal matrix. An implementation for dense matrices is described whose computational cost is roughly twice the cost of an LU factorization with complete pivoting. Because of its flexibility in choosing pivot elements, the algorithm is amenable to implementation with blocked memory access and for sparse matrices.


title: Orbital Motion of Young Binaries in Ophiuchus and Upper Centaurus-Lupus

abstract: We present measurements of the orbital positions and flux ratios of 17 binary and triple systems in the Ophiuchus star forming region and the Upper Centaurus-Lupus cluster based on adaptive optics imaging at the Keck Observatory. We report the detection of visual companions in MML 50 and MML 53 for the first time, as well as the possible detection of a third component in WSB 21. For six systems in our sample, our measurements provide a second orbital position following their initial discoveries over a decade ago. For eight systems with sufficient orbital coverage, we analyze the range of orbital solutions that fit the data. Ultimately, these observations will help provide the groundwork toward measuring precise masses for these pre-main sequence stars and understanding the distribution of orbital parameters in young multiple systems.


title: Serve the shortest queue and Walsh Brownian motion

abstract: We study a single-server Markovian queueing model with $N$ customer classes in which priority is given to the shortest queue. Under a critical load condition, we establish the diffusion limit of the workload and queue length processes in the form of a Walsh Brownian motion (WBM) living in the union of the $N$ nonnegative coordinate axes in $\mathbb{R}^N$ and a linear transformation thereof. This reveals the following asymptotic behavior. Each time that queues begin to build starting from an empty system, one of them becomes dominant in the sense that it contains nearly all the workload in the system, and it remains so until the system becomes (nearly) empty again. The radial part of the WBM, given as a reflected Brownian motion (RBM) on the half-line, captures the total workload asymptotics, whereas its angular distribution expresses how likely it is for each class to become dominant on excursions. As a heavy traffic result it is nonstandard in three ways: (i) In the terminology of Harrison (1995) it is unconventional, in that the limit is not an RBM. (ii) It does not constitute an invariance principle, in that the limit law (specifically, the angular distribution) is not determined solely by the first two moments of the data, and is sensitive even to tie breaking rules. (iii) The proof method does not fully characterize the limit law (specifically, it gives no information on the angular distribution).


title: Ultraviolet laser spectroscopy of aluminum atoms in hollow-cathode lamp

abstract: We report precision measurement of aluminum atoms ${^{2}P_{1/2}}-{^{2}S_{1/2}}$ transition at 394 nm and ${^{2}P_{3/2}}-{^{2}S_{1/2}}$ transition at 396 nm in a hollow-cathode lamp (HCL). Both absorption spectroscopy and saturated absorption spectroscopy (SAS) are performed. From the absorption spectroscopy the Doppler linewidth is estimated to be 2.6 GHz. The SAS spectroscopy is analyzed based on the velocity-changing-effect model. With a frequency comb calibrated wavemeter, the frequencies of ${^{2}P_{1/2}},{F=3}-{^{2}S_{1/2}},{F=2}$ transition and ${^{2}P_{3/2}},{F=4}-{^{2}S_{1/2}},{F=3}$ transition are measured to be 759.905401(10) THz and 756.547403(10) THz, respectively. The hyperfine structure constants of aluminum atoms are determined and compared with previously reported measurement results and theoretical calculation. Reasonable agreement is found for the magnetic dipole constant (A constant), while the electric quadrupole constant (B constant) has a large deviation.


title: Analytic plasma wakefield limits for active plasma lenses

abstract: Active plasma lensing is a promising technology for compact focusing of particle beams that has seen a recent surge of interest. While these lenses can provide strong focusing gradients of order kT/m and focusing in both transverse planes, there are limitations from nonlinear aberrations, causing emittance growth in the beams being focused. One cause of such aberrations is beam-driven plasma wakefields, present if the beam density is sufficiently high. We develop simple, but powerful analytic formulas for the effective focusing gradient from these wakefields, and use this to set limits on which parts of the beam and plasma parameter space permits distortion-free use of active plasma lenses. It is concluded that the application of active plasma lenses to conventional and plasma-based linear colliders may prove very challenging, except perhaps in the final focus system, unless the typical discharge currents used are dramatically increased, and that in general these lenses are better suited for accelerator applications with lower beam intensities.


title: Monopoly pricing with buyer search

abstract: In many shopping scenarios, e.g., in online shopping, customers have a large menu of options to choose from. However, most of the buyers do not browse all the options and make decision after considering only a small part of the menu. To study such buyer's behavior we consider the standard Bayesian monopoly problem for a unit-demand buyer, where the monopolist displays the menu dynamically page after a page to the buyer. The seller aims to maximize the expected revenue over the distribution of buyer's values which we assume are i.i.d. The buyer incurs a fixed cost for browsing through one menu page and would stop if that cost exceeds the increase in her utility. We observe that the optimal posted price mechanism in our dynamic setting may have quite different structure than in the classic static scenario. We find a (relatively) simple and approximately optimal mechanism, that uses part of the items as a "bait" to keep the buyer interested for multiple rounds with low prices, while at the same time showing many other expensive items.


title: Degree bound of P\'olya Positivstellenstaz

abstract: P\'olya's Positivstellensatz on the $1$-simplex says that if $P(x)$ is a real polynomial such that $P(x)>0$ whenever $x \ge 0$, then all the coefficients of $(1+x)^mP(x)$ are positive whenever $m$ is large. Powers-Reznick gave a complexity estimate for P\'olya's Positivstellensatz. Namely, they proved that, for such $P(x)$ of degree $d$, all the coefficients of $(1+x)^mP(x)$ are positive whenever $m > \frac{1}{2} (d^2 -d) \frac{L(P)}{\lambda(P)} - d$. where $\frac{L(P)}{\lambda(P)}$ is an invariant of $P(x)$. For $d=3$ and $d=4$ specifically, we improve Powers-Reznick's bound by showing $m > \frac{3}{2} \frac{L(P)}{\lambda(P)} - 1$ for $d=3$ and $ m > \frac{4232}{2505} \frac{L(P)}{\lambda(P)} - 1$ for $d=4$.