title: Deep matrix factorizations

abstract: Constrained low-rank matrix approximations have been known for decades as powerful linear dimensionality reduction techniques to be able to extract the information contained in large data sets in a relevant way. However, such low-rank approaches are unable to mine complex, interleaved features that underlie hierarchical semantics. Recently, deep matrix factorization (deep MF) was introduced to deal with the extraction of several layers of features and has been shown to reach outstanding performances on unsupervised tasks. Deep MF was motivated by the success of deep learning, as it is conceptually close to some neural networks paradigms. In this paper, we present the main models, algorithms, and applications of deep MF through a comprehensive literature review. We also discuss theoretical questions and perspectives of research.


title: Student-Initiated Action Advising via Advice Novelty

abstract: Action advising is a budget-constrained knowledge exchange mechanism between teacher-student peers that can help tackle exploration and sample inefficiency problems in deep reinforcement learning (RL). Most recently, student-initiated techniques that utilise state novelty and uncertainty estimations have obtained promising results. However, the approaches built on these estimations have some potential weaknesses. First, they assume that the convergence of the student's RL model implies less need for advice. This can be misleading in scenarios with teacher absence early on where the student is likely to learn suboptimally by itself; yet also ignore the teacher's assistance later. Secondly, the delays between encountering states and having them to take effect in the RL model updates in presence of the experience replay dynamics cause a feedback lag in what the student actually needs advice for. We propose a student-initiated algorithm that alleviates these by employing Random Network Distillation (RND) to measure the novelty of a piece of advice. Furthermore, we perform RND updates only for the advised states to ensure that the student's own learning does not impair its ability to leverage the teacher. Experiments in GridWorld and MinAtar show that our approach performs on par with the state-of-the-art and demonstrates significant advantages in the scenarios where the existing methods are prone to fail.


title: A Deep Learning Framework for COVID Outbreak Prediction

abstract: The outbreak of COVID-19 i.e. a variation of coronavirus, also known as novel corona virus causing respiratory disease is a big concern worldwide since the end of December 2019. As of September 12, 2020, it has turned into an epidemic outbreak with more than 29 million confirmed cases and around 1 million reported deaths worldwide. It has created an urgent need to monitor and forecast COVID-19 spread behavior to better control this spread. Among all the popular models for COVID-19 forecasting, statistical models are receiving much attention in media. However, statistical models are showing less accuracy for long term forecasting, as there is high level of uncertainty and required data is also not sufficiently available. In this paper, we propose a comparative analysis of deep learning models to forecast the COVID-19 outbreak as an alternative to statistical models. We propose a new Attention-based encoder-decoder model, named Attention-Long Short Term Memory (AttentionLSTM). LSTM based neural network layer architecture incorporates the idea of fine-grained attention mechanism i.e., attention on hidden state dimensions instead of hidden state vector itself, which is capable of highlighting the importance and contribution of each hidden state dimension. It helps in detection on crucial temporal information, resulting in a highly interpretable network. Additionally, we implement a learnable vector embedding for time. As, time in a vector representation can be easily added with many architectures. This vector representation is called Time2Vec. We have used COVID-19 data repository by the Center for Systems Science and Engineering (CSSE) at Johns Hopkins University to assess the proposed model's performance. The proposed model give superior forecasting accuracy compared to other existing methods.


title: Multi-photon enhancement of the Schwinger pair production mechanism under strong FEL radiation

abstract: In this work we study the production of electron-positron pairs from vacuum in presence of a field resulting from the collision of two counter-propagating FEL beams that undergo Gaussian amplitude fluctuations. Our work aims to the extension of previous works based on the standing wave hypothesis, by including the inherent stochastic amplitude fluctuations of the individual FEL beams. As shown, depending on the order of the process, a large non-linear enhancement in the number of created pairs can be expected over a big intensity window in the multi-photon regime. Vacuum pair creation in view of future plans on the production of ultra-strong and high energy radiation in FEL facilities is also discussed.


title: Minimum Dissipation Theorem for Microswimmers

abstract: We derive a theorem for the lower bound on the energy dissipation rate by a rigid surface-driven active microswimmer of arbitrary shape in a fluid at low Reynolds number. We show that, for any swimmer, the minimum dissipation at a given velocity can be expressed in terms of the resistance tensors of two passive bodies of the same shape with a no-slip and perfect-slip boundary. To achieve the absolute minimum dissipation, the optimal swimmer needs a surface velocity profile that corresponds to the flow around the perfect-slip body, and a propulsive force density that corresponds to the no-slip body. Using this theorem, we propose an alternative definition of the energetic efficiency of microswimmers that, unlike the commonly-used Lighthill efficiency, can never exceed unity. We validate the theory by calculating the efficiency limits of spheroidal swimmers.


title: Computational Approaches for Grocery Home Delivery Services

abstract: Grocery home delivery services require customers to be present when their deliveries arrive. Hence, the grocery retailer and the customer must mutually agree on a time window during which the delivery can be guaranteed. This concept is referred to as the Attended Home Delivery (AHD) problem. The phase during which customers place orders, usually through a web service, constitutes the computationally most challenging part of the logistical processes behind such services. The system must determine potential delivery time windows that can be offered to incoming customers and incrementally build the delivery schedule as new orders are placed. Typically, a Vehicle Routing Problem with Time Windows forms the underlying optimization problem. This work is concerned with the use case given by an international grocery retailer's online shopping service. We present an analysis of efficient solution methods that can be employed to AHD services. We provide several heuristic approaches for tackling the steps mentioned above. However, the basic framework can be easily be adapted to be used for many similar Vehicle Routing applications. We provide a comprehensive computational study comparing several algorithmic strategies, combining heuristics utilizing Local Search operations and Mixed-Integer Linear Programs, tackling the booking process. Finally, we analyze the scalability and suitability of the approaches.


title: Herding stochastic autonomous agents via local control rules and online global target selection strategies

abstract: In this Paper we propose a simple yet effective set of local control rules to make a group of "herder agents" collect and contain in a desired region an ensemble of non-cooperative stochastic "target agents" in the plane. We investigate the robustness of the proposed strategies to variations of the number of target agents and the strength of the repulsive force they feel when in proximity of the herders. Extensive numerical simulations confirm the effectiveness of the approach and are complemented by a more realistic validation on commercially available robotic agents via ROS.


title: Incorporating network based protein complex discovery into automated model construction

abstract: We propose a method for gene expression based analysis of cancer phenotypes incorporating network biology knowledge through unsupervised construction of computational graphs. The structural construction of the computational graphs is driven by the use of topological clustering algorithms on protein-protein networks which incorporate inductive biases stemming from network biology research in protein complex discovery. This structurally constrains the hypothesis space over the possible computational graph factorisation whose parameters can then be learned through supervised or unsupervised task settings. The sparse construction of the computational graph enables the differential protein complex activity analysis whilst also interpreting the individual contributions of genes/proteins involved in each individual protein complex. In our experiments analysing a variety of cancer phenotypes, we show that the proposed methods outperform SVM, Fully-Connected MLP, and Randomly-Connected MLPs in all tasks. Our work introduces a scalable method for incorporating large interaction networks as prior knowledge to drive the construction of powerful computational models amenable to introspective study.


title: Slope plate of sticky soil granular slope instability based on complex network

abstract: The particle discrete element simulation of the instability and failure process of the granular slope accumulator model when the metal plate continues downward is obtained, and the two-dimensional total velocity vector of soil particle velocity and slope slip during the instability and failure of the slope accumulator are obtained. Macro-response processes such as removing the angle of the crack surface and the average velocity in the y-direction of the slope top of the slope accumulation body. Construct a normal force chain undirected network model of the slope accumulation body particles under natural accumulation, and study the location of its slip surface, and The results are compared with the experimental results. Finally, the complex network method is used to analyze the topological characteristics of the contact force chain network of the particles on the slope top of the slope accumulation body, and the average degree, clustering coefficient and average shortest path are obtained during the slope instability of the slope accumulation body. The evolutionary rule of the method is used to verify its accuracy in combination with the strength reduction method. The research results show that the average shortest path can provide a more effective early warning of the instability and failure of slope deposits. A complex network theory is used to study the macro response of the slope deposits and its force chain. The interrelationship between the macroscopic structure of the network provides a new mathematical analysis method for the study of slope instability.


title: A Survey on Explainability in Machine Reading Comprehension

abstract: This paper presents a systematic review of benchmarks and approaches for explainability in Machine Reading Comprehension (MRC). We present how the representation and inference challenges evolved and the steps which were taken to tackle these challenges. We also present the evaluation methodologies to assess the performance of explainable systems. In addition, we identify persisting open research questions and highlight critical directions for future work.


title: Sustained Coherence Characteristics of Failure Process of Shallow Buried Tunnel under Multiple Explosive Loads based on Persistent Homology

abstract: The damage characteristics of a shallow buried tunnel under multiple explosive loads is an important research issue in the design and evaluation of protective engineering. It is of great significance to develop a method for early warning of the safety of the shallow buried features. The discrete element method is used to establish a mechanical model of the shallow buried tunnel. The South Load Equivalent Principle treats blast loads as a series of dynamic forces acting uniformly on the surface. Based on the discrete element method, the dynamic response after each blast load and the damage evolution process of the surrounding rock of the tunnel are obtained. The strength reduction method is used to obtain the surrounding rock of the tunnel. Introduce the theory of continuous homology, and use the mathematical method of continuous homology to quantitatively and qualitatively analyze the failure characteristics of the discrete element model under multiple explosive loads. The results show that the method of continuous homology can accurately reflect the topological characteristics of the surrounding rock of the tunnel The maximum one-dimensional bar code connection radius can effectively warn tunnel instability. This provides a new mathematical method for tunnel safety design and disaster prediction research.


title: Topological feature study of slope failure process via persistent homology-based machine learning

