title: Physics of Prodigious Lyman Continuum Leakers

abstract: An analysis of the dynamics of a star formation event is performed. It is shown that galaxies able to drive leftover gas to sufficient altitudes in a few million years are characterized by two basic properties: small sizes (<1kpc) and high star formation rate surface densities (Sigma_SFR > 10 Msun/yr/kpc2). For the parameter space of relevance, the outflow is primarily driven by supernovae with radiation pressure being significant but subdominant. Our analysis provides the unifying physical origin for a diverse set of observed LyC leakers, including the green-peas galaxies, [SII]-weak galaxies, Lyman-alpha emitters, with these two characteristics as the common denominator. Among verifiable physical properties of LyC leakers, we predict that (1) the newly formed stellar masses are are typically in the range of 1e8-1e10 Msun, except perhaps ULIRGs, (2) the outflow velocities are typically in the range typically of 100-600km/s, but may exceed 1e3 km/s in ULIRGs, with a strong positive correlation between the stellar masses formed and the outflow velocities, (3) the overall escape fraction of galaxies is expected to increase with increasing redshift, given the cosmological trend that galaxies become denser and more compact with increasing redshift. In addition, two interesting by-product predictions are also borne out. First, ULIRGs appear to be in a parameter region where they should be prodigious LyC leakers, unless there is a large ram-pressure. Second, Lyman break galaxies (LBGs) are not supposed to be prodigious LyC leakers in our model, given their claimed effective radii exceeding 1kpc.


title: Hypertoric Hitchin systems and Kirchhoff polynomials

abstract: We define a formal algebraic analogue of hypertoric Hitchin systems, whose complex-analytic counterparts were defined by Hausel-Proudfoot. These are algebraic completely integrable systems associated to a graph. We study the variation of the Tamagawa number of the resulting family of abelian varieties, and show that it is described by the Kirchhoff polynomial of the graph. In particular, this allows us to compute their p-adic volumes. We conclude the article by remarking that these spaces admit a volume preserving tropicalisation.


title: Deep Channel Learning For Large Intelligent Surfaces Aided mm-Wave Massive MIMO Systems

abstract: This letter presents the first work introducing a deep learning (DL) framework for channel estimation in large intelligent surface (LIS) assisted massive MIMO (multiple-input multiple-output) systems. A twin convolutional neural network (CNN) architecture is designed and it is fed with the received pilot signals to estimate both direct and cascaded channels. In a multi-user scenario, each user has access to the CNN to estimate its own channel. The performance of the proposed DL approach is evaluated and compared with state-of-the-art DL-based techniques and its superior performance is demonstrated.


title: Physics-Guided Machine Learning for Scientific Discovery: An Application in Simulating Lake Temperature Profiles

abstract: Physics-based models of dynamical systems are often used to study engineering and environmental systems. Despite their extensive use, these models have several well-known limitations due to simplified representations of the physical processes being modeled or challenges in selecting appropriate parameters. While-state-of-the-art machine learning models can sometimes outperform physics-based models given ample amount of training data, they can produce results that are physically inconsistent. This paper proposes a physics-guided recurrent neural network model (PGRNN) that combines RNNs and physics-based models to leverage their complementary strengths and improves the modeling of physical processes. Specifically, we show that a PGRNN can improve prediction accuracy over that of physics-based models, while generating outputs consistent with physical laws. An important aspect of our PGRNN approach lies in its ability to incorporate the knowledge encoded in physics-based models. This allows training the PGRNN model using very few true observed data while also ensuring high prediction accuracy. Although we present and evaluate this methodology in the context of modeling the dynamics of temperature in lakes, it is applicable more widely to a range of scientific and engineering disciplines where physics-based (also known as mechanistic) models are used, e.g., climate science, materials science, computational chemistry, and biomedicine.


title: Crash Themes in Automated Vehicles: A Topic Modeling Analysis of the California Department of Motor Vehicles Automated Vehicle Crash Database

abstract: Automated vehicle technology promises to reduce the societal impact of traffic crashes. Early investigations of this technology suggest that significant safety issues remain during control transfers between the automation and human drivers and automation interactions with the transportation system. In order to address these issues, it is critical to understand both the behavior of human drivers during these events and the environments where they occur. This article analyzes automated vehicle crash narratives from the California Department of Motor Vehicles automated vehicle crash database to identify safety concerns and gaps between crash types and current areas of focus in the current research. The database was analyzed using probabilistic topic modeling of open-ended crash narratives. Topic modeling analysis identified five themes in the database: driver-initiated transition crashes, sideswipe crashes during left-side overtakes, and rear-end collisions while the vehicle was stopped at an intersection, in a turn lane, and when the crash involved oncoming traffic. Many crashes represented by the driver-initiated transitions topic were also associated with the side-swipe collisions. A substantial portion of the side-swipe collisions also involved motorcycles. These findings highlight previously raised safety concerns with transitions of control and interactions between vehicles in automated mode and the transportation social network. In response to these findings, future empirical work should focus on driver-initiated transitions, overtakes, silent failures, complex traffic situations, and adverse driving environments. Beyond this future work, the topic modeling analysis method may be used as a tool to monitor emergent safety issues.


title: Low power switching of magnetization using enhanced magnetic anisotropy with short-voltage-pulse application

abstract: A low power magnetization switching scheme based on the voltage control of magnetic anisotropy (VCMA) is proposed. In contrast to the conventional switching scheme using VCMA, where the magnetic anisotropy is eliminated during the voltage pulse, the magnetic anisotropy is enhanced to induce precession around the axis close to the easy axis. After turning off the voltage proximately at a half of precession period, the magnetization relaxes to the opposite equilibrium direction. We perform numerical simulations and show that the pulse duration of the proposed switching scheme is as short as a few tens of pico seconds. Such a short pulse duration is beneficial for low power consumption because of the reduction of energy loss by Joule heating.


title: Further Combinatorics and Applications of Two-Toned Tilings

abstract: Integer compositions, integer partitions, Fibonacci numbers, and generalizations of these have recently been shown to be interconnected via two-toned tilings of horizontal grids. In this article, we present refinements of two-toned tilings, describe functions which analyze them, and apply these to generalizations of integer compositions and partitions which interpolate between the two.


title: An investigation of global radial basis function collocation methods applied to Helmholtz problems

abstract: Global radial basis function (RBF) collocation methods with inifinitely smooth basis functions for partial differential equations (PDEs) work in general geometries, and can have exponential convergence properties for smooth solution functions. At the same time, the linear systems that arise are dense and severely ill-conditioned for large numbers of unknowns and small values of the shape parameter that determines how flat the basis functions are. We use Helmholtz equation as an application problem for the theoretical analysis and numerical experiments. We analyse and characterise the convergence properties as a function of the number of unknowns and for different shape parameter ranges. We provide theoretical results for the flat limit of the PDE solutions and investigate when the non-symmetric collocation matrices become singular. We also provide practical strategies for choosing the method parameters and evaluate the results on Helmholtz problems in a curved waveguide geometry.


title: The benefits of synthetic data for action categorization

abstract: In this paper, we study the value of using synthetically produced videos as training data for neural networks used for action categorization. Motivated by the fact that texture and background of a video play little to no significant roles in optical flow, we generated simplified texture-less and background-less videos and utilized the synthetic data to train a Temporal Segment Network (TSN). The results demonstrated that augmenting TSN with simplified synthetic data improved the original network accuracy (68.5%), achieving 71.8% on HMDB-51 when adding 4,000 videos and 72.4% when adding 8,000 videos. Also, training using simplified synthetic videos alone on 25 classes of UCF-101 achieved 30.71% when trained on 2500 videos and 52.7% when trained on 5000 videos. Finally, results showed that when reducing the number of real videos of UCF-25 to 10% and combining them with synthetic videos, the accuracy drops to only 85.41%, compared to a drop to 77.4% when no synthetic data is added.


title: FIS-Nets: Full-image Supervised Networks for Monocular Depth Estimation

abstract: This paper addresses the importance of full-image supervision for monocular depth estimation. We propose a semi-supervised architecture, which combines both unsupervised framework of using image consistency and supervised framework of dense depth completion. The latter provides full-image depth as supervision for the former. Ego-motion from navigation system is also embedded into the unsupervised framework as output supervision of an inner temporal transform network, making monocular depth estimation better. In the evaluation, we show that our proposed model outperforms other approaches on depth estimation.


title: SLO-ML: A Language for Service Level Objective Modelling in Multi-cloud Applications

abstract: Cloud modelling languages (CMLs) are designed to assist customers in tackling the diversity of services in the cloud market. While many CMLs have been proposed in the literature, they lack practical support for automating the selection of services based on the specific service level objectives of a customer's application. We put forward SLO-ML, a novel and generative CML to capture service level requirements and, subsequently, to select the services to honour customer requirements and generate the deployment code appropriate to these services. We present the architectural design of SLO-ML and the associated broker that realises the deployment operations. We rigorously evaluate SLO-ML using a mixed methods approach. First, we exploit an experimental case study with a group of researchers and developers using a real-world cloud application. We also assess overheads through an exhaustive set of empirical scalability tests. Through expressing the levels of gained productivity and experienced usability, we highlight SLO-ML's profound potential in enabling user-centric cloud brokers. We also discuss limitations as application requirements grow.


title: Inflationary predictions of Geometric Inflation

abstract: In the framework of gravitational models obtained from the Geometric Inflation's proposal, where an infinite tower of curvature scalars are included into the action, we compute the slow-roll parameters by the Hubble slow-roll approach. We test the viability of such models as inflationary scenarios, focusing on the tensor-to-scalar ratio, $r$, and spectral scalar index, $n_s$, relation. We find that all models considered here produce inflation and, most of them coincide, some better than others, with the marginalized 95$\%$ CL region given by Planck's data collaboration.


title: Doubly periodic lozenge tilings of a hexagon and matrix valued orthogonal polynomials

