title: Energy Transparency for Deeply Embedded Programs

abstract: Energy transparency is a concept that makes a program's energy consumption visible, from hardware up to software, through the different system layers. Such transparency can enable energy optimizations at each layer and between layers, and help both programmers and operating systems make energy-aware decisions. In this paper, we focus on deeply embedded devices, typically used for Internet of Things (IoT) applications, and demonstrate how to enable energy transparency through existing Static Resource Analysis (SRA) techniques and a new target-agnostic profiling technique, without hardware energy measurements. Our novel mapping technique enables software energy consumption estimations at a higher level than the Instruction Set Architecture (ISA), namely the LLVM Intermediate Representation (IR) level, and therefore introduces energy transparency directly to the LLVM optimizer. We apply our energy estimation techniques to a comprehensive set of benchmarks, including single- and also multi-threaded embedded programs from two commonly used concurrency patterns, task farms and pipelines. Using SRA, our LLVM IR results demonstrate a high accuracy with a deviation in the range of 1% from the ISA SRA. Our profiling technique captures the actual energy consumption at the LLVM IR level with an average error of 3%.


title: Corrective Control to Handle Forecast Uncertainty: A Chance Constrained Optimal Power Flow

abstract: Higher shares of electricity generation from renewable energy sources and market liberalization is increasing uncertainty in power systems operation. At the same time, operation is becoming more flexible with improved control systems and new technology such as phase shifting transformers (PSTs) and high voltage direct current connections (HVDC). Previous studies have shown that the use of corrective control in response to outages contributes to a reduction in operating cost, while maintaining N-1 security. In this work, we propose a method to extend the use of corrective control of PSTs and HVDCs to react to uncertainty. We characterize the uncertainty as continuous random variables, and define the corrective control actions through affine control policies. This allows us to efficiently model control reactions to a large number of uncertainty sources. The control policies are then included in a chance constrained optimal power flow formulation, which guarantees that the system constraints are enforced with a desired probability. By applying an analytical reformulation of the chance constraints, we obtain a second-order cone problem for which we develop an efficient solution algorithm. In a case study for the IEEE 118 bus system, we show that corrective control for uncertainty leads to a decrease in operational cost, while maintaining system security. Further, we demonstrate the scalability of the method by solving the problem for the IEEE 300 bus and the Polish system test cases.


title: Einstein's $R^{\hat{0} \hat{0}}$ equation for non-relativistic sources derived from Einstein's inertial motion and the Newtonian law for relative acceleration

abstract: With Einstein's inertial motion (free-falling and non-rotating relative to gyroscopes), geodesics for non-relativistic particles can intersect repeatedly, allowing one to compute the space-time curvature $R^{\hat{0} \hat{0}}$ exactly. Einstein's $R^{\hat{0} \hat{0}}$ for strong gravitational fields and for relativistic source-matter is identical with the Newtonian expression for the relative radial acceleration of neighboring free-falling test-particles, spherically averaged.--- Einstein's field equations follow from Newtonian experiments, local Lorentz-covariance, and energy-momentum conservation combined with the Bianchi identity.


title: Conformal covariance and the split property

abstract: We show that for a conformal local net of observables on the circle, the split property is automatic. Both full conformal covariance (i.e. diffeomorphism covariance) and the circle-setting play essential roles in this fact, while by previously constructed examples it was already known that even on the circle, M\"obius covariance does not imply the split property. On the other hand, here we also provide an example of a local conformal net living on the two-dimensional Minkowski space, which - although being diffeomorphism covariant - does not have the split property.


title: Pilot Contamination Attack Detection by Key-Confirmation in Secure MIMO Systems

abstract: Many security techniques working at the physical layer need a correct channel state information (CSI) at the transmitter, especially when devices are equipped with multiple antennas. Therefore such techniques are vulnerable to pilot contamination attacks (PCAs) by which an attacker aims at inducing false CSI. In this paper we provide a solution to some PCA methods, by letting two legitimate parties to compare their channel estimates. The comparison is made in order to minimize the information leakage on the channel to a possible attacker. By reasonable assumptions on both the channel knowledge by the attacker and the correlation properties of the attacker and legitimate channels we show the validity of our solution. An accurate analysis of possible attacks and countermeasures is provided, together with a numerical evaluation of the attainable secrecy outage probability when our solution is used in conjunction with beamforming for secret communications.


title: Sequential Linear Quadratic Optimal Control for Nonlinear Switched Systems

abstract: In this contribution, we introduce an efficient method for solving the optimal control problem for an unconstrained nonlinear switched system with an arbitrary cost function. We assume that the sequence of the switching modes are given but the switching time in between consecutive modes remains to be optimized. The proposed method uses a two-stage approach as introduced by Xu and Antsaklis (2004) where the original optimal control problem is transcribed into an equivalent problem parametrized by the switching times and the optimal control policy is obtained based on the solution of a two-point boundary value differential equation. The main contribution of this paper is to use a Sequential Linear Quadratic approach to synthesize the optimal controller instead of solving a boundary value problem. The proposed method is numerically more efficient and scales very well to the high dimensional problems. In order to evaluate its performance, we use two numerical examples as benchmarks to compare against the baseline algorithm. In the third numerical example, we apply the proposed algorithm to the Center of Mass control problem in a quadruped robot locomotion task.


title: Localization and Fractionalization in a Chain of Rotating Atomic Gases

abstract: We describe the phase diagram and thermodynamic properties of a chain of axially-tunnel-coupled fractional quantum Hall systems realized by rotating a series of optical dipole traps about their center. We demonstrate not only a experimentally feasible pathway to a state describable as a Mott-insulator of composite bosons, but also describe the nature of the coherent states at higher tunnel coupling strength, and identify a series of new superfluid phases with rich behavior. The phase diagram directly reveals not only characteristic features of the few-body systems, including the effective mass of composite particle- and hole-like excitations and their interactions, but emergent properties of the chain also reveal a fundamental mapping between the adiabatic dynamics of two dimensional systems governed by particle braiding and the hydrostatic response of the gas in the conducting phases.


title: Discrete Variational Autoencoders

abstract: Probabilistic models with discrete latent variables naturally capture datasets composed of discrete classes. However, they are difficult to train efficiently, since backpropagation through discrete variables is generally not possible. We present a novel method to train a class of probabilistic models with discrete latent variables using the variational autoencoder framework, including backpropagation through the discrete latent variables. The associated class of probabilistic models comprises an undirected discrete component and a directed hierarchical continuous component. The discrete component captures the distribution over the disconnected smooth manifolds induced by the continuous component. As a result, this class of models efficiently learns both the class of objects in an image, and their specific realization in pixels, from unsupervised data, and outperforms state-of-the-art methods on the permutation-invariant MNIST, Omniglot, and Caltech-101 Silhouettes datasets.


title: Chemical abundances for A-and F-type supergiant stars

abstract: We present the stellar parameters and elemental abundances of a set of A--F-type supergiant stars HD\,45674, HD\,180028, HD\,194951 and HD\,224893 using high resolution ($R$\,$\sim$\,42,000) spectra taken from ELODIE library. We present the first results of the abundance analysis for HD\,45674 and HD\,224893. We reaffirm the abundances for HD\,180028 and HD\,194951 studied previously by Luck (2014) respectively. Alpha-elements indicates that objects belong to the thin disc population. From their abundances and its location on the Hertzsprung-Russell diagram seems point out that HD\,45675, HD\,194951 and HD\,224893 are in the post-first dredge-up (post-1DUP) phase and they are moving in the red-blue loop region. HD~180028, on the contary, shows typical abundances of the population I but its evolutionary status could not be satisfactorily defined.


title: Deconstructing Odorant Identity via Primacy in Dual Networks

abstract: In the olfactory system, odor percepts retain their identity despite substantial variations in concentration, timing, and background. We propose a novel strategy for encoding intensity-invariant stimuli identity that is based on representing relative rather than absolute values of the stimulus features. Because, in this scheme, stimulus identity depends on relative amplitudes of stimulus features, identity becomes invariant with respect to variations in intensity and monotonous non-linearities of neuronal responses. In the olfactory system, stimulus identity can be represented by the identities of the p strongest responding odorant receptor types out of a species dependent complement. We show that this information is sufficient to recover sparse stimuli (odorants) via elastic net loss minimization. Such a minimization has to be performed under constraints imposed by the relationships between stimulus features. We map this problem onto the dual problem of minimizing a functional of Lagrange multipliers. The dual problem, in turn, can be solved by a neural network whose Lyapunov function represents the dual Lagrangian. We thus propose that networks in the piriform cortex compute odorant identity and implement dual computations with the sparse activities of individual neurons representing the Lagrange multipliers.


title: First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