abstract: Using software UDEC to simulate the instability failure process of slope under seismic load, studing the dynamic response of slope failure, obtaining the deformation characteristics and displacement cloud map of slope, then analyzing the instability state of slope by using the theory of persistent homology, generates bar code map and extracts the topological characteristics of slope from bar code map. The topological characteristics corresponding to the critical state of slope instability are found, and the relationship between topological characteristics and instability evolution is established. Finally, it provides a topological research tool for slope failure prediction. The results show that the change of the longest Betti 1 bar code reflects the evolution process of the slope and the law of instability failure. Using discrete element method and persistent homology theory to study the failure characteristics of slope under external load can better understand the failure mechanism of slope, provide theoretical basis for engineering protection, and also provide a new mathematical method for slope safety design and disaster prediction research.


title: Phase retrieval with Bregman divergences and application to audio signal recovery

abstract: Phase retrieval (PR) aims to recover a signal from the magnitudes of a set of inner products. This problem arises in many audio signal processing applications which operate on a short-time Fourier transform magnitude or power spectrogram, and discard the phase information. Recovering the missing phase from the resulting modified spectrogram is indeed necessary in order to synthesize time-domain signals. PR is commonly addressed by considering a minimization problem involving a quadratic loss function. In this paper, we adopt a different standpoint. Indeed, the quadratic loss does not properly account for some perceptual properties of audio, and alternative discrepancy measures such as beta-divergences have been preferred in many settings. Therefore, we formulate PR as a new minimization problem involving Bregman divergences. Since these divergences are not symmetric with respect to their two input arguments in general, they lead to two different formulations of the problem. To optimize the resulting objective, we derive two algorithms based on accelerated gradient descent and alternating direction method of multipliers. Experiments conducted on audio signal recovery from spectrograms that are either exact or estimated from noisy observations highlight the potential of our proposed methods for audio restoration. In particular, leveraging some of these Bregman divergences induce better performance than the quadratic loss when performing PR from spectrograms under very noisy conditions.


title: Length scales and scale-free dynamics of dislocations in dense solid solutions

abstract: The fundamental interactions between an edge dislocation and a random solid solution are studied by analyzing dislocation line roughness profiles obtained from molecular dynamics simulations of Fe0.70Ni0.11 Cr0.19 over a range of stresses and temperatures. These roughness profiles reveal the hallmark features of a depinning transition. Namely, below a temperature-dependent critical stress, the dislocation line exhibits roughness in two different length scale regimes which are divided by a so-called correlation length. This correlation length increases with applied stress and at the critical stress (depinning transition or yield stress) formally goes to infinity. Above the critical stress, the line roughness profile converges to that of a random noise field. Motivated by these results, a physical model is developed based on the notion of coherent line bowing over all length scales below the correlation length. Above the correlation length, the solute field prohibits such coherent line bow outs. Using this model, we identify potential gaps in existing theories of solid solution strengthening and show that recent observations of length-dependent dislocation mobilities can be rationalized.


title: Apparent convergence of Pad\'e approximants for the crossover line in finite density QCD

abstract: We propose a novel Bayesian method to analytically continue observables to real baryochemical potential $\mu_B$ in finite density QCD. Taylor coefficients at $\mu_B=0$ and data at imaginary chemical potential $\mu_B^I$ are treated on equal footing. We consider two different constructions for the Pad\'e approximants, the classical multipoint Pad\'e approximation and a mixed approximation that is a slight generalization of a recent idea in Pad\'e approximation theory. Approximants with spurious poles are excluded from the analysis. As an application, we perform a joint analysis of the available continuum extrapolated lattice data for both pseudocritical temperature $T_c$ at $\mu_B^I$ from the Wuppertal-Budapest Collaboration and Taylor coefficients $\kappa_2$ and $\kappa_4$ from the HotQCD Collaboration. An apparent convergence of $[p/p]$ and $[p/p+1]$ sequences of rational functions is observed with increasing $p.$ We present our extrapolation up to $\mu_B\approx 600$ MeV.


title: The Minkowski problem in the Gaussian probability space

abstract: The Minkowski problem in Gaussian probability space is studied in this paper. In addition to providing an existence result on a Gaussian-volume-normalized version of this problem, the main goal of the current work is to provide uniqueness and existence results on the Gaussian Minkowski problem (with no normalization required).


title: High energy neutrino and gamma-ray emission in the jets of the microquasar M33 X-7

abstract: In this work, after testing the reliability of our algorithms through numerical simulations on the well-studied SS 433 Galactic microquasar, we focus on neutrino and $\gamma$-ray emissions from the extragalactic M33 X-7 system. This is a recently discovered X-ray binary system located in the neighbouring galaxy Messier 33 which has not yet been modelled in detail. The neutrino and $\gamma$-ray energy spectra, produced from the magnetized astrophysical jet of M33 X-7, in the context of our method are assumed to originate from the decay (and scattering) processes taking place among the secondary particles produced assuming that, first, hot (relativistic) protons of the jet scatter on thermal ones (p-p interaction mechanism).


title: Upscaling the porosity-permeability relationship of a microporous carbonate to the Darcy scale with machine learning

abstract: The permeability of a pore structure is typically described by stochastic representations of its geometrical attributes. Database-driven numerical solvers for large model domains can only accurately predict large-scale flow behaviour when they incorporate upscaled descriptions of that structure. The upscaling is particularly challenging for rocks with multimodal porosity structures such as carbonates, where several different types of structures are interacting. It is the connectivity both within and between these different structures that controls the porosity-permeability relationship at the larger length scales. Recent advances in machine learning combined with numerical modelling and structural analysis have allowed us to probe the relationship between structure and permeability more deeply. We have used this integrated approach to tackle the challenge of upscaling multimodal and multiscale porous media. We present a novel method for upscaling multimodal porosity-permeability relationships using machine learning based multivariate structural regression. A m-CT image of limestone was divided into sub-volumes and permeability was computed using the DBS model. The porosity-permeability relationship from Menke et al. was used to assign permeability values to the microporosity. Structural attributes of each sub-volume were extracted and then regressed against the solved permeability using an Extra-Trees regression model to derive an upscaled porosity-permeability relationship. Ten upscaled test cases were then modelled at the Darcy scale using the regression and benchmarked against full DBS simulations, a numerically upscaled Darcy model, and a K-C fit. We found good agreement between the full DBS simulations and both the numerical and machine learning upscaled models while the K-C model was a poor predictor in all cases.


title: Centrality-Based Traffic Restriction in Delayed Epidemic Networks

abstract: During an epidemic, infectious individuals might not be detectable until some time after becoming infected. The studies show that carriers with mild or no symptoms are the main contributors to the transmission of a virus within the population. The average time it takes to develop the symptoms causes a delay in the spread dynamics of the disease. When considering the influence of delay on the disease propagation in epidemic networks, depending on the value of the time-delay and the network topology, the peak of epidemic could be considerably different in time, duration, and intensity. Motivated by the recent worldwide outbreak of the COVID-19 virus and the topological extent in which this virus has spread over the course of a few months, this study aims to highlight the effect of time-delay in the progress of such infectious diseases in the meta-population networks rather than individuals or a single population. In this regard, the notions of epidemic network centrality in terms of the underlying interaction graph of the network, structure of the uncertainties, and symptom development duration are investigated to establish a centrality-based analysis of the disease evolution. A convex traffic volume optimization method is then developed to control the outbreak. The control process is done by identifying the sub-populations with the highest centrality and then isolating them while maintaining the same overall traffic volume (motivated by economic considerations) in the meta-population level. The numerical results, along with the theoretical expectations, highlight the impact of time-delay as well as the importance of considering the worst-case scenarios in investigating the most effective methods of epidemic containment.


title: Analyzing Koopman approaches to physics-informed machine learning for long-term sea-surface temperature forecasting

abstract: Accurately predicting sea-surface temperature weeks to months into the future is an important step toward long term weather forecasting. Standard atmosphere-ocean coupled numerical models provide accurate sea-surface forecasts on the scale of a few days to a few weeks, but many important weather systems require greater foresight. In this paper we propose machine-learning approaches sea-surface temperature forecasting that are accurate on the scale of dozens of weeks. Our approach is based in Koopman operator theory, a useful tool for dynamical systems modelling. With this approach, we predict sea surface temperature in the Gulf of Mexico up to 180 days into the future based on a present image of thermal conditions and three years of historical training data. We evaluate the combination of a basic Koopman method with a convolutional autoencoder, and a newly proposed "consistent Koopman" method, in various permutations. We show that the Koopman approach consistently outperforms baselines, and we discuss the utility of our additional assumptions and methods in this sea-surface temperature domain.


title: DeepFakesON-Phys: DeepFakes Detection based on Heart Rate Estimation

abstract: This work introduces a novel DeepFake detection framework based on physiological measurement. In particular, we consider information related to the heart rate using remote photoplethysmography (rPPG). rPPG methods analyze video sequences looking for subtle color changes in the human skin, revealing the presence of human blood under the tissues. In this work we investigate to what extent rPPG is useful for the detection of DeepFake videos. The proposed fake detector named DeepFakesON-Phys uses a Convolutional Attention Network (CAN), which extracts spatial and temporal information from video frames, analyzing and combining both sources to better detect fake videos. This detection approach has been experimentally evaluated using the latest public databases in the field: Celeb-DF and DFDC. The results achieved, above 98% AUC (Area Under the Curve) on both databases, outperform the state of the art and prove the success of fake detectors based on physiological measurement to detect the latest DeepFake videos.


title: Quasar Detection using Linear Support Vector Machine with Learning From Mistakes Methodology

