title: Collapsing Bandits and Their Application to Public Health Interventions

abstract: We propose and study Collpasing Bandits, a new restless multi-armed bandit (RMAB) setting in which each arm follows a binary-state Markovian process with a special structure: when an arm is played, the state is fully observed, thus "collapsing" any uncertainty, but when an arm is passive, no observation is made, thus allowing uncertainty to evolve. The goal is to keep as many arms in the "good" state as possible by planning a limited budget of actions per round. Such Collapsing Bandits are natural models for many healthcare domains in which workers must simultaneously monitor patients and deliver interventions in a way that maximizes the health of their patient cohort. Our main contributions are as follows: (i) Building on the Whittle index technique for RMABs, we derive conditions under which the Collapsing Bandits problem is indexable. Our derivation hinges on novel conditions that characterize when the optimal policies may take the form of either "forward" or "reverse" threshold policies. (ii) We exploit the optimality of threshold policies to build fast algorithms for computing the Whittle index, including a closed-form. (iii) We evaluate our algorithm on several data distributions including data from a real-world healthcare task in which a worker must monitor and deliver interventions to maximize their patients' adherence to tuberculosis medication. Our algorithm achieves a 3-order-of-magnitude speedup compared to state-of-the-art RMAB techniques while achieving similar performance.


title: Error Estimation and Correction from within Neural Network Differential Equation Solvers

abstract: Neural Network Differential Equation (NN DE) solvers have surged in popularity due to a combination of factors: computational advances making their optimization more tractable, their capacity to handle high dimensional problems, easy interpret-ability of their models, etc. However, almost all NN DE solvers suffer from a fundamental limitation: they are trained using loss functions that depend only implicitly on the error associated with the estimate. As such, validation and error analysis of solution estimates requires knowledge of the true solution. Indeed, if the true solution is unknown, we are often reduced to simply hoping that a "low enough" loss implies "small enough" errors, since explicit relationships between the two are not available/well defined. In this work, we describe a general strategy for efficiently constructing error estimates and corrections for Neural Network Differential Equation solvers. Our methods do not require advance knowledge of the true solutions and obtain explicit relationships between loss functions and the error associated with solution estimates. In turn, these explicit relationships directly allow us to estimate and correct for the errors.


title: Optical and electrical feedback cooling of a silica nanoparticle levitated in a Paul trap

abstract: All three motional modes of a charged dielectric nanoparticle in a Paul trap are cooled by direct feedback to temperatures of a few mK. We test two methods, one based on electrical forces and the other on optical forces; for both methods, we find similar cooling efficiencies. Cooling is characterized for both feedback forces as a function of feedback parameters, background pressure, and the particle's position.


title: Nice guys don't always finish last: succeeding in hierarchical organizations

abstract: What are the chances of an ethical individual rising through the ranks of a political party or a corporation in the presence of unethical peers? To answer this question, I consider a four-player two-stage elimination tournament, in which players are partitioned into those willing to be involved in sabotage behavior and those who are not. I show that, under certain conditions, the latter are more likely to win the tournament.


title: Dynamics and Clogging of colloidal monolayers magnetically driven through an heterogeneous landscape

abstract: We combine experiments and numerical simulations to investigate the emergence of clogging in a system of interacting paramagnetic colloidal particles driven against a disordered landscape of larger obstacles. We consider a single aperture in a landscape of immobile silica particles which are irreversibly attached to the substrate. We use an external rotating magnetic field to generate a traveling wave potential which drives the magnetic particles against these obstacles at a constant and frequency tunable speed. Experimentally we find that the particles display an intermittent dynamics with power law distributions at high frequencies. We reproduce these results by using numerical simulations and shows that clogging in our system arises at large frequency, when the particles desynchronize with the moving landscape. Further, we use the model to explore the hidden role of flexibility in the obstacle displacements and the effect of hydrodynamic interactions between the particles. We also consider numerically the situation of a straight wall and investigate the range of parameters where clogging emerges in such case. Our work provides a soft matter test-bed system to investigate the effect of clogging in driven microscale matter.


title: A nonexistence result for CMC surfaces in hyperbolic 3-manifolds

abstract: We prove that a complete hyperbolic 3-manifold of finite volume does not admit a properly embedded noncompact surface of finite topology with constant mean curvature greater than or equal to 1.


title: Systematic investigation of the effect of 56Ni mixing in the early photospheric velocity evolution of stripped-envelope supernovae

abstract: Mixing of 56Ni, whose nuclear decay energy is a major luminosity source in stripped-envelope supernovae, is known to affect the observational properties of stripped-envelope supernovae such as light-curve and color evolution. Here we systematically investigate the effect of 56Ni mixing on the photospheric velocity evolution in stripped-envelope supernovae. We show that 56Ni mixing significantly affects the early photospheric velocity evolution. The photospheric velocity, which is often used to constrain the ejecta mass and explosion energy, significantly varies by just changing the degree of 56Ni mixing. In addition, the models with a small degree of 56Ni mixing show a flattening in the early photospheric velocity evolution, while the fully mixed models show a monotonic decrease. The velocity flattening appears in both helium and carbon+oxygen progenitor explosions with a variety of ejecta mass, explosion energy, and 56Ni mass. Some stripped-envelope supernovae with early photospheric velocity information do show such a flattening. We find that Type Ib SN 2007Y, which has early photospheric velocity information, has a signature of a moderate degree of 56Ni mixing in the photospheric velocity evolution and about a half of the ejecta is mixed in it. The immediate spectroscopic follow-up observations of stripped-envelope supernovae shortly after the explosion providing the early photospheric evolution give an important clue to constrain 56Ni mixing in the ejecta.


title: Combining Differentiable PDE Solvers and Graph Neural Networks for Fluid Flow Prediction

abstract: Solving large complex partial differential equations (PDEs), such as those that arise in computational fluid dynamics (CFD), is a computationally expensive process. This has motivated the use of deep learning approaches to approximate the PDE solutions, yet the simulation results predicted from these approaches typically do not generalize well to truly novel scenarios. In this work, we develop a hybrid (graph) neural network that combines a traditional graph convolutional network with an embedded differentiable fluid dynamics simulator inside the network itself. By combining an actual CFD simulator (run on a much coarser resolution representation of the problem) with the graph network, we show that we can both generalize well to new situations and benefit from the substantial speedup of neural network CFD predictions, while also substantially outperforming the coarse CFD simulation alone.


title: On the relationship between class selectivity, dimensionality, and robustness

abstract: While the relative trade-offs between sparse and distributed representations in deep neural networks (DNNs) are well-studied, less is known about how these trade-offs apply to representations of semantically-meaningful information. Class selectivity, the variability of a unit's responses across data classes or dimensions, is one way of quantifying the sparsity of semantic representations. Given recent evidence showing that class selectivity can impair generalization, we sought to investigate whether it also confers robustness (or vulnerability) to perturbations of input data. We found that mean class selectivity predicts vulnerability to naturalistic corruptions; networks regularized to have lower levels of class selectivity are more robust to corruption, while networks with higher class selectivity are more vulnerable to corruption, as measured using Tiny ImageNetC and CIFAR10C. In contrast, we found that class selectivity increases robustness to multiple types of gradient-based adversarial attacks. To examine this difference, we studied the dimensionality of the change in the representation due to perturbation, finding that decreasing class selectivity increases the dimensionality of this change for both corruption types, but with a notably larger increase for adversarial attacks. These results demonstrate the causal relationship between selectivity and robustness and provide new insights into the mechanisms of this relationship.


title: Sparse Regression for Extreme Values

abstract: We study the problem of selecting features associated with extreme values in high dimensional linear regression. Normally, in linear modeling problems, the presence of abnormal extreme values or outliers is considered an anomaly which should either be removed from the data or remedied using robust regression methods. In many situations, however, the extreme values in regression modeling are not outliers but rather the signals of interest; consider traces from spiking neurons, volatility in finance, or extreme events in climate science, for example. In this paper, we propose a new method for sparse high-dimensional linear regression for extreme values which is motivated by the Subbotin, or generalized normal distribution, which we call the extreme value linear regression model. For our method, we utilize an $\ell_p$ norm loss where $p$ is an even integer greater than two; we demonstrate that this loss increases the weight on extreme values. We prove consistency and variable selection consistency for the extreme value linear regression with a Lasso penalty, which we term the Extreme Lasso, and we also analyze the theoretical impact of extreme value observations on the model parameter estimates using the concept of influence functions. Through simulation studies and a real-world data example, we show that the Extreme Lasso outperforms other methods currently used in the literature for selecting features of interest associated with extreme values in high-dimensional regression.


title: Neutron Star Mergers as the Main Source of R-process: Natal Kicks And Inside-Out Evolution to The Rescue

abstract: Binary neutron star mergers (BNSMs) is currently the most promising source of \textsl{r}-process thanks to the detection of GW170817. The estimated occurring frequency and the amount of mass ejected per merger indicate that BNSMs by itself can account for all the \textsl{r}-process enrichment in the Galaxy. However, the decreasing trend of [Eu/Fe] vs [Fe/H] of disk stars for [Fe/H]$\gtrsim -1$ in the solar neighborhood is inconsistent with the flat trend expected from BNSMs with a standard delay time distribution (DTD) $\propto t^{-1}$. This has led to the suggestion that either additional sources or modification to the DTD of BNSMs is required to match the observations. We investigate the effects of natal kicks received during the birth of neutron star binaries on the chemical evolution of \textsl{r}-process element Eu in the Milky Way by combining the results from the galactic dynamics code \textsc{galpy} with a one-zone Galactic chemical evolution model \textsc{omega}. We show that when key inputs from simulations of the inside-out disk evolution are combined with natal kicks, BNSMs can naturally reproduce the observed decreasing trend of [Eu/Fe] with [Fe/H] in the solar neighborhood without the need for modification to the DTD or additional \textsl{r}-process sources.


