title: A multi-strategy optimizer for arbitrary generic functions in multidimensional space

abstract: An algorithm capable of finding a likely global optimum (minimum) and a set of sub-optimal points for arbitrary generic functions of several variables is presented. The algorithm is designed to deal even with functions of complex behavior, irregular and noisy, with steep variations and exhibiting a lot of local sub-optimal points. The complications of having to deal with a finite domain, as this is usually the case, are taken into account. The method is composed of a number of cascaded stages, each employing a different strategy to improve over the results of the previous stage. Many ideas and concepts employed in known methods are re-elaborated in a coherent scheme, plus several new ideas are introduced. Line minimization plays an important role in most stages, for this purpose a new and powerful algorithm for line minimization is used as well.


title: An algorithm for minimization of arbitrary generic functions in one dimension over a finite domain

abstract: A new algorithm for one-dimensional minimization is described in detail and the results of some tests on practical cases are reported and illustrated. The method requires only punctual computation of the function, and is suitable to be applied in "difficult" cases, that is when the function is highly irregular and has multiple sub-optimal local minima. The algorithm uses quadratic or cubic interpolation and subdivision of intervals in golden ratio as a last resort. It improves over Brent's method and similar ones in several aspects. It manages multiple local minima, takes into account the complications of having to deal with a finite domain, rather than an unlimited one, and has a slightly faster convergence in most cases.


title: Mathematical renormalization in quantum electrodynamics via noncommutative generating series

abstract: In this work, we focus on on the approach by noncommutative formal power series to study the combinatorial aspects of the renormalization at the singularities in $\{0,1,+\infty\}$ of the solutions of nonlinear differential equations involved in quantum electrodynamics.


title: Inconsistency of Measure-Theoretic Probability and Random Behavior of Microscopic Systems

abstract: We report an inconsistency found in probability theory (also referred to as measure-theoretic probability). For probability measures induced by real-valued random variables, we deduce an "equality" such that one side of the "equality" is a probability, but the other side is not. For probability measures induced by extended random variables, we deduce an "equality" such that its two sides are unequal probabilities. The deduced expressions are erroneous only when it can be proved that measure-theoretic probability is a theory free from contradiction. However, such a proof does not exist. The inconsistency appears only in the theory rather than in the physical world, and will not affect practical applications as long as ideal events in the theory (which will not occur physically) are not mistaken for observable events in the real world. Nevertheless, unlike known paradoxes in mathematics, the inconsistency cannot be explained away and hence must be resolved. The assumption of infinite additivity in the theory is relevant to the inconsistency, and may cause confusion of ideal events and real events. As illustrated by an example in this article, since abstract properties of mathematical entities in theoretical thinking are not necessarily properties of physical quantities observed in the real world, mistaking the former for the latter may lead to misinterpreting random phenomena observed in experiments with microscopic systems. Actually the inconsistency is due to the notion of "numbers" adopted in conventional mathematics. A possible way to resolve the inconsistency is to treat "numbers" from the viewpoint of constructive mathematics.


title: Analysis of Low Excitation HDO Transitions Toward the High-Mass Star-forming Regions G34.26+0.15, W51e$_{1}$/e$_{2}$, and W49N

abstract: We present observations of the ground state 1$_{0,1}$-0$_{0,0}$ rotational transition of HDO at 464.925 GHz and the 1$_{1,0}$-1$_{0,1}$ transition at 509.292 GHz toward the three high-mass star forming regions: G34.26+0.15, W49N, and W51e$_{1}$/e$_{2}$, carried out with the Caltech Submillimeter Observatory. The latter transition is observed for the first time from the ground. The spectra are modeled, together with observations of higher-energy HDO transitions, as well as submillimeter dust continuum fluxes from the literature, using a spherically symmetric radiative transfer model to derive the radial distribution of the HDO abundance in the target sources. The abundance profile is divided into an inner hot core region, with kinetic temperatures higher than 100~K, and a cold outer envelope with lower kinetic temperatures. The derived HDO abundance with respect to H$_2$ is (0.3-3.7)$\times 10^{-8}$ in the hot inner region ($T > 100 \mathrm{K}$) and (7.0-10.0)$\times 10^{-11}$ in the cold outer envelope. We also used two H$_{2}^{18}$O fundamental transitions to constrain the H$_{2}$O abundances in the outer envelopes. The HDO/H$_{2}$O ratios in these cold regions are found to be (1.8-3.1)$\times 10^{-3}$ and are consequently higher than in the hot inner regions of these sources.


title: Diameter-Based Active Learning

abstract: To date, the tightest upper and lower-bounds for the active learning of general concept classes have been in terms of a parameter of the learning problem called the splitting index. We provide, for the first time, an efficient algorithm that is able to realize this upper bound, and we empirically demonstrate its good performance.


title: The Warm Circum-Galactic Medium: 10^5-6 K Gas Associated with a Single Galaxy Halo or with an Entire Group of Galaxies?

abstract: In preparation for a Hubble Space Telescope (HST) observing project using the Cosmic Origins Spectrograph (COS), the positions of all AGN targets having high-S/N far-UV G130M spectra were cross-correlated with a large catalog of low-redshift galaxy groups homogeneously selected from the spectroscopic sample of the Sloan Digital Sky Survey (SDSS). Searching for targets behind only those groups at z = 0.1-0.2 (which places the OVI doublet in the wavelength region of peak COS sensitivity) we identified only one potential S/N = 15-20 target, FBQS 1010+3003. An OVI-only absorber was found in its G130M spectrum at z = 0.11326, close to the redshift of a foreground small group of luminous galaxies at z = 0.11685. Because there is no associated Lyalpha absorption, any characterization of this absorber is necessarily minimal; however, the OVI detection likely traces "warm" gas in collisional ionization equilibrium at T ~ 300,000 K. While this discovery is consistent with being interface gas between cooler, photoionized clouds and a hotter intra-group medium, it could also be warm, interface gas associated with the circum-galactic medium (CGM) of the single closest galaxy. In this case a detailed analysis of the galaxy distribution (complete to 0.2 L*) strongly favors the individual galaxy association. This analysis highlights the necessity of both high-S/N > 20 COS data and a deep galaxy redshift survey of the region in order to test more rigorously the association of OVI-absorbing gas with a galaxy group. A Cycle 23 HST/COS program currently is targeting 10 UV-bright AGN behind 12 low-redshift galaxy groups to test the warm, group gas hypothesis.


title: Associated Legendre Functions and Spherical Harmonics of Fractional Degree and Order

