title: Parameterized Hypercomplex Graph Neural Networks for Graph Classification

abstract: Despite recent advances in representation learning in hypercomplex (HC) space, this subject is still vastly unexplored in the context of graphs. Motivated by the complex and quaternion algebras, which have been found in several contexts to enable effective representation learning that inherently incorporates a weight-sharing mechanism, we develop graph neural networks that leverage the properties of hypercomplex feature transformation. In particular, in our proposed class of models, the multiplication rule specifying the algebra itself is inferred from the data during training. Given a fixed model architecture, we present empirical evidence that our proposed model incorporates a regularization effect, alleviating the risk of overfitting. We also show that for fixed model capacity, our proposed method outperforms its corresponding real-formulated GNN, providing additional confirmation for the enhanced expressivity of HC embeddings. Finally, we test our proposed hypercomplex GNN on several open graph benchmark datasets and show that our models reach state-of-the-art performance while consuming a much lower memory footprint with 70& fewer parameters. Our implementations are available at https://github.com/bayer-science-for-a-better-life/phc-gnn.


title: Brought to Light II: Revealing the Origins of Cloaked Spiral Features in Cluster Passive Dwarf Galaxies

abstract: In our companion paper (Brought to Light I: Michea et al. 2021), we reveal spectacular spiral galaxy-like features in deep optical imaging of nine Virgo early-type dwarf galaxies, hidden beneath a dominating smooth stellar disk. Using a new combination of approaches, we find that bar- and spiral-like features contribute 2.2-6.4% of the total flux within 2 R$_{\rm{eff}}$. In this study, we conduct high resolution simulations of cluster harassment of passive dwarf galaxies. Following close pericenter passages of the cluster core, tidal triggering generates features in our model disks that bear a striking resemblance to the observed features. However, we find the disks must be highly rotationally supported (V$_{\rm{peak}}/\sigma_0 \sim 3$), much higher than typically observed. We propose that some early-type dwarfs may contain a few percent of their mass in a cold, thin disk which is buried in the light of a hot, diffuse disk, and only revealed when they undergo tidal triggering. The red optical colors of our sample do not indicate any recent significant star formation, and our simulations show that very plunging pericenter passages (r$_{\rm{peri}}<0.25$r$_{\rm{vir}}$) are required for tidal triggering. Thus, many cluster early-type dwarfs with less plunging orbits may host a yet undetected cold stellar disk component. We discuss possible origin scenarios and consider why similar mass star-forming galaxies in the field are significantly more thin disk dominated than in our cluster sample.


title: Competing correlated states around the zero field Wigner crystallization transition of electrons in two-dimensions

abstract: The competition between kinetic energy and Coulomb interactions in electronic systems can lead to complex many-body ground states with competing superconducting, charge density wave, and magnetic orders. Here we study the low temperature phases of a strongly interacting zinc-oxide-based high mobility two dimensional electron system that displays a tunable metal-insulator transition. Through a comprehensive analysis of the dependence of electronic transport on temperature, carrier density, in-plane and perpendicular magnetic fields, and voltage bias, we provide evidence for the existence of competing correlated metallic and insulating states with varying degrees of spin polarization. Our system features an unprecedented level of agreement with the state-of-the-art Quantum Monte Carlo phase diagram of the ideal jellium model, including a Wigner crystallization transition at a value of the interaction parameter $r_s\sim 30$ and the absence of a pure Stoner transition. In-plane field dependence of transport reveals a new low temperature state with partial spin polarization separating the spin unpolarized metal and the Wigner crystal, which we examine against possible theoretical scenarios such as an anti-ferromagnetic crystal, Coulomb induced micro-emulsions, and disorder driven puddle formation.


title: Does Double Biofeedback Affect Functional Hemispheric Asymmetry and Activity? A Pilot Study

abstract: In the current pilot study, we attempt to find out how double neurofeedback influences functional hemispheric asymmetry and activity. We examined 30 healthy participants (8 males; 22 females, mean age = 29; SD = 8). To measure functional hemispheric asymmetry and activity, we used computer laterometry in the "two-source" lead-lag dichotic paradigm. Double biofeedback included 8 min of EEG oscillation recording with five minutes of basic mode. During the basic mode, the current amplitude of the EEG oscillator gets transformed into feedback sounds while the current amplitude of alpha EEG oscillator is used to modulate the intensity of light signals. Double neurofeedback did not directly influence the asymmetry itself but accelerated individual sound perception characteristics during dichotic listening in the preceding effect paradigm. Further research is needed to investigate the effect of double neurofeedback training on functional brain activity and asymmetry, taking into account participants' age, gender, and motivation.


title: Implementation of Monte-Carlo transport in the general relativistic SpEC code

abstract: Neutrino transport and neutrino-matter interactions are known to play an important role in the evolution of neutron star mergers, and of their post-merger remnants. Neutrinos cool remnants, drive post-merger winds, and deposit energy in the low-density polar regions where relativistic jets may eventually form. Neutrinos also modify the composition of the ejected material, impacting the outcome of nucleosynthesis in merger outflows and the properties of the optical/infrared transients that they power (kilonovae). So far, merger simulations have largely relied on approximate treatments of the neutrinos (leakage, moments) that simplify the equations of radiation transport in a way that makes simulations more affordable, but also introduces unquantifiable errors in the results. To improve on these methods, we recently published a first simulation of neutron star mergers using a low-cost Monte-Carlo algorithm for neutrino radiation transport. Our transport code limits costs in optically thick regions by placing a hard ceiling on the value of the absorption opacity of the fluid, yet all approximations made within the code are designed to vanish in the limit of infinite numerical resolution. We provide here an in-depth description of this algorithm, of its implementation in the SpEC merger code, and of the expected impact of our approximations in optically thick regions. We argue that the latter is a subdominant source of error at the accuracy reached by current simulations, and for the interactions currently included in our code. We also provide tests of the most important features of this code.


title: Composably secure data processing for Gaussian-modulated continuous variable quantum key distribution

abstract: Continuous-variable (CV) quantum key distribution (QKD) employs the quadratures of a bosonic mode to establish a secret key between two remote parties, and this is usually achieved via a Gaussian modulation of coherent states. The resulting secret key rate depends not only on the loss and noise in the communication channel, but also on a series of data processing steps that are needed for transforming shared correlations into a final string of secret bits. Here we consider a Gaussian-modulated coherent-state protocol with homodyne detection in the general setting of composable finite-size security. After simulating the process of quantum communication, the output classical data is post-processed via procedures of parameter estimation, error correction, and privacy amplification. In particular, we analyze the high signal-to-noise regime which requires the use of high-rate (non-binary) low-density parity check codes. We implement all these steps in a Python-based library that allows one to investigate and optimize the protocol parameters to be used in practical experimental implementations of short-range CV-QKD.

categories: quant-ph physics.app-ph physics.comp-ph physics.data-an


title: Evaluating the Morphosyntactic Well-formedness of Generated Texts

abstract: Text generation systems are ubiquitous in natural language processing applications. However, evaluation of these systems remains a challenge, especially in multilingual settings. In this paper, we propose L'AMBRE -- a metric to evaluate the morphosyntactic well-formedness of text using its dependency parse and morphosyntactic rules of the language. We present a way to automatically extract various rules governing morphosyntax directly from dependency treebanks. To tackle the noisy outputs from text generation systems, we propose a simple methodology to train robust parsers. We show the effectiveness of our metric on the task of machine translation through a diachronic study of systems translating into morphologically-rich languages.


title: Continuous Weight Balancing