title: Minimax Efficient Finite-Difference Stochastic Gradient Estimators Using Black-Box Function Evaluations

abstract: Standard approaches to stochastic gradient estimation, with only noisy black-box function evaluations, use the finite-difference method or its variants. While natural, it is open to our knowledge whether their statistical accuracy is the best possible. This paper argues so by showing that central finite-difference is a nearly minimax optimal zeroth-order gradient estimator for a suitable class of objective functions and mean squared risk, among both the class of linear estimators and the much larger class of all (nonlinear) estimators.


title: Are PETs (Privacy Enhancing Technologies) Giving Protection for Smartphones? -- A Case Study

abstract: With smartphone technologies enhanced way of interacting with the world around us, it has also been paving the way for easier access to our private and personal information. This has been amplified by the existence of numerous embedded sensors utilized by millions of apps to users. While mobile apps have positively transformed many aspects of our lives with new functionalities, many of these applications are taking advantage of vast amounts of data, privacy apps, a form of Privacy Enhancing Technology can be an effective privacy management tool for smartphones. To protect against vulnerabilities related to the collection, storage, and sharing of sensitive data, developers are building numerous privacy apps. However, there has been a lack of discretion in this particular area which calls for a proper assessment to understand the far-reaching utilization of these apps among users. During this process we have conducted an evaluation of the most popular privacy apps from our total collection of five hundred and twelve to demonstrate their functionality specific data protections they are claiming to offer, both technologically and conventionally, measuring up to standards. Taking their offered security functionalities as a scale, we conducted forensic experiments to indicate where they are failing to be consistent in maintaining protection. For legitimate validation of security gaps in assessed privacy apps, we have also utilized NIST and OWASP guidelines. We believe this study will be efficacious for continuous improvement and can be considered as a foundation towards a common standard for privacy and security measures for an app's development stage.


title: Estimation and inference on high-dimensional individualized treatment rule in observational data using split-and-pooled de-correlated score

abstract: With the increasing adoption of electronic health records, there is an increasing interest in developing individualized treatment rules, which recommend treatments according to patients' characteristics, from large observational data. However, there is a lack of valid inference procedures for such rules developed from this type of data in the presence of high-dimensional covariates. In this work, we develop a penalized doubly robust method to estimate the optimal individualized treatment rule from high-dimensional data. We propose a split-and-pooled de-correlated score to construct hypothesis tests and confidence intervals. Our proposal utilizes the data splitting to conquer the slow convergence rate of nuisance parameter estimations, such as non-parametric methods for outcome regression or propensity models. We establish the limiting distributions of the split-and-pooled de-correlated score test and the corresponding one-step estimator in high-dimensional setting. Simulation and real data analysis are conducted to demonstrate the superiority of the proposed method.


title: StructureBoost: Efficient Gradient Boosting for Structured Categorical Variables

abstract: Gradient boosting methods based on Structured Categorical Decision Trees (SCDT) have been demonstrated to outperform numerical and one-hot-encodings on problems where the categorical variable has a known underlying structure. However, the enumeration procedure in the SCDT is infeasible except for categorical variables with low or moderate cardinality. We propose and implement two methods to overcome the computational obstacles and efficiently perform Gradient Boosting on complex structured categorical variables. The resulting package, called StructureBoost, is shown to outperform established packages such as CatBoost and LightGBM on problems with categorical predictors that contain sophisticated structure. Moreover, we demonstrate that StructureBoost can make accurate predictions on unseen categorical values due to its knowledge of the underlying structure.


title: An Algorithm for Quantum Computation of Particle Decays

abstract: A quantum algorithm is developed to calculate decay rates and cross sections using quantum resources that scale polynomially in the system size assuming similar scaling for state preparation and time evolution. This is done by computing finite-volume one- and two-particle Green's functions on the quantum hardware. Particle decay rates and two particle scattering cross sections are extracted from the imaginary parts of the Green's function. A $0+1$ dimensional implementation of this method is demonstrated on IBM's superconducting quantum hardware for the decay of a heavy scalar particle to a pair of light scalars.


title: Emergence of Hierarchy in Networked Endorsement Dynamics

abstract: Many social and biological systems are characterized by enduring hierarchies, including those organized around prestige in academia, dominance in animal groups, and desirability in online dating. Despite their ubiquity, the general mechanisms that explain the creation and endurance of such hierarchies are not well understood. We introduce a generative model for the dynamics of hierarchies using time-varying networks in which new links are formed based on the preferences of nodes in the current network and old links are forgotten over time. The model produces a range of hierarchical structures, ranging from egalitarianism to bistable hierarchies, and we derive critical points that separate these regimes in the limit of long system memory. Importantly, our model supports statistical inference, allowing for a principled comparison of generative mechanisms using data. We apply the model to study hierarchical structures in empirical data on hiring patterns among mathematicians, dominance relations among parakeets, and friendships among members of a fraternity, observing several persistent patterns as well as interpretable differences in the generative mechanisms favored by each. Our work contributes to the growing literature on statistically grounded models of time-varying networks.


title: Searching for Efficient Architecture for Instrument Segmentation in Robotic Surgery

abstract: Segmentation of surgical instruments is an important problem in robot-assisted surgery: it is a crucial step towards full instrument pose estimation and is directly used for masking of augmented reality overlays during surgical procedures. Most applications rely on accurate real-time segmentation of high-resolution surgical images. While previous research focused primarily on methods that deliver high accuracy segmentation masks, majority of them can not be used for real-time applications due to their computational cost. In this work, we design a light-weight and highly-efficient deep residual architecture which is tuned to perform real-time inference of high-resolution images. To account for reduced accuracy of the discovered light-weight deep residual network and avoid adding any additional computational burden, we perform a differentiable search over dilation rates for residual units of our network. We test our discovered architecture on the EndoVis 2017 Robotic Instruments dataset and verify that our model is the state-of-the-art in terms of speed and accuracy tradeoff with a speed of up to 125 FPS on high resolution images.


title: T-REx: Table Repair Explanations

abstract: Data repair is a common and crucial step in many frameworks today, as applications may use data from different sources and of different levels of credibility. Thus, this step has been the focus of many works, proposing diverse approaches. To assist users in understanding the output of such data repair algorithms, we propose T-REx, a system for providing data repair explanations through Shapley values. The system is generic and not specific to a given repair algorithm or approach: it treats the algorithm as a black box. Given a specific table cell selected by the user, T-REx employs Shapley values to explain the significance of each constraint and each table cell in the repair of the cell of interest. T-REx then ranks the constraints and table cells according to their importance in the repair of this cell. This explanation allows users to understand the repair process, as well as to act based on this knowledge, to modify the most influencing constraints or the original database.


title: Online probabilistic label trees

abstract: We introduce online probabilistic label trees (OPLTs), an algorithm that trains a label tree classifier in a fully online manner without any prior knowledge about the number of training instances, their features and labels. OPLTs are characterized by low time and space complexity as well as strong theoretical guarantees. They can be used for online multi-label and multi-class classification, including the very challenging scenarios of one- or few-shot learning. We demonstrate the attractiveness of OPLTs in a wide empirical study on several instances of the tasks mentioned above.


title: NASGEM: Neural Architecture Search via Graph Embedding Method

abstract: Neural Architecture Search (NAS) automates and prospers the design of neural networks. Estimator-based NAS has been proposed recently to model the relationship between architectures and their performance to enable scalable and flexible search. However, existing estimator-based methods encode the architecture into a latent space without considering graph similarity. Ignoring graph similarity in node-based search space may induce a large inconsistency between similar graphs and their distance in the continuous encoding space, leading to inaccurate encoding representation and/or reduced representation capacity that can yield sub-optimal search results. To preserve graph correlation information in encoding, we propose NASGEM which stands for Neural Architecture Search via Graph Embedding Method. NASGEM is driven by a novel graph embedding method equipped with similarity measures to capture the graph topology information. By precisely estimating the graph distance and using an auxiliary Weisfeiler-Lehman kernel to guide the encoding, NASGEM can utilize additional structural information to get more accurate graph representation to improve the search efficiency. GEMNet, a set of networks discovered by NASGEM, consistently outperforms networks crafted by existing search methods in classification tasks, i.e., with 0.4%-3.6% higher accuracy while having 11%- 21% fewer Multiply-Accumulates. We further transfer GEMNet for COCO object detection. In both one-stage and twostage detectors, our GEMNet surpasses its manually-crafted and automatically-searched counterparts.


title: Just in Time: Personal Temporal Insights for Altering Model Decisions

abstract: The interpretability of complex Machine Learning models is coming to be a critical social concern, as they are increasingly used in human-related decision-making processes such as resume filtering or loan applications. Individuals receiving an undesired classification are likely to call for an explanation -- preferably one that specifies what they should do in order to alter that decision when they reapply in the future. Existing work focuses on a single ML model and a single point in time, whereas in practice, both models and data evolve over time: an explanation for an application rejection in 2018 may be irrelevant in 2019 since in the meantime both the model and the applicant's data can change. To this end, we propose a novel framework that provides users with insights and plans for changing their classification in particular future time points. The solution is based on combining state-of-the-art algorithms for (single) model explanations, ones for predicting future models, and database-style querying of the obtained explanations. We propose to demonstrate the usefulness of our solution in the context of loan applications, and interactively engage the audience in computing and viewing suggestions tailored for applicants based on their unique characteristic.


title: Explaining Natural Language Query Results

abstract: Multiple lines of research have developed Natural Language (NL) interfaces for formulating database queries. We build upon this work, but focus on presenting a highly detailed form of the answers in NL. The answers that we present are importantly based on the provenance of tuples in the query result, detailing not only the results but also their explanations. We develop a novel method for transforming provenance information to NL, by leveraging the original NL query structure. Furthermore, since provenance information is typically large and complex, we present two solutions for its effective presentation as NL text: one that is based on provenance factorization, with novel desiderata relevant to the NL case, and one that is based on summarization. We have implemented our solution in an end-to-end system supporting questions, answers and provenance, all expressed in NL. Our experiments, including a user study, indicate the quality of our solution and its scalability.


