title: Findings of the 2nd Photonics and Electronics Technology for Extreme Scale Computing Workgroup (REPETE): Design Challenges for Socket Level Photonic I/O

abstract: To inform research activities in HPC interconnect of strategic importance to the USG beyond 2018, in January of 2018, the DoD sponsored the 2nd Photonics and Electronics Technology for Extreme-scale Computing (REPETE) workgroup. REPETE investigated new challenges in the area of HPC interconnect inspired by technology trends and challenges of vital importance to USG stakeholders. The REPETE Working group investigated two focus areas of interest to the USG: Socket Level Photonic IO and Cryogenic Photonic IO. The working group spanned industry, academia, and government, in research, development, product, and technology investment areas. The workgroup team focusing on current and future design challenges for socket level photonic IO began discussing technical challenges and current state of applying photonics to off-chip IO in April of 2019 through biweekly meetings that concluded in early September of 2019. The focus of these meetings was to discuss what technology exists for moving the conversion of electrical signaling to photonic signaling from the node (as is currently done) down to within the socket for off-chip IO. Areas discussed include: I/O Requirements and Trends at the Compute Socket for 2025 and Beyond; Current and Near-Term Copper Solutions for Off-Chip, Socket-Level Interconnect; Current Photonic Solutions for Off-Chip or Socket-Level Interconnect; Light Generation for Off-Chip, Socket-Level IO; Fabrication and Packaging of Photonic Integrated Circuits; All-Photonic Switching Technology; and Simulation of Photonic Interconnects.


title: Distributed Simulation and Visualization of The ALICE Detector Magnetic Field

abstract: The ALICE detector at CERN uses properties of the magnetic field acting on charged particles as part of the particle tracking and identification system -- via measuring the strength of bending of charged particles in a magnetic field supplied by electromagnets. During the calibration of the detector the magnetic field was measured by the scientific team. The measurement points were then fitted using Chebyshev polynomials to create a field model. That model is used extensively in particle trajectory reconstruction, but for maximum compatibility it is handled by CPU-only software. We propose multiple approaches to the implementation of the model in a shader language, which allows it to run in a graphical GPU environment. This makes the model run more efficiently, as well as enables it to be used for field visualisation, which was not previously done for ALICE.


title: Spherically symmetric 't Hooft-Polyakov monopoles

abstract: A general analytic spherically symmetric solution of the Bogomol'nyi equations is found. It depends on two constants and one arbitrary function on radius and contains the Bogomol'nyi-Prasad-Sommerfield and Singleton solutions as particular cases. Thus all spherically symmetric 't Hooft-Polyakov monopoles with massless scalar field and minimal energy are derived.


title: A simple model of the strange-metal phase in cuprates

abstract: Hole doping of superconducting cuprates generates lattice defects of O atoms. At and beyond closing of the pseudogap, p =< p*, they form a highly symmetric superlattice. Umklapp processes involving reciprocal lattice vectors associated with both the host lattice and the O superlattice could account for the linear temperature dependence of resistivity.


title: Game Redesign in No-regret Game Playing

abstract: We study the game redesign problem in which an external designer has the ability to change the payoff function in each round, but incurs a design cost for deviating from the original game. The players apply no-regret learning algorithms to repeatedly play the changed games with limited feedback. The goals of the designer are to (i) incentivize all players to take a specific target action profile frequently; and (ii) incur small cumulative design cost. We present game redesign algorithms with the guarantee that the target action profile is played in T-o(T) rounds while incurring only o(T) cumulative design cost. Game redesign describes both positive and negative applications: a benevolent designer who incentivizes players to take a target action profile with better social welfare compared to the solution of the original game, or a malicious attacker whose target action profile benefits themselves but not the players. Simulations on four classic games confirm the effectiveness of our proposed redesign algorithms.


title: WannSymm: A symmetry analysis code for Wannier orbitals

abstract: We derived explicit expressions of symmetry operators on Wannier basis, and implemented these operators in WannSymm software. Based on this implementation, WannSymm can i) symmetrize the real-space Hamiltonian output from Wannier90 code, ii) generate symmetry operators of the little group at a specific k-point, and iii) perform symmetry analysis for Wannier band structure. In general, symmetrized Hamiltonians yield improved results compared with the original ones when they are employed for nodal structure searching, surface Green's function calculations, and other model calculations.


title: MOTENS: A Pedagogical Design Model for Serious Cyber Games

