title: Turbulent channel flow over an anisotropic porous wall - Drag increase and reduction

abstract: The effect of the variations of the permeability tensor on the close-to-the-wall behaviour of a turbulent channel flow bounded by porous walls is explored using a set of direct numerical simulations. It is found that the total drag can be either reduced or increased by more than $20\%$ by adjusting the permeability directional properties. Drag reduction is achieved for the case of materials with permeability in the vertical direction lower than the one in the wall-parallel planes. This configuration limits the wall normal velocity at the interface while promoting an increase of the tangential slip velocity leading to an almost "one-component" turbulence where the low- and high-speed streaks coherence is strongly enhanced. On the other hand, strong drag increase is found when a high wall-normal and low wall-parallel permeabilities are prescribed. In this condition, the enhancement of the wall-normal fluctuations due to the reduced wall-blocking effect triggers the onset of structures which are strongly correlated in the spanwise direction, a phenomenon observed by other authors in flows over isotropic porous layers or over ribletted walls with large protrusion heights. The use of anisotropic porous walls for drag reduction is particularly attractive since equal gains can be achieved at different Reynolds numbers by rescaling the magnitude of the permeability only.


title: A van Benthem Theorem for Fuzzy Modal Logic

abstract: We present a fuzzy (or quantitative) version of the van Benthem theorem, which characterizes propositional modal logic as the bisimulation-invariant fragment of first-order logic. Specifically, we consider a first-order fuzzy predicate logic along with its modal fragment, and show that the fuzzy first-order formulas that are non-expansive w.r.t. the natural notion of bisimulation distance are exactly those that can be approximated by fuzzy modal formulas.


title: Exploring hyperon structure with electromagnetic transverse densities

abstract: We explore the structure of the spin-1/2 flavor-octet baryons (hyperons) through their electromagnetic transverse densities. The transverse densities describe the distribution of charge and magnetization at fixed light-front time and enable a spatial representation of the baryons as relativistic systems. At peripheral distances b~1/M_pi the transverse densities are computed using a new method that combines chiral effective field theory and dispersion analysis. The peripheral isovector densities arise from two-pion exchange, which includes the rho-meson resonance through elastic unitarity. The isoscalar densities are estimated from vector meson exchange (omega, phi). We find that the "pion cloud" in the charged Sigma hyperons is comparable to the nucleon, while in the Xi it is suppressed. The Lambda-Sigma^0 transition density is pure isovector and represents a clear manifestiation of peripheral two-pion dynamics.


title: Manifestation of Superposition and Coherence in $\cal PT$-symmetry through the $\eta$-inner Product

abstract: Through the $\eta$-inner product, we investigate $\cal PT$-symmetric quantum mechanics from the viewpoint of superposition and coherence theory. It is argued that $\cal PT$-symmetric quantum systems are endowed with stationary superposition or superposition-free properties. A physical interpretation of $\eta$-inner product in $\mathbb C^2$ is given through the Stokes parameters, showing the difference between broken and unbroken $\cal PT$-symmetric quantum systems.


title: Presqu'un immeuble pour le groupe des automorphismes mod\'er\'es

abstract: Inspired by the Bruhat-Tits building of SL$_n$($\mathbb Q_p$), we construct a complete metric space X with an action of the tame automorphism group of the affine space Tame($K^n$). The points in X are certain monomial valuations, and X admits a natural structure of Euclidean CW-complex of dimension n-1. When n = 3, and for K of characteristic zero, we prove that X has non-positive curvature and is simply connected, hence is a CAT(0) space. As an application we obtain the linearizability of finite subgroups in Tame($K^3$).


title: MALS-NOT: Identifying Radio-Bright Quasars for the MeerKAT Absorption Line Survey

abstract: We present a preparatory spectroscopic survey to identify radio-bright, high-redshift quasars for the MeerKAT Absorption Line Survey (MALS). The candidates have been selected on the basis of a single flux density limit at 1.4 GHz (>200 mJy) together with mid-infrared color criteria from the Wide-field Infrared Survey Explorer (WISE). Through spectroscopic observations using the Nordic Optical Telescope, we identify 72 quasars out of 99 candidates targeted. We measure the spectroscopic redshifts based on characteristic, broad emission lines present in the spectra. Of these 72 quasars, 64 and 48 objects are at sufficiently high redshift (z>0.6 and z>1.4) to be used for the L-band and UHF-band spectroscopic follow-up with the Square Kilometre Array (SKA) precursor in South Africa: the MeerKAT.


title: Completing the Structural Analysis of the 2x4 Permutation Classes

abstract: We study the structure and enumeration of the final two 2x4 permutation classes, completing a research program that has spanned almost two decades. For both classes, careful structural analysis produces a complicated functional equation. One of these equations is solved with the guess-and-check paradigm, while the other is solved with kernel method-like techniques and Gr\"obner basis calculations.


title: Alternative Spreadsheet Model Designs for an Operations Management Model Embedded in a Periodic Business Process

abstract: We present a widely-used operations management model used in supply and distribution planning, that is typically embedded in a periodic business process that necessitates model modification and reuse. We consider three alternative spreadsheet implementations, a data-driven design, a canonical (textbook) design, and a novel (table-driven) technical design. We evaluate each regarding suitability for accuracy, modification, analysis, and transfer. We consider the degree of training and technical sophistication required to utilize each design. The data-driven design provides insight into poor spreadsheet practices by na\"ive modelers. The technical design can be modified for new data and new structural elements without manual writing or editing of cell formulas, thus speeding modification and reducing risk of error. The technical design has potential for use with other classes of models. We identify opportunities for future research.


title: Investigation of the Electronical and Optical Properties of Quantum Well Lasers with Slightly Doped Tunnel Junction

abstract: We experimentally investigate and analyze the electrical and optical characteristics of GaAs-based conventional quantum well laser diodes and the quantum well laser diodes with slightly-doped tunnel junction. It was found that TJ LD show a nonlinear S-shape I-V characteristic.It was also found that the internal quantum efficiency measured by 21% and 87.3% for the TJ LD and the conventional LD, respectively. Furthermore, compared with the conventional LD, it was found that we could achieve 15 nm broadband spectrum from the TJ LD due to lasing dynamics reflects the current dynamics. The results may also lead to the realization of more applications.


title: ALMA [CI] observations toward the central region of a Seyfert galaxy NGC 613

abstract: We report ALMA observations of [CI]($^3P_1-^3P_0$), C$^{13}$O, and CO$^{18}$ ($J=1-0$) toward the central region of a nearby Seyfert galaxy NGC 613. The very high resolutions of $0.26"\times0.23"(=22\times20$ pc) for [CI] and $0.42"\times0.35"(=36\times30$ pc) for C$^{13}$O, and CO$^{18}$ resolve the circum-nuclear disk (CND) and star-forming ring. The distribution of [CI] in the ring resembles that of the CO emission, although [CI] is prominent in the CND. This can be caused by the low intensities of the CO isotopes due to the low optical depths under the high temperature in the CND. We found that the intensity ratios of [CI] to C$^{12}$O(3-2) ($R_{\rm CI/CO}$) and to C$^{13}$O(1-0) ($R_{\rm CI/C^{13}O}$) are high at several positions around the edge of the ring. The spectral profiles of CO lines mostly correspond each other in the spots of the ring and high $R_{\rm CI/CO}$, but those of [CI] at spots of high $R_{\rm CI/CO}$ are different from CO. These results indicate that [CI] at the high $R_{\rm CI/CO}$ traces different gas from that traced by the CO lines. The [CI] kinematics along the minor axis of NGC 613 could be interpreted as a bubbly molecular outflow. The outflow rate of molecular gas is higher than star formation rate in the CND. The flow could be mainly boosted by the AGN through its radio jets.


title: Krasovskii-Subbotin approach to mean field type differential games

abstract: A mean field type differential game is a mathematical model of a large system of identical agents under mean-field interaction controlled by two players with opposite purposes. We study the case when the dynamics of each agent is given by ODE and the players can observe the distribution of the agents. We construct suboptimal strategies and prove the existence of the value function.


title: Prime spectra of abelian 2-categories and categorifications of Richardson varieties

abstract: We describe a general framework for prime, completely prime, semiprime, and primitive ideals of an abelian 2-category. This provides a noncommutative version of Balmer's prime spectrum of a tensor triangulated category. These notions are based on containment conditions in terms of thick subcategories of an abelian category and thick ideals of an abelian 2-category. We prove categorical analogs of the main properties of noncommutative prime spectra. Similar notions, starting with Serre subcategories of an abelian category and Serre ideals of an abelian 2-category, are developed. They are linked to Serre prime spectra of $\mathbb{Z}_+$-rings. As an application, we construct a categorification of the quantized coordinate rings of open Richardson varieties for symmetric Kac-Moody groups, by constructing Serre completely prime ideals of monoidal categories of modules of the KLR algebras, and by taking Serre quotients with respect to them.


title: Optimizing the signal-to-noise ratio of biphoton distribution measurements

abstract: Single-photon-sensitive cameras can now be used as massively parallel coincidence counters for entangled photon pairs. This enables measurement of biphoton joint probability distributions with orders-of-magnitude greater dimensionality and faster acquisition speeds than traditional raster scanning of point detectors; to date, however, there has been no general formula available to optimize data collection. Here we analyze the dependence of such measurements on count rate, detector noise properties, and threshold levels. We derive expressions for the biphoton joint probability distribution and its signal-to-noise ratio (SNR), valid beyond the low-count regime up to detector saturation. The analysis gives operating parameters for global optimum SNR that may be specified prior to measurement. We find excellent agreement with experimental measurements within the range of validity, and discuss discrepancies with the theoretical model for high thresholds. This work enables optimized measurement of the biphoton joint probability distribution in high-dimensional joint Hilbert spaces.


title: Realization problems for diffeomorphism groups

