title: Evaluating the performance of personal, social, health-related, biomarker and genetic data for predicting an individuals future health using machine learning: A longitudinal analysis

abstract: As we gain access to a greater depth and range of health-related information about individuals, three questions arise: (1) Can we build better models to predict individual-level risk of ill health? (2) How much data do we need to effectively predict ill health? (3) Are new methods required to process the added complexity that new forms of data bring? The aim of the study is to apply a machine learning approach to identify the relative contribution of personal, social, health-related, biomarker and genetic data as predictors of future health in individuals. Using longitudinal data from 6830 individuals in the UK from Understanding Society (2010-12 to 2015-17), the study compares the predictive performance of five types of measures: personal (e.g. age, sex), social (e.g. occupation, education), health-related (e.g. body weight, grip strength), biomarker (e.g. cholesterol, hormones) and genetic single nucleotide polymorphisms (SNPs). The predicted outcome variable was limiting long-term illness one and five years from baseline. Two machine learning approaches were used to build predictive models: deep learning via neural networks and XGBoost (gradient boosting decision trees). Model fit was compared to traditional logistic regression models. Results found that health-related measures had the strongest prediction of future health status, with genetic data performing poorly. Machine learning models only offered marginal improvements in model accuracy when compared to logistic regression models, but also performed well on other metrics e.g. neural networks were best on AUC and XGBoost on precision. The study suggests that increasing complexity of data and methods does not necessarily translate to improved understanding of the determinants of health or performance of predictive models of ill health.


title: Four-point interfacial correlation functions in two dimensions. Exact results from field theory and numerical simulations

abstract: We derive exact analytic results for several four-point correlation functions for statistical models exhibiting phase separation in two-dimensions. Our theoretical results are then specialized to the Ising model on the two-dimensional strip and found to be in excellent agreement with high-precision Monte Carlo simulations.


title: Unified Spatio-Temporal Modeling for Traffic Forecasting using Graph Neural Network

abstract: Research in deep learning models to forecast traffic intensities has gained great attention in recent years due to their capability to capture the complex spatio-temporal relationships within the traffic data. However, most state-of-the-art approaches have designed spatial-only (e.g. Graph Neural Networks) and temporal-only (e.g. Recurrent Neural Networks) modules to separately extract spatial and temporal features. However, we argue that it is less effective to extract the complex spatio-temporal relationship with such factorized modules. Besides, most existing works predict the traffic intensity of a particular time interval only based on the traffic data of the previous one hour of that day. And thereby ignores the repetitive daily/weekly pattern that may exist in the last hour of data. Therefore, we propose a Unified Spatio-Temporal Graph Convolution Network (USTGCN) for traffic forecasting that performs both spatial and temporal aggregation through direct information propagation across different timestamp nodes with the help of spectral graph convolution on a spatio-temporal graph. Furthermore, it captures historical daily patterns in previous days and current-day patterns in current-day traffic data. Finally, we validate our work's effectiveness through experimental analysis, which shows that our model USTGCN can outperform state-of-the-art performances in three popular benchmark datasets from the Performance Measurement System (PeMS). Moreover, the training time is reduced significantly with our proposed USTGCN model.


title: Carrier concentrations and optical conductivity of a band-inverted semimetal in two and three dimensions

abstract: Here we study the single-particle, electronic transport and optical properties of a gapped system described by a simple two-band Hamiltonian with inverted valence bands. We analyze its properties in the three-dimensional (3D) and the two-dimensional (2D) case. The insulating phase changes into a metallic phase when the band gap is set to zero. The metallic phase in the 3D case is characterized by a nodal surface. This nodal surface is equivalent to a nodal ring in two dimensions. Within a simple theoretical framework, we calculate the density of state, the total and effective charge carrier concentration, the Hall concentration and the Hall coefficient, for both 2D and 3D cases. The main result is that the three concentrations always differ from one another in the present model. These concentrations can then be used to resolve the nature of the electronic ground state. Similarly, the optical conductivity is calculated and discussed for the insulating phase. We show that there are no optical excitations in the metallic phase. Finally, we compare the calculated optical conductivity with the rule-of-thumb derivation using the joint density of states.


title: Radio stars of the SKA

abstract: Radio emission from stars can be used, e.g., to study ionized winds or stellar flares. The radio emission is faint and studies have been limited to few objects. The Square Kilometer Array (SKA) will bring a survey ability to the topic of radio stars. In this paper we investigate what the SKA can detect, and what sensitivity will be required for deep surveys of the stellar Milky Way. We focus on the radio emission from OB stars, Be stars, flares from M dwarfs, and Ultra Compact HII regions. The stellar distribution in the Milky Way is simulated using the Besancon model, and various relations are used to predict their radio flux. We find that the full SKA will easily detect all UltraCompact HII regions. At the limit of 10 nJy at 5 GHz, the SKA can detect 1500 Be stars and 50 OB stars per square degree, out to several kpc. It can also detect flares from 4500 M dwarfs per square degree. At 100 nJy the numbers become about 8 times smaller. SKA surveys of the Galactic plane should be designed for high sensitivity. Deep imaging should consider the significant number of faint flares in the field, even outside the plane of the Milky Way.


title: Semantic Data Augmentation for End-to-End Mandarin Speech Recognition

abstract: End-to-end models have gradually become the preferred option for automatic speech recognition (ASR) applications. During the training of end-to-end ASR, data augmentation is a quite effective technique for regularizing the neural networks. This paper proposes a novel data augmentation technique based on semantic transposition of the transcriptions via syntax rules for end-to-end Mandarin ASR. Specifically, we first segment the transcriptions based on part-of-speech tags. Then transposition strategies, such as placing the object in front of the subject or swapping the subject and the object, are applied on the segmented sentences. Finally, the acoustic features corresponding to the transposed transcription are reassembled based on the audio-to-text forced-alignment produced by a pre-trained ASR system. The combination of original data and augmented one is used for training a new ASR system. The experiments are conducted on the Transformer[2] and Conformer[3] based ASR. The results show that the proposed method can give consistent performance gain to the system. Augmentation related issues, such as comparison of different strategies and ratios for data combination are also investigated.


title: An Open-Source Bayesian Atmospheric Radiative Transfer (BART) Code: I. Design, Tests, and Application to Exoplanet HD 189733 b

abstract: We present the open-source Bayesian Atmospheric Radiative Transfer (BART) retrieval package, which produces estimates and uncertainties for an atmosphere's thermal profile and chemical abundances from observations. Several BART components are also stand-alone packages, including the parallel Multi-Core Markov chain Monte Carlo (MC3), which implements several Bayesian samplers; a line-by-line radiative-transfer model, transit; a code that calculates Thermochemical Equilibrium Abundances, TEA; and a test suite for verifying radiative-transfer and retrieval codes, BARTTest. The codes are in Python and C. BART and TEA are under a Reproducible Research (RR) license, which requires reviewed-paper authors to publish a compendium of all inputs, codes, and outputs supporting the paper's scientific claims. BART and TEA produce the compendium's content. Otherwise, these codes are under permissive open-source terms, as are MC3 and BARTTest, for any purpose. This paper presents an overview of the code, BARTTest, and an application to eclipse data for exoplanet HD 189733 b. Appendices address RR methodology for accelerating science, a reporting checklist for retrieval papers, the spectral resolution required for synthetic tests, and a derivation of the effective sample size required to estimate any Bayesian posterior distribution to a given precision, which determines how many iterations to run. Paper II, by Cubillos et al., presents the underlying radiative-transfer scheme and an application to transit data for exoplanet HAT-P-11b. Paper III, by Blecic et al., discusses the initialization and post-processing routines, with an application to eclipse data for exoplanet WASP-43b. We invite the community to use and improve BART and its components at http://GitHub.com/ExOSPORTS/BART/.


title: Enumeration of Far-Apart Pairs by Decreasing Distance for Faster Hyperbolicity Computation

abstract: Hyperbolicity is a graph parameter which indicates how much the shortest-path distance metric of a graph deviates from a tree metric. It is used in various fields such as networking, security, and bioinformatics for the classification of complex networks, the design of routing schemes, and the analysis of graph algorithms. Despite recent progress, computing the hyperbolicity of a graph remains challenging. Indeed, the best known algorithm has time complexity $O(n^{3.69})$, which is prohibitive for large graphs, and the most efficient algorithms in practice have space complexity $O(n^2)$. Thus, time as well as space are bottlenecks for computing hyperbolicity. In this paper, we design a tool for enumerating all far-apart pairs of a graph by decreasing distances. A node pair $(u, v)$ of a graph is far-apart if both $v$ is a leaf of all shortest-path trees rooted at $u$ and $u$ is a leaf of all shortest-path trees rooted at $v$. This notion was previously used to drastically reduce the computation time for hyperbolicity in practice. However, it required the computation of the distance matrix to sort all pairs of nodes by decreasing distance, which requires an infeasible amount of memory already for medium-sized graphs. We present a new data structure that avoids this memory bottleneck in practice and for the first time enables computing the hyperbolicity of several large graphs that were far out-of-reach using previous algorithms. For some instances, we reduce the memory consumption by at least two orders of magnitude. Furthermore, we show that for many graphs, only a very small fraction of far-apart pairs have to be considered for the hyperbolicity computation, explaining this drastic reduction of memory. As iterating over far-apart pairs in decreasing order without storing them explicitly is a very general tool, we believe that our approach might also be relevant to other problems.


title: An Open-source Bayesian Atmospheric Radiative Transfer (BART) Code: II. The Transit Radiative-transfer Module and Retrieval of HAT-P-11b

