title: Rayleigh-B\'enard convection of a model emulsion: anomalous heat-flux fluctuations and finite-size droplets effects

abstract: We present mesoscale numerical simulations of Rayleigh-B\'enard (RB) convection in a two-dimensional model emulsion. The systems under study are constituted of finite-size droplets, whose concentration Phi_0 is systematically varied from small (Newtonian emulsions) to large values (non-Newtonian emulsions). We focus on the characterisation of the heat transfer properties close to the transition from conductive to convective states, where it is known that a homogeneous Newtonian system exhibits a steady flow and a time-independent heat flux. In marked contrast, emulsions exhibit a non-steady dynamics with fluctuations in the heat flux. In this paper, we aim at the characterisation of such non-steady dynamics via detailed studies on the time-averaged heat flux and its fluctuations. To understand the time-averaged heat flux, we propose a side-by-side comparison between the emulsion system and a single-phase (SP) system, whose viscosity is constructed from the shear rheology of the emulsion. We show that such local closure works well only when a suitable degree of coarse-graining (at the droplet scale) is introduced in the local viscosity. To delve deeper into the fluctuations in the heat flux, we propose a side-by-side comparison between a Newtonian emulsion and a non-Newtonian emulsion, at fixed time-averaged heat flux. This comparison elucidates that finite-size droplets and the non-Newtonian rheology cooperate to trigger enhanced heat-flux fluctuations at the droplet scales. These enhanced fluctuations are rooted in the emergence of space correlations among distant droplets, which we highlight via direct measurements of the droplets displacement and the characterisation of the associated correlation function. The observed findings offer insights on heat transfer properties for confined systems possessing finite-size constituents.


title: A new route for enantio-sensitive structure determination by photoelectron scattering on molecules in the gas phase

abstract: X-ray as well as electron diffraction are powerful tools for structure determination of molecules. Studies on randomly oriented molecules in the gas-phase address cases in which molecular crystals cannot be generated or the interaction-free molecular structure is to be addressed. Such studies usually yield partial geometrical information, such as interatomic distances. Here, we present a complementary approach, which allows obtaining insight to the structure, handedness and even detailed geometrical features of molecules in the gas phase. Our approach combines Coulomb explosion imaging, the information that is encoded in the molecular frame diffraction pattern of core-shell photoelectrons and ab initio computations. Using a loop-like analysis scheme we are able to deduce specific molecular coordinates with sensitivity even to the handedness of chiral molecules and the positions of individual atoms, as, e.g., protons.


title: Higher-derivative Heterotic Double Field Theory and Classical Double Copy

abstract: The generalized Kerr-Schild ansatz (GKSA) is a powerful tool for constructing exact solutions in Double Field Theory (DFT). In this paper we focus in the heterotic formulation of DFT, considering up to four-derivative terms in the action principle, while the field content is perturbed by the GKSA. We study the inclusion of the generalized version of the Green-Schwarz mechanism to this setup, in order to reproduce the low energy effective heterotic supergravity upon parametrization. This formalism reproduces higher-derivative heterotic background solutions where the metric tensor and Kalb-Ramond field are perturbed by a pair of null vectors. Next we study higher-derivative contributions to the classical double copy structure. After a suitable identification of the null vectors with a pair of $U(1)$ gauge fields, the dynamics is given by a pair of Maxwell equations plus higher derivative corrections in agreement with the KLT relation.


title: Iterative Power Algorithm for Global Optimization with Quantics Tensor Trains

abstract: Optimization algorithms play a central role in chemistry since optimization is the computational keystone of most molecular and electronic structure calculations. Herein, we introduce the iterative power algorithm (IPA) for global optimization and a formal proof of convergence for both discrete and continuous global search problems, which is essential for applications in chemistry such as molecular geometry optimization. IPA implements the power iteration method in quantics tensor train (QTT) representations. Analogous to the imaginary time propagation method with infinite mass, IPA starts with an initial probability distribution and iteratively applies the recurrence relation $\rho_{k+1}({\bf x})=U({\bf x})\rho_k({\bf x})/\|U\rho_k\|_{L^1}$, where $U({\bf x})=e^{-V({\bf x})}$ is defined in terms of the potential energy surface (PES) $V({\bf x})$. Upon convergence, the probability distribution becomes a delta function, so the global minimum can be obtained as the position expectation value. QTT representations are generated by fast adaptive interpolation of multidimensional arrays to bypass the curse of dimensionality and the need to evaluate $V({\bf x})$ for all possible ${\bf x}$. We illustrate the capabilities of IPA for global search optimization of two multidimensional PESs, including a differentiable model PES of a DNA chain and a discrete non-differentiable PES, $V(p)=\text{mod}(N,p)$, that resolves the prime factors of an integer $N$, with $p$ in the space of prime numbers folded as a $d$-dimensional $2_1 \times 2_2 \times \cdots \times 2_d$ tensor. We find that IPA resolves multiple degenerate global minima even when separated by large energy barriers in the highly rugged landscape of the potentials. Therefore, IPA should be of great interest for a wide range of other optimization problems ubiquitous in molecular and electronic structure calculations.


title: Designing optimal masks for a multi-object spectrometer

abstract: This paper concerns a new optimization problem arising in the management of a multi-object spectrometer with a configurable slit unit. The field of view of the spectrograph is divided into contiguous and parallel spatial bands, each one associated with two opposite sliding metal bars that can be positioned to observe one astronomical object. Thus several objects can be analyzed simultaneously within a configuration of the bars called a mask. Due to the high demand from astronomers, pointing the spectrograph's field of view to the sky, rotating it, and selecting the objects to conform a mask is a crucial optimization problem for the efficient use of the spectrometer. The paper describes this problem, presents a Mixed Integer Linear Programming formulation for the case where the rotation angle is fixed, presents a non-convex formulation for the case where the rotation angle is unfixed, describes a heuristic approach for the general problem, and discusses computational results on real-world and randomly-generated instances.

categories: astro-ph.IM astro-ph.EP astro-ph.GA math.CO


title: Quaternionic quantum harmonic oscillator

abstract: In this article we obtained the harmonic oscillator solution for quaternionic quantum mechanics ($\mathbbm{H}$QM) in the real Hilbert space, both in the analytic method and in the algebraic method. The quaternionic solutions have many additional possibilities if compared to complex quantum mechanics ($\mathbbm{C}$QM), and thus there are many possible applications to these results in future research.


title: On the Implementation Complexity of Digital Full-Duplex Self-Interference Cancellation

abstract: In-band full-duplex systems promise to further increase the throughput of wireless systems, by simultaneously transmitting and receiving on the same frequency band. However, concurrent transmission generates a strong self-interference signal at the receiver, which requires the use of cancellation techniques. A wide range of techniques for analog and digital self-interference cancellation have already been presented in the literature. However, their evaluation focuses on cases where the underlying physical parameters of the full-duplex system do not vary significantly. In this paper, we focus on adaptive digital cancellation, motivated by the fact that physical systems change over time. We examine some of the different cancellation methods in terms of their performance and implementation complexity, considering the cost of both cancellation and training. We then present a comparative analysis of all these methods to determine which perform better under different system performance requirements. We demonstrate that with a neural network approach, the reduction in arithmetic complexity for the same cancellation performance relative to a state-of-the-art polynomial model is several orders of magnitude.


title: Optimal operating policies for organic Rankine cycles for waste heat recovery under transient conditions

abstract: Waste heat recovery for trucks via organic Rankine cycle is a promising technology to reduce fuel consumption and emissions. As the vehicles are operated in street traffic, the heat source is subject to strong fluctuations. Consequently, such disturbances have to be considered to enable safe and efficient operation. Herein, we find optimal operating policies for several representative scenarios by means of dynamic optimization and discuss the implications on control strategy design. First, we optimize operation of a typical driving cycle with data from a test rig. Results indicate that operating the cycle at minimal superheat is an appropriate operating policy. Second, we consider a scenario where the permissible expander power is temporarily limited, which is realistic in street traffic. In this case, an operating policy with flexible superheat can reduce the losses associated with operation at minimal superheat by up to 53 % in the considered scenario. As the duration of power limitation increases, other constraints might become active which results in part of the exhaust gas being bypassed, hence reduced savings.


title: Task Adaptive Pretraining of Transformers for Hostility Detection

abstract: Identifying adverse and hostile content on the web and more particularly, on social media, has become a problem of paramount interest in recent years. With their ever increasing popularity, fine-tuning of pretrained Transformer-based encoder models with a classifier head are gradually becoming the new baseline for natural language classification tasks. In our work, we explore the gains attributed to Task Adaptive Pretraining (TAPT) prior to fine-tuning of Transformer-based architectures. We specifically study two problems, namely, (a) Coarse binary classification of Hindi Tweets into Hostile or Not, and (b) Fine-grained multi-label classification of Tweets into four categories: hate, fake, offensive, and defamation. Building up on an architecture which takes emojis and segmented hashtags into consideration for classification, we are able to experimentally showcase the performance upgrades due to TAPT. Our system (with team name 'iREL IIIT') ranked first in the 'Hostile Post Detection in Hindi' shared task with an F1 score of 97.16% for coarse-grained detection and a weighted F1 score of 62.96% for fine-grained multi-label classification on the provided blind test corpora.


title: Topological susceptibility of pure gauge theory using Density of States

abstract: The topological susceptibility of the SU(3) pure gauge theory is calculated in the deconfined phase at temperatures up to $10T_c$. At such large temperatures the susceptibility is suppressed, topologically non-trivial configurations are extremely rare. Thus, direct lattice simulations are not feasible. The density of states (DoS) method is designed to simulate rare events, we present an application of the DoS method to the problem of high temperature topological susceptibility. We reconstruct the histogram of the charge sectors that one could have obtained in a naive importance sampling. Our findings are perfectly consistent with a free instanton gas.


title: On inclusions between quantified provability logics

abstract: We investigate several consequences of inclusion relations between quantified provability logics. Moreover, we give a necessary and sufficient condition for the inclusion relation between quantified provability logics with respect to $\Sigma_1$ arithmetical interpretations.