abstract: The field of Astronomy requires the collection and assimilation of vast volumes of data. The data handling and processing problem has become severe as the sheer volume of data produced by scientific instruments each night grows exponentially. This problem becomes extensive for conventional methods of processing the data, which was mostly manual, but is the perfect setting for the use of Machine Learning approaches. While building classifiers for Astronomy, the cost of losing a rare object like supernovae or quasars to detection losses is far more severe than having many false positives, given the rarity and scientific value of these objects. In this paper, a Linear Support Vector Machine (LSVM) is explored to detect Quasars, which are extremely bright objects in which a supermassive black hole is surrounded by a luminous accretion disk. In Astronomy, it is vital to correctly identify quasars, as they are very rare in nature. Their rarity creates a class-imbalance problem that needs to be taken into consideration. The class-imbalance problem and high cost of misclassification are taken into account while designing the classifier. To achieve this detection, a novel classifier is explored, and its performance is evaluated. It was observed that LSVM along with Ensemble Bagged Trees (EBT) achieved a 10x reduction in the False Negative Rate, using the Learning from Mistakes methodology.


title: From Trees to Continuous Embeddings and Back: Hyperbolic Hierarchical Clustering

abstract: Similarity-based Hierarchical Clustering (HC) is a classical unsupervised machine learning algorithm that has traditionally been solved with heuristic algorithms like Average-Linkage. Recently, Dasgupta reframed HC as a discrete optimization problem by introducing a global cost function measuring the quality of a given tree. In this work, we provide the first continuous relaxation of Dasgupta's discrete optimization problem with provable quality guarantees. The key idea of our method, HypHC, is showing a direct correspondence from discrete trees to continuous representations (via the hyperbolic embeddings of their leaf nodes) and back (via a decoding algorithm that maps leaf embeddings to a dendrogram), allowing us to search the space of discrete binary trees with continuous optimization. Building on analogies between trees and hyperbolic space, we derive a continuous analogue for the notion of lowest common ancestor, which leads to a continuous relaxation of Dasgupta's discrete objective. We can show that after decoding, the global minimizer of our continuous relaxation yields a discrete tree with a (1 + epsilon)-factor approximation for Dasgupta's optimal tree, where epsilon can be made arbitrarily small and controls optimization challenges. We experimentally evaluate HypHC on a variety of HC benchmarks and find that even approximate solutions found with gradient descent have superior clustering quality than agglomerative heuristics or other gradient based algorithms. Finally, we highlight the flexibility of HypHC using end-to-end training in a downstream classification task.


title: Mediating Artificial Intelligence Developments through Negative and Positive Incentives

abstract: The field of Artificial Intelligence (AI) is going through a period of great expectations, introducing a certain level of anxiety in research, business and also policy. This anxiety is further energised by an AI race narrative that makes people believe they might be missing out. Whether real or not, a belief in this narrative may be detrimental as some stake-holders will feel obliged to cut corners on safety precautions, or ignore societal consequences just to "win". Starting from a baseline model that describes a broad class of technology races where winners draw a significant benefit compared to others (such as AI advances, patent race, pharmaceutical technologies), we investigate here how positive (rewards) and negative (punishments) incentives may beneficially influence the outcomes. We uncover conditions in which punishment is either capable of reducing the development speed of unsafe participants or has the capacity to reduce innovation through over-regulation. Alternatively, we show that, in several scenarios, rewarding those that follow safety measures may increase the development speed while ensuring safe choices. Moreover, in {the latter} regimes, rewards do not suffer from the issue of over-regulation as is the case for punishment. Overall, our findings provide valuable insights into the nature and kinds of regulatory actions most suitable to improve safety compliance in the contexts of both smooth and sudden technological shifts.


title: Optical-microwave pump-probe studies of electronic properties in novel materials

abstract: Combined microwave-optical pump-probe methods are emerging to study the quantum state of spin qubit centers and the charge dynamics in semiconductors. A major hindrance is the limited bandwidth of microwave irradiation/detection circuitry which could be overcome with the use of broadband coplanar waveguides (CPW). We present the development and performance characterization of two spectrometers: an optically detected magnetic resonance spectrometer (ODMR) and a microwave detected photoconductivity measurement. In the first method light serves as detection and microwaves excite the investigated medium, while in the second the roles are interchanged. The performance is demonstrated by measuring ODMR maps on the nitrogen-vacancy center in diamond and time resolved photoconductivity in p-doped silicon. The results demonstrate both an efficient coupling of the microwave irradiation to the samples as well as an excellent sensitivity for minute changes in sample conductivity.


title: Krieger's type of nonsingular Poisson suspensions and IDPFT systems

abstract: Given an infinite countable discrete amenable group $\Gamma$, we construct explicitly sharply weak mixing nonsingular Poisson $\Gamma$-actions of each Krieger's type: $III_\lambda$, for $\lambda\in[0,1]$, and $II_\infty$. The result is new even for $\Gamma=\Bbb Z$. As these Poisson suspension actions are over very special dissipative base, we obtain also new examples of sharply weak mixing nonsingular Bernoulli $\Gamma$-actions and IDPFT systems of each possible Krieger's type.


title: Momentum via Primal Averaging: Theoretical Insights and Learning Rate Schedules for Non-Convex Optimization

abstract: Momentum methods are now used pervasively within the machine learning community for training non-convex models such as deep neural networks. Empirically, they out perform traditional stochastic gradient descent (SGD) approaches. In this work we develop a Lyapunov analysis of SGD with momentum (SGD+M), by utilizing a equivalent rewriting of the method known as the stochastic primal averaging (SPA) form. This analysis is much tighter than previous theory in the non-convex case, and due to this we are able to give precise insights into when SGD+M may out-perform SGD, and what hyper-parameter schedules will work and why.


title: A note on the large charge expansion in 4d CFT

abstract: In this letter, we discuss certain universal predictions of the large charge expansion in conformal field theories with $U(1)$ symmetry, mainly focusing on four-dimensional theories. We show that, while in three dimensions quantum fluctuations are responsible for the existence of a theory-independent $Q^0$ term in the scaling dimension $\Delta_Q$ of the lightest operator with fixed charge $Q\gg 1$, in four dimensions the same mechanism provides a universal $Q^0\log Q$ correction to $\Delta_Q$. Previous works discussing four-dimensional theories failed in identifying this term. We also compute the first subleading correction to the OPE coefficient corresponding to the insertion of an arbitrary primary operator with small charge $q\ll Q$ in between the minimal energy states with charge $Q$ and $Q+q$, both in three and four dimensions. This contribution does not depend on the operator insertion and, similarly to the quantum effects in $\Delta_Q$, in four dimensions it scales logarithmically with $Q$.


title: Estimation of copulas via Maximum Mean Discrepancy

abstract: This paper deals with robust inference for parametric copula models. Estimation using Canonical Maximum Likelihood might be unstable, especially in the presence of outliers. We propose to use a procedure based on the Maximum Mean Discrepancy (MMD) principle. We derive non-asymptotic oracle inequalities, consistency and asymptotic normality of this new estimator. In particular, the oracle inequality holds without any assumption on the copula family, and can be applied in the presence of outliers or under misspecification. Moreover, in our MMD framework, the statistical inference of copula models for which there exists no density with respect to the Lebesgue measure on $[0,1]^d$, as the Marshall-Olkin copula, becomes feasible. A simulation study shows the robustness of our new procedures, especially compared to pseudo-maximum likelihood estimation. An R package implementing the MMD estimator for copula models is available.


title: Optimized computation of tight focusing of short pulses using mapping to periodic space

abstract: When a pulsed, few-cycle electromagnetic wave is focused by optics with f-number smaller than two, the frequency components it contains are focused to different regions of space, building up a complex electromagnetic field structure. Accurate numerical computation of this structure is essential for many applications such as the analysis, diagnostics, and control of high-intensity laser-matter interactions. However, straightforward use of finite-difference methods can impose unacceptably high demands on computational resources, owing to the necessity of resolving far-field and near-field zones at sufficiently high resolution to overcome numerical dispersion effects. Here, we present a procedure for fast computation of tight focusing by mapping a spherically curved far-field region to periodic space, where the field can be advanced by a dispersion-free spectral solver. In many cases of interest, the mapping reduces both run time and memory requirements by a factor of order 10, making it possible to carry out simulations on a desktop machine or a single node of a supercomputer. We provide an open-source C++ implementation with Python bindings and demonstrate its use for a desktop machine, where the routine provides the opportunity to use the resolution sufficient for handling the pulses with spectra spanning over several octaves. The described approach can facilitate the stability analysis of theoretical proposals, the studies based on statistical inferences, as well as the overall development and analysis of experiments with tightly-focused short laser pulses.


title: Chiral Triclinic Metamaterial Crystals Supporting Isotropic Acoustical Activity and Isotropic Chiral Phonons

abstract: Recent work predicted the existence of isotropic chiral phonon dispersion relations of the lowest bands connected to isotropic acoustical activity in cubic crystalline approximants of 3D chiral icosahedral metamaterial quasicrystals. While these architectures are fairly broadband and presumably robust against fabrication tolerances due to orientation averaging, they are extremely complex, very hard to manufacture experimentally, and they show effects which are about an order of magnitude smaller compared to those of ordinary highly anisotropic chiral cubic metamaterial crystals. Here, we propose and analyze a chiral triclinic metamaterial crystal exhibiting broadband isotropic acoustical activity. These 3D truss lattices are much less complex and exhibit substantially larger effects than the 3D quasicrystals at the price of being somewhat more susceptible to fabrication tolerances. This susceptibility originates from the fact that we have tailored the lowest two transverse phonon bands to exhibit an "accidental" degeneracy in momentum space.


title: A Feasibility-Driven Approach to Control-Limited DDP