abstract: We propose a simple method by which to choose sample weights for problems with highly imbalanced or skewed traits. Rather than naively discretizing regression labels to find binned weights, we take a more principled approach -- we derive sample weights from the transfer function between an estimated source and specified target distributions. Our method outperforms both unweighted and discretely-weighted models on both regression and classification tasks. We also open-source our implementation of this method (https://github.com/Daniel-Wu/Continuous-Weight-Balancing) to the scientific community.


title: Inductive and Coinductive Topological Generation with Church's thesis and the Axiom of Choice

abstract: In this work we consider an extension MFcind of the Minimalist Foundation MF for predicative constructive mathematics with the addition of inductive and coinductive definitions sufficient to generate Sambin's Positive topologies, namely Martin-L\"of-Sambin formal topologies equipped with a Positivity relation (used to describe pointfree formal closed subsets). In particular the intensional level of MFcind, called mTTcind, is defined by extending with coinductive definitions another theory mTTind extending the intensional level mTT of MF with the sole addition of inductive definitions. In previous work we have shown that mTTind is consistent with Formal Church's Thesis CT and the Axiom of Choice AC via an interpretation in Aczel's CZF+REA. Our aim is to show the expectation that the addition of coinductive definitions to mTTind does not increase its consistency strength by reducing the consistency of mTTcind+CT+AC to the consistency of CZF+REA through various interpretations. We actually reach our goal in two ways. One way consists in first interpreting mTTcind+CT+AC in the theory extending CZF with the Union Regular Extension Axiom, REA_U, a strengthening of REA, and the Axiom of Relativized Dependent Choice, RDC. The theory CZF+REA_U+RDC is then interpreted in MLS*, a version of Martin-L\"of's type theory with Palmgren's superuniverse S. A last step consists in interpreting MLS* back into CZF+REA. The alternative way consists in first interpreting mTTcind+AC+CT directly in a version of Martin-L\"of's type theory with Palmgren's superuniverse extended with CT, which is then interpreted back to CZF+REA. A key benefit of the first way is that the theory CZF+REA_U+RDC also supports the intended set-theoretic interpretation of the extensional level of MFcind. Finally, all the theories considered, except mTTcind+AC+CT, are shown to be of the same proof-theoretic strength.


title: One-loop impact factor for the emission of two gluons

abstract: We consider one-loop five-point QCD amplitudes in next-to-multi-Regge kinematics, and evaluate the one-loop impact factor for the emission of two gluons. This is the last ingredient which is necessary to evaluate the gluon-jet impact factor at NNLO accuracy in $\as$. It is also the first instance in which loop-level QCD amplitudes are evaluated in next-to-multi-Regge kinematics, which requires to apply a different reggeisation ansatz to each colour-ordered amplitude.


title: Numerical Analysis of the Primordial Power Spectrum for (Small Field) Inflationary Potentials

abstract: We study small field models of inflation, which, against previous expectations, yield significant Gravitational Wave (GW) signal while reproducing other measured observable quantities in the Cosmic Microwave Background (CMB). We numerically study these, using previously published analytic works as general guidelines. We first discuss the framework necessary to understand model building, and some of its motivations. We review the slow-roll paradigm, derive the slow-roll parameters, and discuss different formulations thereof. We review the Lyth bound and its theoretical descendants, we outline the small/large field taxonomy and their characterization in the current nomenclature. We then present our models and the methods employed in their building and examination. We employ MCMC simulations to evaluate model likelihood and by process of marginalization extract the most probable coefficients for these inflationary potentials. An additional method employed is a multinomial fit, where we create a correspondence between coefficients and observables. This allows the use of observable values directly to yield the most likely coefficients. We compare the results of the two methods and evaluate the level of tuning required for these models. We discuss an apparent discrepancy between analytical approaches of evaluating Primordial Power Spectrum (PPS) observables and the precise numerical results in our models. We identify some of the sources of this discrepancy and remark on their meaning in the age of precision cosmology. Finally, we present the results of our study, for the most likely inflationary models with polynomial potentials of degree 5, and 6. We demonstrate our ability to produce potentials that yield GW with a tensor-to-scalar ratio r = 0.03. This is a realistic expectation of GW detection sensitivity in the near future.


title: Hidden symmetry between rotational tidal Love numbers of spinning neutron stars

abstract: The coupling between the angular momentum of a compact object and an external tidal field gives rise to the "rotational" tidal Love numbers, which affect the tidal deformability of a spinning self-gravitating body and enter the gravitational waveform of a binary inspiral at high post-Newtonian order. We provide numerical evidence for a surprising "hidden" symmetry among the rotational tidal Love numbers with opposite parities, which are associated to perturbations belonging to separate sectors. This symmetry, whose existence had been suggested on the basis of a Lagrangian description of the tidal interaction in a binary system, holds independently of the equation of state of the star.


title: Benchmarks for Deep Off-Policy Evaluation

abstract: Off-policy evaluation (OPE) holds the promise of being able to leverage large, offline datasets for both evaluating and selecting complex policies for decision making. The ability to learn offline is particularly important in many real-world domains, such as in healthcare, recommender systems, or robotics, where online data collection is an expensive and potentially dangerous process. Being able to accurately evaluate and select high-performing policies without requiring online interaction could yield significant benefits in safety, time, and cost for these applications. While many OPE methods have been proposed in recent years, comparing results between papers is difficult because currently there is a lack of a comprehensive and unified benchmark, and measuring algorithmic progress has been challenging due to the lack of difficult evaluation tasks. In order to address this gap, we present a collection of policies that in conjunction with existing offline datasets can be used for benchmarking off-policy evaluation. Our tasks include a range of challenging high-dimensional continuous control problems, with wide selections of datasets and policies for performing policy selection. The goal of our benchmark is to provide a standardized measure of progress that is motivated from a set of principles designed to challenge and test the limits of existing OPE methods. We perform an evaluation of state-of-the-art algorithms and provide open-source access to our data and code to foster future research in this area.


title: Rectification-based Knowledge Retention for Continual Learning

abstract: Deep learning models suffer from catastrophic forgetting when trained in an incremental learning setting. In this work, we propose a novel approach to address the task incremental learning problem, which involves training a model on new tasks that arrive in an incremental manner. The task incremental learning problem becomes even more challenging when the test set contains classes that are not part of the train set, i.e., a task incremental generalized zero-shot learning problem. Our approach can be used in both the zero-shot and non zero-shot task incremental learning settings. Our proposed method uses weight rectifications and affine transformations in order to adapt the model to different tasks that arrive sequentially. Specifically, we adapt the network weights to work for new tasks by "rectifying" the weights learned from the previous task. We learn these weight rectifications using very few parameters. We additionally learn affine transformations on the outputs generated by the network in order to better adapt them for the new task. We perform experiments on several datasets in both zero-shot and non zero-shot task incremental learning settings and empirically show that our approach achieves state-of-the-art results. Specifically, our approach outperforms the state-of-the-art non zero-shot task incremental learning method by over 5% on the CIFAR-100 dataset. Our approach also significantly outperforms the state-of-the-art task incremental generalized zero-shot learning method by absolute margins of 6.91% and 6.33% for the AWA1 and CUB datasets, respectively. We validate our approach using various ablation studies.


title: Gamma-convergence of fractional Gaussian perimeter

abstract: We prove the $\Gamma$-convergence of the renormalised fractional Gaussian $s$-perimeter to the Gaussian perimeter as $s\to 1^-$. Our definition of fractional perimeter comes from that of the fractional powers of Ornstein-Uhlenbeck operator given via Bochner subordination formula. As a typical feature of the Gaussian setting, the constant appearing in front of the $\Gamma$-limit does not depend on the dimension.


title: Alfven wave heating in partially ionized thin threads of solar prominences

abstract: There is observational evidence of the presence of small-amplitude transverse magnetohydrodynamic (MHD) waves with a wide range of frequencies in the threads of solar prominences. It is believed that the waves are driven at the photosphere and propagate along the magnetic field lines up to prominences suspended in the corona. The dissipation of MHD wave energy in the partially ionized prominence plasma is a heating mechanism whose relevance needs to be explored. Here we consider a simple 1D model for a non-uniform thin thread and investigate the heating associated with dissipation of Alfven waves. The model assumes an ad hoc density profile and a uniform pressure, while the temperature and ionization degree are self-consistently computed considering either LTE or non-LTE approximations for the hydrogen ionization. A broadband driver for Alfven waves is placed at one end of the magnetic field line, representing photospheric excitation. The Alfvenic perturbations along the thread are obtained by solving the linearized MHD equations for a partially ionized plasma in the single-fluid approximation.We find that wave heating in the partially ionized part of the thread is significant enough to compensate for energy losses due to radiative cooling. A greater amount of heating is found in the LTE case because the ionization degree for core prominence temperatures is lower than that in the non-LTE approximation. This results in a greater level of dissipation due to ambipolar diffusion in the LTE case. Conversely, in the hot coronal part of the model, the plasma is fully ionized and wave heating is negligible. The results of this simple model suggest that MHD wave heating can be relevant for the energy balance in prominences. Further studies based on more elaborate models are required.


title: Gap-tunable of Tunneling Time in Graphene Magnetic Barrier

abstract: We study the tunneling time of Dirac fermions in graphene magnetic barrier through an electrostatic potential and a mass term. This latter generates an energy gap in the spectrum and therefore affects the proprieties of tunneling of the system. For clarification, we first start by deriving the eigenspinors solutions of Dirac equation and second connect them to the incident, reflected and transmitted beam waves. This connection allows us to obtain the corresponding phases shifts and consequently compute the group delay time in transmission and reflection. Our numerical results show that the group delay time depends strongly on the energy gap in the tunneling process through single barrier. Moreover, we find that the group approaches unity at some critical value of the energy gap and becomes independent to the strengths of involved physical parameters.


title: Taking the temperature of a pure quantum state

abstract: Temperature is a deceptively simple concept that still raises deep questions at the forefront of quantum physics research. The observation of thermalisation in completely isolated quantum systems, such as cold-atom quantum simulators, implies that a temperature can be assigned even to individual, pure quantum states. Here, we propose a scheme to measure the temperature of such pure states through quantum interference. Our proposal involves interferometry of an auxiliary qubit probe, which is prepared in a superposition state and subsequently decoheres due to weak coupling with a closed, thermalised many-body system. Using only a few basic assumptions about chaotic quantum systems -- namely, the eigenstate thermalisation hypothesis and the emergence of hydrodynamics at long times -- we show that the qubit undergoes pure exponential decoherence at a rate that depends on the temperature of its surroundings. We verify our predictions by numerical experiments on a quantum spin chain that thermalises after absorbing energy from a periodic drive. Our work provides a general method to measure the temperature of isolated, strongly interacting systems under minimal assumptions.

categories: quant-ph cond-mat.mes-hall cond-mat.quant-gas cond-mat.stat-mech hep-th


title: Reducing the conjugacy problem for relatively hyperbolic automorphisms to peripheral components

abstract: We use relative hyperbolicity of mapping tori and Dehn fillings of relatively hyperbolic groups to solve the conjugacy problem between certain outer automorphisms. We reduce this problem to algorithmic problems entirely expressed in terms of the parabolic subgroups of the mapping tori. As an immediate application, we solve the conjugacy problem for outer automorphisms of free groups, whose polynomial part is piecewise inner. This proposes a path toward a full solution to the conjugacy problem for $\mathrm{Out}(F_n)$.


title: A Sublattice Phase-Field Model for Direct CALPHAD Database Coupling

abstract: The phase-field method has been established as a de facto standard for simulating the microstructural evolution of materials. In quantitative modeling the assessment and compilation of thermodynamic/kinetic data is largely dominated by the CALPHAD approach, which has produced a large set of experimentally and computationally generated Gibbs free energy and atomic mobility data in a standardized format: the thermodynamic database (TDB) file format. Harnessing this data for the purpose of phase-field modeling is an ongoing effort encompassing a wide variety of approaches. In this paper, we aim to directly link CALPHAD data to the phase-field method, without intermediate fitting or interpolation steps. We introduce a model based on the Kim-Kim-Suzuki (KKS) approach. This model includes sublattice site fractions and can directly utilize data from TDB files. Using this approach, we demonstrate the model on the U-Zr and Mo-Ni-Re systems.


title: VSS: A Storage System for Video Analytics [Technical Report]

abstract: We present a new video storage system (VSS) designed to decouple high-level video operations from the low-level details required to store and efficiently retrieve video data. VSS is designed to be the storage subsystem of a video data management system (VDBMS) and is responsible for: (1) transparently and automatically arranging the data on disk in an efficient, granular format; (2) caching frequently-retrieved regions in the most useful formats; and (3) eliminating redundancies found in videos captured from multiple cameras with overlapping fields of view. Our results suggest that VSS can improve VDBMS read performance by up to 54%, reduce storage costs by up to 45%, and enable developers to focus on application logic rather than video storage and retrieval.


title: Unsupervised Disentanglement of Linear-Encoded Facial Semantics

abstract: We propose a method to disentangle linear-encoded facial semantics from StyleGAN without external supervision. The method derives from linear regression and sparse representation learning concepts to make the disentangled latent representations easily interpreted as well. We start by coupling StyleGAN with a stabilized 3D deformable facial reconstruction method to decompose single-view GAN generations into multiple semantics. Latent representations are then extracted to capture interpretable facial semantics. In this work, we make it possible to get rid of labels for disentangling meaningful facial semantics. Also, we demonstrate that the guided extrapolation along the disentangled representations can help with data augmentation, which sheds light on handling unbalanced data. Finally, we provide an analysis of our learned localized facial representations and illustrate that the semantic information is encoded, which surprisingly complies with human intuition. The overall unsupervised design brings more flexibility to representation learning in the wild.


title: Functional and spatial rewiring jointly generate convergent-divergent units in self-organizing networks

abstract: Self-organization through adaptive rewiring of random neural networks generates brain-like topologies comprising modular small-world structures with rich club effects, merely as the product of optimizing the network topology. In the nervous system, spatial organization is optimized no less by rewiring, through minimizing wiring distance and maximizing spatially aligned wiring layouts. We show that such spatial organization principles interact constructively with adaptive rewiring, contributing to establish the networks' connectedness and modular structures. We use an evolving neural network model with weighted and directed connections, in which neural traffic flow is based on consensus and advection dynamics, to show that wiring cost minimization supports adaptive rewiring in creating convergent-divergent unit structures. Convergent-divergent units consist of a convergent input-hub, connected to a divergent output-hub via subnetworks of intermediate nodes, which may function as the computational core of the unit. The prominence of minimizing wiring distance in the dynamic evolution of the network determines the extent to which the core is encapsulated from the rest of the network, i.e., the context-sensitivity of its computations. This corresponds to the central role convergent-divergent units play in establishing context-sensitivity in neuronal information processing.


title: Seasonal Contrast: Unsupervised Pre-Training from Uncurated Remote Sensing Data

abstract: Remote sensing and automatic earth monitoring are key to solve global-scale challenges such as disaster prevention, land use monitoring, or tackling climate change. Although there exist vast amounts of remote sensing data, most of it remains unlabeled and thus inaccessible for supervised learning algorithms. Transfer learning approaches can reduce the data requirements of deep learning algorithms. However, most of these methods are pre-trained on ImageNet and their generalization to remote sensing imagery is not guaranteed due to the domain gap. In this work, we propose Seasonal Contrast (SeCo), an effective pipeline to leverage unlabeled data for in-domain pre-training of remote sensing representations. The SeCo pipeline is composed of two parts. First, a principled procedure to gather large-scale, unlabeled and uncurated remote sensing datasets containing images from multiple Earth locations at different timestamps. Second, a self-supervised algorithm that takes advantage of time and position invariance to learn transferable representations for remote sensing applications. We empirically show that models trained with SeCo achieve better performance than their ImageNet pre-trained counterparts and state-of-the-art self-supervised learning methods on multiple downstream tasks. The datasets and models in SeCo will be made public to facilitate transfer learning and enable rapid progress in remote sensing applications.


title: The Intrinsic Communication in Power Systems: A New Perspective to Understand Synchronization Stability

abstract: Synchronization is an essential element in three-phase ac electric power systems. The large-scale integration of converter-interfaced resources leads to the power grid transformation from voltage-source-dominated to voltage-current-source-composite, which also raises new challenges to model and analyze the system synchronization. In this article, we present the intrinsic analogy of a power system to a communication system, which is here called power-communication isomorphism. Based on this isomorphism, we revisit power system synchronization stability from a communication perspective and thereby establish a theory that unifies the synchronization dynamics of heterogeneous power apparatuses. In particular, the proposed theory is used to interpret and analyze, for example, the different power control speeds of power apparatuses; the role of network line dynamics in synchronization analysis; and the synchronization capability of current-source-dominated (e.g., grid-following-converter-dominated) grids. The findings are verified on the IEEE 68-bus test system.


title: Binarized Neural Networks for Resource-Constrained On-Device Gait Identification

abstract: User authentication through gait analysis is a promising application of discriminative neural networks -- particularly due to the ubiquity of the primary sources of gait accelerometry, in-pocket cellphones. However, conventional machine learning models are often too large and computationally expensive to enable inference on low-resource mobile devices. We propose that binarized neural networks can act as robust discriminators, maintaining both an acceptable level of accuracy while also dramatically decreasing memory requirements, thereby enabling on-device inference. To this end, we propose BiPedalNet, a compact CNN that nearly matches the state-of-the-art on the Padova gait dataset, with only 1/32 of the memory overhead.


title: The mean free path of ionizing photons at 5 < z < 6: evidence for rapid evolution near reionization

abstract: The mean free path of ionizing photons, $\lambda_{\rm mfp}$, is a key factor in the photoionization of the intergalactic medium (IGM). At $z \gtrsim 5$, however, $\lambda_{\rm mfp}$ may be short enough that measurements towards QSOs are biased by the QSO proximity effect. We present new direct measurements of $\lambda_{\rm mfp}$ that address this bias and extend up to $z \sim 6$ for the first time. Our measurements at $z \sim 5$ are based on data from the Giant Gemini GMOS survey and new Keck LRIS observations of low-luminosity QSOs. At $z \sim 6$ we use QSO spectra from Keck ESI and VLT X-Shooter. We measure $\lambda_{\rm mfp} = 9.09^{+1.62}_{-1.28}$ proper Mpc and $0.75^{+0.65}_{-0.45}$ proper Mpc (68% confidence) at $z = 5.1$ and 6.0, respectively. The results at $z = 5.1$ are consistent with existing measurements, suggesting that bias from the proximity effect is minor at this redshift. At $z = 6.0$, however, we find that neglecting the proximity effect biases the result high by a factor of two or more. Our measurement at $z = 6.0$ falls well below extrapolations from lower redshifts, indicating rapid evolution in $\lambda_{\rm mfp}$ over $5 < z < 6$. This evolution disfavors models in which reionization ended early enough that the IGM had time to fully relax hydrodynamically by $z = 6$, but is qualitatively consistent with models wherein reionization completed at $z = 6$ or even significantly later. Our mean free path results are most consistent with late reionization models wherein the IGM is still 20% neutral at $z=6$, although our measurement at $z = 6.0$ is even lower than these models prefer.


title: A Robust Stackelberg Game for Cyber-Security Investment in Networked Control Systems

abstract: We present a resource-planning game for cyber-security of networked control systems (NCS). The NCS is assumed to be operating in closed-loop using a linear state-feedback $\mathcal{H}_2$ controller. A zero-sum, two-player Stackelberg game (SG) is developed between an attacker and a defender for this NCS. The attacker aims to disable communication of selected nodes and thereby render the feedback gain matrix to be sparse, leading to degradation of closed-loop performance, while the defender aims to prevent this loss by investing in the protection of targeted nodes. Both players trade their $\mathcal{H}_2$-performance objectives for the costs of their actions. The standard backward induction method is modified to determine a cost-based Stackelberg equilibrium (CBSE) that saves the players' costs without degrading the control performance. We analyze the dependency of a CBSE on the relative budgets of the players as well as on the node "importance" order. Moreover, a robust-defense method is developed for the realistic case when the defender is not informed about the attacker's resources. The proposed algorithms are validated using examples from wide-area control of electric power systems. It is demonstrated that reliable and robust defense is feasible unless the defender's resources are severely limited relative to the attacker's resources. We also show that the proposed methods are robust to time-varying model uncertainties and thus are suitable for long-term security investment in realistic NCSs. Finally, we employ computationally efficient genetic algorithms (GA) to compute the optimal strategies of the attacker and the defender in realistic large power systems.


title: Geodesic deviation equation in Brans-Dicke theory in arbitrary dimensions

abstract: In this paper, we study the geodesic deviation equation (GDE) within the context of the Brans-Dicke (BD) theory in $D$ dimensions. Then, we restrict our attention to the GDE for the fundamental observers and null vector field past directed. Concerning the latter, in order to apply the retrieved GDE for the cosmological models as well as examining exact solutions, we study some cosmological models in general relativity (GR) and in the context of the BD theory. For the BD framework, we provide two appropriate settings for dynamical system, which can be considered as the most extended case scenarios with respect to those studied previously. Then, we investigate a few well-known BD cosmological models, for which we present the modified and corrected exact solutions. Moreover, we show that the retrieved formalism of the GDE can also be properly applicable for the modified BD theory (MBDT), whose matter and potential emerge from the geometry. For all of our herein cosmological models, we investigate the energy conditions and depict the behavior of the deviation vector and observer area distance. We demonstrate that why the MBDT can be considered as a more appropriate candidate to describe the universe in accordance with the observational data as well as from theoretical viewpoint.


title: Tracking Knowledge Propagation Across Wikipedia Languages

abstract: In this paper, we present a dataset of inter-language knowledge propagation in Wikipedia. Covering the entire 309 language editions and 33M articles, the dataset aims to track the full propagation history of Wikipedia concepts, and allow follow up research on building predictive models of them. For this purpose, we align all the Wikipedia articles in a language-agnostic manner according to the concept they cover, which results in 13M propagation instances. To the best of our knowledge, this dataset is the first to explore the full inter-language propagation at a large scale. Together with the dataset, a holistic overview of the propagation and key insights about the underlying structural factors are provided to aid future research. For example, we find that although long cascades are unusual, the propagation tends to continue further once it reaches more than four language editions. We also find that the size of language editions is associated with the speed of propagation. We believe the dataset not only contributes to the prior literature on Wikipedia growth but also enables new use cases such as edit recommendation for addressing knowledge gaps, detection of disinformation, and cultural relationship analysis.


title: B\'ezoutians and the $\mathbb{A}^1$-degree

abstract: We prove that both the local and global $\mathbb{A}^1$-degree of an endomorphism of affine space can be computed in terms of the multivariate B\'ezoutian. In particular, we show that the B\'ezoutian bilinear form, the Scheja--Storch form, and the $\mathbb{A}^1$-degree for complete intersections are isomorphic. Our global theorem generalizes Cazanave's theorem in the univariate case, and our local theorem generalizes Kass--Wickelgren's theorem on EKL forms and the local degree. This result provides an algebraic formula for local and global degrees in motivic homotopy theory.


title: TUSQ: Targeted High-Utility Sequence Querying

abstract: Significant efforts have been expended in the research and development of a database management system (DBMS) that has a wide range of applications for managing an enormous collection of multisource, heterogeneous, complex, or growing data. Besides the primary function (i.e., create, delete, and update), a practical and impeccable DBMS can interact with users through information selection, that is, querying with their targets. Previous querying algorithms, such as frequent itemset querying and sequential pattern querying (SPQ) have focused on the measurement of frequency, which does not involve the concept of utility, which is helpful for users to discover more informative patterns. To apply the querying technology for wider applications, we incorporate utility into target-oriented SPQ and formulate the task of targeted utility-oriented sequence querying. To address the proposed problem, we develop a novel algorithm, namely targeted high-utility sequence querying (TUSQ), based on two novel upper bounds suffix remain utility and terminated descendants utility as well as a vertical Last Instance Table structure. For further efficiency, TUSQ relies on a projection technology utilizing a compact data structure called the targeted chain. An extensive experimental study conducted on several real and synthetic datasets shows that the proposed algorithm outperformed the designed baseline algorithm in terms of runtime, memory consumption, and candidate filtering.


title: Probing the atmospheric precipitable water vapor with SOFIA, Part I, Measurements of the water vapor overburden with FIFI-LS

abstract: We report on the measurements of telluric water vapor made with the instrument FIFI-LS on SOFIA. Since November 2018, FIFI-LS has measured the water vapor overburden with the same measurement setup on each science flight with about 10 data points throughout the flight. This created a large sample of 469 measurements at different locations, flight altitudes and seasons. The paper describes the measurement principle in detail and provides some trend analysis on the 3 parameters. This presents the first systematic analysis with SOFIA based on in situ observations.


title: HAD-Net: A Hierarchical Adversarial Knowledge Distillation Network for Improved Enhanced Tumour Segmentation Without Post-Contrast Images

abstract: Segmentation of enhancing tumours or lesions from MRI is important for detecting new disease activity in many clinical contexts. However, accurate segmentation requires the inclusion of medical images (e.g., T1 post contrast MRI) acquired after injecting patients with a contrast agent (e.g., Gadolinium), a process no longer thought to be safe. Although a number of modality-agnostic segmentation networks have been developed over the past few years, they have been met with limited success in the context of enhancing pathology segmentation. In this work, we present HAD-Net, a novel offline adversarial knowledge distillation (KD) technique, whereby a pre-trained teacher segmentation network, with access to all MRI sequences, teaches a student network, via hierarchical adversarial training, to better overcome the large domain shift presented when crucial images are absent during inference. In particular, we apply HAD-Net to the challenging task of enhancing tumour segmentation when access to post-contrast imaging is not available. The proposed network is trained and tested on the BraTS 2019 brain tumour segmentation challenge dataset, where it achieves performance improvements in the ranges of 16% - 26% over (a) recent modality-agnostic segmentation methods (U-HeMIS, U-HVED), (b) KD-Net adapted to this problem, (c) the pre-trained student network and (d) a non-hierarchical version of the network (AD-Net), in terms of Dice scores for enhancing tumour (ET). The network also shows improvements in tumour core (TC) Dice scores. Finally, the network outperforms both the baseline student network and AD-Net in terms of uncertainty quantification for enhancing tumour segmentation based on the BraTs 2019 uncertainty challenge metrics. Our code is publicly available at: https://github.com/SaverioVad/HAD_Net


title: Fair and Distributed Dynamic Optimal Transport for Resource Allocation over Networks

abstract: Optimal transport is a framework that facilitates the most efficient allocation of a limited amount of resources. However, the most efficient allocation scheme does not necessarily preserve the most fairness. In this paper, we establish a framework which explicitly considers the fairness of dynamic resource allocation over a network with heterogeneous participants. As computing the transport strategy in a centralized fashion requires significant computational resources, it is imperative to develop computationally light algorithm that can be applied to large scale problems. To this end, we develop a fully distributed algorithm for fair and dynamic optimal transport with provable convergence using alternating method of multipliers. In the designed algorithm, each corresponding pair of resource supplier and receiver compute their own solutions and update the transport schemes through negotiation iteratively which do not require a central planner. The distributed algorithm can yield a fair and efficient resource allocation mechanism over a network. We corroborate the obtained results through case studies.


title: On the early-time behavior of quantum subharmonic generation

abstract: A few years ago Avetissian {\it et al.} \cite{Avetissian2014,Avetissian2015} discovered that the exponential growth rate of the stimulated annihilation photons from a singlet positronium Bose-Einstein condensate should be proportional to the square root of the positronium number density, not to the number density itself. In order to elucidate this surprising result obtained via a field-theoretical analysis, we point out that the basic physics involved is the same as that of resonant subharmonic transitions between two quantum oscillators. Using this model, we show that nonlinearities of the type discovered by Avetissian {\it et al.} are not unique to positronium and in fact will be encountered in a wide range of systems that can be modeled as nonlinearly coupled quantum oscillators.


title: Speed Up Zig-Zag

abstract: Zig-Zag is Piecewise Deterministic Markov Process, efficiently used for simulation in an MCMC setting. As we show in this article, it fails to be exponentially ergodic on heavy tailed target distributions. We introduce an extension of the Zig-Zag process by allowing the process to move with a non-constant speed function $s$, depending on the current state of the process. We call this process Speed Up Zig-Zag (SUZZ). We provide conditions that guarantee stability properties for the SUZZ process, including non-explosivity, exponential ergodicity in heavy tailed targets and central limit theorem. Interestingly, we find that using speed functions that induce explosive deterministic dynamics may lead to stable algorithms that can even mix faster. We further discuss the choice of an efficient speed function by providing an efficiency criterion for the one-dimensional process and we support our findings with simulation results.


title: Deconfinement of classical Yang-Mills color fields in a disorder potential

abstract: We study numerically and analytically the behavior of classical Yang-Mills color fields in a random one-dimensional potential described by the Anderson model with disorder. Above a certain threshold the nonlinear interactions of Yang-Mills fields lead to chaos and deconfinement of color wavepackets with their subdiffusive spreading in space. The algebraic exponent of the second moment growth in time is found to be in a range of 0.3 to 0.4. Below the threshold color wavepackets remain confined even if a very slow spreading at very long times is not excluded due to subtle nonlinear effects and the Arnold diffusion for the case when initially color packets are located in a close vicinity. In a case of large initial separation of color wavepackets they remain well confined and localized in space. We also present comparison with the behavior of the one-component field model of discrete Anderson nonlinear Schr\"odinger equation with disorder.


title: Weak-strong uniqueness property for models of compressible viscous fluids near vacuum

abstract: We extend the weak-strong uniqueness principle to general models of compressible viscous fluids near/on the vacuum. In particular, the physically relevant case of positive density with polynomial decay at infinity is considered.


title: Mapping the Electromagnetic Fields of Heavy-Ion Collisions with the Breit-Wheeler Process

abstract: Ultra-relativistic heavy-ion collisions are expected to produce the strongest electromagnetic fields in the known Universe. These highly-Lorentz contracted fields can manifest themselves as linearly polarized quasi-real photons that can interact via the Breit-Wheeler process to produce lepton anti-lepton pairs. The energy and momentum distribution of the produced dileptons carry information about the strength and spatial distribution of the colliding fields. Recently it has been demonstrated that photons from these fields can interact even in heavy-ion collisions with hadronic overlap, providing a purely electromagnetic probe of the produced medium. In this review we discuss the recent theoretical progress and experimental advances for mapping the ultra-strong electromagnetic fields produced in heavy-ion collisions via measurement of the Breit-Wheeler process.


title: Text Classification Using Hybrid Machine Learning Algorithms on Big Data

abstract: Recently, there are unprecedented data growth originating from different online platforms which contribute to big data in terms of volume, velocity, variety and veracity (4Vs). Given this nature of big data which is unstructured, performing analytics to extract meaningful information is currently a great challenge to big data analytics. Collecting and analyzing unstructured textual data allows decision makers to study the escalation of comments/posts on our social media platforms. Hence, there is need for automatic big data analysis to overcome the noise and the non-reliability of these unstructured dataset from the digital media platforms. However, current machine learning algorithms used are performance driven focusing on the classification/prediction accuracy based on known properties learned from the training samples. With the learning task in a large dataset, most machine learning models are known to require high computational cost which eventually leads to computational complexity. In this work, two supervised machine learning algorithms are combined with text mining techniques to produce a hybrid model which consists of Na\"ive Bayes and support vector machines (SVM). This is to increase the efficiency and accuracy of the results obtained and also to reduce the computational cost and complexity. The system also provides an open platform where a group of persons with a common interest can share their comments/messages and these comments classified automatically as legal or illegal. This improves the quality of conversation among users. The hybrid model was developed using WEKA tools and Java programming language. The result shows that the hybrid model gave 96.76% accuracy as against the 61.45% and 69.21% of the Na\"ive Bayes and SVM models respectively.


title: Z2 topological order and first-order quantum phase transitions in systems with combinatorial gauge symmetry

abstract: We study a generalization of the two-dimensional transverse-field Ising model, combining both ferromagnetic and antiferromagnetic two-body interactions, that hosts exact global and local Z2 gauge symmetries. Using exact diagonalization and stochastic series expansion quantum Monte Carlo methods, we confirm the existence of the topological phase in line with previous theoretical predictions. Our simulation results show that the transition between the confined topological phase and the deconfined paramagnetic phase is of first-order, in contrast to the conventional Z2 lattice gauge model in which the transition maps onto that of the standard Ising model and is continuous. We further generalize the model by replacing the transverse field on the gauge spins with a ferromagnetic XX interaction while keeping the local gauge symmetry intact. We find that the Z2 topological phase remains stable, while the paramagnetic phase is replaced by a ferromagnetic phase. The topological-ferromagnetic quantum phase transition is also of first-order. For both models, we discuss the low-energy spinon and vison excitations of the topological phase and their avoided level crossings associated with the first-order quantum phase transitions.


title: Deformed relativistic kinematics on curved spacetime -- a geometric approach

abstract: Deformed relativistic kinematics have been considered as a way to capture residual effects of quantum gravity. It has been shown that they can be understood geometrically in terms of a curved momentum space on a flat spacetime. In this article we present a systematic analysis under which conditions and how deformed relativistic kinematics, encoded in a momentum space metric on flat spacetime, can be lifted to curved spacetimes in terms of a self-consistent cotangent bundle geometry, which leads to purely geometric, geodesic motion of freely falling point particles. We comment how this construction is connected to, and offers a new perspective on, non-commutative spacetimes. From geometric consistency conditions we find that momentum space metrics can be consistently lifted to curved spacetimes if they either lead to a dispersion relation which is homogeneous in the momenta, or, if they satisfy a specific symmetry constraint. The latter is relevant for the momentum space metrics encoding the most studied deformed relativistic kinematics. For these, the constraint can only be satisfied in a momentum space basis in which the momentum space metric is invariant under linear local Lorentz transformations. We discuss how this result can be interpreted and the consequences of relaxing some conditions and principles of the construction from which we started.


title: Purely arithmetic PDE's over a p-adic field I: delta-characters and delta-modular forms

abstract: A formalism of arithmetic partial differential equations (PDEs) is being developed in which one considers several arithmetic differentiations at one fixed prime. In this theory solutions can be defined in algebraically closed p-adic fields. As an application we show that for at least two arithmetic directions every elliptic curve possesses a non-zero arithmetic PDE Manin map of order 1; such maps do not exist in the arithmetic ODE case. Similarly we construct and study "genuinely PDE" differential modular forms. As further applications we derive a Theorem of the Kernel and a Reciprocity Theorem for arithmetic PDE Manin maps and also a finiteness Diophantine result for modular parameterizations. We also prove structure results for the spaces of "PDE differential modular forms defined on the ordinary locus." We also produce a system of differential equations satisfied by our PDE modular forms based on Serre and Euler operators.


title: Gluons, Light and Heavy Quarks in the Instanton Vacuum

abstract: In this review we will concentrate on the nonperturbative effects in the properties of hadrons made from the light-heavy and heavy-heavy quarks in the framework of Instanton Liquid Model (ILM) of QCD vacuum. We briefly discuss the main features of ILM and its applicability in the heavy quark sector. The properties of gluonic systems, light and heavy quark correlators in the instanton background are also analyzed. Consideration of the both, perturbative and nonperturbative, gluon effects in the instanton background for the single heavy quark will lead to the mass shift due to the direct-instanton nonperturbative and ILM modified perturbative contributions, respectively. For the interacting heavy quark-antiquarks, the potential consists the direct instanton induced part and the one-gluon exchange (OGE) perturbative part. OGE interactions are screened at large distances due to the nonperturbative dynamics. We discuss the estimations of instanton contributions in the phenomenological Cornell type potential model. As related to the experimental data we discuss the charmonium properties and the role of instanton effects in their spectra and transitions. We discuss also the main features of heavy-light quarks systems in the ILM. As an example, it is considering the process of pions emission by exited heavy quarkonium states.


title: Learning Lipschitz Feedback Policies from Expert Demonstrations: Closed-Loop Guarantees, Generalization and Robustness

abstract: In this work, we propose a framework to learn feedback control policies with guarantees on closed-loop generalization and adversarial robustness. These policies are learned directly from expert demonstrations, contained in a dataset of state-control input pairs, without any prior knowledge of the task and system model. We use a Lipschitz-constrained loss minimization scheme to learn feedback policies with certified closed-loop robustness, wherein the Lipschitz constraint serves as a mechanism to tune the generalization performance and robustness to adversarial disturbances. Our analysis exploits the Lipschitz property to obtain closed-loop guarantees on generalization and robustness of the learned policies. In particular, we derive a finite sample bound on the policy learning error and establish robust closed-loop stability under the learned control policy. We also derive bounds on the closed-loop regret with respect to the expert policy and the deterioration of closed-loop performance under bounded (adversarial) disturbances to the state measurements. Numerical results validate our analysis and demonstrate the effectiveness of our robust feedback policy learning framework. Finally, our results suggest the existence of a potential tradeoff between nominal closed-loop performance and adversarial robustness, and that improvements in nominal closed-loop performance can only be made at the expense of robustness to adversarial perturbations.


title: Gaussian fluctuation for Gaussian Wishart matrices of overall correlation

abstract: In this note, we study the Gaussian fluctuations for the Wishart matrices $d^{-1}\mathcal{X}_{n, d}\mathcal{X}^{T}_{n, d}$, where $\mathcal{X}_{n, d}$ is a $n\times d$ random matrix whose entries are jointly Gaussian and correlated with row and column covariance functions given by $r$ and $s$ respectively such that $r(0)=s(0)=1$. Under the assumptions $s\in \ell^{4/3}(\mathbb{Z})$ and $\|r\|_{\ell^1(\mathbb{Z})}< \sqrt{6}/2$, we establish the $\sqrt{n^3/d}$ convergence rate for the Wasserstein distance between a normalization of $d^{-1}\mathcal{X}_{n, d}\mathcal{X}^{T}_{n, d}$ and the corresponding Gaussian ensemble. This rate is the same as the optimal one computed in \cite{JL15,BG16,BDER16} for the total variation distance, in the particular case where the Gaussian entries of $\mathcal{X}_{n, d}$ are independent. Similarly, we obtain the $\sqrt{n^{2p-1}/d}$ convergence rate for the Wasserstein distance in the setting of random $p$-tensors of overall correlation. Our analysis is based on the Malliavin-Stein approach.


title: Implicit completeness criterion in three-valued logic in terms of maximal classes

abstract: Implicit expressability was introduced by A.V. Kuznetsov in 1979 as generalization of functional expressability. Set of functions is called implicitly complete if any function has an implicit representation over this set. The system of all implicitly maximal classes in three-valued logic is described. The implicit completeness criterion is stated.


title: Design and Control of a Midair-Reconfigurable Quadcopter using Unactuated Hinges

abstract: A novel quadcopter capable of changing shape mid-flight is presented, allowing for operation in four configurations with the capability of sustained hover in three. This is accomplished without requiring actuators beyond the four motors typical of a quadcopter. Morphing is achieved through freely-rotating hinges that allow the vehicle arms to fold downwards by either reducing or reversing thrust forces. Constraints placed on the control inputs of the vehicle prevent the arms from folding or unfolding unexpectedly. This allows for the use of existing quadcopter controllers and trajectory generation algorithms with only minimal added complexity. For our experimental vehicle at hover, we find that these constraints result in a 36% reduction of the maximum yaw torque the vehicle can produce, but do not result in a reduction of the maximum thrust or roll and pitch torques. Experimental results show that, for a typical maneuver, the added limits have a negligible effect on trajectory tracking performance. Finally, the ability to change configurations is shown to enable the vehicle to traverse small passages, perch on hanging wires, and perform limited grasping tasks.


title: Neutron stars with a crossover equation of state

abstract: The question of whether quark matter exists in neutron stars is a long standing one. Generally one finds that a first order phase transition from baryons to quarks softens the equation of state so much that the star would collapse into a black hole. We consider a crossover equation of state, similar to the crossover that is found in lattice QCD studies at finite temperature and zero or small baryon chemical potentials. We find that with reasonable parameters it may be possible to support neutron stars up to about 2.2 solar masses. In that case 1 to 10% of the pressure would be contributed by quark matter in the central core of the highest mass stars.


title: Exploiting Invariance in Training Deep Neural Networks

abstract: Inspired by two basic mechanisms in animal visual systems, we introduce a feature transform technique that imposes invariance properties in the training of deep neural networks. The resulting algorithm requires less parameter tuning, trains well with an initial learning rate 1.0, and easily generalizes to different tasks. We enforce scale invariance with local statistics in the data to align similar samples at diverse scales. To accelerate convergence, we enforce a GL(n)-invariance property with global statistics extracted from a batch such that the gradient descent solution should remain invariant under basis change. Profiling analysis shows our proposed modifications takes 5% of the computations of the underlying convolution layer. Tested on convolutional networks and transformer networks, our proposed technique requires fewer iterations to train, surpasses all baselines by a large margin, seamlessly works on both small and large batch size training, and applies to different computer vision and language tasks.


title: Fitting the $Z b \bar b$ vertex in the two-Higgs-doublet model and in the three-Higgs-doublet model

abstract: We investigate the new contributions to the parameters $g_L$ and $g_R$ of the $Z b \bar b$ vertex in a multi-Higgs-doublet model (MHDM). We emphasize that those contributions generally worsen the fit of those parameters to the experimental data. We propose a solution to this problem, wherein $g_R$ has the opposite sign from the one predicted by the Standard Model; this solution, though, necessitates light scalars and large Yukawa couplings in the MHDM.


title: Detection of Rotational Variability in Floofy Objects at Optical Wavelengths

abstract: Phase resolved observations of planetary bodies allow us to understand the longitudinal and latitudinal variations that make each one unique. Rotational variations have been detected in several types of astronomical bodies beyond those of planetary mass, including asteroids, brown dwarfs, and stars. Unexpected rotational variations, such as those presented in this work, reminds us that the universe can be complicated, with more mysteries to uncover. In this work we present evidence for a new class of astronomical objects we identify as "floofy" with observational distinctions between several sub-types of these poorly understood objects. Using optical observations contributed by the community, we have identified rotational variation in several of these floofy objects, which suggests that they may have strong differences between their hemispheres, likely caused by differing reflectivity off their surfaces. Additional sub-types show no rotational variability suggesting a uniform distribution of reflective elements on the floofy object. While the work here is a promising step towards the categorization of floofy objects, further observations with more strictly defined limits on background light, illumination angles, and companion objects are necessary to develop a better understanding of the many remaining mysteries of these astronomical objects.


title: Comparison of wide-band vibrotactile and friction modulation surface gratings

abstract: This study seeks to understand conditions under which virtual gratings produced via vibrotaction and friction modulation are perceived as similar and to find physical origins in the results. To accomplish this, we developed two single-axis devices, one based on electroadhesion and one based on out-of-plane vibration. The two devices had identical touch surfaces, and the vibrotactile device used a novel closed-loop controller to achieve precise control of out-of-plane plate displacement under varying load conditions across a wide range of frequencies. A first study measured the perceptual intensity equivalence curve of gratings generated under electroadhesion and vibrotaction across the 20-400Hz frequency range. A second study assessed the perceptual similarity between two forms of skin excitation given the same driving frequency and same perceived intensity. Our results indicate that it is largely the out-of-plane velocity that predicts vibrotactile intensity relative to shear forces generated by friction modulation. A high degree of perceptual similarity between gratings generated through friction modulation and through vibrotaction is apparent and tends to scale with actuation frequency suggesting perceptual indifference to the manner of fingerpad actuation in the upper frequency range.


title: An optimal complexity spectral method for Navier--Stokes simulations in the ball

abstract: We develop a spectral method for solving the incompressible generalized Navier--Stokes equations in the ball with no-flux and prescribed slip boundary conditions. The algorithm achieves an optimal complexity per time step of $\mathcal{O}(N\log^2(N))$, where $N$ is the number of spatial degrees of freedom. The method relies on the poloidal-toroidal decomposition of solenoidal vector fields, the double Fourier sphere method, the Fourier and ultraspherical spectral method, and the spherical harmonics transform to decouple the Navier--Stokes equations and achieve the desired complexity and spectral accuracy.

categories: math.NA cs.NA physics.comp-ph physics.flu-dyn


title: Computational Fermi level engineering and doping-type conversion of Ga2O3 via three-step synthesis process

abstract: Ga2O3 is being actively explored for high-power and high-temperature electronics, deep-ultraviolet optoelectronics, and other applications. Efficient n-type doping of Ga2O3 has been achieved, but p-type doping faces fundamental obstacles due to compensation, deep acceptor levels, and the polaron transport mechanism of free holes. However, aside from achieving p-type conductivity, plenty of opportunity exists to engineer the position of the Fermi level for improved design of Ga2O3 based devices. We use first-principles defect theory and defect equilibrium calculations to simulate a 3-step growth-annealing-quench synthesis protocol for hydrogen assisted Mg doping in beta-Ga2O3, taking into account the gas phase equilibrium between H2, O2 and H2O, which determines the H chemical potential. We predict Ga2O3 doping-type conversion to a net p-type regime after growth under reducing conditions in the presence of H2 followed by O-rich annealing, which is a similar process to the Mg acceptor activation by H removal in GaN. For equilibrium annealing there is an optimal temperature that maximizes the Ga2O3 net acceptor density for a given Mg doping level, which is further increased in the non-equilibrium annealing scenario without re-equilibration. After quenching to operating temperature, the Ga2O3 Fermi level drops below mid-gap down to about +1.5 eV above the valence band maximum, creating a significant number of uncompensated neutral MgGa0 acceptors. The Fermi level reduction down to +1.5 eV and suppression of free electron density in this doping type converted (NA > ND) Ga2O3 material is of significance and impact for the design of Ga2O3 power electronics devices.


title: Frequency Estimation under Local Differential Privacy [Experiments, Analysis and Benchmarks]

abstract: Private collection of statistics from a large distributed population is an important problem, and has led to large scale deployments from several leading technology companies. The dominant approach requires each user to randomly perturb their input, leading to guarantees in the local differential privacy model. In this paper, we place the various approaches that have been suggested into a common framework, and perform an extensive series of experiments to understand the tradeoffs between different implementation choices. Our conclusion is that for the core problems of frequency estimation and heavy hitter identification, careful choice of algorithms can lead to very effective solutions that scale to millions of users


title: Tip-Induced $\beta$-Hydrogen Dissociation in an Alkyl Group Bound on Si(001)

abstract: Atomic-scale chemical modification of surface-adsorbed ethyl groups on Si(001) was induced and studied by means of scanning tunneling microscopy. Tunneling at sample bias > +1.5V leads to tip-induced C-H cleavage of a $\beta$-hydrogen of the covalently bound ethyl configuration. The reaction is characterized by the formation of an additional Si-H and a Si-C bond. The reaction probability shows a linear dependence on the tunneling current at 300 K; the reaction is largely suppressed at 50 K. The observed tip-induced surface reaction at room temperature is thus attributed to a one-electron excitation in combination with thermal activation.


title: The Metal-Poor Metallicity Distribution of the Ancient Milky Way

abstract: We present a low metallicity map of the Milky Way consisting of $\sim$111,000 giants with $-3.5 \lesssim$ [Fe/H] $\lesssim -$0.75, based on public photometry from the second data release of the SkyMapper survey. These stars extend out to $\sim$7kpc from the solar neighborhood and cover the main Galactic stellar populations, including the thick disk and the inner halo. Notably, this map can reliably differentiate metallicities down to [Fe/H] $\sim -3.0$, and thus provides an unprecedented view into the ancient, metal-poor Milky Way. Among the more metal-rich stars in our sample ([Fe/H] $> -2.0$), we recover a clear spatial dependence of decreasing mean metallicity as a function of scale height that maps onto the thick disk component of the Milky Way. When only considering the very metal-poor stars in our sample ([Fe/H] $< -$2), we recover no such spatial dependence in their mean metallicity out to a scale height of $|Z|\sim7$ kpc. We find that the metallicity distribution function (MDF) of the most metal-poor stars in our sample ($-3.0 <$ [Fe/H] $< -2.3$) is well fit with an exponential profile with a slope of $\Delta\log(N)/\Delta$[Fe/H] = 1.52$\pm$0.05, and shifts to $\Delta\log(N)/\Delta$[Fe/H] = 1.53$\pm$0.10 after accounting for target selection effects. For [Fe/H] $< -2.3$, the MDF is largely insensitive to scale height $|Z|$ out to $\sim5$kpc, showing that very and extremely metal-poor stars are in every galactic component.


title: Observations of the Quiet Sun During the Deepest Solar Minimum of the Past Century with Chandrayaan-2 XSM -- Elemental Abundances in the Quiescent Corona

abstract: Elements with low First Ionization Potential (FIP) are known to be three to four times more abundant in active region loops of the solar corona than in the photosphere. There have been observations suggesting that this observed "FIP bias" may be different in other parts of the solar corona and such observations are thus important in understanding the underlying mechanism. The Solar X-ray Monitor (XSM) on board the Chandrayaan-2 mission carried out spectroscopic observations of the Sun in soft X-rays during the 2019-20 solar minimum, considered to be the quietest solar minimum of the past century. These observations provided a unique opportunity to study soft X-ray spectra of the quiescent solar corona in the absence of any active regions. By modelling high resolution broadband X-ray spectra from XSM, we estimate the temperature and emission measure during periods of possibly the lowest solar X-ray intensity. We find that the derived parameters remain nearly constant over time with a temperature around 2 MK, suggesting the emission is dominated by X-ray Bright Points (XBPs). We also obtain the abundances of Mg, Al, and Si relative to H, and find that the FIP bias is ~2, lower than the values observed in active regions.


title: Observations of the Quiet Sun During the Deepest Solar Minimum of the Past Century with Chandrayaan-2 XSM -- Sub-A Class Microflares Outside Active Regions

abstract: Solar flares, with energies ranging over several orders of magnitude, result from impulsive release of energy due to magnetic reconnection in the corona. Barring a handful, almost all microflares observed in X-rays are associated with the solar active regions. Here we present, for the first time, a comprehensive analysis of a large sample of quiet Sun microflares observed in soft X-rays by the Solar X-ray Monitor (XSM) on board the Chandrayaan-2 mission during the 2019-20 solar minimum. A total of 98 microflares having peak flux below GOES A-level were observed by the XSM during observations spanning 76 days. By using the derived plasma temperature and emission measure of these events obtained by fitting the XSM spectra along with volume estimates from concurrent imaging observations in EUV with the Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA), we estimated their thermal energies to be ranging from 3e26 to 6e27 erg. We present the frequency distribution of the quiet Sun microflares with energy and discuss the implications of these observations of small-scale magnetic reconnection events outside active regions on coronal heating.