abstract: We discuss recent results and open questions on the broad theme of (Nielsen) realization problems. Beyond realizing subgroups of mapping class groups, there are many other natural instances where one can ask if a surjection from a group of diffeomorphisms of a manifold to another group admits a section over particular subgroups. This survey includes many open problems, and some short proofs of new results that are illustrative of key techniques; drawing attention to parallels between problems arising in different areas.


title: A New Registration Approach for Dynamic Analysis of Calcium Signals in Organs

abstract: Wing disc pouches of fruit flies are a powerful genetic model for studying physiological intercellular calcium ($Ca^{2+}$) signals for dynamic analysis of cell signaling in organ development and disease studies. A key to analyzing spatial-temporal patterns of $Ca^{2+}$ signal waves is to accurately align the pouches across image sequences. However, pouches in different image frames may exhibit extensive intensity oscillations due to $Ca^{2+}$ signaling dynamics, and commonly used multimodal non-rigid registration methods may fail to achieve satisfactory results. In this paper, we develop a new two-phase non-rigid registration approach to register pouches in image sequences. First, we conduct segmentation of the region of interest. (i.e., pouches) using a deep neural network model. Second, we obtain an optimal transformation and align pouches across the image sequences. Evaluated using both synthetic data and real pouch data, our method considerably outperforms the state-of-the-art non-rigid registration methods.


title: Nonmetricity formulation of general relativity and its scalar-tensor extension

abstract: Einstein's celebrated theory of gravitation can be presented in three forms: general relativity, teleparallel gravity, and the rarely considered before symmetric teleparallel gravity. Extending the latter, we introduce a new class of theories where a scalar field is coupled nonminimally to nonmetricity $Q$, which here encodes the gravitational effects like curvature $R$ in general relativity or torsion $T$ in teleparallel gravity. We point out the similarities and differences with analogous scalar-curvature and scalar-torsion theories by discussing the field equations, role of connection, conformal transformations, relation to $f(Q)$ theory, and cosmology. The equations for spatially flat universe coincide with those of teleparallel dark energy, thus allowing to explain accelerating expansion.


title: 2MASS J13243553+6358281 is an Early T-Type Planetary-Mass Object in the AB Doradus Moving Group

abstract: We present new radial velocity and trigonometric parallax measurements indicating that the unusually red and photometrically variable T2 dwarf 2MASS J13243553+6358281 is a member of the young (~150 Myr) AB Doradus moving group based on its space velocity. We estimate its model-dependent mass in the range 11-12 $M_{\rm Jup}$ at the age of the AB Doradus moving group, and its trigonometric parallax distance of 12.7 $\pm$ 1.5 pc makes it one of the nearest known isolated planetary-mass objects. The unusually red continuum of 2MASS J13243553+6358281 in the near-infrared was previously suspected to be caused by an unresolved L+T brown dwarf binary, although it was never observed with high-spatial resolution imaging. This new evidence of youth suggests that a low surface gravity may be sufficient to explain this peculiar feature. Using the new parallax we find that its absolute $J$-band magnitude is ~0.4 mag fainter than equivalent-type field brown dwarfs, suggesting that the binary hypothesis is unlikely. The fundamental properties of 2MASS J13243553+6358281 follow the spectral type sequence of other known high-likelihood members of the AB Doradus moving group. The effective temperature of 2MASS J13243553+6358281 provides the first precise constraint on the L/T transition at a known young age, and indicates that it happens at a temperature of ~1150 K at ~150 Myr, compared to ~1250 K for field brown dwarfs.


title: A Strange Vertex Condition Coming From Nowhere

abstract: We prove norm-resolvent and spectral convergence in $L^2$ of solutions to the Neumann Poisson problem $-\Delta u_\varepsilon = f$ on a domain $\Omega_\varepsilon$ perforated by Dirichlet-holes and shrinking to a 1-dimensional interval. The limit $u$ satisfies an equation of the type $-u''+\mu u = f$ on the interval $(0,1)$, where $\mu$ is a positive constant. As an application we study the convergence of solutions in perforated graph-like domains. We show that if the scaling between the edge neighbourhood and the vertex neighbourhood is chosen correctly, the constant $\mu$ will appear in the vertex condition of the limit problem. In particular, this implies that the spectrum of the resulting quantum graph is altered in a controlled way by the perforation.


title: Practical Bayesian Modeling and Inference for Massive Spatial Datasets On Modest Computing Environments

abstract: With continued advances in Geographic Information Systems and related computational technologies, statisticians are often required to analyze very large spatial datasets. This has generated substantial interest over the last decade, already too vast to be summarized here, in scalable methodologies for analyzing large spatial datasets. Scalable spatial process models have been found especially attractive due to their richness and flexibility and, particularly so in the Bayesian paradigm, due to their presence in hierarchical model settings. However, the vast majority of research articles present in this domain have been geared toward innovative theory or more complex model development. Very limited attention has been accorded to approaches for easily implementable scalable hierarchical models for the practicing scientist or spatial analyst. This article is submitted to the Practice section of the journal with the aim of developing massively scalable Bayesian approaches that can rapidly deliver Bayesian inference on spatial process that are practically indistinguishable from inference obtained using more expensive alternatives. A key emphasis is on implementation within very standard (modest) computing environments (e.g., a standard desktop or laptop) using easily available statistical software packages without requiring message-parsing interfaces or parallel programming paradigms. Key insights are offered regarding assumptions and approximations concerning practical efficiency.


title: Edu-Edition Spreadsheet Competency Framework

abstract: Based on the Spreadsheet Competency Framework for finance professionals, in the present paper we introduce the Edu-Edition of the Spreadsheet Competency Framework (E2SCF). We claim that building spreadsheet competences should start in education, as early as possible, and this process is a lot more effective if support arrives from expert teachers. The main feature of E2SCF is high mathability computer-supported real world problem solving. This approach is based on - from the very beginning of training - a two-directional knowledge transfer, data and error analysis and handling, and the programming aspect of spreadsheets. Based on these features, E2SCF is set up for basic and general users to build up firm spreadsheet knowledge and to develop transferable problem solving skills and competences.


title: Approximating power by weights

abstract: Determining the power distribution of the members of a shareholder meeting or a legislative committee is a well-known problem for many applications. In some cases it turns out that power is nearly proportional to relative voting weights, which is very beneficial for both theoretical considerations and practical computations with many members. We present quantitative approximation results with precise error bounds for several power indices as well as impossibility results for such approximations between power and weights.


title: Asymptotically safe f(R)-gravity coupled to matter I: the polynomial case

abstract: We use the functional renormalization group equation for the effective average action to study the non-Gaussian renormalization group fixed points (NGFPs) arising within the framework of f(R)-gravity minimally coupled to an arbitrary number of scalar, Dirac, and vector fields. Based on this setting we provide comprehensible estimates which gravity-matter systems give rise to NGFPs suitable for rendering the theory asymptotically safe. The analysis employs an exponential split of the metric fluctuations and retains a 7-parameter family of coarse-graining operators allowing the inclusion of non-trivial endomorphisms in the regularization procedure. For vanishing endomorphisms, it is established that gravity coupled to the matter content of the standard model of particle physics and many beyond the standard model extensions exhibit NGFPs whose properties are strikingly similar to the case of pure gravity: there are two UV-relevant directions and the position and critical exponents converge rapidly when higher powers of the scalar curvature are included. Conversely, none of the phenomenologically interesting gravity-matter systems exhibits a stable NGFP when a Type II coarse graining operator is employed. Our analysis resolves this tension by demonstrating that the NGFPs seen in the two settings belong to different universality classes.


title: High-Speed Photometry of Gaia14aae: An Eclipsing AM CVn That Challenges Formation Models

abstract: AM CVn-type systems are ultra-compact, hydrogen-deficient accreting binaries with degenerate or semi-degenerate donors. The evolutionary history of these systems can be explored by constraining the properties of their donor stars. We present high-speed photometry of Gaia14aae, an AM CVn with a binary period of 49.7 minutes and the first AM CVn in which the central white dwarf is fully eclipsed by the donor star. Modelling of the lightcurves of this system allows for the most precise measurement to date of the donor mass of an AM CVn, and relies only on geometric and well-tested physical assumptions. We find a mass ratio $q = M_2/M_1 = 0.0287 \pm 0.0020$ and masses $M_1 = 0.87 \pm 0.02 M_\odot$ and $M_2 = 0.0250 \pm 0.0013 M_\odot$. We compare these properties to the three proposed channels for AM CVn formation. Our measured donor mass and radius do not fit with the contraction that is predicted for AM CVn donors descended from white dwarfs or helium stars at long orbital periods. The donor properties we measure fall in a region of parameter space in which systems evolved from hydrogen-dominated cataclysmic variables are expected, but such systems should show spectroscopic hydrogen, which is not seen in Gaia14aae. The evolutionary history of this system is therefore not clear. We consider a helium-burning star or an evolved cataclysmic variable to be the most likely progenitors, but both models require additional processes and/or fine-tuning to fit the data. Additionally, we calculate an updated ephemeris which corrects for an anomalous time measurement in the previously published ephemeris.


title: Goal-Oriented Chatbot Dialog Management Bootstrapping with Transfer Learning

abstract: Goal-Oriented (GO) Dialogue Systems, colloquially known as goal oriented chatbots, help users achieve a predefined goal (e.g. book a movie ticket) within a closed domain. A first step is to understand the user's goal by using natural language understanding techniques. Once the goal is known, the bot must manage a dialogue to achieve that goal, which is conducted with respect to a learnt policy. The success of the dialogue system depends on the quality of the policy, which is in turn reliant on the availability of high-quality training data for the policy learning method, for instance Deep Reinforcement Learning. Due to the domain specificity, the amount of available data is typically too low to allow the training of good dialogue policies. In this paper we introduce a transfer learning method to mitigate the effects of the low in-domain data availability. Our transfer learning based approach improves the bot's success rate by 20% in relative terms for distant domains and we more than double it for close domains, compared to the model without transfer learning. Moreover, the transfer learning chatbots learn the policy up to 5 to 10 times faster. Finally, as the transfer learning approach is complementary to additional processing such as warm-starting, we show that their joint application gives the best outcomes.