title: Critical exponent $\eta$ at $O(1/N^3)$ in the chiral XY model using the large $N$ conformal bootstrap

abstract: We compute the $O(1/N^3)$ correction to the critical exponent $\eta$ in the chiral XY or chiral Gross-Neveu model in $d$-dimensions. As the leading order vertex anomalous dimension vanishes, the direct application of the large $N$ conformal bootstrap formalism is not immediately possible. To circumvent this we consider the more general Nambu-Jona-Lasinio model for a general non-abelian Lie group. Taking the abelian limit of the exponents of this model produces those of the chiral XY model. Subsequently we provide improved estimates for $\eta$ in the three dimensional chiral XY model for various values of $N$.


title: Cluster and conquer: The morphodynamics of invasion of a compliant substrate by active rods

abstract: The colonisation of a soft passive material by motile cells such as bacteria is common in biology. The resulting colonies of the invading cells are often observed to exhibit intricate patterns whose morphology and dynamics can depend on a number of factors, particularly the mechanical properties of the substrate and the motility of the individual cells. We use simulations of a minimal 2D model of self-propelled rods moving through with a passive compliant medium consisting of particles that offer elastic resistance before being plastically displaced from their equilibrium positions. It is observed that the motility-induced clustering of active (self-propelled) particles is crucial for understanding the morphodynamics of colonisation. Clustering enables motile colonies to spread faster than they would have as isolated particles. The colonisation rate depends non-monotonically on substrate stiffness with a distinct maximum at a non-zero value of substrate stiffness. This is observed to be due to a change in the morphology of clusters. Furrow networks created by the active particles have a fractal-like structure whose dimension varies systematically with substrate stiffness but is less sensitive to particle activity. The power-law growth exponent of the furrowed area is smaller than unity, suggesting that, to sustain such extensive furrow networks, colonies must regulate their overall growth rate.


title: Connection between inverse engineering and optimal control in shortcuts to adiabaticity

abstract: We consider fast high-fidelity quantum control by using a shortcut to adiabaticity (STA) technique and optimal control theory (OCT). Three specific examples, including expansion of cold atoms from the harmonic trap, atomic transport by moving harmonic trap, and spin dynamics in the presence of dissipation, are explicitly detailed. Using OCT as a qualitative guide, we demonstrate how STA protocols designed from inverse engineering method, can approach with very high precision optimal solutions built about physical constraints, by a proper choice of the interpolation function and with a very reduced number of adjustable parameters.


title: An Optimization Framework for Power Infrastructure Planning

abstract: The ubiquitous expansion and transformation of the energy supply system involves large-scale power infrastructure construction projects. In the view of investments of more than a million dollars per kilometre, planning authorities aim to minimise the resistances posed by multiple stakeholders. Mathematical optimisation research offers efficient algorithms to compute globally optimal routes based on geographic input data. We propose a framework that utilizes a graph model where vertices represent possible locations of transmission towers, and edges are placed according to the feasible distance between neighbouring towers. In order to cope with the specific challenges arising in linear infrastructure layout, we first introduce a variant of the Bellman-Ford algorithm that efficiently computes the minimal-angle shortest path. Secondly, an iterative procedure is proposed that yields a locally optimal path at considerably lower memory requirements and runtime. Third, we discuss and analyse methods to output k diverse path alternatives. Experiments on real data show that compared to previous work, our approach reduces the resistances by more than 10% in feasible time, while at the same time offering much more flexibility and functionality. Our methods are demonstrated in a simple and intuitive graphical user interface, and an open-source package (LION) is available at https://pypi.org/project/lion-sp.


title: Equalized Recovery State Estimators for Linear Systems with Delayed and Missing Observations

abstract: This paper presents a dynamic state observer design for discrete-time linear time-varying systems that robustly achieves equalized recovery despite delayed or missing observations, where the set of all temporal patterns for the missing or delayed data is modeled by a finite-length language. By introducing a mapping of the language onto a reduced event-based language, we design a state estimator that adapts based on the history of available data at each step, and satisfies equalized recovery for all patterns in the reduced language. In contrast to existing equalized recovery estimators, the proposed design considers the equalized recovery level as a decision variable, which enables us to directly obtain the global minimum for the intermediate recovery level, resulting in improved estimation performance. Finally, we demonstrate the effectiveness of the proposed observer when compared to existing approaches using several illustrative examples.


title: On the suboptimality of the p-version discontinuous Galerkin methods for first order hyperbolic problems

abstract: We address the issue of the suboptimality in the p-version discontinuous Galerkin (dG) methods for first order hyperbolic problems. The convergence rate is derived for the upwind dG scheme on tensor product meshes in any dimension. The standard proof in seminal work [14] leads to suboptimal convergence in terms of the polynomial degree by 3/2 order for general convection fields, with the exception of piecewise multi-linear convection fields, which rather yield optimal convergence. Such suboptimality is not observed numerically. Thus, it might be caused by a limitation of the analysis, which we partially overcome: for a special class of convection fields, we shall show that the dG method has a p-convergence rate suboptimal by 1/2 order only.


title: Paradoxes of Probabilistic Programming

abstract: Probabilistic programming languages allow programmers to write down conditional probability distributions that represent statistical and machine learning models as programs that use observe statements. These programs are run by accumulating likelihood at each observe statement, and using the likelihood to steer random choices and weigh results with inference algorithms such as importance sampling or MCMC. We argue that naive likelihood accumulation does not give desirable semantics and leads to paradoxes when an observe statement is used to condition on a measure-zero event, particularly when the observe statement is executed conditionally on random data. We show that the paradoxes disappear if we explicitly model measure-zero events as a limit of positive measure events, and that we can execute these type of probabilistic programs by accumulating infinitesimal probabilities rather than probability densities. Our extension improves probabilistic programming languages as an executable notation for probability distributions by making it more well-behaved and more expressive, by allowing the programmer to be explicit about which limit is intended when conditioning on an event of measure zero.


title: Generate Natural Language Explanations for Recommendation

abstract: Providing personalized explanations for recommendations can help users to understand the underlying insight of the recommendation results, which is helpful to the effectiveness, transparency, persuasiveness and trustworthiness of recommender systems. Current explainable recommendation models mostly generate textual explanations based on pre-defined sentence templates. However, the expressiveness power of template-based explanation sentences are limited to the pre-defined expressions, and manually defining the expressions require significant human efforts. Motivated by this problem, we propose to generate free-text natural language explanations for personalized recommendation. In particular, we propose a hierarchical sequence-to-sequence model (HSS) for personalized explanation generation. Different from conventional sentence generation in NLP research, a great challenge of explanation generation in e-commerce recommendation is that not all sentences in user reviews are of explanation purpose. To solve the problem, we further propose an auto-denoising mechanism based on topical item feature words for sentence generation. Experiments on various e-commerce product domains show that our approach can not only improve the recommendation accuracy, but also the explanation quality in terms of the offline measures and feature words coverage. This research is one of the initial steps to grant intelligent agents with the ability to explain itself based on natural language sentences.


title: Elastic 3D-2D Image Registration

abstract: We propose a method to non-rigidly align a three-dimensional (3D) volumetric image with a two-dimensional (2D) planar image representing a projection of the deformed volume. The application in mind comes from biological studies in which 2D intravital microscopy videos of living tissue are recorded, after which the tissue is excised and a more detailed 3D volume microscopy is performed. Coregistration of both data sets allows to combine the temporal (but 2D) information with more detailed spatial 3D information. Our approach is variational and uses a hyperelastic deformation regularization, as is appropriate for biological material. As a particular feature, the out of plane deformation is estimated based on the out of focus blur inside the 2D microscopy image. The approach becomes computationally feasible through the use of a coarse-to-fine optimization strategy and higher order optimization methods.


title: Context-Aware Target Apps Selection and Recommendation for Enhancing Personal Mobile Assistants

abstract: Users install many apps on their smartphones, raising issues related to information overload for users and resource management for devices. Moreover, the recent increase in the use of personal assistants has made mobile devices even more pervasive in users' lives. This paper addresses two research problems that are vital for developing effective personal mobile assistants: target apps selection and recommendation. The former is the key component of a unified mobile search system: a system that addresses the users' information needs for all the apps installed on their devices with a unified mode of access. The latter, instead, predicts the next apps that the users would want to launch. Here we focus on context-aware models to leverage the rich contextual information available to mobile devices. We design an in situ study to collect thousands of mobile queries enriched with mobile sensor data (now publicly available for research purposes). With the aid of this dataset, we study the user behavior in the context of these tasks and propose a family of context-aware neural models that take into account the sequential, temporal, and personal behavior of users. We study several state-of-the-art models and show that the proposed models significantly outperform the baselines.


title: Stable solution of the log-Minkowski problem in the case of many hyperplane symmetries

abstract: In the case of symmetries with respect to n independent linear hyperplanes, the stability of the solution of the Logarithmic Minkowski problem on S^{n-1} is established.


title: A Data-Taking System for Planetary Radar Applications

abstract: Most planetary radar applications require recording of complex voltages at sampling rates of up to 20 MHz. I describe the design and implementation of a sampling system that has been installed at the Arecibo Observatory, Goldstone Solar System Radar, and Green Bank Telescope. After many years of operation, these data-taking systems have enabled the acquisition of hundreds of data sets, many of which still await publication.


title: Isomonodromic Laplace Transform with Coalescing Eigenvalues and Confluence of Fuchsian Singularities