abstract: Trigonometric formulas are derived for certain families of associated Legendre functions of fractional degree and order, for use in approximation theory. These functions are algebraic, and when viewed as Gauss hypergeometric functions, belong to types classified by Schwarz, with dihedral, tetrahedral, or octahedral monodromy. The dihedral Legendre functions are expressed in terms of Jacobi polynomials. For the last two monodromy types, an underlying `octahedral' polynomial, indexed by the degree and order and having a non-classical kind of orthogonality, is identified, and recurrences for it are worked out. It is a (generalized) Heun polynomial, not a hypergeometric one. For each of these families of algebraic associated Legendre functions, a representation of the rank-2 Lie algebra so(5,C) is generated by the ladder operators that shift the degree and order of the corresponding solid harmonics. All such representations of so(5,C) are shown to have a common value for each of its two Casimir invariants. The Dirac singleton representations of so(3,2) are included.


title: Comb-like Turing patterns embedded in Hopf oscillations: Spatially localized states outside the 2:1 frequency locked region

abstract: A generic distinct mechanism for the emergence of spatially localized states embedded in an oscillatory background is demonstrated by using 2:1 frequency locking oscillatory system. The localization is of Turing type and appears in two space dimensions as a comb-like state in either $\pi$ phase shifted Hopf oscillations or inside a spiral core. Specifically, the localized states appear in absence of the well known flip-flop dynamics (associated with collapsed homoclinic snaking) that is known to arise in the vicinity of Hopf-Turing bifurcation in one space dimension. Derivation and analysis of three Hopf-Turing amplitude equations in two space dimensions reveals a local dynamics pinning mechanism for Hopf fronts, which in turn allows the emergence of perpendicular (to the Hopf front) Turing states. The results are shown to agree well with the comb-like core size that forms inside spiral waves. In the context of 2:1 resonance, these localized states form outside the 2:1 resonance region and thus extend the frequency locking domain for spatially extended media, such as periodically driven Belousov-Zhabotinsky chemical reactions. Implications to chlorite-iodide-malonic-acid and shaken granular media are also addressed.


title: Multimodal Clustering for Community Detection

abstract: Multimodal clustering is an unsupervised technique for mining interesting patterns in $n$-adic binary relations or $n$-mode networks. Among different types of such generalized patterns one can find biclusters and formal concepts (maximal bicliques) for 2-mode case, triclusters and triconcepts for 3-mode case, closed $n$-sets for $n$-mode case, etc. Object-attribute biclustering (OA-biclustering) for mining large binary datatables (formal contexts or 2-mode networks) arose by the end of the last decade due to intractability of computation problems related to formal concepts; this type of patterns was proposed as a meaningful and scalable approximation of formal concepts. In this paper, our aim is to present recent advance in OA-biclustering and its extensions to mining multi-mode communities in SNA setting. We also discuss connection between clustering coefficients known in SNA community for 1-mode and 2-mode networks and OA-bicluster density, the main quality measure of an OA-bicluster. Our experiments with 2-, 3-, and 4-mode large real-world networks show that this type of patterns is suitable for community detection in multi-mode cases within reasonable time even though the number of corresponding $n$-cliques is still unknown due to computation difficulties. An interpretation of OA-biclusters for 1-mode networks is provided as well.


title: DepthSynth: Real-Time Realistic Synthetic Data Generation from CAD Models for 2.5D Recognition

abstract: Recent progress in computer vision has been dominated by deep neural networks trained over large amounts of labeled data. Collecting such datasets is however a tedious, often impossible task; hence a surge in approaches relying solely on synthetic data for their training. For depth images however, discrepancies with real scans still noticeably affect the end performance. We thus propose an end-to-end framework which simulates the whole mechanism of these devices, generating realistic depth data from 3D models by comprehensively modeling vital factors e.g. sensor noise, material reflectance, surface geometry. Not only does our solution cover a wider range of sensors and achieve more realistic results than previous methods, assessed through extended evaluation, but we go further by measuring the impact on the training of neural networks for various recognition tasks; demonstrating how our pipeline seamlessly integrates such architectures and consistently enhances their performance.


title: Inertial manifolds for 1D reaction-diffusion-advection systems. Part II: periodic boundary conditions

abstract: This is the second part of our study of the Inertial Manifolds for 1D systems of reaction-diffusion-advection equations initiated in \cite{KZI} and it is devoted to the case of periodic boundary conditions. It is shown that, in contrast to the case of Dirichlet or Neumann boundary conditions, considered in the first part, Inertial Manifolds may not exist in the case of systems endowed by periodic boundary conditions. However, as also shown, inertial manifolds still exist in the case of scalar reaction-diffusion-advection equations. Thus, the existence or non-existence of inertial manifolds for this class of dissipative systems strongly depend on the choice of boundary conditions.


title: Estimating the reproductive number, total outbreak size, and reporting rates for Zika epidemics in South and Central America

abstract: As South and Central American countries prepare for increased birth defects from Zika virus outbreaks and plan for mitigation strategies to minimize ongoing and future outbreaks, understanding important characteristics of Zika outbreaks and how they vary across regions is a challenging and important problem. We developed a mathematical model for the 2015 Zika virus outbreak dynamics in Colombia, El Salvador, and Suriname. We fit the model to publicly available data provided by the Pan American Health Organization, using Approximate Bayesian Computation to estimate parameter distributions and provide uncertainty quantification. An important model input is the at-risk susceptible population, which can vary with a number of factors including climate, elevation, population density, and socio-economic status. We informed this initial condition using the highest historically reported dengue incidence modified by the probable dengue reporting rates in the chosen countries. The model indicated that a country-level analysis was not appropriate for Colombia. We then estimated the basic reproduction number, or the expected number of new human infections arising from a single infected human, to range between 4 and 6 for El Salvador and Suriname with a median of 4.3 and 5.3, respectively. We estimated the reporting rate to be around 16% in El Salvador and 18% in Suriname with estimated total outbreak sizes of 73,395 and 21,647 people, respectively. The uncertainty in parameter estimates highlights a need for research and data collection that will better constrain parameter ranges.


title: Interplay of chiral and helical states in a Quantum Spin Hall Insulator lateral junction

abstract: We study the electronic transport across an electrostatically-gated lateral junction in a HgTe quantum well, a canonical 2D topological insulator, with and without applied magnetic field. We control carrier density inside and outside a junction region independently and hence tune the number and nature of 1D edge modes propagating in each of those regions. Outside the 2D gap, magnetic field drives the system to the quantum Hall regime, and chiral states propagate at the edge. In this regime, we observe fractional plateaus which reflect the equilibration between 1D chiral modes across the junction. As carrier density approaches zero in the central region and at moderate fields, we observe oscillations in resistance that we attribute to Fabry-Perot interference in the helical states, enabled by the broken time reversal symmetry. At higher fields, those oscillations disappear, in agreement with the expected absence of helical states when band inversion is lifted.


title: A Kennicutt-Schmidt relation at molecular cloud scales and beyond

abstract: Using N-body/gasdynamic simulations of a Milky Way-like galaxy we analyse a Kennicutt-Schmidt relation, $\Sigma_{SFR} \propto \Sigma_{gas}^N$, at different spatial scales. We simulate synthetic observations in CO lines and UV band. We adopt the star formation rate defined in two ways: based on free fall collapse of a molecular cloud - $\Sigma_{SFR, cl}$, and calculated by using a UV flux calibration - $\Sigma_{SFR, UV}$. We study a KS relation for spatially smoothed maps with effective spatial resolution from molecular cloud scales to several hundred parsecs. We find that for spatially and kinematically resolved molecular clouds the $\Sigma_{SFR, cl} \propto \Sigma_{\rm gas}^N$ relation follows the power-law with index $N \approx 1.4$. Using UV flux as SFR calibrator we confirm a systematic offset between the $\Sigma_{\rm UV}$ and $\Sigma_{\rm gas}$ distributions on scales compared to molecular cloud sizes. Degrading resolution of our simulated maps for surface densities of gas and star formation rates we establish that there is no relation $\Sigma_{\rm SFR, UV} - \Sigma_{\rm gas}$ below the resolution $\sim 50$ pc. We find a transition range around scales $\sim 50-120$ pc, where the power-law index $N$ increases from 0 to 1-1.8 and saturates for scales larger $\sim 120$ pc. A value of the index saturated depends on a surface gas density threshold and it becomes steeper for higher $\Sigma_{gas}$ threshold. Averaging over scales with size of $>150$ pc the power-law index $N$ equals 1.3-1.4 for surface gas density threshold $\sim 5 M_\odot$pc$^{-2}$. At scales $>120$ pc surface SFR densities determined by using CO data and UV flux, $\Sigma_{\rm SFR, UV}/\Sigma_{\rm SFR, cl}$, demonstrate a discrepancy about a factor of 3. We argue that this may be originated from overestimating (constant) values of conversion factor, star formation efficiency or UV calibration used in our analysis.


title: Soft Label Memorization-Generalization for Natural Language Inference

abstract: Often when multiple labels are obtained for a training example it is assumed that there is an element of noise that must be accounted for. It has been shown that this disagreement can be considered signal instead of noise. In this work we investigate using soft labels for training data to improve generalization in machine learning models. However, using soft labels for training Deep Neural Networks (DNNs) is not practical due to the costs involved in obtaining multiple labels for large data sets. We propose soft label memorization-generalization (SLMG), a fine-tuning approach to using soft labels for training DNNs. We assume that differences in labels provided by human annotators represent ambiguity about the true label instead of noise. Experiments with SLMG demonstrate improved generalization performance on the Natural Language Inference (NLI) task. Our experiments show that by injecting a small percentage of soft label training data (0.03% of training set size) we can improve generalization performance over several baselines.


title: Non-Abelian Geometric Phases Carried by the Spin Fluctuation Tensor

abstract: The expectation values of the first and second moments of the quantum mechanical spin operator can be used to define a spin vector and spin fluctuation tensor, respectively. The former is a vector inside the unit ball in three space, while the latter is represented by an ellipsoid in three space. They are both experimentally accessible in many physical systems. By considering transport of the spin vector along loops in the unit ball it is shown that the spin fluctuation tensor picks up geometric phase information. For the physically important case of spin one, the geometric phase is formulated in terms of an SO(3) operator. Loops defined in the unit ball fall into two classes: those which do not pass through the origin and those which pass through the origin. The former class of loops subtend a well defined solid angle at the origin while the latter do not and the corresponding geometric phase is non-Abelian. To deal with both classes, a notion of generalized solid angle is introduced, which helps to clarify the interpretation of the geometric phase information. The experimental systems that can be used to observe this geometric phase are also discussed.


title: Nearly Maximally Predictive Features and Their Dimensions

abstract: Scientific explanation often requires inferring maximally predictive features from a given data set. Unfortunately, the collection of minimal maximally predictive features for most stochastic processes is uncountably infinite. In such cases, one compromises and instead seeks nearly maximally predictive features. Here, we derive upper-bounds on the rates at which the number and the coding cost of nearly maximally predictive features scales with desired predictive power. The rates are determined by the fractal dimensions of a process' mixed-state distribution. These results, in turn, show how widely-used finite-order Markov models can fail as predictors and that mixed-state predictive features offer a substantial improvement.

categories: cond-mat.stat-mech cs.IT math.IT nlin.CD stat.ML


title: Superintegrable classical Zernike system

abstract: We consider the differential equation that Zernike proposed to classify aberrations of wavefronts in a circular pupil, as if it were a classical Hamiltonian with a non-standard potential. The trajectories turn out to be closed ellipses. We show that this is due to the existence of higher-order invariants that close into a cubic Higgs algebra. The Zernike classical system thus belongs to the class of superintegrable systems. Its Hamilton-Jacobi action separates in three vertical projections of polar coordinates of a sphere, polar and equidistant coordinates on half-hyperboloids, and also in elliptic coordinates on the sphere.


title: Optimal Experiment Design for Causal Discovery from Fixed Number of Experiments

abstract: We study the problem of causal structure learning over a set of random variables when the experimenter is allowed to perform at most $M$ experiments in a non-adaptive manner. We consider the optimal learning strategy in terms of minimizing the portions of the structure that remains unknown given the limited number of experiments in both Bayesian and minimax setting. We characterize the theoretical optimal solution and propose an algorithm, which designs the experiments efficiently in terms of time complexity. We show that for bounded degree graphs, in the minimax case and in the Bayesian case with uniform priors, our proposed algorithm is a $\rho$-approximation algorithm, where $\rho$ is independent of the order of the underlying graph. Simulations on both synthetic and real data show that the performance of our algorithm is very close to the optimal solution.


title: eXpose: A Character-Level Convolutional Neural Network with Embeddings For Detecting Malicious URLs, File Paths and Registry Keys

abstract: For years security machine learning research has promised to obviate the need for signature based detection by automatically learning to detect indicators of attack. Unfortunately, this vision hasn't come to fruition: in fact, developing and maintaining today's security machine learning systems can require engineering resources that are comparable to that of signature-based detection systems, due in part to the need to develop and continuously tune the "features" these machine learning systems look at as attacks evolve. Deep learning, a subfield of machine learning, promises to change this by operating on raw input signals and automating the process of feature design and extraction. In this paper we propose the eXpose neural network, which uses a deep learning approach we have developed to take generic, raw short character strings as input (a common case for security inputs, which include artifacts like potentially malicious URLs, file paths, named pipes, named mutexes, and registry keys), and learns to simultaneously extract features and classify using character-level embeddings and convolutional neural network. In addition to completely automating the feature design and extraction process, eXpose outperforms manual feature extraction based baselines on all of the intrusion detection problems we tested it on, yielding a 5%-10% detection rate gain at 0.1% false positive rate compared to these baselines.


title: Probing the Quantum States of a Single Atom Transistor at Microwave Frequencies

abstract: The ability to apply GHz frequencies to control the quantum state of a single $P$ atom is an essential requirement for the fast gate pulsing needed for qubit control in donor based silicon quantum computation. Here we demonstrate this with nanosecond accuracy in an all epitaxial single atom transistor by applying excitation signals at frequencies up to $\approx$ 13 GHz to heavily phosphorous doped silicon leads. These measurements allow the differentiation between the excited states of the single atom and the density of states in the one dimensional leads. Our pulse spectroscopy experiments confirm the presence of an excited state at an energy $\approx$ 9 meV consistent with the first excited state of a single $P$ donor in silicon. The relaxation rate of this first excited state to ground is estimated to be larger than 2.5 GHz, consistent with theoretical predictions. These results represent a systematic investigation of how an atomically precise single atom transistor device behaves under rf excitations.


title: Quantum superintegrable Zernike system

abstract: We consider the differential equation that Zernike proposed to classify aberrations of wavefronts in a circular pupil, whose value at the boundary can be nonzero. On this account the quantum Zernike system, where that differential equation is seen as a Schr\"odinger equation with a potential, is special in that it has a potential and boundary condition that are not standard in quantum mechanics. We project the disk on a half-sphere and there we find that, in addition to polar coordinates, this system separates in two additional coordinate systems (non-orthogonal on the pupil disk), which lead to Schr\"odinger-type equations with P\"oschl-Teller potentials, whose eigen-solutions involve Legendre, Gegenbauer and Jacobi polynomials. This provides new expressions for separated polynomial solutions of the original Zernike system that are real. The operators which provide the separation constants are found to participate in a superintegrable cubic Higgs algebra.


title: Replicating and Scaling up Qualitative Analysis using Crowdsourcing: A Github-based Case Study

abstract: Due to the difficulties in replicating and scaling up qualitative studies, such studies are rarely verified. Accordingly, in this paper, we leverage the advantages of crowdsourcing (low costs, fast speed, scalable workforce) to replicate and scale-up one state-of-the-art qualitative study. That qualitative study explored 20 GitHub pull requests to learn factors that influence the fate of pull requests with respect to approval and merging. As a secondary study, using crowdsourcing at a cost of $200, we studied 250 pull requests from 142 GitHub projects. The prior qualitative findings are mapped into questions for crowds workers. Their answers were converted into binary features to build a predictor which predicts whether code would be merged with median F1 scores of 68%. For the same large group of pull requests, the median F1 scores could achieve 90% by a predictor built with additional features defined by prior quantitative results. Based on this case study, we conclude that there is much benefit in combining different kinds of research methods. While qualitative insights are very useful for finding novel insights, they can be hard to scale or replicate. That said, they can guide and define the goals of scalable secondary studies that use (e.g.) crowdsourcing+data mining. On the other hand, while data mining methods are reproducible and scalable to large data sets, their results may be spectacularly wrong since they lack contextual information. That said, they can be used to test the stability and external validity, of the insights gained from a qualitative analysis.


title: Comparison of Confidence Interval Estimators: an Index Approach

abstract: In many statistical problems, several estimators are usually available for interval estimation of a parameter of interest, and hence, the selection of an appropriate estimator is important. The criterion for a good estimator is to have a high coverage probability close to the nominal level and a shorter interval length. However, these two concepts are in opposition to each other: high and low coverages are associated with longer and shorter interval lengths respectively. Some methods, such as bootstrap calibration, modify the nominal level to improve the coverage and thereby allow the selection of intervals based on interval lengths only. Nonetheless, these methods are computationally expensive. In this paper, we propose an index which offers an easy to compute approach of comparing confidence interval estimators based on a compromise between the coverage probability and the confidence interval length. We illustrate that the confidence interval index has range of values within the neighborhood of the range of the coverage probability, [0,1]. In addition, a good confidence interval estimator has an index value approaching 1; and a bad confidence interval has an index value approaching 0. A simulation study was conducted to assess the finite sample performance of the index. The proposed index is illustrated with a practical example from the literature


title: Dynamics of Bound Magnon Pairs in the Quasi-One-Dimensional Frustrated Magnet LiCuVO_4

abstract: We report on the dynamics of the spin-1/2 quasi-one-dimensional frustrated magnet LiCuVO$\mathrm{_4}$ measured by nuclear spin relaxation in high magnetic fields 10--34 T, in which the ground state has spin-density-wave order. The spin fluctuations in the paramagnetic phase exhibit striking anisotropy with respect to the magnetic field. The transverse excitation spectrum probed by $^{51}$V nuclei has an excitation gap, which increases with field. On the other hand, the gapless longitudinal fluctuations sensed by $^7$Li nuclei grow with lowering temperature, but tend to be suppressed with increasing field. Such anisotropic spin dynamics and its field dependence agree with the theoretical predictions and are ascribed to the formation of bound magnon pairs, a remarkable consequence of the frustration between ferromagnetic nearest neighbor and antiferromagnetic next-nearest-neighbor interactions.


title: Millimeter Wave Beam-Selection Using Out-of-Band Spatial Information

abstract: Millimeter wave (mmWave) communication is one feasible solution for high data-rate applications like vehicular-to-everything communication and next generation cellular communication. Configuring mmWave links, which can be done through channel estimation or beam-selection, however, is a source of significant overhead. In this paper, we propose to use spatial information extracted at sub-6 GHz to help establish the mmWave link. First, we review the prior work on frequency dependent channel behavior and outline a simulation strategy to generate multi-band frequency dependent channels. Second, assuming: (i) narrowband channels and a fully digital architecture at sub-6 GHz; and (ii) wideband frequency selective channels, OFDM signaling, and an analog architecture at mmWave, we outline strategies to incorporate sub-6 GHz spatial information in mmWave compressed beam selection. We formulate compressed beam-selection as a weighted sparse signal recovery problem, and obtain the weighting information from sub-6 GHz channels. In addition, we outline a structured precoder/combiner design to tailor the training to out-of-band information. We also extend the proposed out-of-band aided compressed beam-selection approach to leverage information from all active OFDM subcarriers. The simulation results for achievable rate show that out-of-band aided beam-selection can reduce the training overhead of in-band only beam-selection by 4x.


title: Learning Vector Autoregressive Models with Latent Processes

abstract: We study the problem of learning the support of transition matrix between random processes in a Vector Autoregressive (VAR) model from samples when a subset of the processes are latent. It is well known that ignoring the effect of the latent processes may lead to very different estimates of the influences among observed processes, and we are concerned with identifying the influences among the observed processes, those between the latent ones, and those from the latent to the observed ones. We show that the support of transition matrix among the observed processes and lengths of all latent paths between any two observed processes can be identified successfully under some conditions on the VAR model. From the lengths of latent paths, we reconstruct the latent subgraph (representing the influences among the latent processes) with a minimum number of variables uniquely if its topology is a directed tree. Furthermore, we propose an algorithm that finds all possible minimal latent graphs under some conditions on the lengths of latent paths. Our results apply to both non-Gaussian and Gaussian cases, and experimental results on various synthetic and real-world datasets validate our theoretical results.


title: A Tetrahedron-tiling Method for Crystal Structure Prediction

abstract: Reliable and robust methods of predicting the crystal structure of a compound, based only on its chemical composition, is crucial to the study of materials and their applications. Despite considerable ongoing research efforts, crystal structure prediction remains a challenging problem that demands large computational resources. Here we propose an efficient approach for first-principles crystal structure prediction. The new method explores and finds crystal structures by tiling together elementary tetrahedra that are energetically favorable and geometrically matching each other. This approach has three distinguishing features: a favorable building unit, an efficient calculation of local energy, and a stochastic Monte Carlo simulation of crystal growth. By applying the method to the crystal structure prediction of various materials, we demonstrate its validity and potential as a promising alternative to current methods.


title: Massive Galileon Positivity Bounds

abstract: The EFT coefficients in any gapped, scalar, Lorentz invariant field theory must satisfy positivity requirements if there is to exist a local, analytic Wilsonian UV completion. We apply these bounds to the tree level scattering amplitudes for a massive Galileon. The addition of a mass term, which does not spoil the non-renormalization theorem of the Galileon and preserves the Galileon symmetry at loop level, is necessary to satisfy the lowest order positivity bound. We further show that a careful choice of successively higher derivative corrections are necessary to satisfy the higher order positivity bounds. There is then no obstruction to a local UV completion from considerations of tree level 2-to-2 scattering alone. To demonstrate this we give an explicit example of such a UV completion.


title: Projected Hartree-Fock as a Polynomial of Particle-Hole Excitations and Its Combination With Variational Coupled Cluster Theory

abstract: Projected Hartree-Fock theory provides an accurate description of many kinds of strong correlation but does not properly describe weakly correlated systems. Coupled cluster theory, in contrast, does the opposite. It therefore seems natural to combine the two so as to describe both strong and weak correlations with high accuracy in a relatively black-box manner. Combining the two approaches, however, is made more difficult by the fact that the two techniques are formulated very differently. In earlier work, we showed how to write spin-projected Hartree-Fock in a coupled-cluster-like language. Here, we fill in the gaps in that earlier work. Further, we combine projected Hartree-Fock and coupled cluster theory in a variational formulation and show how the combination performs for the description of the Hubbard Hamiltonian and for several small molecular systems.


title: Maximum Size of a Family of Pairwise Graph-Different Permutations

abstract: Two permutations of the vertices of a graph $G$ are called $G$-different if there exists an index $i$ such that $i$-th entry of the two permutations form an edge in $G$. We bound or determine the maximum size of a family of pairwise $G$-different permutations for various graphs $G$. We show that for all balanced bipartite graphs $G$ of order $n$ with minimum degree $n/2 - o(n)$, the maximum number of pairwise $G$-different permutations of the vertices of $G$ is $2^{(1-o(1))n}$. We also present examples of bipartite graphs $G$ with maximum degree $O(\log n)$ that have this property. We explore the problem of bounding the maximum size of a family of pairwise graph-different permutations when an unlimited number of disjoint vertices is added to a given graph. We determine this exact value for the graph of 2 disjoint edges, and present some asymptotic bounds relating to this value for graphs consisting of the union of $n/2$ disjoint edges.


title: Depth Creates No Bad Local Minima

abstract: In deep learning, \textit{depth}, as well as \textit{nonlinearity}, create non-convex loss surfaces. Then, does depth alone create bad local minima? In this paper, we prove that without nonlinearity, depth alone does not create bad local minima, although it induces non-convex loss surface. Using this insight, we greatly simplify a recently proposed proof to show that all of the local minima of feedforward deep linear neural networks are global minima. Our theoretical results generalize previous results with fewer assumptions, and this analysis provides a method to show similar results beyond square loss in deep linear models.


title: A Linear Time Algorithm for Solving #2SAT on Cactus Formulas

abstract: An O(n+m)-time algorithm is presented for counting the number of models of a two Conjunctive Normal Form Formula F that represents a Cactus graph, where n is the number of variables and m is the number of clauses of F. Although, it was already known that this class of formulas could be computed in polynomial time, we compare our proposal algorithm with two state of the art implementations for the same problem, sharpSAT and countAntom. The results of the comparison show that our algorithm outperforms both implementations, and it can be considered as a base case for general counting of two Conjunctive Normal Form Formulas.


title: Private and Secure Coordination of Match-Making for Heavy-Duty Vehicle Platooning

abstract: A secure and private framework for inter-agent communication and coordination is developed. This allows an agent, in our case a fleet owner, to ask questions or submit queries in an encrypted fashion using semi-homomorphic encryption. The submitted query can be about the interest of the other fleet owners for using a road at a specific time of the day, for instance, for the purpose of collaborative vehicle platooning. The other agents can then provide appropriate responses without knowing the content of the questions or the queries. Strong privacy and security guarantees are provided for the agent who is submitting the queries. It is also shown that the amount of the information that this agent can extract from the other agent is bounded. In fact, with submitting one query, a sophisticated agent can at most extract the answer to two queries. This secure communication platform is used subsequently to develop a distributed coordination mechanisms among fleet owners.


title: Multi-epoch, high spatial resolution observations of multiple T Tauri systems

abstract: Context. In multiple pre-main-sequence systems the lifetime of circumstellar disks appears to be shorter than around single stars, and the actual dissipation process may depend on the binary parameters of the systems. Aims. We report high spatial resolution observations of multiple T Tauri systems at optical and infrared wavelengths. We determine if the components are gravitationally bound and orbital motion is visible, derive orbital parameters and investigate possible correlations between the binary parameters and disk states. Methods. We selected 18 T Tau multiple systems (16 binary and two triple systems, yielding $16 + 2\times2=20$ binary pairs) in the Taurus-Auriga star forming region from the survey by Leinert et al. (1993), with spectral types from K1 to M5 and separations from 0.22" (31 AU) to 5.8" (814 AU). We analysed data acquired in 2006-07 at Calar Alto using the AstraLux lucky imaging system, along with data from SPHERE and NACO at the VLT, and from the literature. Results. We found ten pairs to orbit each other, five pairs that may show orbital motion and five likely common proper motion pairs. We found no obvious correlation between the stellar parameters and binary configuration. The 10 $\mu$m infra-red excess varies between 0.1 and 7.2 magnitudes (similar to the distribution in single stars, where it is between 1.7 and 9.1), implying that the presence of the binary star does not greatly influence the emission from the inner disk. Conclusions. We have detected orbital motion in young T Tauri systems over a timescale of $\approx20$ years. Further observations with even longer temporal baseline will provide crucial information on the dynamics of these young stellar systems.


title: Model-based reinforcement learning in differential graphical games

abstract: This paper seeks to combine differential game theory with the actor-critic-identifier architecture to determine forward-in-time, approximate optimal controllers for formation tracking in multi-agent systems, where the agents have uncertain heterogeneous nonlinear dynamics. A continuous control strategy is proposed, using communication feedback from extended neighbors on a communication topology that has a spanning tree. A model-based reinforcement learning technique is developed to cooperatively control a group of agents to track a trajectory in a desired formation. Simulation results are presented to demonstrate the performance of the developed technique.


title: Surge-like oscillations above sunspot light bridges driven by magnetoacoustic shocks

abstract: High-resolution observations of the solar chromosphere and transition region often reveal surge-like oscillatory activities above sunspot light bridges. These oscillations are often interpreted as intermittent plasma jets produced by quasi-periodic magnetic reconnection. We have analyzed the oscillations above a light bridge in a sunspot using data taken by the Interface Region Imaging Spectrograph (IRIS). The chromospheric 2796\AA{}~images show surge-like activities above the entire light bridge at any time, forming an oscillating wall. Within the wall we often see that the Mg~{\sc{ii}}~k 2796.35\AA{}~line core first experiences a large blueshift, and then gradually decreases to zero shift before increasing to a red shift of comparable magnitude. Such a behavior suggests that the oscillations are highly nonlinear and likely related to shocks. In the 1400\AA{}~passband which samples emission mainly from the Si~{\sc{iv}}~ion, the most prominent feature is a bright oscillatory front ahead of the surges. We find a positive correlation between the acceleration and maximum velocity of the moving front, which is consistent with numerical simulations of upward propagating slow-mode shock waves. The Si~{\sc{iv}} 1402.77\AA{}~line profile is generally enhanced and broadened in the bright front, which might be caused by turbulence generated through compression or by the shocks. These results, together with the fact that the oscillation period stays almost unchanged over a long duration, lead us to propose that the surge-like oscillations above light bridges are caused by shocked p-mode waves leaked from the underlying photosphere.


title: Can Boltzmann Machines Discover Cluster Updates ?

abstract: Boltzmann machines are physics informed generative models with wide applications in machine learning. They can learn the probability distribution from an input dataset and generate new samples accordingly. Applying them back to physics, the Boltzmann machines are ideal recommender systems to accelerate Monte Carlo simulation of physical systems due to their flexibility and effectiveness. More intriguingly, we show that the generative sampling of the Boltzmann Machines can even discover unknown cluster Monte Carlo algorithms. The creative power comes from the latent representation of the Boltzmann machines, which learn to mediate complex interactions and identify clusters of the physical system. We demonstrate these findings with concrete examples of the classical Ising model with and without four spin plaquette interactions. Our results endorse a fresh research paradigm where intelligent machines are designed to create or inspire human discovery of innovative algorithms.

categories: physics.comp-ph cond-mat.stat-mech cs.LG stat.ML


title: Universal many-body response of heavy impurities coupled to a Fermi sea

abstract: In this work we discuss the dynamical response of heavy quantum impurities immersed in a Fermi gas at zero and at finite temperature. Studying both the frequency and the time domain allows one to identify interaction regimes that are characterized by distinct many-body dynamics. From this theoretical study a picture emerges in which impurity dynamics is universal on essentially all time scales, and where the high-frequency few-body response is related to the long-time dynamics of the Anderson orthogonality catastrophe by Tan relations. Our theoretical description relies on different and complementary approaches: functional determinants give an exact numerical solution for time- and frequency-resolved responses, bosonization provides analytical expressions at low temperatures, and the theory of Toeplitz determinants allows one to analytically predict response up to high temperatures. Using these approaches we predict the thermal decoherence rate and prove that within the considered model the fastest rate of long-time decoherence is given by $\gamma=\pi k_BT/4$. We show that Feshbach resonances in cold atomic systems give access to new interaction regimes where quantum effects prevail even in the thermal regime of many-body dynamics. The key signature of this phenomenon is a crossover between exponential decay rates of the real-time Ramsey signal. It is shown that the physics of the orthogonality catastrophe is experimentally observable up to temperatures $T/T_F\lesssim 0.2$ where it leaves its fingerprint in a power-law temperature dependence of thermal spectral weight and we review how this phenomenon is related to the physics of heavy ions in liquid $^3$He and the formation of Fermi polarons. The presented results are in excellent agreement with recent experiments on LiK mixtures, and we predict several phenomena that can be tested using currently available experimental technology.

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


title: Effect of doping on lattice instabilities of single-layer 1H-TaS2

abstract: Recent ARPES measurements of single-layer 1H-TaS2 grown on Au(111) suggest strong electron doping from the substrate. In addition, STM/STS measurements on this system show suppression of the charge-density-wave (CDW) instability that occurs in bulk 2H-TaS2. We present results from ab initio DFT calculations of free-standing single-layer 1H-TaS2 to explore the effects of doping on the CDW. In the harmonic approximation, we find that a lattice instability along the Gamma-M line occurs in the undoped monolayer, consistent with the bulk 3x3 CDW ordering vector. Doping removes the CDW instability, in agreement with the experimental findings. The doping and momentum dependence of both the electron-phonon coupling and of the bare phonon energy (unscreened by metallic electrons) determine the stability of lattice vibrations. Electron doping also causes an expansion of the lattice, so strain is a secondary but also relevant effect.


title: Study of even-parity Rydberg and autoionizing states of lutetium by laser resonance ionization spectroscopy

abstract: Multi-step laser resonance ionization spectroscopy of lutetium (Lu) has been performed at TRIUMF's off-line laser ion source test stand. The even-parity Rydberg series $6s^2nd$ $^2D_{3/2}$, $6s^2nd$ $^2D_{5/2}$ and $6s^2ns$ $^2S_{1/2}$ were observed converging to the 6s$^2$ ionization potential. The experimental results has been compared to previous work. 51 levels of Rydberg series $6s^2nd$ $^2D_{5/2}$ and 52 levels of Rydberg series $6s^2ns$ $^2S_{1/2}$ were reported new. Additionally six even-parity autoionization (AI) series converging to Lu ionic states $5d6s$ $^3D_1$ and $5d6s$ $^3D_2$ were observed. The level energies of these AI states were measured. The configurations of the AI states were assigned by relativistic multichannel theory (RMCT) within the framework of multichannel quantum defect theory (MQDT).


title: Verlinde Gravity and AdS/CFT

abstract: Verlinde argued \cite{Verlinde:2016toy} that baryonic matter in de Sitter Universe with a Hubble constant $H$ produces a bulk contribution to the vacuum entanglement entropy $\delta S_{ent}$. Using the de Sitter temperature $T_{dS}=\frac{H}{2\pi}$, he interpreted $\delta M$ in the first law $T_{dS}\ \delta S_{ent}=-\delta M$ with the Dark Matter energy. We use standard tools of the gauge/gravity correspondence to present some evidence for the Verlinde's proposal.


title: The Shattered Gradients Problem: If resnets are the answer, then what is the question?

abstract: A long-standing obstacle to progress in deep learning is the problem of vanishing and exploding gradients. Although, the problem has largely been overcome via carefully constructed initializations and batch normalization, architectures incorporating skip-connections such as highway and resnets perform much better than standard feedforward architectures despite well-chosen initialization and batch normalization. In this paper, we identify the shattered gradients problem. Specifically, we show that the correlation between gradients in standard feedforward networks decays exponentially with depth resulting in gradients that resemble white noise whereas, in contrast, the gradients in architectures with skip-connections are far more resistant to shattering, decaying sublinearly. Detailed empirical evidence is presented in support of the analysis, on both fully-connected networks and convnets. Finally, we present a new "looks linear" (LL) initialization that prevents shattering, with preliminary experiments showing the new initialization allows to train very deep networks without the addition of skip-connections.


title: Convergence of the age structure of general schemes of population processes

abstract: We consider a family of general branching processes with reproduction parameters depending on the age of the individual as well as the population age structure and a parameter $K$, which may represent the carrying capacity. These processes are Markovian in the age structure. In a previous paper the Law of Large Numbers as $K\to\infty$ was derived. Here we prove the Central Limit Theorem, namely the weak convergence of the fluctuation processes in an appropriate Skorokhod space. We also show that the limit is driven by a stochastic partial differential equation.


title: Mind the Gap: A Study in Global Development through Persistent Homology

abstract: The Gapminder project set out to use statistics to dispel simplistic notions about global development. In the same spirit, we use persistent homology, a technique from computational algebraic topology, to explore the relationship between country development and geography. For each country, four indicators, gross domestic product per capita; average life expectancy; infant mortality; and gross national income per capita, were used to quantify the development. Two analyses were performed. The first considers clusters of the countries based on these indicators, and the second uncovers cycles in the data when combined with geographic border structure. Our analysis is a multi-scale approach that reveals similarities and connections among countries at a variety of levels. We discover localized development patterns that are invisible in standard statistical methods.


title: Sparse Bounds for Spherical Maximal Functions

abstract: We consider the averages of a function $ f$ on $ \mathbb R ^{n}$ over spheres of radius $ 0< r< \infty $ given by $ A_{r} f (x) = \int_{\mathbb S ^{n-1}} f (x-r y) \; d \sigma (y)$, where $ \sigma $ is the normalized rotation invariant measure on $ \mathbb S ^{n-1}$. We prove a sharp range of sparse bounds for two maximal functions, the first the lacunary spherical maximal function, and the second the full maximal function. $$ M_{{lac}} f = \sup_{j\in \mathbb Z } A_{2^j} f , \qquad M_{{full}} f = \sup_{ r>0 } A_{r} f . $$ The sparse bounds are very precise variants of the known $L^p$ bounds for these maximal functions. They are derived from known $ L ^{p}$-improving estimates for the localized versions of these maximal functions, and the indices in our sparse bound are sharp. We derive novel weighted inequalities for weights in the intersection of certain Muckenhoupt and reverse H\"older classes.


title: Tunnel transport and interlayer excitons in bilayer fractional quantum Hall systems

abstract: In a bilayer system consisting of a composite-fermion Fermi sea in each layer, the tunnel current is exponentially suppressed at zero bias, followed by a strong peak at a finite bias voltage $V_{\rm max}$. This behavior, which is qualitatively different from that observed for the electron Fermi sea, provides fundamental insight into the strongly correlated non-Fermi liquid nature of the CF Fermi sea and, in particular, offers a window into the short-distance high-energy physics of this state. We identify the exciton responsible for the peak current and provide a quantitative account of the value of $V_{\rm max}$. The excitonic attraction is shown to be quantitatively significant, and its variation accounts for the increase of $V_{\rm max}$ with the application of an in-plane magnetic field. We also estimate the critical Zeeman energy where transition occurs from a fully spin polarized composite fermion Fermi sea to a partially spin polarized one, carefully incorporating corrections due to finite width and Landau level mixing, and find it to be in satisfactory agreement with the Zeeman energy where a qualitative change has been observed for the onset bias voltage [Eisenstein et al., Phys. Rev. B 94, 125409 (2016)]. For fractional quantum Hall states, we predict a substantial discontinuous jump in $V_{\rm max}$ when the system undergoes a transition from a fully spin polarized state to a spin singlet or a partially spin polarized state.


title: The interplay between Dirac systems, Morse families and Interconnection

abstract: We provide a generalization of the notion of Dirac system by using Morse families to intrinsically embrace the dynamics associated with different physical systems such as constrained variational calculus, optimal control, Lagrangian mechanics, etc., from both the Lagrangian and the Hamiltonian viewpoint. Besides, we introduce a geometric construction of Port-Dirac systems and use it to describe the interconnection of Dirac systems. Special emphasis on theoretical and applied examples is given throughout the paper to show the applicability of the proposed framework.


title: Enabling Sparse Winograd Convolution by Native Pruning

abstract: Sparse methods and the use of Winograd convolutions are two orthogonal approaches, each of which significantly accelerates convolution computations in modern CNNs. Sparse Winograd merges these two and thus has the potential to offer a combined performance benefit. Nevertheless, training convolution layers so that the resulting Winograd kernels are sparse has not hitherto been very successful. By introducing a Winograd layer in place of a standard convolution layer, we can learn and prune Winograd coefficients "natively" and obtain sparsity level beyond 90% with only 0.1% accuracy loss with AlexNet on ImageNet dataset. Furthermore, we present a sparse Winograd convolution algorithm and implementation that exploits the sparsity, achieving up to 31.7 effective TFLOP/s in 32-bit precision on a latest Intel Xeon CPU, which corresponds to a 5.4x speedup over a state-of-the-art dense convolution implementation.


title: Market-Driven Energy Storage Planning for Microgrids with Renewable Energy Systems Using Stochastic Programming

abstract: Battery Energy Storage Systems (BESS) can mitigate effects of intermittent energy production from renewable energy sources and play a critical role in peak shaving and demand charge management. To optimally size the BESS from an economic perspective, the trade-off between BESS investment costs, lifetime, and revenue from utility bill savings along with microgrid ancillary services must be taken into account. The optimal size of a BESS is solved via a stochastic optimization problem considering wholesale market pricing. A stochastic model is used to schedule arbitrage services for energy storage based on the forecasted energy market pricing while accounting for BESS cost trends, the variability of renewable energy resources, and demand prediction. The uniqueness of the approach proposed in this paper lies in the convex optimization programming framework that computes a globally optimal solution to the financial trade-off solution. The approach is illustrated by application to various realistic case studies based on pricing and demand data from the California Independent System Operator (CAISO). The case study results give insight in optimal BESS sizing from a cost perspective, based on both yearly scheduling and daily BESS operation.


title: Some Families of Super Congruences Involving Alternating Multiple Harmonic Sums

abstract: Let $p$ be a prime. In this short note we study some families of super congruences involving the following alternating sums \begin{equation*} \sum_{\substack{j_1+j_2+\cdots+j_n=2 p^r p\nmid j_1 j_2 \cdots j_n}} \frac{(-1)^{j_1+\cdots+j_b}}{j_1\cdots j_n} \pmod{p^r}, \end{equation*} which extend similar statements proved by Shen and Cai who treated the cases when $n=4,5$.


title: Casimir energy for acoustic phonons in graphene

abstract: We find the Casimir energy, at finite temperature, for acoustic phonons in a Graphene sheet suspended over a rectangular trench, and the corresponding Casimir forces are interpreted as correction terms to the built-in tensions of the Graphene. We show that these corrections generally break the tensional isotropy of the membrane, and can increase or decrease the membrane tension. We demonstrate that for a narrow rectangular trench with side-lengths in the order of few nanometers and few micrometers, these temperature corrections are expected to be noticeable ($\sim 10^{-4} N/m$) at the room temperature. These corrections would be even more considerable by increasing the temperature, and can be applied for adjusting the built-in tension of the Graphene. Consequently we introduce a corrected version for the fundamental resonance frequency of the Graphene resonator.


title: Parallel Structure from Motion from Local Increment to Global Averaging

abstract: In this paper, we tackle the accurate and consistent Structure from Motion (SfM) problem, in particular camera registration, far exceeding the memory of a single computer in parallel. Different from the previous methods which drastically simplify the parameters of SfM and sacrifice the accuracy of the final reconstruction, we try to preserve the connectivities among cameras by proposing a camera clustering algorithm to divide a large SfM problem into smaller sub-problems in terms of camera clusters with overlapping. We then exploit a hybrid formulation that applies the relative poses from local incremental SfM into a global motion averaging framework and produce accurate and consistent global camera poses. Our scalable formulation in terms of camera clusters is highly applicable to the whole SfM pipeline including track generation, local SfM, 3D point triangulation and bundle adjustment. We are even able to reconstruct the camera poses of a city-scale data-set containing more than one million high-resolution images with superior accuracy and robustness evaluated on benchmark, Internet, and sequential data-sets.


title: Observations of M31 and M33 with the Fermi Large Area Telescope: a galactic center excess in Andromeda?

abstract: The Fermi LAT has opened the way for comparative studies of cosmic rays (CRs) and high-energy objects in the Milky Way (MW) and in other, external, star-forming galaxies. Using 2 yr of observations with the Fermi LAT, local Group galaxy M31 was detected as a marginally extended gamma-ray source, while only an upper limit (UL) has been derived for the other nearby galaxy M33. We revisited the gamma-ray emission in the direction of M31 and M33 using more than 7 yr of LAT Pass 8 data in the energy range 0.1-100 GeV, presenting detailed morphological and spectral analyses. M33 remains undetected and we computed an UL for it. This revised UL remains consistent with the observed correlation between gamma-ray luminosity and star-formation rate tracers and implies an average CR density in M33 that is at most half of that of the MW. M31 is detected with a significance of nearly 10 sigma and to be extended with 4 sigma. Its spectrum is consistent with a power law. The spatial distribution of the emission is consistent with a uniform disk with a radius of 0.4 deg and no offset from the center of M31, but nonuniform intensity distributions cannot be excluded. The flux from M31 appears confined to the inner regions of the galaxy and does not fill the disk or extend far from it. The gamma-ray signal is not correlated with regions rich in gas or star-formation activity suggesting that the emission is not interstellar in origin, unless the energetic particles radiating in gamma rays do not originate in recent star formation. Alternative and nonexclusive interpretations are that the emission results from a population of millisecond pulsars dispersed in the bulge and disk of M31 by disrupted globular clusters or from the decay or annihilation of dark matter particles, similar to what has been proposed to account for the so-called Galactic Center excess found in Fermi-LAT observations of the MW.


title: Approximation by translates of a single function of functions in space induced by the convolution with a given function

abstract: We study approximation by arbitrary linear combinations of $n$ translates of a single function of periodic functions. We construct some methods of this approximation for functions in a class induced by the convolution with a given function, and prove upper bounds of $L_p$-the approximation convergence rate by these methods, when $n \to \infty$, for $1 < p < \infty$, and lower bounds of the quantity of best approximation of this class by arbitrary linear combinations of $n$ translates of arbitrary function, for the particular case $p=2$.


title: Evaluating Optical Classification for {\em Fermi} Blazar Candidates with a Statistical method using Broadband Spectral Indices

abstract: We aim to test if a blazar candidate of uncertain-type (BCU) in the third Fermi active galactic nuclei catalog (3LAC) can be potentially classified as a BL Lac object or a flat spectrum radio quasar (FSRQ) by performing a statistical analysis of its broadband spectral properties. We find that 34% of the radio-selected BCUs (583 BCUs) are BL Lac-like and 20% of them are FSRQ-like, which maybe within 90% level of confidence. Similarly, 77.3% of the X-ray selected BCUs (176 BCUs) are evaluated as BL Lac-like and 6.8% of them may be FSRQ-like sources. And 88.7% of the BL Lac-like BCUs that have synchrotron peak frequencies available are high synchrotron peaked BL Lacs in the X-ray selected BCUs. The percentages are accordingly 62% and 7.3% in the sample of 124 optical-selected BCUs. The high ratio of source numbers of the BL Lac-like to the FSRQ-like BCUs in the X-ray and optically selected BCU samples is due to the selection effect. Examining the consistency between our evaluation and spectroscopic identification case by case with a sample of 78 radio-selected BCUs, it is found that the statistical analysis and its resulting classifications agree with the results of the optical follow-up spectroscopic observations. Our observation campaign for high-|\rho_{s}| BCUs , i.e., |\rho_{s}|>0.8, selected with our method is ongoing.


title: Effective field theory description of halo nuclei

abstract: Nuclear halos emerge as new degrees of freedom near the neutron and proton driplines. They consist of a core and one or a few nucleons which spend most of their time in the classically-forbidden region outside the range of the interaction. Individual nucleons inside the core are thus unresolved in the halo configuration, and the low-energy effective interactions are short-range forces between the core and the valence nucleons. Similar phenomena occur in clusters of $^4$He atoms, cold atomic gases near a Feshbach resonance, and some exotic hadrons. In these weakly-bound quantum systems universal scaling laws for s-wave binding emerge that are independent of the details of the interaction. Effective field theory (EFT) exposes these correlations and permits the calculation of non-universal corrections to them due to short-distance effects, as well as the extension of these ideas to systems involving the Coulomb interaction and/or binding in higher angular-momentum channels. Halo nuclei exhibit all these features. Halo EFT, the EFT for halo nuclei, has been used to compute the properties of single-neutron, two-neutron, and single-proton halos of s-wave and p-wave type. This review summarizes these results for halo binding energies, radii, Coulomb dissociation, and radiative capture, as well as the connection of these properties to scattering parameters, thereby elucidating the universal correlations between all these observables. We also discuss how Halo EFT's encoding of the long-distance physics of halo nuclei can be used to check and extend ab initio calculations that include detailed modeling of their short-distance dynamics.


title: The Active Atlas: Combining 3D Anatomical Models with Texture Detectors

abstract: While modern imaging technologies such as fMRI have opened exciting new possibilities for studying the brain in vivo, histological sections remain the best way to study the anatomy of the brain at the level of single neurons. The histological atlas changed little since 1909 and localizing brain regions is a still a labor intensive process performed only by experienced neuro-anatomists. Existing digital atlases such as the Allen Brain atlas are limited to low resolution images which cannot identify the detailed structure of the neurons. We have developed a digital atlas methodology that combines information about the 3D organization of the brain and the detailed texture of neurons in different structures. Using the methodology we developed an atlas for the mouse brainstem and mid-brain, two regions for which there are currently no good atlases. Our atlas is "active" in that it can be used to automatically align a histological stack to the atlas, thus reducing the work of the neuroanatomist.


title: Faster DB-scan and HDB-scan in Low-Dimensional Euclidean Spaces

abstract: We present a new algorithm for the widely used density-based clustering method DBscan. Our algorithm computes the DBscan-clustering in $O(n\log n)$ time in $\mathbb{R}^2$, irrespective of the scale parameter $\varepsilon$ (and assuming the second parameter MinPts is set to a fixed constant, as is the case in practice). Experiments show that the new algorithm is not only fast in theory, but that a slightly simplified version is competitive in practice and much less sensitive to the choice of $\varepsilon$ than the original DBscan algorithm. We also present an $O(n\log n)$ randomized algorithm for HDBscan in the plane---HDBscan is a hierarchical version of DBscan introduced recently---and we show how to compute an approximate version of HDBscan in near-linear time in any fixed dimension.


title: Towards A Rigorous Science of Interpretable Machine Learning

abstract: As machine learning systems become ubiquitous, there has been a surge of interest in interpretable machine learning: systems that provide explanation for their outputs. These explanations are often used to qualitatively assess other criteria such as safety or non-discrimination. However, despite the interest in interpretability, there is very little consensus on what interpretable machine learning is and how it should be measured. In this position paper, we first define interpretability and describe when interpretability is needed (and when it is not). Next, we suggest a taxonomy for rigorous evaluation and expose open questions towards a more rigorous science of interpretable machine learning.


title: Effects of extended components of \oiii on the correlation between the \oiii luminosity and the power-law continuum luminosity for active galactic nuclei

abstract: In this manuscript, we check the well-known correlation between \oiii luminosity and continuum luminosity (\loiiicon) for AGN by a large sample of 1982 SDSS QSOs with $z<0.8$ and with high quality spectra. The strong correlation of \loiiicon can be found, similar as previous results for AGN. Moreover, among the 1982 QSOs, there are 708 QSOs with the [O~{\sc iii}]$\lambda$5007\AA\ described by two components: one core component plus one extended component. Based on the luminosity from the core components ($L_{\rm [O~\textsc{iii}],~narrow}$) and from the extended components ($L_{\rm [O~\textsc{iii}],~ext}$), we confirm that the correlation of $L_{{\rm [O~\textsc{iii}],~ext}}~-~L_{{\rm 5100\textsc{\AA}}}$ is more stronger and tighter than the correlations on the total \oiii luminosity and on the luminosity of the core components of the \oiii lines. Therefore, the luminosity of the extended components should be better applied to trace AGN intrinsic luminosity. Meanwhile, we have found strong line width correlation and line luminosity correlation between the core components and the extended components, indicating the extended components of the \oiii lines should be not due to commonly considered radial flows in the common \oiii line clouds. And virial effects due to gravity of central black holes naturally lead to the wider extended components from regions more nearer to central black holes. Finally, we can say that the reported correlation of $L_{{\rm [O~\textsc{iii}],~narrow}}~-~L_{{\rm 5100\textsc{\AA}}}$ on the core components of the \oiii lines should be more better to estimate AGN intrinsic luminosity in Type-2 narrow line AGN, because of totally/partly obscured extended components.


title: Cosmology in beyond-generalized Proca theories

abstract: The beyond-generalized Proca theories are the extension of second-order massive vector-tensor theories (dubbed generalized Proca theories) with two transverse vector modes and one longitudinal scalar besides two tensor polarizations. Even with this extension, the propagating degrees of freedom remain unchanged on the isotropic cosmological background without an Ostrogradski instability. We study the cosmology in beyond-generalized Proca theories by paying particular attention to the dynamics of late-time cosmic acceleration and resulting observational consequences. We derive conditions for avoiding ghosts and instabilities of tensor, vector, and scalar perturbations and discuss viable parameter spaces in concrete models allowing the dark energy equation of state smaller than $-1$. The propagation speeds of those perturbations are subject to modifications beyond the domain of generalized Proca theories. There is a mixing between scalar and matter sound speeds, but such a mixing is suppressed during most of the cosmic expansion history without causing a new instability. On the other hand, we find that derivative interactions arising in beyond-generalized Proca theories give rise to important modifications to the cosmic growth history. The growth rate of matter perturbations can be compatible with the redshift-space distortion data due to the realization of gravitational interaction weaker than that in generalized Proca theories. Thus, it is possible to distinguish the dark energy model in beyond-generalized Proca theories from the counterpart in generalized Proca theories as well as from the $\Lambda$CDM model.


title: Titan's Atmosphere and Climate

abstract: Titan is the only moon with a substantial atmosphere, the only other thick N$_{2}$ atmosphere besides Earth's, the site of extraordinarily complex atmospheric chemistry that far surpasses any other solar system atmosphere, and the only other solar system body with stable liquid currently on its surface. The connection between Titan's surface and atmosphere is also unique in our Solar system; atmospheric chemistry produces materials that are deposited on the surface and subsequently altered by surface-atmosphere interactions such as aeolian and fluvial processes resulting in the formation of extensive dune fields and expansive lakes and seas. Titan's atmosphere is favorable for organic haze formation, which combined with the presence of some oxygen bearing molecules indicates that Titan's atmosphere may produce molecules of prebiotic interest. The combination of organics and liquid, in the form of water in a subsurface ocean and methane/ethane in the surface lakes and seas, means that Titan may be the ideal place in the solar system to test ideas about habitability, prebiotic chemistry, and the ubiquity and diversity of life in the Universe. The Cassini-Huygens mission to the Saturn system has provided a wealth of new information allowing for study of Titan as a complex system. Here I review our current understanding of Titan's atmosphere and climate forged from the powerful combination of Earth-based observations, remote sensing and \emph{in situ} spacecraft measurements, laboratory experiments, and models. I conclude with some of our remaining unanswered questions as the incredible era of exploration with Cassini-Huygens comes to an end.


title: Mather measures for space-time periodical nonconvex hamiltonians

abstract: The main goal of this paper is to construct Mather measures for space-time periodical nonconvex Hamiltonians.


title: Extending the LHC Reach for New Physics with Sub-Millimeter Displaced Vertices

abstract: Particles with a sub-millimeter decay length appear in many models of physics beyond the Standard Model. However, their longevity has been often ignored in their LHC searches and they have been regarded as promptly-decaying particles. In this letter, we show that, by requiring displaced vertices on top of the event selection criteria used in the ordinary search strategies for promptly-decaying particles, we can considerably extend the LHC reach for particles with a decay length of $\gtrsim 100~\mu{\rm m}$. We discuss a way of reconstructing sub-millimeter displaced vertices by exploiting the same technique used for the primary vertex reconstruction on the assumption that the metastable particles are always pair-produced and their decay products contain high-$p_{\rm T}$ jets. We show that, by applying a cut based on displaced vertices on top of standard kinematical cuts for the search of new particles, the LHC reach can be significantly extended if the decay length is $\gtrsim 100~\mu{\rm m}$. In addition, we may measure the lifetime of the target particle through the reconstruction of displaced vertices, which plays an important role in understanding the new physics behind the metastable particles.


title: Testing redMaPPer centring probabilities using galaxy clustering and galaxy-galaxy lensing

abstract: Galaxy cluster centring is a key issue for precision cosmology studies using galaxy surveys. Mis-identification of central galaxies causes systematics in various studies such as cluster lensing, satellite kinematics, and galaxy clustering. The red-sequence Matched-filter Probabilistic Percolation (redMaPPer) estimates the probability that each member galaxy is central from photometric information rather than specifying one central galaxy. The redMaPPer estimates can be used for calibrating the off-centring effect, however, the centring algorithm has not previously been well-tested. We test the centring probabilities of redMaPPer cluster catalog using the projected cross correlation between redMaPPer clusters with photometric red galaxies and galaxy-galaxy lensing. We focus on the subsample of redMaPPer clusters in which the redMaPPer central galaxies (RMCGs) are not the brightest member galaxies (BMEM) and both of them have spectroscopic redshift. This subsample represents nearly 10% of the whole cluster sample. We find a clear difference in the cross-correlation measurements between RMCGs and BMEMs, and the estimated centring probability is 74$\pm$10% for RMCGs and 13$\pm$4% for BMEMs in the Gaussian offset model and 78$\pm$9% for RMCGs and 5$\pm$5% for BMEMs in the NFW offset model. These values are in agreement with the centring probability values reported by redMaPPer (75% for RMCG and 10% for BMEMs) within 1$\sigma$. Our analysis provides a strong consistency test of the redMaPPer centring probabilities. Our results suggest that redMaPPer centring probabilities are reliably estimated. We confirm that the brightest galaxy in the cluster is not always the central galaxy as has been shown in previous works.


title: Bridging Finite and Super Population Causal Inference

abstract: There are two general views in causal analysis of experimental data: the super population view that the units are an independent sample from some hypothetical infinite populations, and the finite population view that the potential outcomes of the experimental units are fixed and the randomness comes solely from the physical randomization of the treatment assignment. These two views differs conceptually and mathematically, resulting in different sampling variances of the usual difference-in-means estimator of the average causal effect. Practically, however, these two views result in identical variance estimators. By recalling a variance decomposition and exploiting a completeness-type argument, we establish a connection between these two views in completely randomized experiments. This alternative formulation could serve as a template for bridging finite and super population causal inference in other scenarios.


title: Complete Classification of Two-Dimensional Algebras

abstract: A complete classification of two-dimensional algebras over algebraically closed fields is provided


title: Spontaneous Activity in the Visual Cortex is Organized by Visual Streams

abstract: Large-scale functional networks have been extensively studied using resting state functional magnetic resonance imaging. However, the pattern, organization, and function of fine-scale network activity remain largely unknown. Here we characterized the spontaneously emerging visual cortical activity by applying independent component analysis to resting state fMRI signals exclusively within the visual cortex. In this sub-system scale, we observed about 50 spatially independent components that were reproducible within and across subjects, and analyzed their spatial patterns and temporal relationships to reveal the intrinsic parcellation and organization of the visual cortex. We found that the visual cortical parcels were aligned with the steepest gradient of cortical myelination, and organized into functional modules segregated along the dorsal/ventral pathways and foveal/peripheral early visual areas. In contrast, cortical retinotopy, folding, and cytoarchitecture impose limited constraints to the organization of resting state activity. From these findings, we conclude that spontaneous activity patterns in the visual cortex are primarily organized by visual streams, likely reflecting feedback network interactions.


title: A Volcanic Hydrogen Habitable Zone

abstract: The classical habitable zone is the circular region around a star in which liquid water could exist on the surface of a rocky planet. The outer edge of the traditional N2-CO2-H2O habitable zone (HZ) extends out to nearly 1.7 AU in our Solar System, beyond which condensation and scattering by CO2 outstrips its greenhouse capacity. Here, we show that volcanic outgassing of atmospheric H2 on a planet near the outer edge can extend the habitable zone out to ~2.4 AU in our solar system. This wider volcanic hydrogen habitable zone (N2-CO2-H2O-H2) can be sustained as long as volcanic H2 output offsets its escape from the top of the atmosphere. We use a single-column radiative-convective climate model to compute the HZ limits of this volcanic hydrogen habitable zone for hydrogen concentrations between 1% and 50%, assuming diffusion-limited atmospheric escape. At a hydrogen concentration of 50%, the effective stellar flux required to support the outer edge decreases by ~35% to 60% for M to A stars. The corresponding orbital distances increase by ~30% to 60%. The inner edge of this HZ only moves out by ~0.1 to 4% relative to the classical HZ because H2 warming is reduced in dense H2O atmospheres. The atmospheric scale heights of such volcanic H2 atmospheres near the outer edge of the HZ also increase, facilitating remote detection of atmospheric signatures.


title: Spinning Witten Diagrams

abstract: We develop a systematic framework to compute the conformal partial wave expansions (CPWEs) of tree-level four-point Witten diagrams with totally symmetric external fields of arbitrary mass and integer spin in AdS$_{d+1}$. Employing this framework, we determine the CPWE of a generic exchange Witten diagram with spinning exchanged field. As an intermediate step, we diagonalise the linear map between spinning three-point conformal structures and spinning cubic couplings in AdS. As a concrete application, we compute all exchange diagrams in the type A higher-spin gauge theory on AdS$_{d+1}$, which is conjectured to be dual to the free scalar $O\left(N\right)$ model. Given a CFT$_d$, our results provide the complete holographic reconstruction of all cubic couplings involving totally symmetric fields in the putative dual theory on AdS$_{d+1}$.


title: An efficient approach to suppress the negative role of contrarian oscillators in synchronization

abstract: It has been found that contrarian oscillators usually take a negative role in the collective behaviors formed by conformist oscillators. However, experiments revealed that it is also possible to achieve a strong coherence even when there are contrarians in the system such as neuron networks with both excitable and inhibitory neurons. To understand the underlying mechanism of this abnormal phenomenon, we here consider a complex network of coupled Kuramoto oscillators with mixed positive and negative couplings and present an efficient approach, i.e. tit-for-tat strategy, to suppress the negative role of contrarian oscillators in synchronization and thus increase the order parameter of synchronization. Two classes of contrarian oscillators are numerically studied and a brief theoretical analysis is provided to explain the numerical results.


title: The Bressoud-G\"ollnitz-Gordon Theorem for Overpartitions of even moduli

abstract: We give an overpartition analogue of Bressoud's combinatorial generalization of the G\"ollnitz-Gordon theorem for even moduli in general case. Let $\widetilde{O}_{k,i}(n)$ be the number of overpartitions of $n$ whose parts satisfy certain difference condition and $\widetilde{P}_{k,i}(n)$ be the number of overpartitions of $n$ whose non-overlined parts satisfy certain congruence condition. We show that $\widetilde{O}_{k,i}(n)=\widetilde{P}_{k,i}(n)$ for $1\leq i<k$.


title: Weighted Calder\'{o}n-Zygmund estimates for weak solutions of quasi-linear degenerate elliptic equations

abstract: This paper studies the Sobolev regularity estimates for weak solutions of a class of degenerate, and singular quasi-linear elliptic problems of the form $\text{div}[\mathbf{A}(x,u, \nabla u)]= \text{div}[\mathbf{F}]$ with non-homogeneous Dirichlet boundary conditions over bounded non-smooth domains. The coefficients $\mathbf{A}$ could be be singular, and degenerate or both in $x$ in the sense that they behave like some weight function $\mu$, which is in the $A_2$ class of Muckenhoupt weights. Global and interior weighted $W^{1,p}(\Omega, \omega)$-regularity estimates are established for weak solutions of these equations with some other weight function $\omega$. The results obtained are even new for the case $\mu =1$ because of the dependence on the solution $u$ of $\mathbf{A}$. In case of linear equations, our $W^{1,p}$-regularity estimates can be viewed as the Sobolev's counterpart of the H\"{o}lder's regularity estimates established by B. Fabes, C. E. Kenig, and R. P. Serapioni.


title: Progress Estimation and Phase Detection for Sequential Processes

abstract: Process modeling and understanding are fundamental for advanced human-computer interfaces and automation systems. Most recent research has focused on activity recognition, but little has been done on sensor-based detection of process progress. We introduce a real-time, sensor-based system for modeling, recognizing and estimating the progress of a work process. We implemented a multimodal deep learning structure to extract the relevant spatio-temporal features from multiple sensory inputs and used a novel deep regression structure for overall completeness estimation. Using process completeness estimation with a Gaussian mixture model, our system can predict the phase for sequential processes. The performance speed, calculated using completeness estimation, allows online estimation of the remaining time. To train our system, we introduced a novel rectified hyperbolic tangent (rtanh) activation function and conditional loss. Our system was tested on data obtained from the medical process (trauma resuscitation) and sports events (Olympic swimming competition). Our system outperformed the existing trauma-resuscitation phase detectors with a phase detection accuracy of over 86%, an F1-score of 0.67, a completeness estimation error of under 12.6%, and a remaining-time estimation error of less than 7.5 minutes. For the Olympic swimming dataset, our system achieved an accuracy of 88%, an F1-score of 0.58, a completeness estimation error of 6.3% and a remaining-time estimation error of 2.9 minutes.


title: Calibration of the Large Area X-ray Proportional Counter (LAXPC) instrument on-board AstroSat

abstract: We present the calibration and background model for the Large Area X-ray Proportional Counter (LAXPC) detectors on-board AstroSat. LAXPC instrument has three nominally identical detectors to achieve large collecting area. These detectors are independent of each other and in the event analysis mode, they record the arrival time and energy of each photon that is detected. The detectors have a time-resolution of 10 $\mu$s and a dead-time of about 42 $\mu$s. This makes LAXPC ideal for timing studies. The energy resolution and peak channel to energy mapping were obtained from calibration on ground using radioactive sources coupled with GEANT4 simulations of the detectors. The response matrix was further refined from observations of the Crab X-ray source after launch. At around 20 keV the energy resolution of detector is about 10--15\%, while the combined effective area of the 3 detectors is about 6000 cm$^2$.


title: Fluctuation of Dynamical Robustness in a Networked Oscillators System

abstract: In this work, we study the dynamical robustness in a system consisting of both active and inactive oscillators. We analytically show that the dynamical robustness of such system is determined by the cross link density between active and inactive subpopulations, which depends on the specific process of inactivation. It is the multi-valued dependence of the cross link density on the control parameter, i.e., the ratio of inactive oscillators in the system, that leads to the fluctuation of the critical points. We further investigate how different network topologies and inactivation strategies affect the fluctuation. Our results explain why the fluctuation is more obvious in heterogeneous networks than in homogeneous ones, and why the low-degree nodes are crucial in terms of dynamical robustness. The analytical results are supported by numerical verifications.


title: Show, Attend and Interact: Perceivable Human-Robot Social Interaction through Neural Attention Q-Network

abstract: For a safe, natural and effective human-robot social interaction, it is essential to develop a system that allows a robot to demonstrate the perceivable responsive behaviors to complex human behaviors. We introduce the Multimodal Deep Attention Recurrent Q-Network using which the robot exhibits human-like social interaction skills after 14 days of interacting with people in an uncontrolled real world. Each and every day during the 14 days, the system gathered robot interaction experiences with people through a hit-and-trial method and then trained the MDARQN on these experiences using end-to-end reinforcement learning approach. The results of interaction based learning indicate that the robot has learned to respond to complex human behaviors in a perceivable and socially acceptable manner.


title: An Optimization Framework with Flexible Inexact Inner Iterations for Nonconvex and Nonsmooth Programming

abstract: In recent years, numerous vision and learning tasks have been (re)formulated as nonconvex and nonsmooth programmings(NNPs). Although some algorithms have been proposed for particular problems, designing fast and flexible optimization schemes with theoretical guarantee is a challenging task for general NNPs. It has been investigated that performing inexact inner iterations often benefit to special applications case by case, but their convergence behaviors are still unclear. Motivated by these practical experiences, this paper designs a novel algorithmic framework, named inexact proximal alternating direction method (IPAD) for solving general NNPs. We demonstrate that any numerical algorithms can be incorporated into IPAD for solving subproblems and the convergence of the resulting hybrid schemes can be consistently guaranteed by a series of simple error conditions. Beyond the guarantee in theory, numerical experiments on both synthesized and real-world data further demonstrate the superiority and flexibility of our IPAD framework for practical use.


title: Analysis of Agent Expertise in Ms. Pac-Man using Value-of-Information-based Policies

abstract: Conventional reinforcement learning methods for Markov decision processes rely on weakly-guided, stochastic searches to drive the learning process. It can therefore be difficult to predict what agent behaviors might emerge. In this paper, we consider an information-theoretic cost function for performing constrained stochastic searches that promote the formation of risk-averse to risk-favoring behaviors. This cost function is the value of information, which provides the optimal trade-off between the expected return of a policy and the policy's complexity; policy complexity is measured by number of bits and controlled by a single hyperparameter on the cost function. As the policy complexity is reduced, the agents will increasingly eschew risky actions. This reduces the potential for high accrued rewards. As the policy complexity increases, the agents will take actions, regardless of the risk, that can raise the long-term rewards. The obtainable reward depends on a single, tunable hyperparameter that regulates the degree of policy complexity. We evaluate the performance of value-of-information-based policies on a stochastic version of Ms. Pac-Man. A major component of this paper is the demonstration that ranges of policy complexity values yield different game-play styles and explaining why this occurs. We also show that our reinforcement-learning search mechanism is more efficient than the others we utilize. This result implies that the value of information theory is appropriate for framing the exploitation-exploration trade-off in reinforcement learning.


title: Nontrivial standing wave state in frequency-weighted Kuramoto model

abstract: Synchronization in a frequency-weighted Kuramoto model with a uniform frequency distribution is studied. We plot the bifurcation diagram and identify the asymptotic coherent states. Numerical simulations show that the system undergoes two first-order transitions in both the forward and backward directions. Apart from the trivial phase-locked state, a novel nonstationary coherent state, i.e., a nontrivial standing wave state is observed and characterized. In this state, oscillators inside the coherent clusters are not frequency-locked as they would be in the usual standing wave state. Instead, their average frequencies are locked to a constant. The critical coupling strength from the incoherent state to the nontrivial standing wave state can be obtained by performing linear stability analysis. The theoretical results are supported by the numerical simulations.


title: A quantum degenerate Bose-Fermi mixture of $^{41}$K and $^6$Li

abstract: We report a new apparatus for the study of two-species quantum degenerate mixture of $^{41}$K and $^6$Li atoms. We develop and combine several advanced cooling techniques to achieve both large atom number and high phase space density of the two-species atom clouds. Furthermore, we build a high-efficiency two-species magnetic transport system to transfer atom clouds from the 3D magneto-optical-trap chamber to a full glass science chamber of extreme high vacuum environment and good optical access. We perform a forced radio-frequency evaporative cooling for $^{41}$K atoms while the $^6$Li atoms are sympathetically cooled in an optically-plugged magnetic trap. Finally, we achieve the simultaneous quantum degeneracy for the $^{41}$K and $^6$Li atoms. The Bose-Einstein condensate of $^{41}$K has 1.4$\times$10$^5$ atoms with a condensate fraction of about 62%, while the degenerate Fermi gas of $^6$Li has a total atom number of 5.4$\times$10$^5$ at 0.25 Fermi temperature.


title: Topological recursion with hard edges

abstract: We prove a Givental type decomposition for partition functions that arise out of topological recursion applied to spectral curves. Copies of the Konstevich-Witten KdV tau function arise out of regular spectral curves and copies of the Brezin-Gross-Witten KdV tau function arise out of irregular spectral curves. We present the example of this decomposition for the matrix model with two hard edges and spectral curve $(x^2-4)y^2=1$


title: Critical behaviour in two-dimensional Coulomb Glass at zero temperature

abstract: The lattice model of Coulomb Glass in two dimensions with box-type random field distribution is studied at zero temperature for system size upto $96^{2}$. To obtain the minimum energy state we annealed the system using Monte Carlo simulation followed by further minimization using cluster-flipping. The values of the critical exponents are determined using the standard finite size scaling. We found that the correlation length $\xi$ diverges with an exponent $\nu=1.0$ at the critical disorder $W_{c} = 0.2253$ and that $\chi_{dis} \approx \xi^{4-\bar{\eta}}$ with $\bar{\eta}=2$ for the disconnected susceptibility. The staggered magnetization behaves discontinuously around the transition and the critical exponent of magnetization $\beta=0$. The probability distribution of the staggered magnetization shows a three peak structure which is a characteristic feature for the phase coexistence at first-order phase transition. In addition to this, at the critical disorder we have also studied the properties of the domain for different system sizes. In contradiction with the Imry-Ma arguments, we found pinned and non-compact domains where most of the random field energy was contained in the domain wall. Our results are also inconsistent with Binder's roughening picture.


title: Molecular clouds in the extreme outer Galaxy between $l=34.75^{\circ}$ to $45.25^{\circ}$

abstract: We present the result of an unbiased CO survey in Galactic range of 34.75$^{\circ}\leq$l$\leq$ 45.25$^{\circ}$ and -5.25$^{\circ}\leq$b$\leq$ 5.25$^{\circ}$, and velocity range beyond the Outer arm. A total of 168 molecular clouds (MCs) are identified within the Extreme Outer Galaxy~(EOG) region, and 31 of these MCs are associated with 13CO emission. However, none of them show significant C18O emission under current detection limit. The typical size and mass of these MCs are 5~pc and 3$\times$10$^{3}M_{\sun}$, implying the lack of large and massive MCs in the EOG region. Similar to MCs in the outer Galaxy, the velocity dispersions of EOG clouds are also correlated with their sizes, however, are well displaced below the scaling relationship defined by the inner Galaxy MCs. These MCs with a median Galactocentric radius of 12.6 kpc, show very different distributions from those of the MCs in the Outer arm published in our previous paper, while roughly follow the Outer Scutum-Centaurus arm defined by Dame & Thaddeus 2011. This result may provide a robust evidence for the existence of the Outer Scutum-Centaurus arm. The lower limit of the total mass of this segment is about 2.7$\times$10$^5$ $M_{\sun}$, which is about one magnitude lower than that of the Outer arm. The mean thickness of gaseous disk is about 1.45$^{\circ}$ or 450~pc, and the scale height is about 1.27$^{\circ}$ or 400~pc above the $b=0^{\circ}$ plane. The warp traced by CO emission is very obvious in the EOG region and its amplitude is consistent with the predictions by other warp models using different tracers, such as dust, HI and stellar components of our Galaxy.