title: Evolutionary branching via replicator-mutator equations

abstract: We consider a class of non-local reaction-diffusion problems, referred to as replicator-mutator equations in evolutionary genetics. For a confining fitness function, we prove well-posedness and write the solution explicitly, via some underlying Schr\"odinger spectral elements (for which we provide new and non-standard estimates). As a consequence, the long time behaviour is determined by the principal eigenfunction or ground state. Based on this, we discuss (rigorously and via numerical explorations) the conditions on the fitness function and the mutation rate for evolutionary branching to occur.


title: Unfair and Anomalous Evolutionary Dynamics from Fluctuating Payoffs

abstract: Evolution occurs in populations of reproducing individuals. Reproduction depends on the payoff a strategy receives. The payoff depends on the environment that may change over time, on intrinsic uncertainties, and on other sources of randomness. These temporal variations in the payoffs can affect which traits evolve. Understanding evolutionary game dynamics that are affected by varying payoffs remains difficult. Here we study the impact of arbitrary amplitudes and covariances of temporally varying payoffs on the dynamics. The evolutionary dynamics may be "unfair", meaning that, on average, two coexisting strategies may persistently receive different payoffs. This mechanism can induce an anomalous coexistence of cooperators and defectors in the Prisoner's Dilemma, and an unexpected selection reversal in the Hawk-Dove game.


title: The bipolar jet of the symbiotic star R Aquarii: A study of its morphology using the high-resolution HST WFC3/UVIS camera

abstract: R Aqr is a symbiotic binary system consisting of a Mira variable with a pulsation period of 387 days and a hot companion which is presumably a white dwarf with an accretion disk. This binary system is the source of a prominent bipolar gaseous outflow. We use high spatial resolution and sensitive images from the Hubble Space Telescope to identify and investigate the different structural components that form the complex morphology of the R Aqr jet . Methods. We present new high-resolution HST WFC3/UVIS narrow-band images of the R Aqr jet obtained in 2013/14 in the light of the [OIII] 5007, [OI] 6300, [NII] 6583, and Ha emission lines. These images also allow us to produce detailed maps of the jet flow in several line ratios such as [OIII]/[OI] and [NII]/[OI] which are sensitive to the outflow temperature and its hydrogen ionisation fraction. The new emission maps together with archival HST data are used to derive and analyse the proper motion of prominent emitting features which can be traced over 20 years with the HST observations. The images reveal the fine gas structure of the jet out to distances of a few ten arcseconds from the central region, as well as in the innermost region, within a few arcseconds around the stellar source. They reveal for the first time the straight highly-collimated jet component which can be traced to up to 900 AU in the NE direction. Images in [OIII], [OI], and [NII] clearly show a helical pattern in the jet beams which may derive from the small-scale precession of the jet. The highly-collimated jet is accompanied by a wide opening angle outflow which is filled by low excitation gas. The position angles of the jet structures as well as their opening angles are calculated. Our measurements of the proper motions of some prominent emission knots confirm the scenario of gas acceleration during the propagation of the outflow.


title: Three dimensional spectrometer

abstract: We present a novel design of 3D spectrometer that can retrieve 3D spectral profile in a single measurement. The 3D spectrometer design is built upon the concept of compressed sensing to make it possible to retrieve 3D information from 2D data from a screen/camera. In contrast to common spectrometers, the 3D spectrometer uses a wide slit instead of a narrow slit and retrieve the 3D datacube that consists of 2D spatial and 1D spectral information. Numerical tests were performed to simulate the retrieval of spectral profiles. The results show that the retrieved profiles match well the original profiles. It is also shown that the retrieved signal from the 3D spectrometer is robust enough for a further post-processing analysis.


title: The spectrum of non-centrosymmetrically layered spherical cavity resonator. I.The mode decomposition method

abstract: We develop a theoretical method for solving Maxwell's equations to obtain the frequency spectra of inhomogeneous and asymmetric cavity resonators using a couple of effective Debye-type potentials. The structure we study specifically is the layered spherical cavity resonator with symmetrically or asymmetrically inserted inner dielectric sphere. The comparison of the exact numerical results obtained for the frequency spectrum of layered cavity resonator with centrosymmetrically inserted sphere and the spectrum found from the suggested theory reveals good agreement at the initial part of the frequency axis. The coincidence accuracy depends on the number of trial resonant modes that we use while approving our method numerically.


title: Dual terahertz comb spectroscopy with a single free-running fibre laser

abstract: Dual THz comb spectroscopy has the potential to be used as universal THz spectroscopy with high spectral resolution, high spectral accuracy, and broad spectral coverage; however, the requirement for dual stabilized femtosecond lasers hampers its versatility due to the bulky size, high complexity, and high cost. We here report the first demonstration of dual THz comb spectroscopy using a single free-running fibre laser. By tuning the cavity-loss-dependent gain profile with an intracavity Lyot filter together with precise management of the cavity length and dispersion, dual-wavelength pulsed light beams with slightly detuned repetition frequencies are generated in a single laser cavity. Due to sharing of the same cavity, such pulsed light beams suffer from common-mode fluctuation of the repetition frequency, and hence the corresponding frequency difference between them is passively stable around a few hundred hertz within millihertz fluctuation. This considerably stable frequency difference enables dual THz comb spectroscopy with a single free-running fibre laser. While greatly reducing the size, complexity, and cost of the laser source by use of a single free-running fibre laser, the dual THz comb spectroscopy system maintains a spectral bandwidth and dynamic range of spectral power comparable to a system equipped with dual stabilized fibre lasers, and can be effectively applied to high-precision spectroscopy of acetonitrile gas at atmospheric pressure. The demonstrated results indicate that this system is an attractive solution for practical applications of not only THz spectroscopy but also THz-pulse-based measurements.


title: The effects of anger on automated long-term-spectra based speaker-identification

abstract: Forensic speaker identification has traditionally considered approaches based on long term spectra analysis as especially robust, given that they work well for short recordings, are not sensitive to changes in the intensity of the sample, and continue to function in the presence of noise and limited passband. We find, however, that anger induces a significant distortion of the acoustic signal for long term spectra analysis purposes. Even moderate anger offsets speaker identification results by 33% in the direction of a different speaker altogether. Thus, caution should be exercised when applying this tool.


title: Snort Intrusion Detection System with Intel Software Guard Extension (Intel SGX)

abstract: Network Function Virtualization (NFV) promises the benefits of reduced infrastructure, personnel, and management costs by outsourcing network middleboxes to the public or private cloud. Unfortunately, running network functions in the cloud entails security challenges, especially for complex stateful services. In this paper, we describe our experiences with hardening the king of middleboxes - Intrusion Detection Systems (IDS) - using Intel Software Guard Extensions (Intel SGX) technology. Our IDS secured using Intel SGX, called SEC-IDS, is an unmodified Snort 3 with a DPDK network layer that achieves 10Gbps line rate. SEC-IDS guarantees computational integrity by running all Snort code inside an Intel SGX enclave. At the same time, SEC-IDS achieves near-native performance, with throughput close to 100 percent of vanilla Snort 3, by retaining network I/O outside of the enclave. Our experiments indicate that performance is only constrained by the modest Enclave Page Cache size available on current Intel SGX Skylake based E3 Xeon platforms. Finally, we kept the porting effort minimal by using the Graphene-SGX library OS. Only 27 Lines of Code (LoC) were modified in Snort and 178 LoC in Graphene-SGX itself.


title: Learning Semantic Segmentation with Diverse Supervision

abstract: Models based on deep convolutional neural networks (CNN) have significantly improved the performance of semantic segmentation. However, learning these models requires a large amount of training images with pixel-level labels, which are very costly and time-consuming to collect. In this paper, we propose a method for learning CNN-based semantic segmentation models from images with several types of annotations that are available for various computer vision tasks, including image-level labels for classification, box-level labels for object detection and pixel-level labels for semantic segmentation. The proposed method is flexible and can be used together with any existing CNN-based semantic segmentation networks. Experimental evaluation on the challenging PASCAL VOC 2012 and SIFT-flow benchmarks demonstrate that the proposed method can effectively make use of diverse training data to improve the performance of the learned models.


title: Adaptive Memory Networks

abstract: We present Adaptive Memory Networks (AMN) that processes input-question pairs to dynamically construct a network architecture optimized for lower inference times for Question Answering (QA) tasks. AMN processes the input story to extract entities and stores them in memory banks. Starting from a single bank, as the number of input entities increases, AMN learns to create new banks as the entropy in a single bank becomes too high. Hence, after processing an input-question(s) pair, the resulting network represents a hierarchical structure where entities are stored in different banks, distanced by question relevance. At inference, one or few banks are used, creating a tradeoff between accuracy and performance. AMN is enabled by dynamic networks that allow input dependent network creation and efficiency in dynamic mini-batching as well as our novel bank controller that allows learning discrete decision making with high accuracy. In our results, we demonstrate that AMN learns to create variable depth networks depending on task complexity and reduces inference times for QA tasks.


title: Quantum energy exchange and refrigeration: A full-counting statistics approach

abstract: We formulate a full-counting statistics description to study energy exchange in multi-terminal junctions. Our approach applies to quantum systems that are coupled either additively or non-additively (cooperatively) to multiple reservoirs. We derive a Markovian Redfield-type equation for the counting-field dependent reduced density operator. Under the secular approximation, we confirm that the cumulant generating function satisfies the heat exchange fluctuation theorem. Our treatment thus respects the second law of thermodynamics. We exemplify our formalism on a multi-terminal two-level quantum system, and apply it to realize the smallest quantum absorption refrigerator, operating through engineered reservoirs, and achievable only through a cooperative bath interaction model.