title: Dynamical Quantization of Contact Structures

abstract: We construct a dynamical quantization for contact manifolds in terms of a flat connection acting on a Hilbert tractor bundle. We show that this contact quantization, which is independent of the choice of contact form, can be obtained by quantizing the Reeb dynamics of an ambient strict contact manifold equivariantly with respect to an R+-action. The contact quantization further determines a certain contact tractor connection whose parallel sections determine a distinguished choice of Reeb dynamics and their quantization. This relationship relies on tractor constructions from parabolic geometries and mirrors the tight relationship between Einstein metrics and conformal geometries. Finally, we construct in detail the dynamical quantization of the unique tight contact structure on the 3-sphere, where the Holstein-Primakoff transformation makes a surprising appearance.

categories: math-ph hep-th math.DG math.MP math.SG quant-ph


title: Lifted MDS Codes over Finite Fields

abstract: MDS codes are elegant constructions in coding theory and have mode important applications in cryptography, network coding, distributed data storage, communication systems et. In this study, a method is given which MDS codes are lifted to a higher finite field. The presented method satisfies the protection of the distance and creating the MDS code over the $F_q$ by using MDS code over $F_p.$


title: An outer approximation algorithm for multi-objective mixed-integer linear and non-linear programming

abstract: In this paper, we present the first outer approximation algorithm for multi-objective mixed-integer linear programming problems with any number of objectives. The algorithm also works for certain classes of non-linear programming problems. It produces the non-dominated extreme points as well as the facets of the convex hull of these points. The algorithm relies on an oracle which solves single-objective weighted-sum problems and we show that the required number of oracle calls is polynomial in the number of facets of the convex hull of the non-dominated extreme points in the case of multiobjective mixed-integer programming (MOMILP). Thus, for MOMILP problems for which the weighted-sum problem is solvable in polynomial time, the facets can be computed with incremental-polynomial delay. From a practical perspective, the algorithm starts from a valid lower bound set for the non-dominated extreme points and iteratively improves it. Therefore it can be used in multi-objective branch-and-bound algorithms and still provide a valid bound set at any stage, even if interrupted before converging. Moreover, the oracle produces Pareto optimal solutions, which makes the algorithm also attractive from the primal side in a multi-objective branch-and-bound context. Finally, the oracle can also be called with any relaxation of the primal problem, and the obtained points and facets still provide a valid lower bound set. A computational study on a set of benchmark instances from the literature and new non-linear multi-objective instances is provided.


