title: Stability of Metabolic Networks via Linear-In-Flux-Expressions

abstract: The methodology named LIFE (Linear-in-Flux-Expressions) was developed with the purpose of simulating and analyzing large metabolic systems. With LIFE, the number of model parameters is reduced by accounting for correlations among the parameters of the system. Perturbation analysis on LIFE systems results in less overall variability of the system, leading to results that more closely resemble empirical data. These systems can be associated to graphs, and characteristics of the graph give insight into the dynamics of the system. This work addresses two main problems: 1. for fixed metabolite levels, find all fluxes for which the metabolite levels are an equilibrium, and 2. for fixed fluxes, find all metabolite levels which are equilibria for the system. We characterize the set of solutions for both problems, and show general results relating stability of systems to the structure of the associated graph. We show that there is a structure of the graph necessary for stable dynamics. Along with these general results, we show how stability analysis from the fields of network flows, compartmental systems, control theory and Markov chains apply to LIFE systems.


title: Renormalized oscillation theory for linear Hamiltonian systems on [0,1] via the Maslov index

abstract: Working with a general class of linear Hamiltonian systems on $[0, 1]$, we show that renormalized oscillation results can be obtained in a natural way through consideration of the Maslov index associated with appropriately chosen paths of Lagrangian subspaces of $\mathbb{C}^{2n}$. We verify that our applicability class includes Dirac and Sturm-Liouville systems, as well as a system arising from differential-algebraic equations for which the spectral parameter appears nonlinearly.


title: Towards efficient random metasurfaces

abstract: Random media introduce large degrees of freedom in device design and can thus address challenges in manipulating optical waves. Wave shaping with metasurfaces has mainly utilized periodic or quasi-periodic grids, and, the potential of random arrangement of particles for devices has only come under investigation recently. The main difficulty in pursuing random metasurfaces is the identification of the degrees of freedom that optimize their efficiencies and functions. They can also encode information using the statistics of particles distribution. We propose a phase-map that accounts for the statistical nature of random media. The method takes into account effects of random near-field couplings that introduce phase errors by affecting the phase shift of elements. The proposed approach increases the efficiency of our random metasurface devices by up to ~20%. This work paves the way towards the efficient design of random metasurfaces with potential applications in highly secure optical cryptography and information encoding.


title: A Visual Attention Grounding Neural Model for Multimodal Machine Translation

abstract: We introduce a novel multimodal machine translation model that utilizes parallel visual and textual information. Our model jointly optimizes the learning of a shared visual-language embedding and a translator. The model leverages a visual attention grounding mechanism that links the visual semantics with the corresponding textual semantics. Our approach achieves competitive state-of-the-art results on the Multi30K and the Ambiguous COCO datasets. We also collected a new multilingual multimodal product description dataset to simulate a real-world international online shopping scenario. On this dataset, our visual attention grounding model outperforms other methods by a large margin.


title: Structural symmetry and membrane curvature sensing

abstract: Symmetry is closely intertwined with the function, genetics, and chemical properties of multiprotein complexes. Here, we explore the relation between structural symmetry and the ability of membrane proteins to sense and induce membrane curvature, which is a key factor for modulating the shape and organization of cell membranes. Using coarse-grained simulations and elasticity theory, we show that the potential for direction-dependent membrane curvature sensing is limited to asymmetric proteins, dimers, and tetramers, and argue that one should expect this anisotropy to be strongest for dimers. Odd and higher-order symmetries strongly suppress directional curvature sensing. This classification gives a new perspective on the structure-function relation for membrane proteins, and simplifies the task of translating between molecular sensing mechanisms and their large-scale cellular consequences.


title: Learning Models for Shared Control of Human-Machine Systems with Unknown Dynamics

abstract: We present a novel approach to shared control of human-machine systems. Our method assumes no a priori knowledge of the system dynamics. Instead, we learn both the dynamics and information about the user's interaction from observation through the use of the Koopman operator. Using the learned model, we define an optimization problem to compute the optimal policy for a given task, and compare the user input to the optimal input. We demonstrate the efficacy of our approach with a user study. We also analyze the individual nature of the learned models by comparing the effectiveness of our approach when the demonstration data comes from a user's own interactions, from the interactions of a group of users and from a domain expert. Positive results include statistically significant improvements on task metrics when comparing a user-only control paradigm with our shared control paradigm. Surprising results include findings that suggest that individualizing the model based on a user's own data does not effect the ability to learn a useful dynamic system. We explore this tension as it relates to developing human-in-the-loop systems further in the discussion.


title: The Unexpected Role of Evolving Longitudinal Electric Fields in Generating Energetic Electrons in Relativistically Transparent Plasmas

abstract: Superponderomotive-energy electrons are observed experimentally from the interaction of an intense laser pulse with a relativistically transparent target. For a relativistically transparent target, kinetic modeling shows that the generation of energetic electrons is dominated by energy transfer within the main, classically overdense, plasma volume. The laser pulse produces a narrowing, funnel-like channel inside the plasma volume that generates a field structure responsible for the electron heating. The field structure combines a slowly evolving azimuthal magnetic field, generated by a strong laser-driven longitudinal electron current, and, unexpectedly, a strong propagating longitudinal electric field, generated by reflections off the walls of the funnel-like channel. The magnetic field assists electron heating by the transverse electric field of the laser pulse through deflections, whereas the longitudinal electric field directly accelerates the electrons in the forward direction. The longitudinal electric field produced by reflections is 30 times stronger than that in the incoming laser beam and the resulting direct laser acceleration contributes roughly one third of the energy transferred by the transverse electric field of the laser pulse to electrons of the super-ponderomotive tail.


title: Robust Text Classifier on Test-Time Budgets

abstract: We propose a generic and interpretable learning framework for building robust text classification model that achieves accuracy comparable to full models under test-time budget constraints. Our approach learns a selector to identify words that are relevant to the prediction tasks and passes them to the classifier for processing. The selector is trained jointly with the classifier and directly learns to incorporate with the classifier. We further propose a data aggregation scheme to improve the robustness of the classifier. Our learning framework is general and can be incorporated with any type of text classification model. On real-world data, we show that the proposed approach improves the performance of a given classifier and speeds up the model with a mere loss in accuracy performance.


title: An elementary introduction to information geometry

abstract: In this survey, we describe the fundamental differential-geometric structures of information manifolds, state the fundamental theorem of information geometry, and illustrate some use cases of these information manifolds in information sciences. The exposition is self-contained by concisely introducing the necessary concepts of differential geometry, but proofs are omitted for brevity.


title: Probabilistic Model of Object Detection Based on Convolutional Neural Network

abstract: The combination of a CNN detector and a search framework forms the basis for local object/pattern detection. To handle the waste of regional information and the defective compromise between efficiency and accuracy, this paper proposes a probabilistic model with a powerful search framework. By mapping an image into a probabilistic distribution of objects, this new model gives more informative outputs with less computation. The setting and analytic traits are elaborated in this paper, followed by a series of experiments carried out on FDDB, which show that the proposed model is sound, efficient and analytic.


title: Improving Breast Cancer Detection using Symmetry Information with Deep Learning

abstract: Convolutional Neural Networks (CNN) have had a huge success in many areas of computer vision and medical image analysis. However, there is still an immense potential for performance improvement in mammogram breast cancer detection Computer-Aided Detection (CAD) systems by integrating all the information that the radiologist utilizes, such as symmetry and temporal data. In this work, we proposed a patch based multi-input CNN that learns symmetrical difference to detect breast masses. The network was trained on a large-scale dataset of 28294 mammogram images. The performance was compared to a baseline architecture without symmetry context using Area Under the ROC Curve (AUC) and Competition Performance Metric (CPM). At candidate level, AUC value of 0.933 with 95% confidence interval of [0.920, 0.954] was obtained when symmetry information is incorporated in comparison with baseline architecture which yielded AUC value of 0.929 with [0.919, 0.947] confidence interval. By incorporating symmetrical information, although there was no a significant candidate level performance again (p = 0.111), we have found a compelling result at exam level with CPM value of 0.733 (p = 0.001). We believe that including temporal data, and adding benign class to the dataset could improve the detection performance.


title: Can we leverage rating patterns from traditional users to enhance recommendations for children?

abstract: Recommender algorithms performance is often associated with the availability of sufficient historical rating data. Unfortunately, when it comes to children, this data is seldom available. In this paper, we report on an initial analysis conducted to examine the degree to which data about traditional users, i.e., adults, can be leveraged to enhance the recommendation process for children.


title: Binary Image Features Proposed to Empower Computer Vision

abstract: This literature has proposed three fast and easy computable image features to improve computer vision by offering more human-like vision power. These features are not based on image pixels absolute or relative intensity; neither based on shape or colour. So, no complex pixel by pixel calculation is required. For human eyes, pixel by pixel calculation is like seeing an image with maximum zoom which is done only when a higher level of details is required. Normally, first we look at an image to get an overall idea about it to know whether it deserves further investigation or not. This capacity of getting an idea at a glance is analysed and three basic features are proposed to empower computer vision. Potential of proposed features is tested and established through different medical dataset. Achieved accuracy in classification demonstrates possibilities and potential of the use of the proposed features in image processing.


title: Galaxy formation in the Planck Millennium: the atomic hydrogen content of dark matter halos

abstract: We present recalibrations of the GALFORM semi-analytical model of galaxy formation in a new N-body simulation with the Planck cosmology. The Planck Millennium simulation uses more than 128 billion particles to resolve the matter distribution in a cube of $800$ Mpc on a side, which contains more than 77 million dark matter haloes with mass greater than $2.12 \times 10^{9} h^{-1} {\rm M_{\odot}}$ at the present day. Only minor changes to a very small number of model parameters are required in the recalibration. We present predictions for the atomic hydrogen content (HI) of dark matter halos, which is a key input into the calculation of the HI intensity mapping signal expected from the large-scale structure of the Universe. We find that the HI mass $-$ halo mass relation displays a clear break at the halo mass above which AGN heating suppresses gas cooling, $\approx 3 \times 10^{11} h^{-1} M_{\rm \odot}$. Below this halo mass, the HI content of haloes is dominated by the central galaxy; above this mass it is the combined HI content of satellites that prevails. We find that the HI mass - halo mass relation changes little with redshift up to $z=3$. The bias of HI sources shows a scale dependence that gets more pronounced with increasing redshift.