abstract: We analyze a random lozenge tiling model of a large regular hexagon, whose underlying weight structure is periodic of period $2$ in both the horizontal and vertical directions. This is a determinantal point process whose correlation kernel is expressed in terms of non-Hermitian matrix valued orthogonal polynomials. This model belongs to a class of models for which the existing techniques for studying asymptotics cannot be applied. The novel part of our method consists of establishing a connection between matrix valued and scalar valued orthogonal polynomials. This allows to simplify the double contour formula for the kernel obtained by Duits and Kuijlaars by reducing the size of a Riemann-Hilbert problem. The proof relies on the fact that the matrix valued weight possesses eigenvalues that live on an underlying Riemann surface $\mathcal{M}$ of genus $0$. We consider this connection of independent interest; it is natural to expect that similar ideas can be used for other matrix valued orthogonal polynomials, as long as the corresponding Riemann surface $\mathcal{M}$ is of genus $0$. The rest of the method consists of two parts, and mainly follows the lines of a previous work of Charlier, Duits, Kuijlaars and Lenells. First, we perform a Deift-Zhou steepest descent analysis to obtain asymptotics for the scalar valued orthogonal polynomials. The main difficulty is the study of an equilibrium problem in the complex plane. Second, the asymptotics for the orthogonal polynomials are substituted in the double contour integral and the latter is analyzed using the saddle point method. Our main results are the limiting densities of the lozenges in the disordered flower-shaped region. However, we stress that the method allows in principle to rigorously compute other meaningful probabilistic quantities in the model.


title: Codimension-$1$ Simplices in Divisible Convex Domains

abstract: Properly embedded simplices in a convex divisible domain $\Omega \subset \mathbb{R} \textrm{P}^d$ behave somewhat like flats in Riemannian manifolds, so we call them flats. We show that the set of codimension-$1$ flats has image which is a finite collection of disjoint virtual $(d-1)$-tori in the compact quotient manifold. If this collection of virtual tori is non-empty, then the components of its complement are cusped convex projective manifolds with type $d$ cusps.


title: Plectic Galois action on CM points and connected components of Hilbert modular varieties

abstract: We expand on Nekov\'a\v{r}'s construction of the plectic half transfer to define a plectic Galois action on Hilbert modular varieties. More precisely, we study in a unifying fashion Shimura varieties associated to groups that differ only in the centre from $R_{F/\mathbb{Q}}{\rm GL}_2$. We define plectic Galois actions on the CM points and on the set of connected components of these Shimura varieties, and show that these two actions are compatible. This extends the plectic conjecture of Nekov\'a\v{r}--Scholl.


title: Solution of the logarithmic coefficients conjecture in some families of univalent functions