title: Exponential sums with multiplicative coefficients and applications

abstract: We show that if an exponential sum with multiplicative coefficients is large then the associated multiplicative function is "pretentious". This leads to applications in the circle method, and a natural interpretation of the local-global principle.


title: Revisiting Bayesian Optimization in the light of the COCO benchmark

abstract: It is commonly believed that Bayesian optimization (BO) algorithms are highly efficient for optimizing numerically costly functions. However, BO is not often compared to widely different alternatives, and is mostly tested on narrow sets of problems (multimodal, low-dimensional functions), which makes it difficult to assess where (or if) they actually achieve state-of-the-art performance. Moreover, several aspects in the design of these algorithms vary across implementations without a clear recommendation emerging from current practices, and many of these design choices are not substantiated by authoritative test campaigns. This article reports a large investigation about the effects on the performance of (Gaussian process based) BO of common and less common design choices. The experiments are carried out with the established COCO (COmparing Continuous Optimizers) software. It is found that a small initial budget, a quadratic trend, high-quality optimization of the acquisition criterion bring consistent progress. Using the GP mean as an occasional acquisition contributes to a negligible additional improvement. Warping degrades performance. The Mat\'ern 5/2 kernel is a good default but it may be surpassed by the exponential kernel on irregular functions. Overall, the best EGO variants are competitive or improve over state-of-the-art algorithms in dimensions less or equal to 5 for multimodal functions. The code developed for this study makes the new version (v2.1.1) of the R package DiceOptim available on CRAN. The structure of the experiments by function groups allows to define priorities for future research on Bayesian optimization.