abstract: Differential dynamic programming (DDP) is a direct single shooting method for trajectory optimization. Its efficiency derives from the exploitation of temporal structure (inherent to optimal control problems) and explicit roll-out/integration of the system dynamics. However, it suffers from numerical instability and, when compared to direct multiple shooting methods, it has limited initialization options (allows initialization of controls, but not of states) and lacks proper handling of control constraints. In this work, we tackle these issues with a feasibility-driven approach that regulates the dynamic feasibility during the numerical optimization and ensures control limits. Our feasibility search emulates the numerical resolution of a direct multiple shooting problem with only dynamics constraints. We show that our approach (named BOX-FDDP) has better numerical convergence than BOX-DDP+ (a single shooting method), and that its convergence rate and runtime performance are competitive with state-of-the-art direct transcription formulations solved using the interior point and active set algorithms available in KNITRO. We further show that BOX-FDDP decreases the dynamic feasibility error monotonically--as in state-of-the-art nonlinear programming algorithms. We demonstrate the benefits of our approach by generating complex and athletic motions for quadruped and humanoid robots. Finally, we highlight that BOX-FDDP is suitable for model predictive control in legged robots.


title: Competitive Algorithms for the Online Multiple Knapsack Problem with Application to Electric Vehicle Charging

abstract: We introduce and study a general version of the fractional online knapsack problem with multiple knapsacks, heterogeneous constraints on which items can be assigned to which knapsack, and rate-limiting constraints on the assignment of items to knapsacks. This problem generalizes variations of the knapsack problem and of the one-way trading problem that have previously been treated separately, and additionally finds application to the real-time control of electric vehicle (EV) charging. We introduce a new algorithm that achieves a competitive ratio within an additive factor of one of the best achievable competitive ratios for the general problem and matches or improves upon the best-known competitive ratio for special cases in the knapsack and one-way trading literatures. Moreover, our analysis provides a novel approach to online algorithm design based on an instance-dependent primal-dual analysis that connects the identification of worst-case instances to the design of algorithms. Finally, we illustrate the proposed algorithm via trace-based experiments of EV charging.


title: Resilient Urban Housing Markets: Shocks vs. Fundamentals

abstract: In the face of a pandemic, urban protests, and an affordability crisis, is the desirability of dense urban settings at a turning point? Assessing cities' long term trends remains challenging. The first part of this chapter describes the short-run dynamics of the housing market in 2020. Evidence from prices and price-to-rent ratios suggests expectations of resilience. Zip-level evidence suggests a short-run trend towards suburbanization, and some impacts of urban protests on house prices. The second part of the chapter analyzes the long-run dynamics of urban growth between 1970 and 2010. It analyzes what, in such urban growth, is explained by short-run shocks as opposed to fundamentals such as education, industrial specialization, industrial diversification, urban segregation, and housing supply elasticity. This chapter's original results as well as a large established body of literature suggest that fundamentals are the key drivers of growth. The chapter illustrates this finding with two case studies: the New York City housing market after September 11, 2001; and the San Francisco Bay Area in the aftermath of the 1989 Loma Prieta earthquake. Both areas rebounded strongly after these shocks, suggesting the resilience of the urban metropolis.


title: Response to comment on "Lack of a genuine time crystal in a chiral soliton model" by \"Ohberg and Wright

abstract: In the paper [Phys. Rev. Research 2, 032038] we have analyzed a chiral soliton model and shown that despite the claim of \"Ohberg and Wright [Phys. Rev. Lett. 124, 178902], there is no indication that a genuine quantum time crystal can be observed in the system. Here, we response to the recent comment on our paper written by \"Ohberg and Wright.

categories: cond-mat.quant-gas cond-mat.mes-hall physics.atom-ph quant-ph


title: Dynamic Phase Alignment in Inertial Alfven Turbulence

abstract: In weakly-collisional plasma environments with sufficiently low electron beta, Alfv\'enic turbulence transforms into inertial Alfv\'enic turbulence at scales below the electron skin-depth, $k_\perp d_e\gtrsim 1$. We argue that, in inertial Alfv\'enic turbulence, both energy and generalized kinetic helicity exhibit direct cascades. We demonstrate that the two cascades are compatible due to the existence of a strong scale-dependence of the phase alignment angle between velocity and magnetic field fluctuations, with the phase alignment angle scaling as $\cos\alpha_k\propto k_{\perp}^{-1}$. The kinetic and magnetic energy spectra scale as $\propto k_{\perp}^{-5/3}$ and $\propto k_{\perp}^{-11/3}$, respectively. As a result of the dual direct cascade, the generalized-helicity spectrum scales as $\propto k^{-5/3}_{\perp}$, implying progressive balancing of the turbulence as the cascade proceeds to smaller scales in the $k_{\perp} d_e \gg 1$ range. Turbulent eddies exhibit a phase-space anisotropy $k_{\parallel} \propto k_{\perp}^{5/3}$, consistent with critically-balanced inertial Alfv\'en fluctuations. Our results may be applicable to a variety of geophysical, space, and astrophysical environments, including the Earth's magnetosheath and ionosphere, solar corona, non-relativistic pair plasmas, as well as to strongly rotating non-ionized fluids.


title: PL density invariant for type II degenerating K3 surfaces, Moduli compactification and hyperKahler metrics

abstract: A protagonist here is a new-type invariant for type II degenerations of K3 surfaces, which is explicit PL (piecewise linear) convex function from the interval with at most 18 non-linear points. Forgetting its actual function behaviour, it also classifies the type II degenerations into several combinatorial types, depending on the type of root lattices as appeared in classical examples. From differential geometric viewpoint, the function is obtained as the density function of the limit measure on the collapsing hyperKahler metrics to conjectural segments, as in the work of Honda-Sun-Zhang. On the way, we also reconstruct a moduli compactification of elliptic K3 surfaces in the works of Brunyate, Ascher-Bejleri, Alexeev-Brunyate-Engel in a more elementary manner, analyze the cusps more explicitly. We also interpret the glued hyperKahler fibration of Hein-Sun-Viaclovsky-Zhang as a special case from our viewpoint, discuss other cases, and possible relations with Landau-Ginzburg models in the mirror symmetry context.


title: Learning to be safe, in finite time

abstract: This paper aims to put forward the concept that learning to take safe actions in unknown environments, even with probability one guarantees, can be achieved without the need for an unbounded number of exploratory trials, provided that one is willing to relax its optimality requirements mildly. We focus on the canonical multi-armed bandit problem and seek to study the exploration-preservation trade-off intrinsic within safe learning. More precisely, by defining a handicap metric that counts the number of unsafe actions, we provide an algorithm for discarding unsafe machines (or actions), with probability one, that achieves constant handicap. Our algorithm is rooted in the classical sequential probability ratio test, redefined here for continuing tasks. Under standard assumptions on sufficient exploration, our rule provably detects all unsafe machines in an (expected) finite number of rounds. The analysis also unveils a trade-off between the number of rounds needed to secure the environment and the probability of discarding safe machines. Our decision rule can wrap around any other algorithm to optimize a specific auxiliary goal since it provides a safe environment to search for (approximately) optimal policies. Simulations corroborate our theoretical findings and further illustrate the aforementioned trade-offs.


title: Rigidity and Flexibility of Isometric Extensions

abstract: In this paper we consider the rigidity and flexibility of $C^{1, \theta}$ isometric extensions and we show that the H\"older exponent $\theta_0=\frac12$ is critical in the following sense: if $u\in C^{1,\theta}$ is an isometric extension of a smooth isometric embedding of a codimension one submanifold $\Sigma$ and $\theta> \frac12$, then the tangential connection agrees with the Levi-Civita connection along $\Sigma$. On the other hand, for any $\theta<\frac12$ we can construct $C^{1,\theta}$ isometric extensions via convex integration which violate such property. As a byproduct we get moreover an existence theorem for $C^{1, \theta}$ isometric embeddings, $\theta<\frac12$, of compact Riemannian manifolds with $C^1$ metrics and sharper amount of codimension.


title: Cutpoints of invariant subcontinua of polynomial Julia sets

abstract: We prove fixed point results for branched covering maps $f$ of the plane. For complex polynomials $P$ with Julia set $J_P$ these imply that periodic cutpoints of some invariant subcontinua of $J_P$ are also cutpoints of $J_P$. We deduce that, under certain assumptions on invariant subcontinua $Q$ of $J_P$, every Riemann ray to $Q$ landing at a periodic repelling/parabolic point $x\in Q$ is isotopic to a Riemann ray to $J_P$ relative to $Q$.


title: A quantum moat barrier, realized with a finite square well

abstract: The notion of a double well potential typically involves two regions of space separated by a repulsive potential barrier. The solution is a wave function that is suppressed in the barrier region and localized in the two surrounding regions. Remarkably, we illustrate that similar solutions can be achieved using an attractive "barrier" potential (a "quantum moat") instead of a repulsive one (a "quantum wall"). The reason this works is intimately connected to the concepts of "orthogonalized plane waves" and the pseudopotential method, both originally used to understand electronic band structures in solids. While the main goal of this work is to use a simple model to demonstrate the barrier-like attribute of a quantum moat, we also show how the pseudopotential method is used to greatly improve the efficiency of constructing wave functions for this system using matrix diagonalization.


title: A computationally efficient reconstruction algorithm for circular cone-beam computed tomography using shallow neural networks