abstract: For univalent and normalized functions $f$ the logarithmic coefficients $\gamma_n(f)$ are determined by the formula $\log(f(z)/z)=\sum_{n=1}^{\infty}2\gamma_n(f)z^n$. In the paper \cite{Pon} the authors posed the conjecture that a locally univalent function in the unit disk, satisfying the condition \[ \Re\left\{1+zf''(z)/f'(z)\right\}<1+\lambda/2\quad (z\in \mathbb{D}), \] fulfill also the following inequality: $$|\gamma_n(f)|\le \lambda/(2n(n+1)).$$ Here $\lambda$ is a real number such that $0<\lambda\le 1$. In the paper we confirm that the conjecture is true, and sharp.


title: Motivation to Run in One-Day Cricket

abstract: In this paper we introduce a new methodology to determine an optimal coefficient for a positive finite measure of batting average, strike rate, and bowling average of a player in order to get an optimal score of a team under dynamic modeling using a path integral method. We also introduce new run dynamics modeled as a stochastic differential equation in order to incorporate the average weather conditions at the cricket ground, the weather condition on the day of the match including sudden deterioration which leads to a partial or complete stop of the game, total attendance, and home field advantage.


title: A Scalable Framework for Quality Assessment of RDF Datasets

abstract: Over the last years, Linked Data has grown continuously. Today, we count more than 10,000 datasets being available online following Linked Data standards. These standards allow data to be machine readable and inter-operable. Nevertheless, many applications, such as data integration, search, and interlinking, cannot take full advantage of Linked Data if it is of low quality. There exist a few approaches for the quality assessment of Linked Data, but their performance degrades with the increase in data size and quickly grows beyond the capabilities of a single machine. In this paper, we present DistQualityAssessment -- an open source implementation of quality assessment of large RDF datasets that can scale out to a cluster of machines. This is the first distributed, in-memory approach for computing different quality metrics for large RDF datasets using Apache Spark. We also provide a quality assessment pattern that can be used to generate new scalable metrics that can be applied to big data. The work presented here is integrated with the SANSA framework and has been applied to at least three use cases beyond the SANSA community. The results show that our approach is more generic, efficient, and scalable as compared to previously proposed approaches.


title: Urban2Vec: Incorporating Street View Imagery and POIs for Multi-Modal Urban Neighborhood Embedding

abstract: Understanding intrinsic patterns and predicting spatiotemporal characteristics of cities require a comprehensive representation of urban neighborhoods. Existing works relied on either inter- or intra-region connectivities to generate neighborhood representations but failed to fully utilize the informative yet heterogeneous data within neighborhoods. In this work, we propose Urban2Vec, an unsupervised multi-modal framework which incorporates both street view imagery and point-of-interest (POI) data to learn neighborhood embeddings. Specifically, we use a convolutional neural network to extract visual features from street view images while preserving geospatial similarity. Furthermore, we model each POI as a bag-of-words containing its category, rating, and review information. Analog to document embedding in natural language processing, we establish the semantic similarity between neighborhood ("document") and the words from its surrounding POIs in the vector space. By jointly encoding visual, textual, and geospatial information into the neighborhood representation, Urban2Vec can achieve performances better than baseline models and comparable to fully-supervised methods in downstream prediction tasks. Extensive experiments on three U.S. metropolitan areas also demonstrate the model interpretability, generalization capability, and its value in neighborhood similarity analysis.


title: Ensemble Grammar Induction For Detecting Anomalies in Time Series

abstract: Time series anomaly detection is an important task, with applications in a broad variety of domains. Many approaches have been proposed in recent years, but often they require that the length of the anomalies be known in advance and provided as an input parameter. This limits the practicality of the algorithms, as such information is often unknown in advance, or anomalies with different lengths might co-exist in the data. To address this limitation, previously, a linear time anomaly detection algorithm based on grammar induction has been proposed. While the algorithm can find variable-length patterns, it still requires preselecting values for at least two parameters at the discretization step. How to choose these parameter values properly is still an open problem. In this paper, we introduce a grammar-induction-based anomaly detection method utilizing ensemble learning. Instead of using a particular choice of parameter values for anomaly detection, the method generates the final result based on a set of results obtained using different parameter values. We demonstrate that the proposed ensemble approach can outperform existing grammar-induction-based approaches with different criteria for selection of parameter values. We also show that the proposed approach can achieve performance similar to that of the state-of-the-art distance-based anomaly detection algorithm.


title: Simulation of electron-proton scattering events by a Feature-Augmented and Transformed Generative Adversarial Network (FAT-GAN)

abstract: We apply generative adversarial network (GAN) technology to build an event generator that simulates particle production in electron-proton scattering that is free of theoretical assumptions about underlying particle dynamics. The difficulty of efficiently training a GAN event simulator lies in learning the complicated patterns of the distributions of the particles physical properties. We develop a GAN that selects a set of transformed features from particle momenta that can be generated easily by the generator, and uses these to produce a set of augmented features that improve the sensitivity of the discriminator. The new Feature-Augmented and Transformed GAN (FAT-GAN) is able to faithfully reproduce the distribution of final state electron momenta in inclusive electron scattering, without the need for input derived from domain-based theoretical assumptions. The developed technology can play a significant role in boosting the science of existing and future accelerator facilities, such as the Electron-Ion Collider.


title: A practical guide to using the AngCor package

abstract: This guide contains brief, practically oriented instructions for the use of the AngCor package for computing angular correlations following nuclear reactions. This guide serves as a lasting reference for the AngCor package and complements the GitHub repository for the AngCor code.


title: Exploring the dispersion measure of the Milky Way halo

abstract: Fast radio bursts offer the opportunity to place new constraints on the mass and density profile of hot and ionized gas in galactic haloes. We test here the X-ray emission and dispersion measure predicted by different gas profiles for the halo of the Milky Way. We examine a range of models, including entropy stability conditions and external pressure continuity. We find that incorporating constraints from X-ray observations leads to favouring dispersion measures on the lower end of the range given by these models. We show that the dispersion measure of the Milky Way halo could be less than 10 cm$^{-3}$ pc in the most extreme model we consider, which is based on constraints from O VII absorption lines. However, the models allowed by the soft X-ray constraints span more than an order of magnitude in dispersion measures. Additional information on the distribution of gas in the Milky Way halo could be obtained from the signature of a dipole in the dispersion measure of fast radio bursts across the sky, but this will be a small effect for most cases.


title: The order of the product of two elements in finite nilpotent groups

abstract: An old problem in group theory is that of describing how the order of an element behaves under multiplication. To generalize some classical bounds concerning the order $\mathrm o(ab)$ of two elements $a, b$ in a finite abelian group to the non-commutative case, we replace $\mathrm o(ab)$ with a notion of mutual order $\mathrm o(a, b)$, defined as the least positive integer $n$ such that $a^nb^n = 1$. Motivated by this, we then compare $\mathrm o(ab)$ and $\mathrm o(a, b)$ in finite nilpotent groups, and show that in a group of class $\gamma$, the ratio $\mathrm o(ab)/\mathrm o(a, b)$ lies in some fixed finite set $S(\gamma) \subset \mathbb Q$, whose elements do not involve prime factors exceeding $\gamma$. In particular, we generalize a result of P. Hall, which asserts that $\mathrm o(ab) = \mathrm o(a, b)$ in $p$-groups with $p > \gamma$. We end with a more detailed analysis for groups of class 2, which allows one to give a more explicit description of $\mathrm o(ab)/\mathrm o(a, b)$.


title: Hamiltonian neural networks for solving equations of motion

abstract: There has been a wave of interest in applying machine learning to study dynamical systems. We present a Hamiltonian neural network that solves the differential equations that govern dynamical systems. This is an equation-driven machine learning method where the optimization process of the network depends solely on the predicted functions without using any ground truth data. The model learns solutions that satisfy, up to an arbitrarily small error, Hamilton's equations and, therefore, conserve the Hamiltonian invariants. The choice of an appropriate activation function drastically improves the predictability of the network. Moreover, an error analysis is derived and states that the numerical errors depend on the overall network performance. The Hamiltonian network is then employed to solve the equations for the nonlinear oscillator and the chaotic Henon-Heiles dynamical system. In both systems, a symplectic Euler integrator requires two orders more evaluation points than the Hamiltonian network in order to achieve the same order of the numerical error in the predicted phase space trajectories.


title: D2M: Dynamic Defense and Modeling of Adversarial Movement in Networks

abstract: Given a large enterprise network of devices and their authentication history (e.g., device logons), how can we quantify network vulnerability to lateral attack and identify at-risk devices? We systematically address these problems through D2M, the first framework that models lateral attacks on enterprise networks using multiple attack strategies developed with researchers, engineers, and threat hunters in the Microsoft Defender Advanced Threat Protection group. These strategies integrate real-world adversarial actions (e.g., privilege escalation) to generate attack paths: a series of compromised machines. Leveraging these attack paths and a novel Monte-Carlo method, we formulate network vulnerability as a probabilistic function of the network topology, distribution of access credentials and initial penetration point. To identify machines at risk to lateral attack, we propose a suite of five fast graph mining techniques, including a novel technique called AnomalyShield inspired by node immunization research. Using three real-world authentication graphs from Microsoft and Los Alamos National Laboratory (up to 223,399 authentications), we report the first experimental results on network vulnerability to lateral attack, demonstrating D2M's unique potential to empower IT admins to develop robust user access credential policies.


title: Singular integral equations with applications to travelling waves for doubly nonlinear diffusion

abstract: We consider a family of singular Volterra integral equations that appear in the study of monotone travelling-wave solutions for a family of diffusion-convection-reaction equations involving the $p$-Laplacian operator. Our results extend the ones due to B.\,Gilding for the case $p=2$. The fact that $p\neq2$ modifies the nature of the singularity in the integral equation, and introduces the need to develop some new tools for the analysis. The results for the integral equation are then used to study the existence and properties of travelling-wave solutions for doubly nonlinear diffusion-reaction equations in terms of the constitutive functions of the problem.


title: Electromagnetic radiation from the pre-equilibrium/pre-hydro stage of the quark-gluon plasma

abstract: Photons produced in the pre-equilibrium/pre-hydro stage of the quark-gluon plasma produced in relativistic heavy-ion collisions were computed using parton distribution functions obtained from solutions of the Boltzmann equation. The effect of the initial gluon momentum anisotropy $\xi$ and the dependence on the saturation momentum $Q_s$ was investigated. We see that small $Q_s$ results in a photon yield enhancement, whereas a larger $Q_s$ results in a pre-equilibrium photon suppression, owing to the strict constraint of matching to experimental energy density.


title: Asymptotics of Cross-Validation

abstract: Cross validation is a central tool in evaluating the performance of machine learning and statistical models. However, despite its ubiquitous role, its theoretical properties are still not well understood. We study the asymptotic properties of the cross validated-risk for a large class of models. Under stability conditions, we establish a central limit theorem and Berry-Esseen bounds, which enable us to compute asymptotically accurate confidence intervals. Using our results, we paint a big picture for the statistical speed-up of cross validation compared to a train-test split procedure. A corollary of our results is that parametric M-estimators (or empirical risk minimizers) benefit from the "full" speed-up when performing cross-validation under the training loss. In other common cases, such as when the training is done using a surrogate loss or a regularizer, we show that the behavior of the cross-validated risk is complex with a variance reduction which may be smaller or larger than the "full" speed-up, depending on the model and the underlying distribution. We allow the number of folds to grow with the number of observations at any rate.


title: TraMoS V. Updated ephemeris and multi-epoch monitoring of the hot Jupiters WASP-18Ab, WASP-19b, and WASP-77Ab

abstract: We present 22 new transit observations of the exoplanets WASP-18Ab, WASP-19b, and WASP-77Ab, from the Transit Monitoring in the South (TraMoS) project. We simultaneously model our newly collected transit light curves with archival photometry and radial velocity data, to obtain refined physical and orbital parameters. We include TESS light curves of the three exoplanets to perform an extended analysis of the variations in their transit mid-time (TTV) and to refine their planetary orbital ephemeris. We did not find significant $\rm TTV_{RMS}$ variations larger than 47, 65, and 86 seconds for WASP-18Ab, WASP-19b, and WASP-77Ab, respectively. Dynamical simulations were carried out to constrain the masses of a possible perturber. The observed RMS could be produced by a perturber body with an upper limit mass of 9, 2.5, 11 and $4~M_{\oplus}$ in 1:2, 1:3, 2:1, and 3:1 resonances in the WASP-18Ab system. In the case of WASP-19b, companions with masses up to 0.26, 0.65, 1 and $2.8~M_{\oplus}$, in 1:2, 2:1, 3:1, and 5:3 resonances, produce the RMS. And for the WASP-77Ab system, a planet with mass between $1.5-9~M_{\oplus}$ in 1:2, 1:3, 2:1, 2:3, 3:1, 3:5, 5:3 resonances. Comparing our results with RV variations, we discard massive companions with $350~M_{\oplus}$ in 17:5 resonance for WASP-18Ab, $95~M_{\oplus}$ in 4:1 resonance for WASP-19b and $105~M_{\oplus}$ in 5:2 resonance for WASP-77Ab. Finally, using a Lomb-Scargle period search we find no evidence of a periodic trend on our TTV data for the three exoplanets.


title: GradientDICE: Rethinking Generalized Offline Estimation of Stationary Values

abstract: We present GradientDICE for estimating the density ratio between the state distribution of the target policy and the sampling distribution in off-policy reinforcement learning. GradientDICE fixes several problems of GenDICE (Zhang et al., 2020), the state-of-the-art for estimating such density ratios. Namely, the optimization problem in GenDICE is not a convex-concave saddle-point problem once nonlinearity in optimization variable parameterization is introduced to ensure positivity, so any primal-dual algorithm is not guaranteed to converge or find the desired solution. However, such nonlinearity is essential to ensure the consistency of GenDICE even with a tabular representation. This is a fundamental contradiction, resulting from GenDICE's original formulation of the optimization problem. In GradientDICE, we optimize a different objective from GenDICE by using the Perron-Frobenius theorem and eliminating GenDICE's use of divergence. Consequently, nonlinearity in parameterization is not necessary for GradientDICE, which is provably convergent under linear function approximation.


title: A Family of Pairwise Multi-Marginal Optimal Transports that Define a Generalized Metric

abstract: The Optimal transport (OT) problem is rapidly finding its way into machine learning. Favoring its use are its metric properties. Many problems admit solutions with guarantees only for objects embedded in metric spaces, and the use of non-metrics can complicate solving them. Multi-marginal OT (MMOT) generalizes OT to simultaneously transporting multiple distributions. It captures important relations that are missed if the transport only involves two distributions. Research on MMOT, however, has been focused on its existence, uniqueness, practical algorithms, and the choice of cost functions. There is a lack of discussion on the metric properties of MMOT, which limits its theoretical and practical use. Here, we prove new generalized metric properties for a family of pairwise MMOTs. We first explain the difficulty of proving this via two negative results. Afterward, we prove the MMOTs' metric properties. Finally, we show that the generalized triangle inequality of this family of MMOTs cannot be improved. We illustrate the superiority of our MMOTs over other generalized metrics, and over non-metrics in both synthetic and real tasks.


title: Multichannel ALOHA with Exploration Phase

abstract: In this paper, we consider exploration for multichannel ALOHA by transmitting preambles before transmitting data packets and show that the maximum throughput can be improved by a factor of 2 - exp(-1) = 1.632, which can be seen as the gain of exploration. In the proposed approach, a base station (BS) needs to send the feedback information to active users to inform the numbers of transmitted preambles in multiple channels, which can be reliably estimated as in compressive random access. Simulation results also confirm the results from analysis.


title: Counterfactual Programming for Optimal Control

abstract: In recent years, considerable work has been done to tackle the issue of designing control laws based on observations to allow unknown dynamical systems to perform pre-specified tasks. At least as important for autonomy, however, is the issue of learning which tasks can be performed in the first place. This is particularly critical in situations where multiple (possibly conflicting) tasks and requirements are demanded from the agent, resulting in infeasible specifications. Such situations arise due to over-specification or dynamic operating conditions and are only aggravated when the dynamical system model is learned through simulations. Often, these issues are tackled using regularization and penalties tuned based on application-specific expert knowledge. Nevertheless, this solution becomes impractical for large-scale systems, unknown operating conditions, and/or in online settings where expert input would be needed during the system operation. Instead, this work enables agents to autonomously pose, tune, and solve optimal control problems by compromising between performance and specification costs. Leveraging duality theory, it puts forward a counterfactual optimization algorithm that directly determines the specification trade-off while solving the optimal control problem.


title: Summary of the 14th IACHEC Meeting

abstract: We summarize the 14th meeting of the International Astronomical Consortium for High Energy Calibration (IACHEC) held at \textit{Shonan Village} (Kanagawa, Japan) in May 2019. Sixty scientists directly involved in the calibration of operational and future high-energy missions gathered during 3.5 days to discuss the status of the cross-calibration between the current international complement of X-ray observatories, and the possibilities to improve it. This summary consists of reports from the various WGs with topics ranging from the identification and characterization of standard calibration sources, multi-observatory cross-calibration campaigns, appropriate and new statistical techniques, calibration of instruments and characterization of background, communication and preservation of knowledge, and results for the benefit of the astronomical community.


title: Self-consistent Semi-analytic Modeling of Feedback During Primordial Star Formation and Reionization

abstract: We present a new semi-analytic model of the formation of the first stars. Our method takes dark matter halo merger trees (including 3-dimensional spatial information) from cosmological N-body simulations as input and applies analytic prescriptions to compute both the Population III (Pop III) and metal-enriched star formation histories. We have developed a novel method to accurately compute the major feedback processes affecting Pop III star formation: H$_2$ photodissociation from Lyman-Werner (LW) radiation, suppression of star formation due to inhomogeneous reionization, and metal enrichment via supernovae winds. Our method utilizes a grid-based approach relying on fast Fourier transforms (FFTs) to rapidly track the LW intensity, ionization fraction, and metallicity in 3-dimensions throughout the simulation box. We present simulations for a wide range of astrophysical model parameters from $z\approx 30-6$. Initially long-range LW feedback and local metal enrichment and reionization feedback dominate. However, for $z \lesssim 15$ we find that the star formation rate density (SFRD) of Pop III stars is impacted by the combination of external metal enrichment (metals from one halo polluting other pristine halos) and inhomogeneous reionization. We find that the interplay of these processes is particularly important for the Pop III SFRD at $z \lesssim 10$. Reionization feedback delays star formation long enough for metal bubbles to reach halos that would otherwise form Pop III stars. Including these effects can lead to more than an order of magnitude decrease in the Pop III SFRD at $z=6$ compared to LW feedback alone.


title: Single Atom Qubits: Acceptors

abstract: Acceptor dopant atoms in silicon have recently been identified as compelling candidates for spin-based quantum technologies. Interest in acceptor qubits ultimately derives from the properties of acceptor bound holes, where spin-orbit coupling quantizes total angular momentum $J=3/2$ rather than spin, endowing them with quadrupolar couplings to electric and elastic fields. This property of single-atom acceptor qubits makes them amenable to electric control and coupling over long distances using phonon, capacitive, or microwave photon mediated interactions. Accepted as a contribution within a roadmap for quantum technologies, this section reviews progress on single-acceptor devices, observation of ultra-long coherence times of acceptors in isotope purified silicon, and comments on future challenges in single qubit measurement, long distance coupling, and scalable acceptor-spin-based technologies.


title: Gun Source and Muzzle Head Detection

abstract: There is a surging need across the world for protection against gun violence. There are three main areas that we have identified as challenging in research that tries to curb gun violence: temporal location of gunshots, gun type prediction and gun source (shooter) detection. Our task is gun source detection and muzzle head detection, where the muzzle head is the round opening of the firing end of the gun. We would like to locate the muzzle head of the gun in the video visually, and identify who has fired the shot. In our formulation, we turn the problem of muzzle head detection into two sub-problems of human object detection and gun smoke detection. Our assumption is that the muzzle head typically lies between the gun smoke caused by the shot and the shooter. We have interesting results both in bounding the shooter as well as detecting the gun smoke. In our experiments, we are successful in detecting the muzzle head by detecting the gun smoke and the shooter.


title: ABSent: Cross-Lingual Sentence Representation Mapping with Bidirectional GANs

abstract: A number of cross-lingual transfer learning approaches based on neural networks have been proposed for the case when large amounts of parallel text are at our disposal. However, in many real-world settings, the size of parallel annotated training data is restricted. Additionally, prior cross-lingual mapping research has mainly focused on the word level. This raises the question of whether such techniques can also be applied to effortlessly obtain cross-lingually aligned sentence representations. To this end, we propose an Adversarial Bi-directional Sentence Embedding Mapping (ABSent) framework, which learns mappings of cross-lingual sentence representations from limited quantities of parallel data.


title: Joint Visual-Temporal Embedding for Unsupervised Learning of Actions in Untrimmed Sequences

abstract: Understanding the structure of complex activities in untrimmed videos is a challenging task in the area of action recognition. One problem here is that this task usually requires a large amount of hand-annotated minute- or even hour-long video data, but annotating such data is very time consuming and can not easily be automated or scaled. To address this problem, this paper proposes an approach for the unsupervised learning of actions in untrimmed video sequences based on a joint visual-temporal embedding space. To this end, we combine a visual embedding based on a predictive U-Net architecture with a temporal continuous function. The resulting representation space allows detecting relevant action clusters based on their visual as well as their temporal appearance. The proposed method is evaluated on three standard benchmark datasets, Breakfast Actions, INRIA YouTube Instructional Videos, and 50 Salads. We show that the proposed approach is able to provide a meaningful visual and temporal embedding out of the visual cues present in contiguous video frames and is suitable for the task of unsupervised temporal segmentation of actions.


title: Tjurina number of a local complete intersection curve

abstract: Differently from the Milnor number, no formula relating the Tjurina number of a reducible algebroid curve to invariants of its branches was known. The aim of this work is to provide such a formula for a complete intersection algebroid curve with several branches in terms of invariants of its branches and the maximal points of the set of values of K\"ahler differentials on components of the curve.


title: Magnon-Phonon Quantum Correlation Thermometry

abstract: A large fraction of quantum science and technology requires low-temperature environments such as those afforded by dilution refrigerators. In these cryogenic environments, accurate thermometry can be difficult to implement, expensive, and often requires calibration to an external reference. Here, we theoretically propose a primary thermometer based on measurement of a hybrid system consisting of phonons coupled via a magnetostrictive interaction to magnons. Thermometry is based on a cross-correlation measurement in which the spectrum of back-action driven motion is used to scale the thermomechanical motion, providing a direct measurement of the phonon temperature independent of experimental parameters. Combined with a simple low-temperature compatible microwave cavity read-out, this primary thermometer is expected to become a popular thermometer for experiments below 1 K.


title: Testing biological network motif significance with exponential random graph models

abstract: Analysis of the structure of biological networks often uses statistical tests to establish the over-representation of motifs, which are thought to be important building blocks of such networks, related to their biological functions. However, there is disagreement as to the statistical significance of these motifs, and there are potential problems with standard methods for estimating this significance. Exponential random graph models (ERGMs) are a class of statistical model that can overcome some of the shortcomings of commonly used methods for testing the statistical significance of motifs. ERGMs were first introduced into the bioinformatics literature over ten years ago but have had limited application to biological networks, possibly due to the practical difficulty of estimating model parameters. Advances in estimation algorithms now afford analysis of much larger networks in practical time. We illustrate the application of ERGM to both an undirected protein-protein interaction (PPI) network and directed gene regulatory networks. ERGM models indicate over-representation of triangles in the PPI network, and confirm results from previous research as to over-representation of transitive triangles (feed-forward loop) in an E. coli and a yeast regulatory network. We also confirm, using ERGMs, previous research showing that under-representation of the cyclic triangle (feedback loop) can be explained as a consequence of other topological features.


title: Strain-induced tuning of optical properties of layered Mo$S_2$

abstract: The sensitive correlation between optical parameters and strain in Mo$S_2$ results in a totally different approach to tune the optical properties. Usually, an external source of strain is employed to monitor the optical and vibrational properties of a material. It is always challenging to have a precise control over the strain and its consequences on material properties. Here, we report the presence of a compressive strain in Mo$S_2$ crystalline powder and nanosheets obtained via the process of ball-milling and probe sonication. The diffraction peaks in the X-ray diffraction pattern shift to higher 2$\theta$ value implying a compressive strain that increases with the processing time. The absorption spectra, photoluminescence and Raman modes are blue-shifted w.r.t the bulk unprocessed sample. The observed blue-shift is attributed to the presence of compressive strain in the samples. Whereas in thin nano-sheets of Mo$S_2$, it is very likely that both quantum confinement as well as strain result in the observed blue-shift. These results indicate that by optimizing the processing conditions and/or time, a strain of desired amount and hence tunable shift in optical properties of material can be achieved.


title: Perpendicularly Magnetized Ni / Pt (001) Epitaxial Superlattice

abstract: A perpendicularly magnetized ferromagnetic layer is an important building block for recent/future highdensity spintronic memory applications. This paper reports on the fabrication of perpendicularly magnetized Ni / Pt superlattices and the characterization of their structures and magnetic properties. The optimization of film growth conditions allowed us to grow epitaxial Ni / Pt (001) superlattices on SrTiO$_{3}$ (001) single crystal substrates. We investigated their structural parameters and magnetic properties as a function of the Ni layer thickness, and obtained a high uniaxial magnetic anisotropy energy of 1.9 x 10$^{6}$ erg/cm$^{3}$ for a [Ni (4.0 nm) / Pt (1.0 nm)] superlattice. In order to elucidate the detailed mechanism on perpendicular magnetic anisotropy for the Ni / Pt (001) superlattices, the experimental results were compared with the first-principles calculations. It has been found that the strain effect is a prime source of the emergence of perpendicular magnetic anisotropy.


title: Learning Robust and Multilingual Speech Representations

abstract: Unsupervised speech representation learning has shown remarkable success at finding representations that correlate with phonetic structures and improve downstream speech recognition performance. However, most research has been focused on evaluating the representations in terms of their ability to improve the performance of speech recognition systems on read English (e.g. Wall Street Journal and LibriSpeech). This evaluation methodology overlooks two important desiderata that speech representations should have: robustness to domain shifts and transferability to other languages. In this paper we learn representations from up to 8000 hours of diverse and noisy speech data and evaluate the representations by looking at their robustness to domain shifts and their ability to improve recognition performance in many languages. We find that our representations confer significant robustness advantages to the resulting recognition systems: we see significant improvements in out-of-domain transfer relative to baseline feature sets and the features likewise provide improvements in 25 phonetically diverse languages including tonal languages and low-resource languages.


title: Time- and frequency-limited H2-optimal model order reduction of bilinear control systems

abstract: In the time- and frequency-limited model order reduction, a reduced-order approximation of the original high-order model is sought to ensure superior accuracy in some desired time and frequency intervals. We first consider the time-limited H2-optimal model order reduction problem for bilinear control systems and derive first-order optimality conditions that a local optimum reduced-order model should satisfy. We then propose a heuristic algorithm that generates a reduced-order model, which tends to achieve these optimality conditions. The frequency-limited and the time-limited H2-pseudo-optimal model reduction problems are also considered wherein we restrict our focus on constructing a reduced-order model that satisfies a subset of the respective optimality conditions for the local optimum. Two new algorithms have been proposed that enforce two out of four optimality conditions on the reduced-order model upon convergence. The algorithms are tested on three numerical examples to validate the theoretical results presented in the paper. The numerical results confirm the efficacy of the proposed algorithms.


title: Blocked Clusterwise Regression

abstract: A recent literature in econometrics models unobserved cross-sectional heterogeneity in panel data by assigning each cross-sectional unit a one-dimensional, discrete latent type. Such models have been shown to allow estimation and inference by regression clustering methods. This paper is motivated by the finding that the clustered heterogeneity models studied in this literature can be badly misspecified, even when the panel has significant discrete cross-sectional structure. To address this issue, we generalize previous approaches to discrete unobserved heterogeneity by allowing each unit to have multiple, imperfectly-correlated latent variables that describe its response-type to different covariates. We give inference results for a k-means style estimator of our model and develop information criteria to jointly select the number clusters for each latent variable. Monte Carlo simulations confirm our theoretical results and give intuition about the finite-sample performance of estimation and model selection. We also contribute to the theory of clustering with an over-specified number of clusters and derive new convergence rates for this setting. Our results suggest that over-fitting can be severe in k-means style estimators when the number of clusters is over-specified.


title: Developing an Augmented Reality Tourism App through User-Centred Design (Extended Version)

abstract: Augmented Reality (AR) bridges the gap between the physical and virtual world. Through overlaying graphics on natural environments, users can immerse themselves in a tailored environment. This offers great benefits to mobile tourism, where points of interest (POIs) can be annotated on a smartphone screen. While a variety of apps currently exist, usability issues can discourage users from embracing AR. Interfaces can become cluttered with icons, with POI occlusion posing further challenges. In this paper, we use user-centred design (UCD) to develop an AR tourism app. We solicit requirements through a synthesis of domain analysis, tourist observation and semi-structured interviews. Whereas previous user-centred work has designed mock-ups, we iteratively develop a full Android app. This includes overhead maps and route navigation, in addition to a detailed AR browser. The final product is evaluated by 20 users, who participate in a tourism task in a UK city. Users regard the system as usable and intuitive, and suggest the addition of further customisation. We finish by critically analysing the challenges of a user-centred methodology.


title: Describing and Predicting Online Items with Reshare Cascades via Dual Mixture Self-exciting Processes

abstract: It is well-known that online behavior is long-tailed, with most cascaded actions being short and a few being very long. A prominent drawback in generative models for online events is the inability to describe unpopular items well. This work addresses these shortcomings by proposing dual mixture self-exciting processes to jointly learn from groups of cascades. We first start from the observation that maximum likelihood estimates for content virality and influence decay are separable in a Hawkes process. Next, our proposed model, which leverages a Borel mixture model and a kernel mixture model, jointly models the unfolding of a heterogeneous set of cascades. When applied to cascades of the same online items, the model directly characterizes their spread dynamics and supplies interpretable quantities, such as content virality and content influence decay, as well as methods for predicting the final content popularities. On two retweet cascade datasets -- one relating to YouTube videos and the second relating to controversial news articles -- we show that our models capture the differences between online items at the granularity of items, publishers and categories. In particular, we are able to distinguish between far-right, conspiracy, controversial and reputable online news articles based on how they diffuse through social media, achieving an F1 score of 0.945. On holdout datasets, we show that the dual mixture model provides, for reshare diffusion cascades especially unpopular ones, better generalization performance and, for online items, accurate item popularity predictions.


title: Universal Algebraic Controllers and System Identification

abstract: In this document, some structured operator approximation theoretical methods for system identification of nearly eventually periodic systems, are presented. Let $\mathbb{C}^{n\times m}$ denote the algebra of $n\times m$ complex matrices. Given $\varepsilon>0$, an arbitrary discrete-time dynamical system $(\Sigma,\mathcal{T})$ with state-space $\Sigma$ contained in the finite dimensional Hilbert space $\mathbb{C}^n$, whose state-transition map $\mathcal{T}:\Sigma\times ([0,\infty)\cap \mathbb{Z})\to \Sigma$ is unknown or partially known, and needs to be determined based on some sampled data in a finite set $\hat{\Sigma}=\{x_t\}_{1\leq t\leq m}\subset \Sigma$ according to the rule $\mathcal{T}(x_t,1)=x_{t+1}$ for each $1\leq t\leq m-1$, and given $x\in \hat{\Sigma}$. We study the solvability of the existence problems for two triples $(p,A,\varphi)$ and $(p,A_\eta,\Phi)$ determined by a polynomial $p\in \mathbb{C}[z]$ with $\deg(p)\leq m$, a matrix root $A\in\mathbb{C}^{m\times m}$ and an approximate matrix root $A_\eta\in\mathbb{C}^{r\times r}$ of $p(z)=0$ with $r\leq m$, two completely positive linear multiplicative maps $\varphi:\mathbb{C}^{m\times m}\to \mathbb{C}^{n\times n}$ and $\Phi:\mathbb{C}^{r\times r}\to \mathbb{C}^{n\times n}$, such that $\|\mathcal{T}(x,t)-\varphi(A^t)x\|\leq\varepsilon$ and $\|\Phi(A_\eta^t)x-\varphi(A^t)x\|\leq\varepsilon$, for each integer $t\geq 1$ such that $\|\mathcal{T}(x,t)-y\|\leq \varepsilon$ for some $y\in \hat{\Sigma}$. Some numerical implementations of these techniques for the reduced-order predictive simulation of dynamical systems in continuum and quantum mechanics, are outlined.

categories: math.NA cs.NA cs.SY eess.SY math.OA math.OC


title: RTO Frequency Regulation Market Clearing Formulations and the Effect of Opportunity Costs on Electricity Prices

abstract: The Regional Transmission Operators each have unique methods to procure frequency regulation reserves used to track real power fluctuations on the grid. Some market clearing practices result in high regulation prices with large spread compared to other markets. We present the frequency regulation market clearing formulations--as derived from operator tariffs--of the ISO New England, PJM Interconnection, and Midcontinent ISO. We offer test case examples to explain the historic market pricing behavior seen within each system operator and conclude this behavior is due to the inclusion of estimated lost opportunity costs in market clearing.


title: Perturbed rank 2 Poisson systems and periodic orbits on Casimir invariant manifolds

abstract: A class of n-dimensional Poisson systems reducible to an unperturbed harmonic oscillator shall be considered. In such case, perturbations leaving invariant a given symplectic leaf shall be investigated. Our purpose will be to analyze the bifurcation phenomena of periodic orbits as a result of these perturbations in the period annulus associated to the unperturbed harmonic oscillator. This is accomplished via the averaging theory up to an arbitrary order in the perturbation parameter e. In that theory we shall also use both branching theory and singularity theory of smooth maps to analyze the bifurcation phenomena at points where the implicit function theorem is not applicable. When the perturbation is given by a polynomial family, the associated Melnikov functions are polynomial and tools of computational algebra based on Gr\"obner basis are employed in order to reduce the generators of some polynomial ideals needed to analyze the bifurcation problem. When the most general perturbation of the harmonic oscillator by a quadratic perturbation field is considered, the complete bifurcation diagram (except at a high codimension subset) in the parameter space is obtained. Examples are given.

categories: math-ph math.CA math.MP math.SG nlin.SI physics.class-ph


title: The Secret is in the Spectra: Predicting Cross-lingual Task Performance with Spectral Similarity Measures

abstract: Performance in cross-lingual NLP tasks is impacted by the (dis)similarity of languages at hand: e.g., previous work has suggested there is a connection between the expected success of bilingual lexicon induction (BLI) and the assumption of (approximate) isomorphism between monolingual embedding spaces. In this work we present a large-scale study focused on the correlations between monolingual embedding space similarity and task performance, covering thousands of language pairs and four different tasks: BLI, parsing, POS tagging and MT. We hypothesize that statistics of the spectrum of each monolingual embedding space indicate how well they can be aligned. We then introduce several isomorphism measures between two embedding spaces, based on the relevant statistics of their individual spectra. We empirically show that 1) language similarity scores derived from such spectral isomorphism measures are strongly associated with performance observed in different cross-lingual tasks, and 2) our spectral-based measures consistently outperform previous standard isomorphism measures, while being computationally more tractable and easier to interpret. Finally, our measures capture complementary information to typologically driven language distance measures, and the combination of measures from the two families yields even higher task performance correlations.