title: Quasi-fibered boundary pseudodifferential operators

abstract: We develop a pseudodifferential calculus for differential operators associated to quasi-fibered boundary metrics (QFB metrics), a class of metrics including the quasi-asymptotically conical metrics (QAC metrics) of Degeratu-Mazzeo and the quasi-asymptotically locally Euclidean metrics (QALE metrics) of Joyce. Introducing various principal symbols, we introduce the notion of fully elliptic QFB operators and show that those are Fredholm when acting on QFB Sobolev spaces. For QAC metrics, we also develop a pseudodifferential calculus for the conformally related class of Qb metrics. We use these calculi to construct a parametrix for the Hodge-deRham operator of certain QFB metrics, allowing us to show that it is Fredholm on suitable Sobolev spaces and that the space of $L^2$ harmonic forms is finite dimensional. Our parametrix is obtained by inverting certain model operators at infinity, inversions that we achieve in part through a fine understanding of the low energy limit of the resolvent of the Hodge-deRham operator. Our parametrix also implies that $L^2$ harmonic forms decay faster at infinity than an arbitrary $L^2$ form, the extra decay being quantified in terms of a small negative power of the distance function. This decay of $L^2$ harmonic forms is used in a companion paper to study the $L^2$ cohomology of some $QFB$ metrics.