abstract: Circular cone-beam (CCB) Computed Tomography (CT) has become an integral part of industrial quality control, materials science and medical imaging. The need to acquire and process each scan in a short time naturally leads to trade-offs between speed and reconstruction quality, creating a need for fast reconstruction algorithms capable of creating accurate reconstructions from limited data. In this paper we introduce the Neural Network Feldkamp-Davis-Kress (NN-FDK) algorithm. This algorithm adds a machine learning component to the FDK algorithm to improve its reconstruction accuracy while maintaining its computational efficiency. Moreover, the NN-FDK algorithm is designed such that it has low training data requirements and is fast to train. This ensures that the proposed algorithm can be used to improve image quality in high throughput CT scanning settings, where FDK is currently used to keep pace with the acquisition speed using readily available computational resources. We compare the NN-FDK algorithm to two standard CT reconstruction algorithms and to two popular deep neural networks trained to remove reconstruction artifacts from the 2D slices of an FDK reconstruction. We show that the NN-FDK reconstruction algorithm is substantially faster in computing a reconstruction than all the tested alternative methods except for the standard FDK algorithm and we show it can compute accurate CCB CT reconstructions in cases of high noise, a low number of projection angles or large cone angles. Moreover, we show that the training time of an NN-FDK network is orders of magnitude lower than the considered deep neural networks, with only a slight reduction in reconstruction accuracy.


title: Supercomputing with MPI meets the Common Workflow Language standards: an experience report

abstract: Use of standards-based workflows is still somewhat unusual by high-performance computing users. In this paper we describe the experience of using the Common Workflow Language (CWL) standards to describe the execution, in parallel, of MPI-parallelised applications. In particular, we motivate and describe the simple extension to the specification which was required, as well as our implementation of this within the CWL reference runner. We discuss some of the unexpected benefits, such as simple use of HPC-oriented performance measurement tools, and CWL software requirements interfacing with HPC module systems. We close with a request for comment from the community on how these features could be adopted within versions of the CWL standards.


title: Importance of long-range channel Sr displacements for the narrow emission in Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$ phosphor

abstract: The recently discovered Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$ red phosphor, candidate for the next generation of eco-efficient white light-emitting diodes, exhibits excellent emission spectral position and exceptionally small linewidth. It belongs to the UCr$_4$C$_4$-structure family of phosphors, containing many potential candidates commercial phosphors, whose small linewidth, tentatively ascribed to the high-symmetry cuboid environment of the doping site, has drawn the attention of researchers in the last five years. We use density functional theory, $\Delta$SCF method and configuration coordinate models (CCM) to provide a complete characterization of this material. Using a multi-dimensional CCM, an accurate description of the coupling of the vibronic structure with the electronic 5d$\rightarrow$4f transition is obtained, including the partial Huang-Rhys factors and frequency of the dominant modes. We show that, in addition to the first-coordination shell cuboid deformation mode, low-frequency phonon modes involving chains of strontium atoms along the tetragonal axis shape the emission linewidth in Sr[Li$_2$Al$_2$O$_2$N$_2$]:Eu$^{2+}$. This finding sheds new light on the emission properties of UCr$_4$C$_4$-structure phosphors, possessing similar Ca/Sr/Ba channel. Our approach provides a robust theoretical framework to systematically study the emission spectra of such Eu-doped phosphors, and predict candidates with expected similar or even sharper linewidth.


title: An invariance principle for gradient flows in the space of probability measures

abstract: We seek to establish qualitative convergence results to a general class of evolution PDEs described by gradient flows in optimal transportation distances. These qualitative convergence results come from dynamical systems under the general name of LaSalle Invariance Principle. By combining some of the basic notions of gradient flow theory and dynamical systems, we are able to reproduce this invariance principle in the setting of evolution PDEs under general assumptions. We apply this abstract theory to a non-exhaustive list of examples that recover, simplify, and even extend the results in their respective literatures.


title: Forced Variational Integrators for the Formation Control of Multi-Agent Systems

abstract: Formation control of autonomous agents can be seen as a physical system of individuals interacting with local potentials, and whose evolution can be described by a Lagrangian function. In this paper, we construct and implement forced variational integrators for the formation control of autonomous agents modeled by double integrators. In particular, we provide an accurate numerical integrator with a lower computational cost than traditional solutions. We find error estimations for the rate of the energy dissipated along with the agents' motion to achieve desired formations. Consequently, this permits to provide sufficient conditions on the simulation's time step for the convergence of discrete formation control systems such as the consensus problem in discrete systems. We present practical applications such as the rapid estimation of regions of attraction to desired shapes in distance-based formation control.


title: Elastically driven Kelvin-Helmholtz-like instability in planar channel flow

abstract: Kelvin-Helmholtz instability (KHI) is widely spread in nature on scales from micrometer up to Galactic one. This instability refers to the growth of perturbation of an interface between two parallel streams of Newtonian fluids with different velocities and densities, destabilized by shear strain and stabilized by density stratification with the heavier fluid at the bottom. Here, we report the discovery of the purely elastic KH-like instability in planar straight channel flow of viscoelastic fluid, which is theoretically considered to be stable. However, despite the remarkable similarity to the Newtonian KHI temporal interface dynamics, the elastic KHI reveals qualitatively different instability mechanism. Indeed, the velocity difference across the interface strongly fluctuates and non-monotonically varies in time due to energy pumping by elastic waves, detected in the flow. A correlation of the elastic wave intensity and efficiency of the elastic KHI in different regimes suggests that the competition between the destabilizing factor of the elastic waves and the stabilizing effect of the elastic stress difference generated by the velocity difference at the interface is the novel instability mechanism of the elastic KHI.


title: Proton path reconstruction for pCT using Neural Networks

abstract: The Most Likely Path formalism (MLP) is widely established as the most statistically precise method for proton path reconstruction in proton computed tomography (pCT). However, while this method accounts for small-angle Multiple Coulomb Scattering (MCS) and energy loss, inelastic nuclear interactions play an influential role in a significant number of proton paths. By applying cuts based on energy and direction, tracks influenced by nuclear interactions are largely discarded from the MLP analysis. In this work we propose a new method to estimate the proton paths based on a Deep Neural Network (DNN). Through this approach, estimates of proton paths equivalent to MLP predictions have been achieved in the case where only MCS occurs, together with an increased accuracy when nuclear interactions are present. Moreover, our tests indicate that the DNN algorithm can be considerably faster than the MLP algorithm.


title: A posteriori Error Estimates for Numerical Solutions to Hyperbolic Conservation Laws

abstract: The paper is concerned with a posteriori error bounds for a wide class of numerical schemes, for $n\times n$ hyperbolic conservation laws in one space dimension. These estimates are achieved by a "post-processing algorithm", checking that the numerical solution retains small total variation, and computing its oscillation on suitable subdomains. The results apply, in particular, to solutions obtained by the Godunov or the Lax-Friedrichs scheme, backward Euler approximations, and the method of periodic smoothing. Some numerical implementations are presented.


title: Injective edge-coloring of graphs with given maximum degree

abstract: A coloring of edges of a graph $G$ is injective if for any two distinct edges $e_1$ and $e_2$, the colors of $e_1$ and $e_2$ are distinct if they are at distance $1$ in $G$ or in a common triangle. Naturally, the injective chromatic index of $G$, $\chi'_{inj}(G)$, is the minimum number of colors needed for an injective edge-coloring of $G$. We study how large can be the injective chromatic index of $G$ in terms of maximum degree of $G$ when we have restrictions on girth and/or chromatic number of $G$. We also compare our bounds with analogous bounds on the strong chromatic index.


title: Circuit-QED with phase-biased Josephson weak links

abstract: By coupling a superconducting weak link to a microwave resonator, recent experiments probed the spectrum and achieved the quantum manipulation of Andreev states in various systems. However, the quantitative understanding of the response of the resonator to changes in the occupancy of the Andreev levels, which are of fermionic nature, is missing. Here, using Bogoliubov-de Gennes formalism to describe the weak link and a general formulation of the coupling to the resonator, we calculate the shift of the resonator frequency as a function of the levels occupancy and describe how transitions are induced by phase or electric field microwave drives. We apply this formalism to analyze recent experimental results obtained using circuit-QED techniques on superconducting atomic contacts and semiconducting nanowire Josephson junctions.


title: Confinement and the Global $SU(3)$ Color Symmetry

abstract: The global $SU(3)$ color symmetry and its physical consequences are discussed. The N\"{o}ther current is actually governed by the conserved matter current of color charges if the color field generated by this charge is properly polarized. The color field strength of a charge can have a uniform part due to the nontrivial QCD vacuum field and the nonzero gluon condensate, which implies that the self-energy of a system with a net color charge is infinite and thereby cannot exist as a free state. This is precisely what the color confinement means. Accordingly, the Cornell type potential with the feature of the Casimir scaling is derived for a color singlet system composed of a static color charge and an anti-charge. The uniform color field also implies that a hadron has a minimal size and a minimal energy. Furthermore, the global $SU(3)$ color symmetry requires that the minimal irreducible color singlet systems can only be $q\bar{q}$, $qqq$, $gg$, $ggg$, $q\bar{q}g$, $qqqg$ and $\bar{q}\bar{q}\bar{q}g$, etc., as such a multi-quark systems can only exist as a molecular configurations if there are no other binding mechanisms.


title: The RFML Ecosystem: A Look at the Unique Challenges of Applying Deep Learning to Radio Frequency Applications

abstract: While deep machine learning technologies are now pervasive in state-of-the-art image recognition and natural language processing applications, only in recent years have these technologies started to sufficiently mature in applications related to wireless communications. In particular, recent research has shown deep machine learning to be an enabling technology for cognitive radio applications as well as a useful tool for supplementing expertly defined algorithms for spectrum sensing applications such as signal detection, estimation, and classification (termed here as Radio Frequency Machine Learning, or RFML). A major driver for the usage of deep machine learning in the context of wireless communications is that little, to no, a priori knowledge of the intended spectral environment is required, given that there is an abundance of representative data to facilitate training and evaluation. However, in addition to this fundamental need for sufficient data, there are other key considerations, such as trust, security, and hardware/software issues, that must be taken into account before deploying deep machine learning systems in real-world wireless communication applications. This paper provides an overview and survey of prior work related to these major research considerations. In particular, we present their unique considerations in the RFML application space, which are not generally present in the image, audio, and/or text application spaces.


title: A note on the amount of information borrowed from external data in hybrid controlled trials with time-to-event outcomes