title: Deep Mask For X-ray Based Heart Disease Classification

abstract: We build a deep learning model to detect and classify heart disease using $X-ray$. We collect data from several hospitals and public datasets. After preprocess we get 3026 images including disease type VSD, ASD, TOF and normal control. The main problem we have to solve is to enable the network to accurately learn the characteristics of the heart, to ensure the reliability of the network while increasing accuracy. By learning the doctor's diagnostic experience, labeling the image and using tools to extract masks of heart region, we train a U-net to generate a mask to give more attention. It forces the model to focus on the characteristics of the heart region and obtain more reliable results.


title: Variational Procedure for Higher-Derivative Mechanical Models in a Fractional Integral Framework

abstract: We present both the Lagrangian and Hamiltonian procedures for treating higher-order equations of motion for mechanical models by adopting the Riemann-Liouville Fractional integral to describe their action. We point out and discuss its efficacy and difficulties. We also present the physical and geometric interpretations for the approach we pursue and present the details of a higher-order harmonic oscillator.


title: Nuclei Detection Using Mixture Density Networks

abstract: Nuclei detection is an important task in the histology domain as it is a main step toward further analysis such as cell counting, cell segmentation, study of cell connections, etc. This is a challenging task due to the complex texture of histology image, variation in shape, and touching cells. To tackle these hurdles, many approaches have been proposed in the literature where deep learning methods stand on top in terms of performance. Hence, in this paper, we propose a novel framework for nuclei detection based on Mixture Density Networks (MDNs). These networks are suitable to map a single input to several possible outputs and we utilize this property to detect multiple seeds in a single image patch. A new modified form of a cost function is proposed for training and handling patches with missing nuclei. The probability maps of the nuclei in the individual patches are next combined to generate the final image-wide result. The experimental results show the state-of-the-art performance on complex colorectal adenocarcinoma dataset.


title: Deep multiscale convolutional feature learning for weakly supervised localization of chest pathologies in X-ray images

abstract: Localization of chest pathologies in chest X-ray images is a challenging task because of their varying sizes and appearances. We propose a novel weakly supervised method to localize chest pathologies using class aware deep multiscale feature learning. Our method leverages intermediate feature maps from CNN layers at different stages of a deep network during the training of a classification model using image level annotations of pathologies. During the training phase, a set of \emph{layer relevance weights} are learned for each pathology class and the CNN is optimized to perform pathology classification by convex combination of feature maps from both shallow and deep layers using the learned weights. During the test phase, to localize the predicted pathology, the multiscale attention map is obtained by convex combination of class activation maps from each stage using the \emph{layer relevance weights} learned during the training phase. We have validated our method using 112000 X-ray images and compared with the state-of-the-art localization methods. We experimentally demonstrate that the proposed weakly supervised method can improve the localization performance of small pathologies such as nodule and mass while giving comparable performance for bigger pathologies e.g., Cardiomegaly


title: High Quality Facial Surface and Texture Synthesis via Generative Adversarial Networks

abstract: In the past several decades, many attempts have been made to model synthetic realistic geometric data. The goal of such models is to generate plausible 3D geometries and textures. Perhaps the best known of its kind is the linear 3D morphable model (3DMM) for faces. Such models can be found at the core of many computer vision applications such as face reconstruction, recognition and authentication to name just a few. Generative adversarial networks (GANs) have shown great promise in imitating high dimensional data distributions. State of the art GANs are capable of performing tasks such as image to image translation as well as auditory and image signal synthesis, producing novel plausible samples from the data distribution at hand. Geometric data is generally more difficult to process due to the inherent lack of an intrinsic parametrization. By bringing geometric data into an aligned space, we are able to map the data onto a 2D plane using a universal parametrization. This alignment process allows for efficient processing of digitally scanned geometric data via image processing tools. Using this methodology, we propose a novel face synthesis model for generation of realistic facial textures together with their corresponding geometry. A GAN is employed in order to imitate the space of parametrized human textures, while corresponding facial geometries are generated by learning the best 3DMM coefficients for each texture. The generated textures are mapped back onto the corresponding geometries to obtain new generated high resolution 3D faces.


title: Controlling Over-generalization and its Effect on Adversarial Examples Generation and Detection

abstract: Convolutional Neural Networks (CNNs) significantly improve the state-of-the-art for many applications, especially in computer vision. However, CNNs still suffer from a tendency to confidently classify out-distribution samples from unknown classes into pre-defined known classes. Further, they are also vulnerable to adversarial examples. We are relating these two issues through the tendency of CNNs to over-generalize for areas of the input space not covered well by the training set. We show that a CNN augmented with an extra output class can act as a simple yet effective end-to-end model for controlling over-generalization. As an appropriate training set for the extra class, we introduce two resources that are computationally efficient to obtain: a representative natural out-distribution set and interpolated in-distribution samples. To help select a representative natural out-distribution set among available ones, we propose a simple measurement to assess an out-distribution set's fitness. We also demonstrate that training such an augmented CNN with representative out-distribution natural datasets and some interpolated samples allows it to better handle a wide range of unseen out-distribution samples and black-box adversarial examples without training it on any adversaries. Finally, we show that generation of white-box adversarial attacks using our proposed augmented CNN can become harder, as the attack algorithms have to get around the rejection regions when generating actual adversaries.


title: Experimental Study of Beam Dynamics In The PIP-II MEBT Prototype

abstract: The Proton Improvement Plan, Stage Two (PIP-II) is a program of upgrades proposed for the Fermilab injection complex, which central part is an 800 MeV, 2 mA CW SRF linac. A prototype of the PIP-II linac front end called PIP-II Injector Test (PIP2IT) is being built at Fermilab. As of now, a 15 mA DC, 30-keV H- ion source, a 2 m-long Low Energy Beam Transport (LEBT), a 2.1 MeV CW RFQ, followed by a 10 m Medium Energy Beam Transport (MEBT) have been assembled and commissioned. The MEBT bunch-by-bunch chopping system and the requirement of a low uncontrolled beam loss put stringent limitations on the beam envelope and its variation. Measurements of transverse and longitudinal beam dynamics in the MEBT were performed in the range of 1-10 mA of the RFQ beam current. Almost all measurements are made with 10 {\mu}s beam pulses in order to avoid damage to the beam line. This report presents measurements of the transverse optics with differential trajectories, reconstruction of the beam envelope with scrapers and an Allison emittance scanner, as well as bunch length measurements with a Fast Faraday Cup.


title: The Strong Slope Conjecture and torus knots

abstract: We observe that the strong slope conjecture implies that the degree of the colored Jones polynomial detects all torus knots. As an application we obtain that an adequate knot that has the same colored Jones polynomial degrees as a torus knot must be a $(2,q)$-torus knot.


title: Sublattice symmetry breaking and Kondo-effect enhancement in strained graphene

abstract: Kondo physics in doped monolayer graphene is predicted to exhibit unusual features due to the linear vanishing of the pristine material's density of states at the Dirac point. Despite several attempts, conclusive experimental observation of the phenomenon remains elusive. One likely obstacle to identification is a very small Kondo temperature scale $T_K$ in situations where the chemical potential lies near the Dirac point. We propose tailored mechanical deformations of monolayer graphene as a means of revealing unique fingerprints of the Kondo effect. Inhomogeneous strains are known to produce specific alternating changes in the local density of states (LDOS) away from the Dirac point that signal sublattice symmetry breaking effects. Small LDOS changes can be amplified in an exponential increase or decrease of $T_K$ for magnetic impurities attached at different locations. We illustrate this behavior in two deformation geometries: a circular 'bubble' and a long fold, both described by Gaussian displacement profiles. We calculate the LDOS changes for modest strains and analyze the relevant Anderson impurity model describing a magnetic atom adsorbed in either a 'top-site' or a 'hollow-site' configuration. Numerical renormalization-group solutions of the impurity model suggest that higher expected $T_K$ values, combined with distinctive spatial patterns under variation of the point of graphene attachment, make the top-site configuration the more promising for experimental observation of signatures of the Kondo effect. The strong strain sensitivity of $T_K$ may lift top-site Kondo physics into the range experimentally accessible using local probes such as scanning tunneling microscopy.


title: Probabilistic Graphical Modeling approach to dynamic PET direct parametric map estimation and image reconstruction

abstract: In the context of dynamic emission tomography, the conventional processing pipeline consists of independent image reconstruction of single time frames, followed by the application of a suitable kinetic model to time activity curves (TACs) at the voxel or region-of-interest level. The relatively new field of 4D PET direct reconstruction, by contrast, seeks to move beyond this scheme and incorporate information from multiple time frames within the reconstruction task. Existing 4D direct models are based on a deterministic description of voxels' TACs, captured by the chosen kinetic model, considering the photon counting process the only source of uncertainty. In this work, we introduce a new probabilistic modeling strategy based on the key assumption that activity time course would be subject to uncertainty even if the parameters of the underlying dynamic process were known. This leads to a hierarchical Bayesian model, which we formulate using the formalism of Probabilistic Graphical Modeling (PGM). The inference of the joint probability density function arising from PGM is addressed using a new gradient-based iterative algorithm, which presents several advantages compared to existing direct methods: it is flexible to an arbitrary choice of linear and nonlinear kinetic model; it enables the inclusion of arbitrary (sub)differentiable priors for parametric maps; it is simpler to implement and suitable to integration in computing frameworks for machine learning. Computer simulations and an application to real patient scan showed how the proposed approach allows us to weight the importance of the kinetic model, providing a bridge between indirect and deterministic direct methods.