title: DAP: Detection-Aware Pre-training with Weak Supervision

abstract: This paper presents a detection-aware pre-training (DAP) approach, which leverages only weakly-labeled classification-style datasets (e.g., ImageNet) for pre-training, but is specifically tailored to benefit object detection tasks. In contrast to the widely used image classification-based pre-training (e.g., on ImageNet), which does not include any location-related training tasks, we transform a classification dataset into a detection dataset through a weakly supervised object localization method based on Class Activation Maps to directly pre-train a detector, making the pre-trained model location-aware and capable of predicting bounding boxes. We show that DAP can outperform the traditional classification pre-training in terms of both sample efficiency and convergence speed in downstream detection tasks including VOC and COCO. In particular, DAP boosts the detection accuracy by a large margin when the number of examples in the downstream task is small.


title: Robustness Certification for Point Cloud Models

abstract: The use of deep 3D point cloud models in safety-critical applications, such as autonomous driving, dictates the need to certify the robustness of these models to real-world transformations. This is technically challenging, as it requires a scalable verifier tailored to point cloud models that handles a wide range of semantic 3D transformations. In this work, we address this challenge and introduce 3DCertify, the first verifier able to certify the robustness of point cloud models. 3DCertify is based on two key insights: (i) a generic relaxation based on first-order Taylor approximations, applicable to any differentiable transformation, and (ii) a precise relaxation for global feature pooling, which is more complex than pointwise activations (e.g., ReLU or sigmoid) but commonly employed in point cloud models. We demonstrate the effectiveness of 3DCertify by performing an extensive evaluation on a wide range of 3D transformations (e.g., rotation, twisting) for both classification and part segmentation tasks. For example, we can certify robustness against rotations by $\pm$60{\deg} for 95.7% of point clouds, and our max pool relaxation increases certification by up to 15.6%.