abstract: In situations where it is difficult to enroll patients in randomized controlled trials, external data can improve efficiency and feasibility. In such cases, adaptive trial designs could be used to decrease enrollment in the control arm of the trial by updating the randomization ratio at the interim analysis. Updating the randomization ratio requires an estimate of the amount of information effectively borrowed from external data, which is typically done with a linear approximation. However, this linear approximation is not always a reliable estimate, which could potentially lead to sub-optimal randomization ratio updates. In this note, we highlight this issue through simulations for exponential time-to-event outcomes, because in this simple setting there is an exact solution available for comparison. We also propose a potential generalization that could complement the linear approximation in more complex settings, discuss challenges for this generalization, and recommend best practices for computing and interpreting estimates of the effective number of events borrowed.


title: Medical imaging data structure extended to multiple modalities and anatomical regions

abstract: Brain Imaging Data Structure (BIDS) allows the user to organise brain imaging data into a clear and easy standard directory structure. BIDS is widely supported by the scientific community and is considered a powerful standard for management. The original BIDS is limited to images or data related to the brain. Medical Imaging Data Structure (MIDS) was therefore conceived with the objective of extending this methodology to other anatomical regions and other types of imaging systems in these areas.


title: Community detection, pattern recognition, and hypergraph-based learning: approaches using metric geometry and persistent homology

abstract: Hypergraph data appear and are hidden in many places in the modern age. They are data structure that can be used to model many real data examples since their structures contain information about higher order relations among data points. One of the main contributions of our paper is to introduce a new topological structure to hypergraph data which bears a resemblance to a usual metric space structure. Using this new topological space structure of hypergraph data, we propose several approaches to study community detection problem, detecting persistent features arising from homological structure of hypergraph data. Also based on the topological space structure of hypergraph data introduced in our paper, we introduce a modified nearest neighbors methods which is a generalization of the classical nearest neighbors methods from machine learning. Our modified nearest neighbors methods have an advantage of being very flexible and applicable even for discrete structures as in hypergraphs. We then apply our modified nearest neighbors methods to study sign prediction problem in hypegraph data constructed using our method.


title: On an extension of the generalized BGW tau-function

abstract: For an arbitrary solution to the Burgers--KdV hierarchy, we define the tau-tuple $(\tau_1,\tau_2)$ of the solution. We show that the product $\tau_1\tau_2$ admits Buryak's residue formula. Therefore, according to Alexandrov's theorem, $\tau_1\tau_2$ is a tau-function of the KP hierarchy. We then derive a formula for the affine coordinates for the point of the Sato Grassmannian corresponding to the tau-function $\tau_1\tau_2$ explicitly in terms of those for $\tau_1$. Applications to the analogous open extension of the generalized BGW tau-function and to the open partition function are given.


title: Interaction of electrons with acoustic phonons in AlN/GaN resonant tunnelling nanostructures at different temperatures

abstract: The theory of the interaction of electrons with acoustic phonons in multilayer nitride-based AlN/GaN nanostructures was developed for the first time at $T\geqslant 0$ using the method of finite-temperature Green's functions and Dyson equation. Components of the Hamiltonian describing the system of electrons with acoustic phonons and the magnitudes of the electron spectrum shifts due to the electron-phonon interaction were obtained. Dependences of the electronic spectrum levels and spectrum of the acoustic phonons were found depending on the position of the internal potential barrier in the studied nanostructure. The temperature shifts of the electronic spectrum and decay rates were calculated for various values of temperature $T$.


title: Analysis of KNN Density Estimation

abstract: We analyze the $\ell_1$ and $\ell_\infty$ convergence rates of k nearest neighbor density estimation method. Our analysis includes two different cases depending on whether the support set is bounded or not. In the first case, the probability density function has a bounded support and is bounded away from zero. We show that kNN density estimation is minimax optimal under both $\ell_1$ and $\ell_\infty$ criteria, if the support set is known. If the support set is unknown, then the convergence rate of $\ell_1$ error is not affected, while $\ell_\infty$ error does not converge. In the second case, the probability density function can approach zero and is smooth everywhere. Moreover, the Hessian is assumed to decay with the density values. For this case, our result shows that the $\ell_\infty$ error of kNN density estimation is nearly minimax optimal. The $\ell_1$ error does not reach the minimax lower bound, but is better than kernel density estimation.


title: Learning Set Functions that are Sparse in Non-Orthogonal Fourier Bases

abstract: Many applications of machine learning on discrete domains, such as learning preference functions in recommender systems or auctions, can be reduced to estimating a set function that is sparse in the Fourier domain. In this work, we present a new family of algorithms for learning Fourier-sparse set functions. They require at most $nk - k \log_2 k + k$ queries (set function evaluations), under mild conditions on the Fourier coefficients, where $n$ is the size of the ground set and $k$ the number of non-zero Fourier coefficients. In contrast to other work that focused on the orthogonal Walsh-Hadamard transform, our novel algorithms operate with recently introduced non-orthogonal Fourier transforms that offer different notions of Fourier-sparsity. These naturally arise when modeling, e.g., sets of items forming substitutes and complements. We demonstrate effectiveness on several real-world applications.


title: Leptonic CP Violation and its Origin

abstract: This PhD thesis deals with the current experimental hint for large CP violation in the leptonic sector. We first quantify its global statistical significance, by combining the results from relevant neutrino oscillation experiments. As the hint is indirect, other new physics apart from three massive neutrinos could be present, masking the results as direct leptonic CP violation has not yet been conclusively observed. Thus, we move on to checking the robustness of the hint against the framework in which the experimental data is interpreted; in other words, its physical origin. To that purpose, we confront with data the scenario that is less bounded by other experiments: neutral current Non-Standard Interactions (NSI) among neutrinos and matter. In the end, CP violation induced by neutrino masses and leptonic mixing is found to be quite robust, due to the large amount of neutrino flavour transition data collected along three decades. Nevertheless, the next generation LBL accelerator experiments are aimed at precision measurements that could be more severely affected. Luckily, in the last years the COHERENT experiment has provided independent constraints on NSI. In the final chapter of the thesis, the data from the COHERENT experiment is analysed and integrated into the global analyses from previous chapters. These first results could be greatly improved by increasing the statistics of the signal and/or by performing measurements with different nuclei sensitive to different NSI models. For that, the European Spallation Source is an ideal future facility, whose prospects for bounding NSI are also explored. In summary, this thesis takes a global approach to rigorously assess whether cutting-edge leptonic flavour measurements are pointing towards a new strong violation of a symmetry of Nature.


title: Regularity of the optimal sets for the second Dirichlet eigenvalue

abstract: This paper is dedicated to the regularity of the optimal sets for the second eigenvalue of the Dirichlet Laplacian. Precisely, we prove that if the set $\Omega$ minimizes the functional \[ \mathcal F_\Lambda(\Omega)=\lambda_2(\Omega)+\Lambda |\Omega|, \] among all subsets of a smooth bounded open set $D\subset \mathbb{R}^d$, where $\lambda_2(\Omega)$ is the second eigenvalue of the Dirichlet Laplacian on $\Omega$ and $\Lambda>0$ is a fixed constant, then $\Omega$ is equivalent to the union of two disjoint open sets $\Omega_+$ and $\Omega_-$, which are $C^{1,\alpha}$-regular up to a (possibly empty) closed set of Hausdorff dimension at most $d-5$, contained in the one-phase free boundaries $D\cap \partial\Omega_+\setminus\partial\Omega_-$ and $D\cap\partial\Omega_-\setminus\partial\Omega_+$.


title: On the apparent molar volume of methanol in water-methanol mixtures. Composition and temperature effects from molecular dynamics study

abstract: We explore the effects of composition and temperature on the apparent molar volumes of species of water-methanol mixtures. Isothermal-isobaric molecular dynamics simulations are used with this purpose. Several combinations of models for water and for methanol are explored. Validity of predictions concerned with a puzzling minimum of apparent molar volume of methanol in water-rich solutions is tested against experimental results.


abstract: A relation between $\frac{1}{2}$-derivations of Lie algebras and transposed Poisson algebras was established. Some non-trivial transposed Poisson algebras with a certain Lie algebra (Witt algebra, algebra $\mathcal{W}(a,-1)$, thin Lie algebra and solvable Lie algebra with abelian nilpotent radical) were constructed. In particular, we constructed an example of the transposed Poisson algebra with associative and Lie parts isomorphic to the Laurent polynomials and the Witt algebra. On the other side, it was proven that there are no non-trivial transposed Poisson algebras with Lie algebra part isomorphic to a semisimple finite-dimensional algebra, a simple finite-dimensional superalgebra, the Virasoro algebra, $N=1$ and $N=2$ superconformal algebras, or a semisimple finite-dimensional $n$-Lie algebra.


title: Study of the Kondo and high-temperature limits of the slave-boson and X-boson methods

abstract: In this Letter we study the periodic Anderson model, employing both the slave-boson and the X-boson approaches in the mean field approximation. We investigate the breakdown of the slave-boson at intermediate temperatures when the total occupation number of particles Nt = Nf + Nc is keep constant, where Nf and Nc are respectively the occupation numbers of the localized and conduction electrons, and we show that the high-temperature limit of the slave-boson is Nf = Nc = Nt /2. We also compare the results of the two approaches in the Kondo limit and we show that at low-temperatures the X-boson exhibits a phase transition, from the Kondo heavy fermion (K-HF) regime to a local moment magnetic regime (LMM).


title: Superconductor-insulator transition induced by pressure within the X-boson approach

abstract: The pressure induced superconducting phase diagram is calculated for an extension of the periodic Anderson model (PAM) in the $ U = \infty $ limit taking into account the effect of a nearest neighbor attractive interaction between f-electrons. We analyze the role of the chemical potential compared to several plots of the f-band density of states and we also found a superconductor-insulator transition induced by pressure when the chemical potential crosses the hybridization gap.