abstract: We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola component with curvature beta = 0.3(+0.2,-0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p = 2.2. Assuming that the power-law electron distribution extends down to the minimum electron Lorentz factor gamma_min = 1 and that there is one proton per electron, an unrealistically high total jet power L_p of roughly 10^47 erg/s is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio of at least 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma_br1 of roughly 100. In either case, the jet composition is expected to be strongly matter-dominated.


title: A comprehensive search for stellar bowshock nebulae in the Milky Way: a catalog of 709 mid-infrared selected candidates

abstract: We identify 709 arc-shaped mid-infrared nebula in 24 micron Spitzer Space Telescope or 22 micron Wide Field Infrared Explorer surveys of the Galactic Plane as probable dusty interstellar bowshocks powered by early-type stars. About 20% are visible at 8 microns or shorter mid-infrared wavelengths as well. The vast majority (660) have no previous identification in the literature. These extended infrared sources are strongly concentrated near Galactic mid-Plane with an angular scale height of ~0.6 degrees. All host a symmetrically placed star implicated as the source of a stellar wind sweeping up interstellar material. These are candidate "runaway" stars potentially having high velocities in the reference frame of the local medium. Among the 286 objects with measured proper motions, we find an unambiguous excess having velocity vectors aligned with the infrared morphology --- kinematic evidence that many of these are "runaway" stars with large peculiar motions responsible for the bowshock signature. We discuss a population of "in-situ" bowshocks (103 objects) that face giant HII regions where the relative motions between the star and ISM may be caused by bulk outflows from an overpressured bubble. We also identify 58 objects that face 8 micron bright-rimmed clouds and apparently constitute a sub-class of in-situ bowshocks where the stellar wind interacts with a photo-evaporative flow from an eroding molecular cloud interface (i.e., "PEF bowshocks"). Orientations of the arcuate nebulae exhibit a correlation over small angular scales, indicating that external influences such as HII regions are responsible for producing some bowshock nebulae. However, the vast majority of this sample appear to be isolated (499 objects) from obvious external influences.


title: On the Determination of Star Bodies from Their Half-Sections

abstract: We obtain explicit inversion formulas for the Radon-like transform that assigns to a function on the unit sphere the integrals of that function over hemispheres lying in lower dimensional central cross-sections. The results are applied to determination of star bodies from the volumes of their central half-sections.


title: A couple of real hyperbolic disc bundles over surfaces

abstract: Applying the techniques developed in [AGG], we construct new real hyperbolic manifolds whose underlying topology is that of a disc bundle over a closed orientable surface. By the Gromov-Lawson-Thurston conjecture [GLT], such bundles $M\to S$ should satisfy the inequality $|eM/\chi S|\leqslant1$, where $eM$ stands for the Euler number of the bundle and $\chi S$, for the Euler characteristic of the surface. In this paper, we construct new examples that provide a maximal value of $|eM/\chi S|=\frac35$ among all known examples. The former maximum, belonging to Feng Luo [Luo], was $|eM/\chi S|=\frac12$.


title: Nuclear Physics and Astrophysics of Neutrino Oscillations

abstract: For a long time very little experimental information was available about neutrino properties, even though a minute neutrino mass has intriguing cosmological and astrophysical implications. This situation has changed in recent decades: intense experimental activity to measure many neutrino properties took place. Some of these developments and their implications for astrophysics and cosmology are briefly reviewed with a particular emphasis on neutrino magnetic moments and collective neutrino oscillations


title: Breaking the Bandwidth Barrier: Geometrical Adaptive Entropy Estimation

abstract: Estimators of information theoretic measures such as entropy and mutual information are a basic workhorse for many downstream applications in modern data science. State of the art approaches have been either geometric (nearest neighbor (NN) based) or kernel based (with a globally chosen bandwidth). In this paper, we combine both these approaches to design new estimators of entropy and mutual information that outperform state of the art methods. Our estimator uses local bandwidth choices of $k$-NN distances with a finite $k$, independent of the sample size. Such a local and data dependent choice improves performance in practice, but the bandwidth is vanishing at a fast rate, leading to a non-vanishing bias. We show that the asymptotic bias of the proposed estimator is universal; it is independent of the underlying distribution. Hence, it can be pre-computed and subtracted from the estimate. As a byproduct, we obtain a unified way of obtaining both kernel and NN estimators. The corresponding theoretical contribution relating the asymptotic geometry of nearest neighbors to order statistics is of independent mathematical interest.


title: A lower bound on the spectrum of unimodular networks

abstract: Unimodular networks are a generalization of finite graphs in a stochastic sense. We prove a lower bound to the spectral radius of the adjacency operator and of the Markov operator of an unimodular network in terms of its average degree. This allows to prove an Alon-Boppana type bound for the largest eigenvalues in absolute value of large, connected, bounded degree graphs, which generalizes the Alon-Boppana theorem for regular graphs. A key step is establishing a lower bound to the spectral radius of a unimodular tree in terms of its average degree. Similarly, we provide a lower bound on the volume growth rate of an unimodular tree in terms of its average degree.


title: Enumerating cycles in the graph of overlapping permutations

abstract: The graph of overlapping permutations is a directed graph that is an analogue to the De Bruijn graph. It consists of vertices that are permutations of length $n$ and edges that are permutations of length $n+1$ in which an edge $a_1\cdots a_{n+1}$ would connect the standardization of $a_1\cdots a_n$ to the standardization of $a_2\cdots a_{n+1}$. We examine properties of this graph to determine where directed cycles can exist, to count the number of directed $2$-cycles within the graph, and to enumerate the vertices that are contained within closed walks and directed cycles of more general lengths.


title: Quasi-stability and Exponential Attractors for A Non-Gradient System---Applications to Piston-Theoretic Plates with Internal Damping

abstract: We consider a nonlinear (Berger or Von Karman) clamped plate model with a {\em piston-theoretic} right hand side---which include non-dissipative, non-conservative lower order terms. The model arises in aeroelasticity when a panel is immersed in a high velocity linear potential flow; in this case the effect of the flow can be captured by a dynamic pressure term written in terms of the material derivative of the plate's displacement. The effect of fully-supported internal damping is studied for both Berger and von Karman dynamics. The non-dissipative nature of the dynamics preclude the use of strong tools such as backward-in-time smallness of velocities and finiteness of the dissipation integral. Modern quasi-stability techniques are utilized to show the existence of compact global attractors and generalized fractal exponential attractors. Specific results depending on the size of the damping parameter and the nonlinearity in force. For the Berger plate, in the presence of large damping, the existence of a proper global attractor (whose fractal dimension is finite in the state space) is shown via a decomposition of the nonlinear dynamics. This leads to the construction of a compact set upon which quasi-stability theory can be implemented. Numerical investigations for appropriate 1-D models are presented which explore and support the abstract results presented herein.


title: Explicit symplectic approximation of nonseparable Hamiltonians: algorithm and long time performance

abstract: Explicit symplectic integrators have been important tools for accurate and efficient approximations of mechanical systems with separable Hamiltonians. For the first time, the article proposes for arbitrary Hamiltonians similar integrators, which are explicit, of any even order, symplectic in an extended phase space, and with pleasant long time properties. They are based on a mechanical restraint that binds two copies of phase space together. Using backward error analysis, KAM theory, and additional multiscale analysis, an error bound of $\mathcal{O}(T\delta^l \omega)$ is established for integrable systems, where $T$, $\delta$, $l$ and $\omega$ are respectively the (long) simulation time, step size, integrator order, and some binding constant. For non-integrable systems with positive Lyapunov exponents, such an error bound is generally impossible, but satisfactory statistical behaviors were observed in a numerical experiment with a nonlinear Schr\"{o}dinger equation.

categories: math.NA math.DS physics.class-ph physics.comp-ph


title: Unsupervised cryo-EM data clustering through adaptively constrained K-means algorithm

abstract: In single-particle cryo-electron microscopy (cryo-EM), K-means clustering algorithm is widely used in unsupervised 2D classification of projection images of biological macromolecules. 3D ab initio reconstruction requires accurate unsupervised classification in order to separate molecular projections of distinct orientations. Due to background noise in single-particle images and uncertainty of molecular orientations, traditional K-means clustering algorithm may classify images into wrong classes and produce classes with a large variation in membership. Overcoming these limitations requires further development on clustering algorithms for cryo-EM data analysis. We propose a novel unsupervised data clustering method building upon the traditional K-means algorithm. By introducing an adaptive constraint term in the objective function, our algorithm not only avoids a large variation in class sizes but also produces more accurate data clustering. Applications of this approach to both simulated and experimental cryo-EM data demonstrate that our algorithm is a significantly improved alterative to the traditional K-means algorithm in single-particle cryo-EM analysis.


title: Automated Segmentation of Retinal Layers from Optical Coherent Tomography Images Using Geodesic Distance

abstract: Optical coherence tomography (OCT) is a non-invasive imaging technique that can produce images of the eye at the microscopic level. OCT image segmentation to localise retinal layer boundaries is a fundamental procedure for diagnosing and monitoring the progression of retinal and optical nerve disorders. In this paper, we introduce a novel and accurate geodesic distance method (GDM) for OCT segmentation of both healthy and pathological images in either two- or three-dimensional spaces. The method uses a weighted geodesic distance by an exponential function, taking into account both horizontal and vertical intensity variations. The weighted geodesic distance is efficiently calculated from an Eikonal equation via the fast sweeping method. The segmentation is then realised by solving an ordinary differential equation with the geodesic distance. The results of the GDM are compared with manually segmented retinal layer boundaries/surfaces. Extensive experiments demonstrate that the proposed GDM is robust to complex retinal structures with large curvatures and irregularities and it outperforms the parametric active contour algorithm as well as the graph theoretic based approaches for delineating the retinal layers in both healthy and pathological images.


title: A Trip to the Cataclysmic Binary Zoo: Detailed Follow-Up of 35 Recently-Discovered Systems

abstract: We report follow-up studies of 35 recently-discovered cataclysmic variables (CVs), 32 of which were found in large, automated synoptic sky surveys. The objects were selected for observational tractability. For 34 of the objects we present mean spectra and spectroscopic orbital periods, and for one more we give an eclipse-based period. Thirty-two of the period determinations are new, and three of these refine published estimates based on superhump periods. The remaining three of our determinations confirm previously published periods. Twenty of the stars are confirmed or suspected dwarf novae with periods shorter than 3 hours, but we also find three apparent polars (AM Her stars), and six systems with P > 5 h, five of which have secondary stars visible in their spectra, from which we estimate distances when possible. The orbital period distribution of this sample is very similar to that of previously discovered CVs.


title: A Galerkin BEM for high-frequency scattering problems based on frequency dependent changes of variables

abstract: In this paper we develop a class of efficient Galerkin boundary element methods for the solution of two-dimensional exterior single-scattering problems. Our approach is based upon construction of Galerkin approximation spaces confined to the asymptotic behaviour of the solution through a certain direct sum of appropriate function spaces weighted by the oscillations in the incident field of radiation. Specifically, the function spaces in the illuminated/shadow regions and the shadow boundaries are simply algebraic polynomials whereas those in the transition regions are generated utilizing novel, yet simple, \emph{frequency dependent changes of variables perfectly matched with the boundary layers of the amplitude} in these regions. While, on the one hand, we rigorously verify for smooth convex obstacles that these methods require only an $\mathcal{O}\left( k^{\epsilon} \right)$ increase in the number of degrees of freedom to maintain any given accuracy independent of frequency, and on the other hand, remaining in the realm of smooth obstacles they are applicable in more general single-scattering configurations. The most distinctive property of our algorithms is their \emph{remarkable success} in approximating the solution in the shadow region when compared with the algorithms available in the literature.


title: Magnonic Interferometric Switch for Multi-Valued Logic Circuits

abstract: We investigated a possible use of the magnonic interferometric switches in multi-valued logic circuits. The switch is a three-terminal device consisting of two spin channels where input, control, and output signals are spin waves. Signal modulation is achieved via the interference between the source and gate spin waves. We report experimental data on a micrometer scale prototype based on Y3Fe2(FeO4)3 structure. The output characteristics are measured at different angles of the bias magnetic field. The On/Off ratio of the prototype exceeds 36 dB at room temperature. Experimental data is complemented by the theoretical analysis and the results of micro magnetic simulations showing spin wave propagation in a micrometer size magnetic junction. We also present the results of numerical modeling illustrating the operation of a nanometer-size switch consisting of just 20 spins in the source-drain channel. The utilization of spin wave interference as a switching mechanism makes it possible to build nanometer-scale logic gates, and minimize energy per operation, which is limited only by the noise margin. The utilization of phase in addition to amplitude for information encoding offers an innovative route towards multi-state logic circuits. We describe possible implementation of the three-value logic circuits based on the magnonic interferometric switches. The advantages and shortcomings inherent in interferometric switches are also discussed.


title: Measurement of the Yb I $^1S_0 - ^1P_1$ transition frequency at 399 nm using an optical frequency comb

abstract: We determine the frequency of the Yb I $^1S_0 - ^1P_1$ transition at 399 nm using an optical frequency comb. Although this transition was measured previously using an optical transfer cavity [D. Das et al., Phys. Rev. A 72, 032506 (2005)], recent work has uncovered significant errors in that method. We compare our result of 751 526 533.49 $\pm$ 0.33 MHz for the Yb-174 isotope with those from the literature and discuss observed differences. We verify the correctness of our method by measuring the frequencies of well-known transitions in Rb and Cs, and by demonstrating proper control of systematic errors in both laser metrology and atomic spectroscopy. We also demonstrate the effect of quantum interference due to hyperfine structure in a divalent atomic system and present isotope shift measurements for all stable isotopes.


title: Loop expansion and the bosonic representation of loop quantum gravity

abstract: We introduce a new loop expansion that provides a resolution of the identity in the Hilbert space of loop quantum gravity on a fixed graph. We work in the bosonic representation obtained by the canonical quantization of the spinorial formalism. The resolution of the identity gives a tool for implementing the projection of states in the full bosonic representation onto the space of solutions to the Gauss and area matching constraints of loop quantum gravity. This procedure is particularly efficient in the semiclassical regime, leading to explicit expressions for the loop expansions of coherent, heat kernel and squeezed states.


title: Heat Kernel Renormalization on Manifolds with Boundary

abstract: In the monograph Renormalization and Effective Field Theory, Costello made two major advances towards the mathematical formulation of quantum field theory. Firstly, he developed an inductive position space renormalization procedure for constructing effective field theories that is based on heat kernel regularization of the propagator. Secondly, he gave a rigorous formulation of quantum gauge theory within effective field theory that makes use of the BV formalism. In this work, we extend Costello's inductive renormalization procedure from manifolds without boundary to a class of manifolds with boundary. In addition, we reorganize the presentation of the preexisting material, filling in details and strengthening the results.


title: Systems of quasi-variational inequalities related to the switching problem

abstract: We prove the existence of weak solution for a system of quasi-variational inequalities related to a switching problem with dynamic driven by operator associated with a semi-Dirichlet form and with measure data. We give a stochastic representation of solutions in terms of solutions of a system of reflected BSDEs with oblique reflection. As a by-product, we prove the existence of an optimal strategy in the switching problem and show regularity of the payoff function.


title: Stellar nuclei and inner polar disks in lenticular galaxies

abstract: I analyze statistics of the stellar population properties for stellar nuclei and bulges of nearby lenticular galaxies in different environments by using panoramic spectral data of the integral-field spectrograph SAURON retrieved from the open archive of Isaac Newton Group. I estimate also the fraction of nearby lenticular galaxies having inner polar gaseous disks by exploring the volume-limited sample of early-type galaxies of the ATLAS-3D survey. By inspecting the two-dimensional velocity fields of the stellar and gaseous components with running tilted-ring technique, I have found 7 new cases of the inner polar disks. Together with those, the frequency of inner polar disks in nearby S0 galaxies reaches 10% that is much higher than the frequency of large-scale polar rings. Interestingly, the properties of the nuclear stellar populations in the inner polar ring hosts are statistically the same as those in the whole S0 sample implying similar histories of multiple gas accretion events from various directions.


title: Entanglement Entropy in the $\sigma$-Model with the de Sitter Target Space

abstract: We derive the formula of the entanglement entropy between the left and right oscillating modes of the $\sigma$-model with the de Sitter target space. To this end, we study the theory in the \emph{cosmological gauge} in which the non-vanishing components of the metric on the two-dimensional base space are functions of the expansion parameter of the de Sitter space. The model is embedded in the causal north pole diamond of the Penrose diagram. We argue that the cosmological gauge is natural to the $\sigma$-model as it is compatible with the canonical quantization relations. In this gauge, we obtain a new general solution to the equations of motion in terms of time-independent oscillating modes. The constraint structure is adequate for quantization in the Gupta-Bleuler formalism. We construct the space of states as a one-parameter family of Hilbert spaces and give the Bargmann-Fock and Jordan-Schwinger representations of it. Also, we give a simple description of the physical subspace as an infinite product of $\mathcal{D}^{+}_{\frac{1}{2}}$ in the positive discreet series irreducible representations of the $SU(1,1)$ group. We construct the map generated by the Hamiltonian between states at two different values of time and show how it produces the entanglement of left and right excitations. Next, we derive the formula of the entanglement entropy of the reduced density matrix for the ground state acted upon by the Hamiltonian map. Finally, we determine the asymptotic form of the entanglement entropy of a single mode bi-oscillator in the limit of large values of time.


title: Fundamental Limitation on Cooling under Classical Noise

abstract: We prove a general theorem that the action of arbitrary classical noise or random unitary channels can not increase the maximum population of any eigenstate of an open quantum system, assuming initial system-environment factorization. Such factorization is the conventional starting point for descriptions of open system dynamics. In particular, our theorem implies that a system can not be ideally cooled down unless it is initially prepared as a pure state. The resultant inequality rigorously constrains the possibility of cooling the system solely through temporal manipulation, i.e., dynamical control over the system Hamiltonian without resorting to measurement based cooling methods.


title: Quantum lattice gas algorithmic representation of gauge field theory

abstract: Presented is a quantum lattice gas algorithm to efficiently model a system of Dirac particles interacting through an intermediary gauge field. The algorithm uses a fixed qubit array to represent both the spacetime and the particles contained in the spacetime. Despite being a lattice based algorithm, Lorentz invariance is preserved down to the grid scale, with the continuum Dirac Hamiltonian generating the local unitary evolution even at that scale: there is nonlinear scaling between the smallest observable time and that time measured in the quantum field theory limit, a kind of time dilation effect that emerges on small scales but has no effect on large scales. The quantum lattice gas algorithm correctly accounts for the anticommutative braiding of indistinguishable fermions---it does not suffer the Fermi-sign problem. It provides a highly convergent numerical simulation for strongly-correlated fermions equal to a covariant path integral, presented here for the case when a Dirac particle's Compton wavelength is large compared to the grid scale of the qubit array.


title: Fitted Learning: Models with Awareness of their Limits

abstract: Though deep learning has pushed the boundaries of classification forward, in recent years hints of the limits of standard classification have begun to emerge. Problems such as fooling, adding new classes over time, and the need to retrain learning models only for small changes to the original problem all point to a potential shortcoming in the classic classification regime, where a comprehensive a priori knowledge of the possible classes or concepts is critical. Without such knowledge, classifiers misjudge the limits of their knowledge and overgeneralization therefore becomes a serious obstacle to consistent performance. In response to these challenges, this paper extends the classic regime by reframing classification instead with the assumption that concepts present in the training set are only a sample of the hypothetical final set of concepts. To bring learning models into this new paradigm, a novel elaboration of standard architectures called the competitive overcomplete output layer (COOL) neural network is introduced. Experiments demonstrate the effectiveness of COOL by applying it to fooling, separable concept learning, one-class neural networks, and standard classification benchmarks. The results suggest that, unlike conventional classifiers, the amount of generalization in COOL networks can be tuned to match the problem.


title: Optimal ALOHA-like Random Access with Heterogeneous QoS Guarantees for Multi-Packet Reception Aided Visible Light Communications

abstract: There is a paucity of random access protocols designed for alleviating collisions in visible light communication (VLC) systems, where carrier sensing is hard to be achieved due to the directionality of light. To resolve the problem of collisions, we adopt the successive interference cancellation (SIC) algorithm to enable the coordinator to simultaneously communicate with multiple devices, which is referred to as the multi-packet reception (MPR) capability. However, the MPR capability could be fully utilized only when random access algorithms are properly designed. Considering the characteristics of the SIC aided random access VLC system, we propose a novel effective capacity (EC)-based ALOHA-like distributed random access algorithm for MPR-aided uplink VLC systems having heterogeneous quality-of-service (QoS) guarantees. Firstly, we model the VLC network as a conflict graph and derive the EC for each device. Then, we formulate the VLC QoS-guaranteed random access problem as a saturation throughput maximization problem subject to multiple statistical QoS constraints. Finally, the resultant non-concave optimization problem (OP) is solved by a memetic search algorithm relying on invasive weed optimization and differential evolution (IWO-DE). We demonstrate that our derived EC expression matches the Monte Carlo simulation results accurately, and the performance of our proposed algorithms is competitive.


title: Learning to learn with backpropagation of Hebbian plasticity

abstract: Hebbian plasticity is a powerful principle that allows biological brains to learn from their lifetime experience. By contrast, artificial neural networks trained with backpropagation generally have fixed connection weights that do not change once training is complete. While recent methods can endow neural networks with long-term memories, Hebbian plasticity is currently not amenable to gradient descent. Here we derive analytical expressions for activity gradients in neural networks with Hebbian plastic connections. Using these expressions, we can use backpropagation to train not just the baseline weights of the connections, but also their plasticity. As a result, the networks "learn how to learn" in order to solve the problem at hand: the trained networks automatically perform fast learning of unpredictable environmental features during their lifetime, expanding the range of solvable problems. We test the algorithm on various on-line learning tasks, including pattern completion, one-shot learning, and reversal learning. The algorithm successfully learns how to learn the relevant associations from one-shot instruction, and fine-tunes the temporal dynamics of plasticity to allow for continual learning in response to changing environmental parameters. We conclude that backpropagation of Hebbian plasticity offers a powerful model for lifelong learning.


title: Quantum lattice gas model of spin-2 Bose-Einstein condensates and closed-form analytical continuation of nonlinear interactions in spin-2 superfluids

abstract: Presented is an unitary operator splitting method for handling the spin-density interaction in spinor Bose-Einstein condensates. The zero temperature behavior of a spinor BEC is given by mean field theory, where the Hamiltonian includes a nonlinear hyperfine spin interaction. This hyperfine interaction has a diagonal probability-density term (leading to the usual Gross-Pitaevskii type equation of motion) but also has a nondiagonal spin-density term. Since the F=2 spinor BEC (spin-2 BEC) has a non-Abelian superfluid phase (nonperturbative cyclic phase in the strong spin-density coupling regime), an infinite-order expansion of the quantum evolution operator is needed for quantum simulation applications. An infinite-order expansion, obtained by analytical continuation and expressed in analytically closed form, for the spin-2 BEC is presented.


title: Decomposition of radiation energy into work and heat

abstract: We investigate energy transfer by the radiation from a cavity quantum electrodynamics (QED) system in the context of quantum thermodynamics. We propose a method of decomposing it into work and heat within the framework of quantum master equations. We find that the work and heat correspond respectively to the coherent and incoherent parts of the radiation. In the derivation of the method, it is crucial to investigate the dynamics of the system that receives the radiation from the cavity.


title: Cohomological invariants of genus three hyperelliptic curves

abstract: We compute the cohomological invariants with coefficients in $\mathbb{Z}/p\mathbb{Z}$ of the stack $\mathscr{H}_3$ of hyperelliptic curves of genus $3$ over an algebraically closed field.


title: A millisecond pulsar candidate in a 21-hr orbit: 3FGL J0212.1+5320

abstract: We present the discovery of a variable optical counterpart to the unidentified gamma-ray source 3FGL J0212.1+5320, and argue this is a new compact binary millisecond pulsar (MSP) candidate. We show 3FGL J0212.1+5320 hosts a semi-detached binary with a 0.86955$\pm$0.00015 d orbital period and a F6-type companion star at an estimated distance of D=1.1$\pm$0.2 kpc, with a radial velocity curve semi-amplitude K$_2$=214.1$\pm$5.0 km s$^{-1}$ and a projected rotational velocity of Vsin(i)=73.2$\pm$1.6 km s$^{-1}$. We find a hard X-ray source at the same location with a 0.5$-$10 keV luminosity L$_\mathrm{X}$=2.6$\times$10$^{32}$ (D/1.1 kpc)$^2$ erg s$^{-1}$, which strengthens the MSP identification. Our results imply a mass ratio q=M$_2$/M$_1$=0.26$^{+0.02}_{-0.03}$ if the companion star fills its Roche lobe, and q$\gtrsim$0.26 in any case. This classifies 3FGL J0212.1+5320 as a "redback" binary MSP; if its MSP nature is confirmed, this will be the brightest compact binary MSP in the optical band (r'$\simeq$14.3 mag) and will have the longest orbital period among Galactic field systems (nearly 21 hr). Based on the light curve peak-to-peak amplitude ($\Delta$r=0.19 mag), we further suggest that the orbital inclination is high and the putative pulsar mass is close to canonical (M$_1$$\simeq$1.3$-$1.6 M$_\odot$). Finally, we discuss the lack of heating signatures and asymmetric optical light curves in the context of other redback MSPs.


title: Optimization methods for frame conditioning and application to graph Laplacian scaling

abstract: A frame is scalable if each of its vectors can be rescaled in such a way that the resulting set becomes a Parseval frame. In this paper, we consider four different optimization problems for determining if a frame is scalable. We offer some algorithms to solve these problems. We then apply and extend our methods to the problem of reweighing (finite) graph so as to minimize the condition number of the resulting Laplacian.


title: From Physical to Cyber: Escalating Protection for Personalized Auto Insurance

abstract: Nowadays, auto insurance companies set personalized insurance rate based on data gathered directly from their customers' cars. In this paper, we show such a personalized insurance mechanism -- wildly adopted by many auto insurance companies -- is vulnerable to exploit. In particular, we demonstrate that an adversary can leverage off-the-shelf hardware to manipulate the data to the device that collects drivers' habits for insurance rate customization and obtain a fraudulent insurance discount. In response to this type of attack, we also propose a defense mechanism that escalates the protection for insurers' data collection. The main idea of this mechanism is to augment the insurer's data collection device with the ability to gather unforgeable data acquired from the physical world, and then leverage these data to identify manipulated data points. Our defense mechanism leveraged a statistical model built on unmanipulated data and is robust to manipulation methods that are not foreseen previously. We have implemented this defense mechanism as a proof-of-concept prototype and tested its effectiveness in the real world. Our evaluation shows that our defense mechanism exhibits a false positive rate of 0.032 and a false negative rate of 0.013.


title: Explicit Schoen surfaces

abstract: We give an explicit construction for the $4$-dimensional family of Schoen surfaces by computing equations for their canonical images, which are $40$-nodal complete intersections of a quadric and the Igusa quartic in $\mathbb P^4$. We then study a particularly interesting example, with $240$ automorphisms and maximal Picard number.


title: Latent Dependency Forest Models

abstract: Probabilistic modeling is one of the foundations of modern machine learning and artificial intelligence. In this paper, we propose a novel type of probabilistic models named latent dependency forest models (LDFMs). A LDFM models the dependencies between random variables with a forest structure that can change dynamically based on the variable values. It is therefore capable of modeling context-specific independence. We parameterize a LDFM using a first-order non-projective dependency grammar. Learning LDFMs from data can be formulated purely as a parameter learning problem, and hence the difficult problem of model structure learning is circumvented. Our experimental results show that LDFMs are competitive with existing probabilistic models.


title: The unrelaxed dynamical structure of the galaxy cluster Abell 85

abstract: For the first time, we explore the dynamics of the central region of a galaxy cluster within $r_{500}\sim 600h^{-1}$~kpc from its center by combining optical and X-ray spectroscopy. We use (1) the caustic technique that identifies the cluster substructures and their galaxy members with optical spectroscopic data, and (2) the X-ray redshift fitting procedure that estimates the redshift distribution of the intracluster medium (ICM). We use the spatial and redshift distributions of the galaxies and of the X-ray emitting gas to associate the optical substructures to the X-ray regions. When we apply this approach to Abell 85 (A85), a complex dynamical structure of A85 emerges from our analysis: a galaxy group, with redshift $z=0.0509 \pm 0.0021$ is passing through the cluster center along the line of sight dragging part of the ICM present in the cluster core; two additional groups, at redshift $z=0.0547 \pm 0.0022$ and $z=0.0570 \pm 0.0020$, are going through the cluster in opposite directions, almost perpendicularly to the line of sight, and have substantially perturbed the dynamics of the ICM. An additional group in the outskirts of A85, at redshift $z=0.0561 \pm 0.0023$, is associated to a secondary peak of the X-ray emission, at redshift $z=0.0583^{+0.0039}_{-0.0047}$. Although our analysis and results on A85 need to be confirmed by high-resolution spectroscopy, they demonstrate how our new approach can be a powerful tool to constrain the formation history of galaxy clusters by unveiling their central and surrounding structures.


title: Robinson Stability of Parametric Constraint Systems via Variational Analysis

abstract: This paper investigates a well-posedness property of parametric constraint systems named here Robinson stability. Based on advanced tools of variational analysis and generalized differentiation, we derive first-order and second-order conditions for this property under minimal constraint qualifications and establish relationships of Robinson stability with other well-posedness properties in variational analysis and optimization. The results obtained are applied to robust Lipschitzian stability of parametric variational systems.


title: Evaluation of bipartite entanglement between two optical multi-mode systems using mode translation symmetry

abstract: Optical multi-mode systems provide large scale Hilbert spaces that can be accessed and controlled using single photon sources, linear optics and photon detection. Here, we consider the bipartite entanglement generated by coherently distributing $M$ photons in $M$ modes to two separate locations, where linear optics and photon detection is used to verify the non-classical correlations between the two $M$-mode systems. We show that the entangled state is symmetric under mode shift operations performed in the two systems and use this symmetry to derive correlations between photon number distributions detected after a discrete Fourier transform (DFT) of the modes. The experimentally observable correlations can be explained by a simple and intuitive rule that relates the sum of the output mode indices to the eigenvalue of the input state under the mode shift operation. Since the photon number operators after the DFT do not commute with the initial photon number operators, entanglement is necessary to achieve strong correlations in both the initial mode photon numbers and the photon numbers observed after the DFT. We can therefore derive entanglement witnesses based on the experimentally observable correlations in both photon number distributions, providing a practical criterion for the evaluation of large scale entanglement in optical multi-mode systems. Our method thus demonstrates how non-classical signatures in large scale optical quantum circuits can be accessed experimentally by choosing an appropriate combination of modes in which to detect the photon number distributions that characterize the quantum coherences of the state.


title: Slices of Parameter Space of Cubic Polynomials

abstract: In this paper, we study slices of the parameter space of cubic polynomials, up to affine conjugacy, given by a fixed value of the multiplier at a non-repelling fixed point. In particular, we study the location of the $main\, cubioid$ in this parameter space. The $main\, cubioid$ is the set of affine conjugacy classes of complex cubic polynomials that have certain dynamical properties generalizing those of polynomials $z^2+c$ for $c$ in the filled main cardioid.


title: Nodal Variational Principle for Excited States

abstract: It is proven that the exact excited-state wave function and energy may be obtained by minimizing the energy expectation value of trial wave functions that are constrained only to have the correct nodes of the state of interest. This excited-state nodal minimum principle has the advantage that it requires neither minimization with the constraint of wave-function orthogonality to all lower eigenstates nor the antisymmetry of the trial wave functions. It is also found that the minimization over the entire space can be partitioned into several interconnected minimizations within the individual nodal regions, and the exact excited-state energy may be obtained by a minimization in just one or several of these nodal regions. For the proofs of the theorem, it is observed that the many-electron eigenfunction (excited state as well as ground state), restricted to a nodal region, is equivalent to a ground-state wave function of one electron in a higher-dimensional space; and, alternatively, an explicit excited-state energy variational expression is utilized by generalizing the Jacobi method of multiplicative variation. In corollaries, error functions are constructed for cases for which the nodes are not necessarily exact. The exact nodes minimize the energy error functions with respect to nodal variations. Simple numerical illustrations of the error functions are presented.


title: Spin accumulation in disordered topological insulator thin films

abstract: Topological insulator (TI) thin films differ from the more commonly studied semi-infinite bulk TIs in that the former possesses both top and bottom surfaces where the surface states localized at different surfaces can couple to one another due to the finite thickness of the film. In the presence of an in-plane mangnetization TI thin films display two distinct phases depending on which of the inter-surface coupling or the magnetization is stronger. In this work, we consider a TI thin film system with an in-plane magnetization and calculate numerically the resulting spin accumulation on both surfaces of the film due to an in-plane electric field to linear order. We describe a numerical scheme for performing the Kubo calculation calculation in which we include impurity scattering and vertex corrections. We find that the sums of the spin accumulation over the two surfaces in the in-plane direction perpendicular to the magnetization, and in the out of plane direction, are antisymmetric in Fermi energy about the charge neutrality point and are non-vanishing only when the symmetry between the top and bottom TI surfaces is broken. The impurity scattering, in general, diminishes the magnitude of the spin accumulation and can also change the sign of the spin accumulation at some Fermi energies where the accumulation is small.


title: Determination of Pedestrian Flow Performance Based on Video Tracking and Microscopic Simulations

abstract: One of the objectives of understanding pedestrian behavior is to predict the effect of proposed changes in the design or evaluation of pedestrian facilities. We want to know the impact to the user of the facilities, as the design of the facilities change. That impact was traditionally evaluated by level of service standards. Another design criterion to measure the impact of design change is measured by the pedestrian flow performance index. This paper describes the determination of pedestrian flow performance based video tracking or any microscopic pedestrian simulation models. Most of pedestrian researches have been done on a macroscopic level, which is an aggregation of all pedestrian movement in pedestrian areas into flow, average speed and area module. Macroscopic level, however, does not consider the interaction between pedestrians. It is also not well suited for prediction of pedestrian flow performance in pedestrian areas or in buildings with some obstruction, that reduces the effective width of the walkways. On the other hand, the microscopic level has a more general usage and considers detail in the design. More efficient pedestrian flow can even be reached with less space. Those results have rejected the linearity assumption of space and flow in the macroscopic level.


title: On minimality of two-bridge knots

abstract: A knot is called minimal if its knot group admits epimorphisms onto the knot groups of only the trivial knot and itself. In this paper, we determine which two-bridge knot $\mathfrak{b}(p,q)$ is minimal where $q \leq 6$ or $p \leq 100$.


title: Supermassive black holes formed by direct collapse of inflationary perturbations

abstract: We propose a mechanism of producing a new type of primordial perturbations that collapse to primordial black holes whose mass can be as large as necessary for them to grow to the supermassive black holes observed at high redshifts, without contradicting COBE/FIRAS upper limits on cosmic microwave background (CMB) spectral distortions. In our model, the observable Universe consists of two kinds of many small patches which experienced different expansion histories during inflation. Primordial perturbations large enough to form primordial black holes are realized on patches that experienced more Hubble expansion than the others. By making these patches the minor component, the rarity of supermassive black holes can be explained. On the other hand, most regions of the Universe experienced the standard history and, hence, only have standard almost-scale-invariant adiabatic perturbations confirmed by observations of CMB or large-scale structures of the Universe. Thus, our mechanism can evade the constraint from the nondetection of the CMB distortion set by the COBE/FIRAS measurement. Our model predicts the existence of supermassive black holes even at redshifts much higher than those observed. Hence, our model can be tested by future observations peeking into the higher-redshift Universe.


title: A Unified Hamiltonian Solution to Maxwell-Schrodinger Equations for Modeling Electromagnetic Field-Particle Interaction

abstract: A novel unified Hamiltonian approach is proposed to solve Maxwell-Schrodinger equation for modeling the interaction between classical electromagnetic (EM) fields and particles. Based on the Hamiltonian of electromagnetics and quantum mechanics, a unified Maxwell-Schrodinger system is derived by the variational principle. The coupled system is well-posed and symplectic, which ensures energy conserving property during the time evolution. However, due to the disparity of wavelengths of EM waves and that of electron waves, a numerical implementation of the finite-difference time-domain (FDTD) method to the multiscale coupled system is extremely challenging. To overcome this difficulty, a reduced eigenmode expansion technique is first applied to represent the wave function of the particle. Then, a set of ordinary differential equations (ODEs) governing the time evolution of the slowly-varying expansion coefficients are derived to replace the original Schrodinger equation. Finally, Maxwell's equations represented by the vector potential with a Coulomb gauge, together with the ODEs, are solved self-consistently. For numerical examples, the interaction between EM fields and a particle is investigated for both the closed, open and inhomogeneous electromagnetic systems. The proposed approach not only captures the Rabi oscillation phenomenon in the closed cavity but also captures the effects of radiative decay and shift in the open free space. After comparing with the existing theoretical approximate models, it is found that the approximate models break down in certain cases where a rigorous self-consistent approach is needed. This work is helpful for the EM simulation of emerging nanodevices or next-generation quantum electrodynamic systems.

categories: physics.optics math.NA math.SG physics.comp-ph quant-ph


title: Demixing Sines and Spikes: Robust Spectral Super-resolution in the Presence of Outliers

abstract: We consider the problem of super-resolving the line spectrum of a multisinusoidal signal from a finite number of samples, some of which may be completely corrupted. Measurements of this form can be modeled as an additive mixture of a sinusoidal and a sparse component. We propose to demix the two components and super-resolve the spectrum of the multisinusoidal signal by solving a convex program. Our main theoretical result is that-- up to logarithmic factors-- this approach is guaranteed to be successful with high probability for a number of spectral lines that is linear in the number of measurements, even if a constant fraction of the data are outliers. The result holds under the assumption that the phases of the sinusoidal and sparse components are random and the line spectrum satisfies a minimum-separation condition. We show that the method can be implemented via semidefinite programming and explain how to adapt it in the presence of dense perturbations, as well as exploring its connection to atomic-norm denoising. In addition, we propose a fast greedy demixing method which provides good empirical results when coupled with a local nonconvex-optimization step.


title: Investigating Magnetic Activity in the Galactic Centre by Global MHD Simulation

abstract: By performing a global magnetohydrodynamical (MHD) simulation for the Milky Way with an axisymmetric gravitational potential, we propose that spatially dependent amplification of magnetic fields possibly explains the observed noncircular motion of the gas in the Galactic centre (GC) region. The radial distribution of the rotation frequency in the bulge region is not monotonic in general. The amplification of the magnetic field is enhanced in regions with stronger differential rotation, because magnetorotational instability and field-line stretching are more effective. The strength of the amplified magnetic field reaches >~ 0.5 mG, and radial flows of the gas are excited by the inhomogeneous transport of angular momentum through turbulent magnetic field that is amplified in a spatially dependent manner. As a result, the simulated position-velocity diagram exhibits a time-dependent asymmetric parallelogram-shape owing to the intermittency of the magnetic turbulence; the present model provides a viable alternative to the bar-potential-driven model for the parallelogram shape of the central molecular zone. In addition, Parker instability (magnetic buoyancy) creates vertical magnetic structure, which would correspond to observed molecular loops, and frequently excited vertical flows. Furthermore, the time-averaged net gas flow is directed outward, whereas the flows are highly time dependent, which would contribute to the outflow from the bulge.


title: Clustering sequence data with mixture Markov chains with covariates using multiple simplex constrained optimization routine (MSiCOR)

abstract: Mixture Markov Model (MMM) is a widely used tool to cluster sequences of events coming from a finite state-space. However the MMM likelihood being multi-modal, the challenge remains in its maximization. Although Expectation-Maximization (EM) algorithm remains one of the most popular ways to estimate the MMM parameters, however convergence of EM algorithm is not always guaranteed. Given the computational challenges in maximizing the mixture likelihood on the constrained parameter space, we develop a pattern search-based global optimization technique which can optimize any objective function on a collection of simplexes, which is eventually used to maximize MMM likelihood. This is shown to outperform other related global optimization techniques. In simulation experiments, the proposed method is shown to outperform the expectation-maximization (EM) algorithm in the context of MMM estimation performance. The proposed method is applied to cluster Multiple sclerosis (MS) patients based on their treatment sequences of disease-modifying therapies (DMTs). We also propose a novel method to cluster people with MS based on DMT prescriptions and associated clinical features (covariates) using MMM with covariates. Based on the analysis, we divided MS patients into 3 clusters. Further cluster-specific summaries of relevant covariates indicate patient differences among the clusters.


title: On Singer's conjecture for the fifth algebraic transfer

abstract: Let $P_k:= \mathbb{F}_2[x_1,x_2,\ldots ,x_k]$ be the polynomial algebra in $k$ variables with the degree of each $x_i$ being $1,$ regarded as a module over the mod-$2$ Steenrod algebra $\mathcal{A},$ and let $GL_k$ be the general linear group over the prime field $\mathbb{F}_2$ which acts naturally on $P_k$. We study the hit problem, set up by Frank Peterson, of finding a minimal set of generators for the polynomial algebra $P_k$ as a module over the mod-2 Steenrod algebra, $\mathcal{A}$. These results are used to study the Singer algebraic transfer which is a homomorphism from the homology of the mod-$2$ Steenrod algebra, $\mbox{Tor}^{\mathcal{A}}_{k, k+n}(\mathbb{F}_2, \mathbb{F}_2),$ to the subspace of $\mathbb{F}_2\otimes_{\mathcal{A}}P_k$ consisting of all the $GL_k$-invariant classes of degree $n.$ In this paper, we explicitly compute the hit problem for $k = 5$ and the degree $7.2^s-5$ with $s$ an arbitrary positive integer. Using this result, we show that Singer's conjecture for the algebraic transfer is true in the case $k=5$ and the above degree.


title: Characterizations and Effective Computation of Supremal Relatively Observable Sublanguages

abstract: Recently we proposed relative observability for supervisory control of discrete-event systems under partial observation. Relative observability is closed under set unions and hence there exists the supremal relatively observable sublanguage of a given language. In this paper we present a new characterization of relative observability, based on which an operator on languages is proposed whose largest fixpoint is the supremal relatively observable sublanguage. Iteratively applying this operator yields a monotone sequence of languages; exploiting the linguistic concept of support based on Nerode equivalence, we prove for regular languages that the sequence converges finitely to the supremal relatively observable sublanguage, and the operator is effectively computable. Moreover, for the purpose of control, we propose a second operator that in the regular case computes the supremal relatively observable and controllable sublanguage. The computational effectiveness of the operator is demonstrated on a case study.


title: On Performance Modeling for MANETs under General Limited Buffer Constraint

abstract: Understanding the real achievable performance of mobile ad hoc networks (MANETs) under practical network constraints is of great importance for their applications in future highly heterogeneous wireless network environments. This paper explores, for the first time, the performance modeling for MANETs under a general limited buffer constraint, where each network node maintains a limited source buffer of size $B_s$ to store its locally generated packets and also a limited shared relay buffer of size $B_r$ to store relay packets for other nodes. Based on the Queuing theory and birth-death chain theory, we first develop a general theoretical framework to fully depict the source/relay buffer occupancy process in such a MANET, which applies to any distributed MAC protocol and any mobility model that leads to the uniform distribution of nodes' locations in steady state. With the help of this framework, we then derive the exact expressions of several key network performance metrics, including achievable throughput, throughput capacity, and expected end-to-end delay. We further conduct case studies under two network scenarios and provide the corresponding theoretical/simulation results to demonstrate the application as well as the efficiency of our theoretical framework. Finally, we present extensive numerical results to illustrate the impacts of buffer constraint on the performance of a buffer-limited MANET.


title: Distributed Nash equilibrium seeking for aggregative games with coupled constraints

abstract: In this paper, we study a distributed continuous-time design for aggregative games with coupled constraints in order to seek the generalized Nash equilibrium by a group of agents via simple local information exchange. To solve the problem, we propose a distributed algorithm based on projected dynamics and non-smooth tracking dynamics, even for the case when the interaction topology of the multi-agent network is time-varying. Moreover, we prove the convergence of the non-smooth algorithm for the distributed game by taking advantage of its special structure and also combining the techniques of the variational inequality and Lyapunov function.


title: Sisyphus Laser Cooling of a Polyatomic Molecule

abstract: We perform magnetically-assisted Sisyphus laser cooling of the triatomic free radical strontium monohydroxide (SrOH). This is achieved with principal optical cycling in the rotationally closed $P\left(N"=1\right)$ branch of either the $\tilde{X}^{2}\Sigma^{+}\left(000\right)\leftrightarrow\tilde{A}^{2}\Pi_{1/2}\left(000\right)$ or the $\tilde{X}^{2}\Sigma^{+}\left(000\right)\leftrightarrow\tilde{B}^{2}\Sigma^{+}\left(000\right)$ vibronic transitions. Molecules lost into the excited vibrational states during the cooling process are repumped back through the $\tilde{B}\left(000\right)$ state for both the $\left(100\right)$ level of the Sr-O stretching mode and the $\left(02^{0}0\right)$ level of the bending mode. The transverse temperature of a SrOH molecular beam is reduced in one dimension by two orders of magnitude to $\sim700\ {\rm \mu K}$. This approach opens a path towards creating a variety of ultracold polyatomic molecules, including much larger ones, by means of direct laser cooling.


title: The geometry properties of parity and time reversal operators in two dimensional spaces

abstract: The parity operator $\cal P$ and time reversal operator $\cal T$ are two important operators in the quantum theory, in particular, in the $\cal PT$-symmetric quantum theory. By using the concrete forms of $\cal P$ and $\cal T$, we discuss their geometrical properties in two dimensional spaces. It is showed that if $\cal T$ is given, then all $\cal P$ links with the quadric surfaces; if $\cal P$ is given, then all $\cal T$ links with the quadric curves. Moreover, we give out the generalized unbroken $\cal PT$-symmetric condition of an operator. The unbroken $\cal PT$-symmetry of a Hermitian operator is also showed in this way.


title: The 26 December 2001 Solar Event Responsible for GLE63. I. Observations of a Major Long-Duration Flare with the Siberian Solar Radio Telescope

abstract: Ground Level Enhancements (GLEs) of cosmic-ray intensity occur, on average, once a year. Due to their rareness, studying the solar sources of GLEs is especially important to approach understanding their origin. The SOL2001-12-26 eruptive-flare event responsible for GLE63 seems to be challenging in some aspects. Deficient observations limited its understanding. Analysis of extra observations found for this event provided new results shading light on the flare. This article addresses the observations of this flare with the Siberian Solar Radio Telescope (SSRT). Taking advantage of its instrumental characteristics, we analyze the detailed SSRT observations of a major long-duration flare at 5.7 GHz without cleaning the images. The analysis confirms that the source of GLE63 was associated with an event in active region 9742 that comprised two flares. The first flare (04:30-05:03 UT) reached a GOES importance of about M1.6. Two microwave sources were observed, whose brightness temperatures at 5.7 GHz exceeded 10 MK. The main flare, up to the M7.1 importance, started at 05:04 UT, and occurred in strong magnetic fields. The observed microwave sources reached about 250 MK. They were not static. Having appeared on the weaker-field periphery of the active region, the microwave sources moved toward each other nearly along the magnetic neutral line, approaching a stronger-field core of the active region, and then moved away from the neutral line like expanding ribbons. These motions rule out an association of the non-thermal microwave sources with a single flaring loop.


title: Spine decomposition and $L\log L$ criterion for superprocesses with non-local branching mechanisms

abstract: In this paper, we provide a pathwise spine decomposition for superprocesses with both local and non-local branching mechanisms under a martingale change of measure. This result complements the related results obtained in Evans (1993), Kyprianou et al. (2012) and Liu, Ren and Song (2009) for superprocesses with purely local branching mechanisms and in Chen, Ren and Song (2016) and Kyprianou and Palau (2016) for multitype superprocesses. As an application of this decomposition, we obtain necessary/sufficient conditions for the limit of the fundamental martingale to be non-degenerate. In particular, we obtain extinction properties of superprocesses with non-local branching mechanisms as well as a Kesten-Stigum $L\log L$ theorem for the fundamental martingale.


title: Tighter bound of Sketched Generalized Matrix Approximation

abstract: Generalized matrix approximation plays a fundamental role in many machine learning problems, such as CUR decomposition, kernel approximation, and matrix low rank approximation. Especially with today's applications involved in larger and larger dataset, more and more efficient generalized matrix approximation algorithems become a crucially important research issue. In this paper, we find new sketching techniques to reduce the size of the original data matrix to develop new matrix approximation algorithms. Our results derive a much tighter bound for the approximation than previous works: we obtain a $(1+\epsilon)$ approximation ratio with small sketched dimensions which implies a more efficient generalized matrix approximation.


title: Self-triggered Model Predictive Control for Continuous-Time Systems: A Multiple Discretizations Approach

abstract: In this paper, we propose a new self-triggered formulation of Model Predictive Control for continuous-time linear networked control systems. Our control approach, which aims at reducing the number of transmitting control samples to the plant, is derived by parallelly solving optimal control problems with different sampling time intervals. The controller then picks up one sampling pattern as a transmission decision, such that a reduction of communication load and the stability will be obtained. The proposed strategy is illustrated through comparative simulation examples.


title: Spectral and scattering theory for Gauss-Bonnet operators on perturbed topological crystals

abstract: In this paper we investigate the spectral and the scattering theory of Gauss--Bonnet operators acting on perturbed periodic combinatorial graphs. Two types of perturbation are considered: either a multiplication operator by a short-range or a long-range potential, or a short-range type modification of the graph. For short-range perturbations, existence and completeness of local wave operators are proved. In addition, similar results are provided for the Laplacian acting on edges.


title: Three-dimensional network of Drosophila brain hemisphere

abstract: The first step to understanding brain function is to determine the brain's network structure. We report a three-dimensional analysis of the brain network of the fruit fly Drosophila melanogaster by synchrotron-radiation tomographic microscopy. A skeletonized wire model of the left half of the brain network was built by tracing the three-dimensional distribution of X-ray absorption coefficients. The obtained models of neuronal processes were classified into groups on the basis of their three-dimensional structures. These classified groups correspond to neuronal tracts that send long-range projections or repeated structures of the optic lobe. The skeletonized model is also composed of neuronal processes that could not be classified into the groups. The distribution of these unclassified structures correlates with the distribution of contacts between neuronal processes. This suggests that neurons that cannot be classified into typical structures should play important roles in brain functions. The quantitative description of the brain network provides a basis for structural and statistical analyses of the Drosophila brain. The challenge is to establish a methodology for reconstructing the brain network in a higher-resolution image, leading to a comprehensive understanding of the brain structure.


title: Kleitman's conjecture about families of given size minimizing the number of $k$-chains

abstract: A central theorem in combinatorics is Sperner's Theorem, which determines the maximum size of a family $\mathcal{F}\subseteq \mathcal{P}(n)$ that does not contain a $2$-chain $F_1\subsetneq F_2$. Erd\H{o}s later extended this result and determined the largest family not containing a $k$-chain $F_1\subsetneq \ldots \subsetneq F_k$. Erd\H{o}s and Katona and later Kleitman asked how many such chains must appear in families whose size is larger than the corresponding extremal result. This question was resolved for $2$-chains by Kleitman in $1966$, who showed that amongst families of size $M$ in $\mathcal{P}(n)$, the number of $2$-chains is minimized by a family whose sets are taken as close to the middle layer as possible. He also conjectured that the same conclusion should hold for all $k$, not just $2$. The best result on this question is due to Das, Gan and Sudakov who showed that Kleitman's conjecture holds for families whose size is at most the size of the $k+1$ middle layers of $\mathcal{P}(n)$, provided $k\leq n-6$. Our main result is that for every fixed $k$ and $\epsilon>0$, if $n$ is sufficiently large then Kleitman's conjecture holds for families of size at most $(1-\epsilon)2^n$, thereby establishing Kleitman's conjecture asymptotically. Our proof is based on ideas of Kleitman and Das, Gan and Sudakov. Several open problems are also given.


title: X-ray microtomography in biology

abstract: Progress in high-resolution x-ray microtomography has provided us with a practical approach to determining three-dimensional (3D) structures of opaque samples at micrometer to submicrometer resolution. In this review, we give an introduction to hard x-ray microtomography and its application to the visualization of 3D structures of biological soft tissues. Practical aspects of sample preparation, handling, data collection, 3D reconstruction, and structure analysis are described. Furthermore, different sample contrasting methods are approached in detail. Examples of microtomographic studies are overviewed to present an outline of biological applications of x-ray microtomography. We also provide perspectives of biological microtomography as the convergence of sciences in x-ray optics, biology, and structural analysis.


title: A dilaton-pion mass relation

abstract: Recently, Golterman and Shamir presented an effective field theory which is supposed to describe the low-energy physics of the pion and the dilaton in an $SU(N_c)$ gauge theory with $N_f$ Dirac fermions in the fundamental representation. By employing this formulation with a slight but important modification, we derive a relation between the dilaton mass squared~$m_\tau^2$, with and without the fermion mass~$m$, and the pion mass squared~$m_\pi^2$ to the leading order of the chiral logarithm. This is analogous to a similar relation obtained by Matsuzaki and~Yamawaki on the basis of a somewhat different low-energy effective field theory. Our relation reads $m_\tau^2=m_\tau^2|_{m=0}+KN_f\hat{f}_\pi^2m_\pi^2/(2\hat{f}_\tau^2)+O(m_\pi^4\ln m_\pi^2)$ with~$K=9$, where $\hat{f}_\pi$ and~$\hat{f}_\tau$ are decay constants of the pion and the dilaton, respectively. This mass relation differs from the one derived by Matsuzaki and~Yamawaki on the points that $K=(3-\gamma_m)(1+\gamma_m)$, where $\gamma_m$ is the mass anomalous dimension, and the leading chiral logarithm correction is~$O(m_\pi^2\ln m_\pi^2)$. For~$\gamma_m\sim1$, the value of the decay constant~$\hat{f}_\tau$ estimated from our mass relation becomes $\sim50\%$ larger than $\hat{f}_\tau$ estimated from the relation of Matsuzaki and~Yamawaki.


title: Defect production in nonequilibrium phase transitions: Experimental investigation of the Kibble-Zurek mechanism in a two-qubit quantum simulator

abstract: Systems passing through quantum critical points at finite rates have a finite probability of undergoing transitions between different eigenstates of the instantaneous Hamiltonian. This mechanism was proposed by Kibble as the underlying mechanism for the formation of topological defects in the early universe and by Zurek for condensed matter systems. Here, we use a system of nuclear spins as an experimental quantum simulator undergoing a non-equilibrium quantum phase transition. The experimental data confirm the validity of the Kibble-Zurek mechanism of defect formation.


title: Inverse Compton emission from a cosmic-ray precursor in RX J1713.7-3946

abstract: Recently, the High Energy Stereoscopic System (H.E.S.S.) reported two new interesting results for a gamma-ray emitting supernova remnant, RX J1713.7-3946 (G347.3-0.5). The first result is the establishment of a broken power-law spectrum of GeV-TeV gamma-rays. The other is a more extended gamma-ray spatial profile than the one in the X-ray band. In this paper, we show both of these results can be explained by inverse Compton emission from accelerated electrons. If the maximum energy of electrons being accelerated decreases with time, the broken power-law spectrum can be generated by accumulation. Furthermore, the extended component of gamma-ray profile can be interpreted as a CR precursor of currently accelerated electrons.


title: Contact Binaries as Viable Distance Indicators: New, Competitive (V)JHKs Period-Luminosity Relations

abstract: Based on the largest catalogs currently available, comprising 6090 contact binaries (CBs) and 2167 open clusters, we determine the near-infrared $JHK_{\rm s}$ CB period--luminosity (PL) relations, for the first time achieving the low levels of intrinsic scatter that make these relations viable as competitive distance calibrators. To firmly establish our distance calibration on the basis of open cluster CBs, we require that (i) the CB of interest must be located inside the core radius of its host cluster; (ii) the CB's proper motion must be located within the $2\sigma$ distribution of that of its host open cluster; and (iii) the CB's age, $t$, must be comparable to that of its host cluster, i.e., $\Delta \log (t\mbox{ yr}^{-1}) <0.3$. We thus select a calibration sample of 66 CBs with either open cluster distances or accurate space-based parallaxes. The resulting near-infrared PL relations, for both late-type (i.e., W Ursae Majoris-type) and---for the first time---early-type CBs, are as accurate as the well-established $JHK_{\rm s}$ Cepheid PL relations, (characterized by single-band statistical uncertainties of $\sigma < 0.10$ mag). We show that CBs can be used as viable distance tracers, yielding distances with uncertainties of better than 5\% for 90\% of the 6090 CBs in our full sample. By combining the full $JHK_{\rm s}$ photometric data set, CBs can trace distances with an accuracy, $\sigma=0.05 \mbox{ (statistical)} \pm0.03 \mbox{ (systematic)}$ mag. The 102 CBs in the Large Magellanic Cloud are used to determine a distance modulus to the galaxy of $(m-M_V)_0^{\rm LMC}=18.41\pm0.20$ mag.


title: A method for estimating spatial resolution of real image in the Fourier domain

abstract: Spatial resolution is a fundamental parameter in structural sciences. In crystallography, the resolution is determined from the detection limit of high-angle diffraction in reciprocal space. In electron microscopy, correlation in the Fourier domain is used for estimating the resolution. In this paper, we report a method for estimating the spatial resolution of real images from a logarithmic intensity plot in the Fourier domain. The logarithmic intensity plots of test images indicated that the full width at half maximum of a Gaussian point-spread function can be estimated from the images. The spatial resolution of imaging X-ray microtomography using Fresnel zone-plate optics was also estimated with this method. A cross section of a test object visualized with the imaging microtomography indicated that square-wave patterns up to 120-nm pitch were resolved. The logarithmic intensity plot was calculated from a tomographic cross section of brain tissue. The full width at half maximum of the point spread function estimated from the plot coincided with the resolution determined from the test object. These results indicated that the logarithmic intensity plot in the Fourier domain provides an alternative measure of the spatial resolution without explicitly defining a noise criterion.


title: Estimation of presampling modulation transfer function in synchrotron radiation microtomography

abstract: The spatial resolution achieved by recent synchrotron radiation microtomographs should be estimated from the modulation transfer function (MTF) on the micrometer scale. Step response functions of a synchrotron radiation microtomograph were determined by the slanted edge method by using high-precision surfaces of diamond crystal and ion-milled aluminum wire. Tilted reconstruction was introduced to enable any edge to be used as the slanted edge by defining the reconstruction pixel matrix in an arbitrary orientation. MTFs were estimated from the step response functions of the slanted edges. The obtained MTFs coincided with MTF values estimated from square-wave patterns milled on the aluminum surface. Although x-ray refraction influences should be taken into account to evaluate MTFs, any flat surfaces with nanometer roughness can be used to determine the spatial resolutions of microtomographs.


title: Submicrometer tomographic resolution examined using a micro-fabricated test object

abstract: To estimate the spatial resolution of microtomographs, a test object on the submicrometer scale was prepared by focused ion beam milling and subjected to microtomographic analysis. Since human tissues are composed of cells and extracellular matrices with micrometer and submicrometer structures, it is important to investigate the three-dimensional spatial resolution of microtomographs used to visualize microstructures of human tissues. The resolutions along the direction within the tomographic slice plane (in-plane resolution) and perpendicular to it (through-plane resolution) were determined from the modulation transfer function of square-wave patterns. The in-plane resolution was estimated to be 1.2 um from the modulation transfer function of the non-zoom image. In contrast, the zoom image gave the in-plane resolution of 0.8 um. This in-plane resolution is comparable to the through-plane resolution, which was estimated to be 0.8 um. Although the two-dimensional radiographs were taken with the pixel width of half the x-ray optics resolution, these three-dimensional resolution analyses indicated that the zoom reconstruction should be performed to achieve the in-plane resolution comparable to the x-ray optics resolution. The submicrometer three-dimensional analysis was applied in the structural study of human cerebral tissue stained with high-Z elements and the obtained tomograms revealed that the microtomographic analysis allows visualization of the subcellular structures of the cerebral tissue.


title: Ashwin: Plug-and-Play System for Machine-Human Image Annotation

abstract: We present an end-to-end machine-human image annotation system where each component can be attached in a plug-and-play fashion. These components include Feature Extraction, Machine Classifier, Task Sampling and Crowd Consensus.


title: Three-dimensional microtomographic imaging of human brain cortex

abstract: This paper describes an x-ray microtomographic technique for imaging the three-dimensional structure of the human cerebral cortex. Neurons in the brain constitute a neural circuit as a three-dimensional network. The brain tissue is composed of light elements that give little contrast in a hard x-ray transmission image. The contrast was enhanced by staining neural cells with metal compounds. The obtained structure revealed the microarchitecture of the gray and white matter regions of the frontal cortex, which is responsible for the higher brain functions.


title: On the structure of etale motivic cohomology

abstract: We discuss the structure of integral etale motivic cohomology groups of smooth and projective schemes over algebraically closed fields, finite fields, local fields, and arithmetic schemes.


title: Pressure-induced topological phase transition in polar semiconductor BiTeBr

abstract: We performed X-ray diffraction and electrical resistivity measurement up to pressures of 5 GPa and the first-principles calculations utilizing experimental structural parameters to investigate the pressure-induced topological phase transition in BiTeBr having a noncentrosymmetric layered structure (space group P3m1). The P3m1 structure remains stable up to pressures of 5 GPa; the ratio of lattice constants, c/a, has a minimum at pressures of 2.5 - 3 GPa. In the same range, the temperature dependence of resistivity changes from metallic to semiconducting at 3 GPa and has a plateau region between 50 and 150 K in the semiconducting state. Meanwhile, the pressure variation of band structure shows that the bulk band-gap energy closes at 2.9 GPa and re-opens at higher pressures. Furthermore, according to the Wilson loop analysis, the topological nature of electronic states in noncentrosymmetric BiTeBr at 0 and 5 GPa are explicitly revealed to be trivial and non-trivial, respectively. These results strongly suggest that pressure-induced topological phase transition in BiTeBr occurs at the pressures of 2.9 GPa.


title: Uniform Regularity and Vanishing Viscosity limit for the chemotaxis-Navier-Stokes system in a 3D bounded domain

abstract: We investigate the uniform regularity and vanishing viscosity limit for the incompressible chemotaxis-Navier-Stokes system in a smooth bounded domain $\Omega\subset\mathbb{R}^3$. It is shown that there exists a unique strong solution of the incompressible chemotaxis-Navier-Stokes system in a finite time interval which is independent of the viscosity coefficient. Moreover, the solution is uniformly bounded in a conormal Sobolev space, which allows us to take the vanishing viscosity limit to obtain the incompressible inviscid chemotaxis-Navier-Stokes system.


title: Exceptional Points in three-dimensional Nanostructures

abstract: Exceptional points (EPs) are degeneracies in open wave systems where at least two energy levels and their corresponding eigenstates coalesce. We report evidence of the existence of EPs in 3D plasmonic nanostructures. The systems are composed of coupled plasmonic nanoresonators and can be judiciously and systematically driven to EPs by controlling symmetry-compatible modes via their near-field and far-field interactions. The proposed platform opens the way to the investigation of EPs for enhanced light-matter interactions and applications in communication, sensing and imaging.


title: A Decoupling Two-Grid Method for the Steady-State Poisson-Nernst-Planck Equations

abstract: Poisson-Nernst-Planck equations are widely used to describe the electrodiffusion of ions in a solvated biomolecular system. Two kinds of two-grid finite element algorithms are proposed to decouple the steady-state Poisson-Nernst-Planck equations by coarse grid finite element approximations. Both theoretical analysis and numerical experiments show the efficiency and effectiveness of the two-grid algorithms for solving Poisson-Nernst-Planck equations.


title: Infrared-Faint Radio Sources in the SERVS deep fields: Pinpointing AGNs at high redshift

abstract: Infrared-Faint Radio Sources (IFRS) represent an unexpected class of objects relatively bright at radio wavelength, but unusually faint at infrared (IR) and optical wavelengths. A recent and extensive campaign on the radio-brightest IFRSs (S_{1.4GHz} >= 10 mJy) has provided evidence that most of them (if not all) contain an AGN. Still uncertain is the nature of the radio-faintest ones (S_{1.4GHz} <= 1 mJy). The scope of this paper is to assess the nature of the radio-faintest IFRSs, testing their classification and improving the knowledge of their IR properties making use of the most sensitive IR survey available so far: the Spitzer Extragalactic Representative Volume Survey (SERVS). We also explore how the criteria of IFRSs can be fine-tuned to pinpoint radio-loud AGNs at very high redshift (z > 4). We analysed a number of IFRS samples identified in SERVS fields, including a new sample (21 sources) extracted from the Lockman Hole. 3.6 and 4.5 mum IR counterparts of the 64 sources located in the SERVS fields were searched for, and, when detected, their IR properties were studied. We compared the radio/IR properties of the IR-detected IFRSs with those expected for a number of known classes of objects. We found that they are mostly consistent with a mixture of high-redshift (z >= 3) radio-loud AGNs. The faintest ones (S_{1.4GHz} ~ 100 muJy), however, could be also associated with nearer (z ~ 2) dust-enshrouded star-burst galaxies. We also argue that, while IFRSs with radio-to-IR ratios > 500 can very efficiently pinpoint radio-loud AGNs at redshift 2 < z < 4, lower radio-to-IR ratios (~ 100--200) are expected for higher redshift radio-loud AGNs.


title: CERES: A Set of Automated Routines for Echelle Spectra

abstract: We present the Collection of Elemental Routines for Echelle Spectra (CERES). These routines were developed for the construction of automated pipelines for the reduction, extraction and analysis of spectra acquired with different instruments, allowing the obtention of homogeneous and standardised results. This modular code includes tools for handling the different steps of the processing: CCD image reductions, identification and tracing of the echelle orders, optimal and rectangular extraction, computation of the wavelength solution, estimation of radial velocities, and rough and fast estimation of the atmospheric parameters. Currently, CERES has been used to develop automated pipelines for thirteen different spectrographs, namely CORALIE, FEROS, HARPS, ESPaDOnS, FIES, PUCHEROS, FIDEOS, CAFE, DuPont/Echelle, Magellan/Mike, Keck/HIRES, Magellan/PFS and APO/ARCES, but the routines can be easily used in order to deal with data coming from other spectrographs. We show the high precision in radial velocity that CERES achieves for some of these instruments and we briefly summarize some results that have already been obtained using the CERES pipelines.


title: Harmonics Effect on Ion-Bulk Waves in CH Plasmas

abstract: The harmonics effect on ion-bulk (IBk) waves has been researched by Vlasov simulation. The condition of excitation of a large-amplitude IBk waves is given to explain the phenomenon of strong short-wavelength electrostatic activity in solar wind. When $k$ is much lower than $k_{lor}/2$ ($k_{lor}$ is the wave number at loss-of-resonance point), the IBk waves will not be excited to a large amplitude, because a large part of energy will be spread to harmonics. The nature of nonlinear IBk waves in the condition of $k<k_{lor}/2$ is undamped Bernstein-Greene-Kruskal-like waves with harmonics superposition. Only when the wave number $k$ of IBk waves satisfies $k_{lor}/2\lesssim k\leq k_{lor}$, can a large-amplitude and mono-frequency IBk wave be excited. These results give a guidance for a novel scattering mechanism related to IBk waves in the field of laser plasma interaction.


title: Deformable Map Matching for Uncertain Loop-Less Maps

abstract: In the classical context of robotic mapping and localization, map matching is typically defined as the task of finding a rigid transformation (i.e., 3DOF rotation/translation on the 2D moving plane) that aligns the query and reference maps built by mobile robots. This definition is valid in loop-rich trajectories that enable a mapper robot to close many loops, for which precise maps can be assumed. The same cannot be said about the newly emerging autonomous navigation and driving systems, which typically operate in loop-less trajectories that have no large loop (e.g., straight paths). In this paper, we propose a solution that overcomes this limitation by merging the two maps. Our study is motivated by the observation that even when there is no large loop in either the query or reference map, many loops can often be obtained in the merged map. We add two new aspects to map matching: (1) image retrieval with discriminative deep convolutional neural network (DCNN) features, which efficiently generates a small number of good initial alignment hypotheses; and (2) map merge, which jointly deforms the two maps to minimize differences in shape between them. To realize practical computation time, we also present a preemption scheme that avoids excessive evaluation of useless map-matching hypotheses. To verify our approach experimentally, we created a novel collection of uncertain loop-less maps by utilizing the recently published North Campus Long-Term (NCLT) dataset and its ground-truth GPS data. The results obtained using these map collections confirm that our approach improves on previous map-matching approaches.


title: Bell inequality of frequency-bin entangled photon pairs with time-resolved detection

abstract: Entanglement, describing the inseparability of a quantum multiparty system, is one of the most intriguing features of quantum mechanics. Violation of Bell inequality, for ruling out the possibility of local hidden variable theories, is commonly used as a strong witness for quantum entanglement. In previous Bell test experiments with photonic entanglement based on two-photon coincidence measurement, the photon temporal wave packets are absorbed completely by the detectors. That is, the photon coherence time is much shorter than the detection time window. Here we demonstrate generation of frequency-bin entangled narrowband biphotons, and for the first time, test the Clauser-Horne-Shimony-Holt (CHSH) Bell inequality |S|<= 2 for their nonlocal temporal correlations with time-resolved detection. We obtain a maximum |S| value of 2.52+/-0.48 that violates the CHSH inequality. Our result will have applications in quantum information processing involving time-frequency entanglement.