abstract: This and companion papers by Harrington et al. and Blecic et al. present the Bayesian Atmospheric Radiative Transfer (BART) code, an open-source, open-development package to characterize extrasolar-planet atmospheres. BART combines a thermochemical equilibrium abundances (TEA), a radiative-transfer (Transit), and a Bayesian statistical (MC3) module to constrain atmospheric temperatures and molecular abundances for given spectroscopic observations. Here, we describe the Transit radiative-transfer package, an efficient line-by-line radiative-transfer C code for one-dimensional atmospheres, developed by P. Rojo and further modified by the UCF exoplanet group. This code produces transmission and hemisphere-integrated emission spectra. Transit handles line-by-line opacities from HITRAN, Partridge \& Schwenke ({\water}), Schwenke (TiO), and Plez (VO); and collision-induced absorption from Borysow, HITRAN, and ExoMol. Transit emission-spectra models agree with models from C. Morley (priv. comm.) within a few percent. We applied BART to the {\Spitzer} and {\Hubble} transit observations of the Neptune-sized planet HAT-P-11b. Our results generally agree with those from previous studies, constraining the {\water} abundance and finding an atmosphere enhanced in heavy elements. Different conclusions start to emerge when we make different assumptions from other studies. The BART source code and documentation are available at https://github.com/exosports/BART.


title: An open-source Bayesian atmospheric radiative transfer (BART) code: III. Initialization, atmospheric profile generator, post-processing routines, and application to exoplanet WASP-43b

abstract: This and companion papers by Harrington et al. 2021, submitted and Cubillos et al. 2021, submitted describe an open-source retrieval framework, Bayesian Atmospheric Radiative Transfer (BART), available to the community under the reproducible-research license via https://github.com/exosports/BART . BART is a radiative-transfer code (transit, https://github.com/exosports/transit , Rojo 2009, 2009ASPC..420..321R), initialized by the Thermochemical Equilibrium Abundances (TEA, https://github.com/dzesmin/TEA , Blecic et al. 2016, arXiv:1505.06392) code, and driven through the parameter phase space by a differential-evolution Markov-chain Monte Carlo (MC3, https://github.com/pcubillos/mc3 , Cubillos et al. 2017, arXiv:1610.01336) sampler. In this paper we give a brief description of the framework, and its modules that can be used separately for other scientific purposes; outline the retrieval analysis flow; present the initialization routines, describing in detail the atmospheric profile generator and the temperature and species parameterizations; and specify the post-processing routines and outputs, concentrating on the spectrum band integrator, the best-fit model selection, and the contribution functions. We also present an atmospheric analysis of WASP-43b secondary eclipse data obtained from space- and ground-based observations. We compare our results with the results from the literature, and investigate how the inclusion of additional opacity sources influence the best-fit model.


title: Alpha-decay spectroscopy of $^{257}$Rf

abstract: The decay properties of the ground state and excited states of $^{257}$Rf have been investigated with the detector array GABRIELA at the FLNR, Dubna. The electromagnetic decay of a new excited state in $^{253}$No has been observed. The state lies 750 keV above the ground state and is favourably populated in the alpha decay of the low-lying spin isomer of $^{257}$Rf. It decays to the 9/2$^-$ ground state by an M1 transition and is assigned the 11/2$^-$[725] Nilsson configuration. The presence of this state suggests a possible reinterpretation of the decay of the high-K isomer in $^{253}$No.


title: Quantifying Unknown Entanglement by Neural Networks

abstract: Quantum entanglement plays a crucial role in quantum information processing tasks and quantum mechanics, hence quantifying unknown entanglement is a fundamental task. However, this is also challenging, as entanglement cannot be measured by any observables directly. In this paper, we train neural networks to quantify unknown entanglement, where the input features of neural networks are the outcome statistics data produced by locally measuring target quantum states, and the training labels are well-chosen quantities. For bipartite quantum states, this quantity is coherent information, which is a lower bound for the entanglement of formation and the entanglement of distillation. For multipartite quantum states, we choose this quantity as the geometric measure of entanglement. It turns out that the neural networks we train have very good performance in quantifying unknown quantum states, and can beat previous approaches like semi-device-independent protocols for this problem easily in both precision and application range. We also observe a surprising phenomenon that on quantum states with stronger quantum nonlocality, the neural networks tend to have better performance, though we do not provide them any knowledge on quantum nonlocality.


title: Spatio-Temporal Pruning and Quantization for Low-latency Spiking Neural Networks

abstract: Spiking Neural Networks (SNNs) are a promising alternative to traditional deep learning methods since they perform event-driven information processing. However, a major drawback of SNNs is high inference latency. The efficiency of SNNs could be enhanced using compression methods such as pruning and quantization. Notably, SNNs, unlike their non-spiking counterparts, consist of a temporal dimension, the compression of which can lead to latency reduction. In this paper, we propose spatial and temporal pruning of SNNs. First, structured spatial pruning is performed by determining the layer-wise significant dimensions using principal component analysis of the average accumulated membrane potential of the neurons. This step leads to 10-14X model compression. Additionally, it enables inference with lower latency and decreases the spike count per inference. To further reduce latency, temporal pruning is performed by gradually reducing the timesteps while training. The networks are trained using surrogate gradient descent based backpropagation and we validate the results on CIFAR10 and CIFAR100, using VGG architectures. The spatiotemporally pruned SNNs achieve 89.04% and 66.4% accuracy on CIFAR10 and CIFAR100, respectively, while performing inference with 3-30X reduced latency compared to state-of-the-art SNNs. Moreover, they require 8-14X lesser compute energy compared to their unpruned standard deep learning counterparts. The energy numbers are obtained by multiplying the number of operations with energy per operation. These SNNs also provide 1-4% higher robustness against Gaussian noise corrupted inputs. Furthermore, we perform weight quantization and find that performance remains reasonably stable up to 5-bit quantization.


title: Generation of Local Unitary Groups

abstract: Let $E$ be a two-dimensional \'etale algebra over a non-Archimedean local field $K$ of characteristic zero. We show that the unitary group of a non-degenerate hermitian lattice over $E$ is generated by symmetries and rescaled Eichler isometries. In the appendix we show that unless $E/K$ is a ramified dyadic field extension and the residue field has two elements, symmetries suffice.


title: A class of new stable, explicit methods to solve the non-stationary heat equation

abstract: We present a class of new explicit and stable numerical algorithms to solve the spatially discretized linear heat or diffusion equation. After discretizing the space and the time variables like conventional finite difference methods, we do not approximate the time derivatives by finite differences, but use constant neighbor and linear neighbour approximations to decouple the ordinary differential equations and solve them analytically. During this process, the timestep-size appears not in polynomial, but in exponential form with negative exponents, which guarantees stability. We compare the performance of the new methods with analytical and numerical solutions. According to our results, the methods are first and second order in time and can be much faster than the commonly used explicit or implicit methods, especially in the case of extremely large stiff systems.

categories: math.NA cs.NA math-ph math.MP physics.comp-ph


title: A development of a compact software package for systematic error studies in muon $g-2$/EDM experiment at J-PARC

abstract: The J-PARC muon g-2/EDM experiment aims to measure the muon magnetic moment anomaly $a_{\mu} = (g -2 ) / 2 ~ $. The target sensitivity for $a_{\mu}$ is a statistical uncertainty of $450 \times 10^{-9}$, and $a_{\mu}$ can be extracted by measuring the spin precession of muons ($\omega_a$). The $\omega_a$ can be extracted from the distribution of muon decay time, and the time is measured by reconstructing tracks of positrons which are decayed from the muons. To extract $\omega_a$ precisely, the systematic effects on $\omega_a$ should be controlled under the aforementioned sensitivity. To study systematic effects, we develop a compact simulation package, and using our developed package, we study several systematic effects and discuss results. We discuss the package details and studies of systematic effects on $\omega_a$ based on our package.


title: Confined massive circumstellar shell in type IIL SN 2008fq

abstract: I explore a phenomenon of the circumstellar H alpha and H beta absorption lines in the spectrum of luminous type IIL SN 2008fq taken on day 6.8 after the discovery. The absorpion is identified with the radiatively accelerated preshock wind up to ~800 km/s. The required initial luminosity is attributed to the earlier circumstellar interaction with the confined dense shell of ~0.08 Msun. The modelling of similar luminous type IIL SN 1998S based on the same approach results in the comparable shell mass of ~0.1 Msun. More than 1 dex larger mass of the confined dense shell of both SNe IIL compared to that of type IIP supernovae is attributed to the larger progenitor mass of type IIL supernovae.


title: Visformer: The Vision-friendly Transformer