abstract: In the last few years, serious games have become popular, with a consensus of the benefits for teaching cyber security awareness and education. However, there is still a lack of pedagogical driven methodologies and tools to support serious games design to ensure they achieve the learning objectives. This paper proposes MOTENS, a pedagogical model, to design serious cyber games based on the gaps we identified in the current games design models and the lessons learnt from creating a serious tabletop game called Riskio, designed to teach cyber security awareness and education. The MOTENS model has six high-level components. Five components are linked to the games/design mechanics, and one component, `Theory', that supports the design's cognitive principles, including players' motivation. The model is used to design serious cyber games and goes through five stages, from identifying and segmenting target players, steps to creating game mechanics linked to pedagogy instruction and then to testing to create a serious game that is designed to achieve the games learning objectives.


title: Semantic Identifiers and DNS Names for IoT

abstract: In this paper, we propose a scheme for representing semantic metadata of IoT devices in compact identifiers and DNS names to enable simple discovery and search with standard DNS servers. Our scheme defines a binary identifier as a sequence of bits: a Context to use and several bits of fields corresponding to semantic properties specific to the Context. The bit string is then encoded as base32 characters and registered in DNS. Furthermore, we use the compact semantic DNS names to offer support for search and discovery. We propose to take advantage of the DNS system as the basic functionality for querying and discovery of semantic properties related to IoT devices. We have defined three specific Contexts for hierarchical semantic properties as well as logical and geographical locations. For this last part, we have developed two prototypes for managing geo-identifiers in LoRa networks, one based on Node and the Redis in-memory database, the other one based on the CoreDNS server.


title: Exploiting Cross-Modal Prediction and Relation Consistency for Semi-Supervised Image Captioning

abstract: The task of image captioning aims to generate captions directly from images via the automatically learned cross-modal generator. To build a well-performing generator, existing approaches usually need a large number of described images, which requires a huge effects on manual labeling. However, in real-world applications, a more general scenario is that we only have limited amount of described images and a large number of undescribed images. Therefore, a resulting challenge is how to effectively combine the undescribed images into the learning of cross-modal generator. To solve this problem, we propose a novel image captioning method by exploiting the Cross-modal Prediction and Relation Consistency (CPRC), which aims to utilize the raw image input to constrain the generated sentence in the commonly semantic space. In detail, considering that the heterogeneous gap between modalities always leads to the supervision difficulty of using the global embedding directly, CPRC turns to transform both the raw image and corresponding generated sentence into the shared semantic space, and measure the generated sentence from two aspects: 1) Prediction consistency. CPRC utilizes the prediction of raw image as soft label to distill useful supervision for the generated sentence, rather than employing the traditional pseudo labeling; 2) Relation consistency. CPRC develops a novel relation consistency between augmented images and corresponding generated sentences to retain the important relational knowledge. In result, CPRC supervises the generated sentence from both the informativeness and representativeness perspectives, and can reasonably use the undescribed images to learn a more effective generator under the semi-supervised scenario.


title: Periodic cosmic evolution in $f(Q)$ gravity formalism

abstract: We study the periodic cosmic transit behavior of accelerated universe in the framework of symmetric teleparallelism. The exact solution of field equations is obtained by employing a well known deceleration parameter (DP) called periodic varying deceleration parameter (PVDP), $q = m cos kt-1$. The viability and physical reliability of the DP are studied by using the observational constraint. The dynamics of periodicity and singularity are addressed in details with respect to time and redshift parameter. Several ECs are discussed in this setting.


title: Clustering of Bank Customers using LSTM-based encoder-decoder and Dynamic Time Warping

abstract: Clustering is an unsupervised data mining technique that can be employed to segment customers. The efficient clustering of customers enables banks to design and make offers based on the features of the target customers. The present study uses a real-world financial dataset (Berka, 2000) to cluster bank customers by an encoder-decoder network and the dynamic time warping (DTW) method. The customer features required for clustering are obtained in four ways: Dynamic Time Warping (DTW), Recency Frequency and Monetary (RFM), LSTM encoder-decoder network, and our proposed hybrid method. Once the LSTM model was trained by customer transaction data, a feature vector of each customer was automatically extracted by the encoder.Moreover, the distance between pairs of sequences of transaction amounts was obtained using DTW. Another vector feature was calculated for customers by RFM scoring. In the hybrid method, the feature vectors are combined from the encoder-decoder output, the DTW distance, and the demographic data (e.g., age and gender). Finally, feature vectors were introduced as input to the k-means clustering algorithm, and we compared clustering results with Silhouette and Davies-Bouldin index. As a result, the clusters obtained from the hybrid approach are more accurate and meaningful than those derived from individual clustering techniques. In addition, the type of neural network layers had a substantial effect on the clusters, and high network error does not necessarily worsen clustering performance.


title: Holographic application in cosmology: Thermodynamics of the Van der Waals cosmic fluid

abstract: This paper is devoted to investigate the thermodynamic stability of a generic cosmological fluid known as Van der Waals fluid in the context of flat FRW universe. It is treated as a perfect fluid that obeys the equation of state $P=\frac{\gamma\rho}{1-\beta\rho}-\alpha\rho^{2}, 0\leq\gamma<1$, where $\rho$ stands for energy density and $P$ stands for pressure of the fluid. In this regard, we discuss the behavior of physical parameters to analyze the evolution of the universe. We investigate whether the cosmological scenario fulfills the third law of thermodynamics using specific heat formalism. Next we discuss the thermal equation of state and by means of adiabatic, specific heat and isothermal conditions from classical thermodynamics we examine the thermal stability.


title: Additive Density-on-Scalar Regression in Bayes Hilbert Spaces with an Application to Gender Economics

abstract: Motivated by research on gender identity norms and the distribution of the woman's share in a couple's total labor income, we consider functional additive regression models for probability density functions as responses with scalar covariates. To preserve nonnegativity and integration to one under vector space operations, we formulate the model for densities in a Bayes Hilbert space, which allows to not only consider continuous densities, but also, e.g., discrete or mixed densities. Mixed ones occur in our application, as the woman's income share is a continuous variable having discrete point masses at zero and one for single-earner couples. Estimation is based on a gradient boosting algorithm, allowing for potentially numerous flexible covariate effects and model selection. We develop properties of Bayes Hilbert spaces related to subcompositional coherence, yielding (odds-ratio) interpretation of effect functions and simplified estimation for mixed densities via an orthogonal decomposition. Applying our approach to data from the German Socio-Economic Panel Study (SOEP) shows a more symmetric distribution in East German than in West German couples after reunification and a smaller child penalty comparing couples with and without minor children. These West-East differences become smaller, but are persistent over time.


title: Grounding force-directed network layouts with latent space models

abstract: Force-directed layout algorithms are ubiquitously-used tools for network visualisation across a multitude of scientific disciplines. However, they lack theoretical grounding which allows to interpret their outcomes rigorously and can guide the choice of specific algorithms for certain data sets. We propose an approach building on latent space models, which assume that the probability of nodes forming a tie depends on their distance in an unobserved latent space. From such latent space models, we derive force equations for a force-directed layout algorithm. Since the forces infer positions which maximise the likelihood of the given network under the latent space model, the force-directed layout becomes interpretable. We implement these forces for unweighted and weighted networks and spatialise different real-world networks. Comparison to existing layout algorithms (not grounded in an interpretable model) reveals that node groups are placed in similar configurations, while said algorithms show a stronger intra-cluster separation of nodes, as well as a tendency to separate clusters more strongly in retweet networks. We also explore the possibility of visualising data traditionally not seen as network data, such as survey data.


title: Sinkformers: Transformers with Doubly Stochastic Attention

abstract: Attention based models such as Transformers involve pairwise interactions between data points, modeled with a learnable attention matrix. Importantly, this attention matrix is normalized with the SoftMax operator, which makes it row-wise stochastic. In this paper, we propose instead to use Sinkhorn's algorithm to make attention matrices doubly stochastic. We call the resulting model a Sinkformer. We show that the row-wise stochastic attention matrices in classical Transformers get close to doubly stochastic matrices as the number of epochs increases, justifying the use of Sinkhorn normalization as an informative prior. On the theoretical side, we show that, unlike the SoftMax operation, this normalization makes it possible to understand the iterations of self-attention modules as a discretized gradient-flow for the Wasserstein metric. We also show in the infinite number of samples limit that, when rescaling both attention matrices and depth, Sinkformers operate a heat diffusion. On the experimental side, we show that Sinkformers enhance model accuracy in vision and natural language processing tasks. In particular, on 3D shapes classification, Sinkformers lead to a significant improvement.


title: Learning Stable Vector Fields on Lie Groups

abstract: Learning robot motions from demonstration requires models able to specify vector fields for the full robot pose when the task is defined in operational space. Recent advances in reactive motion generation have shown that learning adaptive, reactive, smooth, and stable vector fields is possible. However, these approaches define vector fields on a flat Euclidean manifold, while representing vector fields for orientations requires modeling the dynamics in non-Euclidean manifolds, such as Lie Groups. In this paper, we present a novel vector field model that can guarantee most of the properties of previous approaches i.e., stability, smoothness, and reactivity beyond the Euclidean space. In the experimental evaluation, we show the performance of our proposed vector field model to learn stable vector fields for full robot poses as SE(2) and SE(3) in both simulated and real robotics tasks.


title: Federated Learning over Wireless IoT Networks with Optimized Communication and Resources

abstract: To leverage massive distributed data and computation resources, machine learning in the network edge is considered to be a promising technique especially for large-scale model training. Federated learning (FL), as a paradigm of collaborative learning techniques, has obtained increasing research attention with the benefits of communication efficiency and improved data privacy. Due to the lossy communication channels and limited communication resources (e.g., bandwidth and power), it is of interest to investigate fast responding and accurate FL schemes over wireless systems. Hence, we investigate the problem of jointly optimized communication efficiency and resources for FL over wireless Internet of things (IoT) networks. To reduce complexity, we divide the overall optimization problem into two sub-problems, i.e., the client scheduling problem and the resource allocation problem. To reduce the communication costs for FL in wireless IoT networks, a new client scheduling policy is proposed by reusing stale local model parameters. To maximize successful information exchange over networks, a Lagrange multiplier method is first leveraged by decoupling variables including power variables, bandwidth variables and transmission indicators. Then a linear-search based power and bandwidth allocation method is developed. Given appropriate hyper-parameters, we show that the proposed communication-efficient federated learning (CEFL) framework converges at a strong linear rate. Through extensive experiments, it is revealed that the proposed CEFL framework substantially boosts both the communication efficiency and learning performance of both training loss and test accuracy for FL over wireless IoT networks compared to a basic FL approach with uniform resource allocation.


title: Matching-invariant running quark masses in Quantum Chromodynamics

abstract: The conventional quark mass is not continuous at thresholds. In this paper, we derive matchinginvariant quark masses which are continuous everywhere. They are expanded as an obvious function of the logarithmic Lambda scaled energy. The expansion coefficients are related to the original gamma and beta functions, with concretization to four loop level. The results show that the new expressions for the quark masses converge indeed much faster.


title: Modeling extinction and reddening effects by circumstellar dust in the Betelgeuse envelope in the presence of radiative torque disruption

abstract: Circumstellar dust is formed and evolved within the envelope of evolved stars, including Asymptotic Giant Branch (AGB) and Red Supergiant (RSG). The extinction of stellar light by circumstellar dust is vital for interpreting RSG/AGB observations and determining high-mass RSG progenitors of core-collapse supernovae. Nevertheless, circumstellar dust properties are not well understood. Modern understanding of dust evolution suggests that intense stellar radiation can radically change the dust properties across the circumstellar envelope through the RAdiative Torque Disruption (RAT-D) mechanism. In this paper, we study the impacts of RAT-D on the grain size distribution (GSD) of circumstellar dust and model its effects on photometric observations of $\alpha$ Orionis (Betelgeuse). Due to the RAT-D effects, large grains formed in the dust formation zone are disrupted into smaller species of size $a < 0.5\,\rm\mu m$. Using the GSD constrained by the RAT-D effects, we model the visual extinction of background stars and Betelgeuse. We find that the extinction decreases at near-UV, optical, and infrared wavelengths while increasing at far-UV wavelengths. The resulting flux potentially reproduces the observation from the near-UV to near-IR range. Our results can be used to explain dust extinction and photometric observations toward other RSG/AGB stars.


title: Domain Adaptation and Active Learning for Fine-Grained Recognition in the Field of Biodiversity

abstract: Deep-learning methods offer unsurpassed recognition performance in a wide range of domains, including fine-grained recognition tasks. However, in most problem areas there are insufficient annotated training samples. Therefore, the topic of transfer learning respectively domain adaptation is particularly important. In this work, we investigate to what extent unsupervised domain adaptation can be used for fine-grained recognition in a biodiversity context to learn a real-world classifier based on idealized training data, e.g. preserved butterflies and plants. Moreover, we investigate the influence of different normalization layers, such as Group Normalization in combination with Weight Standardization, on the classifier. We discovered that domain adaptation works very well for fine-grained recognition and that the normalization methods have a great influence on the results. Using domain adaptation and Transferable Normalization, the accuracy of the classifier could be increased by up to 12.35 % compared to the baseline. Furthermore, the domain adaptation system is combined with an active learning component to improve the results. We compare different active learning strategies with each other. Surprisingly, we found that more sophisticated strategies provide better results than the random selection baseline for only one of the two datasets. In this case, the distance and diversity strategy performed best. Finally, we present a problem analysis of the datasets.


title: SwiftLane: Towards Fast and Efficient Lane Detection

abstract: Recent work done on lane detection has been able to detect lanes accurately in complex scenarios, yet many fail to deliver real-time performance specifically with limited computational resources. In this work, we propose SwiftLane: a simple and light-weight, end-to-end deep learning based framework, coupled with the row-wise classification formulation for fast and efficient lane detection. This framework is supplemented with a false positive suppression algorithm and a curve fitting technique to further increase the accuracy. Our method achieves an inference speed of 411 frames per second, surpassing state-of-the-art in terms of speed while achieving comparable results in terms of accuracy on the popular CULane benchmark dataset. In addition, our proposed framework together with TensorRT optimization facilitates real-time lane detection on a Nvidia Jetson AGX Xavier as an embedded system while achieving a high inference speed of 56 frames per second.


title: Reconstruction of Sentinel-2 Time Series Using Robust Gaussian Mixture Models -- Application to the Detection of Anomalous Crop Development in wheat and rapeseed crops

abstract: Missing data is a recurrent problem in remote sensing, mainly due to cloud coverage for multispectral images and acquisition problems. This can be a critical issue for crop monitoring, especially for applications relying on machine learning techniques, which generally assume that the feature matrix does not have missing values. This paper proposes a Gaussian Mixture Model (GMM) for the reconstruction of parcel-level features extracted from multispectral images. A robust version of the GMM is also investigated, since datasets can be contaminated by inaccurate samples or features (e.g., wrong crop type reported, inaccurate boundaries, undetected clouds, etc). Additional features extracted from Synthetic Aperture Radar (SAR) images using Sentinel-1 data are also used to provide complementary information and improve the imputations. The robust GMM investigated in this work assigns reduced weights to the outliers during the estimation of the GMM parameters, which improves the final reconstruction. These weights are computed at each step of an Expectation-Maximization (EM) algorithm by using outlier scores provided by the isolation forest algorithm. Experimental validation is conducted on rapeseed and wheat parcels located in the Beauce region (France). Overall, we show that the GMM imputation method outperforms other reconstruction strategies. A mean absolute error (MAE) of 0.013 (resp. 0.019) is obtained for the imputation of the median Normalized Difference Index (NDVI) of the rapeseed (resp. wheat) parcels. Other indicators (e.g., Normalized Difference Water Index) and statistics (for instance the interquartile range, which captures heterogeneity among the parcel indicator) are reconstructed at the same time with good accuracy. In a dataset contaminated by irrelevant samples, using the robust GMM is recommended since the standard GMM imputation can lead to inaccurate imputed values.


title: Forcing axioms via ground model interpretations

abstract: We study principles of the form: if a name $\sigma$ is forced to have a certain property $\varphi$, then there is a ground model filter $g$ such that $\sigma^g$ satisfies $\varphi$. We prove a general correspondence connecting these name principles to forcing axioms. Special cases of the main theorem are: Any forcing axiom can be expressed as a name principle. For instance, $\mathsf{PFA}$ is equivalent to a principle for rank $1$ names (equivalently, nice names) for subsets of $\omega_1$, and a principle for rank $2$ names for sets of reals. Moreover, $\lambda$-bounded forcing axioms are equivalent to name principles. Bagaria's characterisation of $\mathsf{BFA}$ via generic absoluteness is a corollary. We further systematically study name principles where $\varphi$ is a notion of largeness for subsets of $\omega_1$ (such as being unbounded, stationary or in the club filter) and corresponding forcing axioms.


title: Multiple-ionization energy difference of 163Ho and 163Dy atoms

abstract: The multiple-ionization energy difference of 163Ho and 163Dy atoms is evaluated for the ionization degree q = 30, 48, 56. The calculations are performed by means of the large-scale relativistic configuration interaction method combined with the many-body perturbation theory. The quantum electrodynamics, nuclear recoil, and frequency-dependent Breit interaction corrections are taken into account. The obtained theoretical values are within 1 eV uncertainty. These results can help to increase accuracy of the laboratory neutrino mass limit, provided they are accompanied by the corresponding experiment on electron capture in Ho and a precise measurement of the ion mass difference.


title: Generic Poincar\'{e}-Bendixson Theorem for singularly perturbed monotone systems with respect to cones of rank-$2$

abstract: We investigate the singularly perturbed monotone systems with respect to cones of rank $2$ and obtain the so called Generic Poincar\'{e}-Bendixson theorem for such perturbed systems, that is, for a bounded positively invariant set, there exists an open and dense subset $\mathcal{P}$ such that for each $z\in \mathcal{P}$, the $\omega$-limit set $\omega(z)$ that contains no equilibrium points is a closed orbit.


title: Safe rules for the identification of zeros in the solutions of the SLOPE problem

abstract: In this paper we propose a methodology to accelerate the resolution of the so-called "Sorted L-One Penalized Estimation" (SLOPE) problem. Our method leverages the concept of "safe screening", well-studied in the literature for \textit{group-separable} sparsity-inducing norms, and aims at identifying the zeros in the solution of SLOPE. More specifically, we derive a set of \(\tfrac{n(n+1)}{2}\) inequalities for each element of the \(n\)-dimensional primal vector and prove that the latter can be safely screened if some subsets of these inequalities are verified. We propose moreover an efficient algorithm to jointly apply the proposed procedure to all the primal variables. Our procedure has a complexity \(\mathcal{O}(n\log n + LT)\) where \(T\leq n\) is a problem-dependent constant and \(L\) is the number of zeros identified by the tests. Numerical experiments confirm that, for a prescribed computational budget, the proposed methodology leads to significant improvements of the solving precision.


title: Double Proton Transfer in Hydrated Formic Acid Dimer: Interplay of Spatial Symmetry and Solvent-Generated Force on Reactivity

abstract: The double proton transfer (DPT) reaction in hydrated formic acid dimer (FAD) is investigated at molecular-level detail. For this, a global and reactive machine learned (ML) potential energy surface (PES) is developed to run extensive (more than 100 ns) mixed ML/MM molecular dynamics (MD) simulations in explicit molecular mechanics (MM) solvent at MP2-quality for the solute. Simulations with fixed - as in a conventional empirical force field - and conformationally fluctuating - as available from the ML-based PES - charge models for FAD shows significant impact on the competition between DPT and dissociation of FAD into two formic acid monomers. With increasing temperature the barrier height for DPT in solution changes by about 10% ($\sim 1$ kcal/mol) between 300 K and 600 K. The rate for DPT is largest, $\sim 1$ ns$^{-1}$, at 350 K and decreases for higher temperatures due to destabilisation and increased probability for dissociation of FAD. The water solvent is found to promote the first proton transfer by exerting a favourable solvent-induced Coulomb force along the O-H$\cdots$O hydrogen bond whereas the second proton transfer is significantly controlled by the O-O separation and other conformational degrees of freedom. Double proton transfer in hydrated FAD is found to involve a subtle interplay and balance between structural and electrostatic factors.


title: Surjective isometries on Banach sequence spaces: a survey

abstract: In this survey, we present several results related to characterizing the surjective isometries on Banach sequence spaces. Our survey includes full proofs of these characterizations for the classical spaces as well as more recent results for combinatorial Banach spaces and Tsirelson-type spaces. Along the way, we pose many open problems related to the structure of the group of surjective isometries and characterizations of the group of surjective isometries for various Banach spaces.


title: Emergence of a periodically rotating one-point cluster in a thermodynamic Cucker-Smale ensemble

abstract: We study emergent behaviors of thermomechanical Cucker-Smale (TCS) ensemble confined in a harmonic potential field. In the absence of external force field, emergent dynamics of TCS particles has been extensively studied recently under various frameworks formulated in terms of initial configuration, system parameters and network topologies. Moreover, the TCS model does not exhibit rotating motions in the absence of an external force field. In this paper, we show the emergence of periodically rotating one-point cluster for the TCS model in a harmonic potential field using elementary energy estimates and continuity argument. We also provide several numerical simulations and compare them with analytical results.


title: An analysis on metric-driven multi-target sensor management: GOSPA versus OSPA

abstract: This paper presents an analysis on sensor management using a cost function based on a multi-target metric, in particular, the optimal subpattern-assignment (OSPA) metric, the unnormalised OSPA (UOSPA) metric and the generalised OSPA (GOSPA) metric (\alpha=2). We consider the problem of managing an array of sensors, where each sensor is able to observe a region of the surveillance area, not covered by other sensors, with a given sensing cost. We look at the case in which there are far-away, independent potential targets, at maximum one per sensor region. In this set-up, the optimal action using GOSPA is taken for each sensor independently, as we may expect. On the contrary, as a consequence of the spooky effect at a distance in optimal OSPA/UOSPA estimation, the optimal actions for different sensors using OSPA and UOSPA are entangled.


title: From scalar fields on quantum spaces to blobbed topological recursion

abstract: We review the construction of the $\lambda\phi^4$-model on noncommutative geometries via exact solutions of Dyson-Schwinger equations and explain how this construction relates via (blobbed) topological recursion to problems in algebraic and enumerative geometry.


title: Automatic Guide Generation for Stan via NumPyro

abstract: Stan is a very popular probabilistic language with a state-of-the-art HMC sampler but it only offers a limited choice of algorithms for black-box variational inference. In this paper, we show that using our recently proposed compiler from Stan to Pyro, Stan users can easily try the set of algorithms implemented in Pyro for black-box variational inference. We evaluate our approach on PosteriorDB, a database of Stan models with corresponding data and reference posterior samples. Results show that the eight algorithms available in Pyro offer a range of possible compromises between complexity and accuracy. This paper illustrates that compiling Stan to another probabilistic language can be used to leverage new features for Stan users, and give access to a large set of examples for language developers who implement these new features.


title: Algebraic cycles and Fano threefolds of genus 8

abstract: We show that prime Fano threefolds $Y$ of genus 8 have a multiplicative Chow-K\"unneth decomposition, in the sense of Shen-Vial. As a consequence, a certain tautological subring of the Chow ring of powers of $Y$ injects into cohomology.


title: Evaluation of binomial series with harmonic numbers

abstract: We define a special function related to the digamma function and use it to evaluate in closed form various series involving binomial coefficients and harmonic numbers.


title: Optimality conditions for mathematical programs with orthogonality type constraints

abstract: We consider the class of mathematical programs with orthogonality type constraints (MPOC). Orthogonality type constraints appear by reformulating the sparsity constraint via auxiliary binary variables and relaxing them afterwards. For MPOC a necessary optimality condition in terms of T-stationarity is stated. The justification of T-stationarity is threefold. First, it allows to capture the global structure of MPOC in terms of Morse theory, i. e. deformation and cell-attachment results are established. For that, nondegeneracy for the T-stationary points is introduced and shown to hold at a generic MPOC. Second, we prove that Karush-Kuhn-Tucker points of the Scholtes-type regularization converge to T-stationary points of MPOC. This is done under the MPOC-tailored linear independence constraint qualification (LICQ), which turns out to be a generic property too. Third, we show that T-stationarity applied to the relaxation of sparsity constrained nonlinear optimization (SCNO) naturally leads to its M-stationary points. Moreover, we argue that all T-stationary points of this relaxation become degenerate.


title: Federated Unlearning via Class-Discriminative Pruning

abstract: We explore the problem of selectively forgetting categories from trained CNN classification models in the federated learning (FL). Given that the data used for training cannot be accessed globally in FL, our insights probe deep into the internal influence of each channel. Through the visualization of feature maps activated by different channels, we observe that different channels have a varying contribution to different categories in image classification. Inspired by this, we propose a method for scrubbing the model clean of information about particular categories. The method does not require retraining from scratch, nor global access to the data used for training. Instead, we introduce the concept of Term Frequency Inverse Document Frequency (TF-IDF) to quantize the class discrimination of channels. Channels with high TF-IDF scores have more discrimination on the target categories and thus need to be pruned to unlearn. The channel pruning is followed by a fine-tuning process to recover the performance of the pruned model. Evaluated on CIFAR10 dataset, our method accelerates the speed of unlearning by 8.9x for the ResNet model, and 7.9x for the VGG model under no degradation in accuracy, compared to retraining from scratch. For CIFAR100 dataset, the speedups are 9.9x and 8.4x, respectively. We envision this work as a complementary block for FL towards compliance with legal and ethical criteria.


title: HDRVideo-GAN: Deep Generative HDR Video Reconstruction

abstract: High dynamic range (HDR) videos provide a more visually realistic experience than the standard low dynamic range (LDR) videos. Despite having significant progress in HDR imaging, it is still a challenging task to capture high-quality HDR video with a conventional off-the-shelf camera. Existing approaches rely entirely on using dense optical flow between the neighboring LDR sequences to reconstruct an HDR frame. However, they lead to inconsistencies in color and exposure over time when applied to alternating exposures with noisy frames. In this paper, we propose an end-to-end GAN-based framework for HDR video reconstruction from LDR sequences with alternating exposures. We first extract clean LDR frames from noisy LDR video with alternating exposures with a denoising network trained in a self-supervised setting. Using optical flow, we then align the neighboring alternating-exposure frames to a reference frame and then reconstruct high-quality HDR frames in a complete adversarial setting. To further improve the robustness and quality of generated frames, we incorporate temporal stability-based regularization term along with content and style-based losses in the cost function during the training procedure. Experimental results demonstrate that our framework achieves state-of-the-art performance and generates superior quality HDR frames of a video over the existing methods.


title: Connes spectral distance and nonlocality of generalized noncommutative phase spaces

abstract: We study the Connes spectral distance of quantum states and analyse the nonlocality of the $4D$ generalized noncommutative phase space. By virtue of the Hilbert-Schmidt operatorial formulation, we obtain the Dirac operator and construct a spectral triple corresponding to the noncommutative phase space. Based on the ball condition, we obtain some constraint relations about the optimal elements, and then calculate the Connes spectral distance between two Fock states. Due to the noncommutativity, the spectral distances between Fock states in generalized noncommutative phase space are shorter than those in normal phase space. This shortening of distances implies some kind of nonlocality caused by the noncommutativity. We also find that these spectral distances in $4D$ generalized noncommutative phase space are additive and satisfy the normal Pythagoras theorem. When the noncommutative parameters equal zero, the results return to those in normal quantum phase space.


title: Flexible Physical Layer Security for Joint Data and Pilots in Future Wireless Networks

abstract: In this work, novel physical layer security (PLS) schemes are proposed for orthogonal frequency-division multiplexing (OFDM) to secure both data and pilots. The majority of previous studies focus on only securing the data without considering the security of the pilots used for channel estimation. However, the leakage of channel state information (CSI) from a legitimate node to an eavesdropper allows the latter to acquire knowledge about the channel of the legitimate nodes. To this end, we propose adaptive and flexible PLS algorithms which can 1) secure data, 2) secure pilots, and 3) jointly secure both data and pilots. Particularly, minimum-phase all-pass channel decomposition is exploited, where the proposed algorithms use the all-pass component to provide security without harming the performance of the legitimate user. In the analysis for data security, we evaluated the secrecy under correlated and uncorrelated eavesdropping channels via closed-form bit error rate (BER) formulas. For pilot security, we analyzed the normalized mean squared error (NMSE) performance of the estimated channel. The simulation results along with theoretical analysis demonstrate that the proposed algorithms can effectively enhance the communication secrecy of the overall system.


title: Resolution-vs.-Accuracy Dilemma in Machine Learning Modeling of Electronic Excitation Spectra

abstract: In this study, we explore the potential of machine learning for modeling molecular electronic spectral intensities as a continuous function in a given wavelength range. Since presently available chemical space datasets provide excitation energies and corresponding oscillator strengths for only a few valence transitions, here, we present a new dataset -- \bigqm -- with 12,880 molecules containing up to 7 CONF atoms and report ground state and excited state properties. A publicly accessible web-based data-mining platform is presented to facilitate on-the-fly screening of several molecular properties including harmonic vibrational and electronic spectra. We present all singlet electronic transitions from the ground state calculated using the time-dependent density functional theory framework with the $\omega$B97XD exchange-correlation functional and a diffuse-function augmented basis set. The resulting spectra predominantly span the X-ray to deep-UV region (10--120 nm). To compare the target spectra with predictions based on small basis sets, we bin spectral intensities and show good agreement is obtained only at the expense of the resolution. Compared to this, machine learning models with latest structural representations trained directly using $<10 \%$ of the target data recover the spectra of the remaining molecules with better accuracies at a desirable $<1$ nm wavelength resolution.


title: Toward a multidimensional analysis of transmission spectroscopy. Part I: Computation of transmission spectra using a 1D, 2D, or 3D atmosphere structure

abstract: Considering the relatively high precision that will be reached by future observatories, it has recently become clear that one dimensional (1D) atmospheric models, in which the atmospheric temperature and composition of a planet are considered to vary only in the vertical, will be unable to represent exoplanetary transmission spectra with a sufficient accuracy. This is particularly true for warm to (ultra-) hot exoplanets because the atmosphere is unable to redistribute all the energy deposited on the dayside, creating a strong thermal and often compositional dichotomy on the planet. This situation is exacerbated by transmission spectroscopy, which probes the terminator region. This is the most heterogeneous region of the atmosphere. However, if being able to compute transmission spectra from 3D atmospheric structures (from a global climate model, e.g.) is necessary to predict realistic observables, it is too computationally expensive to be used in a data inversion framework. For this reason, there is a need for a medium-complexity 2D approach that captures the most salient features of the 3D model in a sufficiently fast implementation. With this in mind, we present a new open-source documented version of Pytmosph3R that handles the computation of transmission spectra for atmospheres with up to three spatial dimensions and can account for time variability. Taking the example of an ultra hot Jupiter, we illustrate how the changing orientation of the planet during the transit can allow us to probe the horizontal variations in the atmosphere. We further implement our algorithm in TauREx to allow the community to perform 2D retrievals. We describe our extensive cross-validation benchmarks and discuss the accuracy and numerical performance of each model.


title: Monohedral Tilings of a Convex Disc with a Smooth Boundary

abstract: In this paper we give a complete description about normal monohedral tilings of a convex disc with smooth boundary where we have at most three topological discs as tiles. This result is a far-reaching generalization of the results of Kurusa, L\'angi and V\'igh \cite{KLV2020}. Some further partial results are proved for non-normal tilings.


title: Computational methods for t-spread monomial ideals

abstract: Let $K$ be a field and $S=K[x_1,\ldots,x_n]$ a standard polynomial ring over $K$. In this paper, some new optimized algorithms to compute the smallest $t$-spread lexicographic set and the smallest $t$-spread strongly stable set containing a given set of $t$-spread monomials of $S$ are presented. Some technical tools allowing to compute the cardinality of $t$-spread strongly stable sets avoiding their construction are given. Then, a \emph{Macaulay2} package, \texttt{TSpreadIdeals}, providing methods to easily manage $t$-spread monomials and $t$-spread ideals is implemented. Some functions to ease the calculation of well known results about algebraic invariants for $t$-spread ideals are also provided.


title: Deep Convolutional Autoencoders as Generic Feature Extractors in Seismological Applications

abstract: The idea of using a deep autoencoder to encode seismic waveform features and then use them in different seismological applications is appealing. In this paper, we designed tests to evaluate this idea of using autoencoders as feature extractors for different seismological applications, such as event discrimination (i.e., earthquake vs. noise waveforms, earthquake vs. explosion waveforms, and phase picking). These tests involve training an autoencoder, either undercomplete or overcomplete, on a large amount of earthquake waveforms, and then using the trained encoder as a feature extractor with subsequent application layers (either a fully connected layer, or a convolutional layer plus a fully connected layer) to make the decision. By comparing the performance of these newly designed models against the baseline models trained from scratch, we conclude that the autoencoder feature extractor approach may only perform well under certain conditions such as when the target problems require features to be similar to the autoencoder encoded features, when a relatively small amount of training data is available, and when certain model structures and training strategies are utilized. The model structure that works best in all these tests is an overcomplete autoencoder with a convolutional layer and a fully connected layer to make the estimation.


title: Validation of point process predictions with proper scoring rules

abstract: We introduce a class of proper scoring rules for evaluating spatial point process forecasts based on summary statistics. These scoring rules rely on Monte-Carlo approximations of expectations and can therefore easily be evaluated for any point process model that can be simulated. In this regard, they are more flexible than the commonly used logarithmic score and other existing proper scores for point process predictions. The scoring rules allow for evaluating the calibration of a model to specific aspects of a point process, such as its spatial distribution or tendency towards clustering. Using simulations we analyze the sensitivity of our scoring rules to different aspects of the forecasts and compare it to the logarithmic score. Applications to earthquake occurrences in northern California, USA and the spatial distribution of Pacific silver firs in Findley Lake Reserve in Washington, USA highlight the usefulness of our scores for scientific model selection.


title: Probabilistic fine-tuning of pruning masks and PAC-Bayes self-bounded learning

abstract: We study an approach to learning pruning masks by optimizing the expected loss of stochastic pruning masks, i.e., masks which zero out each weight independently with some weight-specific probability. We analyze the training dynamics of the induced stochastic predictor in the setting of linear regression, and observe a data-adaptive L1 regularization term, in contrast to the dataadaptive L2 regularization term known to underlie dropout in linear regression. We also observe a preference to prune weights that are less well-aligned with the data labels. We evaluate probabilistic fine-tuning for optimizing stochastic pruning masks for neural networks, starting from masks produced by several baselines. In each case, we see improvements in test error over baselines, even after we threshold fine-tuned stochastic pruning masks. Finally, since a stochastic pruning mask induces a stochastic neural network, we consider training the weights and/or pruning probabilities simultaneously to minimize a PAC-Bayes bound on generalization error. Using data-dependent priors, we obtain a selfbounded learning algorithm with strong performance and numerically tight bounds. In the linear model, we show that a PAC-Bayes generalization error bound is controlled by the magnitude of the change in feature alignment between the 'prior' and 'posterior' data.


title: Model, sample, and epoch-wise descents: exact solution of gradient flow in the random feature model

abstract: Recent evidence has shown the existence of a so-called double-descent and even triple-descent behavior for the generalization error of deep-learning models. This important phenomenon commonly appears in implemented neural network architectures, and also seems to emerge in epoch-wise curves during the training process. A recent line of research has highlighted that random matrix tools can be used to obtain precise analytical asymptotics of the generalization (and training) errors of the random feature model. In this contribution, we analyze the whole temporal behavior of the generalization and training errors under gradient flow for the random feature model. We show that in the asymptotic limit of large system size the full time-evolution path of both errors can be calculated analytically. This allows us to observe how the double and triple descents develop over time, if and when early stopping is an option, and also observe time-wise descent structures. Our techniques are based on Cauchy complex integral representations of the errors together with recent random matrix methods based on linear pencils.


title: Multi-Horizon Planning of Multi-Energy Systems

abstract: In order to reach EU's goal of zero emissions in 2050, the energy system will go through a significant transition over the next decades. To substitute fossil energy carriers, renewable energy sources will be mainly integrated in the power system. Thereby, sector coupling will play a major role by making flexibility from other sectors such as heat or transport accessible to the power system. Planning the cost optimal transition requires a whole-system view over multiple horizons and across all sectors. This imposes the need for multi-energy system (MES) models coupled with multi-horizon investment models. This paper presents two multi-horizon planning approaches to determine the cost optimal pathway of a MES. As a major contribution, we propose a new method to incorporate technology-dependent learning cost curves in the planning problem and show that the resulting mixed-integer linear programming problem can be solved faster with a Benders decomposition technique as compared to a closed optimization. As a further contribution, we demonstrate the usefulness of our approach by showing the MES expansion pathway for a small German test system.


title: Signal-Envelope: A C++ library with Python bindings for temporal envelope estimation

abstract: Signals can be interpreted as composed of a rapidly varying component modulated by a slower varying envelope. Identifying this envelope is an essential operation in signal processing, with applications in areas ranging from seismology to medicine. Conventional envelope detection approaches based on classic methods tend to lack generality, however, and need to be tailored to each specific application in order to yield reasonable results. Taking inspiration from geometric concepts, most notably the theory of alpha-shapes, we introduce a general-purpose library to efficiently extract the envelope of arbitrary signals.


title: On the design of regularized explicit predictive controllers from input-output data

abstract: On the wave of recent advances in data-driven predictive control, we present an explicit predictive controller that can be constructed from a batch of input/output data only. The proposed explicit law is build upon a regularized implicit data-driven predictive control problem, so as to guarantee the uniqueness of the explicit predictive controller. As a side benefit, the use of regularization is shown to improve the capability of the explicit law in coping with noise on the data. The effectiveness of the retrieved explicit law and the repercussions of regularization on noise handling are analyzed on two benchmark simulation case studies, showing the potential of the proposed regularized explicit controller.


title: PropMix: Hard Sample Filtering and Proportional MixUp for Learning with Noisy Labels

abstract: The most competitive noisy label learning methods rely on an unsupervised classification of clean and noisy samples, where samples classified as noisy are re-labelled and "MixMatched" with the clean samples. These methods have two issues in large noise rate problems: 1) the noisy set is more likely to contain hard samples that are in-correctly re-labelled, and 2) the number of samples produced by MixMatch tends to be reduced because it is constrained by the small clean set size. In this paper, we introduce the learning algorithm PropMix to handle the issues above. PropMix filters out hard noisy samples, with the goal of increasing the likelihood of correctly re-labelling the easy noisy samples. Also, PropMix places clean and re-labelled easy noisy samples in a training set that is augmented with MixUp, removing the clean set size constraint and including a large proportion of correctly re-labelled easy noisy samples. We also include self-supervised pre-training to improve robustness to high noisy label scenarios. Our experiments show that PropMix has state-of-the-art (SOTA) results on CIFAR-10/-100(with symmetric, asymmetric and semantic label noise), Red Mini-ImageNet (from the Controlled Noisy Web Labels), Clothing1M and WebVision. In severe label noise bench-marks, our results are substantially better than other methods. The code is available athttps://github.com/filipe-research/PropMix.


title: IVS3D: An Open Source Framework for Intelligent Video Sampling and Preprocessing to Facilitate 3D Reconstruction

abstract: The creation of detailed 3D models is relevant for a wide range of applications such as navigation in three-dimensional space, construction planning or disaster assessment. However, the complex processing and long execution time for detailed 3D reconstructions require the original database to be reduced in order to obtain a result in reasonable time. In this paper we therefore present our framework iVS3D for intelligent pre-processing of image sequences. Our software is able to down sample entire videos to a specific frame rate, as well as to resize and crop the individual images. Furthermore, thanks to our modular architecture, it is easy to develop and integrate plugins with additional algorithms. We provide three plugins as baseline methods that enable an intelligent selection of suitable images and can enrich them with additional information. To filter out images affected by motion blur, we developed a plugin that detects these frames and also searches the spatial neighbourhood for suitable images as replacements. The second plugin uses optical flow to detect redundant images caused by a temporarily stationary camera. In our experiments, we show how this approach leads to a more balanced image sampling if the camera speed varies, and that excluding such redundant images leads to a time saving of 8.1\percent for our sequences. A third plugin makes it possible to exclude challenging image regions from the 3D reconstruction by performing semantic segmentation. As we think that the community can greatly benefit from such an approach, we will publish our framework and the developed plugins open source using the MIT licence to allow co-development and easy extension.


title: Angular structure factor of the hexatic-B liquid crystals: bridging theory and experiment

abstract: We report results from X-ray scattering studies of the angular structure factor of liquid crystal hexatic-B films. According to the sixfold rotational symmetry of the hexatic-B phase, its characteristic scattering splits into six reflections. The shape of the radial and angular cross-sections of these reflections and their temperature evolution are analyzed. We find that over a wide temperature range of the hexatic-B phase existence the angular profiles of the in-plane X-ray scattering are well fitted by the Voigt function, which is a convolution of the Gaussian and Lorentzian functions. This result is supported by the known theoretical considerations of the hexatic structure factor below the smectic-hexatic phase transition temperture. Similar predictions for the angular shape of the hexatic peak in the vicinity of the smectic-hexatic phase transition temperature follow from the multicritical scaling theory of the hexatic-B phase in three dimensions. We find that the specific shape of the hexatic structure factor can be explained by the interplay of two distinct contributions to the free energy of the system, a liquid-like density term and a coupling term between the bond-orientational order and short-range density fluctuations.


title: Probabilistic ODE Solutions in Millions of Dimensions

abstract: Probabilistic solvers for ordinary differential equations (ODEs) have emerged as an efficient framework for uncertainty quantification and inference on dynamical systems. In this work, we explain the mathematical assumptions and detailed implementation schemes behind solving {high-dimensional} ODEs with a probabilistic numerical algorithm. This has not been possible before due to matrix-matrix operations in each solver step, but is crucial for scientifically relevant problems -- most importantly, the solution of discretised {partial} differential equations. In a nutshell, efficient high-dimensional probabilistic ODE solutions build either on independence assumptions or on Kronecker structure in the prior model. We evaluate the resulting efficiency on a range of problems, including the probabilistic numerical simulation of a differential equation with millions of dimensions.


title: Handling Concurrency in Behavior Trees

abstract: This paper addresses the concurrency issues affecting Behavior Trees (BTs), a popular tool to model the behaviors of autonomous agents in the video game and the robotics industry. BT designers can easily build complex behaviors composing simpler ones, which represents a key advantage of BTs. The parallel composition of BTs expresses a way to combine concurrent behaviors that has high potential, since composing pre-existing BTs in parallel results easier than composing in parallel classical control architectures, as finite state machines or teleo-reactive programs. However, BT designers rarely use such composition due to the underlying concurrency problems similar to the ones faced in concurrent programming. As a result, the parallel composition, despite its potential, finds application only in the composition of simple behaviors or where the designer can guarantee the absence of conflicts by design. In this paper, we define two new BT nodes to tackle the concurrency problems in BTs and we show how to exploit them to create predictable behaviors. In addition, we introduce measures to assess execution performance and show how different design choices affect them. We validate our approach in both simulations and the real world. Simulated experiments provide statistically significant data, whereas real-world experiments show the applicability of our method on real robots. We provided an open-source implementation of the novel BT formulation and published all the source code to reproduce the numerical examples and experiments.


title: Identifiability of local and global features of phylogenetic networks from average distances

abstract: Phylogenetic networks extend phylogenetic trees to model non-vertical inheritance, by which a lineage inherits material from multiple parents. The computational complexity of estimating phylogenetic networks from genome-wide data with likelihood-based methods limits the size of networks that can be handled. Methods based on pairwise distances could offer faster alternatives. We study here the information that average pairwise distances contain on the underlying phylogenetic network, by characterizing local and global features that can or cannot be identified. For general networks, we clarify that the root and edge lengths adjacent to reticulations are not identifiable, and then focus on the class of zipped-up semidirected networks. We provide a criterion to swap subgraphs locally, such as 3-cycles, resulting in indistinguishable networks. We propose the "distance split tree", which can be constructed from pairwise distances, and prove that it is a refinement of the network's tree of blobs, capturing the tree-like features of the network. For level-1 networks, this distance split tree is equal to the tree of blobs refined to separate polytomies from blobs, and we prove that the mixed representation of the network is identifiable. The information loss is localized around 4-cycles, for which the placement of the reticulation is unidentifiable. The mixed representation combines split edges for 4-cycles, regular tree and hybrid edges from the semidirected network, and edge parameters that encode all information identifiable from average pairwise distances.


title: cleanTS: Automated (AutoML) Tool to Clean Univariate Time Series at Microscales

abstract: Data cleaning is one of the most important tasks in data analysis processes. One of the perennial challenges in data analytics is the detection and handling of non-valid data. Failing to do so can result in inaccurate analytics and unreliable decisions. The process of properly cleaning such data takes much time. Errors are prevalent in time series data. It is usually found that real world data is unclean and requires some pre-processing. The analysis of large amounts of data is difficult. This paper is intended to provide an easy to use and reliable system which automates the cleaning process of univariate time series data. Automating the process greatly reduces the time required. Visualizing a large amount of data at once is not very effective. To tackle this issue, an R package cleanTS is proposed. The proposed system provides a way to analyze data on different scales and resolutions. Also, it provides users with tools and a benchmark system for comparing various techniques used in data cleaning.


title: Testing network correlation efficiently via counting trees

abstract: We propose a new procedure for testing whether two networks are edge-correlated through some latent vertex correspondence. The test statistic is based on counting the co-occurrences of signed trees for a family of non-isomorphic trees. When the two networks are Erd\H{o}s-R\'enyi random graphs $\mathcal{G}(n,q)$ that are either independent or correlated with correlation coefficient $\rho$, our test runs in $n^{2+o(1)}$ time and succeeds with high probability as $n\to\infty$, provided that $n\min\{q,1-q\} \ge n^{-o(1)}$ and $\rho^2>\alpha \approx 0.338$, where $\alpha$ is Otter's constant so that the number of unlabeled trees with $K$ edges grows as $(1/\alpha)^K$. This significantly improves the prior work in terms of statistical accuracy, running time, and graph sparsity.


title: Multimessenger Astronomy with Neutrinos

abstract: Multimessenger astronomy is arguably the branch of the astroparticle physics field that has seen the most significant developments in recent years. In this manuscript, we will review the state-of-the-art, the recent observations, and the prospects and challenges for the near future. We will give special emphasis to the observation carried out with neutrino telescopes.


title: Characterizations of Stability of Error Bounds for Convex Inequality Constraint Systems

abstract: In this paper, we mainly study error bounds for a single convex inequality and semi-infinite convex constraint systems, and give characterizations of stability of error bounds via directional derivatives. For a single convex inequality, it is proved that the stability of local error bounds under small perturbations is essentially equivalent to the non-zero minimun of the directional derivative at a reference point over the sphere, and the stability of global error bounds is proved to be equivalent to the strictly positive infimum of the directional derivatives, at all points in the boundary of the solution set, over the sphere as well as some mild constraint qualification. When these results are applied to semi-infinite convex constraint systems, characterizations of stability of local and global error bounds under small perturbations are also provided. In particular such stability of error bounds is proved to only require that all component functions in semi-infinite convex constraint systems have the same linear perturbation. Our work demonstrates that verifying the stability of error bounds for convex inequality constraint systems is, to some degree, equivalent to solving the convex optimization/minimization problems (defined by directional derivatives) over the sphere.


title: A Last Switch Dependent Analysis of Satiation and Seasonality in Bandits

abstract: Motivated by the fact that humans like some level of unpredictability or novelty, and might therefore get quickly bored when interacting with a stationary policy, we introduce a novel non-stationary bandit problem, where the expected reward of an arm is fully determined by the time elapsed since the arm last took part in a switch of actions. Our model generalizes previous notions of delay-dependent rewards, and also relaxes most assumptions on the reward function. This enables the modeling of phenomena such as progressive satiation and periodic behaviours. Building upon the Combinatorial Semi-Bandits (CSB) framework, we design an algorithm and prove a bound on its regret with respect to the optimal non-stationary policy (which is NP-hard to compute). Similarly to previous works, our regret analysis is based on defining and solving an appropriate trade-off between approximation and estimation. Preliminary experiments confirm the superiority of our algorithm over both the oracle greedy approach and a vanilla CSB solver.


title: Endpoint estimates for bilinear pseudodifferential operators with symbol in BS^m_{1,1}

abstract: In this paper we establish some endpoint estimates for bilinearpseudodifferential operators with symbol in the class BS^m_{1,1}, involving the space of functions with local bounded mean oscillation bmo. As a consequence we also obtain an endpoint estimate of Kato-Ponce type.


title: Node-based Generalized Friendship Paradox fails

abstract: The Friendship Paradox--the principle that "your friends have more friends than you do"--is a combinatorial fact about degrees in a graph; but given that many web-based social activities are correlated with a user's degree, this fact has been taken more broadly to suggest the empirical principle that "your friends are also more active than you are." This Generalized Friendship Paradox, the notion that any attribute positively correlated with degree obeys the Friendship Paradox, has been established mathematically in a network-level version that essentially aggregates uniformly over all the edges of a network. Here we show, however, that the natural node-based version of the Generalized Friendship Paradox--which aggregates over nodes, not edges--may fail, even for degree-attribute correlations approaching 1. Whether this version holds depends not only on degree-attribute correlations, but also on the underlying network structure and thus can't be said to be a universal phenomenon. We establish both positive and negative results for this node-based version of the Generalized Friendship Paradox and consider its implications for social-network data.


title: Unraveling the Hidden Environmental Impacts of AI Solutions for Environment

abstract: In the past ten years, artificial intelligence has encountered such dramatic progress that it is now seen as a tool of choice to solve environmental issues and in the first place greenhouse gas emissions (GHG). At the same time the deep learning community began to realize that training models with more and more parameters requires a lot of energy and as a consequence GHG emissions. To our knowledge, questioning the complete net environmental impacts of AI solutions for the environment (AI for Green), and not only GHG, has never been addressed directly. In this article, we propose to study the possible negative impacts of AI for Green. First, we review the different types of AI impacts, then we present the different methodologies used to assess those impacts, and show how to apply life cycle assessment to AI services. Finally, we discuss how to assess the environmental usefulness of a general AI service, and point out the limitations of existing work in AI for Green.


title: Tunable spin-charge conversion in class-I topological Dirac semimetals

abstract: We theoretically demonstrate that class-I topological Dirac semimetals (TDSMs) can provide a platform for realizing both electrically and magnetically tunable spin-charge conversion. With time-reversal symmetry, the spin component along the uniaxial rotation axis ($z$-axis) is approximately conserved, which leads to an anisotropic spin Hall effect -- the resulting spin Hall current relies on the relative orientation between the external electric field and the $z$-axis. The application of a magnetic field, on the other hand, breaks time-reversal symmetry, driving the TDSM into a Weyl semimetal phase and, consequently, partially converting the spin current to a charge Hall current. Using the Kubo formulas, we numerically evaluate the spin and charge Hall conductivities based on a low-energy TDSM Hamiltonian together with the Zeeman coupling. Besides the conventional tensor element of the spin Hall conductivity $\sigma_{xy}^z$, we find that unconventional components, such as $\sigma_{xy}^x$ and $\sigma_{xy}^y$, also exist and vary as the magnetic field is rotated. Likewise, the charge Hall conductivity also exhibits appreciable tunability upon variation of the magnetic field. We show that such tunability -- as well as large spin-charge conversion efficiency -- arises from the interplay of symmetry and band topology of the TDSMs.


title: The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

abstract: Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope era high resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constrains in the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models (Marley et al 2021, Morley et al in prep.). In this paper we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide separation giant planets with key molecular species such as CH4, H2O, CO and NH3 in disequilibrium. Our grid covers atmospheres with Teff~[500 K,1300 K], logg~[3.0,5.5] (cgs) and an eddy diffusion parameter of logKzz=2, 4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH4, H2O, CO and NH3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.


title: Annihilating Entanglement Between Cones

abstract: Every multipartite entangled quantum state becomes fully separable after an entanglement breaking quantum channel acted locally on each of its subsystems. Whether there are other quantum channels with this property has been an open problem with important implications for entanglement theory (e.g., for the distillation problem and the PPT squared conjecture). We cast this problem in the general setting of proper convex cones in finite-dimensional vector spaces. The entanglement annihilating maps transform the $k$-fold maximal tensor product of a cone $C_1$ into the $k$-fold minimal tensor product of a cone $C_2$, and the pair $(C_1,C_2)$ is called resilient if all entanglement annihilating maps are entanglement breaking. Our main result is that $(C_1,C_2)$ is resilient if either $C_1$ or $C_2$ is a Lorentz cone. Our proof exploits the symmetries of the Lorentz cones and applies two constructions resembling protocols for entanglement distillation: As a warm-up, we use the multiplication tensors of real composition algebras to construct a finite family of generalized distillation protocols for Lorentz cones, containing the distillation protocol for entangled qubit states by Bennett et al. as a special case. Then, we construct an infinite family of protocols using solutions to the Hurwitz matrix equations. After proving these results, we focus on maps between cones of positive semidefinite matrices, where we derive necessary conditions for entanglement annihilation similar to the reduction criterion in entanglement distillation. Finally, we apply results from the theory of Banach space tensor norms to show that the Lorentz cones are the only cones with a symmetric base for which a certain stronger version of the resilience property is satisfied.


title: Predictive machine learning for prescriptive applications: a coupled training-validating approach

abstract: In this research we propose a new method for training predictive machine learning models for prescriptive applications. This approach, which we refer to as coupled validation, is based on tweaking the validation step in the standard training-validating-testing scheme. Specifically, the coupled method considers the prescription loss as the objective for hyper-parameter calibration. This method allows for intelligent introduction of bias in the prediction stage to improve decision making at the prescriptive stage, and is generally applicable to most machine learning methods, including recently proposed hybrid prediction-stochastic-optimization techniques, and can be easily implemented without model-specific mathematical modeling. Several experiments with synthetic and real data demonstrate promising results in reducing the prescription costs in both deterministic and stochastic models.


title: Uniquely Decodable Multi-Amplitude Sequence for Grant-Free Multiple-Access Adder Channels

abstract: Grant-free multiple-access (GFMA) is a valuable research topic, since it can support multiuser transmission with low latency. This paper constructs novel uniquely-decodable multi-amplitude sequence (UDAS) sets for GFMA systems, which can provide high spectrum efficiency (SE) with low-complexity active user detection (AUD) algorithm. First of all, we propose an UDAS-based multi-dimensional bit interleaving coded modulation (MD-BICM) transmitter; then introduce the definition of UDAS and construct two kinds of UDAS sets based on cyclic and quasi-cyclic matrix modes. Besides, we present a statistic of UDAS feature based AUD algorithm (SoF-AUD), and a joint multiuser detection and improved message passing algorithm for the proposed system. Finally, the active user error rate (AUER) and Shannon limits of the proposed system are deduced in details. Simulation results show that our proposed system can simultaneously support four users without additional redundancy, and the AUER can reach an extremely low value $10^{-5}$ when $E_b/N_0$ is $0$ dB and the length of transmit block is larger than a given value, i.e., 784, verifying the validity and flexibility of the proposed UDAS sets.


title: Magnetization dynamics fingerprints of an excitonic condensate $t_{2g}^{4}$ magnet

abstract: The competition between spin-orbit coupling $\lambda$ and electron-electron interaction $U$ leads to a plethora of novel states of matter, extensively studied in the context of $t_{2g}^4$ and $t_{2g}^5$ materials, such as ruthenates and iridates. Excitonic magnets -- the antiferromagnetic state of bounded electron-hole pairs -- is a prominent example of phenomena driven by those competing energy scales. Interestingly, recent theoretical studies predicted that excitonic magnets can be found in the ground-state of spin-orbit-coupled $t_{2g}^4$ Hubbard models. Here, we present a detailed computational study of the magnetic excitations in that excitonic magnet, employing one-dimensional chains (via density matrix renormalization group) and small two-dimensional clusters (via Lanczos). Specifically, first we show that the low-energy spectrum is dominated by a dispersive (acoustic) magnonic mode, with extra features arising from the $\lambda=0$ state in the phase diagram. Second, and more importantly, we found a novel magnetic excitation forming a high-energy optical mode with the highest intensity at wavevector $q\to 0$. In the excitonic condensation regime at large $U$, we also have found a novel high-energy $\pi$-mode composed solely of orbital excitations. These unique fingerprints of the $t_{2g}^4$ excitonic magnet are important in the analysis of neutron and RIXS experiments.


title: CNN-based Omnidirectional Object Detection for HermesBot Autonomous Delivery Robot with Preliminary Frame Classification

abstract: Mobile autonomous robots include numerous sensors for environment perception. Cameras are an essential tool for robot's localization, navigation, and obstacle avoidance. To process a large flow of data from the sensors, it is necessary to optimize algorithms, or to utilize substantial computational power. In our work, we propose an algorithm for optimizing a neural network for object detection using preliminary binary frame classification. An autonomous outdoor mobile robot with 6 rolling-shutter cameras on the perimeter providing a 360-degree field of view was used as the experimental setup. The obtained experimental results revealed that the proposed optimization accelerates the inference time of the neural network in the cases with up to 5 out of 6 cameras containing target objects.


title: Multi-attribute Pizza Generator: Cross-domain Attribute Control with Conditional StyleGAN

abstract: Multi-attribute conditional image generation is a challenging problem in computervision. We propose Multi-attribute Pizza Generator (MPG), a conditional Generative Neural Network (GAN) framework for synthesizing images from a trichotomy of attributes: content, view-geometry, and implicit visual style. We design MPG by extending the state-of-the-art StyleGAN2, using a new conditioning technique that guides the intermediate feature maps to learn multi-scale multi-attribute entangled representationsof controlling attributes. Because of the complex nature of the multi-attribute image generation problem, we regularize the image generation by predicting the explicit conditioning attributes (ingredients and view). To synthesize a pizza image with view attributesoutside the range of natural training images, we design a CGI pizza dataset PizzaView using 3D pizza models and employ it to train a view attribute regressor to regularize the generation process, bridging the real and CGI training datasets. To verify the efficacy of MPG, we test it on Pizza10, a carefully annotated multi-ingredient pizza image dataset. MPG can successfully generate photo-realistic pizza images with desired ingredients and view attributes, beyond the range of those observed in real-world training data.


title: Characterization of the measurement uncertainty dynamics in an open system

abstract: Uncertainty principle plays a crucial role in quantum mechanics, because it captures the essence of the inevitable randomness associated with the outcomes of two incompatible quantum measurements. Information entropy can perfectly describe this type of randomness in information theory, because entropy can measure the degree of chaos in a given quantum system. However, the quantum-assisted uncertainty of entropy eventually inflate inevitably as the quantum correlations of the system are progressively corrupted by noise from the surrounding environment. In this paper, we investigate the dynamical features of the von Neumann entropic uncertainty in the presence of quantum memory, exploring the time evolution of entropic uncertainty suffer noise from the surrounding. Noteworthily, how the environmental noises affect the uncertainty of entropy is revealed, and specifically we verify how two types of noise environments, AD channel and BPF channel, influence the uncertainty. Meanwhile, we put forward some effective operation strategies to reduce the magnitude of the measurement uncertainty under the open systems. Furthermore, we explore the applications of the uncertainty relation investigated on entanglement witness and quantum channel capacity. Therefore , in open quantum correlation systems, our investigations could provide an insight into quantum measurement estimation


title: Global Tracking and Quantification of Oil and Gas Methane Emissions from Recurrent Sentinel-2 Imagery

abstract: Methane (CH4) emissions estimates from top-down studies over oil and gas basins have revealed systematic under-estimation of CH4 emissions in current national inventories. Sparse but extremely large amounts of CH4 from oil and gas production activities have been detected across the globe, resulting in a significant increase of the overall O&G contribution. However, attribution to specific facilities remains a major challenge unless high-resolution images provide the sufficient granularity within O&G basin. In this paper, we monitor known oil-and-gas infrastructures across the globe using recurrent Sentinel-2 imagery to detect and quantify more than 800 CH4 emissions. In combination with emissions estimates from airborne and Sentinel-5P measurements, we demonstrate the robustness of the fit to a power law from 0.1 tCH4/hr to 600 tCH4/hr. We conclude here that the prevalence of ultra-emitters (> 25tCH4/hr) detected globally by Sentinel-5P directly relates to emission occurrences below its detection threshold. Similar power law coefficients arise from several major oil and gas producers but noticeable differences in emissions magnitudes suggest large differences in maintenance practices and infrastructures across countries.


title: Refined decay bounds on the entries of spectral projectors associated with sparse Hermitian matrices

abstract: Spectral projectors of Hermitian matrices play a key role in many applications, and especially in electronic structure computations. Linear scaling methods for gapped systems are based on the fact that these special matrix functions are localized, which means that the entries decay exponentially away from the main diagonal or with respect to more general sparsity patterns. The relation with the sign function together with an integral representation is used to obtain new decay bounds, which turn out to be optimal in an asymptotic sense. The influence of isolated eigenvalues in the spectrum on the decay properties is also investigated and a superexponential behaviour is predicted.


title: Deflection of MEMS Based Sandwiched Cantilever Beam for Piezoelectric Actuation: Analysis of Lead Zirconate Titanate Piezoceramic Material

abstract: This paper investigates the deflection of micro- electro-mechanical system (MEMS) based sandwiched cantilever beam for piezoelectric actuation using eight of the popular lead zirconatetitanate (PZT)piezoceramic material. The simulation for the investigation was carried away in COMSOL Multiphysics software environment. For the investigation of a certain model, the lowest tip deflection is found to be 38.074 nm for PZT-8 and the highest tip deflection is found to be 82.965 nm for PZT-5H. A 2 mm thick flexible foam core is sandwiched between two 8 mm thick aluminum layers in this model of sandwiched cantilever beam, which is 100 mm long. The gadget is replaced with a foam core with a 10 mm long piezoelectric actuator that is located between x=55 mm and x=65 mm. Along the global x-axis, the beam is aligned. The beam changes substantially when the piezoceramic material is changed. The tip deflection of the cantilever beam varies between 37 and 83 nanometers for different lead zirconate titanate piezoceramic materials. PZT-8 has the lowest deflection of 37.074 nm, whereas PZT-7A has the second lowest deflection of 40.423 nm. The maximum tip deflection discovered is more than twice as high as the lowest tip deflection discovered. The deflection of a sandwich cantilever beam for piezoelectric actuation is clearly influenced by piezoceramic materials which can be used in increasing the efficiency of any mechanical device.


title: Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

abstract: With the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections are combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyse them with a set of four different waveform models. Based on our analysis with about 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to $\pm 250\rm m$ ($2\%$ at $90\%$ credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.


title: The Geometry of Tree-Based Sorting

abstract: We study the connections between sorting and the binary search tree (BST) model, with an aim towards showing that the fields are connected more deeply than is currently appreciated. While any BST can be used to sort by inserting the keys one-by-one, this is a very limited relationship and importantly says nothing about parallel sorting. We show what we believe to be the first formal relationship between the BST model and sorting. Namely, we show that a large class of sorting algorithms, which includes mergesort, quicksort, insertion sort, and almost every instance-optimal sorting algorithm, are equivalent in cost to offline BST algorithms. Our main theoretical tool is the geometric interpretation of the BST model introduced by Demaine et al., which finds an equivalence between searches on a BST and point sets in the plane satisfying a certain property. To give an example of the utility of our approach, we introduce the log-interleave bound, a measure of the information-theoretic complexity of a permutation $\pi$, which is within a $\lg \lg n$ multiplicative factor of a known lower bound in the BST model; we also devise a parallel sorting algorithm with polylogarithmic span that sorts a permutation $\pi$ using comparisons proportional to its log-interleave bound. Our aforementioned result on sorting and offline BST algorithms can be used to show existence of an offline BST algorithm whose cost is within a constant factor of the log-interleave bound of any permutation $\pi$.


title: Physically-motivated basis functions for temperature maps of exoplanets

abstract: Thermal phase curves of exoplanet atmospheres have revealed temperature maps as a function of planetary longitude, often by sinusoidal decomposition of the phase curve. We construct a framework for describing two-dimensional temperature maps of exoplanets with mathematical basis functions derived for a fluid layer on a rotating, heated sphere with drag/friction, which are generalizations of spherical harmonics. These basis functions naturally produce physically-motivated temperature maps for exoplanets with few free parameters. We investigate best practices for applying this framework to temperature maps of hot Jupiters by splitting the problem into two parts: (1) we constrain the temperature map as a function of latitude by tuning the basis functions to reproduce general circulation model (GCM) outputs, since disk-integrated phase curve observations do not constrain this dimension; and (2) we infer the temperature maps of real hot Jupiters using original reductions of several Spitzer phase curves, which directly constrain the temperature variations with longitude. The resulting phase curves can be described with only three free parameters per bandpass -- an efficiency improvement over the usual five or so used to describe sinusoidal decompositions of phase curves. Upon obtaining the hemispherically averaged dayside and nightside temperatures, the standard approach would be to use zero-dimensional box models to infer the Bond albedo and redistribution efficiency. We elucidate the limitation of these box models by demonstrating that negative Bond albedos may be obtained due to a choice of boundary condition on the nightside temperature. We propose generalized definitions for the Bond albedo and heat redistribution efficiency for use with two-dimensional (2D) temperature maps. Open-source software called kelp is provided to efficiently compute these phase curves.


title: Experimental and numerical study of the effect of surface patterning on the frictional properties of polymer surfaces

abstract: We describe benchmark experiments to evaluate the frictional properties of laser patterned low-density polyethylene as a function of sliding velocity, normal force and humidity. The pattern is a square lattice of square cavities with sub-mm spacing. We find that dynamic friction decreases compared to non-patterned surfaces, since stress concentrations lead to anticipated detachment, and that stick-slip behavior is also affected. Friction increases with humidity, and the onset of stick-slip events occurs in the high humidity regime. Experimental results are compared with numerical simulations of a simplified 2-D spring-block model. A good qualitative agreement can be obtained by introducing a deviation from the linear behavior of the Amontons-Coulomb law with the load, due to a saturation in the effective contact area with pressure. This also leads also to the improvement of the quantitative results of the spring-block model by reducing the discrepancy with the experimental results, indicating the robustness of the adopted simplified approach, which could be adopted to design patterned surfaces with controlled friction properties.


title: $\texttt{GRASS}$: Distinguishing Planet-induced Doppler Signatures from Granulation with a Synthetic Spectra Generator

abstract: Owing to recent advances in radial-velocity instrumentation and observation techniques, the detection of Earth-mass planets around Sun-like stars may soon be primarily limited by intrinsic stellar variability. Several processes contribute to this variability, including starspots, pulsations, and granulation. Although many previous studies have focused on techniques to mitigate signals from pulsations and other types of magnetic activity, granulation noise has to date only been partially addressed by empirically-motivated observation strategies and magnetohydrodynamic simulations. To address this deficit, we present the GRanulation And Spectrum Simulator ($\texttt{GRASS}$), a new tool designed to create time-series synthetic spectra with granulation-driven variability from spatially- and temporally-resolved observations of solar absorption lines. In this work, we present $\texttt{GRASS}$, detail its methodology, and validate its model against disk-integrated solar observations. As a first-of-its-kind empirical model for spectral variability due to granulation in a star with perfectly known center-of-mass radial-velocity behavior, $\texttt{GRASS}$ is an important tool for testing new methods of disentangling granular line-shape changes from true Doppler shifts.


title: Variational Bayesian Approximation of Inverse Problems using Sparse Precision Matrices

abstract: Inverse problems involving partial differential equations (PDEs) are widely used in science and engineering. Although such problems are generally ill-posed, different regularisation approaches have been developed to ameliorate this problem. Among them is the Bayesian formulation, where a prior probability measure is placed on the quantity of interest. The resulting posterior probability measure is usually analytically intractable. The Markov Chain Monte Carlo (MCMC) method has been the go-to method for sampling from those posterior measures. MCMC is computationally infeasible for large-scale problems that arise in engineering practice. Lately, Variational Bayes (VB) has been recognised as a more computationally tractable method for Bayesian inference, approximating a Bayesian posterior distribution with a simpler trial distribution by solving an optimisation problem. In this work, we argue, through an empirical assessment, that VB methods are a flexible and efficient alternative to MCMC for this class of problems. We propose a natural choice of a family of Gaussian trial distributions parametrised by precision matrices, thus taking advantage of the inherent sparsity of the inverse problem encoded in its finite element discretisation. We utilise stochastic optimisation to efficiently estimate the variational objective and assess not only the error in the solution mean but also the ability to quantify the uncertainty of the estimate. We test this on PDEs based on the Poisson equation in 1D and 2D. A Tensorflow implementation is made publicly available on GitHub.


title: What drives the scatter of local star-forming galaxies in the BPT diagrams? A Machine Learning based analysis

abstract: We investigate which physical properties are most predictive of the position of local star forming galaxies on the BPT diagrams, by means of different Machine Learning (ML) algorithms. Exploiting the large statistics from the Sloan Digital Sky Survey (SDSS), we define a framework in which the deviation of star-forming galaxies from their median sequence can be described in terms of the relative variations in a variety of observational parameters. We train artificial neural networks (ANN) and random forest (RF) trees to predict whether galaxies are offset above or below the sequence (via classification), and to estimate the exact magnitude of the offset itself (via regression). We find, with high significance, that parameters primarily associated to variations in the nitrogen-over-oxygen abundance ratio (N/O) are the most predictive for the [N II]-BPT diagram, whereas properties related to star formation (like variations in SFR or EW[H$\alpha$]) perform better in the [S II]-BPT diagram. We interpret the former as a reflection of the N/O-O/H relationship for local galaxies, while the latter as primarily tracing the variation in the effective size of the S$^{+}$ emitting region, which directly impacts the [S II]emission lines. This analysis paves the way to assess to what extent the physics shaping local BPT diagrams is also responsible for the offsets seen in high redshift galaxies or, instead, whether a different framework or even different mechanisms need to be invoked.


title: Multi-view Contrastive Graph Clustering

abstract: With the explosive growth of information technology, multi-view graph data have become increasingly prevalent and valuable. Most existing multi-view clustering techniques either focus on the scenario of multiple graphs or multi-view attributes. In this paper, we propose a generic framework to cluster multi-view attributed graph data. Specifically, inspired by the success of contrastive learning, we propose multi-view contrastive graph clustering (MCGC) method to learn a consensus graph since the original graph could be noisy or incomplete and is not directly applicable. Our method composes of two key steps: we first filter out the undesirable high-frequency noise while preserving the graph geometric features via graph filtering and obtain a smooth representation of nodes; we then learn a consensus graph regularized by graph contrastive loss. Results on several benchmark datasets show the superiority of our method with respect to state-of-the-art approaches. In particular, our simple approach outperforms existing deep learning-based methods.