title: Overcharging dilaton black holes in 2+1 dimensions to extremality and beyond

abstract: We test whether static charged dilaton black holes in 2+1 dimensions can be turned into naked singularities, by sending in test particles from infinity. We derive that overcharging is possible and generic for both extremal and nearly extremal black holes. Our analysis also imply that nearly extremal charged dilaton black holes can be continuously driven to extremality and beyond, unlike nearly extremal Ba\u{n}ados-Teitelboim-Zanelli, Kerr and Reissner-Nordstr\"{o}m black holes which are overspun or overcharged by a discrete jump. Thus the weak form of the cosmic censorship conjecture and the third law of black hole thermodynamics are both violated in the interaction of charged dilaton black holes in 2+1 dimensions, with test particles. We also derive that there exists no points where the heat capacity vanishes or diverges in the transition from black holes to naked singularities. The phase transitions that could potentially prevent the formation of naked singularities do not occur.


title: Dynamical fission of the quasiprojectile and isospin equilibration for the system 80Kr+ 48 Ca at 35 MeV/nucleon

abstract: Experimental results concerning the dynamical fission of quasiprojectiles in semiperipheral collisions for the system 80 Kr+ 48 Ca at 35 MeV/nucleon are presented. Data have been collected with four blocks of the FAZIA setup in the first physics experiment of the FAZIA Collaboration. The degree of isospin equilibration between the two fission fragments and its dependence on their charge asymmetry is investigated. The data are compared with the results of the AMD model coupled to GEMINI as an afterburner, in order to get hints about the timescale of the process.


title: Functions of perturbed commuting dissipative operators

abstract: The main objective of the paper is to obtain sharp Lipschitz type estimates for the norm of operator differences $f(L_1,M_1)-f(L_2,M_2)$ for pairs $(L_1,M_1)$ and $(L_2,M_2)$ of commuting maximal dissipative operators. To obtain such estimates, we use double operator integrals with respect to semi-spectral measures associated with the pairs $(L_1,M_1)$ and $(L_2,M_2)$. Note that the situation is considerably more complicated than in the case of functions of two commuting contractions and to overcome difficulties we had to elaborate new techniques. We deduce from the main result H\"older type estimates for operator differences as well as their estimates in Schatten--von Neumann norms.


title: On Parameter Selection of Nonsingular Predefined-time Terminal Sliding Mode With the Fixed-time Convergence Guarantee

abstract: This paper study the parameter selection of predefined-time sliding mode and try to design a general nonsingular predefined-time terminal sliding mode. 1). On parameter selection: Some existing predefined-time sliding modes are designed to focus on the reaching time and ignore the characterization of the equilibrium point i.e. x_e= 0. In actual engineer, the system can only converge to a small neighborhood near the equilibrium point because of the existence of uncertainty i.e. x_e \to 0. The actual equilibrium point should be deduced by taking the limit but not directly solve the equilibrium point. Hence, the actual selection of exponential term is suggested to be not 1 to make system convergence to the equilibrium point within predefined time. 2). On singularity-avoidance: Based on the exponential feature of predefined-time stability systems, a mathematic concept of switching sliding mode \cite{su1994adaptive} is applied in the design of nonsingular predefined-time terminal sliding mode. However, this class of sliding mode switching method will cause a singular problem, which can be described as control input u \to \infty if the system state x \to 0 and any-order time derivative of x is not 0. Hence, a non-singular selection condition for each recursive sliding mode surface is explored in this paper to realize the global finiteness of control input. In addition, an unnecessary chattering of control input will also be caused by the sliding mode switching frequency, and a novel switching condition will be applied in the controller to reduce this chattering. Simulations will be carried out to validate the effect of the novel nonsingular predefined-time sliding mode control method to some extent.


title: X-Fields: Implicit Neural View-, Light- and Time-Image Interpolation

abstract: We suggest to represent an X-Field -a set of 2D images taken across different view, time or illumination conditions, i.e., video, light field, reflectance fields or combinations thereof-by learning a neural network (NN) to map their view, time or light coordinates to 2D images. Executing this NN at new coordinates results in joint view, time or light interpolation. The key idea to make this workable is a NN that already knows the "basic tricks" of graphics (lighting, 3D projection, occlusion) in a hard-coded and differentiable form. The NN represents the input to that rendering as an implicit map, that for any view, time, or light coordinate and for any pixel can quantify how it will move if view, time or light coordinates change (Jacobian of pixel position with respect to view, time, illumination, etc.). Our X-Field representation is trained for one scene within minutes, leading to a compact set of trainable parameters and hence real-time navigation in view, time and illumination.


title: Gravitational Waves in Higher Order Teleparallel Gravity

abstract: The teleparallel equivalent of higher order Lagrangians like $L_{\Box R}=-R+a_{0}R^{2}+a_{1}R\Box R$ can be obtained by means of the boundary term $B=2\nabla_{\mu}T^{\mu}$. In this perspective, we derive the field equations in presence of matter for higher-order teleparallel gravity considering, in particular, sixth-order theories where the $\Box$ operator is linearly included. In the weak field approximation, gravitational wave solutions for these theories are derived. Three states of polarization are found: the two standard $+$ and $\times$ polarizations, namely 2-helicity massless transverse tensor polarizations, and a 0-helicity massive, with partly transverse and partly longitudinal scalar polarization. Moreover, these gravitational waves exhibit four oscillation modes related to four degrees of freedom: the two classical $+$ and $\times$ tensor modes of frequency $\omega_{1}$, related to the standard Einstein waves with $k^{2}_{1}=0$; two mixed longitudinal-transverse scalar modes for each frequencies $\omega_{2}$ and $\omega_{3}$, related to two different 4-wave vectors, $k^{2}_{2}=M_{2}^{2}$ and $k^{2}_{3}=M^{2}_{3}$. The four degrees of freedom are the amplitudes of each individual mode, i.e. $\hat{\epsilon}^{(+)}\left(\omega_{1}\right)$, $\hat{\epsilon}^{(\times)}\left(\omega_{1}\right)$, $\hat{B}_{2}\left(\mathbf{k}\right)$, and $\hat{B}_{3}\left(\mathbf{k}\right)$.


title: The Extreme Space Weather Event in February/March 1941

abstract: Given the infrequency of extreme geomagnetic storms, it is significant to note the concentration of three extreme geomagnetic storms in 1941, whose intensities ranked fourth, twelfth, and fifth within the aa index between 1868-2010. Among them, the geomagnetic storm on 1 March 1941 was so intense that three of the four Dst station magnetograms went off scale. Herein, we reconstruct its time series and measure the storm intensity with an alternative Dst estimate (Dst*). The source solar eruption at 09:29 - 09:38 GMT on 28 February was located at RGO AR 13814 and its significant intensity is confirmed by large magnetic crochets of 35 nT measured at Abinger. This solar eruption most likely released a fast interplanetary coronal mass ejection with estimated speed 2260 km/s. After its impact at 03:57 - 03:59 GMT on 1 March, an extreme magnetic storm was recorded worldwide. Comparative analyses on the contemporary magnetograms show the storm peak intensity of minimum Dst* < -464 nT at 16 GMT, comparable to the most and the second most extreme magnetic storms within the standard Dst index since 1957. This storm triggered significant low-latitude aurorae in the East Asian sector and their equatorward boundary has been reconstructed as 38.5{\deg} in invariant latitude. This result agrees with British magnetograms which indicate auroral oval moving above Abinger at 53.0{\deg} in magnetic latitude. The storm amplitude was even more enhanced in equatorial stations and consequently casts caveats on their usage for measurements of the storm intensity in Dst estimates.


title: Water is a radiation protection agent for ionised pyrrole

abstract: Radiation-induced damage of biological matter is an ubiquitous problem in nature. The influence of the hydration environment is widely discussed, but its exact role remains elusive. Utilising well defined solvated-molecule aggregates, we experimentally observed a hydrogen-bonded water molecule acting as a radiation protection agent for ionised pyrrole, a prototypical aromatic biomolecule. Pure samples of pyrrole and pyrrole(H$_2$O) were outer-valence ionised and the subsequent damage and relaxation processes were studied. Bare pyrrole ions fragmented through the breaking of C-C or N-C covalent bonds. However, for pyrrole(H$_2$O)$^+$, we observed a strong protection of the pyrrole ring through the dissociative release of neutral water or by transferring an electron or proton across the hydrogen bond. Overall, a single water molecule strongly reduces the fragmentation probability and thus the persistent radiation damage of singly-ionised pyrrole.


title: Evaluating Multilingual BERT for Estonian

abstract: Recently, large pre-trained language models, such as BERT, have reached state-of-the-art performance in many natural language processing tasks, but for many languages, including Estonian, BERT models are not yet available. However, there exist several multilingual BERT models that can handle multiple languages simultaneously and that have been trained also on Estonian data. In this paper, we evaluate four multilingual models -- multilingual BERT, multilingual distilled BERT, XLM and XLM-RoBERTa -- on several NLP tasks including POS and morphological tagging, NER and text classification. Our aim is to establish a comparison between these multilingual BERT models and the existing baseline neural models for these tasks. Our results show that multilingual BERT models can generalise well on different Estonian NLP tasks outperforming all baselines models for POS and morphological tagging and text classification, and reaching the comparable level with the best baseline for NER, with XLM-RoBERTa achieving the highest results compared with other multilingual models.


title: On the theta representations of finite inverse monoids

abstract: (I) We study Clifford-Mackey-Rieffel's theory for finite monoid; (II) We prove some results of Theta Representations of finite inverse monoids.


title: From Handcrafted to Deep Features for Pedestrian Detection: A Survey