title: Bulk-boundary quantum oscillations in inhomogeneous Weyl semimetals

abstract: Weyl fermions in an external magnetic field exhibit the chiral anomaly, a non-conservation of chiral fermions. In a Weyl semimetal, a spatially inhomogeneous Weyl node separation causes similar effect by creating an intrinsic pseudo-magnetic field with an opposite sign for nodes of opposite chirality. In the present work we study the interplay of external and intrinsic fields. In particular, we focus on quantum oscillations due to bulk-boundary trajectories. When caused by an external field, such oscillations are a proven experimental technique to detect Weyl semimetals. We show that the intrinsic field leaves hallmarks on such oscillations by decreasing the period of the oscillations in an analytically traceable manner. The oscillations can thus be used to test the effect of an intrinsic field and to extract its strength.


title: Note: Distance-Based Network Partitioning

abstract: A new method for identifying soft communities in networks is proposed. Reference nodes, either selected using a priori information about the network or according to relevant node measurements, are obtained. Distance vectors between each network node and the reference nodes are then used for defining a multidimensional coordinate system representing the community structure of the network at many different scales. For modular networks, the distribution of nodes in this space often results in a well-separated clustered structure, with each cluster corresponding to a community. The potential of the method is illustrated with respect to a spatial network model and the Zachary's karate club network.


title: Accuracy of Shack-Hartmann wavefront sensor using a coherent wound fibre image bundle

abstract: Shack-Hartmann wavefront sensors using wound fibre bundles are desired for multi-object adaptive optical systems to provide large multiplex positioned by Starbugs. The use of the large-sized wound fibre bundle provides the exibility to use more sub-apertures wavefront sensor for ELTs. These compact wavefront sensors take advantage of large focal surfaces such as the Giant Magellan Telescope. The focus of this paper is to study the wound fibre image bundle structure defects effect on the centroid measurement accuracy of a Shack-Hartmann wavefront sensor. We use the first moment centroid method to estimate the centroid of a focused Gaussian beam sampled by a simulated bundle. Spot estimation accuracy with wound fibre image bundle and its structure impact on wavefront measurement accuracy statistics are addressed. Our results show that when the measurement signal to noise ratio is high, the centroid measurement accuracy is dominated by the wound fibre image bundle structure, e.g. tile angle and gap spacing. For the measurement with low signal to noise ratio, its accuracy is influenced by the read noise of the detector instead of the wound fibre image bundle structure defects. We demonstrate this both with simulation and experimentally. We provide a statistical model of the centroid and wavefront error of a wound fibre image bundle found through experiment.


title: Dimension Reduction via Gaussian Ridge Functions

abstract: Ridge functions have recently emerged as a powerful set of ideas for subspace-based dimension reduction. In this paper we begin by drawing parallels between ridge subspaces, sufficient dimension reduction and active subspaces, contrasting between techniques rooted in statistical regression and those rooted in approximation theory. This sets the stage for our new algorithm that approximates what we call a Gaussian ridge function---the posterior mean of a Gaussian process on a dimension-reducing subspace---suitable for both regression and approximation problems. To compute this subspace we develop an iterative algorithm that optimizes over the Stiefel manifold to compute the subspace, followed by an optimization of the hyperparameters of the Gaussian process. We demonstrate the utility of the algorithm on two analytical functions, where we obtain near exact ridge recovery, and a turbomachinery case study, where we compare the efficacy of our approach with three well-known sufficient dimension reduction methods: SIR, SAVE, CR. The comparisons motivate the use of the posterior variance as a heuristic for identifying the suitability of a dimension-reducing subspace.


title: Existence of solitary waves in one dimensional peridynamics

abstract: We give a rigorous proof of existence for solitary waves of a peridynamics model in one space dimension recently investigated by Silling (J. Mech. Phys. Solids 96:121--132, 2016). We adapt the variational framework developed by Friesecke and Wattis (Comm. Math Phys. 161:391--418, 1994) for the Fermi-Pasta-Ulam-Tsingou lattice equations to treat a truncated problem which cuts off short-range interactions, then pass to the limit.


title: Zero-adjusted Birnbaum-Saunders regression model

abstract: In this paper we introduce the zero-adjusted Birnbaum-Saunders regression model. This new model generalizes at least seven Birnbaum-Saunders regression models. The idea of this modeling is mixing a degenerate distribution at zero with a Birnbaum-Saunders distribution. Besides the capacity to account for excess zeros, the zero-adjusted Birnbaum-Saunders distribution additionally produces an attractive modeling structure to right-skewed data. In this model, the mean and precision parameter of the Birnbaum-Saunders distribution and the probability of zeros can be related to linear and/or non-linear predictors through link functions. We derive a type of residual to perform diagnostic analysis and a perturbation scheme for identifying those observations that exert unusual influence on the estimation process. Finally, two applications to real data show the potential of the model.


title: Analysis of Fast Alternating Minimization for Structured Dictionary Learning

abstract: Methods exploiting sparsity have been popular in imaging and signal processing applications including compression, denoising, and imaging inverse problems. Data-driven approaches such as dictionary learning and transform learning enable one to discover complex image features from datasets and provide promising performance over analytical models. Alternating minimization algorithms have been particularly popular in dictionary or transform learning. In this work, we study the properties of alternating minimization for structured (unitary) sparsifying operator learning. While the algorithm converges to the stationary points of the non-convex problem in general, we prove rapid local linear convergence to the underlying generative model under mild assumptions. Our experiments show that the unitary operator learning algorithm is robust to initialization.


title: Numerical solution of variable order fractional differential equations

abstract: A method for the numerical solution of variable order (VO) fractional differential equations (FDE) is presented. The method applies to linear as well as to nonlinear VO-FDEs. The Caputo type VO fractional derivative is employed. First, an simple expression, which approximates the VO fractional derivative, is established and then a procedure based on this approximation is developed to solve VO-FDEs linear and nonlinear, both explicit and implicit. VO-FDEs with variable coefficients are also treated. The method is illustrated by solving the second order VO-FDE describing the response of the VO fractional oscillator, linear and nonlinear (Duffing). However, it can be straightforwardly extended to higher order VO-FDEs. The presented method, in addition to its effectiveness, is simple to implement and program on a computer. The obtained results validate the efficiency and accuracy of the developed method


title: Real-world Multi-object, Multi-grasp Detection

abstract: A deep learning architecture is proposed to predict graspable locations for robotic manipulation. It considers situations where no, one, or multiple object(s) are seen. By defining the learning problem to be classification with null hypothesis competition instead of regression, the deep neural network with RGB-D image input predicts multiple grasp candidates for a single object or multiple objects, in a single shot. The method outperforms state-of-the-art approaches on the Cornell dataset with 96.0% and 96.1% accuracy on image-wise and object- wise splits, respectively. Evaluation on a multi-object dataset illustrates the generalization capability of the architecture. Grasping experiments achieve 96.0% grasp localization and 88.0% grasping success rates on a test set of household objects. The real-time process takes less than .25 s from image to plan.


title: Multipath Communication with Finite Sliding Window Network Coding for Ultra-Reliability and Low Latency

abstract: We use random linear network coding (RLNC) based scheme for multipath communication in the presence of lossy links with different delay characteristics to obtain ultra-reliability and low latency. A sliding window version of RLNC is proposed where the coded packets are generated using packets in a window size and are inserted among systematic packets in different paths. The packets are scheduled in the paths in a round robin fashion proportional to the data rates. We use finite encoding and decoding window size and do not rely on feedback for closing the sliding window, unlike the previous work. Our implementation of two paths with LTE and WiFi characteristics shows that the proposed sliding window scheme achieves better latency compared to the block RLNC code. It is also shown that the proposed scheme achieves low latency communication through multiple paths compared to the individual paths for bursty traffic by translating the throughput on both the paths into latency gain.


title: Nanomechanical characterization of the Kondo charge dynamics in a carbon nanotube

abstract: Using the transversal vibration resonance of a suspended carbon nanotube as charge detector for its embedded quantum dot, we investigate the case of strong Kondo correlations between a quantum dot and its leads. We demonstrate that even when large Kondo conductance is carried at odd electron number, the charging behaviour remains similar between odd and even quantum dot occupation. While the Kondo conductance is caused by higher order processes, a sequential tunneling only model can describe the time-averaged charge. The gate potentials of maximum current and fastest charge increase display a characteristic relative shift, which is suppressed at increased temperature. These observations agree very well with models for Kondo-correlated quantum dots.


title: The Satisfiability of Extended Word Equations: The Boundary Between Decidability and Undecidability

abstract: The study of word equations (or the existential theory of equations over free monoids) is a central topic in mathematics and theoretical computer science. The problem of deciding whether a given word equation has a solution was shown to be decidable by Makanin in the late 1970s, and since then considerable work has been done on this topic. In recent years, this decidability question has gained critical importance in the context of string SMT solvers for security analysis. Further, many extensions (e.g., quantifier-free word equations with linear arithmetic over the length function) and fragments (e.g., restrictions on the number of variables) of this theory are important from a theoretical point of view, as well as for program analysis applications. Motivated by these considerations, we prove several new results and thus shed light on the boundary between decidability and undecidability for many fragments and extensions of the first order theory of word equations.


title: Simulations of relativistic-quantum plasmas using real-time lattice scalar QED