title: Adversarial Attacks on Convolutional Neural Networks in Facial Recognition Domain

abstract: Numerous recent studies have demonstrated how Deep Neural Network (DNN) classifiers can be fooled by adversarial examples, in which an attacker adds perturbations to an original sample, causing the classifier to misclassify the sample. Adversarial attacks that render DNNs vulnerable in real life represent a serious threat in autonomous vehicles, malware filters, or biometric authentication systems. In this paper, we apply Fast Gradient Sign Method to introduce perturbations to a facial image dataset and then test the output on a different classifier that we trained ourselves, to analyze transferability of this method. Next, we craft a variety of different black-box attack algorithms on a facial image dataset assuming minimal adversarial knowledge, to further assess the robustness of DNNs in facial recognition. While experimenting with different image distortion techniques, we focus on modifying single optimal pixels by a large amount, or modifying all pixels by a smaller amount, or combining these two attack approaches. While our single-pixel attacks achieved about a 15% average decrease in classifier confidence level for the actual class, the all-pixel attacks were more successful and achieved up to an 84% average decrease in confidence, along with an 81.6% misclassification rate, in the case of the attack that we tested with the highest levels of perturbation. Even with these high levels of perturbation, the face images remained identifiable to a human. Understanding how these noised and perturbed images baffle the classification algorithms can yield valuable advances in the training of DNNs against defense-aware adversarial attacks, as well as adaptive noise reduction techniques. We hope our research may help to advance the study of adversarial attacks on DNNs and defensive mechanisms to counteract them, particularly in the facial recognition domain.