abstract: Pedestrian detection is an important but challenging problem in computer vision, especially in human-centric tasks. Over the past decade, significant improvement has been witnessed with the help of handcrafted features and deep features. Here we present a comprehensive survey on recent advances in pedestrian detection. First, we provide a detailed review of single-spectral pedestrian detection that includes handcrafted features based methods and deep features based approaches. For handcrafted features based methods, we present an extensive review of approaches and find that handcrafted features with large freedom degrees in shape and space have better performance. In the case of deep features based approaches, we split them into pure CNN based methods and those employing both handcrafted and CNN based features. We give the statistical analysis and tendency of these methods, where feature enhanced, part-aware, and post-processing methods have attracted main attention. In addition to single-spectral pedestrian detection, we also review multi-spectral pedestrian detection, which provides more robust features for illumination variance. Furthermore, we introduce some related datasets and evaluation metrics, and compare some representative methods. We conclude this survey by emphasizing open problems that need to be addressed and highlighting various future directions. Researchers can track an up-to-date list at https://github.com/JialeCao001/PedSurvey.


title: A new ionisation network and radiation transport module in PLUTO

abstract: We introduce a new general-purpose time-dependent ionisation network (IN) and a radiation transport (RT) module in the magneto-hydrodynamic (MHD) code PLUTO. Our ionisation network is reliable for temperatures ranging from 5e3 to 3e8 K, and includes all ionisation states of H, He, C, N, O, Ne, Mg, Si, S and Fe making it suitable for studying a variety of astrophysical scenarios. Radiation loss for each ion-electron pair is calculated using CLOUDY-17 data on-the-fly. Photo-ionisation and charge exchange are the chemical heating mechanisms. The IN is fully coupled to the radiation transport module over a very large range of opacities at different frequencies. The RT module employs a method of short characteristics assuming spherical symmetry. The radiation module requires the assumption of spherical symmetry, while the IN is compatible with full 3D. We also include a simple prescription for dust opacity, grain destruction, and the dust contribution to radiation pressure. We present numerical tests to show the reliability and limitations of the new modules. We also present a post-processing tool to calculate projected column densities and emission spectra.


title: Characters and chromatic symmetric functions

abstract: Let $P$ be a poset, $inc(P)$ its incomparability graph, and $X_{inc(P)}$ the corresponding chromatic symmetric function, as defined by Stanley in {\em Adv. Math.}, {\bf 111} (1995) pp.~166--194. Certain conditions on $P$ imply that the expansions of $X_{inc(P)}$ in standard symmetric function bases yield coefficients which have simple combinatorial interpretations. By expressing these coefficients as character evaluations, we extend several of these interpretations to {\em all} posets $P$. Consequences include new combinatorial interpretations of the permanent and other immanants of totally nonnegative matrices, and of the sum of elementary coefficients in the Shareshian-Wachs chromatic quasisymmetric function $X_{inc(P),q}$ when $P$ is a unit interval order.


title: A baseline Model for the Relationships between Network Operators and Tower Companies

abstract: The introduction of virtualization techniques in radio cellular networks allows the emergence of a business based on the outsourcing of towers hosting antennas and operated by the so-called Tower Companies (TowerCos). In this paper, we develop a baseline business model for studying the potential relationships between network operators and TowerCos. It turns out that the gain in operational costs achieved when network operators outsource the management of towers can be gracefully utilized to reduce prices so as to attract more customers. The price drop has however to be carefully realized so as not to break the market share between operators and to preserve competition. To prove this claim, we adopt in a first step a centralized optimization formulation. In a second step, we develop a game theoretic framework.


title: The quaternionic Goos-Haenchen shift

abstract: We investigate the lateral displacement of electronic waves when partially and totally reflected by a quaternionic potential. Following the analogy between Quantum Mechanics and Optics we introduce a refractive index for the complex and the pure quaternionic case. For incidence greater than the critical one (total reflection) the quaternionic potentials amplify the lateral displacement found in the complex case. For incidence below the critical angle, we find, in the quaternionic case, an additional shift. The analytical formula found for the quaternionic lateral displacement gives the possibility to observe quantitative and qualitative differences between the complex and quaternionic Goos-Haenchen shift.


title: Radio telescope total power mode: improving observation efficiency

abstract: Aims. Radio observing efficiency can be improved by calibrating and reducing the observations in total power mode rather than in frequency, beam, or position-switching modes. Methods. We selected a sample of spectra obtained from the Institut de Radio-Astronomie Millim\'etrique (IRAM) 30-m telescope and the Green Bank Telescope (GBT) to test the feasibility of the method. Given that modern front-end amplifiers for the GBT and direct Local Oscillator injection for the 30 m telescope provide smooth pass bands that are a few tens of megahertz in width, the spectra from standard observations can be cleaned (baseline removal) separately and then co-added directly when the lines are narrow enough (a few km/s), instead of performing the traditional ON minus OFF data reduction. This technique works for frequency-switched observations as well as for position- and beam-switched observations when the ON and OFF data are saved separately. Results. The method works best when the lines are narrow enough and not too numerous so that a secure baseline removal can be achieved. A signal-to-noise ratio improvement of a factor of sqrt(2) is found in most cases, consistent with theoretical expectations. Conclusions. By keeping the traditional observing mode, the fallback solution of the standard reduction technique is still available in cases of suboptimal baseline behavior, sky instability, or wide lines, and to confirm the line intensities. These techniques of total-power-mode reduction can be applied to any radio telescope with stable baselines as long as they record and deliver the ONs and OFFs separately, as is the case for the GBT.


title: A survey on natural language processing (nlp) and applications in insurance

abstract: Text is the most widely used means of communication today. This data is abundant but nevertheless complex to exploit within algorithms. For years, scientists have been trying to implement different techniques that enable computers to replicate some mechanisms of human reading. During the past five years, research disrupted the capacity of the algorithms to unleash the value of text data. It brings today, many opportunities for the insurance industry.Understanding those methods and, above all, knowing how to apply them is a major challenge and key to unleash the value of text data that have been stored for many years. Processing language with computer brings many new opportunities especially in the insurance sector where reports are central in the information used by insurers. SCOR's Data Analytics team has been working on the implementation of innovative tools or products that enable the use of the latest research on text analysis. Understanding text mining techniques in insurance enhances the monitoring of the underwritten risks and many processes that finally benefit policyholders.This article proposes to explain opportunities that Natural Language Processing (NLP) are providing to insurance. It details different methods used today in practice traces back the story of them. We also illustrate the implementation of certain methods using open source libraries and python codes that we have developed to facilitate the use of these techniques.After giving a general overview on the evolution of text mining during the past few years,we share about how to conduct a full study with text mining and share some examples to serve those models into insurance products or services. Finally, we explained in more details every step that composes a Natural Language Processing study to ensure the reader can have a deep understanding on the implementation.


title: A Finite Memory Interacting P\'{o}lya Contagion Network and its Approximating Dynamical Systems

abstract: We introduce a new model for contagion spread using a network of interacting finite memory two-color P\'{o}lya urns, which we refer to as the finite memory interacting P\'{o}lya contagion network. The urns interact in the sense that the probability of drawing a red ball (which represents an infection state) for a given urn, not only depends on the ratio of red balls in that urn but also on the ratio of red balls in the other urns in the network, hence accounting for the effect of spatial contagion. The resulting network-wide contagion process is a discrete-time finite-memory ($M$th order) Markov process, whose transition probability matrix is determined. The stochastic properties of the network contagion Markov process are analytically examined, and for homogeneous system parameters, we characterize the limiting state of infection in each urn. For the non-homogeneous case, given the complexity of the stochastic process, and in the same spirit as the well-studied SIS models, we use a mean-field type approximation to obtain a discrete-time dynamical system for the finite memory interacting P\'{o}lya contagion network. Interestingly, for $M=1$, we obtain a linear dynamical system which exactly represents the corresponding Markov process. For $M>1$, we use mean-field approximation to obtain a nonlinear dynamical system. Furthermore, noting that the latter dynamical system admits a linear variant (realized by retaining its leading linear terms), we study the asymptotic behavior of the linear systems for both memory modes and characterize their equilibrium. Finally, we present simulation studies to assess the quality of the approximation purveyed by the linear and non-linear dynamical systems.


title: Reactive atomistic simulations of Diels-Alder-type reactions: Conformational and dynamic effects in the polar cycloaddition of 2,3-dibromobutadiene radical ions with maleic anhydride

abstract: The kinetics, dynamics and conformational specificities for the ionic Diels-Alder reaction (polar cycloaddition) of maleic anhydride with 2,3-dibromobutadiene radical ions have been studied theoretically using multisurface adiabatic reactive molecular dynamics. A competition of concerted and stepwise reaction pathways was found and both the s-cis and s-trans conformers of of the diene are reactive. The analysis of the minimum dynamic path of the reaction indicates that both, rotations and vibrations of the reactant molecules are important for driving the system towards the transition state. The rates were computed as $k = 5.1 \times 10^{-14}$ s$^{-1}$ for the s-cis and $k = 3.8 \times 10^{-14}$ s$^{-1}$ for the s-trans conformer of 2,3-dibromobutadiene at an internal temperature of 300 K. The present results are to be contrasted with the neutral variant of the title system in which only the gauche conformer of the diene was found to undergo a considerably slower, concerted and mostly synchronous reaction driven by the excitation of rotations. The results presented here inform detailed experimental studies of the dynamics of polar cycloadditions under single-collision conditions in the gas phase.


title: Computing the matrix sine and cosine simultaneously with a reduced number of products

abstract: A new procedure is presented for computing the matrix cosine and sine simultaneously by means of Taylor polynomial approximations. These are factorized so as to reduce the number of matrix products involved. Two versions are developed to be used in single and double precision arithmetic. The resulting algorithms are more efficient than schemes based on Pad\'e approximations for a wide range of norm matrices.