title: FEDReD I: 3D extinction and stellar maps by Bayesian deconvolution

abstract: Context. While Gaia enables to probe in great detail the extended local neighbourhood, the thin disk structure at larger distances remains sparsely explored. Aims. We aim here to build a non-parametric 3D model of the thin disc structures handling both the extinction and the stellar density simultaneously. Methods. We developed a Bayesian deconvolution method in two dimensions: extinction and distance. It uses a reference catalogue which completeness information defines the selection function. It is designed so that any complementary information from other catalogues can be added. It has also been designed to be robust to outliers, frequent in crowded fields, and differential extinction. The prior information is designed to be minimal: only a reference H-R diagram. We derived for this an empirical H-R diagram of the thin disk using Gaia DR2 data and synthetic isochrone-based H-R diagrams can also be used. Results. We validated the method on simulations and real fields using 2MASS and UKIDSS data complemented by Gaia DR2 photometry and parallaxes. We detail the results of two test fields: a 2MASS field centred around the NGC 4815 open cluster which shows an over-density of both extinction and stellar density at the cluster distance, and a UKIDSS field at $l=10^{\circ}$ where we recover the position of the Galactic bar.


title: An Efficient Data Imputation Technique for Human Activity Recognition

abstract: The tremendous applications of human activity recognition are surging its span from health monitoring systems to virtual reality applications. Thus, the automatic recognition of daily life activities has become significant for numerous applications. In recent years, many datasets have been proposed to train the machine learning models for efficient monitoring and recognition of human daily living activities. However, the performance of machine learning models in activity recognition is crucially affected when there are incomplete activities in a dataset, i.e., having missing samples in dataset captures. Therefore, in this work, we propose a methodology for extrapolating the missing samples of a dataset to better recognize the human daily living activities. The proposed method efficiently pre-processes the data captures and utilizes the k-Nearest Neighbors (KNN) imputation technique to extrapolate the missing samples in dataset captures. The proposed methodology elegantly extrapolated a similar pattern of activities as they were in the real dataset.


title: Accelerating Multigrid-based Hierarchical Scientific Data Refactoring on GPUs

abstract: Rapid growth in scientific data and a widening gap between computational speed and I/O bandwidth make it increasingly infeasible to store and share all data produced by scientific simulations. Instead, we need methods for reducing data volumes: ideally, methods that can scale data volumes adaptively so as to enable negotiation of performance and fidelity tradeoffs in different situations. Multigrid-based hierarchical data representations hold promise as a solution to this problem, allowing for flexible conversion between different fidelities so that, for example, data can be created at high fidelity and then transferred or stored at lower fidelity via logically simple and mathematically sound operations. However, the effective use of such representations has been hindered until now by the relatively high costs of creating, accessing, reducing, and otherwise operating on such representations. We describe here highly optimized data refactoring kernels for GPU accelerators that enable efficient creation and manipulation of data in multigrid-based hierarchical forms. We demonstrate that our optimized design can achieve up to 250 TB/s aggregated data refactoring throughput -- 83% of theoretical peak -- on 1024 nodes of the Summit supercomputer. We showcase our optimized design by applying it to a large-scale scientific visualization workflow and the MGARD lossy compression software.


title: Robust Bayesian Classification Using an Optimistic Score Ratio

abstract: We build a Bayesian contextual classification model using an optimistic score ratio for robust binary classification when there is limited information on the class-conditional, or contextual, distribution. The optimistic score searches for the distribution that is most plausible to explain the observed outcomes in the testing sample among all distributions belonging to the contextual ambiguity set which is prescribed using a limited structural constraint on the mean vector and the covariance matrix of the underlying contextual distribution. We show that the Bayesian classifier using the optimistic score ratio is conceptually attractive, delivers solid statistical guarantees and is computationally tractable. We showcase the power of the proposed optimistic score ratio classifier on both synthetic and empirical data.


title: Meta-Learning for One-Class Classification with Few Examples using Order-Equivariant Network