abstract: We consider a Pfaffian system expressing isomonodromy of an irregular system of Okubo type, depending on complex deformation parameters u=(u_1,...,u_n), which are eigenvalues of the leading matrix at the irregular singuilarity. At the same time, we consider a Pfaffian system of non-normalized Schlesinger type expressing isomonodromy of a Fuchsian system, whose poles are the deformation parameters u_1,...,u_n. The parameters vary in a polydisc containing a coalescence locus for the eigenvalues of the leading matrix of the irregular system, corresponding to confluence of the Fuchsian singularities. We construct isomonodromic selected and singular vector solutions of the Fuchsian Pfaffian system together with their isomonodromic connection coefficients, so extending a result of references [4] and [20] to the isomonodromic case, including confluence of singularities. Then, we introduce an isomonodromic Laplace transform of the selected and singular vector solutions, allowing to obtain isomonodromic fundamental solutions for the irregular system, and their Stokes matrices expressed in terms of connection coefficients. These facts, in addition to extending [4] and [20] to the isomonodromic case with coalescences/confluences, allow to prove by means of Laplace transform the main result of reference [11], which is the analytic theory of non-generic isomonodromic deformations of the irregular system with coalescing eigenvalues.


title: Electrically non-neutral ground states of stars

abstract: To approximately compute the non-relativistic ground state of an electrically non-neutral star, an exactly solvable model was recently introduced, and partly solved, by Krivoruchenko, Nadyozhin, and Yudin. The model generalizes the well-known Lane--Emden equation of a polytropic gas ball of index $n=1$ to a two-fluid setting. Here its complete solution is presented in terms of simple elementary functions; it is also generalized to a more-than-two-fluid setting where it remains exactly solvable. It is shown that, given the number of nuclei, a maximal negatively and a maximal positively charged solution exists, plus a continuous family of solutions which interpolates between these extremes. Numerical comparisons show that this exactly solvable model captures the qualitative behavior of the more physical model it is supposed to approximate. Furthermore, it correctly answers the question: how non-neutral can the star be? The answer is independent of the speed of light $c$ and the Planck quantum $h$. It supports Penrose's weak cosmic censorship hypothesis, in the sense that the bounds on the excess charge are compatible with the bound on the charge of a Reissner--Weyl--Nordstr\"om black hole.


title: Surface state at $BaSnO_3$ evidenced by angle-resolved photoemission spectroscopy and ab initio calculations

abstract: Perovskite alkaline earth stannates, such as $BaSnO_3$ and $SrSnO_3$, showing light transparency and high electrical conductivity (when doped), have become promising candidates for novel optoelectrical devices. Such devices are mostly based on hetero-structures and understanding of their electronic structure, which usually deviates from the bulk, is mandatory for exploring a full application potential. Employing angle-resolved photoemission spectroscopy and ab initio calculations we reveal the existence of a 2-dimensional metallic state at the $SnO_2$-terminated surface of a 1\% La-doped $BaSnO_3$ thin film. The observed surface state is characterized by distinct carrier density and a smaller effective mass in comparison with the corresponding bulk values. The small surface effective mass of about $0.12m_e$ can cause an improvement of the electrical conductivity of BSO based heterostructures.


title: Asymptotic analysis of a biphase tumor fluid flow. The weak coupling case

abstract: The aim of this paper is to investigate the asymptotic behavior of a biphase tumor fluid flow derived by 2-scale homogenisation techniques in recent works. This biphase fluid flow model accounts for the capillary wall permeability, and the interstitial avascular phase, both being mixed in the limit homogenised problem. When the vessel walls become more permeable, we show that the biphase fluid flow exhibits a boundary layer that makes the computation of the full problem costly and unstable. In the limit, both capillary and interstitial pressures coincide except in the vicinity of the boundary where different boundary conditions are applied. Thanks to a rigorous asymptotic analysis, we prove that the solution to the full problem can be approached at any order of approximation by a monophasic model with appropriate boundary conditions on the tumor boundary and appropriate correcting terms near the boundary are given. Numerical simulations in spherical geometry illustrate the theoretical results.


title: Analysis of 3GM Callisto Gravity Experiment of the JUICE Mission

abstract: The ESA's JUICE mission will provide a thorough investigation of the Jupiter system and the Galilean moons during its nominal tour, comprising flybys of Europa and Callisto, and an orbital phase about Ganymede at the end of the mission. The 3GM experiment will exploit accurate Doppler and range measurements to determine the moons' orbits and gravity fields (both static and tidal) and infer their interior structure. This paper presents the attainable accuracies of the Callisto geodesy experiment and addresses the effect of different flybys mean anomaly distribution and geometry on the estimation of the tidal Love number k2.

categories: astro-ph.EP astro-ph.IM math.OC physics.space-ph


title: Summability characterizations of positive sequences

abstract: In this paper, we propose extensions for the classical Kummer test, which is a very far-reaching criterion that provides sufficient and necessary conditions for convergence and divergence of series of positive terms. Furthermore, we present and discuss some interesting consequences and examples such as extensions of the Olivier's theorem and Raabe, Bertrand and Gauss's test.


title: CryptoEmu: An Instruction Set Emulator for Computation Over Ciphers

abstract: Fully homomorphic encryption (FHE) allows computations over encrypted data. This technique makes privacy-preserving cloud computing a reality. Users can send their encrypted sensitive data to a cloud server, get encrypted results returned and decrypt them, without worrying about data breaches. This project report presents a homomorphic instruction set emulator, CryptoEmu, that enables fully homomorphic computation over encrypted data. The software-based instruction set emulator is built upon an open-source, state-of-the-art homomorphic encryption library that supports gate-level homomorphic evaluation. The instruction set architecture supports multiple instructions that belong to the subset of ARMv8 instruction set architecture. The instruction set emulator utilizes parallel computing techniques to emulate every functional unit for minimum latency. This project report includes details on design considerations, instruction set emulator architecture, and datapath and control unit implementation. We evaluated and demonstrated the instruction set emulator's performance and scalability on a 48-core workstation. CryptoEmu has shown a significant speedup in homomorphic computation performance when compared with HELib, a state-of-the-art homomorphic encryption library.


title: Profinite groups with an automorphism of prime order whose fixed points have finite Engel sinks

abstract: A right Engel sink of an element $g$ of a group $G$ is a set ${\mathscr R}(g)$ such that for every $x\in G$ all sufficiently long commutators $[...[[g,x],x],\dots ,x]$ belong to ${\mathscr R}(g)$. (Thus, $g$ is a right Engel element precisely when we can choose ${\mathscr R}(g)=\{ 1\}$.) We prove that if a profinite group $G$ admits a coprime automorphism $\varphi $ of prime order such that every fixed point of $\varphi $ has a finite right Engel sink, then $G$ has an open locally nilpotent subgroup. A left Engel sink of an element $g$ of a group $G$ is a set ${\mathscr E}(g)$ such that for every $x\in G$ all sufficiently long commutators $[...[[x,g],g],\dots ,g]$ belong to ${\mathscr E}(g)$. (Thus, $g$ is a left Engel element precisely when we can choose ${\mathscr E}(g)=\{ 1\}$.) We prove that if a profinite group $G$ admits a coprime automorphism $\varphi $ of prime order such that every fixed point of $\varphi $ has a finite left Engel sink, then $G$ has an open pronilpotent-by-nilpotent subgroup.


title: Delay Minimization in Sliced Multi-Cell Mobile Edge Computing (MEC) Systems

abstract: We consider the problem of jointly optimizing users' offloading decisions, communication and computing resource allocation in a sliced multi-cell mobile edge computing (MEC) network. We minimize the weighted sum of the gap between the observed delay at each slice and its corresponding delay requirement, where weights set the priority of each slice. Fractional form of the objective function, discrete subchannel allocation, considered partial offloading, and the interference incorporated in the rate function, make the considered problem a complex mixed integer non-linear programming problem. Thus, we decompose the original problem into two sub-problems: (i) offloading decision-making and (ii) joint computation resource, subchannel, and power allocation. We solve the first sub-problem optimally and for the second sub-problem, leveraging on novel tools from fractional programming and Augmented Lagrangian method, we propose an efficient algorithm whose computational complexity is proved to be polynomial. Using alternating optimization, we solve these two sub-problems iteratively until convergence is obtained. Simulation results demonstrate the convergence of our proposed algorithm and its effectiveness compared to existing schemes.


title: Angular Momentum of Fission Fragments from Microscopic Theory

abstract: During nuclear fission, a heavy nucleus splits into two rotating fragments. The associated angular momentum is large, yet the mechanism of its generation and its dependence on the mass of fragments remain poorly understood. In this Letter, we provide the first microscopic calculations of angular momentum distributions in fission fragments for a wide range of fragment masses. For the benchmark case of $^{239}$Pu($n_{\text{th}}$,f), we find that the angular momentum of the fragments is largely determined by the nuclear shell structure and deformation, and that the heavy fragments therefore typically carry less angular momentum than their light partners. We use the fission model $\tt{FREYA}$ to simulate the emission of neutrons and photons from the fragments. The dependence of the angular momenta on fragment mass after the emission of neutrons and statistical photons is linear for the heavy fragments and either constant or weakly linear for the light fragments, consistent with the universal sawtooth pattern suggested by recent experimental data. Finally, we observe that using microscopic angular momentum distributions modifies the number of emitted photons significantly.


title: The $4$-rank of class groups of $K(\sqrt{n})$

abstract: Let $K/\mathbb{Q}$ be a quadratic extension. In this paper we study the $4$-rank of the class group $\text{Cl}(K(\sqrt{n}))$, where $n$ varies over squarefree rational integers. We show that for $100\%$ of squarefree $n$, the $4$-rank is given by an explicit formula involving the $2$-rank of $\text{Cl}(K)$ and the number of prime factors of $n$ which are inert in $K/\mathbb{Q}$.


title: Hierarchical Dynamic Modeling for Individualized Bayesian Forecasting