title: Narrowband Ca Photometry for Dwarf Spheroidal Galaxies. I. Chemostructural Study on Draco, Sextans, and Canes Venatici I

abstract: A few dozen dwarf satellite galaxies of the Milky Way have been discovered, which are often viewed as the remaining building blocks of our Galaxy. The follow-up spectroscopy showed that dwarf galaxies have a sizeable spread in their metallicities. Several scenarios were suggested to explain the metallicity spread, which can be tested by the structural patterns of stellar subpopulations with distinct metallicities. However, such chemical plus structural examination, to which we refer to as "chemostructural study" is hindered by the lack of stars with spectroscopic metallicity. Here we propose the Ca$-$$by$ photometry as an alternative way to secure metallicities for a 2$-$3 orders of magnitude larger stellar sample than the spectroscopic sample and thus enable us to perform a chemostructural study on dwarf galaxies. In particular, we use the $hk$ index [$\equiv($Ca$-b)-(b-y)$], whose validity as a photometric metallicity indicator (and crass insensitivity to age) for red-giant-branch stars was upheld via Galactic globular clusters, and observe three dwarf spheroidal galaxies$-$Draco, Sextans, and Canes Venatici I (CVnI)$-$with Subaru/Suprime-Cam. We find that in all the galaxies the metal-rich stellar populations are more centrally concentrated than the metal-poor counterparts, suggesting that the central regions of the galaxies underwent extended star formation. Such a negative radial metallicity gradient for Sextans and CVnI opposes to the traditional spectroscopic results. We also find that their metallicity distribution functions (MDFs) can be characterized by a unimodal, skewed Gaussian shape with a metal-rich peak and a metal-poor tail. We discuss their features in the chemostructure and MDF in terms of dwarf galaxy formation theories.