abstract: Real-time lattice quantum electrodynamics (QED) provides a unique tool for simulating plasmas in the strong-field regime, where collective plasma scales are not well-separated from relativistic-quantum scales. As a toy model, we study scalar QED, which describes self-consistent interactions between charged bosons and electromagnetic fields. To solve this model on a computer, we first discretize the scalar-QED action on a lattice, in a way that respects geometric structures of exterior calculus and U(1)-gauge symmetry. The lattice scalar QED can then be solved, in the classical-statistics regime, by advancing an ensemble of statistically equivalent initial conditions in time, using classical field equations obtained by extremizing the discrete action. To demonstrate the capability of our numerical scheme, we apply it to two example problems. The first example is the propagation of linear waves, where we recover analytic wave dispersion relations using numerical spectrum. The second example is an intense laser interacting with a 1D plasma slab, where we demonstrate natural transition from wakefield acceleration to pair production when the wave amplitude exceeds the Schwinger threshold. Our real-time lattice scheme is fully explicit and respects local conservation laws, making it reliable for long-time dynamics. The algorithm is readily parallelized using domain decomposition, and the ensemble may be computed using quantum parallelism in the future.


title: Diophantine approximation on the parabola with non-monotonic approximation functions

abstract: We show that the parabola is of strong Khintchine type for convergence, which is the first result of its kind for curves. Moreover, Jarnik type theorems are established in both the simultaneous and the dual settings, without monotonicity on the approximation function. To achieve the above, we prove a new counting result for the number of rational points with fixed denominators lying close to the parabola, which uses Burgess's bound on short character sums.


title: Toward Optimal Coupon Allocation in Social Networks: An Approximate Submodular Optimization Approach

abstract: CMO Council reports that 71\% of internet users in the U.S. were influenced by coupons and discounts when making their purchase decisions. It has also been shown that offering coupons to a small fraction of users (called seed users) may affect the purchase decisions of many other users in a social network. This motivates us to study the optimal coupon allocation problem, and our objective is to allocate coupons to a set of users so as to maximize the expected cascade. Different from existing studies on influence maximizaton (IM), our framework allows a general utility function and a more complex set of constraints. In particular, we formulate our problem as an approximate submodular maximization problem subject to matroid and knapsack constraints. Existing techniques relying on the submodularity of the utility function, such as greedy algorithm, can not work directly on a non-submodular function. We use $\epsilon$ to measure the difference between our function and its closest submodular function and propose a novel approximate algorithm with approximation ratio $\beta(\epsilon)$ with $\lim_{\epsilon\rightarrow 0}\beta(\epsilon)=1-1/e$. This is the best approximation guarantee for approximate submodular maximization subject to a partition matroid and knapsack constraints, our results apply to a broad range of optimization problems that can be formulated as an approximate submodular maximization problem.


title: Predicting outcomes for games of skill by redefining what it means to win

abstract: The Elo rating system is a highly successful ranking algorithm for games of skill where, by construction, one team wins and the other loses. A primary limitation of the original Elo algorithm is its inability to predict information beyond a match's win-loss probability. Specifically, the victor is awarded the same point bounty if he beats a team by 1 point or 10 points; only the rating difference between the team and its opponent affects the match bounty. In this work, we explain that Elo ratings and predictions can be naturally extended to include margin-of-victory information by simply redefining "what it means to win." We create ratings for each value of the margin-of-victory and use these ratings to predict the full distribution of point spread outcomes for matches which have not yet been played.


title: Implicative algebras: a new foundation for realizability and forcing

abstract: We introduce the notion of implicative algebra, a simple algebraic structure intended to factorize the model constructions underlying forcing and realizability (both in intuitionistic and classical logic). The salient feature of this structure is that its elements can be seen both as truth values and as (generalized) realizers, thus blurring the frontier between proofs and types. We show that each implicative algebra induces a (Set-based) tripos, using a construction that is reminiscent from the construction of a realizability tripos from a partial combinatory algebra. Relating this construction with the corresponding constructions in forcing and realizability, we conclude that the class of implicative triposes encompass all forcing triposes (both intuitionistic and classical), all classical realizability triposes (in the sense of Krivine) and all intuitionistic realizability triposes built from partial combinatory algebras.


title: Microwave-assisted Rydberg Electromagnetically induced transparency

abstract: We demonstrate electromagnetically induced transparency (EIT) in a four-level cascade-like system, where the two upper levels are Rydberg states coupled by a microwave field. A two-photon transition consisting of an off-resonant microwave field and an off-resonant optical field forms an effective coupling field to induce transparency of the probe light. We characterize the Rabi frequency of the effective coupling field, as well as the EIT microwave spectra. The results show that microwave assisted EIT allows us to efficiently access Rydberg states with relatively high orbital angular momentum $\ell=3$, which is promising for the study of exotic Rydberg molecular states.


title: Scalable L\'evy Process Priors for Spectral Kernel Learning

abstract: Gaussian processes are rich distributions over functions, with generalization properties determined by a kernel function. When used for long-range extrapolation, predictions are particularly sensitive to the choice of kernel parameters. It is therefore critical to account for kernel uncertainty in our predictive distributions. We propose a distribution over kernels formed by modelling a spectral mixture density with a L\'evy process. The resulting distribution has support for all stationary covariances--including the popular RBF, periodic, and Mat\'ern kernels--combined with inductive biases which enable automatic and data efficient learning, long-range extrapolation, and state of the art predictive performance. The proposed model also presents an approach to spectral regularization, as the L\'evy process introduces a sparsity-inducing prior over mixture components, allowing automatic selection over model order and pruning of extraneous components. We exploit the algebraic structure of the proposed process for $\mathcal{O}(n)$ training and $\mathcal{O}(1)$ predictions. We perform extrapolations having reasonable uncertainty estimates on several benchmarks, show that the proposed model can recover flexible ground truth covariances and that it is robust to errors in initialization.


title: A generalization of Menon's identity with Dirichlet characters

abstract: The classical Menon's identity [7] states that \begin{equation*}\label{oldbegin1} \sum_{\substack{a\in\Bbb Z_n^\ast }}\gcd(a -1,n)=\varphi(n) \sigma_{0} (n), \end{equation*} where for a positive integer $n$, $\Bbb Z_n^\ast$ is the group of units of the ring $\Bbb Z_n=\Bbb Z/n\Bbb Z$, $\gcd(\ ,\ )$ represents the greatest common divisor, $\varphi(n)$ is the Euler's totient function and $\sigma_{k} (n) =\sum_{d|n } d^{k}$ is the divisor function. In this paper, we generalize Menon's identity with Dirichlet characters in the following way: \begin{equation*} \sum_{\substack{a\in\Bbb Z_n^\ast b_1, ..., b_k\in\Bbb Z_n}} \gcd(a-1,b_1, ..., b_k, n)\chi(a)=\varphi(n)\sigma_k\left(\frac{n}{d}\right), \end{equation*} where $k$ is a non-negative integer and $\chi$ is a Dirichlet character modulo $n$ whose conductor is $d$. Our result can be viewed as an extension of Zhao and Cao's result [16] to $k>0$. It can also be viewed as an extension of Sury's result [12] to Dirichlet characters.


title: Stability for Representations of Hecke Algebras of type $A$

abstract: In this paper we introduce the notion of the stability of a sequence of modules over Hecke algebras. We prove that a finitely generated consistent sequence associated with Hecke algebras is representation stable.


title: Persistent Homology and the Upper Box Dimension

abstract: We introduce a fractal dimension for a metric space defined in terms of the persistent homology of extremal subsets of that space. We exhibit hypotheses under which this dimension is comparable to the upper box dimension; in particular, the dimensions coincide for subsets of $\mathbb{R}^2$ whose upper box dimension exceeds $1.5.$ These results are related to extremal questions about the number of persistent homology intervals of a set of $n$ points in a metric space.


title: How light absorption modifies the radiative force on a microparticle in optical tweezers

abstract: Reflection and refraction of light can be used to trap small dielectric particles in the geometrical optics regime. Absorption of light is usually neglected in theoretical calculations, but it is known that it occurs in the optical trapping of metallic particles. Also, recent experiments with semi-transparent microparticles have shown that absorption of light is important to explain their optically induced oscillations. Here, we propose a generalization of Ashkin's model for the radiative force exerted on a spherical bead, including the contribution due to attenuation/absorption of light in the bulk of the particle. We discuss in detail the balance between refraction, reflection and absorption for different optical parameters and particle sizes. A detailed example is provided in order to clarify how the model can be applied, and it is obtained that the radiative force can either increase or decrease with absorption, depending on the particle size. Our findings contribute to the understanding of optical trapping of light-absorbing particles, and may be used to predict whenever absorption is important in real experiments.

categories: cond-mat.mes-hall physics.app-ph physics.bio-ph physics.optics


title: Scalable Preprocessing of High Volume Bird Acoustic Data

abstract: In this work, we examine the problem of efficiently preprocessing high volume bird acoustic data. We combine several existing preprocessing steps including noise reduction approaches into a single efficient pipeline by examining each process individually. We then utilise a distributed computing architecture to improve execution time. Using a master-slave model with data parallelisation, we developed a near-linear automated scalable system, capable of preprocessing bird acoustic recordings 21.76 times faster with 32 cores over 8 virtual machines, compared to a serial process. This work contributes to the research area of bioacoustic analysis, which is currently very active because of its potential to monitor animals quickly at low cost. Overcoming noise interference is a significant challenge in many bioacoustic studies, and the volume of data in these studies is increasing. Our work makes large scale bird acoustic analyses more feasible by parallelising important bird acoustic processing tasks to significantly reduce execution times.


title: A Kernel Based High Order "Explicit" Unconditionally Stable Scheme for Time Dependent Hamilton-Jacobi Equations

abstract: In this paper, a class of high order numerical schemes is proposed for solving Hamilton-Jacobi (H-J) equations. This work is regarded as an extension of our previous work for nonlinear degenerate parabolic equations, see Christlieb et al. \emph{arXiv preprint arXiv:1707.09294},, which relies on a special kernel-based formulation of the solutions and successive convolution. When applied to the H-J equations, the newly proposed scheme attains genuinely high order accuracy in both space and time, and more importantly, it is unconditionally stable, hence allowing for much larger time step evolution compared with other explicit schemes and saving computational cost. A high order weighted essentially non-oscillatory methodology and a novel nonlinear filter are further incorporated to capture the correct viscosity solution. Furthermore, by coupling the recently proposed inverse Lax-Wendroff boundary treatment technique, this method is very flexible in handing complex geometry as well as general boundary conditions. We perform numerical experiments on a collection of numerical examples, including H-J equations with linear, nonlinear, convex or non-convex Hamiltonians. The efficacy and efficiency of the proposed scheme in approximating the viscosity solution of general H-J equations is verified.