abstract: We present a case study and methodological developments in large-scale hierarchical dynamic modeling for personalized prediction in commerce. The context is supermarket sales, where improved forecasting of customer/household-specific purchasing behavior informs decisions about personalized pricing and promotions on a continuing basis. This is a big data, big modeling and forecasting setting involving many thousands of customers and items on sale, requiring sequential analysis, addressing information flows at multiple levels over time, and with heterogeneity of customer profiles and item categories. Models developed are fully Bayesian, interpretable and multi-scale, with hierarchical forms overlaid on the inherent structure of the retail setting. Customer behavior is modeled at several levels of aggregation, and information flows from aggregate to individual levels. Forecasting at an individual household level infers price sensitivity to inform personalized pricing and promotion decisions. Methodological innovations include extensions of Bayesian dynamic mixture models, their integration into multi-scale systems, and forecast evaluation with context-specific metrics. The use of simultaneous predictors from multiple hierarchical levels improves forecasts at the customer-item level of main interest. This is evidenced across many different households and items, indicating the utility of the modeling framework for this and other individualized forecasting applications.


title: Active Fire Detection in Landsat-8 Imagery: a Large-Scale Dataset and a Deep-Learning Study

abstract: Active fire detection in satellite imagery is of critical importance to the management of environmental conservation policies, supporting decision-making and law enforcement. This is a well established field, with many techniques being proposed over the years, usually based on pixel or region-level comparisons involving sensor-specific thresholds and neighborhood statistics. In this paper, we address the problem of active fire detection using deep learning techniques. In recent years, deep learning techniques have been enjoying an enormous success in many fields, but their use for active fire detection is relatively new, with open questions and demand for datasets and architectures for evaluation. This paper addresses these issues by introducing a new large-scale dataset for active fire detection, with over 150,000 image patches (more than 200 GB of data) extracted from Landsat-8 images captured around the world in August and September 2020, containing wildfires in several locations. The dataset was split in two parts, and contains 10-band spectral images with associated outputs, produced by three well known handcrafted algorithms for active fire detection in the first part, and manually annotated masks in the second part. We also present a study on how different convolutional neural network architectures can be used to approximate these handcrafted algorithms, and how models trained on automatically segmented patches can be combined to achieve better performance than the original algorithms - with the best combination having 87.2% precision and 92.4% recall on our manually annotated dataset. The proposed dataset, source codes and trained models are available on Github (https://github.com/pereira-gha/activefire), creating opportunities for further advances in the field


title: The IITM Earth System Model (IITM ESM)

abstract: Earth System Models (ESM) are important tools that allow us to understand and quantify the physical, chemical & biological mechanisms governing the rates of change of elements of the Earth System, comprising of the atmosphere, ocean, land, cryosphere and biosphere (terrestrial and marine) and related components. ESMs are essentially coupled numerical models which incorporate processes within and across the different Earth system components and are expressed as set of mathematical equations. ESMs are useful for enhancing our fundamental understanding of the climate system, its multi-scale variability, global and regional climatic phenomena and making projections of future climate change. In this chapter, we briefly describe the salient aspects of the Indian Institute of Tropical Meteorology ESM (IITM ESM), that has been developed recently at the IITM, Pune, India, for investigating long-term climate variability and change with focus on the South Asian monsoon.


title: On higher dimensional extremal varieties of general type

abstract: Relations among fundamental invariants play an important role in algebraic geometry. It is known that an $n$-dimensional variety of general type with nef canonical divisor and canonical singularities, whose image $Y$ under the canonical map is of maximal dimension, satisfies $K_X^n \ge 2 (p_g-n)$. We investigate the very interesting extremal situation $K_X^n=2(p_g-n)$, which appears in a number of geometric situations. Since these extremal varieties are natural higher dimensional analogues of Horikawa surfaces, we name them Horikawa varieties. These varieties have been previously dealt with inthe works of Fujita and Kobayashi. We carry out further studies of Horikawa varieties, proving new results on various geometric and topological issues concerning them. In particular, we prove that the geometric genus of those Horikawa varieties whose image under the canonical map is singular is bounded. We give an analogous result for polarized hyperelliptic subcanonical varieties, in particular, for polarized Calabi-Yau and Fano varieties. The pleasing numerology that emerges puts Horikawa's result on surfaces in a broader perspective. We obtain a structure theorem for Horikawa varieties and explore their pluriregularity. We use this to prove optimal results on projective normality of pluricanonical linear systems. We study the fundamental groups of Horikawa varieties, showing that they are simply connected, even if $Y$ is singular. We also prove results on deformations of Horikawa varieties, whose implications on the moduli space make them the higher dimensional analogue of curves of genus $2$.


title: Exciton-photon interactions in semiconductor nanocrystals: {radiative transitions, non-radiative processes,} and environment effects

abstract: In this review we discuss several fundamental processes taking place in semiconductor nanocrystals (quantum dots, QDs) when their electron subsystem interacts with electromagnetic (EM) radiation. The physical phenomena of light emission and EM energy transfer from a QD exciton to other electronic systems such as neighbouring nanocrystals and polarisable 3D (semi-infinite dielectric or metal) and 2D (graphene) materials are considered. The cases of direct (II-VI) and indirect (silicon) band gap semiconductors are compared. We also cover the relevant non-radiative mechanisms such as the Auger process, electron capture on dangling bonds and interaction with phonons. The emphasis is on explaining the underlying physics and illustrating it with calculated and experimental results in a comprehensive, tutorial manner.


title: Facial Biometric System for Recognition using Extended LGHP Algorithm on Raspberry Pi

abstract: In todays world, where the need for security is paramount and biometric access control systems are gaining mass acceptance due to their increased reliability, research in this area is quite relevant. Also with the advent of IOT devices and increased community support for cheap and small computers like Raspberry Pi its convenient than ever to design a complete standalone system for any purpose. This paper proposes a Facial Biometric System built on the client-server paradigm using Raspberry Pi 3 model B running a novel local descriptor based parallel algorithm. This paper also proposes an extended version of Local Gradient Hexa Pattern with improved accuracy. The proposed extended version of LGHP improved performance as shown in performance analysis. Extended LGHP shows improvement over other state-of-the-art descriptors namely LDP, LTrP, MLBP and LVP on the most challenging benchmark facial image databases, i.e. Cropped Extended Yale-B, CMU-PIE, color-FERET, LFW, and Ghallager database. Proposed system is also compared with various patents having similar system design and intent to emphasize the difference and novelty of the system proposed.


title: Conformation-dependent sequence design of polymer chains in melts

abstract: Conformation-dependent design of polymer sequences can be considered as a tool to control macromolecular self-assembly. We consider the monomer unit sequences created via the modification of polymers in a homogeneous melt in accordance with the spatial positions of the monomer units. The geometrical patterns of lamellae, hexagonally packed cylinders, and balls arranged in a body-centered cubic lattice are considered as typical microphase-separated morphologies of block copolymers. Random trajectories of polymer chains are described by the diffusion-type equations and, in parallel, simulated in the computer modeling. The probability distributions of block length $k$, which are analogous to the first-passage probabilities, are calculated analytically and determined from the computer simulations. In any domain, the probability distribution can be described by the asymptote $~k^{-3/2}$ at moderate values of $k$ if the spatial size of the block is less than the smallest characteristic size of the domain. For large blocks, the exponential asymptote $exp(-const \, k a^2/d_{as}^2)$ is valid, $d_{as}$ being the asymptotic domain length (a is the monomer unit size). The number average block lengths and their dispersities change linearly with the block length for lamellae, cylinders, and balls, when the domain is characterized by a single characteristic size. If the domain is described by more than one size, the number average block length can grow nonlinearly with the domain sizes and the length das can depend on all of them.


title: Dynamic optimal transport on networks

abstract: We study a dynamic optimal transport problem on a network. Despite the cost for transport along the edges, an additional cost, scaled with a parameter $\kappa$, has to be paid for interchanging mass between edges and vertices. We show existence f minimisers using duality and discuss the relationship of the model to other metrics such as Fisher-Rao and the classical Wasserstein metric. Finally, we examine the limiting behaviour of the model in terms of the parameter $\kappa$.


title: $L^p$-$L^q$ boundedness of $(k, a)$-Fourier multipliers with applications to Nonlinear equations

abstract: The $(k,a)$-generalised Fourier transform is the unitary operator defined using the $a$-deformed Dunkl harmonic oscillator.The main aim of this paper is to prove $L^p$-$L^q$ boundedness of $(k, a)$-generalised Fourier multipliers. To show the boundedness we first establish Paley inequality and Hausdorff-Young-Paley inequality for $(k, a)$-generalised Fourier transform. We also demonstrate applications of obtained results to study the well-posedness of nonlinear partial differential equations.


title: Analysis of an abstract mixed formulation for viscoelastic problems

abstract: This study provides an abstract framework to analyze mixed formulations in viscoelasticity, in the classic saddle point form. Standard hypothesis for mixed methods are adapted to the Volterra type equations in order to obtain stability of the proposed problem. Error estimates are derived for suitable finite element spaces. We apply the developed theory to a bending moment formulation for a linear viscoelastic Timoshenko beam and for the Laplace operator with memory terms. For both problems we report numerical results to asses the performance of the methods.


title: Deep Reinforcement Learning with Function Properties in Mean Reversion Strategies

abstract: Over the past decades, researchers have been pushing the limits of Deep Reinforcement Learning (DRL). Although DRL has attracted substantial interest from practitioners, many are blocked by having to search through a plethora of available methodologies that are seemingly alike, while others are still building RL agents from scratch based on classical theories. To address the aforementioned gaps in adopting the latest DRL methods, I am particularly interested in testing out if any of the recent technology developed by the leads in the field can be readily applied to a class of optimal trading problems. Unsurprisingly, many prominent breakthroughs in DRL are investigated and tested on strategic games: from AlphaGo to AlphaStar and at about the same time, OpenAI Five. Thus, in this writing, I want to show precisely how to use a DRL library that is initially built for games in a fundamental trading problem; mean reversion. And by introducing a framework that incorporates economically-motivated function properties, I also demonstrate, through the library, a highly-performant and convergent DRL solution to decision-making financial problems in general.


title: Training Deep Architectures Without End-to-End Backpropagation: A Survey on the Provably Optimal Methods

abstract: This tutorial paper surveys provably optimal alternatives to end-to-end backpropagation (E2EBP) -- the de facto standard for training deep architectures. Modular training refers to strictly local training without both the forward and the backward pass, i.e., dividing a deep architecture into several nonoverlapping modules and training them separately without any end-to-end operation. Between the fully global E2EBP and the strictly local modular training, there are weakly modular hybrids performing training without the backward pass only. These alternatives can match or surpass the performance of E2EBP on challenging datasets such as ImageNet, and are gaining increasing attention primarily because they offer practical advantages over E2EBP, which will be enumerated herein. In particular, they allow for greater modularity and transparency in deep learning workflows, aligning deep learning with the mainstream computer science engineering that heavily exploits modularization for scalability. Modular training has also revealed novel insights about learning and has further implications on other important research domains. Specifically, it induces natural and effective solutions to some important practical problems such as data efficiency and transferability estimation.


title: Improving ALMA's data processing effciency using a holistic approach

abstract: ALMA (Atacama Large Millimeter/submillimeter Array) is the world's largest ground-based facility for observations in the millimeter/submillimeter regime. One of ALMA's outstanding characteristics is the large effort dedicated to the quality assurance (QA) of the calibrated and imaged data products offered to the astronomical community. The Data Management Group (DMG), in charge of the data processing, review, and delivery of the ALMA data, consists of approximately 60 experts in data reduction, from the ALMA Regional Centers (ARCs) and the Joint ALMA Observatory (JAO), distributed in fourteen countries. With a throughput of more than 3,000 datasets per year, meeting the goal of delivering the pipeline-able data products within 30 days after data acquisition is a huge challenge. This paper presents (a) the history of data processing at ALMA, (b) the challenges our team had and is still facing, (c) the methodology followed to mitigate the operational risks, (d) the ongoing optimization initiatives, (e) the current data processing status, (f) the strategy which is being followed so that, in a few Cycles from now, a team of approximately 10 data reducers (DRs) at JAO can process and review some 80% of the datasets collected during an observing cycle, and, finally, (g) the important role of the ARCs for processing the remaining datasets.


title: On Some Integral Representation Of $\zeta(n)$ Involving Nielsen's Generalized Polylogarithms And The Related Partition Problem

abstract: In this paper, we study a family of single variable integral representations for some products of $\zeta(2n+1)$, where $\zeta(z)$ is Riemann zeta function and $n$ is positive integer. Such representation involves the integral $Lz(a,b):=\frac{1}{(a-1)!b!}\int_{0}^{1}\log^a (t)\log^b (1-t)dt/t$ with positive integers $a,b$, which is related to Nielsen's generalized polylogarithms. By analyzing the related partition problem, we discuss the structure of such integral representation, especially the condition of expressing products of $\zeta(2n+1)$ by finite $\mathbb{Q}(\pi)$-linear combination of $Lz(a,b)$.


title: Many-particle Quantum Hydrodynamics of Spin-1 Bose-Einstein Condensates

abstract: We develop a novel model of the magnetized spin-1 Bose-Einstein condensate (BEC) of neutral atoms, using the method of many-particle quantum hydrodynamic (QHD) and propose an original derivation of the system of continual equations. We consider bosons with a spin-spin interaction and a short range interaction in the first order in the interaction radius, on the of basis of the self-consistent field approximation of the QHD equations. We demonstrate that the dynamics of the fluid velocity and magnetization is determined by a nontrivial modification of the Euler and Landau-Lifshitz equation, and show that a nontrivial modification of the spin density evolution equation contains the spin torque effect that arises from the self-interactions between spins of the bosons. The properties of the dispersion spectrum of collective excitations are described. We obtain the new contribution of the self-interaction of spins in the spin wave spectrum together with the influence of an external magnetic field and spin-spin interactions between polarized particles. The anisotropic spin wave instability is predicted.


title: DeepFilter: an ECG baseline wander removal filter using deep learning techniques

abstract: According to the World Health Organization, around 36% of the annual deaths are associated with cardiovascular diseases and 90% of heart attacks are preventable. Electrocardiogram signal analysis in ambulatory electrocardiography, during an exercise stress test, and in resting conditions allows cardiovascular disease diagnosis. However, during the acquisition, there is a variety of noises that may damage the signal quality thereby compromising their diagnostic potential. The baseline wander is one of the most undesirable noises. In this work, we propose a novel algorithm for BLW noise filtering using deep learning techniques. The model performance was validated using the QT Database and the MIT-BIH Noise Stress Test Database from Physionet. In addition, several comparative experiments were performed against state-of-the-art methods using traditional filtering procedures as well as deep learning techniques. The proposed approach yields the best results on four similarity metrics: the sum of squared distance, maximum absolute square, percentage of root distance, and cosine similarity with 4.29 (6.35) au, 0.34 (0.25) au, 45.35 (29.69) au and, 91.46 (8.61) au, respectively. The source code of this work, containing our method and related implementations, is freely available on Github.


title: On Farkas' Lemma and Related Propositions in BISH

abstract: In this paper we analyse in the framework of constructive mathematics (BISH) the validity of Farkas' lemma and related propositions, namely the Fredholm alternative for solvability of systems of linear equations, optimality criteria in linear programming, Stiemke's lemma and the Superhedging Duality from mathematical finance, and von Neumann's minimax theorem with application to constructive game theory.


title: Parallax Systematics and Photocenter Motions of Benchmark Eclipsing Binaries in Gaia EDR3

abstract: Previous analyses of various standard candles observed by the Gaia satellite have reported statistically significant systematics in the parallaxes that have improved from $\sim$250 $\mu$as in the first data release (DR1) to 50--80 $\mu$as in the second data release (DR2). Here we examine the parallaxes newly reported in the Gaia early third data release (EDR3) using the same sample of benchmark eclipsing binaries (EBs) we used to assess the DR1 and DR2 parallaxes. We find a mean offset of $-37\pm20$ $\mu$as (Gaia$-$EB), which decreases to $-15\pm18$ $\mu$as after applying the corrections recommended by the Gaia Mission team; global systematics in the Gaia parallaxes have clearly improved and are no longer statistically significant for the EB sample, which spans $5\lesssim{G}\lesssim12$ in brightness and 0.03--3 kpc in distance. We also find that the RUWE goodness-of-fit statistic reported in Gaia EDR3 is highly sensitive to unresolved companions (tertiaries in the case of our EB sample) as well as to photocenter motion of the binaries themselves. RUWE is nearly perfectly correlated ($r^2=0.82$) with photocenter motions down to $\lesssim$0.1 mas, and surprisingly this correlation exists entirely within the nominal "good" RUWE range of 1.0--1.4. This suggests that RUWE values even slightly greater than 1.0 may signify unresolved binaries in Gaia, and that the RUWE value can serve as a quantitative predictor of the photocenter motion.


title: Thermomechanical properties of polypropylene and styrene-ethylene-butylene-styrene blends: a molecular simulation and experimental study

abstract: Polymer blends consisting of two or more polymers are important for a wide variety of industries and processes, but, the precise mechanism of their thermomechanical behaviour is incompletely understood. In order to understand clearly, it is essential to determine the miscibility and interactions between the components in polymer blend and its macroscopic thermomechanical properties. In this study, we performed experiments on SEBS and isotactic PP blends (SP) as well as molecular dynamics simulations, aiming to know the role played by molecular interactions on the thermomechanical properties. To investigate the glass transition temperature (Tg) of SEBS, PP and their blends at different ratio, the unit cell of the polymer molecular structure of each was established. The LAMMPS molecular dynamics method was used to predict the density, specific volume, free volume, enthalpy, kinetic energy, potential energy and bond energy. The (Tg) s of the SEBS, PP and SP blends were predicted by analysing these properties. Interestingly, the simulated values of the Tg of SEBS, PP and their blends showed good agreement with our experimental results obtained from dynamic mechanical analysis (DMA). This technique used in this work can be used in studying glass transition of other complex polymer blends.


title: Towards the processing, review, and delivery of 80% of the ALMA data by the Joint ALMA Observatory (JAO)

abstract: After eight observing Cycles, the Atacama Large Millimeter-submillimeter Array (ALMA) is capable of observing in eight different bands (covering a frequency range from 84 to 950 GHz), with 66 antennas and two correlators. For the current Cycle (7), ALMA offers up to 4300 hours for the 12-m array, and 3000 hours on both the 7-m of the Atacama Compact Array (ACA) and TP Array plus 750 hours in a supplemental call. From the customer perspective (i.e., the astronomical community), ALMA is an integrated product service provider, i.e. it observes in service mode, processes and delivers the data obtained. The Data Management Group (DMG) is in charge of the processing, reviewing, and delivery of the ALMA data and consists of approximately 60 experts in data reduction, from the ALMA Regional Centers (ARCs) and the Joint ALMA Observatory (JAO), distributed in fourteen countries. Prior to their delivery, the ALMA data products go through a thorough quality assurance (QA) process, so that the astronomers can work on their science without the need of significant additional calibration re-processing. Currently, around 90% of the acquired data is processed with the ALMA pipeline (the so called pipeline-able data), while the remaining 10% is processed completely manually. The Level-1 Key Performance Indicator set by the Observatory to DMG is that 90% of the pipeline-able data sets (i.e. some 80% of the data sets observed during an observing cycle) must be processed, reviewed and delivered within 30 days of data acquisition. This paper describes the methodology followed by the JAO in order to process near 80% of the total data observed during Cycle 7, a giant leap with respect to approximately 30% in Cycle 4 (October 2016 - September 2017).


title: A simple counterexample for the permanent-on-top conjecture

abstract: The permanent-on-top conjecture states that the largest eigenvalue of the Schur power matrix of a positive semi-definite Hermitian matrix H is per(H). A counterexample has been found with the help of computers, but here, I present another counterexample that can be checked by hand. My method is to use linear representations of groups to connect the spectrum of the Schur power matrix with the spectra of the entrywise(Hadamard) product matrices of permanental compound matrices. By that, we are able to study the properties of the spectrum of the Schur power matrix through the entrywise(Hadamard) product matrices of permanental compound matrices. The counterexample we find is in fact also a counterexample to a weaker conjecture related to permanental compound matrices. This conjecture was also known to be false, but the new counterexample is smaller than the known one.


title: Impact of a moon on the evolution of a planet's obliquity: a non-resonant case

abstract: We investigate how the variation of the obliquity (the axial tilt) of a hypothetical exo-Earth is effected by the presence of a satellite, an exo-Moon. Namely, we study analytically and numerically how the range of obliquity of the exo-Earth changes if an exo-Moon is added to a system comprised of exo-Sun, exo-Earth and exo-planets. We say that the impact of the exo-Moon is stabilising if upon the addition of the exo-Moon the range of obliquity decreases, while we call the impact destabilising if the range increases as the exo-Moon is added to the system. The problem is considered in a general setup. The exo-Earth is assumed to be rigid, axially symmetric and almost spherical, the difference between the largest and the smallest principal moments of inertia being a small parameter of the problem. Assuming the orbits of the celestial bodies to be quasiperiodic, we apply time averaging to study rotation of the exo-Earth at times large compared to the respective periods. Non-resonant frequencies are assumed. We identify a class of systems for which we prove analytically that the impact of the exo-Moon is stabilising and a class where it is destabilising. We also investigate numerically how the impact of the exo-Moon in a particular system comprised of a star and two planets varies on modifying the geometry of the orbits of the exo-Moon and the second planet and the initial obliquity.


title: Transitions between polarization and radicalization in a temporal bilayer echo chambers model

abstract: Echo chambers and polarisation dynamics are as of late a very prominent topic in scientific communities around the world. As these phenomena directly affect our lives and seemingly more and more as our societies and communication channels evolve it becomes ever so important for us to understand the intricacies of opinion dynamics in the modern era. Here we extend an existing echo chambers model with activity driven agents onto a bi-layer topology and study the dynamics of the polarised state as a function of interlayer couplings. Different cases of such couplings are presented - unidirectional coupling that can be reduced to a mono-layer facing an external bias, symmetric and non-symmetric couplings. We have assumed that initial conditions impose system polarisation and agent opinions are different for both layers. Such a pre-conditioned polarised state can sustain without explicit homophilic interactions provided the coupling strength between agents belonging to different layers is weak enough. For a strong unidirectional or attractive coupling between two layers a discontinuous transition to a radicalised state takes place when mean opinions in both layers are the same. When coupling constants between the layers are of different signs the system exhibits sustained or decaying oscillations. Transitions between these states are analysed using a mean field approximation and classified in the framework of bifurcation theory.


title: Are We There Yet? Learning to Localize in Embodied Instruction Following

abstract: Embodied instruction following is a challenging problem requiring an agent to infer a sequence of primitive actions to achieve a goal environment state from complex language and visual inputs. Action Learning From Realistic Environments and Directives (ALFRED) is a recently proposed benchmark for this problem consisting of step-by-step natural language instructions to achieve subgoals which compose to an ultimate high-level goal. Key challenges for this task include localizing target locations and navigating to them through visual inputs, and grounding language instructions to visual appearance of objects. To address these challenges, in this study, we augment the agent's field of view during navigation subgoals with multiple viewing angles, and train the agent to predict its relative spatial relation to the target location at each timestep. We also improve language grounding by introducing a pre-trained object detection module to the model pipeline. Empirical studies show that our approach exceeds the baseline model performance.


title: Regularising Transformations for Complex Differential Equations with Movable Algebraic Singularities

abstract: In a 1979 paper, K. Okamoto introduced the space of initial values for the six Painlev\'e equations and their associated Hamiltonian systems, showing that these define regular initial value problems at every point of an augmented phase space, a rational surface with certain exceptional divisors removed. We show that the construction of the space of initial values remains meaningful for certain classes of second-order complex differential equations, and more generally, Hamiltonian systems, where all movable singularities of all their solutions are algebraic poles (by some authors denoted the quasi-Painlev\'e property), which is a generalisation of the Painlev\'e property. The difference here is that the initial value problems obtained in the extended phase space become regular only after an additional change of dependent and independent variables. Constructing the analogue of space of initial values for these equations in this way also serves as an algorithm to single out, from a given class of equations or system of equations, those equations which are free from movable logarithmic branch points.


title: TIC 168789840: A Sextuply-Eclipsing Sextuple Star System

abstract: We report the discovery of a sextuply-eclipsing sextuple star system from TESS data, TIC 168789840, also known as TYC 7037-89-1, the first known sextuple system consisting of three eclipsing binaries. The target was observed in Sectors 4 and 5 during Cycle 1, with lightcurves extracted from TESS Full Frame Image data. It was also previously observed by the WASP survey and ASAS-SN. The system consists of three gravitationally-bound eclipsing binaries in a hierarchical structure of an inner quadruple system with an outer binary subsystem. Follow-up observations from several different observatories were conducted as a means of determining additional parameters. The system was resolved by speckle interferometry with a 0."42 separation between the inner quadruple and outer binary, inferring an estimated outer period of ~2 kyr. It was determined that the fainter of the two resolved components is an 8.217 day eclipsing binary, which orbits the inner quadruple that contains two eclipsing binaries with periods of 1.570 days and 1.306 days. MCMC analysis of the stellar parameters has shown that the three binaries of TIC 168789840 are "triplets", as each binary is quite similar to the others in terms of mass, radius, and Teff. As a consequence of its rare composition, structure, and orientation, this object can provide important new insight into the formation, dynamics, and evolution of multiple star systems. Future observations could reveal if the intermediate and outer orbital planes are all aligned with the planes of the three inner eclipsing binaries.


title: Noetherian spaces in choiceless set theory

abstract: I prove several independence results in the choiceless ZF+DC theory which separate algebraic and non-algebraic consequences of the axiom of choice.


title: The maximization of the p-Laplacian energy for a two-phase material

abstract: We consider the optimal arrangement of two diffusion materials in a bounded open set $\Omega\subset \mathbb{R}^N$ in order to maximize the energy. The diffusion problem is modeled by the $p$-Laplacian operator. It is well known that this type of problems has no solution in general and then that it is necessary to work with a relaxed formulation. In the present paper we obtain such relaxed formulation using the homogenization theory, i.e. we replace both materials by microscopic mixtures of them. Then we get some uniqueness results and a system of optimality conditions. As a consequence we prove some regularity properties for the optimal solutions of the relaxed problem. Namely, we show that the flux is in the Sobolev space $H^1(\Omega)^N$ and that the optimal proportion of the materials is derivable in the orthogonal direction to the flux. This will imply that the unrelaxed problem has no solution in general. Our results extend those obtained by the first author for the Laplace operator.


title: The structure of global conservation laws in Galerkin plasma models

abstract: Plasmas are highly nonlinear and multi-scale, motivating a hierarchy of models to understand and describe their behavior. However, there is a scarcity of plasma models of lower fidelity than magnetohydrodynamics (MHD). Galerkin models, obtained by projection of the MHD equations onto a truncated modal basis, can furnish this gap in the lower levels of the model hierarchy. In the present work, we develop low-dimensional Galerkin plasma models which preserve global conservation laws by construction. This additional model structure enables physics-constrained machine learning algorithms that can discover these types of low-dimensional plasma models directly from data. This formulation relies on an energy-based inner product which takes into account all of the dynamic variables. The theoretical results here build a bridge to the extensive Galerkin literature in fluid mechanics, and facilitate the development of physics-constrained reduced-order models from plasma data.


title: Coloring redundant algebraic hypergraphs

abstract: We prove several consistency results in choiceless set theory ZF+DC regarding countable chromatic numbers of various algebraic hypergraphs on Euclidean spaces.


title: Quantum Generative Models for Small Molecule Drug Discovery

abstract: Existing drug discovery pipelines take 5-10 years and cost billions of dollars. Computational approaches aim to sample from regions of the whole molecular and solid-state compounds called chemical space which could be on the order of 1060 . Deep generative models can model the underlying probability distribution of both the physical structures and property of drugs and relate them nonlinearly. By exploiting patterns in massive datasets, these models can distill salient features that characterize the molecules. Generative Adversarial Networks (GANs) discover drug candidates by generating molecular structures that obey chemical and physical properties and show affinity towards binding with the receptor for a target disease. However, classical GANs cannot explore certain regions of the chemical space and suffer from curse-of-dimensionality. A full quantum GAN may require more than 90 qubits even to generate QM9-like small molecules. We propose a qubit-efficient quantum GAN with a hybrid generator (QGAN-HG) to learn richer representation of molecules via searching exponentially large chemical space with few qubits more efficiently than classical GAN. The QGANHG model is composed of a hybrid quantum generator that supports various number of qubits and quantum circuit layers, and, a classical discriminator. QGAN-HG with only 14.93% retained parameters can learn molecular distribution as efficiently as classical counterpart. The QGAN-HG variation with patched circuits considerably accelerates our standard QGANHG training process and avoids potential gradient vanishing issue of deep neural networks. Code is available on GitHub https://github.com/jundeli/quantum-gan.


title: Study of Nonlinear Interaction between Waves and Ocean Currents Using High-Fidelity Simulation and Machine Learning

abstract: Modeling ocean surface waves under complex ocean current conditions is of crucial importance to many naval applications. For example, traveling ships and underwater vehicles generate spatially heterogeneous currents behind them through their drag and propeller motions. The strong currents can influence the surface wave pattern in the ship wake. In this study, the nonlinear interactions between waves and complex wake currents are investigated using numerical simulations. An in-house code is developed for high-fidelity simulations of a nonlinear phase-resolved ocean wavefield interacting with subsurface currents. Several typical wake patterns are simulated using the present numerical method, and the influence of complex currents on the waves is analyzed quantitatively using theoretical solutions of wave-current interactions. We also present a method for solving the inverse problem of deducing the current field based on surface-wave data using machine-learning techniques. A deep neural network is designed for processing spatial-temporal surface wave data. Detailed analyses on the distributions of regression errors and the training dataset-dependency show that the proposed neural network can effectively deduce the current field.


title: A Hardy-Ramanujan type inequality for shifted primes and sifted sets

abstract: We establish an analog of the Hardy-Ramanujan inequality for counting members of sifted sets with a given number of distinct prime factors. In particular, we establish a bound for the number of shifted primes p+a below x with k distinct prime factors, uniformly for all positive integers k.


title: Rate Allocation and Content Placement in Cache Networks

abstract: We introduce the problem of optimal congestion control in cache networks, whereby \emph{both} rate allocations and content placements are optimized \emph{jointly}. We formulate this as a maximization problem with non-convex constraints, and propose solving this problem via (a) a Lagrangian barrier algorithm and (b) a convex relaxation. We prove different optimality guarantees for each of these two algorithms; our proofs exploit the fact that the non-convex constraints of our problem involve DR-submodular functions.


title: Transience of symmetric non-local Dirichlet forms

abstract: We establish transience criteria for symmetric non-local Dirichlet forms on $L^2({\mathbb R}^d)$ in terms of the coefficient growth rates at infinity. Applying these criteria, we find a necessary and sufficient condition for recurrence of Dirichlet forms of symmetric stable-like with unbounded/degenerate coefficients. This condition indicates that both of the coefficient growth rates of small and big jump parts affect the sample path properties of the associated symmetric jump processes.


title: Complexity for link complement States in Chern Simons Theory

abstract: We study notions of complexity for link complement states in Chern Simons theory with compact gauge group $G$. Such states are obtained by the Euclidean path integral on the complement of $n$-component links inside a 3-manifold $M_3$. For the Abelian theory at level $k$ we find that a natural set of fundamental gates exists and one can identify the complexity as differences of linking numbers modulo $k$. Such linking numbers can be viewed as coordinates which embeds all link complement states into $\mathbb{Z}_k ^{\otimes n(n-1)/2}$ and the complexity is identified as the distance with respect to a particular norm. For non-Abelian Chern Simons theories, the situation is much more complicated. We focus here on torus link states and show that the problem can be reduced to defining complexity for a single knot complement state. We suggest a systematic way to choose a set of minimal universal generators for single knot complement states and then evaluate the complexity using such generators. A detailed illustration is shown for $SU(2)_k$ Chern Simons theory and the results can be extended to general compact gauge group.


title: Initial operation and data processing on a system for real-time evaluation of Thomson scattering signals on the Large Helical Device

abstract: A scalable system for real-time analysis of electron temperature and density based on signals from the Thomson scattering diagnostic, initially developed for and installed on the NSTX-U experiment, was recently adapted for the Large Helical Device (LHD) and operated for the first time during plasma discharges. During its initial operation run, it routinely recorded and processed signals for four spatial points at the laser repetition rate of 30 Hz, well within the system's rated capability for 60 Hz. We present examples of data collected from this initial run and describe subsequent adaptations to the analysis code to improve the fidelity of the temperature calculations.


title: The Pauli Exclusion Operator: example of Hooke's atom

abstract: The Pauli Exclusion Operator (PEO) which ensures proper symmetry of the eigenstates of multi-electron systems with respect to exchange of each pair of electrons is introduced. Once PEO is added to the Hamiltonian, no additional constraints on multi-electron wave function due to the Pauli exclusion principle are needed. For two-electron states in two dimensions ($2D$) the PEO can be expressed in a closed form in terms of momentum operators, while in the position representation PEO is a non-local operator. Generalizations of PEO for multi-electron systems is introduced. Several approximations to PEO are discussed. Examples of analytical and numerical calculations of PEO are given for isotropic and anisotropic Hooke's atom in~$2D$. Application of approximate and kernel forms of PEO for calculations of energies and states in~$2D$ Hooke's atom are analyzed. Relation of PEO to standard variational calculations with the use of Slater determinant is discussed.


title: The shifted ODE method for underdamped Langevin MCMC

abstract: In this paper, we consider the underdamped Langevin diffusion (ULD) and propose a numerical approximation using its associated ordinary differential equation (ODE). When used as a Markov Chain Monte Carlo (MCMC) algorithm, we show that the ODE approximation achieves a $2$-Wasserstein error of $\varepsilon$ in $\mathcal{O}\big(d^{\frac{1}{3}}/\varepsilon^{\frac{2}{3}}\big)$ steps under the standard smoothness and strong convexity assumptions on the target distribution. This matches the complexity of the randomized midpoint method proposed by Shen and Lee [NeurIPS 2019] which was shown to be order optimal by Cao, Lu and Wang. However, the main feature of the proposed numerical method is that it can utilize additional smoothness of the target log-density $f$. More concretely, we show that the ODE approximation achieves a $2$-Wasserstein error of $\varepsilon$ in $\mathcal{O}\big(d^{\frac{2}{5}}/\varepsilon^{\frac{2}{5}}\big)$ and $\mathcal{O}\big(\sqrt{d}/\varepsilon^{\frac{1}{3}}\big)$ steps when Lipschitz continuity is assumed for the Hessian and third derivative of $f$. By discretizing this ODE using a third order Runge-Kutta method, we can obtain a practical MCMC method that uses just two additional gradient evaluations per step. In our experiment, where the target comes from a logistic regression, this method shows faster convergence compared to other unadjusted Langevin MCMC algorithms.


title: An intuitive parametric model for 3D compressible hydrodynamical and MHD turbulence

abstract: An analytical model for three-dimensional incompressible turbulence was recently introduced in the hydrodynamics community which, with only a few parameters, shares many properties of experimental and numerical turbulence, notably intermittency (non-Gaussianity), the energy cascade (skewness), and vorticity alignment properties. In view of modeling astrophysical environments, we introduce a manner to extend to compressible fluids the three-dimensional turbulent velocity field model of Chevillard et al. (2010), as well as the three 3D turbulent magnetic field models of Durrive et al. (2020), following the same procedure based on the concept of multiplicative chaos. Our model provides a complementary tool to numerical simulations, as it enables us to generate very quickly fairly realistic velocity fields and magnetic fields, the statistics of which are controllable with intuitive parameters. Therefore our model will also provide a useful tool for observers in astrophysics. Finally, maybe even more than the model itself, it is the very procedure that matters the most: our model is modular, in the sense that it is constructed gradually, with intuitive and physically motivated steps, so that it is prone to many further improvements.

categories: astro-ph.GA astro-ph.SR physics.flu-dyn physics.plasm-ph


title: A Thermodynamic Core using Voltage-Controlled Spin-Orbit-Torque Magnetic Tunnel Junctions

abstract: We present a magnetic implementation of a thermodynamic computing fabric. Magnetic devices within computing cores harness thermodynamics through its voltage-controlled thermal stability; while the evolution of network states is guided by the spin-orbit-torque effect. We theoretically derive the dynamics of the cores and show that the computing fabric can successfully compute ground states of a Boltzmann Machine. Subsequently, we demonstrate the physical realization of these devices based on a CoFeB-MgO magnetic tunnel junction structure. The results of this work pave the path towards the realization of highly efficient, high-performance thermodynamic computing hardware. Finally, this paper will also give a perspective of computing beyond thermodynamic computing.


title: Tuning Dzyaloshinskii-Moriya Interaction in Ferrimagnetic GdCo: A First Principles Approach

abstract: We present a systematic analysis of our ability to tune chiral Dzyaloshinskii-Moriya Interactions (DMI) in compensated ferrimagnetic Pt/GdCo/Pt1-xWx trilayers by cap layer composition. Using first principles calculations, we show that the DMI increases rapidly for only ~ 10% W and saturates thereafter, in agreement with experiments. The calculated DMI shows a spread in values around the experimental mean, depending on the atomic configuration of the cap layer interface. The saturation is attributed to the vanishing of spin orbit coupling energy at the cap layer and the simultaneous constancy at the bottom interface. Additionally, we predict the DMI in Pt/GdCo/X (X=Ta, W, Ir) and find that W in the cap layer favors a higher DMI than Ta and Ir that can be attributed to the difference in d-band alignment around the Fermi level. Our results open up exciting combinatorial possibilities for controlling the DMI in ferrimagnets towards nucleating and manipulating ultrasmall high-speed skyrmions.


title: Asymptotic interplay of states and adapted coupling gains in the Lohe hermitian sphere model

abstract: We study emergent dynamics of the Lohe hermitian sphere (LHS) model with the same free flows under the dynamic interplay between state evolution and adaptive couplings. The LHS model is a complex counterpart of the Lohe sphere (LS) model on the unit sphere in Euclidean space, and when particles lie in the Euclidean unit sphere embedded in $\bbc^{d+1}$, it reduces to the Lohe sphere model. In the absence of interactions between states and coupling gains, emergent dynamics have been addressed in [22]. In this paper, we further extend earlier results in the aforementioned work to the setting in which the state and coupling gains are dynamically interrelated via two types of coupling laws, namely anti-Hebbian and Hebbian coupling laws. In each case, we present two sufficient frameworks leading to complete aggregation depending on the coupling laws, when the corresponding free flow is the same for all particles.


title: Emergent behaviors of homogeneous Lohe Hermitian sphere particles under time-delayed interactions

abstract: We study emergent behaviors of the Lohe hermitian sphere(LHS) model with a time-delay for a homogeneous ensemble. The LHS model is a complex counterpart of the Lohe sphere(LS) aggregation model on the unit sphere in Euclidean space, and it describes the aggregation of particles on the unit hermitian sphere in $\mathbb{C}^d$ with $d \geq 2$, Recently it has been introduced by two authors of this work as a special case of the Lohe tensor model [23]. When the coupling gain pair satisfies a specific linear relation, namely the Stuart-Landau(SL) coupling gain pair, it can be embedded into the LS model on $\mathbb{R}^{2d}$. In this work, we show that if the coupling gain pair is close to the SL coupling pair case, the dynamics of the LHS model exhibits an emergent aggregate phenomenon via the interplay between time-delayed interactions and nonlinear coupling between states. For this, we present several frameworks for complete aggregation and practical aggregation in terms of initial data and system parameters using the Lyapunov functional approach.


title: Tail inequalities for restricted classes of discrete random variables

abstract: Let $X$ be an integrable discrete random variable over $\{0, 1, 2, \ldots\}$ with $\mathbb{P}(X = i + 1) \leq \mathbb{P}(X = i)$ for all $i$. Then for any integer $a \geq 1$, $\mathbb{P}(X \leq a) \leq \mathbb{E}[X] / (2a - 1)$. Let $W$ be an discrete random variable over $\{\ldots, -2, -1, 0, 1, 2, \ldots\}$ with finite second moment where the $\mathbb{P}(W = i)$ values are unimodal. Then $\mathbb{P}(|W - \mathbb{E}[W]| \geq a) \leq (\mathbb{V}(W) + 1 / 12) / (2(a - 1 / 2)^2)$.


title: BERT & Family Eat Word Salad: Experiments with Text Understanding

abstract: In this paper, we study the response of large models from the BERT family to incoherent inputs that should confuse any model that claims to understand natural language. We define simple heuristics to construct such examples. Our experiments show that state-of-the-art models consistently fail to recognize them as ill-formed, and instead produce high confidence predictions on them. As a consequence of this phenomenon, models trained on sentences with randomly permuted word order perform close to state-of-the-art models. To alleviate these issues, we show that if models are explicitly trained to recognize invalid inputs, they can be robust to such attacks without a drop in performance.


title: Visualizing adverse events in clinical trials using correspondence analysis with R-package visae

abstract: We propose to apply stacked CA using contribution biplots as a tool to explore differences in AE data among treatments in clinical trials. We defined five levels of refinement for the analysis based on data derived from the Common Terminology Criteria for Adverse Events (CTCAE) grades, domains, terms and their combinations. In addition, we developed a Shiny app built in an R-package, publicly available on Comprehensive R Archive Network (CRAN), to interactively investigate CA configurations. Data from two randomized controlled trials (RCT) were used to illustrate the proposed methods: NSABP R-04, a neoadjuvant rectal 2x2 factorial trial comparing radiation therapy with either capecitabine (Cape) or 5-fluorouracil (5-FU) alone with or without oxaliplatin (Oxa), and NSABP B-35, a double-blind RCT comparing tamoxifen to anastrozole in postmenopausal women with hormone-positive ductal carcinoma in situ. In the R04 trial (n=1308), CA biplots displayed the discrepancies between single agent treatments and their combinations with Oxa at all levels of AE classes, such that these discrepancies were responsible for the largest portion of the explained variability among treatments. In addition, an interaction effect when adding Oxa to Cape/5-FU was identified when the distance between Cape+Oxa and 5-FU+Oxa was observed to be larger than the distance between 5-FU and Cape, with Cape+Oxa and 5-FU+Oxa in different quadrants of the CA biplots. In the B35 trial (n=3009), CA biplots showed different patterns for non-adherent Anastrozole and Tamoxifen compared with their adherent counterparts. CA with contribution biplot is an effective tool that can be used to summarize AE data in a two-dimensional display while minimizing the loss of information and interpretation.


title: Approximately Strategyproof Tournament Rules in the Probabilistic Setting

abstract: We consider the manipulability of tournament rules which map the results of $\binom{n}{2}$ pairwise matches and select a winner. Prior work designs simple tournament rules such that no pair of teams can manipulate the outcome of their match to improve their probability of winning by more than $1/3$, and this is the best possible among any Condorcet-consistent tournament rule (which selects an undefeated team whenever one exists) [Schneider et al., 2017, Schvartzman et al., 2020]. These lower bounds require the manipulators to know precisely the outcome of all future matches. We take a beyond worst-case view and instead consider tournaments which are "close to uniform": the outcome of all matches are independent, and no team is believed to win any match with probability exceeding $1/2+\varepsilon$. We show that Randomized Single Elimination Bracket [Schneider et al., 2017] and a new tournament rule we term Randomized Death Match have the property that no pair of teams can manipulate the outcome of their match to improve their probability of winning by more than $\varepsilon/3 + 2\varepsilon^2/3$, for all $\varepsilon$, and this is the best possible among any Condorcet-consistent tournament rule. Our main technical contribution is a recursive framework to analyze the manipulability of certain forms of tournament rules. In addition to our main results, this view helps streamline previous analysis of Randomized Single Elimination Bracket, and may be of independent interest.


title: Implementation of Polygonal Mesh Refinement in MATLAB

abstract: We present a simple and efficient MATLAB implementation of the local refinement for polygonal meshes. The purpose of this implementation is primarily educational, especially on the teaching of adaptive virtual element methods.


title: Deep Ensemble Learning-based Approach to Real-time Power System State Estimation

abstract: Power system state estimation (PSSE) is commonly formulated as weighted least-square (WLS) algorithm and solved using iterative methods such as Gauss-Newton methods. However, iterative methods have become more sensitive to system operating conditions than ever before due to the deployment of intermittent renewable energy sources, low carbon technologies (e.g., electric vehicles), and demand response programs. Appropriate PSSE approaches are required to avoid pitfalls of the WLS-based PSSE computations for accurate prediction of operating conditions. This paper proposes a data-driven real-time PSSE using a deep ensemble learning algorithm. In the proposed approach, the ensemble learning setup is formulated with dense residual neural networks as base-learners and multivariate-linear regressor as meta-learner. Historical measurements and states are utilised to train and test the model. The trained model can be used in real-time to estimate power system states (voltage magnitudes and phase angles) using real-time measurements. Most of current data-driven PSSE methods assume the availability of a complete set of measurements, which may not be the case in real power system data-acquisition. This paper adopts multivariate linear regression to forecast system states for instants of missing measurements to assist the proposed PSSE technique. Case studies are performed on various IEEE standard benchmark systems to validate the proposed approach. The results show that the proposed approach outperforms existing data-driven PSSE methods techniques.


title: Unraveling the Dynamic Importance of County-level Features in Trajectory of COVID-19

abstract: The objective of this study was to investigate the importance of multiple county-level features in the trajectory of COVID-19. We examined feature importance across 2,787 counties in the United States using a data-driven machine learning model. We trained random forest models using 23 features representing six key influencing factors affecting pandemic spread: social demographics of counties, population activities, mobility within the counties, movement across counties, disease attributes, and social network structure. Also, we categorized counties into multiple groups according to their population densities, and we divided the trajectory of COVID-19 into three stages: the outbreak stage, the social distancing stage, and the reopening stage. The study aims to answer two research questions: (1) The extent to which the importance of heterogeneous features evolves in different stages; (2) The extent to which the importance of heterogeneous features varies across counties with different characteristics. We fitted a set of random forest models to determine weekly feature importance. The results showed that: (1) Social demographic features, such as gross domestic product, population density, and minority status maintained high-importance features throughout stages of COVID-19 across the 2787 studied counties; (2) Within-county mobility features had the highest importance in county clusters with higher population densities; (3) The feature reflecting the social network structure (Facebook, social connectedness index), had higher importance in the models for counties with higher population densities. The results show that the data-driven machine learning models could provide important insights to inform policymakers regarding feature importance for counties with various population densities and in different stages of a pandemic life cycle.


title: Prediction of Stable Ground-State Binary Sodium-Potassium Interalkalis under High Pressures

abstract: The complex structures and electronic properties of alkali metals and their alloys provide a natural laboratory for studying the interelectronic interactions of metals under compression. A recent theoretical study (J. Phys. Chem. Lett. 2019, 10, 3006) predicted an interesting pressure-induced decomposition-recombination behavior of the Na2K compound over a pressure range of 10 - 500 GPa. However, a subsequent experiment (Phys. Rev. B 2020, 101, 224108) reported the formation of NaK rather than Na2K at pressures above 5.9 GPa. To address this discordance, we study the chemical stability of different stoichiometries of NaxK (x = 1/4, 1/3, 1/2, 2/3, 3/4, 4/3, 3/2 and 1 - 4) by effective structure searching method combined with first-principles calculations. Na2K is calculated to be unstable at 5 - 35 GPa due to the decomposition reaction Na2K-> NaK + Na, coinciding well with the experiment. NaK undergoes a combination-decomposition-recombination process accompanied by an opposite charge-transfer behavior between Na and K with pressure. Besides NaK, two hitherto unknown compounds NaK3 and Na3K2 are uncovered. NaK3 is a typical metallic alloy, while Na3K2 is an electride with strong interstitial electron localization.


title: Quantum Random Number Generation with Uncharacterized Laser and Sunlight

abstract: The entropy or randomness source is an essential ingredient in random number generation. Quantum random number generators generally require well modeled and calibrated light sources, such as a laser, to generate randomness. With uncharacterized light sources, such as sunlight or an uncharacterized laser, genuine randomness is practically hard to be quantified or extracted owing to its unknown or complicated structure. By exploiting a recently proposed source-independent randomness generation protocol, we theoretically modify it by considering practical issues and experimentally realize the modified scheme with an uncharacterized laser and a sunlight source. The extracted randomness is guaranteed to be secure independent of its source and the randomness generation speed reaches 1 Mbps, three orders of magnitude higher than the original realization. Our result signifies the power of quantum technology in randomness generation and paves the way to high-speed semi-self-testing quantum random number generators with practical light sources.


title: Torsors in super-symmetry

abstract: Torsors under affine groups are generalized in the super context by super-torsors under affine super-groups. We investigate those super-torsors by using Hopf-algebra language and techniques. It is explicitly shown, under suitable assumptions, that every super-torsor arises from an ordinary torsor. Especially, the objects with affinity restriction, or namely, the affine super-torsors and the affine ordinary torsors are shown to be precisely in one-to-one correspondence. The results play substantial roles in ongoing construction of super-symmetric Picard-Vessiot theory.


title: Braided Brin-Thompson Groups

abstract: We construct braided versions $sV_{br}$ of the Brin-Thompson groups $sV$ and prove that they are of type $F_\infty$. The proof involves showing that the matching complexes of colored arcs on surfaces are highly connected.