categories: stat.AP physics.data-an physics.med-ph stat.ML


title: Convergence of the Augmented Decomposition Algorithm

abstract: We study the convergence of the Augmented Decomposition Algorithm (ADA) proposed in [32] for solving multi-block separable convex minimization problems subject to linear constraints. We show that the global convergence rate of the exact ADA is o(1/k) under the assumption that there exists a saddle point. We consider the inexact Augmented Decomposition Algorithm (iADA) and establish global and local convergence results under some mild assumptions, by providing a stability result for the maximal monotone operator T associated with the perturbation from both primal and dual perspectives. This result implies the local linear convergence of the inexact ADA for many applications such as the lasso, total variation reconstruction, exchange problem and many other problems from statistics, machine learning and engineering with l1 regularization.


title: Defining the problem of Observation Learning

abstract: This article defines and formulates the problem of observation learning in robotic systems.


title: Performance evaluation of job schedulers on Hadoop YARN

abstract: To solve the limitation of Hadoop on scalability, resource sharing, and application support, the open-source community proposes the next generation of Hadoop's compute platform called Yet Another Resource Negotiator (YARN) by separating resource management functions from the programming model. This separation enables various application types to run on YARN in parallel. To achieve fair resource sharing and high resource utilization, YARN provides the capacity scheduler and the fair scheduler. However, the performance impacts of the two schedulers are not clear when mixed applications run on a YARN cluster. Therefore, in this paper, we study four scheduling-policy combinations (SPCs for short) derived from the two schedulers and then evaluate the four SPCs in extensive scenarios, which consider not only four application types, but also three different queue structures for organizing applications. The experimental results enable YARN managers to comprehend the influences of different SPCs and different queue structures on mixed applications. The results also help them to select a proper SPC and an appropriate queue structure to achieve better application execution performance.


title: Generalized exponential basis for efficient solving of homogeneous diffusion free boundary problems: Russian option pricing

abstract: This paper develops a method for solving free boundary problems for time-homogeneous diffusions. We combine the complete exponential system of solutions for the heat equation, transmutation operators and recently discovered Neumann series of Bessel functions representation for solutions of Sturm-Liouville equations to construct a complete system of solutions for the considered partial differential equations. The conceptual algorithm for the application of the method is presented. The valuation of Russian options with finite horizon is used as a numerical illustration. The solution under different horizons is computed and compared to the results that appear in the literature.


title: Fibrations of AU-contexts beget fibrations of toposes

abstract: Suppose an extension map $U\colon \mathbb{T}_1 \to \mathbb{T}_0$ in the 2-category $\mathfrak{Con}$ of contexts for arithmetic universes satisfies a Chevalley criterion for being an (op)fibration in $\mathfrak{Con}$. If $M$ is a model of $\mathbb{T}_0$ in an elementary topos $\mathcal{S}$ with nno, then the classifier $p\colon\mathcal{S}[\mathbb{T}_1/M]\to\mathcal{S}$ satisfies Johnstone's criterion for being an (op)fibration in the 2-category $\mathcal{E}\mathfrak{Top}$ of elementary toposes (with nno) and geometric morphisms. Along the way, we provide a convenient reformulation of Johnstone's criterion.


title: Corrigendum: Similarity degree of Fourier algebras

abstract: We address two errors made in our paper arXiv:1511.03423. The most significant error is in Theorem 1.1. We repair this error, and show that the main result, Theorem 2.5 of arXiv:1511.03423, is true. The second error is in one of our examples, Remark 2.4 (iv), and we partially resolve it.


title: Delay-Induced Switched States in a slow-fast system

abstract: We consider the two-component delay system $\varepsilon x^{\prime}(t)=-x(t)-y(t)+f(x(t-1)),$ $y^{\prime}(t)=\eta x(t)$ with small parameters $\varepsilon,\eta$, and positive feedback function $f$. Previously, such systems have been reported to model switching in optoelectronic experiments, where each switching induces another one after approximately one delay time, related to one round trip of the signal. In this paper, we study these delay-induced switched states. We provide conditions for their existence and show how the formal limits $\varepsilon\to0$ and/or $\eta\to0$ facilitate our understanding of this phenomenon.


title: Unknown Examples & Machine Learning Model Generalization

abstract: Over the past decades, researchers and ML practitioners have come up with better and better ways to build, understand and improve the quality of ML models, but mostly under the key assumption that the training data is distributed identically to the testing data. In many real-world applications, however, some potential training examples are unknown to the modeler, due to sample selection bias or, more generally, covariate shift, i.e., a distribution shift between the training and deployment stage. The resulting discrepancy between training and testing distributions leads to poor generalization performance of the ML model and hence biased predictions. We provide novel algorithms that estimate the number and properties of these unknown training examples---unknown unknowns. This information can then be used to correct the training set, prior to seeing any test data. The key idea is to combine species-estimation techniques with data-driven methods for estimating the feature values for the unknown unknowns. Experiments on a variety of ML models and datasets indicate that taking the unknown examples into account can yield a more robust ML model that generalizes better.


title: A Phonon Boltzmann Study of Microscale Thermal Transport in $\alpha$-RDX Cook-Off

abstract: The microscale thermal transport properties of $\alpha$RDX are believed to be major factors in the initiation process. In this study we present a thorough examination of phonon properties which dominate energy storage and transport in $\alpha$RDX. The phonon lifetimes are determined for all phonon branches, revealing the characteristic time scale of energy transfer amongst phonon modes. The phonon parameters also serve as inputs to a full Brillouin zone three dimensional phonon transport simulation in the presence of a hotspot. In addition to identifying the phonon mode contributions to thermal transport, and as N-N bond breaking is integral to disassociation, we identify phonon modes corresponding to large N-N bond stretch analyzing the manner in which these modes store and transfer energy.


title: Brain Biomarker Interpretation in ASD Using Deep Learning and fMRI

abstract: Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder. Finding the biomarkers associated with ASD is extremely helpful to understand the underlying roots of the disorder and can lead to earlier diagnosis and more targeted treatment. Although Deep Neural Networks (DNNs) have been applied in functional magnetic resonance imaging (fMRI) to identify ASD, understanding the data-driven computational decision making procedure has not been previously explored. Therefore, in this work, we address the problem of interpreting reliable biomarkers associated with identifying ASD; specifically, we propose a 2-stage method that classifies ASD and control subjects using fMRI images and interprets the saliency features activated by the classifier. First, we trained an accurate DNN classifier. Then, for detecting the biomarkers, different from the DNN visualization works in computer vision, we take advantage of the anatomical structure of brain fMRI and develop a frequency-normalized sampling method to corrupt images. Furthermore, in the ASD vs. control subjects classification scenario, we provide a new approach to detect and characterize important brain features into three categories. The biomarkers we found by the proposed method are robust and consistent with previous findings in the literature. We also validate the detected biomarkers by neurological function decoding and comparing with the DNN activation maps.


title: AuthPDB: Query Authentication for Outsourced Probabilistic Databases

abstract: Spurred by developments such as cloud computing, there are increasing efforts for outsourcing of data management. A company (data owner) who lacks expertise and comptational resources can outsource his data to a third-party service provider (server), who provides storage and query evaluation on the outsourced data as the services. One of the security concerns of the outsourcing paradigm is the integrity of the returned query results on the outsourced data. In this paper, we consider the outsourcing of probabilistic databases, on which query evaluation is of high complexity. A dishonest server may return cheap (and incorrect) query answers, hoping that the client who has weak computational power cannot catch the incorrect results. To address this issue, we design efficient integrity verification methods for both all-answer and top-k query evaluation on outsourced probabilistic databases. Our empirical results demonstrate the effectiveness and efficiency of our verification methods.


title: Direct Splitting Method for the Baum-Connes Conjecture

abstract: We introduce a new method for studying the Baum-Connes conjecture, which we call the direct splitting method. The method can simplify and clarify proofs of some of the known cases of the conjecture. In a separate paper, with J. Brodzki, E. Guentner and N. Higson, a similar idea will be used to give a finite-dimensional proof of the Baum-Connes conjecture for groups which act properly and co-compactly on a finite-dimensional CAT(0)-cubical space.


title: Harnessing Infant Cry for swift, cost-effective Diagnosis of Perinatal Asphyxia in low-resource settings