title: The UV Emission of Stars in LAMOST Survey I. Catalogs

abstract: We present the ultraviolet magnitudes for over three million stars in the LAMOST survey, in which 2,202,116 stars are detected by $GALEX$. For 889,235 undetected stars, we develop a method to estimate their upper limit magnitudes. The distribution of (FUV $-$ NUV) shows that the color declines with increasing effective temperature for stars hotter than 7000 K in our sample, while the trend disappears for the cooler stars due to upper atmosphere emission from the regions higher than their photospheres. For stars with valid stellar parameters, we calculate the UV excesses with synthetic model spectra, and find that the (FUV $-$ NUV) vs. $R'_{\mathrm{FUV}}$ can be fitted with a linear relation and late-type dwarfs tend to have high UV excesses. There are 87,178 and 1,498,103 stars detected more than once in the visit exposures of $GALEX$ in the FUV and NUV, respectively. We make use of the quantified photometric errors to determine statistical properties of the UV variation, including intrinsic variability and the structure function on the timescale of days. The overall occurrence of possible false positives is below 1.3\% in our sample. UV absolute magnitudes are calculated for stars with valid parallaxes, which could serve as a possible reference frame in the NUV. We conclude that the colors related to UV provide good criteria to distinguish between M giants and M dwarfs, and the variability of RR Lyrae stars in our sample is stronger than that of other A and F stars.


title: Decentralized Control of Stochastically Switched Linear System with Unreliable Communication

abstract: We consider a networked control system (NCS) consisting of two plants, a global plant and a local plant, and two controllers, a global controller and a local controller. The global (resp. local) plant follows discrete-time stochastically switched linear dynamics with a continuous global (resp. local) state and a discrete global (resp. local) mode. We assume that the state and mode of the global plant are observed by both controllers while the state and mode of the local plant are only observed by the local controller. The local controller can inform the global controller of the local plant's state and mode through an unreliable TCP-like communication channel where successful transmissions are acknowledged. The objective of the controllers is to cooperatively minimize a modes-dependent quadratic cost over a finite time horizon. Following the method developed in [1] and [2], we construct a dynamic program based on common information and a decomposition of strategies, and use it to obtain explicit optimal strategies for the controllers. In the optimal strategies, both controllers compute a common estimate of the local plant's state. The global controller's action is linear in the state of the global plant and the common estimated state, and the local controller's action is linear in the actual states of both plants and the common estimated state. Furthermore, the gain matrices for the global controller depend on the global mode and its observation about the local mode, while the gain matrices for the local controller depend on the actual modes of both plants and the global controller's observation about the local mode.


title: Complex Network Classification with Convolutional Neural Network

abstract: Classifying large scale networks into several categories and distinguishing them according to their fine structures is of great importance with several applications in real life. However, most studies of complex networks focus on properties of a single network but seldom on classification, clustering, and comparison between different networks, in which the network is treated as a whole. Due to the non-Euclidean properties of the data, conventional methods can hardly be applied on networks directly. In this paper, we propose a novel framework of complex network classifier (CNC) by integrating network embedding and convolutional neural network to tackle the problem of network classification. By training the classifiers on synthetic complex network data and real international trade network data, we show CNC can not only classify networks in a high accuracy and robustness, it can also extract the features of the networks automatically.


title: Dense Instanton-Dyon Liquid Model: Diagrammatics

abstract: We revisit the instanton-dyon liquid model in the confined phase by using a non-linear Debye-Huckel (DH) resummation for the Coulomb interactions induced by the moduli, followed by a cluster expansion. The organization is shown to rapidly converge and yields center symmetry at high density. The dependence of these results on a finite vacuum angle are also discussed. We also formulate the hypernetted chain (HCN) resummation for the dense instanton-dyon liquid and use it to estimate the liquid pair correlation functions in the DH limit. At very low temperature, the dense limit interpolates between chains and rings of instanton-anti-instanton-dyons and a bcc crystal, with strong topological and magnetic correlations.


title: Causal Learning and Explanation of Deep Neural Networks via Autoencoded Activations

abstract: Deep neural networks are complex and opaque. As they enter application in a variety of important and safety critical domains, users seek methods to explain their output predictions. We develop an approach to explaining deep neural networks by constructing causal models on salient concepts contained in a CNN. We develop methods to extract salient concepts throughout a target network by using autoencoders trained to extract human-understandable representations of network activations. We then build a bayesian causal model using these extracted concepts as variables in order to explain image classification. Finally, we use this causal model to identify and visualize features with significant causal influence on final classification.


title: ExpNet: Landmark-Free, Deep, 3D Facial Expressions

abstract: We describe a deep learning based method for estimating 3D facial expression coefficients. Unlike previous work, our process does not relay on facial landmark detection methods as a proxy step. Recent methods have shown that a CNN can be trained to regress accurate and discriminative 3D morphable model (3DMM) representations, directly from image intensities. By foregoing facial landmark detection, these methods were able to estimate shapes for occluded faces appearing in unprecedented in-the-wild viewing conditions. We build on those methods by showing that facial expressions can also be estimated by a robust, deep, landmark-free approach. Our ExpNet CNN is applied directly to the intensities of a face image and regresses a 29D vector of 3D expression coefficients. We propose a unique method for collecting data to train this network, leveraging on the robustness of deep networks to training label noise. We further offer a novel means of evaluating the accuracy of estimated expression coefficients: by measuring how well they capture facial emotions on the CK+ and EmotiW-17 emotion recognition benchmarks. We show that our ExpNet produces expression coefficients which better discriminate between facial emotions than those obtained using state of the art, facial landmark detection techniques. Moreover, this advantage grows as image scales drop, demonstrating that our ExpNet is more robust to scale changes than landmark detection methods. Finally, at the same level of accuracy, our ExpNet is orders of magnitude faster than its alternatives.


title: Modeling polypharmacy side effects with graph convolutional networks

abstract: The use of drug combinations, termed polypharmacy, is common to treat patients with complex diseases and co-existing conditions. However, a major consequence of polypharmacy is a much higher risk of adverse side effects for the patient. Polypharmacy side effects emerge because of drug-drug interactions, in which activity of one drug may change if taken with another drug. The knowledge of drug interactions is limited because these complex relationships are rare, and are usually not observed in relatively small clinical testing. Discovering polypharmacy side effects thus remains an important challenge with significant implications for patient mortality. Here, we present Decagon, an approach for modeling polypharmacy side effects. The approach constructs a multimodal graph of protein-protein interactions, drug-protein target interactions, and the polypharmacy side effects, which are represented as drug-drug interactions, where each side effect is an edge of a different type. Decagon is developed specifically to handle such multimodal graphs with a large number of edge types. Our approach develops a new graph convolutional neural network for multirelational link prediction in multimodal networks. Decagon predicts the exact side effect, if any, through which a given drug combination manifests clinically. Decagon accurately predicts polypharmacy side effects, outperforming baselines by up to 69%. We find that it automatically learns representations of side effects indicative of co-occurrence of polypharmacy in patients. Furthermore, Decagon models particularly well side effects with a strong molecular basis, while on predominantly non-molecular side effects, it achieves good performance because of effective sharing of model parameters across edge types. Decagon creates opportunities to use large pharmacogenomic and patient data to flag and prioritize side effects for follow-up analysis.


title: Neutrino-mass differences studied in a model with one basic neutrino

abstract: In this paper, we propose and study a model, in which there exists only one basic neutrino. Basically, the Lagrangian of the model is obtained from the Lagrangian in the Glashow-Weinberg-Salam electroweak theory by reducing the three flavor neutrinos there to one basic neutrino. In this model, neutrino states with a fixed flavor are interpreted as certain superpositions of states of the basic neutrino and states of the related charged lepton. Neutrino mass states are associated with low-lying eigenstates of the total Hamiltonian. We derive an approximate expression for a ratio of neutrino-mass differences, which gives a value $17$ for the ratio, about half of the experimental result of $33$.


title: Effective model with strong Kitaev interactions for $\alpha$-${\rm RuCl_3}$

abstract: We use an exact numerical diagonalization method to calculate the dynamical spin structure factors (DSFs) of three ab-initio models and one ab-initio-guided model for a honeycomb-lattice magnet $\alpha$-RuCl$_3$. We also use thermal pure quantum states to calculate the temperature dependence of the heat capacity, the nearest-neighbor (NN) spin-spin correlation function, and the static spin structure factor. From the results obtained from these four effective models, we find that, even when the magnetic order is stabilized at low temperature, the intensity at the $\Gamma$ point in the DSFs increases with increasing NN spin correlation. In addition, we find that the four models fail to explain heat-capacity measurements whereas two of the four models succeed in explaining inelastic-neutron-scattering (INS) experiments. In the four models, when temperature decreases, the heat capacity shows a prominent peak at a high temperature where the NN spin-spin correlation function increases. However, the peak temperature in heat capacity is too low in comparison with that observed experimentally. To address these discrepancies, we propose an effective model that includes strong ferromagnetic Kitaev coupling, and we show that this model quantitatively reproduces both INS experiments and heat-capacity measurements. To further examine the adequacy of the proposed model, we calculate the field dependence of the polarized terahertz spectra, which reproduces the experimental results: the spin-gapped excitation survives up to an onset field where the magnetic order disappears and the response in the high-field region is almost linear. Based on these numerical results, we argue that the low-energy magnetic excitation in $\alpha$-RuCl$_3$ is mainly characterized by interactions such as off-diagonal interactions and weak Heisenberg interactions between NN pairs, rather than by the strong Kitaev interactions.