title: New Algorithms for Discrete-Time Parameter Estimation

abstract: We propose two algorithms for discrete-time parameter estimation, one for time-varying parameters under persistent excitation (PE) condition, another for constant parameters under no PE condition. For the first algorithm, we show that in the presence of time-varying unknown parameters, the parameter estimation error converges uniformly to a compact set under conditions of persistent excitation, with the size of the compact set proportional to the time-variation of unknown parameters. Leveraging a projection operator, the second algorithm is shown to result in boundedness guarantees when the plant has constant unknown parameters. Simulations show better convergence results compared to recursive least squares (RLS) and comparable results to RLS with forgetting factor.


title: A Systematic Literature Review on Process-Aware Recommender Systems

abstract: Considering processes of a business in a recommender system is highly advantageous. Although most studies in the business process analysis domain are of descriptive and predictive nature, the feasibility of constructing a process-aware recommender system is assessed in a few works. One reason can be the lack of knowledge on process mining potential for recommendation problems. Therefore, this paper aims to identify and analyze the published studies on process-aware recommender system techniques in business process management and process mining domain. A systematic review was conducted on 33 academic articles published between 2008 and 2020 according to several aspects. In this regard, we provide a state-of-the-art review with critical details and researchers with a better perception of which path to pursue in this field. Moreover, based on a knowledge base and holistic perspective, we discuss some research gaps and open challenges in this field.


abstract: We develop new techniques to compute the weighted $L^2$-cohomology of quasi-fibered boundary metrics (QFB-metrics). Combined with the decay of $L^2$-harmonic forms obtained in a companion paper, this allows us to compute the reduced $L^2$-cohomology for various classes of QFB-metrics. Our results applies in particular to the Nakajima metric on the Hilbert scheme of $n$ points on $\mathbb{C}^2$, for which we can show that the Vafa-Witten conjecture holds. Using the compactification of the monopole moduli space announced by Fritzsch, the first author and Singer, we can also give a proof of the Sen conjecture for the monopole moduli space of magnetic charge 3.


title: Noisy Coherent Population Trapping: Applications to Noise Estimation and Qubit State Preparation

abstract: Coherent population trapping is a well-known quantum phenomenon in a driven $\Lambda$ system, with many applications across quantum optics. However, when a stochastic bath is present in addition to vacuum noise, the observed trapping is no longer perfect. Here we derive a time-convolutionless master equation describing the equilibration of the $\Lambda$ system in the presence of additional temporally correlated classical noise, with an unknown decay parameter. Our simulations show a one-to-one correspondence between the decay parameter and the depth of the characteristic dip in the photoluminescence spectrum, thereby enabling the unknown parameter to be estimated from the observed spectra. We apply our analysis to the problem of qubit state initialization in a $\Lambda$ system via dark states and show how the stochastic bath affects the fidelity of such initialization as a function of the desired dark-state amplitudes. We show that an optimum choice of Rabi frequencies is possible.


title: Typical periods or typical time steps? A multi-model analysis to determine the optimal temporal aggregation for energy system models