abstract: This paper presents a meta-learning framework for few-shots One-Class Classification (OCC) at test-time, a setting where labeled examples are only available for the positive class, and no supervision is given for the negative example. We consider that we have a set of `one-class classification' objective-tasks with only a small set of positive examples available for each task, and a set of training tasks with full supervision (i.e. highly imbalanced classification). We propose an approach using order-equivariant networks to learn a 'meta' binary-classifier. The model will take as input an example to classify from a given task, as well as the corresponding supervised set of positive examples for this OCC task. Thus, the output of the model will be 'conditioned' on the available positive example of a given task, allowing to predict on new tasks and new examples without labeled negative examples. In this paper, we are motivated by an astronomy application. Our goal is to identify if stars belong to a specific stellar group (the 'one-class' for a given task), called \textit{stellar streams}, where each stellar stream is a different OCC-task. We show that our method transfers well on unseen (test) synthetic streams, and outperforms the baselines even though it is not retrained and accesses a much smaller part of the data per task to predict (only positive supervision). We see however that it doesn't transfer as well on the real stream GD-1. This could come from intrinsic differences from the synthetic and real stream, highlighting the need for consistency in the 'nature' of the task for this method. However, light fine-tuning improve performances and outperform our baselines. Our experiments show encouraging results to further explore meta-learning methods for OCC tasks.


title: Covariant vs Contravariant Methods in Differential Geometry

abstract: This is a short essay about some fundamental results on scalar curvature and the two key methods that are used to establish them.


title: Phenomenological model of motility by spatiotemporal modulation of active interactions

abstract: Transport at microscopic length scales is essential in biological systems and various technologies, including microfluidics. Recent experiments achieved self-organized transport phenomena in microtubule active matter using light to modulate motor-protein activity in time and space. Here, we introduce a novel phenomenological model to explain such experiments. Our model, based on spatially modulated particle interactions, reveals a possible mechanism for emergent transport phenomena in light-controlled active matter, including motility and contraction. In particular, the model's analytic treatment elucidates the conservation of the center of mass of activated particles as a fundamental mechanism of material transport and demonstrates the necessity of memory for sustained motility. Furthermore, we generalize the model to explain other phenomena, like microtubule aster-aster interactions induced by more complicated activation geometries. Our results demonstrate that the model provides a possible foundation for the phenomenological understanding of light-controlled active matter, and it will enable the design and optimization of transport protocols for active matter devices.


title: Scalable Computations of Wasserstein Barycenter via Input Convex Neural Networks

abstract: Wasserstein Barycenter is a principled approach to represent the weighted mean of a given set of probability distributions, utilizing the geometry induced by optimal transport. In this work, we present a novel scalable algorithm to approximate the Wasserstein Barycenters aiming at high-dimensional applications in machine learning. Our proposed algorithm is based on the Kantorovich dual formulation of the Wasserstein-2 distance as well as a recent neural network architecture, input convex neural network, that is known to parametrize convex functions. The distinguishing features of our method are: i) it only requires samples from the marginal distributions; ii) unlike the existing approaches, it represents the Barycenter with a generative model and can thus generate infinite samples from the barycenter without querying the marginal distributions; iii) it works similar to Generative Adversarial Model in one marginal case. We demonstrate the efficacy of our algorithm by comparing it with the state-of-art methods in multiple experiments.


title: Charging a quantum battery via non equilibrium heat current

abstract: When a quantum system is subject to a thermal gradient it may sustain a steady non-equilibrium heat current, by entering into a so-called non equilibrium steady state (NESS). Here we show that NESS constitute a thermodynamic resource that can be exploited to fuel a quantum heat engine. This adds to the list of recently reported sources available at the nano-scale, such as coherence, entanglement and quantum measurements. We elucidate this concept by showing analytic and numerical studies of a two-qubits quantum battery that is alternatively charged by a thermal gradient and discharged by application of a properly chosen unitary gate. The presence of a NESS for the charging step guarantees steady operation with positive power output. Decreasing the duration of the charging step results in a time periodic steady state accompanied by increased efficiency and output power. The device is amenable to implementation with different nanotechnology platforms.


title: Deform, Cut and Tear a skinned model using Conformal Geometric Algebra

abstract: In this work, we present a novel, integrated rigged character simulation framework in Conformal Geometric Algebra (CGA) that supports, for the first time, real-time cuts and tears, before and/or after the animation, while maintaining deformation topology. The purpose of using CGA is to lift several restrictions posed by current state-of-the-art character animation & deformation methods. Previous implementations originally required weighted matrices to perform deformations, whereas, in the current state-of-the-art, dual-quaternions handle both rotations and translations, but cannot handle dilations. CGA is a suitable extension of dual-quaternion algebra that amends these two major previous shortcomings: the need to constantly transmute between matrices and dual-quaternions as well as the inability to properly dilate a model during animation. Our CGA algorithm also provides easy interpolation and application of all deformations in each intermediate steps, all within the same geometric framework. Furthermore we also present two novel algorithms that enable cutting and tearing of the input rigged, animated model, while the output model can be further re-deformed. These interactive, real-time cut and tear operations can enable a new suite of applications, especially under the scope of a medical surgical simulation.


title: Discovery of a transient X-ray source Suzaku J1305-4930 in NGC 4945

abstract: We report the serendipitous discovery of a transient X-ray source, Suzaku J1305$-$4930, $\sim$3 kpc southwest of the nucleus of the Seyfert 2 galaxy NGC 4945. Among the seven Suzaku observations of NGC 4945 from 2005 to 2011, Suzaku J1305$-$4930 was detected four times in July and August in 2010. The X-ray spectra are better approximated with a multi-color disk model than a power-law model. At the first detection on 2010 July 4--5, its X-ray luminosity was $(8.9^{+0.2}_{-0.4}) \times 10^{38}$ erg s$^{-1}$ and the temperature at the inner-disk radius ($kT_{\rm in}$) was $1.12\pm0.04$ keV. At the last detection with Suzaku on 2010 August 4--5, the luminosity decreased to $(2.2^{+0.3}_{-0.8}) \times10^{38}$ erg s$^{-1}$ and $kT_{\rm in}$ was $0.62\pm0.07$ keV. The source was not detected on 2011 January 29, about six months after the first detection, with a luminosity upper limit of $2.4\times10^{38}$ erg s$^{-1}$. We also find an absorption feature which is similar to that reported in Cyg X-1. Assuming the standard disk, we suggest that Suzaku J1305$-$4930 consists of a black hole with a mass of $\sim$10 $M_{\odot}$. The relation between the disk luminosity and $kT_{\rm in}$ is not reproduced with the standard model of a constant inner radius but is better approximated with a slim-disk model.


title: URSABench: Comprehensive Benchmarking of Approximate Bayesian Inference Methods for Deep Neural Networks

abstract: While deep learning methods continue to improve in predictive accuracy on a wide range of application domains, significant issues remain with other aspects of their performance including their ability to quantify uncertainty and their robustness. Recent advances in approximate Bayesian inference hold significant promise for addressing these concerns, but the computational scalability of these methods can be problematic when applied to large-scale models. In this paper, we describe initial work on the development ofURSABench(the Uncertainty, Robustness, Scalability, and Accu-racy Benchmark), an open-source suite of bench-marking tools for comprehensive assessment of approximate Bayesian inference methods with a focus on deep learning-based classification tasks


title: A new entropy-variable-based discretization method for minimum entropy moment approximations of linear kinetic equations

abstract: In this contribution we derive and analyze a new numerical method for kinetic equations based on a variable transformation of the moment approximation. Classical minimum-entropy moment closures are a class of reduced models for kinetic equations that conserve many of the fundamental physical properties of solutions. However, their practical use is limited by their high computational cost, as an optimization problem has to be solved for every cell in the space-time grid. In addition, implementation of numerical solvers for these models is hampered by the fact that the optimization problems are only well-defined if the moment vectors stay within the realizable set. For the same reason, further reducing these models by, e.g., reduced-basis methods is not a simple task. Our new method overcomes these disadvantages of classical approaches. The transformation is performed on the semi-discretized level which makes them applicable to a wide range of kinetic schemes and replaces the nonlinear optimization problems by inversion of the positive-definite Hessian matrix. As a result, the new scheme gets rid of the realizability-related problems. Moreover, a discrete entropy law can be enforced by modifying the time stepping scheme. Our numerical experiments demonstrate that our new method is often several times faster than the standard optimization-based scheme.


title: FPT and kernelization algorithms for the k-in-a-tree problem

abstract: The three-in-a-tree problem asks for an induced tree of the input graph containing three mandatory vertices. In 2006, Chudnovsky and Seymour [Combinatorica, 2010] presented the first polynomial time algorithm for this problem, which has become a critical subroutine in many algorithms for detecting induced subgraphs, such as beetles, pyramids, thetas, and even and odd-holes. In 2007, Derhy and Picouleau [Discrete Applied Mathematics, 2009] considered the natural generalization to $k$ mandatory vertices, proving that, when $k$ is part of the input, the problem is $\mathsf{NP}$-complete, and ask what is the complexity of four-in-a-tree. Motivated by this question and the relevance of the original problem, we study the parameterized complexity of $k$-in-a-tree. We begin by showing that the problem is $\mathsf{W[1]}$-hard when jointly parameterized by the size of the solution and minimum clique cover and, under the Exponential Time Hypothesis, does not admit an $n^{o(k)}$ time algorithm. Afterwards, we use Courcelle's Theorem to prove fixed-parameter tractability under cliquewidth, which prompts our investigation into which parameterizations admit single exponential algorithms; we show that such algorithms exist for the unrelated parameterizations treewidth, distance to cluster, and distance to co-cluster. In terms of kernelization, we present a linear kernel under feedback edge set, and show that no polynomial kernel exists under vertex cover nor distance to clique unless $\mathsf{NP} \subseteq \mathsf{coNP}/\mathsf{poly}$. Along with other remarks and previous work, our tractability and kernelization results cover many of the most commonly employed parameters in the graph parameter hierarchy.


title: Recognising Graphic and Matroidal Connectivity Functions

abstract: A {\em connectivity function} on a set $E$ is a function $\lambda:2^E\rightarrow \mathbb R$ such that $\lambda(\emptyset)=0$, that $\lambda(X)=\lambda(E-X)$ for all $X\subseteq E$, and that $\lambda(X\cap Y)+\lambda(X\cup Y)\leq \lambda(X)+\lambda(Y)$ for all $X,Y \subseteq E$. Graphs, matroids and, more generally, polymatroids have associated connectivity functions. In this paper we give a method for identifying when a connectivity function comes from a graph. This method uses no more than a polynomial number of evaluations of the connectivity function. In contrast, we show that the problem of identifying when a connectivity function comes from a matroid cannot be solved in polynomial time. We also show that the problem of identifying when a connectivity function is not that of a matroid cannot be solved in polynomial time.


title: Finite mixture models do not reliably learn the number of components

abstract: Scientists and engineers are often interested in learning the number of subpopulations (or components) present in a data set. A common suggestion is to use a finite mixture model (FMM) with a prior on the number of components. Past work has shown the resulting FMM component-count posterior is consistent; that is, the posterior concentrates on the true, generating number of components. But consistency requires the assumption that the component likelihoods are perfectly specified, which is unrealistic in practice. In this paper, we add rigor to data-analysis folk wisdom by proving that under even the slightest model misspecification, the FMM component-count posterior diverges: the posterior probability of any particular finite number of components converges to 0 in the limit of infinite data. Contrary to intuition, posterior-density consistency is not sufficient to establish this result. We develop novel sufficient conditions that are more realistic and easily checkable than those common in the asymptotics literature. We illustrate practical consequences of our theory on simulated and real data.


title: Properties of fractional integral operators involving the three-parameters Mittag-Leffler function in the kernels with respect to another function

abstract: This paper aims to investigate properties associated with fractional integral operators involving the three-parameters Mittag-Leffler function in the kernels with respect to another function. We prove that the Cauchy problem and the Volterra integral equation are equivalent. We find a closed-form to the solution of the Cauchy problem using successive approximations method and $\psi$-Caputo fractional derivative.


title: Evaluation of Adversarial Training on Different Types of Neural Networks in Deep Learning-based IDSs

abstract: Network security applications, including intrusion detection systems of deep neural networks, are increasing rapidly to make detection task of anomaly activities more accurate and robust. With the rapid increase of using DNN and the volume of data traveling through systems, different growing types of adversarial attacks to defeat them create a severe challenge. In this paper, we focus on investigating the effectiveness of different evasion attacks and how to train a resilience deep learning-based IDS using different Neural networks, e.g., convolutional neural networks (CNN) and recurrent neural networks (RNN). We use the min-max approach to formulate the problem of training robust IDS against adversarial examples using two benchmark datasets. Our experiments on different deep learning algorithms and different benchmark datasets demonstrate that defense using an adversarial training-based min-max approach improves the robustness against the five well-known adversarial attack methods.


title: Myth of scattering in finite volume

abstract: In this notes, we illustrate why the infinite volume scattering amplitude is in fact dispensable when it comes to formulating few-body quantization condition in finite volume. Only subprocess interactions or interactions associated subprocess amplitudes are essential and fundamental ingredients of quantization conditions. After these ingredients are determined, infinite volume scattering amplitude can be computed separately. The underlying reasons are rooted in facts that (1) the final physical process is generated by all subprocess or interactions among particles; (2) the ultimate goal of quantization condition in finite volume is to find stationary solutions of few-body system. That is to say, in the end, it all comes down to the solving of eigenvalue problem, $\hat{H} | n \rangle = E_n | n \rangle $.


title: Instantons and Bows for the Classical Groups

abstract: The construction of Atiyah, Drinfeld, Hitchin, and Manin [ADHM78] provided complete description of all instantons on Euclidean four-space. It was extended by Kronheimer and Nakajima to instantons on ALE spaces, resolutions of orbifolds $\mathbb{R}^4/\Gamma$ by a finite subgroup $\Gamma\subset SU(2).$ We consider a similar classification, in the holomorphic context, of instantons on some of the next spaces in the hierarchy, the ALF multi-Taub-NUT manifolds, showing how they tie in to the bow solutions to Nahm's equations [Che09] via the Nahm correspondence. Recently in [Nak18a] and [NT17], based on [Nak03], Nakajima and Takayama constructed the Coulomb branch of the moduli space of vacua of a quiver gauge theory, tying them to the same space of bow solutions. One can view our construction as describing the same manifold as the Higgs branch of the mirror gauge theory [COS11]. Our construction also yields the monad construction of holomorphic instanton bundles on the multi-Taub-NUT space for any classical compact Lie structure group.


title: Mapping oscillating magnetic fields around rechargeable batteries

abstract: Power storage devices such as batteries are a crucial part of modern technology. The development and use of batteries has accelerated in the past decades, yet there are only a few techniques that allow gathering vital information from battery cells in a nonivasive fashion. A widely used technique to investigate batteries is electrical impedance spectroscopy (EIS), which provides information on how the impedance of a cell changes as a function of the frequency of applied alternating currents. Building on recent developments of inside-out MRI (ioMRI), a technique is presented here which produces spatially-resolved maps of the oscillating magnetic fields originating from the alternating electrical currents distributed within a cell. The technique works by using an MRI pulse sequence synchronized with a gated alternating current applied to the cell terminals. The approach is benchmarked with a current-carrying wire coil, and demonstrated with commercial and prototype lithium-ion cells. Marked changes in the fields are observed for different cell types.


title: $\eta_{c}$-nucleus bound states

abstract: $\eta_c$-nucleus bound state energies are calculated for various nuclei. Essential input for the calculations, namely the medium-modified $D$ and $D^{*}$ meson masses, as well as the density distributions in nuclei, are calculated within the quark-meson coupling (QMC) model. The attractive potentials for the $\eta_c$ meson in the nuclear medium originate from the in-medium enhanced $DD^{*}$ loops in the $\eta_c$ self-energy. Our results suggest that the $\eta_c$ meson should form bound states with all the nuclei considered.


title: Good AI for the Present of Humanity Democratizing AI Governance

abstract: What do Cyberpunk and AI Ethics have to do with each other? Cyberpunk is a sub-genre of science fiction that explores the post-human relationships between human experience and technology. One similarity between AI Ethics and Cyberpunk literature is that both seek to explore future social and ethical problems that our technological advances may bring upon society. In recent years, an increasing number of ethical matters involving AI have been pointed and debated, and several ethical principles and guides have been suggested as governance policies for the tech industry. However, would this be the role of AI Ethics? To serve as a soft and ambiguous version of the law? We would like to advocate in this article for a more Cyberpunk way of doing AI Ethics, with a more democratic way of governance. In this study, we will seek to expose some of the deficits of the underlying power structures of the AI industry, and suggest that AI governance be subject to public opinion, so that good AI can become good AI for all.


title: Closing the Random Graph Gap in Tuza's Conjecture Through the Online Triangle Packing Process

abstract: A long-standing conjecture of Zsolt Tuza asserts that the triangle covering number $\tau(G)$ is at most twice the triangle packing number $\nu(G)$, where the triangle packing number $\nu(G)$ is the maximum size of a set of edge-disjoint triangles in $G$ and the triangle covering number $\tau(G)$ is the minimal size of a set of edges intersecting all triangles. In this paper, we prove that Tuza's conjecture holds in the Erd\H{o}s-R\'enyi random graph $G(n,m)$ for all range of $m$, closing the gap in what was previously known. (Recently, this result was also independently proved by Jeff Kahn and Jinyoung Park.) We employ a random greedy process called the online triangle packing process to produce a triangle packing in $G(n,m)$ and analyze this process by using the differential equations method.


title: Signless Laplacian spectral radius and matching in graphs

abstract: The signless Laplacian matrix of a graph $G$ is given by $Q(G)=D(G)+A(G)$, where $D(G)$ is a diagonal matrix of vertex degrees and $A(G)$ is the adjacency matrix. The largest eigenvalue of $Q(G)$ is called the signless Laplacian spectral radius, denoted by $q_1=q_1(G)$. In this paper, some properties between the signless Laplacian spectral radius and perfect matching in graphs are establish. Let $r(n)$ be the largest root of equation $x^3-(3n-7)x^2+n(2n-7)x-2(n^2-7n+12)=0$. We show that $G$ has a perfect matching for $n=4$ or $n\geq10$, if $q_1(G)>r(n)$, and for $n=6$ or $n=8$, if $q_1(G)>4+2\sqrt{3}$ or $q_1(G)>6+2\sqrt{6}$ respectively, where $n$ is a positive even integer number. Moreover, there exists graphs $K_{n-3}\vee K_1 \vee \overline{K_2}$ such that $q_1(K_{n-3}\vee K_1 \vee \overline{K_2})=r(n)$ if $n\geq4$, a graph $K_2\vee\overline{K_4}$ such that $q_1(K_2\vee\overline{K_4})=4+2\sqrt{3}$ and a graph $K_3\vee\overline{K_5}$ such that $q_1(K_3\vee\overline{K_5})=6+2\sqrt{6}$. These graphs all have no prefect matching.


title: Automatic Probe Movement Guidance for Freehand Obstetric Ultrasound

abstract: We present the first system that provides real-time probe movement guidance for acquiring standard planes in routine freehand obstetric ultrasound scanning. Such a system can contribute to the worldwide deployment of obstetric ultrasound scanning by lowering the required level of operator expertise. The system employs an artificial neural network that receives the ultrasound video signal and the motion signal of an inertial measurement unit (IMU) that is attached to the probe, and predicts a guidance signal. The network termed US-GuideNet predicts either the movement towards the standard plane position (goal prediction), or the next movement that an expert sonographer would perform (action prediction). While existing models for other ultrasound applications are trained with simulations or phantoms, we train our model with real-world ultrasound video and probe motion data from 464 routine clinical scans by 17 accredited sonographers. Evaluations for 3 standard plane types show that the model provides a useful guidance signal with an accuracy of 88.8% for goal prediction and 90.9% for action prediction.


title: Multi-dimensional backward stochastic differential equations of diagonally quadratic generators: the general result

abstract: This paper is devoted to a general solvability of a multi-dimensional backward stochastic differential equation (BSDE) of a diagonally quadratic generator $g(t,y,z)$, by relaxing the assumptions of \citet{HuTang2016SPA} on the generator and terminal value. More precisely, the generator $g(t,y,z)$ can have more general growth and continuity in $y$ in the local solution; while in the global solution, the generator $g(t,y,z)$ can have a skew sub-quadratic but in addition "strictly and diagonally" quadratic growth in the second unknown variable $z$, or the terminal value can be unbounded but the generator $g(t,y,z)$ is "diagonally dependent" on the second unknown variable $z$ (i.e., the $i$-th component $g^i$ of the generator $g$ only depends on the $i$-th row $z^i$ of the variable $z$ for each $i=1,\cdots,n$ ). Three new results are established on the local and global solutions when the terminal value is bounded and the generator $g$ is subject to some general assumptions. When the terminal value is unbounded but is of exponential moments of arbitrary order, an existence and uniqueness result is given under the assumptions that the generator $g(t,y,z)$ is Lipschitz continuous in the first unknown variable $y$, and varies with the second unknown variable $z$ in a "diagonal" , "component-wisely convex or concave", and "quadratically growing" way, which seems to be the first general solvability of systems of quadratic BSDEs with unbounded terminal values. This generalizes and strengthens some existing results via some new ideas.


title: Nanoscale Real-Time Detection of Quantum Vortices at Millikelvin Temperatures

abstract: Since we still lack a theory of classical turbulence, attention has focused on the conceptually simpler turbulence in quantum fluids. Can such systems of identical singly-quantized vortices provide a physically accessible "toy model" of the classical counterpart? That said, we have hitherto lacked detectors capable of the real-time, non-invasive probing of the wide range of length scales involved in quantum turbulence. However, we demonstrate here the real-time detection of quantum vortices by a nanoscale resonant beam in superfluid $^4$He at 10 mK. The basic idea is that we can trap a single vortex along the length of a nanobeam and observe the transitions as a vortex is either trapped or released, which we observe through the shift in the resonant frequency of the beam. With a tuning fork source, we can control the ambient vorticity density and follow its influence on the vortex capture and release rates. But, most important, we show that these devices are capable of probing turbulence on the micron scale.

categories: cond-mat.mes-hall cond-mat.quant-gas physics.flu-dyn physics.ins-det


title: Classification of simple Harish-Chandra modules over the N=1 Ramond algebra

abstract: In this paper, we give a new approach to classify all simple Harish-Chandra modules for the N=1 Ramond algebra based on the so called A-cover theory developed in \cite{BF}


title: Transparency Tools for Fairness in AI (Luskin)

abstract: We propose new tools for policy-makers to use when assessing and correcting fairness and bias in AI algorithms. The three tools are: - A new definition of fairness called "controlled fairness" with respect to choices of protected features and filters. The definition provides a simple test of fairness of an algorithm with respect to a dataset. This notion of fairness is suitable in cases where fairness is prioritized over accuracy, such as in cases where there is no "ground truth" data, only data labeled with past decisions (which may have been biased). - Algorithms for retraining a given classifier to achieve "controlled fairness" with respect to a choice of features and filters. Two algorithms are presented, implemented and tested. These algorithms require training two different models in two stages. We experiment with combinations of various types of models for the first and second stage and report on which combinations perform best in terms of fairness and accuracy. - Algorithms for adjusting model parameters to achieve a notion of fairness called "classification parity". This notion of fairness is suitable in cases where accuracy is prioritized. Two algorithms are presented, one which assumes that protected features are accessible to the model during testing, and one which assumes protected features are not accessible during testing. We evaluate our tools on three different publicly available datasets. We find that the tools are useful for understanding various dimensions of bias, and that in practice the algorithms are effective in starkly reducing a given observed bias when tested on new data.


title: Computing homological residue fields in algebra and topology

abstract: We determine the homological residue fields, in the sense of tensor-triangular geometry, in a series of concrete examples ranging from topological stable homotopy theory to modular representation theory of finite groups.


title: Making learning more transparent using conformalized performance prediction

abstract: In this work, we study some novel applications of conformal inference techniques to the problem of providing machine learning procedures with more transparent, accurate, and practical performance guarantees. We provide a natural extension of the traditional conformal prediction framework, done in such a way that we can make valid and well-calibrated predictive statements about the future performance of arbitrary learning algorithms, when passed an as-yet unseen training set. In addition, we include some nascent empirical examples to illustrate potential applications.


title: Invertibility of Multi-Energy X-ray Transform

abstract: Purpose: The goal is to provide a sufficient condition on the invertibility of a multi-energy (ME) X-ray transform. The energy-dependent X-ray attenuation profiles can be represented by a set of coefficients using the Alvarez-Macovski (AM) method. An ME X-ray transform is a mapping from $N$ AM coefficients to $N$ noise-free energy-weighted measurements, where $N\geq2$. Methods: We apply a general invertibility theorem which tests whether the Jacobian of the mapping $J(\mathbf A)$ has zero values over the support of the mapping. The Jacobian of an arbitrary ME X-ray transform is an integration over all spectral measurements. A sufficient condition of $J(\mathbf A)\neq0$ for all $\mathbf A$ is that the integrand of $J(\mathbf A)$ is $\geq0$ (or $\leq0$) everywhere. Note that the trivial case of the integrand equals to zero everywhere is ignored. With symmetry, we simplified the integrand of the Jacobian into three factors that are determined by the total attenuation, the basis functions, and the energy-weighting functions, respectively. The factor related to total attenuation is always positive, hence the invertibility of the X-ray transform can be determined by testing the signs of the other two factors. Furthermore, we use the Cramer-Rao lower bound (CRLB) to characterize the noise-induced estimation uncertainty and provide a maximum-likelihood (ML) estimator. Conclusions: We have provided a framework to study the invertibility of an arbitrary ME X-ray transform and proved the global invertibility for four types of systems.


title: Two approaches to quantum gravity and M-(atrix) theory at large number of dimensions

abstract: A Gaussian approximation to the bosonic part of M-(atrix) theory with mass deformation is considered at large values of the dimension $d$. From the perspective of the gauge/gravity duality this action reproduces with great accuracy the stringy Hagedorn phase transition from a confinement (black string) phase to a deconfinement (black hole) phase whereas from the perspective of the matrix/geometry approach this action only captures a remnant of the geometric Yang-Mills-to-fuzzy-sphere phase where the fuzzy sphere solution is only manifested as a three-cut configuration termed the "baby fuzzy sphere" configuration. The Yang-Mills phase retains most of its characteristics with two exceptions: i) the uniform distribution inside a solid ball suffers a crossover at very small values of the gauge coupling constant to a Wigner's semi-circle law, and ii) the uniform distribution at small values of the temperatures is non-existent.


title: Neutral monism, perspectivism and the quantum dualism: An essay

abstract: Quantum mechanics in the Wigner-von Neumann interpretation is presented. This is characterized by 1) a quantum dualism between matter and consciousness unified within an informational neutral monism, 2) a quantum perspectivism which is extended to a complementarity between the Copenhagen interpretation and the many-worlds formalism, 3) a psychophysical causal closure akin to Leibniz parallelism and 4) a quantum solipsism, i.e. a reality in which classical states are only potentially-existing until a conscious observation is made.


title: Artificial Intelligence and Machine Learning in 5G Network Security: Opportunities, advantages, and future research trends

abstract: Recent technological and architectural advancements in 5G networks have proven their worth as the deployment has started over the world. Key performance elevating factor from access to core network are softwareization, cloudification and virtualization of key enabling network functions. Along with the rapid evolution comes the risks, threats and vulnerabilities in the system for those who plan to exploit it. Therefore, ensuring fool proof end-to-end (E2E) security becomes a vital concern. Artificial intelligence (AI) and machine learning (ML) can play vital role in design, modelling and automation of efficient security protocols against diverse and wide range of threats. AI and ML has already proven their effectiveness in different fields for classification, identification and automation with higher accuracy. As 5G networks' primary selling point has been higher data rates and speed, it will be difficult to tackle wide range of threats from different points using typical/traditional protective measures. Therefore, AI and ML can play central role in protecting highly data-driven softwareized and virtualized network components. This article presents AI and ML driven applications for 5G network security, their implications and possible research directions. Also, an overview of key data collection points in 5G architecture for threat classification and anomaly detection are discussed.


title: Decay estimates for nonlinear Schr\"odinger equations

abstract: In this short note, we present some decay estimates for nonlinear solutions of 3d quintic, 3d cubic and 2d quintic NLS (nonlinear Schr\"odinger equations).


title: The SAMI -- Fornax Dwarfs Survey I: Sample, observations and the specific stellar angular momentum of dwarf elliptical galaxies

abstract: Dwarf ellipticals are the most common galaxy type in cluster environments, however the challenges associated with their observation mean their formation mechanisms are still poorly understood. To address this, we present deep integral field observations of a sample of 31 low-mass ($10^{7.5} <$ M$_\star < 10^{9.5}$ M$_\odot$) early-type galaxies in the Fornax cluster with the SAMI instrument. For 21 galaxies our observations are sufficiently deep to construct spatially resolved maps of the stellar velocity and velocity dispersion - for the remaining galaxies we extract global velocities and dispersions from aperture spectra only. From the kinematic maps we measure the specific stellar angular momentum $\lambda_R$ of the lowest mass dE galaxies to date. Combining our observations with early-type galaxy data from the literature spanning a large range in stellar mass, we find that $\lambda_R$ decreases towards lower stellar mass, with a corresponding increase in the proportion of slowly rotating galaxies in this regime. The decrease of $\lambda_R$ with mass in our sample dE galaxies is consistent with a similar trend seen in somewhat more massive spiral galaxies from the CALIFA survey. This suggests that the degree of dynamical heating required to produce dEs from low-mass starforming progenitors may be relatively modest, and consistent with a broad range of formation mechanisms.


title: The prescribed curvature problem for entire hypersurfaces in Minkowski space

abstract: We prove three results in this paper. First, we prove for a wide class of functions $\varphi\in C^2(\mathbb{S}^{n-1})$ and $\psi(X, \nu)\in C^2(\mathbb{R}^{n+1}\times\mathbb{H}^n),$ there exists a unique, entire, strictly convex, spacelike hypersurface $M_u$ satisfying $\sigma_k(\kappa[M_u])=\psi(X, \nu)$ and $u(x)\rightarrow |x|+\varphi\left(\frac{x}{|x|}\right)$ as $|x|\rightarrow\infty.$ Second, when $k=n-1, n-2,$ we show the existence and uniqueness of entire, $k$-convex, spacelike hypersurface $M_u$ satisfying $\sigma_k(\kappa[M_u])=\psi(x, u(x))$ and $u(x)\rightarrow |x|+\varphi\left(\frac{x}{|x|}\right)$ as $|x|\rightarrow\infty.$ Last, we obtain the existence and uniqueness of entire, strictly convex, downward translating solitons $M_u$ with prescribed asymptotic behavior at infinity for $\sigma_k$ curvature flow equations. Moreover, we prove that the downward translating solitons $M_u$ have bounded principal curvatures.


title: Density dependent quark mass model revisited: Thermodynamic consistency, stability windows and stellar properties

abstract: In this work a density dependent quark model is revisited, its thermodynamic consistency checked and the stability window for absolutely stable quark matter obtained. The hypotheses of both pure quark matter with equal quark chemical potentials and stellar matter subject to chemical stability and charge neutrality are investigated. The parameters that appear in the density dependent mass and satisfy the Bodmer-Witten conjecture are then used to compute the masses and radii of strange stars. We show that the obtained values are compatible with the recently observed massive stars.


title: Hack.VR: A Programming Game in Virtual Reality

abstract: In this article we describe Hack.VR, an object-oriented programming game in virtual reality. Hack.VR uses a VR programming language in which nodes represent functions and node connections represent data flow. Using this programming framework, players reprogram VR objects such as elevators, robots, and switches. Hack.VR has been designed to be highly interactable both physically and semantically.


title: Fast neural Poincar\'e maps for toroidal magnetic fields

abstract: Poincar\'e maps for toroidal magnetic fields are routinely employed to study gross confinement properties in devices built to contain hot plasmas. In most practical applications, evaluating a Poincar\'e map requires numerical integration of a magnetic field line, a process that can be slow and that cannot be easily accelerated using parallel computations. We show that a novel neural network architecture, the H\'enonNet, is capable of accurately learning realistic Poincar\'e maps from observations of a conventional field-line-following algorithm. After training, such learned Poincar\'e maps evaluate much faster than the field-line integration method. Moreover, the H\'enonNet architecture exactly reproduces the primary physics constraint imposed on field-line Poincar\'e maps: flux preservation. This structure-preserving property is the consequence of each layer in a H\'enonNet being a symplectic map. We demonstrate empirically that a H\'enonNet can learn to mock the confinement properties of a large magnetic island by using coiled hyperbolic invariant manifolds to produce a sticky chaotic region at the desired island location. This suggests a novel approach to designing magnetic fields with good confinement properties that may be more flexible than ensuring confinement using KAM tori.


title: Fast Molecular Compression by a Hyperthermal Collision Gives Bond-Selective Mechanochemistry

abstract: Using electrospray ion beam deposition, we collide the complex molecule Reichardt's Dye (C41H30NO+) at low, hyperthermal translational energy (2-50 eV) with a Cu(100) surface and image the outcome at single-molecule level by scanning tunneling microscopy. We observe bond-selective reaction induced by the translational kinetic energy. The collision impulse compresses the molecule and bends specific bonds, prompting them to react selectively. This dynamics drives the system to seek thermally inaccessible reactive pathways, since the compression timescale (sub-ps) is much shorter than the thermalization timescale (ns), thereby yielding reaction products that are unobtainable thermally.


title: Resonator-assisted single molecule quantum state detection

abstract: We propose a state-sensitive scheme to optically detect a single molecule without a closed transition, through strong coupling to a high-Q whispering-gallery mode high-Q resonator. A background-free signal can be obtained by detecting a molecule-induced transparency in a photon bus waveguide that is critically coupled to the resonator, with a suppressed depumping rate to other molecular states by the cooperativity parameter $C$. We numerically calculate the dynamics of the molecule-resonator coupled system using Lindblad master equations, and develop analytical solutions through the evolution of quasi-steady states in the weak-driving regime. Using Rb$_2$ triplet ground state molecules as an example, we show that high fidelity state readout can be achieved using realistic resonator parameters. We further discuss the case of multiple molecules collectively coupled to a resonator, demonstrating near-unity detection fidelity and negligible population loss.


title: Robotic Grasping using Deep Reinforcement Learning

abstract: In this work, we present a deep reinforcement learning based method to solve the problem of robotic grasping using visio-motor feedback. The use of a deep learning based approach reduces the complexity caused by the use of hand-designed features. Our method uses an off-policy reinforcement learning framework to learn the grasping policy. We use the double deep Q-learning framework along with a novel Grasp-Q-Network to output grasp probabilities used to learn grasps that maximize the pick success. We propose a visual servoing mechanism that uses a multi-view camera setup that observes the scene which contains the objects of interest. We performed experiments using a Baxter Gazebo simulated environment as well as on the actual robot. The results show that our proposed method outperforms the baseline Q-learning framework and increases grasping accuracy by adapting a multi-view model in comparison to a single-view model.


title: Dynamics in cortical activity revealed by resting-state MEG rhythms

abstract: The brain may be thought of as a many-body architecture with a spatio-temporal dynamics described by neuronal structures. The oscillatory nature of brain activity allows these structures (nodes) to be described as a set of coupled oscillators forming a network where the node dynamics, and that of the network topology can be studied. Quantifying its dynamics at various scales is an issue that claims to be explored for several brain activities, e.g., activity at rest. The resting-state associates the underlying brain dynamics of healthy subjects that are not actively compromised with sensory or cognitive processes. Studying its dynamics is highly non-trivial but opens the door to understand the general principles of brain functioning. We hypothesize about how could be the spatio-temporal dynamics of cortical fluctuations for healthy subjects at resting-state. We retrieve the alphabet that reconstructs the dynamics (entropy/complexity) of magnetoencephalograpy signals. We assemble the cortical connectivity to elicit the network's dynamics. We depict an order relation between entropy/complexity for frequency bands. We unveiled that the posterior cortex conglomerates nodes with both stronger dynamics and high clustering for {\alpha} band. The existence of these order relations suggests an emergent phenomenon of each band. Interestingly, we find that the posterior cortex plays a cardinal role in both the dynamics and structure regarding the resting-state. To the best of our knowledge, this is the first study with magnetoencephalograpy involving information theory and network science to better understand the dynamics and structure of brain activity at rest for different bands and scales.


title: Non-uniform dependence on initial data for the Camassa-Holm equation in the critical Besov space

abstract: Whether or not the data-to-solution map of the Cauchy problem for the Camassa-Holm equation and Novikov equation in the critical Besov space $B_{2,1}^{3/2}(\R)$ is not uniformly continuous remains open. In the paper, we aim at solving the open question left the previous works in \cite{Li3,Li4} and give a positive answer to this problem.


title: Detecting Band Profiles of Devices with Conductive Atomic Force Microscopy

abstract: Band profiles of electronic devices are of fundamental importance in determining their properties. A technique that can map the band profile of both the interior and edges of a device at the nanometer scale is highly demanded. Conventional scanning tunneling spectroscopy (STS) can map band structure at the atomic scale, but is limited to the interior of large and conductive samples. Here we develop a contact-mode STS based on conductive atomic force microscope that can remove these constraints. With this technique, we map the band profile of MoS$_2$ transistors with nanometer resolution at room temperature. A band bending of 0.6 eV within 18 nm of the edges of MoS$_2$ on insulating substrate is discovered. This technique will be of great use for both fundamental and applied studies of various electronic devices.


title: Efficient Trajectory Compression and Range Query Processing

abstract: Nowadays, there are ubiquitousness of GPS sensors in various devices collecting, transmitting and storing tremendous trajectory data. However, such an unprecedented scale of GPS data has posed an urgent demand for not only an effective storage mechanism but also an efficient query mechanism. Line simplification in online mode, searving as a mainstream trajectory compression method, plays an important role to attack this issue. But for the existing algorithms, either their time cost is extremely high, or the accuracy loss after the compression is completely unacceptable. To attack this issue, we propose $\epsilon \_$Region based Online trajectory Compression with Error bounded (ROCE for short), which makes the best balance among the accuracy loss, the time cost and the compression rate. The range query serves as a primitive, yet quite essential operation on analyzing trajectories. Each trajectory is usually seen as a sequence of discrete points, and in most previous work, a trajectory is judged to be overlapped with the query region R iff there is at least one point in this trajectory falling in R. But this traditional criteria is not suitable when the queried trajectories are compressed, because there may be hundreds of points discarded between each two adjacent points and the points in each compressed trajectory are quite sparse. And many trajectories could be missing in the result set. To address this, in this paper, a new criteria based on the probability and an efficient Range Query processing algorithm on Compressed trajectories RQC are proposed. In addition, an efficient index \emph{ASP\_tree} and lots of novel techniques are also presented to accelerate the processing of trajectory compression and range queries obviously. Extensive experiments have been done on multiple real datasets, and the results demonstrate superior performance of our methods.


title: Learning Differential Equations that are Easy to Solve

abstract: Differential equations parameterized by neural networks become expensive to solve numerically as training progresses. We propose a remedy that encourages learned dynamics to be easier to solve. Specifically, we introduce a differentiable surrogate for the time cost of standard numerical solvers, using higher-order derivatives of solution trajectories. These derivatives are efficient to compute with Taylor-mode automatic differentiation. Optimizing this additional objective trades model performance against the time cost of solving the learned dynamics. We demonstrate our approach by training substantially faster, while nearly as accurate, models in supervised classification, density estimation, and time-series modelling tasks.


title: Towards Unsupervised Learning for Instrument Segmentation in Robotic Surgery with Cycle-Consistent Adversarial Networks

abstract: Surgical tool segmentation in endoscopic images is an important problem: it is a crucial step towards full instrument pose estimation and it is used for integration of pre- and intra-operative images into the endoscopic view. While many recent approaches based on convolutional neural networks have shown great results, a key barrier to progress lies in the acquisition of a large number of manually-annotated images which is necessary for an algorithm to generalize and work well in diverse surgical scenarios. Unlike the surgical image data itself, annotations are difficult to acquire and may be of variable quality. On the other hand, synthetic annotations can be automatically generated by using forward kinematic model of the robot and CAD models of tools by projecting them onto an image plane. Unfortunately, this model is very inaccurate and cannot be used for supervised learning of image segmentation models. Since generated annotations will not directly correspond to endoscopic images due to errors, we formulate the problem as an unpaired image-to-image translation where the goal is to learn the mapping between an input endoscopic image and a corresponding annotation using an adversarial model. Our approach allows to train image segmentation models without the need to acquire expensive annotations and can potentially exploit large unlabeled endoscopic image collection outside the annotated distributions of image/annotation data. We test our proposed method on Endovis 2017 challenge dataset and show that it is competitive with supervised segmentation methods.


title: Bayesian Neural Networks for Fast SUSY Predictions

abstract: One of the goals of current particle physics research is to obtain evidence for new physics, that is, physics beyond the Standard Model (BSM), at accelerators such as the Large Hadron Collider (LHC) at CERN. The searches for new physics are often guided by BSM theories that depend on many unknown parameters, which, in some cases, makes testing their predictions difficult. In this paper, machine learning is used to model the mapping from the parameter space of the phenomenological Minimal Supersymmetric Standard Model (pMSSM), a BSM theory with 19 free parameters, to some of its predictions. Bayesian neural networks are used to predict cross sections for arbitrary pMSSM parameter points, the mass of the associated lightest neutral Higgs boson, and the theoretical viability of the parameter points. All three quantities are modeled with average percent errors of 3.34% or less and in a time significantly shorter than is possible with the supersymmetry codes from which the results are derived. These results are a further demonstration of the potential for machine learning to model accurately the mapping from the high dimensional spaces of BSM theories to their predictions.


title: Constraints and cosmography of $\Lambda$CDM in presence of viscosity

abstract: In this work, we study two scenarios of the Universe filled by a perfect fluid following the traditional dark energy and a viscous fluid as dark matter. In this sense, we explore the simplest case for the viscosity in the Eckart formalism, a constant, and then, a polynomial function of the redshift. We constrain the phase-space of the model parameters by performing a Bayesian analysis based on Markov Chain Monte Carlo method and using the latest data of the Hubble parameter (OHD), Type Ia Supernovae (SNIa) and Strong Lensing Systems. The first two samples cover the region $0.01<z<2.36$. Based on AIC, we find equally support of these viscous models over Lambda-Cold Dark Matter (LCDM) taking into account OHD or SNIa. On the other hand, we reconstruct the cosmographic parameters ($q,j,s,l$) and find good agreement to LCDM within up to $3\sigma$ CL. Additionally, we find that the cosmographic parameters and the acceleration-deceleration transition are sensible to the parameters related to the viscosity coefficient, making of the viscosity an interesting physical mechanism to modified them.


title: Cultural Cartography with Word Embeddings

abstract: Using the frequency of keywords is a classic approach in the formal analysis of text, but has the drawback of glossing over the relationality of word meanings. Word embedding models overcome this problem by constructing a standardized and continuous "meaning space" where words are assigned a location based on relations of similarity to other words based on how they are used in natural language samples. We show how word embeddings are commensurate with prevailing theories of meaning in sociology and can be put to the task of interpretation via two kinds of navigation. First, one can hold terms constant and measure how the embedding space moves around them--much like astronomers measured the changing of celestial bodies with the seasons. Second, one can also hold the embedding space constant and see how documents or authors move relative to it--just as ships use the stars on a given night to determine their location. Using the empirical case of immigration discourse in the United States, we demonstrate the merits of these two broad strategies for advancing important topics in cultural theory, including social marking, media fields, echo chambers, and cultural diffusion and change more broadly.


title: Supervised Robust Profile Clustering

abstract: In many studies, dimension reduction methods are used to profile participant characteristics. For example, nutrition epidemiologists often use latent class models to characterize dietary patterns. One challenge with such approaches is understanding subtle variations in patterns across subpopulations. Robust Profile Clustering (RPC) provides a dual flexible clustering model, where participants may cluster at two levels: (1) globally, where participants are clustered according to behaviors shared across an overall population, and (2) locally, where individual behaviors can deviate and cluster in subpopulations. We link clusters to a health outcome using a joint model. This model is used to derive dietary patterns in the United States and evaluate case proportion of orofacial clefts. Using dietary consumption data from the 1997-2009 National Birth Defects Prevention Study, a population-based case-control study, we determine how maternal dietary profiles are associated with an orofacial cleft among offspring. Results indicated that mothers who consumed a high proportion of fruits and vegetables compared to meats, such as chicken and beef, had lower odds delivering a child with an orofacial cleft defect.


title: Symbolic Reachability Analysis of High Dimensional Max-Plus Linear Systems

abstract: This work discusses the reachability analysis (RA) of Max-Plus Linear (MPL) systems, a class of continuous-space, discrete-event models defined over the max-plus algebra. Given the initial and target sets, we develop algorithms to verify whether there exist trajectories of the MPL system that, starting from the initial set, eventually reach the target set. We show that RA can be solved symbolically by encoding the MPL system, as well as initial and target sets into difference logic, and then checking the satisfaction of the resulting logical formula via an off-the-shelf satisfiability modulo theories (SMT) solver. The performance and scalability of the developed SMT-based algorithms are shown to clearly outperform state-of-the-art RA algorithms for MPL systems, newly allowing to investigate RA of high-dimensional MPL systems: the verification of models with more than 100 continuous variables shows the applicability of these techniques to MPL systems of industrial relevance.


title: Causal Effects in Twin Studies: the Role of Interference

abstract: The use of twins designs to address causal questions is becoming increasingly popular. A standard assumption is that there is no interference between twins---that is, no twin's exposure has a causal impact on their co-twin's outcome. However, there may be settings in which this assumption would not hold, and this would (1) impact the causal interpretation of parameters obtained by commonly used existing methods; (2) change which effects are of greatest interest; and (3) impact the conditions under which we may estimate these effects. We explore these issues, and we derive semi-parametric efficient estimators for causal effects in the presence of interference between twins. Using data from the Minnesota Twin Family Study, we apply our estimators to assess whether twins' consumption of alcohol in early adolescence may have a causal impact on their co-twins' substance use later in life.


title: Electronic and transport properties of anisotropic semiconductor quantum wires

abstract: Within the effective-mass approximation, we theoretically investigated the electronic and transport properties of 2D semiconductor quantum wires (QWs) with anisotropic effective masses and different orientations with respect to the anisotropic axis. The energy levels in the absence and presence of an external magnetic field are analytically calculated, showing: (i) a strong dependence on the spacing of energy levels related to the alignment QW angle and the anisotropy axis; and (ii) for non-null magnetic field, the quantum Hall edge states are significantly affected by the edge orientation. Moreover, by means of the split-operator technique, we analyzed the time evolution of wavepackets in straight and V-shaped anisotropic QWs and compared the transmission probabilities with those of isotropic systems. In the anisotropic case we found damped oscillations in the average values of velocity in both x and y directions for a symmetric Gaussian wavepacket propagating along a straight wide QW, with the oscillation being more evident as the non-collinearity between the group velocity and momentum vectors increases.


title: Computing the Largest Bond and the Maximum Connected Cut of a Graph

abstract: The cut-set $\partial(S)$ of a graph $G=(V,E)$ is the set of edges that have one endpoint in $S\subset V$ and the other endpoint in $V\setminus S$, and whenever $G[S]$ is connected, the cut $[S,V\setminus S]$ of $G$ is called a connected cut. A bond of a graph $G$ is an inclusion-wise minimal disconnecting set of $G$, i.e., bonds are cut-sets that determine cuts $[S,V\setminus S]$ of $G$ such that $G[S]$ and $G[V\setminus S]$ are both connected. Contrasting with a large number of studies related to maximum cuts, there exist very few results regarding the largest bond of general graphs. In this paper, we aim to reduce this gap on the complexity of computing the largest bond, and the maximum connected cut of a graph. Although cuts and bonds are similar, we remark that computing the largest bond and the maximum connected cut of a graph tends to be harder than computing its maximum cut. We show that it does not exist a constant-factor approximation algorithm to compute the largest bond, unless P = NP. Also, we show that {\sc Largest Bond} and {\sc Maximum Connected Cut} are NP-hard even for planar bipartite graphs, whereas \textsc{Maximum Cut} is trivial on bipartite graphs and polynomial-time solvable on planar graphs. In addition, we show that {\sc Largest Bond} and {\sc Maximum Connected Cut} are NP-hard on split graphs, and restricted to graphs of clique-width $w$ they can not be solved in time $f(w)\times n^{o(w)}$ unless the Exponential Time Hypothesis fails, but they can be solved in time $f(w)\times n^{O(w)}$. Finally, we show that both problems are fixed-parameter tractable when parameterized by the size of the solution, the treewidth, and the twin-cover number.


title: Deep Multi-task Learning for Facial Expression Recognition and Synthesis Based on Selective Feature Sharing

abstract: Multi-task learning is an effective learning strategy for deep-learning-based facial expression recognition tasks. However, most existing methods take into limited consideration the feature selection, when transferring information between different tasks, which may lead to task interference when training the multi-task networks. To address this problem, we propose a novel selective feature-sharing method, and establish a multi-task network for facial expression recognition and facial expression synthesis. The proposed method can effectively transfer beneficial features between different tasks, while filtering out useless and harmful information. Moreover, we employ the facial expression synthesis task to enlarge and balance the training dataset to further enhance the generalization ability of the proposed method. Experimental results show that the proposed method achieves state-of-the-art performance on those commonly used facial expression recognition benchmarks, which makes it a potential solution to real-world facial expression recognition problems.


title: Aligning Videos in Space and Time

abstract: In this paper, we focus on the task of extracting visual correspondences across videos. Given a query video clip from an action class, we aim to align it with training videos in space and time. Obtaining training data for such a fine-grained alignment task is challenging and often ambiguous. Hence, we propose a novel alignment procedure that learns such correspondence in space and time via cross video cycle-consistency. During training, given a pair of videos, we compute cycles that connect patches in a given frame in the first video by matching through frames in the second video. Cycles that connect overlapping patches together are encouraged to score higher than cycles that connect non-overlapping patches. Our experiments on the Penn Action and Pouring datasets demonstrate that the proposed method can successfully learn to correspond semantically similar patches across videos, and learns representations that are sensitive to object and action states.


title: Characterizations of the sphere by means of visual cones: an alternative proof of Matsuura's theorem

abstract: In this work we prove the following: let $K$ be a convex body in the Euclidean space $\mathbb{R}^n$, $n\geq 3$, contained in the interior of the unit ball of $\mathbb{R}^n$, and let $p\in \mathbb{R}^n$ be a point such that, from each point of $\mathbb{S}^{n-1}$, $K$ looks centrally symmetric and $p$ appears as the center, then $K$ is a ball.


title: Distributed Energy Trading and Scheduling among Microgrids via Multiagent Reinforcement Learning

abstract: The development of renewable energy generation empowers microgrids to generate electricity to supply itself and to trade the surplus on energy markets. To minimize the overall cost, a microgrid must determine how to schedule its energy resources and electrical loads and how to trade with others. The control decisions are influenced by various factors, such as energy storage, renewable energy yield, electrical load, and competition from other microgrids. Making the optimal control decision is challenging, due to the complexity of the interconnected microgrids, the uncertainty of renewable energy generation and consumption, and the interplay among microgrids. The previous works mainly adopted the modeling-based approaches for deriving the control decision, yet they relied on the precise information of future system dynamics, which can be hard to obtain in a complex environment. This work provides a new perspective of obtaining the optimal control policy for distributed energy trading and scheduling by directly interacting with the environment, and proposes a multiagent deep reinforcement learning approach for learning the optimal control policy. Each microgrid is modeled as an agent, and different agents learn collaboratively for maximizing their rewards. The agent of each microgrid can make the local scheduling decision without knowing others' information, which can well maintain the autonomy of each microgrid. We evaluate the performances of our proposed method using real-world datasets. The experimental results show that our method can significantly reduce the cost of the microgrids compared with the baseline methods.


title: Robust Geodesic Regression

abstract: This paper studies robust regression for data on Riemannian manifolds. Geodesic regression is the generalization of linear regression to a setting with a manifold-valued dependent variable and one or more real-valued independent variables. The existing work on geodesic regression uses the sum-of-squared errors to find the solution, but as in the classical Euclidean case, the least-squares method is highly sensitive to outliers. In this paper, we use M-type estimators, including the $L_1$, Huber and Tukey biweight estimators, to perform robust geodesic regression, and describe how to calculate the tuning parameters for the latter two. We also show that, on compact symmetric spaces, all M-type estimators are maximum likelihood estimators, and argue for the overall superiority of the $L_1$ estimator over the $L_2$ and Huber estimators on high-dimensional manifolds and over the Tukey biweight estimator on compact high-dimensional manifolds. Results from numerical examples, including analysis of real neuroimaging data, demonstrate the promising empirical properties of the proposed approach.


title: Broadband X-ray characteristics of the transient pulsar GRO J2058+42

abstract: The Be X-ray binary GRO J2058+42 recently went through a Type-II outburst during March-April 2019 lasting for about 50 days. This outburst was detected with the operating all sky X-ray monitors like the Fermi-GBM, Swift-BAT and MAXI-GSC. Two NuSTAR observations were also made, one during the rise and other during the decay of the outburst. It gave us the unique opportunity to analyze the broadband characteristics of the pulsar for the first time and accretion torque characteristics of the pulsar over a range of X-ray luminosity. The pulse profiles are strongly energy dependent, with at least four different pulse components at low energy (< 20 keV) which evolves to a single-peaked profile at high energy (> 30 keV). In each of the narrow energy bands, the pulse profiles are nearly identical in the two NuSTAR observations. The spectra from both the observations are fitted well to a power-law with a Fermi-Dirac type high energy cutoff. We ruled out presence of a cyclotron line in the pulse phase averaged X-ray spectrum in the NuSTAR band with an optical depth greater than 0.15. An iron emission line is detected in both the NuSTAR spectra with an equivalent width of about 125 eV. We looked at the dependence of the spin-up rate on the luminosity and estimated the magnetic field strength from that, which came out to be much higher compared to other known BeXRB pulsars. Lastly, we discussed the inadequacy of the torque-luminosity relation for determination of magnetic field strength of neutron stars.