title: Real Time Collision Detection and Identification for Robotic Manipulators

abstract: The majority of everyday tasks involve interacting with unstructured environments. This implies that, in order for robots to be truly useful they must be able to handle contacts. This paper explores how a particle filter can be used to localize a contact point and estimate the external force. We demonstrate the capability of the particle filter on a simulated 4DoF planar robotic arm, and compare it to a well-established analytical approach.


title: Double phase transition of the Ising model in core-periphery networks

abstract: We study the phase transition of the Ising model in networks with core-periphery structures. By Monte Carlo simulations, we show that prior to the order-disorder phase transition the system organizes into an inhomogeneous intermediate phase in which core nodes are much more ordered than peripheral nodes. Interestingly, the susceptibility shows double peaks at two distinct temperatures. We find that, if the connections between core and periphery increase linearly with network size, the first peak does not exhibit any size-dependent effect, and the second one diverges in the limit of infinite network size. Otherwise, if the connections between core and periphery scale sub-linearly with the network size, both peaks of the susceptibility diverge as power laws in the thermodynamic limit. This suggests the appearance of a double transition phenomenon in the Ising model for the latter case. Moreover, we develop a mean-field theory that agrees well with the simulations.


title: Asymptotic behavior of lifetime sums for random simplicial complex processes

abstract: We study the homological properties of random simplicial complexes. In particular, we obtain the asymptotic behavior of lifetime sums for a class of increasing random simplicial complexes; this result is a higher-dimensional counterpart of Frieze's $\zeta(3)$-limit theorem for the Erd\H{o}s-R\'{e}nyi graph process. The main results include solutions to questions posed in an earlier study by Hiraoka and Shirai about the Linial-Meshulam complex process and the random clique complex process. One of the key elements of the arguments is a new upper bound on the Betti numbers of general simplicial complexes in terms of the number of small eigenvalues of Laplacians on links. This bound can be regarded as a quantitative version of the cohomology vanishing theorem.


title: Ultra-high Q terahertz whispering-gallery modes in a silicon resonator

abstract: We report on the first experimental demonstration of terahertz (THz) whispering-gallery modes (WGMs) with an ultra high quality (Q) factor of $1.5 \times {10}^{4}$ at 0.62THz. The WGMs are observed in a high resistivity float zone silicon (HRFZ-Si) spherical resonator coupled to a sub-wavelength silica waveguide. A detailed analysis of the coherent continuous wave (CW) THz spectroscopy measurements combined with a numerical model based on Mie-Debye-Aden-Kerker (MDAK) theory allows to unambiguously identify the observed higher order radial THz WGMs.


title: The morphological evolution, AGN fractions, dust content, environments, and downsizing of massive green valley galaxies at 0.5<z<2.5 in 3D-HST/CANDELS

abstract: To explore the evolutionary connection among red, green, and blue galaxy populations, based on a sample of massive ($M_* > 10^{10} M_{\odot} $) galaxies at 0.5<z<2.5 in five 3D-HST/CANDELS fields, we investigate the dust content, morphologies, structures, AGN fractions, and environments of these three galaxy populations. Green valley galaxies are found to have intermediate dust attenuation, and reside in the middle of the regions occupied by quiescent and star-forming galaxies in the UVJ diagram. Compared with blue and red galaxy populations at z<2, green galaxies have intermediate compactness and morphological parameters such as Sersic index, concentration, Gini coefficient, and the second order moment of the 20% brightest pixels of a galaxy. Above findings seem to favor the scenario that green galaxies are at transitional phase when star-forming galaxies are being quenched into quiescent status. The green galaxies at z<2 show the highest AGN fraction, suggesting that AGN feedback may have played an important role in star formation quenching. For the massive galaxies at 2<z<2.5, both red and green galaxies are found to have a similarly higher AGN fraction than the blue ones, which implies that AGN feedback may help to keep quiescence of red galaxies at z>2. A significant environmental difference is found between green and red galaxies at z<1.5. Green and blue galaxies at z>0.5 seem to have similar local density distributions, suggesting that environment quenching is not the major mechanism to cease star formation at z>0.5. The fractions of three populations as functions of mass support a "downsizing" quenching picture that the bulk of star formation in more massive galaxies is completed earlier than that of lower mass galaxies.


title: Analysis of stationary points and their bifurcations in the ABC flow

abstract: Analytical expressions for coordinates of stationary points and conditions for their existence in the ABC flow are received. The type of the stationary points is shown analytically to be saddle-node. Exact expressions for eigenvalues and eigenvectors of the stability matrix are given. Behavior of the stationary points along the bifurcation lines is described.


title: Best Practices for a Future Open Code Policy: Experiences and Vision of the Astrophysics Source Code Library

abstract: We are members of the Astrophysics Source Code Library's Advisory Committee and its editor-in-chief. The Astrophysics Source Code Library (ASCL, ascl.net) is a successful initiative that advocates for open research software and provides an infrastructure for registering, discovering, sharing, and citing this software. Started in 1999, the ASCL has been expanding in recent years, with an average of over 200 codes added each year, and now houses over 1,600 code entries.


title: Correlated Insulator Behaviour at Half-Filling in Magic Angle Graphene Superlattices

abstract: Van der Waals (vdW) heterostructures are an emergent class of metamaterials comprised of vertically stacked two-dimensional (2D) building blocks, which provide us with a vast tool set to engineer their properties on top of the already rich tunability of 2D materials. One of the knobs, the twist angle between different layers, plays a crucial role in the ultimate electronic properties of a vdW heterostructure and does not have a direct analog in other systems such as MBE-grown semiconductor heterostructures. For small twist angles, the moir\'e pattern produced by the lattice misorientation creates a long-range modulation. So far, the study of the effect of twist angles in vdW heterostructures has been mostly concentrated in graphene/hexagonal boron nitride (h-BN) twisted structures, which exhibit relatively weak interlayer interaction due to the presence of a large bandgap in h-BN. Here we show that when two graphene sheets are twisted by an angle close to the theoretically predicted 'magic angle', the resulting flat band structure near charge neutrality gives rise to a strongly-correlated electronic system. These flat bands exhibit half-filling insulating phases at zero magnetic field, which we show to be a Mott-like insulator arising from electrons localized in the moir\'e superlattice. These unique properties of magic-angle twisted bilayer graphene (TwBLG) open up a new playground for exotic many-body quantum phases in a 2D platform made of pure carbon and without magnetic field. The easy accessibility of the flat bands, the electrical tunability, and the bandwidth tunability though twist angle may pave the way towards more exotic correlated systems, such as unconventional superconductors or quantum spin liquids.


title: Generating Redundant Features with Unsupervised Multi-Tree Genetic Programming

abstract: Recently, feature selection has become an increasingly important area of research due to the surge in high-dimensional datasets in all areas of modern life. A plethora of feature selection algorithms have been proposed, but it is difficult to truly analyse the quality of a given algorithm. Ideally, an algorithm would be evaluated by measuring how well it removes known bad features. Acquiring datasets with such features is inherently difficult, and so a common technique is to add synthetic bad features to an existing dataset. While adding noisy features is an easy task, it is very difficult to automatically add complex, redundant features. This work proposes one of the first approaches to generating redundant features, using a novel genetic programming approach. Initial experiments show that our proposed method can automatically create difficult, redundant features which have the potential to be used for creating high-quality feature selection benchmark datasets.


title: On the Predictive Risk in Misspecified Quantile Regression

abstract: In the present paper we investigate the predictive risk of possibly misspecified quantile regression functions. The in-sample risk is well-known to be an overly optimistic estimate of the predictive risk and we provide two relatively simple (asymptotic) characterizations of the associated bias, also called expected optimism. We propose estimates for the expected optimism and the predictive risk, and establish their uniform consistency under mild conditions. Our results hold for models of moderately growing size and allow the quantile function to be incorrectly specified. Empirical evidence from our estimates is encouraging as it compares favorably with cross-validation.


title: Goethals--Seidel difference families with symmetric or skew base blocks

abstract: We single out a class of difference families which is widely used in some constructions of Hadamard matrices and which we call Goethals--Seidel (GS) difference families. They consist of four subsets (base blocks) of a finite abelian group of order $v$, which can be used to construct Hadamard matrices via the well-known Goethals--Seidel array. We consider the special class of these families in cyclic groups, where each base block is either symmetric or skew. We omit the well-known case where all four blocks are symmetric. By extending previous computations by several authors, we complete the classification of GS-difference families of this type for odd $v<50$. In particular, we have constructed the first examples of so called good matrices, G-matrices and best matrices of order 43, and good matrices and G-matrices of order 45. We also point out some errors in one of the cited references.


title: Nonlinear resonances in the $ABC$-flow

abstract: In this paper we study resonances of the $ABC$-flow in the near integrable case ($C\ll 1$). This is an interesting example of a Hamiltonian system with 3/2 degrees of freedom in which simultaneous existence of two resonances of the same order is possible. Analytical conditions of the resonance existence are received. It is shown numerically that the largest $n:1$ ($n=1,2,3$) resonances exist, and their energies are equal to theoretical energies in the near integrable case. We provide analytical and numerical evidences for existence of two branches of the two largest $n:1$ ($n=1,2$) resonances in the region of finite motion.


title: Biot's parameters estimation in ultrasound propagation through cancellous bone

abstract: Of interest is the characterization of a cancellous bone immersed in an acoustic fluid. The bone is placed between an ultrasonic point source and a receiver. Cancellous bone is regarded as a porous medium saturated with fluid according to Biot's theory. This model is coupled with the fluid in an open pore configuration and solved by means of the Finite Volume Method. Characterization is posed as a Bayesian parameter estimation problem in Biot's model given pressure data collected at the receiver. As a first step we present numerical results in 2D for signal recovery. It is shown that as point estimators, the Conditional Mean outperforms the classical PDE-constrained minimization solution.