title: Error estimation and uncertainty quantification for first time to a threshold value

abstract: Classical a posteriori error analysis for differential equations quantifies the error in a Quantity of Interest (QoI) which is represented as a bounded linear functional of the solution. In this work we consider a posteriori error estimates of a quantity of interest that cannot be represented in this fashion, namely the time at which a threshold is crossed for the first time. We derive two representations for such errors and use an adjoint-based a posteriori approach to estimate unknown terms that appear in our representation. The first representation is based on linearizations using Taylor's Theorem. The second representation is obtained by implementing standard root-finding techniques. We provide several examples which demonstrate the accuracy of the methods. We then embed these error estimates within a framework to provide error bounds on a cumulative distribution function when parameters of the differential equations are uncertain.


title: The Baghdad Atlas: A relational database of inelastic neutron-scattering $(n,n'\gamma)$ data

abstract: A relational database has been developed based on the original ($n,n'\gamma$) work carried out by A. M. Demidov $et$ $al$., at the Nuclear Research Institute in Baghdad, Iraq [$"Atlas$ $of$ $Gamma$-$Ray$ $Spectra$ $from$ $the$ $Inelastic$ $Scattering$ $of$ $Reactor$ $Fast$ $Neutrons"$, Nuclear Research Institute, Baghdad, Iraq (Moscow, Atomizdat 1978)] for 105 independent measurements comprising 76 elemental samples of natural composition and 29 isotopically-enriched samples. The information from this Atlas includes: $\gamma$-ray energies and relative intensities; nuclide and level data corresponding to the residual nucleus and meta data associated with the target sample that allows for the extraction of the flux-weighted ($n,n'\gamma$) cross sections for a given transition relative to a defined value. The optimized angular-distribution-corrected fast-neutron flux-weighted partial $\gamma$-ray cross section for the production of the 846.8-keV $2^{+}_{1} \rightarrow 0^{+}_{\rm gs}$ $\gamma$-ray transition in $^{56}$Fe, determined to be $\langle \sigma_{\gamma} \rangle = 143(29)$ mb, is used for this purpose. However, different values for the adopted cross section can be readily implemented to accommodate user preference based on revised determinations of this quantity. The Atlas ($n,n'\gamma$) data has been compiled into a series of CSV-style ASCII data sets and a suite of Python scripts have been developed to build and install the database locally. The database can then be accessed directly through the SQLite engine, or using alternative methods such as the Jupyter Notebook Python-browser interface. Several examples exploiting different interaction methodologies are distributed with the complete software package.


title: The Algorithm of Factor 12: Generating the Information Carved on the Plimpton 322 Tablet

abstract: The right triangles represented in the Plimpton 322 data have integer sides $(a, b, d)$ with $a < b < d$, by inspecting the data we propose $b=M\,Q_{M}$ with $M,Q_{M}$ positive integers. In this way we present a new data-driven family of algorithms in order to generate the whole content of the Plimpton-322 tablet. Each algorithm of the family corresponds to one given value of the bundling factor or mak\c{s}arum $M$. The algorithm with $M=12$ arise naturally from data and is developed entirely through this work. We could generate the complete Plimpton 322 table by using three different schemes in the framework of the twelve factor algorithm. Finally we found by varying $M$ in the general algorithm, we can be able to find all the Pythagorean Triples.


title: VERONICA: Expressive and Precise Concurrent Information Flow Security (Extended Version with Technical Appendices)

abstract: Methods for proving that concurrent software does not leak its secrets has remained an active topic of research for at least the past four decades. Despite an impressive array of work, the present situation remains highly unsatisfactory. With contemporary compositional proof methods one is forced to choose between expressiveness (the ability to reason about a wide variety of security policies), on the one hand, and precision (the ability to reason about complex thread interactions and program behaviours), on the other. Achieving both is essential and, we argue, requires a new style of compositional reasoning. We present VERONICA, the first program logic for proving concurrent programs information flow secure that supports compositional, high-precision reasoning about a wide range of security policies and program behaviours (e.g. expressive de-classification, value-dependent classification, secret-dependent branching). Just as importantly, VERONICA embodies a new approach for engineering such logics that can be re-used elsewhere, called decoupled functional correctness (DFC). DFC leads to a simple and clean logic, even while achieving this unprecedented combination of features. We demonstrate the virtues and versatility of VERONICA by verifying a range of example programs, beyond the reach of prior methods.


title: A Graph-Based Approach for Active Learning in Regression

abstract: Active learning aims to reduce labeling efforts by selectively asking humans to annotate the most important data points from an unlabeled pool and is an example of human-machine interaction. Though active learning has been extensively researched for classification and ranking problems, it is relatively understudied for regression problems. Most existing active learning for regression methods use the regression function learned at each active learning iteration to select the next informative point to query. This introduces several challenges such as handling noisy labels, parameter uncertainty and overcoming initially biased training data. Instead, we propose a feature-focused approach that formulates both sequential and batch-mode active regression as a novel bipartite graph optimization problem. We conduct experiments on both noise-free and noisy settings. Our experimental results on benchmark data sets demonstrate the effectiveness of our proposed approach.


title: A Review of Cybersecurity Incidents in the Water Sector

abstract: This study presents a critical review of disclosed, documented, and malicious cybersecurity incidents in the water sector to inform safeguarding efforts against cybersecurity threats. The review is presented within a technical context of industrial control system architectures, attack-defense models, and security solutions. Fifteen incidents were selected and analyzed through a search strategy that included a variety of public information sources ranging from federal investigation reports to scientific papers. For each individual incident, the situation, response, remediation, and lessons learned were compiled and described. The findings of this review indicate an increase in the frequency, diversity, and complexity of cyberthreats to the water sector. Although the emergence of new threats, such as ransomware or cryptojacking, was found, a recurrence of similar vulnerabilities and threats, such as insider threats, was also evident, emphasizing the need for an adaptive, cooperative, and comprehensive approach to water cyberdefense.


title: Computing real time correlation functions on a hybrid classical/quantum computer

abstract: Quantum devices may overcome limitations of classical computers in studies of nuclear structure functions and parton Wigner distributions of protons and nuclei. In this talk, we discuss a worldline approach to compute nuclear structure functions in the high energy Regge limit of QCD using a hybrid quantum computer, by expressing the fermion determinant in the QCD path integral as a quantum mechanical path integral over $0+1$-dimensional fermionic and bosonic world-lines in background gauge fields. Our simplest example of computing the well-known dipole model result for the structure function $F_2$ in the high energy Regge limit is feasible with NISQ era technology using few qubits and shallow circuits. This example can be scaled up in complexity and extended in scope to compute structure functions, scattering amplitudes and other real-time correlation functions in QCD, relevant for example to describe non-equilibrium transport of quarks and gluons in a Quark-Gluon-Plasma.


title: ART$^2$: A 3D Parallel Multi-wavelength Radiative Transfer Code for Continuum and Atomic and Molecular Lines

abstract: ART$^2$ is a 3D multi-wavelength Monte Carlo radiative transfer (RT) code that couples continuum and emission lines to track the propagation of photons and their interactions with the interstellar medium (ISM). The previous version of ART$^2$, which included continuum and Ly$\alpha$ line, has been extensively applied to hydrodynamics simulations to study multi-band properties of galaxies and ISM. Here, we describe new implementations of non-local thermodynamic equilibrium RT of molecular and atomic fine structure emission lines, and the parallelization of the code using a number of novel methods. The new ART$^2$ can efficiently and self-consistently produce a full spectrum that includes both continuum and lines such as [CII], [NII], [OIII], Ly$\alpha$, and CO. These essential features, together with the multi-phase ISM model and the adaptive grid, make ART$^2$ a multi-purpose code to study multi-wavelength properties of a wide range of astrophysical systems from planetary disks to large-scale structures. To demonstrate the capability of the new ART$^2$, we applied it to two hydrodynamics simulations: the zoom-in Milky Way Simulation to obtain panchromatic properties of individual galaxies, and the large-scale IllustrisTNG100 Simulation to obtain global properties such as the line intensity mappings. These products are essential for a broad array of studies such as the correlations between physical and panchromatic properties and their evolution. By enabling direct comparison between numerical simulations and multi-band observations, ART$^2$ provides a crucial theoretical framework for the interpretations of existing observations, the plan for future surveys, and the synergy between multi-band galaxy surveys and line intensity mappings. Therefore, ART$^2$ is a powerful and versatile tool to bridge the gap between theories and observations of cosmic structures.


title: Mitigation of Backgrounds from Cosmogenic $^{137}$Xe in Xenon Gas Experiments using $^{3}$He Neutron Capture

abstract: \Xe{136} is used as the target medium for many experiments searching for \bbnonu. Despite underground operation, cosmic muons that reach the laboratory can produce spallation neutrons causing activation of detector materials. A potential background that is difficult to veto using muon tagging comes in the form of \Xe{137} created by the capture of neutrons on \Xe{136}. This isotope decays via beta decay with a half-life of 3.8 minutes and a \Qb\ of $\sim$4.16 MeV. This work proposes and explores the concept of adding a small percentage of \He{3} to xenon as a means to capture thermal neutrons and reduce the number of activations in the detector volume. When using this technique we find the contamination from \Xe{137} activation can be reduced to negligible levels in tonne and multi-tonne scale high pressure gas xenon neutrinoless double beta decay experiments running at any depth in an underground laboratory.


title: Neutrino transport with Monte Carlo method: I. Towards fully consistent implementation of nucleon recoils in core-collapse supernova simulations

abstract: The small energy exchange via nucleon recoils in neutrino-nucleon scattering is now supposed to be one of the important factors for successful explosion of core-collapse supernovae (CCSNe) as they can change neutrino spectra through accumulation of a large number of scatterings. In finite-difference methods employed for neutrino transport in CCSN simulations, we normally can not afford to deploy a large enough number of energy bins needed to resolve this small energy exchange and sub-grid techniques are employed one way or another. In this paper we study quantitatively with the Monte Carlo (MC) method how well such a treatment performs. We first investigate the effects of nucleon recoils on the neutrino spectra and confirm that the average energy is reduced by $\sim$15% for heavy-lepton neutrinos and by much smaller quantities for other types of neutrinos in a typical post-bounce situation. It is also observed that the nucleon scattering dominates the electron scattering in the thermalization of neutrino spectra in all flavors. We then study possible artifacts that the coarse energy grid may produce in the finite-difference methods. In order to mimic the latter calculation, we re-distribute MC particles in each energy bin after a certain interval in a couple of ways and study how the results are affected and depend on the energy-resolution. We also discuss possible implications of our results for the finite-difference methods.


title: Updated constraints on superconducting cosmic strings from the astronomy of fast radio bursts

abstract: In this article we update constraints on superconducting cosmic strings (SCSs) in light of the recent observational developments of fast radio bursts (FRBs) astronomy. Assuming strings follow an exponential distribution characterized by current, we show that two parameters in our context, which are the characteristic tension ($G\mu$) and a parameter which describes the aforementioned exponential distribution ($I_c$), can be constrained by FRB experiments. Particularly, we investigate data sets from Parkes and ASKAP. We looked through a parameter space where $G\mu \sim [10^{-17}, 10^{-12}]$ and $I_c\sim [10^{-1}, 10^{2}]$ GeV, and found our results show that Parkes jointly with ASKAP can constrain the parameter space for SCSs.


title: Analysis of Y00 Protocol under Quantum Generalization of a Fast Correlation Attack: Toward Information-Theoretic Security

abstract: In our previous work, it was demonstrated that the attacker could not pin-down the correct keys to start the Y00 protocol with a probability of one under the assistance of unlimitedly long known-plaintext attacks and optimal quantum measurements on the attacker's quantum memory. However, there were several assumptions that the Y00 system utilized linear-feedback shift registers as pseudo-random-number generators, and a fast correlation attack was disabled by irregular mapping on the Y00 systems. This study generalizes such an attack to remove the assumptions of the previous work. The framework of the security analyses of this study reiterates two well-known results from the past: (1) Y00 systems would be cryptanalyzed when the system is not designed well; (2) the system is possibly information-theoretically secure when the system is designed well, although the attacker's confidence in the correct key increases over time but the success probability of key recovery does not reach unity in finite time; (3) the breach time of the shared keys is increasingly threatened with time. Hence, a key-refreshment procedure for the Y00 protocol is provided herein. Such security analyses are important not only in key refreshment but also in initial key agreement situations.


title: Steady-state diffusion approximations of Markov chains: error analysis via the discrete Poisson equation

abstract: This paper uses the generator approach of Stein's method to analyze the gap between steady-state distributions of Markov chains and diffusion processes. Until now, the standard way to invoke Stein's method for this problem was to use the Poisson equation for the diffusion as a starting point. The main technical difficulty with this approach is obtaining bounds on the derivatives of the solution to the Poisson equation, also known as Stein factor bounds. In this paper we propose starting with the discrete Poisson equation of the Markov chain. An important step in our approach is extending the discrete Poisson equation to be defined on the continuum where the diffusion is defined, and we achieve this by using interpolation. Although there are still Stein factor bounds to prove, these now correspond to the finite differences of the discrete Poisson equation solution, as opposed to the derivatives of the solution to the continuous one. Discrete Stein factor bounds can be easier to obtain, for instance when the drift is not everywhere differentiable, when the diffusion has a state-dependent diffusion coefficient, or in the presence of a reflecting boundary condition. We use the join the shortest queue model in the Halfin-Whitt regime as a working example to illustrate the methodology. We show that the steady-state approximation error of the diffusion limit converges to zero at a rate of $1/\sqrt{n}$, where $n$ is the number of servers in the system.


title: Adversarial Incremental Learning

abstract: Although deep learning performs really well in a wide variety of tasks, it still suffers from catastrophic forgetting -- the tendency of neural networks to forget previously learned information upon learning new tasks where previous data is not available. Earlier methods of incremental learning tackle this problem by either using a part of the old dataset, by generating exemplars or by using memory networks. Although, these methods have shown good results but using exemplars or generating them, increases memory and computation requirements. To solve these problems we propose an adversarial discriminator based method that does not make use of old data at all while training on new tasks. We particularly tackle the class incremental learning problem in image classification, where data is provided in a class-based sequential manner. For this problem, the network is trained using an adversarial loss along with the traditional cross-entropy loss. The cross-entropy loss helps the network progressively learn new classes while the adversarial loss helps in preserving information about the existing classes. Using this approach, we are able to outperform other state-of-the-art methods on CIFAR-100, SVHN, and MNIST datasets.


title: Characterizing Transition-Metal Dichalcogenide Thin-Films using Hyperspectral Imaging and Machine Learning

abstract: Atomically thin polycrystalline transition-metal dichalcogenides (TMDs) are relevant to both fundamental science investigation and applications. TMD thin-films present uniquely difficult challenges to effective nanoscale crystalline characterization. Here we present a method to quickly characterize the nanocrystalline grain structure and texture of monolayer WS2 films using scanning nanobeam electron diffraction coupled with multivariate statistical analysis of the resulting data. Our analysis pipeline is highly generalizable and is a useful alternative to the time consuming, complex, and system-dependent methodology traditionally used to analyze spatially resolved electron diffraction measurements.


title: Multi-Participant Multi-Class Vertical Federated Learning

abstract: Federated learning (FL) is a privacy-preserving paradigm for training collective machine learning models with locally stored data from multiple participants. Vertical federated learning (VFL) deals with the case where participants sharing the same sample ID space but having different feature spaces, while label information is owned by one participant. Current studies of VFL only support two participants, and mostly focus on binaryclass logistic regression problems. In this paper, we propose the Multi-participant Multi-class Vertical Federated Learning (MMVFL) framework for multi-class VFL problems involving multiple parties. Extending the idea of multi-view learning (MVL), MMVFL enables label sharing from its owner to other VFL participants in a privacypreserving manner. To demonstrate the effectiveness of MMVFL, a feature selection scheme is incorporated into MMVFL to compare its performance against supervised feature selection and MVL-based approaches. Experiment results on real-world datasets show that MMVFL can effectively share label information among multiple VFL participants and match multi-class classification performance of existing approaches.


title: An automatic deep learning-based workflow for glioblastoma survival prediction using pre-operative multimodal MR images

abstract: We proposed a fully automatic workflow for glioblastoma (GBM) survival prediction using deep learning (DL) methods. 285 glioma (210 GBM, 75 low-grade glioma) patients were included. 163 of the GBM patients had overall survival (OS) data. Every patient had four pre-operative MR scans and manually drawn tumor contours. For automatic tumor segmentation, a 3D convolutional neural network (CNN) was trained and validated using 122 glioma patients. The trained model was applied to the remaining 163 GBM patients to generate tumor contours. The handcrafted and DL-based radiomic features were extracted from auto-contours using explicitly designed algorithms and a pre-trained CNN respectively. 163 GBM patients were randomly split into training (n=122) and testing (n=41) sets for survival analysis. Cox regression models with regularization techniques were trained to construct the handcrafted and DL-based signatures. The prognostic power of the two signatures was evaluated and compared. The 3D CNN achieved an average Dice coefficient of 0.85 across 163 GBM patients for tumor segmentation. The handcrafted signature achieved a C-index of 0.64 (95% CI: 0.55-0.73), while the DL-based signature achieved a C-index of 0.67 (95% CI: 0.57-0.77). Unlike the handcrafted signature, the DL-based signature successfully stratified testing patients into two prognostically distinct groups (p-value<0.01, HR=2.80, 95% CI: 1.26-6.24). The proposed 3D CNN generated accurate GBM tumor contours from four MR images. The DL-based signature resulted in better GBM survival prediction, in terms of higher C-index and significant patient stratification, than the handcrafted signature. The proposed automatic radiomic workflow demonstrated the potential of improving patient stratification and survival prediction in GBM patients.


title: Mass inflation and the $C^2$-inextendibility of spherically symmetric charged scalar field dynamical black holes

abstract: It has long been suggested that the Cauchy horizon of dynamical black holes is subject to a weak null singularity, under the mass inflation scenario. We study in spherical symmetry the Einstein-Maxwell-Klein-Gordon equations and \textit{while we do not directly show mass inflation}, we obtain a "mass inflation/ridigity" dichotomy. More precisely, we prove assuming (sufficiently slow) decay of the charged scalar field on the event horizon, that the Cauchy horizon emanating from time-like infinity is $\mathcal{CH}_{i^+}= \mathcal{D} \cup \mathcal{S}$ for two (possibly empty) disjoint connected sets $\mathcal{D}$ and $\mathcal{S}$ such that: _$\mathcal{D}$ (the dynamical set) is a past set on which the Hawking mass blows up (mass inflation scenario). _$\mathcal{S}$ (the static set) is a future set isometric to a Reissner--Nordstr\"{o}m Cauchy horizon i.e.\ the radiation is zero on $\mathcal{S}$. As a consequence, we establish a novel classification of Cauchy horizons into three types: dynamical ($\mathcal{S}=\emptyset$), static ($\mathcal{D}=\emptyset$) or mixed, and prove that $\mathcal{CH}_{i^+}$ is globally $C^2$-inextendible. Our main motivation is the $C^2$ Strong Cosmic Censorship Conjecture for a realistic model of spherical collapse in which charged matter emulates the repulsive role of angular momentum: in our case the Einstein-Maxwell-Klein-Gordon system on one-ended space-times. As a result, we prove in spherical symmetry that: - two-ended asymptotically flat space-times are $C^2$-future-inextendible i.e. $C^2$ Strong Cosmic Censorship is true for Einstein-Maxwell-Klein-Gordon, assuming the decay of the scalar field on the event horizon at the expected rate. - In the one-ended case, the Cauchy horizon emanating from time-like infinity is $C^2$-inextendible. This result suppresses the main obstruction to $C^2$ Strong Cosmic Censorship in spherical collapse.


title: On Duality in $\mathcal{N}=2$ supersymmetric Liouville Theory

abstract: Similarly to the bosonic Liouville theory, the $\mathcal{N}=2$ supersymmetric Liouville theory was conjectured to be equipped with the duality that exchanges the superpotential and the K\"ahler potential. The conjectured duality, however, seems to suffer from a mismatch of the preserved symmetries. More than fifteen years ago, when I was a student, my supervisor Tohru Eguchi gave a beautiful resolution of the puzzle when the supersymmetry is enhanced to $\mathcal{N}=4$ based on his insight into the underlying geometric structure of the $A_1$ singularity. I will review his unpublished but insightful idea and present our attempts to extend it to more general cases.


title: AVATAR -- Machine Learning Pipeline Evaluation Using Surrogate Model

abstract: The evaluation of machine learning (ML) pipelines is essential during automatic ML pipeline composition and optimisation. The previous methods such as Bayesian-based and genetic-based optimisation, which are implemented in Auto-Weka, Auto-sklearn and TPOT, evaluate pipelines by executing them. Therefore, the pipeline composition and optimisation of these methods requires a tremendous amount of time that prevents them from exploring complex pipelines to find better predictive models. To further explore this research challenge, we have conducted experiments showing that many of the generated pipelines are invalid, and it is unnecessary to execute them to find out whether they are good pipelines. To address this issue, we propose a novel method to evaluate the validity of ML pipelines using a surrogate model (AVATAR). The AVATAR enables to accelerate automatic ML pipeline composition and optimisation by quickly ignoring invalid pipelines. Our experiments show that the AVATAR is more efficient in evaluating complex pipelines in comparison with the traditional evaluation approaches requiring their execution.


title: Universally Safe Swerve Manoeuvres for Autonomous Driving

abstract: This paper characterizes safe following distances for on-road driving when vehicles can avoid collisions by either braking or by swerving into an adjacent lane. In particular, we focus on safety as defined in the Responsibility-Sensitive Safety (RSS) framework. We extend RSS by introducing swerve manoeuvres as a valid response in addition to the already present brake manoeuvre. These swerve manoeuvres use the more realistic kinematic bicycle model rather than the double integrator model of RSS. When vehicles are able to swerve and brake, it is shown that their required safe following distance at higher speeds is less than that required through braking alone. In addition, when all vehicles follow this new distance, they are provably safe. The use of the kinematic bicycle model is then validated by comparing these swerve manoeuvres to that of a dynamic single-track model.


title: Ability of Markovian Master Equations to Model Quantum Computers and Other Systems Under Broadband Control

abstract: Most future quantum devices, including quantum computers, require control that is broadband, meaning that the rate of change of the time-dependent Hamiltonian is as fast or faster than the dynamics it generates. In many areas of quantum physics, including quantum technology, one must include dissipation and decoherence induced by the environment. While Markovian master equations provide the only really efficient way to model these effects, these master equations are derived for constant Hamiltonians (or those with a discrete set of well-defined frequencies). In 2006, Alicky, Lidar, and Zanardi [Phys. Rev. A 73, 052311 (2006)] provided detailed qualitative arguments that Markovian master equations could not describe systems under broadband control. Despite apparently broad acceptance of these arguments, such master equations are routinely used to model precisely these systems. This odd state of affairs is likely due to a lack of quantitative results. Here we perform exact simulations of two- and three-level systems coupled to an oscillator bath to obtain quantitative results. Although we confirm that in general Markovian master equations cannot predict the effects of damping under broadband control, we find that there is a widely applicable regime in which they can. Master equations are accurate for weak damping if both the Rabi frequencies and bandwidth of the control are significantly smaller than the system's transition frequencies. They also remain accurate if the bandwidth of control is as large as the frequency of the driven transition so long as this bandwidth does not overlap other transitions. Master equations are thus able to provide accurate descriptions of many quantum information processing protocols for atomic systems.


title: On the Chemical and Kinematic Consistency Between N-rich Metal-poor Field Stars and Enriched Populations in Globular Clusters

abstract: Interesting chemically peculiar field stars may reflect their stellar evolution history and their possible origin in a different environment from where they are found now, which is one of the most important research fields in Galactic archaeology. To explore this further, we have used the CN-CH bands around 4000 A to identify N-rich metal-poor field stars in LAMOST DR3. Here we expand our N-rich metal-poor field star sample to ~100 stars in LAMOST DR5, where 53 of them are newly found in this work. We investigate light elements of the common stars between our sample and APOGEE DR14. While Mg, Al, and Si abundances generally agree with the hypothesis that N-rich metal-poor field stars come from enriched populations in globular clusters, it is still inconclusive for C, N, and O. After integrating the orbits of our N-rich field stars and a control sample of normal metal-poor field stars, we find that N-rich field stars have different orbital parameter distributions compared to the control sample, specifically, apocentric distances, maximum vertical amplitude (Zmax), orbital energy, and z direction angular momentum (Lz). The orbital parameters of N-rich field stars indicate that most of them are inner-halo stars. The kinematics of N-rich field stars support their possible GC origin. The spatial and velocity distributions of our bona fide N-rich field star sample are important observational evidence to constrain simulations of the origin of these interesting objects.


title: On the Age of Information in Internet of Things Systems with Correlated Devices

abstract: In this paper, a real-time Internet of Things (IoT) monitoring system is considered in which multiple IoT devices must transmit timely updates on the status information of a common underlying physical process to a common destination. In particular, a real-world IoT scenario is considered in which multiple (partially) observed status information by different IoT devices are required at the destination, so that the real-time status of the physical process can be properly re-constructed. By taking into account such correlated status information at the IoT devices, the problem of IoT device scheduling is studied in order to jointly minimize the average age of information (AoI) at the destination and the average energy cost at the IoT devices. Particularly, two types of IoT devices are considered: Type-I devices whose status updates randomly arrive and type-II devices whose status updates can be generated-at-will with an associated sampling cost. This stochastic problem is formulated as an infinite horizon average cost Markov decision process (MDP). The optimal scheduling policy is shown to be threshold-based with respect to the AoI at the destination, and the threshold is non-increasing with the channel condition of each device. For a special case in which all devices are type-II, the original MDP can be reduced to an MDP with much smaller state and action spaces. The optimal policy is further shown to have a similar threshold-based structure and the threshold is non-decreasing with an energy cost function of the devices. Simulation results illustrate the structure of the optimal policy and show the effectiveness of the optimal policy compared with a myopic baseline policy.


title: Visual Exploration of Movement Relatedness for Multi-species Ecology Analysis

abstract: Advances in GPS telemetry technology have enabled analysis of animal movement in open areas. Ecologists today are utilizing modern analytic tools to study animal behaviors from large quantity of GPS coordinates. Analytic tools with automatic event extraction functionality can be used to investigate potential interactions between animals by locating relevant segments in movement trajectories. However, such automation can easily overlook the spatial, temporal, social context as well as potential data problems. To this end, this paper explores the visual presentations that also clarify the spatial-temporal contexts, social surroudings, as well as underlying data uncertainties of multi-species animal interactions. The outcome system presents the proximity-based, time-varying relatedness between animal entities through pairwise (PW) or individual-to-group (i-G) perspectives. Focusing on the relational aspects, we employ both static depictions and animations to communicate the travelling of individuals. Our contributions are a novel visualization system that helps investigate the subtle variations of long term spatial-temporal relatedness while considering small group patterns. Our evaluation with movement ecologists shows that the system gives them quick access to valuable clues in discovering insights into multi-species movements and signs of potential interactions.