abstract: The past year has witnessed the rapid development of applying the Transformer module to vision problems. While some researchers have demonstrated that Transformer-based models enjoy a favorable ability of fitting data, there are still growing number of evidences showing that these models suffer over-fitting especially when the training data is limited. This paper offers an empirical study by performing step-by-step operations to gradually transit a Transformer-based model to a convolution-based model. The results we obtain during the transition process deliver useful messages for improving visual recognition. Based on these observations, we propose a new architecture named Visformer, which is abbreviated from the `Vision-friendly Transformer'. With the same computational complexity, Visformer outperforms both the Transformer-based and convolution-based models in terms of ImageNet classification accuracy, and the advantage becomes more significant when the model complexity is lower or the training set is smaller. The code is available at https://github.com/danczs/Visformer.


title: Hydrodynamic Attractor in a Hubble Expansion

abstract: We analytically investigate hydrodynamic attractor solutions in both M\"{u}ller-Israel-Stewart (MIS) and kinetic theories in a viscous fluid system undergoing a Hubble expansion with a fixed expansion rate. We show that the gradient expansion for the MIS theory and the Chapman-Enskog expansion for the Boltzmann equation within the relaxation time approximation are factorially divergent and obtain hydrodynamic attractor solutions by applying the Borel resummation technique to those asymptotic divergent series. In both theories, we find that the hydrodynamic attractor solutions are globally attractive and only the first-order non-hydrodynamic mode exists. We also find that the hydrodynamic attractor solutions in the two theories disagree with each other when gradients become large, and that the speed of the attraction is different. Similarities and differences from hydrodynamic attractors in the Bjorken and Gubser flows are also discussed. Our results push the idea of far-from-equilibrium hydrodynamics in systems undergoing a Hubble expansion.


title: Representations of Extended Carroll Group

abstract: Carroll's group is presented as a group of transformations in a 5-dimensional space ($\mathcal{C}$) obtained by embedding the Euclidean space into a (4; 1)-de Sitter space. Three of the five dimensions of $\mathcal{C}$ are related to $\mathcal{R}^3$, and the other two to mass and time. A covariant formulation of Caroll's group, analogous as introduced by Takahashi to Galilei's group, is deduced. Unit representations are studied.


title: CSI-free Rotary Antenna Beamforming for Massive RF Wireless Energy Transfer

abstract: Radio frequency (RF) wireless energy transfer (WET) is a key technology that may allow seamlessly powering future massive low-energy Internet of Things (IoT) networks. To enable efficient massive WET, channel state information (CSI)-limited/free multi-antenna transmit schemes have been recently proposed in the literature. The idea is to reduce/null the energy costs to be paid by energy harvesting (EH) IoT nodes from participating in large-scale time/power-consuming CSI training, but still enable some transmit spatial gains. In this paper, we take another step forward by proposing a novel CSI-free rotary antenna beamforming (RAB) WET scheme that outperforms all state-of-the-art CSI-free schemes in a scenario where a power beacon (PB) equipped with a uniform linear array (ULA) powers a large set of surrounding EH IoT devices. RAB uses a properly designed CSI-free beamformer combined with a continuous or periodic rotation of the ULA at the PB to provide average EH gains that scale as $0.85\sqrt{M}$, where $M$ is the number of PB's antenna elements. Moreover, a rotation-specific power control mechanism was proposed to i) fairly optimize the WET process if devices' positioning information is available, and/or ii) to avoid hazards to human health in terms of specific absorption rate (SAR). We show that RAB performance even approaches quickly (or surpasses, for scenarios with sufficiently large number of EH devices, or when using the proposed power control) the performance of a traditional full-CSI based transmit scheme, and it is also less sensitive to SAR constraints. Finally, we discuss important practicalities related to RAB such as its robustness against non line-of-sight conditions compared to other CSI-free WET schemes, and its generalizability to scenarios where the PB uses other than a ULA topology.


title: Computer vision in automated parking systems: Design, implementation and challenges

abstract: Automated driving is an active area of research in both industry and academia. Automated Parking, which is automated driving in a restricted scenario of parking with low speed manoeuvring, is a key enabling product for fully autonomous driving systems. It is also an important milestone from the perspective of a higher end system built from the previous generation driver assistance systems comprising of collision warning, pedestrian detection, etc. In this paper, we discuss the design and implementation of an automated parking system from the perspective of computer vision algorithms. Designing a low-cost system with functional safety is challenging and leads to a large gap between the prototype and the end product, in order to handle all the corner cases. We demonstrate how camera systems are crucial for addressing a range of automated parking use cases and also, to add robustness to systems based on active distance measuring sensors, such as ultrasonics and radar. The key vision modules which realize the parking use cases are 3D reconstruction, parking slot marking recognition, freespace and vehicle/pedestrian detection. We detail the important parking use cases and demonstrate how to combine the vision modules to form a robust parking system. To the best of the authors' knowledge, this is the first detailed discussion of a systemic view of a commercial automated parking system.


title: Exponential Competence of Computer Science and Software Engineering Undergraduate Students

abstract: We live in exceptional times in which the entire world is witnessing the exponential spread of a pandemic, which requires to adopt new habits of mind and behaviors. In this paper, I introduce the term exponential competence, which encompasses these cognitive and social skills, and describe a course for computer science and software engineering students in which emphasis is placed on exponential competence. I argue that exponential competence is especially important for computer science and software engineering students, since many of them will, most likely, be required to deal with exponential phenomena in their future professional development.


title: Exploratory Learning Environments for Responsible Management Education Using Lego Serious Play

abstract: Research into responsible management education has largely focused on the merits, attributes, and transformation opportunities to enhance responsible business school education aims. As such, a prominent part of the literature has occupied itself with examining if responsible management modules are inherently considered a non-crucial element of the curriculum and determining the extent to which business schools have introduced such learning content into their curriculum. However, there has been scant research into how to apply novel teaching approaches to engage students and promote responsible management education endeavours. As such, this paper seeks to address this gap through the development of a teaching framework to support educators in designing effective learning environments focused on responsible management education. We will draw on constructivist learning theories and Lego Serious Play (LSP) as a learning enhancement approach to develop a pedagogical framework. LSP is selected due to its increasing application in learning environments to help promote critical discourse, and engage with highly complex problems, whether these are social, economic, environmental, or organisational.


title: Detecting new edge types in a temporal network model

abstract: Networks representing complex systems in nature and society usually involve multiple interaction types. These types suggest essential information on the interactions between components, but not all of the existing types are usually discovered. Therefore, detecting the undiscovered edge types is crucial for deepening our understanding of the network structure. Although previous studies have discussed the edge label detection problem, we still lack effective methods for uncovering previously-undetected edge types. Here, we develop an effective technique to detect undiscovered new edge types in networks by leveraging a novel temporal network model. Both analytical and numerical results show that the prediction accuracy of our method is perfect when the model networks' time parameter approaches infinity. Furthermore, we find that when time is finite, our method is still significantly more accurate than the baseline.


title: Feasibility Study of Aerocapture at Mars with an Innovative Deployable Heat Shield

abstract: Performing orbital insertion around Mars using aerocapture instead of a propulsive orbit insertion manoeuvre allows to save resources and/or increase the payload mass fraction. Aerocapture has never been employed to date because of the high uncertainties in the parameters from which it depends, mainly related to atmospheric density modeling and navigation errors. The purpose of this work is to investigate the feasibility of aerocapture at Mars with an innovative deployable drag device, whose aperture can be modulated in flight, and assess the effects of the main uncertainties on the success of the manoeuvre. This paper starts with the presentation of a parametric bi-dimensional analysis of the effectiveness of aerocapture, for which a wide range of uncertainty levels in the atmospheric density and the ballistic coefficient are considered. Then, an application to a real mission scenario is carried out including the error of the targeting manoeuvre performed at the limit of the sphere of influence of the planet. The analyses show the strong influence of the uncertainties in the atmospheric density and the ballistic coefficient, which significantly narrow the solution space and limit its continuity. However, viable solutions for aerocapture can still be identified even in the worst conditions.


title: Catalogs of C and Python Antipatterns by CS1 Students

abstract: Understanding students' programming misconceptions is critical. Doing so depends on identifying the reasons why students make errors when learning a new programming language. Knowing the misconceptions can help students to improve their reflection about their mistakes and also help instructors to design better teaching strategies. In this technical report, we propose catalogs of antipatterns for two programming languages: C and Python. To accomplish this, we analyzed the codes of 166 CS1 engineering students when they were coding solutions to programming exercises. In our results, we catalog 41 CS1 antipatterns from 95 cataloged misconceptions in C and Python. These antipatterns were separated into three catalogs: C, Python, and antipatterns found in code using both programming languages. For each antipattern, we present code examples, students' solutions (if they are present), a possible solution to avoid the antipattern, among other information.


title: Enhanced dynamic homogenization of hexagonally packed granular materials with elastic interfaces

abstract: It is well known that the classical energetically consistent micropolar model has limits in simulating the frequency band structure of packed granular materials (see Merkel et al., 2011). It is here shown that if a standard continualization of the difference equation of motion of the discrete model is carried out, an equivalent micropolar continuum is obtained which is able to accurately simulate the optical branches of the discrete model. Nevertheless, this homogenized continuum presents non-positive defined elastic potential energy, a deficiency that limits the reliability of the model and implies instability phenomena (destabilizing effects) in the acoustic branches. This drawback is circumvented here through an high-frequency dynamic homogenization scheme which is based on an enhanced continualization of the discrete governing equations into pseudo-differential equations. Through a formal Taylor expansion of the pseudo-differential operators a higher order differential equation corresponding to the governing equation of a non-local continuum thermodynamically consistent are obtained. The resulting approach allows obtaining an equivalent micropolar continuum characterized by inertial non-locality. Moreover, higher order continua with non-local constitutive and inertial terms may be derived. The proposed continuum models are proved to be able to accurately describe both the static and dynamic behavior of the discrete granular model. Finally, the convergence to the response of the discrete system is shown when increasing the order of the higher order continuum.


title: The Myth of Complete AI-Fairness

abstract: The idea of fairness and justice has long and deep roots in Western civilization, and is strongly linked to ethics. It is therefore not strange that it is core to the current discussion about the ethics of development and use of AI systems. In this short paper, I wish to further motivate my position in this matter: ``I will never be completely fair. Nothing ever is. The point is not complete fairness, but the need to establish metrics and thresholds for fairness that ensure trust in AI systems".


title: Agile (data) science: a (draft) manifesto

abstract: Science has a data management problem, as well as a project management problem. While industrial-grade data science teams have embraced the agile mindset, and adopted or created all kind of tools to create reproducible workflows, academia-based science is still (mostly) mired in a mindset that is focused on a single final product (a paper), without focusing on incremental improvement, on any specific problem or customer, or, paying any attention reproducibility. In this report we argue towards the adoption of the agile mindset and agile data science tools in academia, to make a more responsible, and over all, reproducible science.


title: The Effects of Air Quality on the Spread of the COVID-19 Pandemic in Italy: An Artificial Intelligence Approach

abstract: The COVID-19 pandemic considerably affects public health systems around the world. The lack of knowledge about the virus, the extension of this phenomenon, and the speed of the evolution of the infection are all factors that highlight the necessity of employing new approaches to study these events. Artificial intelligence techniques may be useful in analyzing data related to areas affected by the virus. The aim of this work is to investigate any possible relationships between air quality and confirmed cases of COVID-19 in Italian districts. Specifically, we report an analysis of the correlation between daily COVID-19 cases and environmental factors, such as temperature, relative humidity, and atmospheric pollutants. Our analysis confirms a significant association of some environmental parameters with the spread of the virus. This suggests that machine learning models trained on the environmental parameters to predict the number of future infected cases may be accurate. Predictive models may be useful for helping institutions in making decisions for protecting the population and contrasting the pandemic.


title: A Framework for Ethical AI at the United Nations

abstract: This paper aims to provide an overview of the ethical concerns in artificial intelligence (AI) and the framework that is needed to mitigate those risks, and to suggest a practical path to ensure the development and use of AI at the United Nations (UN) aligns with our ethical values. The overview discusses how AI is an increasingly powerful tool with potential for good, albeit one with a high risk of negative side-effects that go against fundamental human rights and UN values. It explains the need for ethical principles for AI aligned with principles for data governance, as data and AI are tightly interwoven. It explores different ethical frameworks that exist and tools such as assessment lists. It recommends that the UN develop a framework consisting of ethical principles, architectural standards, assessment methods, tools and methodologies, and a policy to govern the implementation and adherence to this framework, accompanied by an education program for staff.


title: The Cardan grille approach to the Voynich MS taken to the next level

abstract: The Voynich MS is an illustrated 15th century manuscript, whose text is written in an unknown alphabet, which has not been translated until today. In 2004 Gordon Rugg published a paper in which he proposed that this text is likely to be meaningless, and could have been composed by an alternative application of a so-called Cardan Grille, namely by moving a piece of cardboard with holes over a large table of word fragments, and writing down the words that thus appear. This paper caused considerable discussion in the circles of people interested in the Voynich MS text, but it has not found many followers, even until today. The present paper takes a closer look at the mechanics of this method. Based on this, a more generic method is proposed, which is considerably simpler, both to set up and to execute. It is shown that the unusual word length distribution of the Voynich MS, which is very close to binomial, could be a consequence of the application of such a method. Furthermore, it is argued that this method could not only be used to create meaningless text, but also to encode meaningful text. A first high-level analysis looks at whether such a method could indeed have been applied to create the Voynich MS text, but this is certainly far from conclusive. The main aim of this paper is to inspire further research of the Voynich MS text into a new direction that has not yet been explored in great detail.


title: Bulk High-Entropy Hexaborides

abstract: For the first time, a group of CaB6-typed cubic rare earth high-entropy hexaborides have been successfully fabricated into dense bulk pellets (>98.5% in relative densities). The specimens are prepared from elemental precursors via in-situ metal-boron reactive spark plasma sintering. The sintered bulk pellets are determined to be single-phase without any detectable oxides or other secondary phases. The homogenous elemental distributions have been confirmed at both microscale and nanoscale. The Vickers microhardness are measured to be 16-18 GPa at a standard indentation load of 9.8 N. The nanoindentation hardness and Young's moduli have been measured to be 19-22 GPa and 190-250 GPa, respectively, by nanoindentation test using a maximum load of 500 mN. The material work functions are determined to be 3.7-4.0 eV by ultraviolet photoelectron spectroscopy characterizations, which are significantly higher than that of LaB6.


title: Short-Range Order and Origin of the Low Thermal Conductivity in Compositionally Complex Rare-Earth Niobates and Tantalates

abstract: Rare-earth niobates and tantalates possess low thermal conductivities, which can be further reduced in high-entropy compositions. Here, a large number of 40 compositions are synthesized to investigate the origin of low thermal conductivity. Amongst, 29 possess single (nominally cubic) fluorite phases and most of them are new compositionally complex (medium- or high-entropy) compositions. One new finding is that doping 2 % of light element cations can further reduce thermal conductivity. This large data set enables the discovery of a negative correlation between the thermal conductivity and averaged radius ratio of the 3+/5+ cations. While this ratio is still below the threshold for forming long-range ordered weberite phases, this correlation suggests the reduced thermal conductivity is related to short-range weberite order, which is indeed revealed by diffuse scattering in X-ray and neutron diffraction. Specifically, neutron diffraction is used characterize five selected specimens. A better fit to a weberite structure is found at nanoscale (~1 nm). The characteristic length (domain size) is smaller but with stronger short-range ordering in more insulative materials. As it approaches the Ioffe-Regel limit, the phonon limit breaks down and "diffusons" give rise to the observed amorphous-like thermal conductivity. Disordered oxygen sublattices are also confirmed by neutron diffraction.


title: Zinbiel 2-algebras

abstract: The notions of Zinbiel 2-algebras and 2-term $Z_{\infty}$-algebras are introduced. It is proved that the category of Zinbiel 2-algebras and the category of $2$-term $Z_{\infty}$-algebras are equivalent to each other. Crossed module extensions of Zinbiel algebras are studied. Relationship between Zinbiel 2-algebras and dendriform 2-algebras are found.


title: Generalized heterogeneous hypergeometric functions and the distribution of the largest eigenvalue of an elliptical Wishart matrix

abstract: In this study, we derive the exact distributions of eigenvalues of a singular Wishart matrix under an elliptical model. We define generalized heterogeneous hypergeometric functions with two matrix arguments and provide convergence conditions for these functions. The joint density of eigenvalues and the distribution function of the largest eigenvalue for a singular elliptical Wishart matrix are represented by these functions. Numerical computations for the distribution of the largest eigenvalue were conducted under the matrix-variate $t$ and Kotz-type models.


title: Avoiding bias when inferring race using name-based approaches

abstract: Racial disparity in academia is a widely acknowledged problem. The quantitative understanding of racial based systemic inequalities is an important step towards a more equitable research system. However, because of the lack of robust information on authors' race, few large scale analyses have been performed on this topic. Algorithmic approaches offer one solution, using known information about authors, such as their names, to infer their perceived race. As with any other algorithm, the process of racial inference can generate biases if it is not carefully considered. The goal of this article is to assess the extent to which algorithmic bias is introduced using different approaches for name based racial inference. We use information from the U.S. Census and mortgage applications to infer the race of U.S. affiliated authors in the Web of Science. We estimate the effects of using given and family names, thresholds or continuous distributions, and imputation. Our results demonstrate that the validity of name based inference varies by race/ethnicity and that threshold approaches underestimate Black authors and overestimate White authors. We conclude with recommendations to avoid potential biases. This article lays the foundation for more systematic and less biased investigations into racial disparities in science.


title: A stability preserved time-integration method for nonlinear advection-diffusion models

abstract: A new implicit-explicit local differential transform method (IELDTM) is derived here for time integration of the nonlinear advection-diffusion processes represented by (2+1)-dimensional Burgers equation. The IELDTM is adaptively constructed as stability preserved and high order time integrator for spatially discretized Burgers equation. For spatial discretization of the model equation, the Chebyshev spectral collocation method (ChCM) is utilized. A robust stability analysis and global error analysis of the IELDTM are presented with respect to the direction parameter \theta. With the help of the global error analysis, adaptivity equations are derived to minimize the computational costs of the algorithms. The produced method is shown to eliminate the accuracy disadvantage of the classical \theta-method and the stability disadvantages of the existing DTM-based methods. Two examples of the Burgers equation in one and two dimensions have been solved via the ChCM-IELDTM hybridization, and the produced results are compared with the literature. The present time integrator has been proven to produce more accurate numerical results than the MATLAB solvers, ode45 and ode15s.


title: Linking open-source code commits and MOOC grades to evaluate massive online open peer review

abstract: Massive Open Online Courses (MOOCs) have been used by students as a low-cost and low-touch educational credential in a variety of fields. Understanding the grading mechanisms behind these course assignments is important for evaluating MOOC credentials. A common approach to grading free-response assignments is massive scale peer-review, especially used for assignments that are not easy to grade programmatically. It is difficult to assess these approaches since the responses typically require human evaluation. Here we link data from public code repositories on GitHub and course grades for a large massive-online open course to study the dynamics of massive scale peer review. This has important implications for understanding the dynamics of difficult to grade assignments. Since the research was not hypothesis-driven, we described the results in an exploratory framework. We find three distinct clusters of repeated peer-review submissions and use these clusters to study how grades change in response to changes in code submissions. Our exploration also leads to an important observation that massive scale peer-review scores are highly variable, increase, on average, with repeated submissions, and changes in scores are not closely tied to the code changes that form the basis for the re-submissions.


title: COVID-19 Modeling: A Review

abstract: The SARS-CoV-2 virus and COVID-19 disease have posed unprecedented and overwhelming demand, challenges and opportunities to domain, model and data driven modeling. This paper provides a comprehensive review of the challenges, tasks, methods, progress, gaps and opportunities in relation to modeling COVID-19 problems, data and objectives. It constructs a research landscape of COVID-19 modeling tasks and methods, and further categorizes, summarizes, compares and discusses the related methods and progress of modeling COVID-19 epidemic transmission processes and dynamics, case identification and tracing, infection diagnosis and medical treatments, non-pharmaceutical interventions and their effects, drug and vaccine development, psychological, economic and social influence and impact, and misinformation, etc. The modeling methods involve mathematical and statistical models, domain-driven modeling by epidemiological compartmental models, medical and biomedical analysis, AI and data science in particular shallow and deep machine learning, simulation modeling, social science methods, and hybrid modeling.


title: Flexible Educational Software Architecture

abstract: EAs.LiT is an e-assessment management and analysis software for which contextual requirements and usage scenarios changed over time. Based on these factors and further development activities, the decision was made to adopt a microservice architecture for EAs.LiT version 2 in order to increase its flexibility to adapt to new and changed circumstances. This architectural style and a few adopted technologies, like RDF as a data format, enabled an eased implementation of various use cases. Thus we consider the microservice architecture productive and recommend it for usage in other educational projects. The specific architecture of EAs.LiT version 2 is presented within this article, targeting to enable other educational projects to adopt it by using our work as a foundation or template.


title: EduPal leaves no professor behind: Supporting faculty via a peer-powered recommender system

abstract: The swift transitions in higher education after the COVID-19 outbreak identified a gap in the pedagogical support available to faculty. We propose a smart, knowledge-based chatbot that addresses issues of knowledge distillation and provides faculty with personalized recommendations. Our collaborative system crowdsources useful pedagogical practices and continuously filters recommendations based on theory and user feedback, thus enhancing the experiences of subsequent peers. We build a prototype for our local STEM faculty as a proof concept and receive favorable feedback that encourages us to extend our development and outreach, especially to underresourced faculty.


title: Whether the Health Care Practices For the Patients With Comorbidities Have Changed After the Outbreak of COVID-19; Big Data Public Sentiment Analysis

abstract: After the pandemic of SARS-CoV-2, it has influenced the health care practices around the world. Initial investigations indicate that patients with comorbidities are more fragile to this SARS-CoV-2 infection. They suggested postponing the routine treatment of cancer patients. However, few meta-analyses suggested evidences are not sufficient to hold the claim of the frailty of cancer patients to COVID-19, and they are not in favour of shelving the scheduled procedures. There are recent studies in which medical professionals, according to their competence, are referring to change the routine practices on how to manage the applicable therapeutic resources judiciously to combat this vital infection. This is a different study that reveals the cancer patients' viewpoint about how health care practices have been changed in their opinion during this pandemic year? Are they satisfied with their treatment or not? To serve the purpose, we gathered more than 60000 relevant tweets from Twitter to analyse the sentiment of cancer patients around the world. Our findings demonstrate that there is a surge in argument about cancer and its treatment after the outbreak of COVID-19. Most of the tweets are reasonable (52.6%) compared to the negative ones (24.3). We developed polarity and subjectivity distribution to better recognise the positivity/negativity in the sentiment. Results reveal that the polarity range of positive tweets is within the range of 0 to 0.5. Which means the tendency in the tweets is not so much positive but surely not negative. It is a piece of modest statistical evidence in support of how natural language processing (NLP) can be accepted to better understand the patient's behaviour in real-time, and it may facilitate the medical professional to make better decision to organise the routine management of cancer patients.


title: Computing Grain Boundary Diagrams of Thermodynamic and Mechanical Properties

abstract: Computing the grain boundary (GB) counterparts to bulk phase diagrams represents an emerging research direction. Using a classical embrittlement model system Ga-doped Al alloy, this study demonstrates the feasibility of computing temperature- and composition-dependent GB diagrams to represent not only equilibrium thermodynamic and structural characters, but also mechanical properties. Specifically, hybrid Monte Carlo and molecular dynamics (MC/MD) simulations are used to obtain the equilibrium GB structure as a function of temperature and composition. Simulated GB structures are validated by aberration-corrected scanning transmission electron microscopy. Subsequently, MD tensile tests are performed on the simulated equilibrium GB structures. GB diagrams are computed for not only GB adsorption and structural disorder, but also interfacial structural and chemical widths, MD ultimate tensile strength, and MD tensile toughness. This study suggests a research direction to investigate GB composition-structure-property relationships via computing GB diagrams of thermodynamic, structural, and mechanical (or potentially other) properties.


title: Artificial Intelligence Based Prognostic Maintenance of Renewable Energy Systems: A Review of Techniques, Challenges, and Future Research Directions

abstract: Since the depletion of fossil fuels, the world has started to rely heavily on renewable sources of energy. With every passing year, our dependency on the renewable sources of energy is increasing exponentially. As a result, complex and hybrid generation systems are being designed and developed to meet the energy demands and ensure energy security in a country. The continual improvement in the technology and an effort towards the provision of uninterrupted power to the end-users is strongly dependent on an effective and fault resilient Operation and Maintenance (O&M) system. Ingenious algorithms and techniques are hence been introduced aiming to minimize equipment and plant downtime. Efforts are being made to develop robust Prognostic Maintenance systems that can identify the faults before they occur. To this aim, complex Data Analytics and Machine Learning (ML) techniques are being used to increase the overall efficiency of these prognostic maintenance systems. This paper provides an overview of the predictive/prognostic maintenance frameworks reported in the literature. We pay a particular focus to the approaches, challenges including data-related issues, such as the availability and quality of the data and data auditing, feature engineering, interpretability, and security issues. Being a key aspect of ML-based solutions, we also discuss some of the commonly used publicly available datasets in the domain. The paper also identifies key future research directions. We believe such detailed analysis will provide a baseline for future research in the domain.


title: On the p-biharmonic submanifolds and stress p-bienergy tensors

abstract: In this paper, we consider p-biharmonic submanifolds of aspace form. We give the necessary and sufficient conditions for a submanifold to be p-biharmonic in a space form. We present some new properties for the stress p-bienergy tensor.


title: Fermi-Surface Modeling of Light-Rare-Earth Hexaborides with 2D-ACAR Spectroscopy

abstract: Two dimensional angular correlation of the positron annihilation radiation (2D-ACAR) spectra are measured for $\mathrm{LaB}_6$ along high symmetry directions and compared with first principle calculations based on density functional theory (DFT). This allows the modeling of the Fermi surface in terms of ellipsoid electron pockets centered at $X$-points elongated along the $\Sigma$ axis (${\Gamma-M}$ direction). The obtained structure is in agreement with quantum oscillation measurements and previous band structure calculations. For the isostructural topologically not-trivial $\mathrm{SmB}_6$ the similar ellipsoids are connected through necks that have significantly smaller radii in the case of $\mathrm{LaB}_6$. A theoretical analysis of the 2D-ACAR spectra is also performed for $\mathrm{CeB}_6$ including the on-site repulsion $U$ correction to the local-density approximation (LDA+$U$) of the DFT. The similarities of 2D-ACAR spectra and the Fermi-surface projections of these two compounds allow to infer that both $\mathrm{LaB}_6$ and $\mathrm{CeB}_6$ are topologically trivial correlated metals.


title: Synthetic 3D Data Generation Pipeline for Geometric Deep Learning in Architecture

abstract: With the growing interest in deep learning algorithms and computational design in the architectural field, the need for large, accessible and diverse architectural datasets increases. We decided to tackle this problem by constructing a field-specific synthetic data generation pipeline that generates an arbitrary amount of 3D data along with the associated 2D and 3D annotations. The variety of annotations, the flexibility to customize the generated building and dataset parameters make this framework suitable for multiple deep learning tasks, including geometric deep learning that requires direct 3D supervision. Creating our building data generation pipeline we leveraged architectural knowledge from experts in order to construct a framework that would be modular, extendable and would provide a sufficient amount of class-balanced data samples. Moreover, we purposefully involve the researcher in the dataset customization allowing the introduction of additional building components, material textures, building classes, number and type of annotations as well as the number of views per 3D model sample. In this way, the framework would satisfy different research requirements and would be adaptable to a large variety of tasks. All code and data are made publicly available.


title: Mutual Contrastive Learning for Visual Representation Learning

abstract: We present a collaborative learning method called Mutual Contrastive Learning (MCL) for general visual representation learning. The core idea of MCL is to perform mutual interaction and transfer of contrastive distributions among a cohort of networks. A crucial component of MCL is Interactive Contrastive Learning (ICL). Compared with vanilla contrastive learning, ICL can aggregate cross-network embedding information and maximize the lower bound to the mutual information between two networks. This enables each network to learn extra contrastive knowledge from others, leading to better feature representations for visual recognition tasks. We emphasize that the resulting MCL is conceptually simple yet empirically powerful. It is a generic framework that can be applied to both supervised and self-supervised representation learning. Experimental results on image classification and transfer learning to object detection show that MCL can lead to consistent performance gains, demonstrating that MCL can guide the network to generate better feature representations. Code is available at https://github.com/winycg/MCL.


title: Plectic p-adic invariants

abstract: For modular elliptic curves over number fields of narrow class number one, and with multiplicative reduction at a collection of p-adic primes, we define new p-adic invariants. Inspired by Nekovar and Scholl's plectic conjectures, we believe these invariants control the Mordell--Weil group of higher rank elliptic curves and we support our expectations with numerical experiments.


title: Evaluating the Values of Sources in Transfer Learning

abstract: Transfer learning that adapts a model trained on data-rich sources to low-resource targets has been widely applied in natural language processing (NLP). However, when training a transfer model over multiple sources, not every source is equally useful for the target. To better transfer a model, it is essential to understand the values of the sources. In this paper, we develop SEAL-Shap, an efficient source valuation framework for quantifying the usefulness of the sources (e.g., domains/languages) in transfer learning based on the Shapley value method. Experiments and comprehensive analyses on both cross-domain and cross-lingual transfers demonstrate that our framework is not only effective in choosing useful transfer sources but also the source values match the intuitive source-target similarity.


title: An action for nonlinear dislocation dynamics

abstract: An action functional is developed for nonlinear dislocation dynamics. This serves as a first step towards the application of effective field theory in physics to evaluate its potential in obtaining a macroscopic description of dislocation dynamics describing the plasticity of crystalline solids. Connections arise between the continuum mechanics and material science of defects in solids, effective field theory techniques in physics, and fracton tensor gauge theories. The scheme that emerges from this work for generating a variational principle for a nonlinear pde system is general, as is demonstrated by doing so for nonlinear elastostatics involving a stress response function that is not necessarily hyperelastic.


title: Sommerfeld enhancement factor in two-dimensional Dirac materials

abstract: In this work the above-band gap absorption spectrum in two-dimensional Dirac materials is calculated with account for the interaction between the photocarriers. Both the screened Rytova-Keldysh and pure Coulomb attraction potentials between the electron and hole are used in the study. We find that, in the materials under consideration, the interaction enhances the absorbance in the narrow interband edge region, in a sharp contrast to the band model with the parabolic free-carrier energy dispersion. We develop an approximation of the weak interaction which allows us to reproduce the main features of the exactly calculated Sommerfeld factor. We show a substantial reduction of this factor at higher frequencies due to the single-particle energy renormalization.


title: Algorithmic Solution for Systems of Linear Equations, in $\mathcal{O}(mn)$ time

abstract: We present a novel algorithm attaining excessively fast, the sought solution of linear systems of equations. The algorithm is short in its basic formulation and, by definition, vectorized, while the memory allocation demands are trivial, because, for each iteration, only one dimension of the given input matrix $\mathbf X$ is utilized. The execution time is very short compared with state-of-the-art methods, exhibiting $> \times 10^2$ speed-up and low memory allocation demands, especially for non-square Systems of Linear Equations, with ratio of equations versus features high (tall systems), or low (wide systems) accordingly. The accuracy is high and straightforwardly controlled, and the numerical results highlight the efficiency of the proposed algorithm, in terms of computation time, solution accuracy and memory demands. The paper also comprises a theoretical proof for the algorithmic convergence, and we extend the implementation of the proposed algorithmic rationale to feature selection tasks.


title: Statistical Modeling of Airborne Virus Transmission Through Imperfectly Fitted Face Masks

abstract: The rapid emergence and the disastrous impact of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) pandemic on public health, societies, and economies around the world has created an urgent need for understanding the pathways critical for virus transmission. Airborne virus transmission by asymptomatic SARS-CoV-2-infected individuals is considered to be a major contributor to the spread of SARS-CoV-2 and social distancing and wearing of face masks in public have been implemented as countermeasures in many countries. However, a comprehensive risk assessment framework for the airborne transmission of SARS-CoV-2 incorporating realistic assumptions on the filtration of infectious aerosols (IAs) by face masks is not available yet. In particular, in most end-to-end models for airborne virus transmission, it is neglected that the stochastic spread of IAs through imperfectly fitted face masks depends on the dynamics of the breathing of the wearer. In this paper, we consider airborne virus transmission from an infected but asymptomatic person to a healthy person, both wearing imperfectly fitted face masks, in an indoor environment. By framing the end-to-end virus transmission as a Molecular Communications (MC) system, we obtain a statistical description of the number of IAs inhaled by the healthy person subject to the respective configurations of the face masks of both persons. We demonstrate that the exhalation and inhalation air flow dynamics have a significant impact on the stochastic filtering of IAs by the face masks. Furthermore, we show that the fit of the face mask of the infected person can highly impact the infection probability. We conclude that the proposed MC model may contribute a valuable assessment tool to fight the spread of SARS-CoV-2 as it encompasses the randomness of the transmission process and enables comprehensive risk analysis beyond statistical averages.


title: Multi-scale PIIFD for Registration of Multi-source Remote Sensing Images

abstract: This paper aims at providing multi-source remote sensing images registered in geometric space for image fusion. Focusing on the characteristics and differences of multi-source remote sensing images, a feature-based registration algorithm is implemented. The key technologies include image scale-space for implementing multi-scale properties, Harris corner detection for keypoints extraction, and partial intensity invariant feature descriptor (PIIFD) for keypoints description. Eventually, a multi-scale Harris-PIIFD image registration algorithm framework is proposed. The experimental results of four sets of representative real data show that the algorithm has excellent, stable performance in multi-source remote sensing image registration, and can achieve accurate spatial alignment, which has strong practical application value and certain generalization ability.


title: Robust decision-making under risk and ambiguity

abstract: Economists often estimate economic models on data and use the point estimates as a stand-in for the truth when studying the model's implications for optimal decision-making. This practice ignores model ambiguity, exposes the decision problem to misspecification, and ultimately leads to post-decision disappointment. Using statistical decision theory, we develop a framework to explore, evaluate, and optimize robust decision rules that explicitly account for estimation uncertainty. We show how to operationalize our analysis by studying robust decisions in a stochastic dynamic investment model in which a decision-maker directly accounts for uncertainty in the model's transition dynamics.


title: Detecting and Matching Related Objects with One Proposal Multiple Predictions

abstract: Tracking players in sports videos is commonly done in a tracking-by-detection framework, first detecting players in each frame, and then performing association over time. While for some sports tracking players is sufficient for game analysis, sports like hockey, tennis and polo may require additional detections, that include the object the player is holding (e.g. racket, stick). The baseline solution for this problem involves detecting these objects as separate classes, and matching them to player detections based on the intersection over union (IoU). This approach, however, leads to poor matching performance in crowded situations, as it does not model the relationship between players and objects. In this paper, we propose a simple yet efficient way to detect and match players and related objects at once without extra cost, by considering an implicit association for prediction of multiple objects through the same proposal box. We evaluate the method on a dataset of broadcast ice hockey videos, and also a new public dataset we introduce called COCO +Torso. On the ice hockey dataset, the proposed method boosts matching performance from 57.1% to 81.4%, while also improving the meanAP of player+stick detections from 68.4% to 88.3%. On the COCO +Torso dataset, we see matching improving from 47.9% to 65.2%. The COCO +Torso dataset, code and pre-trained models will be released at https://github.com/foreverYoungGitHub/detect-and-match-related-objects.


title: Plectic Stark-Heegner points

abstract: We propose a conjectural construction of global points on modular elliptic curves over arbitrary number fields, generalizing both the p-adic construction of Heegner points via Cerednik-Drinfeld uniformization and the definition of classical Stark-Heegner points. In alignment with Nekovar and Scholl's plectic conjectures, we expect the non-triviality of these plectic Stark-Heegner points to control the Mordell-Weil group of higher rank elliptic curves. We provide some indirect evidence for our conjectures by showing that higher order derivatives of anticyclotomic p-adic L-functions compute plectic invariants.


title: A Splicing Approach to Best Subset of Groups Selection

abstract: Best subset of groups selection (BSGS) is the process of selecting a small part of non-overlapping groups to achieve the best interpretability on the response variable. It has attracted increasing attention and has far-reaching applications in practice. However, due to the computational intractability of BSGS in high-dimensional settings, developing efficient algorithms for solving BSGS remains a research hotspot. In this paper,we propose a group-splicing algorithm that iteratively detects the relevant groups and excludes the irrelevant ones. Moreover, coupled with a novel group information criterion, we develop an adaptive algorithm to determine the optimal model size. Under mild conditions, it is certifiable that our algorithm can identify the optimal subset of groups in polynomial time with high probability. Finally, we demonstrate the efficiency and accuracy of our methods by comparing them with several state-of-the-art algorithms on both synthetic and real-world datasets.


title: Optimal controller synthesis for timed systems

abstract: Weighted timed games are zero-sum games played by two players on a timed automaton equipped with weights, where one player wants to minimise the cumulative weight while reaching a target. Used in a reactive synthesis perspective, this quantitative extension of timed games allows one to measure the quality of controllers in real-time systems. Weighted timed games are notoriously difficult and quickly undecidable, even when restricted to non-negative weights. For non-negative weights, the largest class that can be analysed has been introduced by Bouyer, Jaziri and Markey in 2015. Though the value problem is undecidable, the authors show how to approximate the value by considering regions with a refined granularity. In this work, we extend this class to incorporate negative weights, allowing one to model energy for instance, and prove that the value can still be approximated, with the same complexity. A small restriction also allows us to obtain a class of decidable weighted timed games with negative weights and an arbitrary number of clocks. In addition, we show that a symbolic algorithm, relying on the paradigm of value iteration, can be used as an approximation/computation schema over these classes. We also consider the special case of untimed weighted games, where the same fragments are solvable in polynomial time: this contrasts with the pseudo-polynomial complexity, known so far, for weighted games without restrictions.


title: Dissipation enhancement for a degenerated parabolic equation

abstract: In this paper, we quantitatively consider the enhanced-dissipation effect of the advection term to the parabolic $p$-Laplacian equations. More precisely, we show the mixing property of flow for the passive scalar enhances the dissipation process of the $p$-Laplacian in the sense of $L^2$ decay, that is, the $L^2$ decay can be arbitrarily fast. The main ingredient of our argument is to understand the underlying iteration structure inherited from the parabolic $p$-Laplacian equations. This extends the dissipation enhancement result of the advection diffusion equation by Yuanyuan Feng and Gautam Iyer into a non-linear setting.


title: Learning from Event Cameras with Sparse Spiking Convolutional Neural Networks

abstract: Convolutional neural networks (CNNs) are now the de facto solution for computer vision problems thanks to their impressive results and ease of learning. These networks are composed of layers of connected units called artificial neurons, loosely modeling the neurons in a biological brain. However, their implementation on conventional hardware (CPU/GPU) results in high power consumption, making their integration on embedded systems difficult. In a car for example, embedded algorithms have very high constraints in term of energy, latency and accuracy. To design more efficient computer vision algorithms, we propose to follow an end-to-end biologically inspired approach using event cameras and spiking neural networks (SNNs). Event cameras output asynchronous and sparse events, providing an incredibly efficient data source, but processing these events with synchronous and dense algorithms such as CNNs does not yield any significant benefits. To address this limitation, we use spiking neural networks (SNNs), which are more biologically realistic neural networks where units communicate using discrete spikes. Due to the nature of their operations, they are hardware friendly and energy-efficient, but training them still remains a challenge. Our method enables the training of sparse spiking convolutional neural networks directly on event data, using the popular deep learning framework PyTorch. The performances in terms of accuracy, sparsity and training time on the popular DVS128 Gesture Dataset make it possible to use this bio-inspired approach for the future embedding of real-time applications on low-power neuromorphic hardware.


title: Ge thin-films with tantalum diffusion-barriers for use in Nb-based superconductor technology

abstract: Germanium thin films are an excellent candidate for use as a low-loss dielectric in superconducting microwave resonators, a low-loss inter-layer metal wiring dielectric, and passivation layers in microwave and Josephson junction devices. In Ge/Nb structures deposited at 400 {\deg}C, we observe intermixing over as much as 20 nm. The addition of a 10 nm Ta diffusion barrier layer reduces the superconductor/dielectric intermixing to less than 5 nm and enhances the structural properties of deposited a-Ge layers based on Raman spectroscopy. Additionally, superconducting microwave resonators fabricated at room-temperature on crystalline Ge substrates with a Ta barrier layer show marked improvement in total and power-dependent two-level system microwave losses.


title: FedDPGAN: Federated Differentially Private Generative Adversarial Networks Framework for the Detection of COVID-19 Pneumonia

abstract: Existing deep learning technologies generally learn the features of chest X-ray data generated by Generative Adversarial Networks (GAN) to diagnose COVID-19 pneumonia. However, the above methods have a critical challenge: data privacy. GAN will leak the semantic information of the training data which can be used to reconstruct the training samples by attackers, thereby this method will leak the privacy of the patient. Furthermore, for this reason that is the limitation of the training data sample, different hospitals jointly train the model through data sharing, which will also cause the privacy leakage. To solve this problem, we adopt the Federated Learning (FL) frame-work which is a new technique being used to protect the data privacy. Under the FL framework and Differentially Private thinking, we propose a FederatedDifferentially Private Generative Adversarial Network (FedDPGAN) to detectCOVID-19 pneumonia for sustainable smart cities. Specifically, we use DP-GAN to privately generate diverse patient data in which differential privacy technology is introduced to make sure the privacy protection of the semantic information of training dataset. Furthermore, we leverage FL to allow hospitals to collaboratively train COVID-19 models without sharing the original data. Under Independent and Identically Distributed (IID) and non-IID settings, The evaluation of the proposed model is on three types of chest X-ray (CXR) images dataset (COVID-19, normal, and normal pneumonia). A large number of the truthful reports make the verification of our model can effectively diagnose COVID-19 without compromising privacy.


title: Understanding and Avoiding AI Failures: A Practical Guide

abstract: As AI technologies increase in capability and ubiquity, AI accidents are becoming more common. Based on normal accident theory, high reliability theory, and open systems theory, we create a framework for understanding the risks associated with AI applications. In addition, we also use AI safety principles to quantify the unique risks of increased intelligence and human-like qualities in AI. Together, these two fields give a more complete picture of the risks of contemporary AI. By focusing on system properties near accidents instead of seeking a root cause of accidents, we identify where attention should be paid to safety for current generation AI systems.


title: Vision-based Driver Assistance Systems: Survey, Taxonomy and Advances

abstract: Vision-based driver assistance systems is one of the rapidly growing research areas of ITS, due to various factors such as the increased level of safety requirements in automotive, computational power in embedded systems, and desire to get closer to autonomous driving. It is a cross disciplinary area encompassing specialised fields like computer vision, machine learning, robotic navigation, embedded systems, automotive electronics and safety critical software. In this paper, we survey the list of vision based advanced driver assistance systems with a consistent terminology and propose a taxonomy. We also propose an abstract model in an attempt to formalize a top-down view of application development to scale towards autonomous driving system.


title: Localization formulas of cohomology intersection numbers

abstract: We revisit the localization formulas of cohomology intersection numbers associated to a logarithmic connection. The main contribution of this paper is threefold: we prove the localization formula of the cohomology intersection number of logarithmic forms in terms of residue of a connection; we prove that the leading term of the Laurent expansion of the cohomology intersection number is Grothendieck residue when the connection is hypergeometric; and we prove that the leading term of stringy integral discussed by Arkani-Hamed, He and Lam is nothing but the self-cohomology intersection number of the canonical form.


title: Effect of collisions on non-adiabatic electron dynamics in ITG-driven microturbulence

abstract: Non-adiabatic electron response leads to significant changes in Ion Temperature Gradient (ITG) eigenmodes, leading in particular to fine-structures that are significantly extended along the magnetic field lines at corresponding Mode Rational Surfaces (MRSs). These eigenmodes can nonlinearly interact with themselves to drive zonal flows via the so-called self-interaction mechanism. In this paper, the effect of collisions on these processes are studied. In presence of non-adiabatic electrons, the linear growth rate of ITG eigenmodes decreases with increasing collisionality. Detailed velocity space analysis of the distribution function shows that this results from collisions leading to a more adiabatic-like response of electrons away from MRSs. In linear simulations, collisions are furthermore found to broaden the radial width of the fine-structures, which translates to narrower tails of the eigenmode in extended ballooning space. The characteristic parallel scale length associated to these tails is shown to scale with the mean free path of electrons. In nonlinear turbulence simulations accounting for physically relevant values of collisionality, the fine structures located at MRSs, together with the associated drive of zonal flows via self-interaction, are shown to persist and play a significant role.


title: Consistency issues in Gaussian Mixture Models reduction algorithms

abstract: In many contexts Gaussian Mixtures (GM) are used to approximate probability distributions, possibly time-varying. In some applications the number of GM components exponentially increases over time, and reduction procedures are required to keep them reasonably limited. The GM reduction (GMR) problem can be formulated by choosing different measures of the dissimilarity of GMs before and after reduction, like the Kullback-Leibler Divergence (KLD) and the Integral Squared Error (ISE). Since in no case the solution is obtained in closed form, many approximate GMR algorithms have been proposed in the past three decades, although none of them provides optimality guarantees. In this work we discuss the importance of the choice of the dissimilarity measure and the issue of consistency of all steps of a reduction algorithm with the chosen measure. Indeed, most of the existing GMR algorithms are composed by several steps which are not consistent with a unique measure, and for this reason may produce reduced GMs far from optimality. In particular, the use of the KLD, of the ISE and normalized ISE is discussed and compared in this perspective.


title: Bayesian Numerical Methods for Nonlinear Partial Differential Equations

abstract: The numerical solution of differential equations can be formulated as an inference problem to which formal statistical approaches can be applied. However, nonlinear partial differential equations (PDEs) pose substantial challenges from an inferential perspective, most notably the absence of explicit conditioning formula. This paper extends earlier work on linear PDEs to a general class of initial value problems specified by nonlinear PDEs, motivated by problems for which evaluations of the right-hand-side, initial conditions, or boundary conditions of the PDE have a high computational cost. The proposed method can be viewed as exact Bayesian inference under an approximate likelihood, which is based on discretisation of the nonlinear differential operator. Proof-of-concept experimental results demonstrate that meaningful probabilistic uncertainty quantification for the unknown solution of the PDE can be performed, while controlling the number of times the right-hand-side, initial and boundary conditions are evaluated. A suitable prior model for the solution of the PDE is identified using novel theoretical analysis of the sample path properties of Mat\'{e}rn processes, which may be of independent interest.


title: Algorithms for finding generalized minimum aberration designs

abstract: Statistical design of experiments is widely used in scientific and industrial investigations. A generalized minimum aberration (GMA) orthogonal array is optimum under the well-established, so-called GMA criterion, and such an array can extract as much information as possible at a fixed cost. Finding GMA arrays is an open (yet fundamental) problem in design of experiments because constructing such arrays becomes intractable as the number of runs and factors increase. We develop two directed enumeration algorithms that call the integer programming with isomorphism pruning algorithm of Margot (2007) for the purpose of finding GMA arrays. Our results include 16 GMA arrays that were not previously in the literature, along with documentation of the efficiencies that made the required calculations possible within a reasonable budget of computer time. We also validate heuristic algorithms against a GMA array catalog, by showing that they quickly output near GMA arrays, and then use the heuristics to find near GMA arrays when enumeration is computationally burdensome.


title: Exploiting Transitivity Constraints for Entity Matching in Knowledge Graphs

abstract: The goal of entity matching in knowledge graphs is to identify entities that refer to the same real-world objects using some similarity metric. The result of entity matching can be seen as a set of entity pairs interpreted as the same-as relation. However, the identified set of pairs may fail to satisfy some structural properties, in particular transitivity, that are expected from the same-as relation. In this work, we show that an ad-hoc enforcement of transitivity, i.e. taking the transitive closure, on the identified set of entity pairs may decrease precision dramatically. We therefore propose a methodology that starts with a given similarity measure, generates a set of entity pairs that are identified as referring to the same real-world objects, and applies the cluster editing algorithm to enforce transitivity without adding many spurious links, leading to overall improved performance.


title: High-field specific heat and entropy obtained from adiabatic temperature change

abstract: Specific heat and entropy are relevant thermodynamic properties, which may be used as macroscopic probes to microscopic properties of materials under ambient conditions and under high applied fields. However, the measurement of specific heat under intense external fields can be a challenging task, as well as to obtain the entropy in the same conditions. Here, we describe a method to obtain high-field specific heat and entropy from measurements of specific heat under ambient conditions and direct temperature change induced by adiabatic field changes. We derive straightforward thermodynamic equations to calculate the specific heat and entropy and our results agree satisfactorily with experimental data of specific heat under magnetic field, electric field, and pressure.


title: Predictive analytics using Social Big Data and machine learning

abstract: The ever-increase in the quality and quantity of data generated from day-to-day businesses operations in conjunction with the continuously imported related social data have made the traditional statistical approaches inadequate to tackle such data floods. This has dictated researchers to design and develop advance and sophisticated analytics that can be incorporated to gain valuable insights that benefit the business domain. This chapter sheds the light on core aspects that lay the foundations for social big data analytics. In particular, the significance of predictive analytics in the context of SBD is discussed fortified with presenting a framework for SBD predictive analytics. Then, various predictive analytical algorithms are introduced with their usage in several important application and top-tier tools and APIs. A case study on using predictive analytics to social data is provided supported with experiments to substantiate the significance and utility of predictive analytics.


title: Understanding and Accelerating EM Algorithm's Convergence by Fair Competition Principle and Rate-Verisimilitude Function

abstract: Why can the Expectation-Maximization (EM) algorithm for mixture models converge? Why can different initial parameters cause various convergence difficulties? The Q-L synchronization theory explains that the observed data log-likelihood L and the complete data log-likelihood Q are positively correlated; we can achieve maximum L by maximizing Q. According to this theory, the Deterministic Annealing EM (DAEM) algorithm's authors make great efforts to eliminate locally maximal Q for avoiding L's local convergence. However, this paper proves that in some cases, Q may and should decrease for L to increase; slow or local convergence exists only because of small samples and unfair competition. This paper uses marriage competition to explain different convergence difficulties and proposes the Fair Competition Principle (FCP) with an initialization map for improving initializations. It uses the rate-verisimilitude function, extended from the rate-distortion function, to explain the convergence of the EM and improved EM algorithms. This convergence proof adopts variational and iterative methods that Shannon et al. used for analyzing rate-distortion functions. The initialization map can vastly save both algorithms' running times for binary Gaussian mixtures. The FCP and the initialization map are useful for complicated mixtures but not sufficient; we need further studies for specific methods.


title: On sums of Fibonacci numbers with few binary digits

abstract: In this paper we completely solve the Diophantine equation $F_n+F_m=2^{a_1}+2^{a_2}+2^{a_3}+2^{a_4}+2^{a_5}$, where $F_k$ denotes the $k$-th Fibonacci number. In addition to complex linear forms in logarithms and the Baker-Davenport reduction method, we use $p$-adic versions of both tools.


title: MAQ-CaF: A Modular Air Quality Calibration and Forecasting method for cross-sensitive pollutants

abstract: The climatic challenges are rising across the globe in general and in worst hit under-developed countries in particular. The need for accurate measurements and forecasting of pollutants with low-cost deployment is more pertinent today than ever before. Low-cost air quality monitoring sensors are prone to erroneous measurements, frequent downtimes, and uncertain operational conditions. Such a situation demands a prudent approach to ensure an effective and flexible calibration scheme. We propose MAQ-CaF, a modular air quality calibration, and forecasting methodology, that side-steps the challenges of unreliability through its modular machine learning-based design which leverages the potential of IoT framework. It stores the calibrated data both locally and remotely with an added feature of future predictions. Our specially designed validation process helps to establish the proposed solution's applicability and flexibility without compromising accuracy. CO, SO2, NO2, O3, PM1.0, PM2.5 and PM10 were calibrated and monitored with reasonable accuracy. Such an attempt is a step toward addressing climate change's global challenge through appropriate monitoring and air quality tracking across a wider geographical region via affordable monitoring.


title: Topics in global real analytic geometry

abstract: In the first half of twentieth century the theory of complex analytic functions and of their zerosets was fully developed. The definition of holomorphic function has a local nature. Germs of holomorphic functions form a distinguished subring of the ring of germs of continuous functions. This way came out the notion of analytic space. The definition of a complex analytic set is by local models as in the case of complex manifolds. But while local models for manifolds are open sets of ${\mathbb C}^n$, a local model of an analytic space is the zeroset of finitely many analytic functions on an open set of ${\mathbb C}^n$ together with a sheaf of continuous function to be called holomorphic. Towards the years $50$ of the last century, Cartan, Whitney, Bruhat and others tried to formulate the notion of analytic space over ${\mathbb R}$. Immediately they realize that the real sets verifying a definition similar to the complex one form a cathegory whose elements can get unpleasant behaviour. In particular this cathegory does not get the good properties of complex analytic spaces, as for instance coherence of their structural sheaves and Theorems A and B do not hold in general.


title: Non-Gaussian Quantum States and Where to Find Them

abstract: Gaussian states have played on important role in the physics of continuous-variable quantum systems. They are appealing for the experimental ease with which they can be produced, and for their compact and elegant mathematical description. Nevertheless, many proposed quantum technologies require us to go beyond the realm of Gaussian states and introduce non-Gaussian elements. In this Tutorial, we provide a roadmap for the physics of non-Gaussian quantum states. We introduce the phase-space representations as a framework to describe the different properties of quantum states in continuous-variable systems. We then use this framework in various ways to explore the structure of the state space. We explain how non-Gaussian states can be characterised not only through the negative values of their Wigner function, but also via other properties such as quantum non-Gaussianity and the related stellar rank. For multimode systems, we are naturally confronted with the question of how non-Gaussian properties behave with respect to quantum correlations. To answer this question, we first show how non-Gaussian states can be created by performing measurements on a subset of modes in a Gaussian state. Then, we highlight that these measured modes must be correlated via specific quantum correlations to the remainder of the system to create quantum non-Gaussian or Wigner-negative states. On the other hand, non-Gaussian operations are also shown to enhance or even create quantum correlations. Finally, we will demonstrate that Wigner negativity is a requirement to violate Bell inequalities and to achieve a quantum computational advantage. At the end of the Tutorial, we also provide an overview of several experimental realisations of non-Gaussian quantum states in quantum optics and beyond.


title: Valid Heteroskedasticity Robust Testing

abstract: Tests based on heteroskedasticity robust standard errors are an important technique in econometric practice. Choosing the right critical value, however, is not simple at all: conventional critical values based on asymptotics often lead to severe size distortions; and so do existing adjustments including the bootstrap. To avoid these issues, we suggest to use smallest size-controlling critical values, the generic existence of which we prove in this article for the commonly used test statistics. Furthermore, sufficient and often also necessary conditions for their existence are given that are easy to check. Granted their existence, these critical values are the canonical choice: larger critical values result in unnecessary power loss, whereas smaller critical values lead to over-rejections under the null hypothesis, make spurious discoveries more likely, and thus are invalid. We suggest algorithms to numerically determine the proposed critical values and provide implementations in accompanying software. Finally, we numerically study the behavior of the proposed testing procedures, including their power properties.


title: Gaussian complex zeroes are not always normal: limit theorems on the disc

abstract: We study the zeroes of a family of random holomorphic functions on the unit disc, distinguished by their invariance with respect to the hyperbolic geometry. Our main finding is a transition in the limiting behaviour of the number of zeroes in a large hyperbolic disc. We find a normal distribution if the covariance decays faster than a certain critical value. In contrast, in the regime of 'long-range dependence' when the covariance decays slowly, the limiting distribution is skewed. For a closely related model we emphasise a link with Gaussian multiplicative chaos.


title: Radiation damage effects on helium diffusion in zircon

abstract: We report the effects of radiation damage on helium diffusion in zircon using data from molecular dynamics simulations. We observe an increase in activation energy for helium diffusion as a result of radiation damage and increasing structural disorder. The activation energy in a heavily damaged region is smaller than in a completely amorphous system which is correlated with remaining order in the cation sublattices of the damaged structure not present in the fully amorphized system. The reduction of activation energy is related to the disappearance of fast diffusion pathways that present in the crystal. Consistent with the change in activation energy, we observe the accumulation of helium atoms in the damaged structure and discuss the implications of this effect for the formation of helium bubbles and zircon's performance as an encapsulation material for nuclear waste.


title: Direct numerical simulation of two-phase pipe flow: influence of the domain length on the flow regime

abstract: We perform direct numerical simulations of a kerosene-water mixture in pipe flow under realistic experimental conditions by solving the Cahn-Hilliard-Navier-Stokes equations. We compute the linear stability of core-annular flow of kerosene and water in a vertical pipe and find that it is highly unstable. By performing DNS initialized with a slightly perturbed core-annular flow, we show that the system transitions to turbulence and finally relaxes into a turbulent slug flow regime provided that the pipe is sufficiently long. This configuration presents mild turbulence and large scale three-dimensional recirculation patterns inside the slugs. Our work highlights the need for applying nonlinear-dynamics approaches and carefully selecting the domain length to investigate the patterns observed in two-phase pipe flows and demonstrates the capabilities of phase field methods to reliably simulate experimental conditions.