abstract: Energy system models are challenged by the need for high temporal and spatial resolutions in or-der to appropriately depict the increasing share of intermittent renewable energy sources, storage technologies, and the growing interconnectivity across energy sectors. This study evaluates methods for maintaining the computational viability of these models by ana-lyzing different temporal aggregation techniques that reduce the number of time steps in their in-put time series. Two commonly-employed approaches are the representation of time series by a subset of single (typical) time steps, or by groups of consecutive time steps (typical periods). We test these techniques for two different energy system models that are implemented using the Frame-work for Integrated Energy System Assessment (FINE) by benchmarking the optimization results based on aggregation to those of the fully resolved models and investigating whether the optimal aggregation method can, a priori, be determined based on the clustering indicators. The results reveal that typical time steps consistenly outperform typical days with respect to cluster-ing indicators, but do not lead to more accurate optimization results when applied to a model that takes numerous storage technologies into account. Although both aggregation techniques are ca-pable of coupling the aggregated time steps, typical days offer more options to depict storage oper-ations, whereas typical time steps are more effective for models that neglect time-linking con-straints. In summary, the adequate choice of aggregation methods strongly depends on the mathematical structure of the considered energy system optimization model, and a priori decisions of a sufficient temporal aggregation are only possible with good knowledge of this mathematical structure.


title: Boundary of the boundary for random walks on groups

abstract: We study fine structure related to finitely supported random walks on infinite finitely generated discrete groups, largely motivated by dimension group techniques. The unfaithful extreme harmonic functions (defined only on proper space-time cones), aka unfaithful pure traces, can be represented on systems of finite support, avoiding dead ends. This motivates properties of the random walk (WC) and of the group (SWC) which become of interest in their own right. While all abelian groups satisfy WC, the do not satisfy SWC; however some abelian by finite groups do satisfy the latter, and we characterize when this occurs. In general, we determine the maximal order ideals, aka, maximal proper space-time subcones of that generated by the group element $1$ at time zero), and show that the corresponding quotients are stationary simple dimension groups, and that all such can occur for the free group on two generators. We conclude with a case study of the discrete Heisenberg group, determining among other things, the pure traces (these are the unfaithful ones, not arising from characters).


title: Scalable adaptive cubic regularization methods

abstract: Adaptive cubic regularization (ARC) methods for unconstrained optimization compute steps from linear systems involving a shifted Hessian in the spirit of the Levenberg-Marquardt and trust-region methods. The standard approach consists in performing an iterative search for the shift akin to solving the secular equation in trust-region methods. Such search requires computing the Cholesky factorization of a tentative shifted Hessian at each iteration, which limits the size of problems that can be reasonably considered. We propose a scalable implementation of ARC named ARCqK in which we solve a set of shifted systems concurrently by way of an appropriate modification of the Lanczos formulation of the conjugate gradient (CG) method. At each iteration of ARCqK to solve a problem with n variables, a range of m << n shift parameters is selected. The computational overhead in CG beyond the Lanczos process is thirteen scalar operations to update five vectors of length m and two n-vector updates for each value of the shift. The CG variant only requires one Hessian-vector product and one dot product per iteration, independently of the number of shift parameters. Solves corresponding to inadequate shift parameters are interrupted early. All shifted systems are solved inexactly. Such modest cost makes our implementation scalable and appropriate for large-scale problems. We provide a new analysis of the inexact ARC method including its worst case evaluation complexity, global and asymptotic convergence. We describe our implementation and provide preliminary numerical observations that confirm that for problems of size at least 100, our implementation of ARCqK is more efficient than a classic Steihaug-Toint trust region method. Finally, we generalize our convergence results to inexact Hessians and nonlinear least-squares problems.


title: Stellar Metallicities from SkyMapper Photometry II: Precise photometric metallicities of $\sim$280,000 giant stars with [Fe/H] $< -0.75$ in the Milky Way

abstract: The Milky Way's metal-poor stars are nearby ancient objects that are used to study early chemical evolution and the assembly and structure of the Milky Way. Here we present reliable metallicities of $\sim280,000$ stars with $-3.75 \lesssim$ [Fe/H] $\lesssim -0.75$ down to $g=17$ derived using metallicity-sensitive photometry from the second data release (DR2) of the SkyMapper Southern Survey. We use the dependency of the flux through the SkyMapper $v$ filter on the strength of the Ca II K absorption features, in tandem with SkyMapper $u,g,i$ photometry, to derive photometric metallicities for these stars. We find that metallicities derived in this way compare well to metallicities derived in large-scale spectroscopic surveys, and use such comparisons to calibrate and quantify systematics as a function of location, reddening, and color. We find good agreement with metallicities from the APOGEE, LAMOST, and GALAH surveys, based on a standard deviation of $\sigma\sim0.25$dex of the residuals of our photometric metallicities with respect to metallicities from those surveys. We also compare our derived photometric metallicities to metallicities presented in a number of high-resolution spectroscopic studies to validate the low metallicity end ([Fe/H] $< -2.5$) of our photometric metallicity determinations. In such comparisons, we find the metallicities of stars with photometric [Fe/H] $< -2.5$ in our catalog show no significant offset and a scatter of $\sigma\sim$0.31dex level relative to those in high-resolution work when considering the cooler stars ($g-i > 0.65$) in our sample. We also present an expanded catalog containing photometric metallicities of $\sim720,000$ stars as a data table for further exploration of the metal-poor Milky Way.


title: Non-Existence and Existence of Shock Profiles in the Bemfica-Disconzi-Noronha Model

abstract: This note studies a four-field hyperbolic PDE model that was recently introduced by Bemfica, Disconzi, and Noronha for the pure radiation fluid with viscosity, and asks whether shock waves admit continuous profiles in this description. The model containing two free parameters mu, nu and being causal whenever one chooses (mu,nu) from a certain range C subset R2, this paper shows that for any choice of (mu,nu) in the interior of C, there is a dichotomy in so far as (i) shocks of sufficiently small amplitude admit profiles and (ii) certain other, thus necessarily non-small, shocks do not. This finding does not preclude the possibility that if one chooses (mu,nu) from a specific part S of the boundary of C, the dichotomy disappears and all shocks have profiles; the parameter set S corresponds to the "sharply causal" case, in which one of the characteristic speeds of the dissipation operator is the speed of light.


title: A genuinely natural information measure

abstract: The theoretical measuring of information was famously initiated by Shannon in his mathematical theory of communication, in which he proposed a now widely used quantity, the entropy, measured in bits. Yet, in the same paper, Shannon also chose to measure the information in continuous systems in nats, which differ from bits by the use of the natural rather than the binary logarithm. We point out that there is nothing natural about the choice of logarithm basis, rather it is arbitrary. We remedy this problematic state of affairs by proposing a genuinely natural measure of information, which we dub gnats. We show that gnats have many advantages in information theory, and propose to adopt the underlying methodology throughout science, arts and everyday life.


title: Rapidly-Declining Hostless Type Ia Supernova KSP-OT-201509b from the KMTNet Supernova Program: Transitional Nature and Constraint on $^{56}$Ni Distribution and Progenitor Type

abstract: We report the early discovery and multi-color ($BVI$) high-cadence light curve analyses of a rapidly-declining sub-Chandrasekhar Type Ia supernova KSP-OT-201509b (= AT2015cx) from the KMTNet Supernova Program. The Phillips parameter and color stretch parameter of KSP-OT-201509b (= AT2015cx) are $\Delta$$M_{B,15}$ $\simeq$ 1.62 mag and $s_{BV}$ $\simeq$ 0.54, respectively, at an inferred redshift of 0.072. These, together with other measured parameters (such as the strength of the secondary $I$-band peak, colors and luminosity), identify the source to be a rapidly-declining Type Ia of transitional nature that is closer to Branch Normal than 91bg-like. Its early light curve evolution and bolometric luminosity are consistent with those of homologously expanding ejecta powered by radioactive decay and a Type Ia SN explosion with 0.32 $\pm$ 0.01 $M_{\odot}$ of synthesized $^{56}$Ni mass, 0.84 $\pm$ 0.12 $M_{\odot}$ of ejecta mass and (0.61 $\pm$ 0.14) $\times$ 10$^{51}$ erg of ejecta kinetic energy. While its $B-V$ and $V-I$ colors evolve largely synchronously with the changes in the $I$-band light curve as found in other supernovae, we also find the presence of an early redward evolution in $V-I$ prior to --10 days since peak. The bolometric light curve of the source is compatible with a stratified $^{56}$Ni distribution extended to shallow layers of the exploding progenitor. Comparisons between the observed light curves and those predicted from ejecta-companion interactions clearly disfavor Roche Lobe-filling companion stars at large separation distances, thus supporting a double-degenerate scenario for its origin. The lack of any apparent host galaxy in our deep stack images reaching a sensitivity limit of $\sim$ 28 $\rm mag\;arcsec^{-2}$ makes KSP-OT-201509b a hostless Type Ia supernova and offers new insights into supernova host galaxy environments.


title: A probabilistic deep learning approach to automate the interpretation of multi-phase diffraction spectra

abstract: Autonomous synthesis and characterization of inorganic materials requires the automatic and accurate analysis of X-ray diffraction spectra. For this task, we designed a probabilistic deep learning algorithm to identify complex multi-phase mixtures. At the core of this algorithm lies an ensemble convolutional neural network trained on simulated diffraction spectra, which are systematically augmented with physics-informed perturbations to account for artifacts that can arise during experimental sample preparation and synthesis. Larger perturbations associated with off-stoichiometry are also captured by supplementing the training set with hypothetical solid solutions. Spectra containing mixtures of materials are analyzed with a newly developed branching algorithm that utilizes the probabilistic nature of the neural network to explore suspected mixtures and identify the set of phases that maximize confidence in the prediction. Our model is benchmarked on simulated and experimentally measured diffraction spectra, showing exceptional performance with accuracies exceeding those given by previously reported methods based on profile matching and deep learning. We envision that the algorithm presented here may be integrated in experimental workflows to facilitate the high-throughput and autonomous discovery of inorganic materials.