abstract: Perinatal Asphyxia is one of the top three causes of infant mortality in developing countries, resulting to the death of about 1.2 million newborns every year. At its early stages, the presence of asphyxia cannot be conclusively determined visually or via physical examination, but by medical diagnosis. In resource-poor settings, where skilled attendance at birth is a luxury, most cases only get detected when the damaging consequences begin to manifest or worse still, after death of the affected infant. In this project, we explored the approach of machine learning in developing a low-cost diagnostic solution. We designed a support vector machine-based pattern recognition system that models patterns in the cries of known asphyxiating infants (and normal infants) and then uses the developed model for classification of `new' infants as having asphyxia or not. Our prototype has been tested in a laboratory setting to give prediction accuracy of up to 88.85%. If higher accuracies can be obtained, this research may be a key contributor to the 4th Millennium Development Goal (MDG) of reducing mortality in under-five children.


title: Image Charges Re-Imagined

abstract: We discuss the grounded, equipotential ellipse in two-dimensional electrostatics to illustrate different ways of extending the domain of the potential and placing image charges such that homogeneous boundary conditions are satisfied. In particular, we compare and contrast the Kelvin and Sommerfeld image methods.

categories: physics.class-ph hep-th math-ph math.MP physics.hist-ph


title: Surface structures of ZrO$_2$ films on Rh(111): From two layers to bulk termination

abstract: We have studied zirconia films on a Rh(111) substrate with thicknesses in the range of 2-10 monolayers (ML) using scanning tunneling microscopy (STM) and low-energy electron diffraction (LEED). Zirconia was deposited using a UHV-compatible sputter source, resulting in layer-by-layer growth and good uniformity of the films. For thicknesses of 2-4 ML, a layer-dependent influence of the substrate on the structure of the thin films is observed. Beyond this thickness, films show a (2 $\times$ 1) or a distorted (2 $\times$ 2) surface structure with respect to cubic ZrO$_2$(111); these structures correspond to tetragonal and monoclinic zirconia, respectively. The tetragonal phase occurs for annealing temperatures of up to 730 {\deg}C; transformation to the thermodynamically stable monoclinic phase occurs after annealing at 850 {\deg}C or above. High-temperature annealing also breaks up the films and exposes the Rh(111) substrate. We argue that the tetragonal films are stabilized by oxygen deficiency, while the monoclinic films are only weakly defective and show band bending at defects and grain boundaries. This observation is in agreement with positive charge being responsible for the grain-boundary blocking effect in zirconia-based solid electrolytes. Our work introduces the tetragonal and monoclinic 5 ML-thick ZrO$_2$ films on Rh(111) as well-suited model system for surface-science studies on ZrO$_2$ as they do not exhibit the charging problems of thicker films or the bulk material and show better homogeneity and stability than the previously-studied ZrO$_2$/Pt(111) system.


title: Detailed Balance and Exact Results for Density Fluctuations in Supersonic Turbulence

abstract: The probabilistic approach to turbulence is applied to investigate density fluctuations in supersonic turbulence. We derive kinetic equations for the probability distribution function (PDF) of the logarithm of the density field, $s$, in compressible turbulence in two forms: a first-order partial differential equation involving the average divergence conditioned on the flow density, $\langle \nabla \cdot {\bs u} | s\rangle$, and a Fokker-Planck equation with the drift and diffusion coefficients equal to $-\langle {\bs u} \cdot \nabla s | s\rangle$ and $\langle {\bs u} \cdot \nabla s | s\rangle$, respectively. Assuming statistical homogeneity only, the detailed balance at steady state leads to two exact results, $\langle \nabla \cdot {\bs u} | s \rangle =0$, and $\langle {\bs u} \cdot \nabla s | s\rangle=0$. The former indicates a balance of the flow divergence over all expanding and contracting regions at each given density. The exact results provide an objective criterion to judge the accuracy of numerical codes with respect to the density statistics in supersonic turbulence. We also present a method to estimate the effective numerical diffusion as a function of the flow density and discuss its effects on the shape of the density PDF.


title: Thermally Oxidized Two-dimensional TaS2 as a High-\k{appa} Gate Dielectric for MoS2 Field-Effect Transistors

abstract: We report a new approach to integrating high-\k{appa} dielectrics in both bottom- and top-gated MoS2 field-effect transistors (FETs) through thermal oxidation and mechanical assembly of layered twodimensional (2D) TaS2. Combined X-ray photoelectron spectroscopy (XPS), optical microscopy, atomic force microscopy (AFM), and capacitance-voltage (C-V) measurements confirm that multilayer TaS2 flakes can be uniformly transformed to Ta2O5 with a high dielectric constant of ~ 15.5 via thermal oxidation, while preserving the geometry and ultra-smooth surfaces of 2D TMDs. Top-gated MoS2 FETs fabricated using the thermally oxidized Ta2O5 as gate dielectric demonstrate a high current on/off ratio approaching 106, a subthreshold swing (SS) down to 61 mV/dec, and a field-effect mobility exceeding 60 cm2V-1 s-1 at room temperature, indicating high dielectric quality and low interface trap density.


title: GlymphVIS: Visualizing Glymphatic Transport Pathways Using Regularized Optimal Transport

abstract: The glymphatic system (GS) is a transit passage that facilitates brain metabolic waste removal and its dysfunction has been associated with neurodegenerative diseases such as Alzheimer's disease. The GS has been studied by acquiring temporal contrast enhanced magnetic resonance imaging (MRI) sequences of a rodent brain, and tracking the cerebrospinal fluid injected contrast agent as it flows through the GS. We present here a novel visualization framework, GlymphVIS, which uses regularized optimal transport (OT) to study the flow behavior between time points at which the images are taken. Using this regularized OT approach, we can incorporate diffusion, handle noise, and accurately capture and visualize the time varying dynamics in GS transport. Moreover, we are able to reduce the registration mean-squared and infinity-norm error across time points by up to a factor of 5 as compared to the current state-of-the-art method. Our visualization pipeline yields flow patterns that align well with experts' current findings of the glymphatic system.


title: Bound on the rate of entropy change in open quantum systems

abstract: We study the temporal rate of variations of the von Neumann entropy in an open quantum system which interacts with a bath. We show that for almost all initial states of the bath and the system, the time-average of the rate of entropy change is bounded by a function which depends on various properties of the system and environment, and is mostly relatively small. This result holds true under fairly general conditions in almost any arbitrary quantum system.


title: Can a composite heart rate variability biomarker shed new insights about autism spectrum disorder in school-aged children?

abstract: High-frequency heart rate variability (HRV) has identified parasympathetic nervous system alterations in autism spectrum disorder (ASD). In a cohort of school-aged children with and without ASD, we test a set of alternative linear and nonlinear HRV measures, including phase rectified signal averaging, applied to a segment of resting ECG, for associations with ASD vs. other psychiatric conditions. Using machine learning, we identify HRV measures derived from time, frequency, and geometric signal-analytical domains that (1) identify children with ASD relative to peers with receiver operating curve area of .89, and (2) differentiate such children from those with conduct problems or depression. Despite the small cohort and lack of prospective external validation, these preliminary results warrant larger prospective validation studies.


title: Reproducible and Interpretable Spiculation Quantification for Lung Cancer Screening

abstract: Spiculations are important predictors of lung cancer malignancy, which are spikes on the surface of the pulmonary nodules. In this study, we proposed an interpretable and parameter-free technique to quantify the spiculation using area distortion metric obtained by the conformal (angle-preserving) spherical parameterization. We exploit the insight that for an angle-preserved spherical mapping of a given nodule, the corresponding negative area distortion precisely characterizes the spiculations on that nodule. We introduced novel spiculation scores based on the area distortion metric and spiculation measures. We also semi-automatically segment lung nodule (for reproducibility) as well as vessel and wall attachment to differentiate the real spiculations from lobulation and attachment. A simple pathological malignancy prediction model is also introduced. We used the publicly-available LIDC-IDRI dataset pathologists (strong-label) and radiologists (weak-label) ratings to train and test radiomics models containing this feature, and then externally validate the models. We achieved AUC$=$0.80 and 0.76, respectively, with the models trained on the 811 weakly-labeled LIDC datasets and tested on the 72 strongly-labeled LIDC and 73 LUNGx datasets; the previous best model for LUNGx had AUC$=$0.68. The number-of-spiculations feature was found to be highly correlated (Spearman's rank correlation coefficient $\rho = 0.44$) with the radiologists' spiculation score. We developed a reproducible and interpretable, parameter-free technique for quantifying spiculations on nodules. The spiculation quantification measures was then applied to the radiomics framework for pathological malignancy prediction with reproducible semi-automatic segmentation of nodule. Using our interpretable features (size, attachment, spiculation, lobulation), we were able to achieve higher performance than previous models.


title: ParaNet - Using Dense Blocks for Early Inference

abstract: DenseNets have been shown to be a competitive model among recent convolutional network architectures. These networks utilize Dense Blocks, which are groups of densely connected layers where the output of a hidden layer is fed in as the input of every other layer following it. In this paper, we aim to improve certain aspects of DenseNet, especially when it comes to practicality. We introduce ParaNet, a new architecture that constructs three pipelines which allow for early inference. We additionally introduce a cascading mechanism such that different pipelines are able to share parameters, as well as logit matching between the outputs of the pipelines. We separately evaluate each of the newly introduced mechanisms of ParaNet, then evaluate our proposed architecture on CIFAR-100.


title: Trajectory Tracking Control of a Flexible Spine Robot, With and Without a Reference Input

abstract: The Underactuated Lightweight Tensegrity Robotic Assistive Spine (ULTRA Spine) project is an ongoing effort to develop a flexible, actuated backbone for quadruped robots. In this work, model-predictive control is used to track a trajectory in the robot's state space, in simulation. The state trajectory used here corresponds to a bending motion of the spine, with translations and rotations of the moving vertebrae. Two different controllers are presented in this work: one that does not use a reference input but includes smoothing constrants, and a second one that uses a reference input without smoothing. For the smoothing controller, without reference inputs, the error converges to zero, while the simpler-to-tune controller with an input reference shows small errors but not complete convergence. It is expected that this controller will converge as it is improved further.


title: Continuous Generation of Volumetric Images During Stereotactic Body Radiation Therapy Using Periodic kV Imaging and an External Respiratory Surrogate

abstract: We present a technique for continuous generation of volumetric images during SBRT using periodic kV imaging and an external respiratory surrogate signal to drive a patient-specific PCA motion model. Using the on-board imager, kV radiographs are acquired every 3 seconds and used to fit the parameters of a motion model so that it matches observed changes in internal patient anatomy. A multi-dimensional correlation model is established between the motion model parameters and the external surrogate position and velocity, enabling volumetric image reconstruction between kV imaging time points. Performance of the algorithm was evaluated using 10 realistic eXtended CArdiac-Torso (XCAT) digital phantoms including 3D anatomical respiratory deformation programmed with 3D tumor positions measured with orthogonal kV imaging of implanted fiducial gold markers. The clinically measured ground truth 3D tumor positions provided a dataset with realistic breathing irregularities, and the combination of periodic on-board kV imaging with recorded external respiratory surrogate signal was used for correlation modeling to account for any changes in internal-external correlation. The three-dimensional tumor positions are reconstructed with an average root mean square error (RMSE) of 1.47 mm, and an average 95th percentile 3D positional error of 2.80 mm compared with the clinically measured ground truth 3D tumor positions. This technique enables continuous 3D anatomical image generation based on periodic kV imaging of internal anatomy without the additional dose of continuous kV imaging. The 3D anatomical images produced using this method can be used for treatment verification and delivered dose computation in the presence of irregular respiratory motion.


title: Voice Conversion with Conditional SampleRNN

abstract: Here we present a novel approach to conditioning the SampleRNN generative model for voice conversion (VC). Conventional methods for VC modify the perceived speaker identity by converting between source and target acoustic features. Our approach focuses on preserving voice content and depends on the generative network to learn voice style. We first train a multi-speaker SampleRNN model conditioned on linguistic features, pitch contour, and speaker identity using a multi-speaker speech corpus. Voice-converted speech is generated using linguistic features and pitch contour extracted from the source speaker, and the target speaker identity. We demonstrate that our system is capable of many-to-many voice conversion without requiring parallel data, enabling broad applications. Subjective evaluation demonstrates that our approach outperforms conventional VC methods.


title: Quantification of Local Metabolic Tumor Volume Changes by Registering Blended PET-CT Images for Prediction of Pathologic Tumor Response

abstract: Quantification of local metabolic tumor volume (MTV) chan-ges after Chemo-radiotherapy would allow accurate tumor response evaluation. Currently, local MTV changes in esophageal (soft-tissue) cancer are measured by registering follow-up PET to baseline PET using the same transformation obtained by deformable registration of follow-up CT to baseline CT. Such approach is suboptimal because PET and CT capture fundamentally different properties (metabolic vs. anatomy) of a tumor. In this work we combined PET and CT images into a single blended PET-CT image and registered follow-up blended PET-CT image to baseline blended PET-CT image. B-spline regularized diffeomorphic registration was used to characterize the large MTV shrinkage. Jacobian of the resulting transformation was computed to measure the local MTV changes. Radiomic features (intensity and texture) were then extracted from the Jacobian map to predict pathologic tumor response. Local MTV changes calculated using blended PET-CT registration achieved the highest correlation with ground truth segmentation (R=0.88) compared to PET-PET (R=0.80) and CT-CT (R=0.67) registrations. Moreover, using blended PET-CT registration, the multivariate prediction model achieved the highest accuracy with only one Jacobian co-occurrence texture feature (accuracy=82.3%). This novel framework can replace the conventional approach that applies CT-CT transformation to the PET data for longitudinal evaluation of tumor response.


title: Integrity Authentication for SQL Query Evaluation on Outsourced Databases: A Survey

abstract: Spurred by the development of cloud computing, there has been considerable recent interest in the Database-as-a-Service (DaaS) paradigm. Users lacking in expertise or computational resources can outsource their data and database management needs to a third-party service provider. Outsourcing, however, raises an important issue of result integrity: how can the client verify with lightweight overhead that the query results returned by the service provider are correct (i.e., the same as the results of query execution locally)? This survey focuses on categorizing and reviewing the progress on the current approaches for result integrity of SQL query evaluation in the DaaS model. The survey also includes some potential future research directions for result integrity verification of the outsourced computations.


title: Shape of the zeroth Landau level in graphene with non-diagonal disorder

abstract: Non-diagonal (bond) disorder in graphene broadens Landau levels (LLs) in the same way as random potential. The exception is the zeroth LL, $n=0$, which is robust to the bond disorder, since it does not mix different $n=0$ states within a given valley. The mechanism of broadening of the $n=0$ LL is the inter-valley scattering. Several numerical simulations of graphene with bond disorder had established that $n=0$ LL is not only anomalously narrow but also that its shape is very peculiar with three maxima, one at zero energy, $E=0$, and two others at finite energies $\pm E$. We study theoretically the structure of the states in $n=0$ LL in the presence of bond disorder. Adopting the assumption that the bond disorder is strongly anisotropic, namely, that one type of bonds is perturbed much stronger than other two, allowed us to get an analytic expression for the density of states which agrees with numerical simulations remarkably well. On the qualitative level, our key finding is that delocalization of $E=0$ state has a dramatic back effect on the density of states near $E=0$. The origin of this unusual behavior is the strong correlation of eigenstates in different valleys.


title: A Deterministic Self-Organizing Map Approach and its Application on Satellite Data based Cloud Type Classification

abstract: A self-organizing map (SOM) is a type of competitive artificial neural network, which projects the high-dimensional input space of the training samples into a low-dimensional space with the topology relations preserved. This makes SOMs supportive of organizing and visualizing complex data sets and have been pervasively used among numerous disciplines with different applications. Notwithstanding its wide applications, the self-organizing map is perplexed by its inherent randomness, which produces dissimilar SOM patterns even when being trained on identical training samples with the same parameters every time, and thus causes usability concerns for other domain practitioners and precludes more potential users from exploring SOM based applications in a broader spectrum. Motivated by this practical concern, we propose a deterministic approach as a supplement to the standard self-organizing map. In accordance with the theoretical design, the experimental results with satellite cloud data demonstrate the effective and efficient organization as well as simplification capabilities of the proposed approach.


title: A Trio Neural Model for Dynamic Entity Relatedness Ranking

abstract: Measuring entity relatedness is a fundamental task for many natural language processing and information retrieval applications. Prior work often studies entity relatedness in static settings and an unsupervised manner. However, entities in real-world are often involved in many different relationships, consequently entity-relations are very dynamic over time. In this work, we propose a neural networkbased approach for dynamic entity relatedness, leveraging the collective attention as supervision. Our model is capable of learning rich and different entity representations in a joint framework. Through extensive experiments on large-scale datasets, we demonstrate that our method achieves better results than competitive baselines.


title: To Cluster, or Not to Cluster: An Analysis of Clusterability Methods

abstract: Clustering is an essential data mining tool that aims to discover inherent cluster structure in data. For most applications, applying clustering is only appropriate when cluster structure is present. As such, the study of clusterability, which evaluates whether data possesses such structure, is an integral part of cluster analysis. However, methods for evaluating clusterability vary radically, making it challenging to select a suitable measure. In this paper, we perform an extensive comparison of measures of clusterability and provide guidelines that clustering users can reference to select suitable measures for their applications.


title: A Molecular Dynamics Study of The Translation and Rotation of Amphiphilic Janus Nanoparticles at a Vapor-Liquid Surface

abstract: We study the effects of heterogeneity on interfacial pinning and hydrodynamic drag using molecular dynamics (MD) simulations of Janus nanospheres at a liquid/vapor interface. We construct the free energy landscape for this system, both in the continuum approximation using surfaces tensions and the flat-interface approximation and atomistically using MD and thermodynamic integration. The results of the two methods differ in detail due to interfacial distortion and finite width, as well as thermal fluctuations, and only the MD landscape is consistent with simulations of a nanosphere approaching the interface from the liquid or vapor side. When dragged along an interface, these Janus particles exhibit a velocity-dependent tilt accompanied by a weak variation in drag force, but never an enhancement of the drag force beyond the value when fully immersed. This velocity dependence arises when the interface is pinned at heterogeneities and prevents the particle from rotating, and similar behavior is observed for homogeneous but non-spherical particles. The occurrence of different particle orientations having different drag coefficients may lead to an apparent violation of the Stokes-Einstein relation.


title: BOP: Benchmark for 6D Object Pose Estimation

abstract: We propose a benchmark for 6D pose estimation of a rigid object from a single RGB-D input image. The training data consists of a texture-mapped 3D object model or images of the object in known 6D poses. The benchmark comprises of: i) eight datasets in a unified format that cover different practical scenarios, including two new datasets focusing on varying lighting conditions, ii) an evaluation methodology with a pose-error function that deals with pose ambiguities, iii) a comprehensive evaluation of 15 diverse recent methods that captures the status quo of the field, and iv) an online evaluation system that is open for continuous submission of new results. The evaluation shows that methods based on point-pair features currently perform best, outperforming template matching methods, learning-based methods and methods based on 3D local features. The project website is available at bop.felk.cvut.cz.


title: $L_p$ and almost sure convergence of estimation on heavy tail index under random censoring

abstract: In this paper, we prove $L_p,\ p\geq 2$ and almost sure convergence of tail index estimator mentioned in \cite{grama2008} under random censoring and several assumptions. $p$th moment of the error of the estimator is proved to be of order $O\left(\frac{1}{\log^{m\kappa/2}n}\right)$ with given assumptions. We also perform several finite sample simulations to quantify performance of this estimator. Finite sample results show that the proposed estimator is effective in finding underlying tail index even when censor rate is high.


title: Aperiodic Array Synthesis for Multi-User MIMO Applications

abstract: This paper demonstrates the advantages of aperiodic arrays in multi-user multiple-input multiple-output systems for future mobile communication applications. We propose a novel aperiodic array synthesis method which account for the statistics of the propagation channel and the adaptive beamforming algorithm. Clear performance gains in line-of-sight dominated propagation environments are achieved in terms of the signal-to-interference-plus-noise ratio, the sum rate capacity, as well as the spread of the amplifier output power as compared to their regular counterparts. We also show that the performance is not sacrificed in rich scattering environments. Hence, aperiodic array layouts can provide performance gains in millimeter-wave applications with a dominating line-of-sight component.


title: Dynamics of a spatially developing liquid jet with slower coaxial gas flow

abstract: A three-dimensional round liquid jet within a low-speed coaxial gas flow is numerically simulated and explained via vortex dynamics ($\lambda_2$ method). The instabilities on the liquid-gas interface reflect well the vortex interactions around the interface. Certain key features are identified for the first time. Two types of surface deformations are distinguished, which are separated by a large indentation on the jet stem: First, those near the jet start-up cap are encapsulated inside the recirculation zone behind the cap. These deformations are directly related to the dynamics of the growing cap and well explained by the vortices generated there. Second, deformations occurring farther upstream of the cap are mainly driven by the Kelvin-Helmholtz (KH) instability at the interface. Three-dimensional deformations occur in the vortex structures first, and the initially axisymmetric KH vortices deform and lead to several liquid lobes, which stretch first as thinning sheets and then either continue stretching directly into elongated ligaments - at lower relative velocity - or perforate to create liquid bridges and holes - at higher relative velocity. The different scenarios depend on Weber and Reynolds numbers based on the relative gas-liquid velocity as was found in the temporal studies. The deformations in the upstream region are well portrayed in a frame moving with the convective velocity of the liquid jet. The usefulness of the temporal analyses are now established.


title: Groupoids and Relative Internality

abstract: In a stable theory, a stationary type $q \in S(A)$ internal to a family of partial types $\mathcal{P}$ over $A$ gives rise to a type-definable group, called its binding group. This group is isomorphic to the group $\mathrm{Aut}(q/\mathcal{P},A)$ of permutations of the set of realizations of $q$, induced by automorphisms of the monster model, fixing $\mathcal{P} \cup A$ pointwise. In this paper, we investigate families of internal types varying uniformly, what we will call relative internality. We prove that the binding groups also vary uniformly, and are the isotropy groups of a natural type-definable groupoid (and even more). We then investigate how properties of this groupoid are related to properties of the type. In particular, we obtain internality criteria for certain 2-analysable types, and a sufficient condition for a type to preserve internality.


title: Phase factors of periodically driven two-level systems

abstract: Using a perturbative solution for a periodically driven two-level quantum system, we show how to obtain phase factors for both a two-level quantum system and two two-level quantum systems non-interacting and interacting. The method is easily implemented by numerical routines and presents the advantage of being stable for long-time periods. We furthermore explore the possibility of implementing a quantum phase gate using the perturbative solution.


title: Rate-Splitting for Multi-Antenna Non-Orthogonal Unicast and Multicast Transmission: Spectral and Energy Efficiency Analysis

abstract: In a Non-Orthogonal Unicast and Multicast (NOUM) transmission system, a multicast stream intended to all the receivers is superimposed in the power domain on the unicast streams. One layer of Successive Interference Cancellation (SIC) is required at each receiver to remove the multicast stream before decoding its intended unicast stream. In this paper, we first show that a linearly-precoded 1-layer Rate-Splitting (RS) strategy at the transmitter can efficiently exploit this existing SIC receiver architecture. We further propose multi-layer transmission strategies based on the generalized RS and power-domain Non-Orthogonal Multiple Access (NOMA). Two different objectives are studied for the design of the precoders, namely, maximizing the Weighted Sum Rate (WSR) of the unicast messages and maximizing the system Energy Efficiency (EE), both subject to Quality of Service (QoS) rate requirements of all the messages and a sum power constraint. A Weighted Minimum Mean Square Error (WMMSE)-based algorithm and a Successive Convex Approximation (SCA)-based algorithm are proposed to solve the WSR and EE problems, respectively. Numerical results show that the proposed RS-assisted NOUM transmission strategies are more spectrally and energy efficient than the conventional Multi-User Linear-Precoding (MU-LP), Orthogonal Multiple Access (OMA) and power-domain NOMA in a wide range of user deployments (with a diversity of channel directions, channel strengths and qualities of channel state information at the transmitter) and network loads (underloaded and overloaded regimes). It is superior for the downlink multi-antenna NOUM transmission.


title: A Bayesian Approach to Restricted Latent Class Models for Scientifically-Structured Clustering of Multivariate Binary Outcomes

abstract: In this paper, we propose a general framework for combining evidence of varying quality to estimate underlying binary latent variables in the presence of restrictions imposed to respect the scientific context. The resulting algorithms cluster the multivariate binary data in a manner partly guided by prior knowledge. The primary model assumptions are that 1) subjects belong to classes defined by unobserved binary states, such as the true presence or absence of pathogens in epidemiology, or of antibodies in medicine, or the "ability" to correctly answer test questions in psychology, 2) a binary design matrix $\Gamma$ specifies relevant features in each class, and 3) measurements are independent given the latent class but can have different error rates. Conditions ensuring parameter identifiability from the likelihood function are discussed and inform the design of a novel posterior inference algorithm that simultaneously estimates the number of clusters, design matrix $\Gamma$, and model parameters. In finite samples and dimensions, we propose prior assumptions so that the posterior distribution of the number of clusters and the patterns of latent states tend to concentrate on smaller values and sparser patterns, respectively. The model readily extends to studies where some subjects' latent classes are known or important prior knowledge about differential measurement accuracy is available from external sources. The methods are illustrated with an analysis of protein data to detect clusters representing auto-antibody classes among scleroderma patients.


title: Big bounce and closed universe from spin and torsion

abstract: We analyze the dynamics of a homogeneous and isotropic universe in the Einstein$--$Cartan theory of gravity. The coupling between the spin and torsion prevents gravitational singularities and replaces the Big Bang with a nonsingular big bounce, at which the universe transitions from contraction to expansion. We show that a closed universe exists only when the product of the scale factor and temperature is higher than a particular threshold, contrary to a flat universe and an open universe, which are not restricted. During inflation, this product must increase to another threshold, so that the universe can reach dark-energy acceleration.


title: Composable block solvers for the four-field double porosity/permeability model

abstract: The objective of this paper is twofold. First, we propose two composable block solver methodologies to solve the discrete systems that arise from finite element discretizations of the double porosity/permeability (DPP) model. The DPP model, which is a four-field mathematical model, describes the flow of a single-phase incompressible fluid in a porous medium with two distinct pore-networks and with a possibility of mass transfer between them. Using the composable solvers feature available in PETSc and the finite element libraries available under the Firedrake Project, we illustrate two different ways by which one can effectively precondition these large systems of equations. Second, we employ the recently developed performance model called the Time-Accuracy-Size (TAS) spectrum to demonstrate that the proposed composable block solvers are scalable in both the parallel and algorithmic sense. Moreover, we utilize this spectrum analysis to compare the performance of three different finite element discretizations (classical mixed formulation with H(div) elements, stabilized continuous Galerkin mixed formulation, and stabilized discontinuous Galerkin mixed formulation) for the DPP model. Our performance spectrum analysis demonstrates that the composable block solvers are fine choices for any of these three finite element discretizations. Sample computer codes are provided to illustrate how one can easily implement the proposed block solver methodologies through PETSc command line options.


title: When You Should Use Lists in Haskell (Mostly, You Should Not)

abstract: We comment on the over-use of lists in functional programming. With this respect, we review history of Haskell and some of its libraries, and hint at current developments.


title: Code Generation for Higher Inductive Types

abstract: Higher inductive types are inductive types that include nontrivial higher-dimensional structure, represented as identifications that are not reflexivity. While work proceeds on type theories with a computational interpretation of univalence and higher inductive types, it is convenient to encode these structures in more traditional type theories with mature implementations. However, these encodings involve a great deal of error-prone additional syntax. We present a library that uses Agda's metaprogramming facilities to automate this process, allowing higher inductive types to be specified with minimal additional syntax.


title: Theory-Driven Automated Content Analysis of Suicidal Tweets : Using Typicality-Based Classification for LDA Dataset

abstract: This study provides a methodological framework for the computer to classify tweets according to variables of the Theory of Planned Behavior. We present a sequential process of automated text analysis which combined supervised approach and unsupervised approach in order to make the computer to detect one of TPB variables in each tweet. We conducted Latent Dirichlet Allocation (LDA), Nearest Neighbor, and then assessed "typicality" of newly labeled tweets in order to predict classification boundary. Furthermore, this study reports findings from a content analysis of suicide-related tweets which identify traits of information environment in Twitter. Consistent with extant literature about suicide coverage, the findings demonstrate that tweets often contain information which prompt perceived behavior control of committing suicide, while rarely provided deterring information on suicide. We conclude by highlighting implications for methodological advances and empirical theory studies.


title: Reducing model bias in a deep learning classifier using domain adversarial neural networks in the MINERvA experiment

abstract: We present a simulation-based study using deep convolutional neural networks (DCNNs) to identify neutrino interaction vertices in the MINERvA passive targets region, and illustrate the application of domain adversarial neural networks (DANNs) in this context. DANNs are designed to be trained in one domain (simulated data) but tested in a second domain (physics data) and utilize unlabeled data from the second domain so that during training only features which are unable to discriminate between the domains are promoted. MINERvA is a neutrino-nucleus scattering experiment using the NuMI beamline at Fermilab. $A$-dependent cross sections are an important part of the physics program, and these measurements require vertex finding in complicated events. To illustrate the impact of the DANN we used a modified set of simulation in place of physics data during the training of the DANN and then used the label of the modified simulation during the evaluation of the DANN. We find that deep learning based methods offer significant advantages over our prior track-based reconstruction for the task of vertex finding, and that DANNs are able to improve the performance of deep networks by leveraging available unlabeled data and by mitigating network performance degradation rooted in biases in the physics models used for training.


title: $p$-adic Mellin Amplitudes

abstract: In this paper, we propose a $p$-adic analog of Mellin amplitudes for scalar operators, and present the computation of the general contact amplitude as well as arbitrary-point tree-level amplitudes for bulk diagrams involving up to three internal lines, and along the way obtain the $p$-adic version of the split representation formula. These amplitudes share noteworthy similarities with the usual (real) Mellin amplitudes for scalars, but are also significantly simpler, admitting closed-form expressions where none are available over the reals. The dramatic simplicity can be attributed to the absence of descendant fields in the $p$-adic formulation.


title: Directed flow in an extended multiphase transport model

abstract: We have studied the rapidity-odd directed flow in $^{197}$Au+$^{197}$Au collisions in the beam energy range from $\sqrt{s_{NN}}$ = 7.7 to 39 GeV within the framework of an extended multiphase transport model with both partonic and hadronic mean-field potentials incorporated. Effects of the partonic scatterings, mean-field potentials, hadronization, and hadronic evolution on the directed flow are investigated, and it is found that the final directed flow is mostly sensitive to the partonic scatterings and the hadronization mechanism. Our study shows that a negative slope of the proton directed flow does not necessarily need the equation of state with a first-order phase transition.


title: Simple Proofs for the Derivative Estimates of the Holomorphic Motion near Two Boundary Points of the Mandelbrot Set

abstract: For the complex quadratic family $q_c:z\mapsto z^2+c$, it is known that every point in the Julia set $J(q_c)$ moves holomorphically on $c$ except at the boundary points of the Mandelbrot set. In this note, we present short proofs of the following derivative estimates of the motions near the boundary points $1/4$ and $-2$: for each $z = z(c)$ in the Julia set, the derivative $dz(c)/dc$ is uniformly $O(1/\sqrt{1/4-c})$ when real $c\nearrow 1/4$; and is uniformly $O(1/\sqrt{-2-c})$ when real $c\nearrow -2$. These estimates of the derivative imply Hausdorff convergence of the Julia set $J(q_c)$ when $c$ approaches these boundary points. In particular, the Hausdorff distance between $J(q_c)$ with $0\le c<1/4$ and $J(q_{1/4})$ is exactly $\sqrt{1/4-c}$.


title: Backward Stochastic Riccati Equation with Jumps associated with Stochastic Linear Quadratic Optimal Control with Jumps and Random Coefficients

abstract: In this paper, we investigate the solvability of matrix valued Backward stochastic Riccati equations with jumps (BSREJ), which is associated with a stochastic linear quadratic (SLQ) optimal control problem with random coefficients and driven by both Brownian motion and Poisson jumps. By dynamic programming principle, Doob-Meyer decomposition and inverse flow technique, the existence and uniqueness of the solution for the BSREJ is established. The difficulties addressed to this issue not only are brought from the high nonlinearity of the generator of the BSREJ like the case driven only by Brownian motion, but also from that i) the inverse flow of the controlled linear stochastic differential equation driven by Poisson jumps may not exist without additional technical condition, and ii) how to show the inverse matrix term involving jump process in the generator is well-defined. Utilizing the structure of the optimal problem, we overcome these difficulties and establish the existence of the solution. In additional, a verification theorem for BSREJ is given which implies the uniqueness of the solution.


title: Multi-fold contour integrals of certain ratios of Euler gamma functions from Feynman diagrams: orthogonality of triangles

abstract: We observe a property of orthogonality of the Mellin-Barnes transformation of the triangle one-loop diagrams, which follows from our previous papers [JHEP {\bf 0808} (2008) 106, JHEP {\bf 1003} (2010) 051, JMP {\bf 51} (2010) 052304]. In those papers it has been established that Usyukina-Davydychev functions are invariant with respect to Fourier transformation. This has been proved at the level of graphs and also via the Mellin-Barnes transformation. We partially apply to one-loop massless scalar diagram the same trick in which the Mellin-Barnes transformation was involved and obtain the property of orthogonality of the corresponding MB transforms under integration over contours in two complex planes with certain weight. This property is valid in an arbitrary number of dimensions.


title: On the motive of Kapustka-Rampazzo's Calabi-Yau threefolds

abstract: Kapustka and Rampazzo have exhibited pairs of Calabi-Yau threefolds $X$ and $Y$ that are L-equivalent and derived equivalent, without being birational. We complete the picture by showing that $X$ and $Y$ have isomorphic Chow motives.


title: On the motive of intersections of two Grassmannians in ${\mathbb{P}}^9$

abstract: Using intersections of two Grassmannians in ${\mathbb{P}}^9$, Ottem-Rennemo and Borisov-C\u{a}ld\u{a}raru-Perry have exhibited pairs of Calabi-Yau threefolds $X$ and $Y$ that are deformation equivalent, L-equivalent and derived equivalent, but not birational. To complete the picture, we show that $X$ and $Y$ have isomorphic Chow motives.


title: Invariant Sets in Quasiperiodically Forced Dynamical Systems

abstract: This paper addresses structures of state space in quasiperiodically forced dynamical systems. We develop a theory of ergodic partition of state space in a class of measure-preserving and dissipative flows, which is a natural extension of the existing theory for measure-preserving maps. The ergodic partition result is based on eigenspace at eigenvalue 0 of the associated Koopman operator, which is realized via time-averages of observables, and provides a constructive way to visualize a low-dimensional slice through a high-dimensional invariant set. We apply the result to the systems with a finite number of attractors and show that the time-average of a continuous observable is well-defined and reveals the invariant sets, namely, a finite number of basins of attraction. We provide a characterization of invariant sets in the quasiperiodically forced systems. A theoretical result on uniform boundedness of the invariant sets is presented. The series of theoretical results enables numerical analysis of invariant sets in the quasiperiodically forced systems based on the ergodic partition and time-averages. Using this, we analyze a nonlinear model of complex power grids that represents the short-term swing instability, named the coherent swing instability. We show that our theoretical results can be used to understand stability regions in such complex systems.


title: Predictive Model of Hydrogen Trapping and Bubbling in Nanovoids in BCC Metals

abstract: Interplay between hydrogen and nanovoids, despite long-recognized as a central aspect in hydrogen-induced damages in structural materials, remains poorly understood. Focusing on tungsten as a model BCC system, the present study, for the first time, explicitly demonstrated sequential adsorption of hydrogen adatoms on Wigner-Seitz squares of nanovoids with distinct energy levels. Interaction between hydrogen adatoms on the nanovoid surface is shown to be dominated by pairwise power law repulsion. A predictive model was established for quantitative prediction of configurations and energetics of hydrogen adatoms in nanovoids. This model, further combined with equation of states of hydrogen gas, enables prediction of hydrogen molecule formation in nanovoids. Multiscale simulations based on the predictive model were performed, showing excellent agreement with experiments. This work clarifies fundamental physics and provides full-scale predictive model for hydrogen trapping and bubbling in nanovoids, offering long-sought mechanistic insights crucial for understanding hydrogen-induced damages in structural materials.


title: Some inequalities for interpolational operator means

abstract: Using the properties of geometric mean, we shall show for any $0\le \alpha ,\beta \le 1$, \[f\left( A{{\nabla }_{\alpha }}B \right)\le f\left( \left( A{{\nabla }_{\alpha }}B \right){{\nabla }_{\beta }}A \right){{\sharp}_{\alpha }}f\left( \left( A{{\nabla }_{\alpha }}B \right){{\nabla }_{\beta }}B \right)\le f\left( A \right){{\sharp}_{\alpha }}f\left( B \right)\] whenever $f$ is a non-negative operator log-convex function, $A,B\in \mathcal{B}\left( \mathcal{H} \right)$ are positive operators, and $0\le \alpha ,\beta \le 1$. As an application of this operator mean inequality, we present several refinements of the Aujla subadditive inequality for operator monotone decreasing functions. Also, in a similar way, we consider some inequalities of Ando's type. Among other things, it is shown that if $\Phi $ is a positive linear map, then \[\Phi \left( A{{\sharp}_{\alpha }}B \right)\le \Phi \left( \left( A{{\sharp}_{\alpha }}B \right){{\sharp}_{\beta }}A \right){{\sharp}_{\alpha }}\Phi \left( \left( A{{\sharp}_{\alpha }}B \right){{\sharp}_{\beta }}B \right)\le \Phi \left( A \right){{\sharp}_{\alpha }}\Phi \left( B \right).\]


title: Spontaneous symmetry breaking induced by quantum monitoring

abstract: Spontaneous symmetry breaking (SSB) is responsible for structure formation in scenarios ranging from condensed matter to cosmology. SSB is broadly understood in terms of perturbations to the Hamiltonian governing the dynamics or to the state of the system. We study SSB due to quantum monitoring of a system via continuous quantum measurements. The acquisition of information during the measurement process induces a measurement back-action that seeds SSB. In this setting, by monitoring different observables, an observer can tailor the topology of the vacuum manifold, the pattern of symmetry breaking and the nature of the resulting domains and topological defects.


title: Multiobjective Optimization Training of PLDA for Speaker Verification

abstract: Most current state-of-the-art text-independent speaker verification systems take probabilistic linear discriminant analysis (PLDA) as their backend classifiers. The parameters of PLDA are often estimated by maximizing the objective function, which focuses on increasing the value of log-likelihood function, but ignoring the distinction between speakers. In order to better distinguish speakers, we propose a multi-objective optimization training for PLDA. Experiment results show that the proposed method has more than 10% relative performance improvement in both EER and MinDCF on the NIST SRE14 i-vector challenge dataset, and about 20% relative performance improvement in EER on the MCE18 dataset.


title: Heterotic Type I Strings at High Temperature

abstract: We show that the high temperature limits of the heterotic E8xE8 and Spin(32)/Z2 strings and their Type I A/B superstring duals are finite and convergent. The Hagedorn growth of the degeneracies in the string mass level expansion is suppressed by an exponential linear in the mass level number for both heterotic strings, and suppressed by the exponential of the negative square root of the mass level number for the Type IB superstring. However, in the massless gauge field theoretic limit of the Type IB open and closed superstring, we find clear evidence for the thermal deconfinement phase transition at the self-dual temperature by examining the annulus graph alone. Above the self-dual temperature, there is a discontinuity in the first derivative with respect to temperature of both the free energy, and the heavy quark potential, leading to a deconfined thermal gluon ensemble, with universal 1/r potential, and temperature dependent corrections, as predicted by Luscher and Weisz. A number of essential aspects of the worldsheet formalism of the heterotic strings are derived in an appendix, deducing thereby the O8-D0-D8brane Type IA duals of all of the heterotic CHL island universe moduli spaces.


title: A non-linear fluid suspension model for blood flow

abstract: Motivated by the complex rheological behaviors observed in small/micro scale blood vessels, such as the Fahraeus effect, plasma-skimming, shear-thinning, etc., we develop a non-linear suspension model for blood. The viscosity is assumed to depend on the volume fraction (hematocrit) and the shear rate. The migration of the red blood cells (RBCs) is studied using a concentration flux equation. A parametric study with two representative problems, namely simple shear flow and a pressure driven flow demonstrate the ability of this reduced-order model to reproduce several key features of the two-fluid model (mixture theory approach), with much lower computational cost.


title: A Comparison of the Taguchi Method and Evolutionary Optimization in Multivariate Testing

abstract: Multivariate testing has recently emerged as a promising technique in web interface design. In contrast to the standard A/B testing, multivariate approach aims at evaluating a large number of values in a few key variables systematically. The Taguchi method is a practical implementation of this idea, focusing on orthogonal combinations of values. This paper evaluates an alternative method: population-based search, i.e. evolutionary optimization. Its performance is compared to that of the Taguchi method in several simulated conditions, including an orthogonal one designed to favor the Taguchi method, and two realistic conditions with dependences between variables. Evolutionary optimization is found to perform significantly better especially in the realistic conditions, suggesting that it forms a good approach for web interface design in the future.


title: Multi-scale CNN stereo and pattern removal technique for underwater active stereo system

abstract: Demands on capturing dynamic scenes of underwater environments are rapidly growing. Passive stereo is applicable to capture dynamic scenes, however the shape with textureless surfaces or irregular reflections cannot be recovered by the technique. In our system, we add a pattern projector to the stereo camera pair so that artificial textures are augmented on the objects. To use the system at underwater environments, several problems should be compensated, i.e., refraction, disturbance by fluctuation and bubbles. Further, since surface of the objects are interfered by the bubbles, projected patterns, etc., those noises and patterns should be removed from captured images to recover original texture. To solve these problems, we propose three approaches; a depth-dependent calibration, Convolutional Neural Network(CNN)-stereo method and CNN-based texture recovery method. A depth-dependent calibration is our analysis to find the acceptable depth range for approximation by center projection to find the certain target depth for calibration. In terms of CNN stereo, unlike common CNNbased stereo methods which do not consider strong disturbances like refraction or bubbles, we designed a novel CNN architecture for stereo matching using multi-scale information, which is intended to be robust against such disturbances. Finally, we propose a multi-scale method for bubble and a projected-pattern removal method using CNNs to recover original textures. Experimental results are shown to prove the effectiveness of our method compared with the state of the art techniques. Furthermore, reconstruction of a live swimming fish is demonstrated to confirm the feasibility of our techniques.


title: Detection and Mitigation of Attacks on Transportation Networks as a Multi-Stage Security Game

abstract: In recent years, state-of-the-art traffic-control devices have evolved from standalone hardware to networked smart devices. Smart traffic control enables operators to decrease traffic congestion and environmental impact by acquiring real-time traffic data and changing traffic signals from fixed to adaptive schedules. However, these capabilities have inadvertently exposed traffic control to a wide range of cyber-attacks, which adversaries can easily mount through wireless networks or even through the Internet. Indeed, recent studies have found that a large number of traffic signals that are deployed in practice suffer from exploitable vulnerabilities, which adversaries may use to take control of the devices. Thanks to the hardware-based failsafes that most devices employ, adversaries cannot cause traffic accidents directly by setting compromised signals to dangerous configurations. Nonetheless, an adversary could cause disastrous traffic congestion by changing the schedule of compromised traffic signals, thereby effectively crippling the transportation network. To provide theoretical foundations for the protection of transportation networks from these attacks, we introduce a game-theoretic model of launching, detecting, and mitigating attacks that tamper with traffic-signal schedules. We show that finding optimal strategies is a computationally challenging problem, and we propose efficient heuristic algorithms for finding near optimal strategies. We also introduce a Gaussian-process based anomaly detector, which can alert operators to ongoing attacks. Finally, we evaluate our algorithms and the proposed detector using numerical experiments based on the SUMO traffic simulator.


title: Long live the scientists: Tracking the scientific fame of great minds in physics

abstract: This study utilizes global digitalized books and articles to examine the scientific fame of the most influential physicists. Our research reveals that the greatest minds are gone but not forgotten. Their scientific impacts on human history have persisted for centuries. We also find evidence in support of own-group fame preference, i.e., that the scientists have greater reputations in their home countries or among scholars sharing the same languages. We argue that, when applied appropriately, Google Books and Ngram Viewer can serve as promising tools for altmetrics, providing a more comprehensive picture of the impacts scholars and their achievements have made beyond academia.


title: A power-law upper bound on the correlations in the 2D random field Ising model

abstract: As first asserted by Y. Imry and S-K Ma, the famed discontinuity of the magnetization as function of the magnetic field in the two dimensional Ising model is eliminated, for all temperatures, through the addition of quenched random magnetic field of uniform variance, even if that is small. This statement is quantified here by a power-law upper bound on the decay rate of the effect of boundary conditions on the magnetization in finite systems, as function of the distance to the boundary. Unlike exponential decay which is only proven for strong disorder or high temperature, the power-law upper bound is established here for all field strengths and at all temperatures, including zero, for the case of independent Gaussian random field. Our analysis proceeds through a streamlined and quantified version of the Aizenman-Wehr proof of the Imry-Ma rounding effect.


title: Random Matrix Theory Model for Mean Notch Depth of the Diagonally Loaded MVDR Beamformer for a Single Interferer Case

abstract: Adaptive beamformers (ABFs) suppress interferers by placing a notch in the beampattern at the interferer direction. This suppres- sion improves detection of a weaker signals in the presence of strong interferers. Hence the notch depth plays a crucial role in determining the adaptive gain obtained from using ABF over conventional beam- forming. This research derives models for the mean notch depth of a diagonally loaded MVDR ABF for a single interferer case. The model describes the mean notch depth as a function of number of snapshots, the number of sensors in the array, the interferer to noise ratio (INR) level, the interferer direction and the diagonal loading level. The derivation uses random matrix theory results on the be- havior of the eigenvectors of sample covariance matrix. The notch depth predicted by the model is shown to be in close agreement with simulation results over a range of INRs and snapshots.


title: Hyperscaling Internet Graph Analysis with D4M on the MIT SuperCloud

abstract: Detecting anomalous behavior in network traffic is a major challenge due to the volume and velocity of network traffic. For example, a 10 Gigabit Ethernet connection can generate over 50 MB/s of packet headers. For global network providers, this challenge can be amplified by many orders of magnitude. Development of novel computer network traffic analytics requires: high level programming environments, massive amount of packet capture (PCAP) data, and diverse data products for "at scale" algorithm pipeline development. D4M (Dynamic Distributed Dimensional Data Model) combines the power of sparse linear algebra, associative arrays, parallel processing, and distributed databases (such as SciDB and Apache Accumulo) to provide a scalable data and computation system that addresses the big data problems associated with network analytics development. Combining D4M with the MIT SuperCloud manycore processors and parallel storage system enables network analysts to interactively process massive amounts of data in minutes. To demonstrate these capabilities, we have implemented a representative analytics pipeline in D4M and benchmarked it on 96 hours of Gigabit PCAP data with MIT SuperCloud. The entire pipeline from uncompressing the raw files to database ingest was implemented in 135 lines of D4M code and achieved speedups of over 20,000.


title: Intraband divergences in third order optical response of 2D systems

abstract: The existence of large nonlinear optical coefficients is one of the preconditions for using nonlinear optical materials in nonlinear optical devices. For a crystal, such large coefficients can be achieved by matching photon energies with resonant energies between different bands, and so the details of the crystal band structure play an important role. Here we demonstrate that large third-order nonlinearities can also be generally obtained by a different strategy: As any of the incident frequencies or the sum of any two or three frequencies approaches zero, the doped or excited populations of electronic states lead to divergent contributions in the induced current density. We refer to these as intraband divergences, by analogy with the behavior of Drude conductivity in linear response. Physically, such resonant processes can be associated with a combination of inraband and interband optical transitions. Current-induced second order nonlinearity, coherent current injection, and jerk currents are all related to such divergences, and we find similar divergences in degenerate four wave mixing and cross-phase modulation under certain conditions. These divergences are limited by intraband relaxation parameters, and lead to a large optical response from a high quality sample; we find they are very robust with respect to variations in the details of the band structure. To clearly track all of these effects, we analyze gapped graphene, describing the electrons as massive Dirac fermions; under the relaxation time approximation, we derive analytic expressions for the third order conductivities, and identify the divergences that arise in describing the associated nonlinear phenomena.


title: Database-Agnostic Workload Management

abstract: We present a system to support generalized SQL workload analysis and management for multi-tenant and multi-database platforms. Workload analysis applications are becoming more sophisticated to support database administration, model user behavior, audit security, and route queries, but the methods rely on specialized feature engineering, and therefore must be carefully implemented and reimplemented for each SQL dialect, database system, and application. Meanwhile, the size and complexity of workloads are increasing as systems centralize in the cloud. We model workload analysis and management tasks as variations on query labeling, and propose a system design that can support general query labeling routines across multiple applications and database backends. The design relies on the use of learned vector embeddings for SQL queries as a replacement for application-specific syntactic features, reducing custom code and allowing the use of off-the-shelf machine learning algorithms for labeling. The key hypothesis, for which we provide evidence in this paper, is that these learned features can outperform conventional feature engineering on representative machine learning tasks. We present the design of a database-agnostic workload management and analytics service, describe potential applications, and show that separating workload representation from labeling tasks affords new capabilities and can outperform existing solutions for representative tasks, including workload sampling for index recommendation and user labeling for security audits.