title: Unlabelled Sensing: A Sparse Bayesian Learning Approach

abstract: We address the recovery of sparse vectors in an overcomplete, linear and noisy multiple measurement framework, where the measurement matrix is known upto a permutation of its rows. We derive sparse Bayesian learning (SBL) based updates for joint recovery of the unknown sparse vector and the sensing order, represented using a permutation matrix. We model the sparse matrix using multiple uncorrelated and correlated vectors, and in particular, we use the first order AR model for the correlated sparse vectors. We propose the Permutation-MSBL and a Kalman filtering based Permutation-KSBL algorithm for low-complex joint recovery of the uncorrelated and correlated sparse vectors, jointly with the permutation matrix. The novelty of this work emerges in providing a simple update step for the permutation matrix using the rearrangement inequality. We demonstrate the mean square error and the permutation recovery performance of the proposed algorithms using Monte Carlo simulations.


title: Interpretable Deep Convolutional Neural Networks via Meta-learning

abstract: Model interpretability is a requirement in many applications in which crucial decisions are made by users relying on a model's outputs. The recent movement for "algorithmic fairness" also stipulates explainability, and therefore interpretability of learning models. And yet the most successful contemporary Machine Learning approaches, the Deep Neural Networks, produce models that are highly non-interpretable. We attempt to address this challenge by proposing a technique called CNN-INTE to interpret deep Convolutional Neural Networks (CNN) via meta-learning. In this work, we interpret a specific hidden layer of the deep CNN model on the MNIST image dataset. We use a clustering algorithm in a two-level structure to find the meta-level training data and Random Forest as base learning algorithms to generate the meta-level test data. The interpretation results are displayed visually via diagrams, which clearly indicates how a specific test instance is classified. Our method achieves global interpretation for all the test instances without sacrificing the accuracy obtained by the original deep CNN model. This means our model is faithful to the deep CNN model, which leads to reliable interpretations.


title: Block4Forensic: An Integrated Lightweight Blockchain Framework for Forensics Applications of Connected Vehicles

abstract: Today's vehicles are becoming cyber-physical systems that do not only communicate with other vehicles but also gather various information from hundreds of sensors within them. These developments help create smart and connected (e.g., self-driving) vehicles that will introduce significant information to drivers, manufacturers, insurance companies and maintenance service providers for various applications. One such application that is becoming crucial with the introduction of self-driving cars is the forensic analysis for traffic accidents. The utilization of vehicle-related data can be instrumental in post-accident scenarios to find out the faulty party, particularly for self-driving vehicles. With the opportunity of being able to access various information on the cars, we propose a permissioned blockchain framework among the various elements involved to manage the collected vehicle-related data. Specifically, we first integrate Vehicular Public Key Management (VPKI) to the proposed blockchain to provide membership establishment and privacy. Next, we design a fragmented ledger that will store detailed data related to vehicle such as maintenance information/history, car diagnosis reports, etc. The proposed forensic framework enables trustless, traceable and privacy-aware post-accident analysis with minimal storage and processing overhead.


title: Optimal interpolation formulas in $W_2^{(m,m-1)}$ space

abstract: In the present paper optimal interpolation formulas are constructed in $W_2^{(m,m-1)}(0,1)$ space. Explicit formulas for coefficients of optimal interpolation formulas are obtained. Some numerical results are presented.


title: A Solar cycle correlation of coronal element abundances in Sun-as-a-star observations

abstract: The elemental composition in the coronae of low-activity solar-like stars appears to be related to fundamental stellar properties such as rotation, surface gravity, and spectral type. Here we use full-Sun observations from the Solar Dynamics Observatory, to show that when the Sun is observed as a star, the variation of coronal composition is highly correlated with a proxy for solar activity, the F10.7 cm radio flux, and therefore with the solar cycle phase. Similar cyclic variations should therefore be detectable spectroscopically in X-ray observations of solar analogs. The plasma composition in full-disk observations of the Sun is related to the evolution of coronal magnetic field activity. Our observations therefore introduce an uncertainty into the nature of any relationship between coronal composition and fixed stellar properties. The results highlight the importance of systematic full-cycle observations for understanding the elemental composition of solar-like stellar coronae.


title: Comparison of permutationally invariant polynomials, neural networks, and Gaussian approximation potentials in representing water interactions through many-body expansions

abstract: The accurate representation of multidimensional potential energy surfaces is a necessary requirement for realistic computer simulations of molecular systems. The continued increase in computer power accompanied by advances in correlated electronic structure methods nowadays enable routine calculations of accurate interaction energies for small systems, which can then be used as references for the development of analytical potential energy functions (PEFs) rigorously derived from many-body expansions. Building on the accuracy of the MB-pol many-body PEF, we investigate here the performance of permutationally invariant polynomials, neural networks, and Gaussian approximation potentials in representing water two-body and three-body interaction energies, denoting the resulting potentials PIP-MB-pol, BPNN-MB-pol, and GAP-MB-pol, respectively. Our analysis shows that all three analytical representations exhibit similar levels of accuracy in reproducing both two-body and three-body reference data as well as interaction energies of small water clusters obtained from calculations carried out at the coupled cluster level of theory, the current gold standard for chemical accuracy. These results demonstrate the synergy between interatomic potentials formulated in terms of a many-body expansion, such as MB-pol, that are physically sound and transferable, and machine-learning techniques that provide a flexible framework to approximate the short-range interaction energy terms.


title: Detecting Zones and Threat on 3D Body for Security in Airports using Deep Machine Learning

abstract: In this research, it was used a segmentation and classification method to identify threat recognition in human scanner images of airport security. The Department of Homeland Security's (DHS) in USA has a higher false alarm, produced from theirs algorithms using today's scanners at the airports. To repair this problem they started a new competition at Kaggle site asking the science community to improve their detection with new algorithms. The dataset used in this research comes from DHS at https://www.kaggle.com/c/passenger-screening-algorithm-challenge/data According to DHS: "This dataset contains a large number of body scans acquired by a new generation of millimeter wave scanner called the High Definition-Advanced Imaging Technology (HD-AIT) system. They are comprised of volunteers wearing different clothing types (from light summer clothes to heavy winter clothes), different body mass indices, different genders, different numbers of threats, and different types of threats". Using Python as a principal language, the preprocessed of the dataset images extracted features from 200 bodies using: intensity, intensity differences and local neighbourhood to detect, to produce segmentation regions and label those regions to be used as a truth in a training and test dataset. The regions are subsequently give to a CNN deep learning classifier to predict 17 classes (that represents the body zones): zone1, zone2, ... zone17 and zones with threat in a total of 34 zones. The analysis showed the results of the classifier an accuracy of 98.2863% and a loss of 0.091319, as well as an average of 100% for recall and precision.


title: Orbital Alignment of Main-Belt Comets

abstract: We examine the orbital element distribution of main-belt comets (MBCs), which are objects that exhibit cometary activity yet orbit in the main asteroid belt, and may be potentially useful as tracers of ice in the inner solar system. We find that the currently known and currently active MBCs have remarkably similar longitudes of perihelion, which are also aligned with that of Jupiter. The clustered objects have significantly higher current osculating eccentricities relative to their proper eccentricities, consistent with their orbits being currently, though only temporarily, secularly excited in osculating eccentricity due to Jupiter's influence. At the moment, most MBCs seem to have current osculating elements that may be particularly favorable for the object becoming active (e.g., maybe because of higher perihelion temperatures or higher impact velocities causing an effective increase in the size of the potential triggering impactor population). At other times, other icy asteroids will have those favorable conditions and might become MBCs at those times as well.


title: Three-dimensional structure of clumpy outflow from supercritical accretion flow onto black holes

abstract: We perform global three-dimensional (3D) radiation-hydrodynamic (RHD) simulations of out- flow from supercritical accretion flow around a 10 Msun black hole. We only solve the outflow part, starting from the axisymmetric 2D simulation data in a nearly steady state but with small perturbations in a sinusoidal form being added in the azimuthal direction. The mass accretion rate onto the black hole is ~10^2 L_E/c^2 in the underlying 2D simulation data and the outflow rate is ~10 L_E/c^2 (with LE and c being the Eddington luminosity and speed of light, respectively). We first confirm the emergence of clumpy outflow, which was discovered by the 2D RHD simulations, above the photosphere located at a few hundreds of Schwarzschild radii (r_S) from the central black hole. As prominent 3D features we find that the clumps have the shape of a torn sheet, rather than a cut string, and that they are rotating around the central black hole with a sub-Keplerian velocity at a distance of ~10^3 r_S from the center. The typical clump size is ~30 r_S or less in the radial direction, and is more elongated in the angular directions, ~hundreds of r_S at most. The sheet separation ranges from 50 to 150 r_S. We expect stochastic time variations when clumps pass across the line of the sight of a distant observer. Variation timescales are estimated to be several seconds for a black hole with mass of ten to several tens of Msun, in rough agreement with the observations of some ultra-luminous X-ray sources.


title: An Instability in Variational Inference for Topic Models

abstract: Topic models are Bayesian models that are frequently used to capture the latent structure of certain corpora of documents or images. Each data element in such a corpus (for instance each item in a collection of scientific articles) is regarded as a convex combination of a small number of vectors corresponding to `topics' or `components'. The weights are assumed to have a Dirichlet prior distribution. The standard approach towards approximating the posterior is to use variational inference algorithms, and in particular a mean field approximation. We show that this approach suffers from an instability that can produce misleading conclusions. Namely, for certain regimes of the model parameters, variational inference outputs a non-trivial decomposition into topics. However --for the same parameter values-- the data contain no actual information about the true decomposition, and hence the output of the algorithm is uncorrelated with the true topic decomposition. Among other consequences, the estimated posterior mean is significantly wrong, and estimated Bayesian credible regions do not achieve the nominal coverage. We discuss how this instability is remedied by more accurate mean field approximations.