title: Digital Contact Tracing: Large-scale Geolocation Data as an Alternative to Bluetooth-based Apps' Failure

abstract: The currently deployed contact-tracing mobile apps have failed as an efficient solution in the context of the COVID-19 pandemic. None of them has managed to attract the number of active users required to achieve an efficient operation. This urges the research community to re-open the debate and explore new avenues that lead to efficient contact-tracing solutions. This paper contributes to this debate with an alternative contact-tracing solution that leverages already available geolocation information owned by BigTech companies with very large penetration rates in most countries adopting contact-tracing mobile apps. Moreover, our solution provides sufficient privacy guarantees to protect the identity of infected users as well as precluding Health Authorities from obtaining the contact graph from individuals.


title: Large Deviation Principles for Block and Step Graphon Random Graph Models

abstract: Borgs, Chayes, Gaudio, Petti and Sen [arXiv:2007.14508] proved a large deviation principle for block model random graphs with rational block ratios. We strengthen their result by allowing any block ratios (and also establish a simpler formula for the rate function). We apply the new result to derive a large deviation principle for graph sampling from any given step graphon.


title: Time-Efficient and High-Quality Graph Partitioning for Graph Dynamic Scaling

abstract: The dynamic scaling of distributed computations plays an important role in the utilization of elastic computational resources, such as the cloud. It enables the provisioning and de-provisioning of resources to match dynamic resource availability and demands. In the case of distributed graph processing, changing the number of the graph partitions while maintaining high partitioning quality imposes serious computational overheads as typically a time-consuming graph partitioning algorithm needs to execute each time repartitioning is required. In this paper, we propose a dynamic scaling method that can efficiently change the number of graph partitions while keeping its quality high. Our idea is based on two techniques: preprocessing and very fast edge partitioning, called graph edge ordering and chunk-based edge partitioning, respectively. The former converts the graph data into an ordered edge list in such a way that edges with high locality are closer to each other. The latter immediately divides the ordered edge list into an arbitrary number of high-quality partitions. The evaluation with the real-world billion-scale graphs demonstrates that our proposed approach significantly reduces the repartitioning time, while the partitioning quality it achieves is on par with that of the best existing static method.


title: Semiclassical analysis on compact nilmanifolds

abstract: In this paper, we develop a semiclassical calculus on compact nil-manifolds. As an application of the former, we obtain asymptotics such as generalised Weyl laws for positive Rockland operators on any graded compact nil-manifolds. We also define and study semiclassical limits and quantum limits in this context.


title: Abnormal topological refraction into free medium at sub-wavelength scale in valley phononic crystal plates

abstract: In this work we propose a topological valley phononic crystal plate and we extensively investigate the refraction of valley modes into the surrounding homogeneous medium. This phononic crystal includes two sublattices of resonators (A and B) modeled by mass-spring systems. We show that two edge states confined at the AB/BA and BA/AB type domain walls exhibit different symmetries in physical space and energy peaks in the Fourier space. As a result, distinct refraction behaviors, especially through an armchair cut edge, are observed. On the other hand, the decay depth of these localized topological modes, which is found to be solely determined by the relative resonant strength between the scatterers, significantly affects the refraction patterns. More interestingly, the outgoing traveling wave through a zigzag interface becomes evanescent when operating at deep sub-wavelength scale. This is realized by tuning the average resonant strength. We show that the evanescent modes only exist along a particular type of outlet edge, and that they can couple with both topological interface states. We also present two designs of topological functional devices, including an elastic one-way transmission waveguide and a near-ideal monopole/dipole emitter, both based on our phononic structure.


title: A deterministic matching method for exact matchings to compare the outcome of different interventions

abstract: Statistical matching methods are widely used in the social and health sciences to estimate causal effects using observational data. Often the objective is to find comparable groups with similar covariate distributions in a dataset, with the aim to reduce bias in a random experiment. We aim to develop a foundation for deterministic methods which provide results with low bias, while retaining interpretability. The proposed method matches on the covariates and calculates all possible maximal exact matchesfor a given dataset without adding numerical errors. Notable advantages of our method over existing matching algorithms are that all available information for exact matches is used, no additional bias is introduced, it can be combined with other matching methods for inexact matching to reduce pruning and that the result is calculated in a fast and deterministic way. For a given dataset the result is therefore provably unique for exact matches in the mathematical sense. We provide proofs, instructions for implementation as well as a numerical example calculated for comparison on a complete survey.


title: X-ray pulsar XTE J1858+034: discovery of the cyclotron line and the revised optical identification

abstract: We present results of a detailed investigation of the poorly studied X-ray pulsar XTE J1858+034 based on the data obtained with the NuSTAR observatory during the outburst of the source in 2019. The spectral analysis resulted in the discovery of a cyclotron absorption feature in the source spectrum at ~48 keV both in the pulse phase averaged and resolved spectra. Accurate X-ray localization of the source using the NuSTAR and Chandra observatories allowed us to accurately determine the position of the X-ray source and identify the optical companion of the pulsar. The analysis of the counterpart properties suggested that the system is likely a symbiotic binary hosting an X-ray pulsar and a late type companion star of K-M classes rather than Be X-ray binary as previously suggested.


title: Blaschke-Santal\'o inequalities for Minkowski and Asplund endomorphisms

abstract: It is shown that each monotone Minkowski endomorphism of convex bodies gives rise to an isoperimetric inequality which directly implies the classical Urysohn inequality. Among this large family of new inequalities, the only affine invariant one - the Blaschke-Santal\'o inequality - turns out to be the strongest one. A further extension of these inequalities to merely weakly monotone Minkowski endomorphisms is proven to be impossible. Moreover, for functional analogues of monotone Minkowski endomorphisms, a family of analytic inequalities for log-concave functions is established which generalizes the functional Blaschke-Santal\'o inequality.


title: Federated Learning Based Proactive Handover in Millimeter-wave Vehicular Networks

abstract: Proactive handover can avoid frequent handovers and reduce handover delay, which plays an important role in maintaining the quality of service (QoS) for mobile users in millimeter-wave vehicular networks. To reduce the communication cost of training the learning model for proactive handover, we propose a federated learning (FL) framework. The proposed FL framework can accommodate the limited storage capacity of each user, increase the number of users who participate in the FL, and adapt to the dynamic mobility pattern. Simulation results validate the effectiveness of the proposed FL framework. Compared to reactive handover schemes, the proposed handover scheme can reduce unnecessary handovers and improve the QoS of users simultaneously.


title: Analysis of key flavors of event-driven predictive maintenance using logs of phenomena described by Weibull distributions

abstract: This work explores two approaches to event-driven predictive maintenance in Industry 4.0 that cast the problem at hand as a classification or a regression one, respectively, using as a starting point two state-of-the-art solutions. For each of the two approaches, we examine different data preprocessing techniques, different prediction algorithms and the impact of ensemble and sampling methods. Through systematic experiments regarding the aspectsmentioned above,we aimto understand the strengths of the alternatives, and more importantly, shed light on how to navigate through the vast number of such alternatives in an informed manner. Our work constitutes a key step towards understanding the true potential of this type of data-driven predictive maintenance as of to date, and assist practitioners in focusing on the aspects that have the greatest impact.


title: AGRNet: Adaptive Graph Representation Learning and Reasoning for Face Parsing

abstract: Face parsing infers a pixel-wise label to each facial component, which has drawn much attention recently. Previous methods have shown their success in face parsing, which however overlook the correlation among facial components. As a matter of fact, the component-wise relationship is a critical clue in discriminating ambiguous pixels in facial area. To address this issue, we propose adaptive graph representation learning and reasoning over facial components, aiming to learn representative vertices that describe each component, exploit the component-wise relationship and thereby produce accurate parsing results against ambiguity. In particular, we devise an adaptive and differentiable graph abstraction method to represent the components on a graph via pixel-to-vertex projection under the initial condition of a predicted parsing map, where pixel features within a certain facial region are aggregated onto a vertex. Further, we explicitly incorporate the image edge as a prior in the model, which helps to discriminate edge and non-edge pixels during the projection, thus leading to refined parsing results along the edges. Then, our model learns and reasons over the relations among components by propagating information across vertices on the graph. Finally, the refined vertex features are projected back to pixel grids for the prediction of the final parsing map. To train our model, we propose a discriminative loss to penalize small distances between vertices in the feature space, which leads to distinct vertices with strong semantics. Experimental results show the superior performance of the proposed model on multiple face parsing datasets, along with the validation on the human parsing task to demonstrate the generalizability of our model.


title: Improved Spin-State Energy Differences of Fe(II) molecular and crystalline complexes via the Hubbard U-corrected Density

abstract: We recently showed that the DFT+U approach with a linear-response U yields adiabatic energy differences biased towards high spin [Mariano et al. J. Chem. Theory Comput. 2020, 16, 6755-6762]. Such bias is removed here by employing a density-corrected DFT approach where the PBE functional is evaluated on the Hubbard U -corrected density. The adiabatic energy differences of six Fe(II) molecular complexes computed using this approach, named here PBE[U], are in excellent agreement with coupled luster-corrected CASPT2 values for both weak- and strong-field ligands resulting in a mean absolute error (MAE) of 0.44 eV, smaller than the recently proposed Hartree-Fock density-corrected DFT (1.22 eV) and any other tested functional, including the best performer TPSSh (0.49 eV). We take advantage of the computational efficiency of this approach and compute the adiabatic energy differences of five molecular crystals using PBE[U] with periodic boundary conditions. The results show, again, an excellent agreement (MAE=0.07 eV) with experimentally-extracted values and a superior performance compared with the best performers TPSSh (MAE=0.08 eV) and M06-L (MAE=0.31 eV) computed on molecular fragments.


title: Iterative Facial Image Inpainting Based on an Encoder-Generator Architecture

abstract: Facial image inpainting is a challenging problem as it requires generating new pixels that include semantic information for masked key components in a face, e.g., eyes and nose. Recently, remarkable methods have been proposed in this field. Most of these approaches use encoder-decoder architectures and have different limitations such as allowing unique results for a given image and a particular mask. Alternatively, some optimization-based approaches generate promising results using different masks with generator networks. However, these approaches are computationally more expensive. In this paper, we propose an efficient solution to the facial image inpainting problem using the Cyclic Reverse Generator (CRG) architecture, which provides an encoder-generator model. We use the encoder to embed a given image to the generator space and incrementally inpaint the masked regions until a plausible image is generated; we trained a discriminator model to assess the quality of the generated images during the iterations and determine the convergence. After the generation process, for the post-processing, we utilize a Unet model that we trained specifically for this task to remedy the artifacts close to the mask boundaries. We empirically observed that only a few iterations are sufficient to generate realistic images with the proposed model. Since the models are not trained for particular mask types, our method allows applying sketch-based inpaintings, using a variety of mask types, and producing multiple and diverse results. We compared our method with the state-of-the-art models both quantitatively and qualitatively, and observed that our method can compete with the other models in all mask types; it is particularly better in images where larger masks are utilized. Our code, dataset and models are available at: https://github.com/yahyadogan72/iterative facial image inpainting.


abstract: SWELTO -- Space WEather Laboratory in Turin Observatory is a conceptual framework where new ideas for the analysis of space-based and ground-based data are developed and tested. The input data are (but not limited to) remote sensing observations (EUV images of the solar disk, Visible Light coronagraphic images, radio dynamic spectra, etc...), in situ plasma measurements (interplanetary plasma density, velocity, magnetic field, etc...), as well as measurements acquired by local sensors and detectors (radio antenna, fluxgate magnetometer, full-sky cameras, located in OATo). The output products are automatic identification, tracking, and monitoring of solar stationary and dynamic features near the Sun (coronal holes, active regions, coronal mass ejections, etc...), and in the interplanetary medium (shocks, plasmoids, corotating interaction regions, etc...), as well as reconstructions of the interplanetary medium where solar disturbances may propagate from the Sun to the Earth and beyond. These are based both on empirical models and numerical MHD simulations. The aim of SWELTO is not only to test new data analysis methods for future application for Space Weather monitoring and prediction purposes, but also to procure, test and deploy new ground-based instrumentation to monitor the ionospheric and geomagnetic responses to solar activity. Moreover, people involved in SWELTO are active in outreach to disseminate the topics related with Space Weather to students and the general public.


title: Computability of Data-Word Transductions over Different Data Domains

abstract: In this paper, we investigate the problem of synthesizing computable functions of infinite words over an infinite alphabet (data $\omega$-words). The notion of computability is defined through Turing machines with infinite inputs which can produce the corresponding infinite outputs in the limit. We use non-deterministic transducers equipped with registers, an extension of register automata with outputs, to describe specifications. Being non-deterministic, such transducers may not define functions but more generally relations of data $\omega$-words. In order to increase the expressive power of these machines, we even allow guessing of arbitrary data values when updating their registers. For functions over data $\omega$-words, we identify a sufficient condition (the possibility of determining the next letter to be outputted, which we call next letter problem) under which computability (resp. uniform computability) and continuity (resp. uniform continuity) coincide. We focus on two kinds of data domains: first, the general setting of oligomorphic data, which encompasses any data domain with equality, as well as the setting of rational numbers with linear order; and second, the set of natural numbers equipped with linear order. For both settings, we prove that functionality, i.e. determining whether the relation recognized by the transducer is actually a function, is decidable. We also show that the so-called next letter problem is decidable, yielding equivalence between (uniform) continuity and (uniform) computability. Last, we provide characterizations of (uniform) continuity, which allow us to prove that these notions, and thus also (uniform) computability, are decidable. We even show that all these decision problems are PSpace-complete for $(\mathbb{N},<)$ and for a large class of oligomorphic data domains, including for instance $(\mathbb{Q},<)$.


title: Optical phonon contribution to the thermal conductivity of a quantum paraelectric

abstract: Motivated by recent experimental findings, we study the contribution of a quantum critical optical phonon branch to the thermal conductivity of a paraelectric system. We consider the proximity of the optical phonon branch to transverse acoustic phonon branch and calculate its contribution to the thermal conductivity within the Kubo formalism. We find a low temperature power law dependence of the thermal conductivity as $T^{\alpha}$, with $1 < \alpha < 2$, (lower than $T^3$ behavior) due to optical phonons near the quantum critical point. This result is in accord with the experimental findings and indicates the importance of quantum fluctuations in the thermal conduction in these materials.


title: The potential for solar-diesel hybrid mini-grids in refugee camps: A case study of Nyabiheke camp, Rwanda

abstract: Electricity access in refugee camps is often limited to critical operations for humanitarian agencies and typically powered by fossil fuel generators. We study the economic and environmental benefits that optimised fully renewable and diesel-hybrid mini-grid designs can provide in humanitarian settings by displacing diesel use. Considering the case study of Nyabiheke camp in Rwanda we found that these benefits are substantial, with savings up to 32% of total costs and 83% of emissions, and cost payback times ranging from 0.9 to 6.2 years. Despite of their different cost structures, we find that all hybridisation levels of the system provide cost and emission savings compared to the incumbent diesel system. We highlight how modelling tools can facilitate the introduction and progressive expansion of systems as well as inform operational considerations on the ground. This study demonstrates how financial resources, environmental objectives and operational timeframes will influence the most appropriate system design for humanitarian actors on a case-by-case basis.


title: Shape programming lines of concentrated Gaussian curvature

abstract: Liquid crystal elastomers (LCEs) can undergo large reversible contractions along their nematic director upon heating or illumination. A spatially patterned director within a flat LCE sheet thus encodes a pattern of contraction on heating, which can morph the sheet into a curved shell, akin to how a pattern of growth sculpts a developing organism. Here we consider, theoretically, numerically and experimentally, patterns constructed from regions of radial and circular director, which, in isolation, would form cones and anticones. The resultant surfaces contain curved ridges with sharp V-shaped cross-sections, associated with the boundaries between regions in the patterns. Such ridges may be created in positively and negatively curved variants and, since they bear Gauss curvature (quantified here via the Gauss-Bonnet theorem), they cannot be flattened without energetically prohibitive stretch. Our experiments and numerics highlight that, although such ridges cannot be flattened isometrically, they can deform isometrically by trading the (singular) curvature of the V angle against the (finite) curvature of the ridge line. Furthermore, in finite thickness sheets, the sharp ridges are inevitably non-isometrically blunted to relieve bend, resulting in a modest smearing out of the encoded singular Gauss curvature. We close by discussing the use of such features as actuating linear features, such as probes, tongues and limbs, and highlighting the similarities between these patterns of shape change and those found during the morphogenesis of several biological systems.


title: CLASTER: Clustering with Reinforcement Learning for Zero-Shot Action Recognition

abstract: Zero-shot action recognition is the task of recognizingaction classes without visual examples, only with a seman-tic embedding which relates unseen to seen classes. Theproblem can be seen as learning a function which general-izes well to instances of unseen classes without losing dis-crimination between classes. Neural networks can modelthe complex boundaries between visual classes, which ex-plains their success as supervised models. However, inzero-shot learning, these highly specialized class bound-aries may not transfer well from seen to unseen classes.In this paper we propose a centroid-based representation,which clusters visual and semantic representation, consid-ers all training samples at once, and in this way generaliz-ing well to instances from unseen classes. We optimize theclustering using Reinforcement Learning which we show iscritical for our approach to work. We call the proposedmethod CLASTER and observe that it consistently outper-forms the state-of-the-art in all standard datasets, includ-ing UCF101, HMDB51 and Olympic Sports; both in thestandard zero-shot evaluation and the generalized zero-shotlearning. Further, we show that our model performs com-petitively in the image domain as well, outperforming thestate-of-the-art in many settings.


title: Online Caching with Optimal Switching Regret

abstract: We consider the classical uncoded caching problem from an online learning point-of-view. A cache of limited storage capacity can hold $C$ files at a time from a large catalog. A user requests an arbitrary file from the catalog at each time slot. Before the file request from the user arrives, a caching policy populates the cache with any $C$ files of its choice. In the case of a cache-hit, the policy receives a unit reward and zero rewards otherwise. In addition to that, there is a cost associated with fetching files to the cache, which we refer to as the switching cost. The objective is to design a caching policy that incurs minimal regret while considering both the rewards due to cache-hits and the switching cost due to the file fetches. The main contribution of this paper is the switching regret analysis of a Follow the Perturbed Leader-based anytime caching policy, which is shown to have an order optimal switching regret. In this pursuit, we improve the best-known switching regret bound for this problem by a factor of $\Theta(\sqrt{C}).$ We conclude the paper by comparing the performance of different popular caching policies using a publicly available trace from a commercial CDN server.


title: Magnetoelastic properties of a spin-1/2 Ising-Heisenberg diamond chain in vicinity of a triple coexistence point

abstract: We study magnetoelastic properties of a spin-1/2 Ising-Heisenberg diamond chain, whose elementary unit cell consists of two decorating Heisenberg spins and one nodal Ising spin. It is assumed that each couple of the decorating atoms including the Heisenberg spins harmonically vibrates perpendicularly to the chain axis, while the nodal atoms involving the Ising spins are placed at rigid positions when ignoring their lattice vibrations. An effect of the magnetoelastic coupling on a ground state and finite-temperature properties is particularly investigated close to a triple coexistence point depending on a spring-stiffness constant ascribed to the Heisenberg interaction. The magnetoelastic nature of the Heisenberg dimers is reflected through a non-null plateau of the entropy emergent in a low-temperature region, whereas the specific heat displays an anomalous peak slightly below the temperature region corresponding to the entropy plateau. The magnetization also exhibits a plateau in the same temperature region at almost saturated value before it gradually tends to zero upon increasing of temperature. The magnetic susceptibility displays within the plateau region an inverse temperature dependence, which slightly drops above this plateau, whereas an inverse temperature dependence is repeatedly recovered at high enough temperatures.


title: The phase diagram of the Hubbard model by Variational Auxiliary Field quantum Monte Carlo

abstract: A systematically improvable wave function is proposed for the numerical solution of strongly correlated systems. With a stochastic optimization method, based on the auxiliary field quantum Monte Carlo technique, an effective temperature Teff is defined, probing the distance of the ground state properties of the model in the thermodynamic limit from the ones of the proposed correlated mean-field ansatz. In this way their uncertainties from the unbiased zero temperature limit may be estimated by simple and stable extrapolations well before the so called sign problem gets prohibitive. At finite Teff the convergence of the energy to the thermodynamic limit is indeed shown to be possible in the Hubbard model already for relatively small square lattices with linear dimension L ~10, thanks to appropriate averages over several twisted boundary conditions. Within the estimated energy accuracy of the proposed variational ansatz, two clear phases are identified, as the energy is lowered by spontaneously breaking some symmetries satisfied by the Hubbard Hamiltonian: a) a stripe phase where both spin and translation symmetries are broken, and b) a strong coupling d-wave superconducting phase when the particle number is not conserved and global U(1) symmetry is broken. On the other hand the symmetric phase is stable in a wide region at large doping and small coupling.


title: Mind the Gap when Conditioning Amortised Inference in Sequential Latent-Variable Models

abstract: Amortised inference enables scalable learning of sequential latent-variable models (LVMs) with the evidence lower bound (ELBO). In this setting, variational posteriors are often only partially conditioned. While the true posteriors depend, e.g., on the entire sequence of observations, approximate posteriors are only informed by past observations. This mimics the Bayesian filter -- a mixture of smoothing posteriors. Yet, we show that the ELBO objective forces partially-conditioned amortised posteriors to approximate products of smoothing posteriors instead. Consequently, the learned generative model is compromised. We demonstrate these theoretical findings in three scenarios: traffic flow, handwritten digits, and aerial vehicle dynamics. Using fully-conditioned approximate posteriors, performance improves in terms of generative modelling and multi-step prediction.


title: In situ AC-hysteresis measurements of SPD-processed Cu20(Fe15Co85)80

abstract: The changes of magnetic properties upon heat treatment of a metastable supersaturated solid solution processed by severe plastic deformation are investigated by in-situ AC-hysteresis measurements. Data are analyzed in the framework of dynamic loss theory, with correlative investigations of the microstructural properties. The evolution of hysteresis upon annealing points out that the single-phase supersaturated solid solution remains stable up to 400{\deg}C, then hindering of domain wall motion sets in at this temperature. At 600{\deg}C, a multi phase microstructure is present, causing a significant increase in coercivity.


title: Deadeye: A Novel Preattentive Visualization Technique Based on Dichoptic Presentation

abstract: Preattentive visual features such as hue or flickering can effectively draw attention to an object of interest -- for instance, an important feature in a scientific visualization. These features appear to pop out and can be recognized by our visual system, independently from the number of distractors. Most cues do not take advantage of the fact that most humans have two eyes. In cases where binocular vision is applied, it is almost exclusively used to convey depth by exposing stereo pairs. We present Deadeye, a novel preattentive visualization technique based on presenting different stimuli to each eye. The target object is rendered for one eye only and is instantly detected by our visual system. In contrast to existing cues, Deadeye does not modify any visual properties of the target and, thus, is particularly suited for visualization applications. Our evaluation confirms that Deadeye is indeed perceived preattentively. We also explore a conjunction search based on our technique and show that, in contrast to 3D depth, the task cannot be processed in parallel.


title: Integrating theory and experiments to link local mechanisms and ecosystem-level consequences of vegetation patterns in drylands

abstract: Self-organized spatial patterns of vegetation are frequent in drylands and, because pattern shape correlates with water availability, they have been suggested as important indicators of ecosystem health. However, the mechanisms underlying pattern emergence remain unclear. Some theories hypothesize that patterns could result from a water-mediated scale-dependent feedback (SDF) whereby interactions favoring plant growth dominate at short distances and growth-inhibitory interactions dominate in the long range. However, we know little about how the presence of a focal plant affects the fitness of its neighbors as a function of the inter-individual distance, which is expected to be highly ecosystem-dependent. This lack of empirical knowledge and system dependency challenge the relevance of SDF as a unifying theory for vegetation pattern formation. Assuming that plant interactions are always inhibitory and only their intensity is scale-dependent, alternative theories also recover the typical vegetation patterns observed in nature. Importantly, although these alternative hypotheses lead to visually indistinguishable patterns, they predict contrasting desertification dynamics, which questions the potential use of vegetation patterns as ecosystem-state indicators. To help resolve this issue, we first review existing theories for vegetation self-organization and their conflicting predictions about desertification dynamics. Second, we discuss potential empirical tests via manipulative experiments to identify pattern-forming mechanisms and link them to specific desertification dynamics. A comprehensive view of models, the mechanisms they intend to capture, and experiments to test them in the field will help to better understand both how patterns emerge and improve predictions on the fate of the ecosystems where they form.


title: A Distributed Chunk Calculation Approach for Self-scheduling of Parallel Applications on Distributed-memory Systems

abstract: Loop scheduling techniques aim to achieve load-balanced executions of scientific applications. Dynamic loop self-scheduling (DLS) libraries for distributed-memory systems are typically MPI-based and employ a centralized chunk calculation approach (CCA) to assign variably-sized chunks of loop iterations. We present a distributed chunk calculation approach (DCA) that supports various types of DLS techniques. Using both CCA and DCA, twelve DLS techniques are implemented and evaluated in different CPU slowdown scenarios. The results show that the DLS techniques implemented using DCA outperform their corresponding ones implemented with CCA, especially in extreme system slowdown scenarios.


title: Electron-hole asymmetry in electron systems with orbital degeneracy and correlated hopping

abstract: Microscopic models of electronic subsystems with orbital degeneracy of energy states and non-diagonal matrix elements of electron interactions (correlated hopping) are considered within the configuration-operator approach. Equations for arbitrary-temperature numerical calculation of the doublon concentration for the integer band filling $n=1$ at different forms of the model density of states are derived. The energy spectra obtained within the Green function method are discussed with special emphasis on the role of correlated hopping in transition from itinerant to localized behavior observed in vanadium Magneli phases V$_n$O$_{2n-1}$, transition-metal dichalcogenides NiS$_{2-x}$Se$_x$, fulleride A$_n$C$_{60}$ systems.


title: Random phase approximation with exchange for an accurate description of crystalline polymorphism

abstract: We determine the correlation energy of BN, SiO$_2$ and ice polymorphs employing a recently developed RPAx (random phase approximation with exchange) approach. The RPAx provides larger and more accurate polarizabilities as compared to the RPA, and captures effects of anisotropy. In turn, the correlation energy, defined as an integral over the density-density response function, gives improved binding energies without the need for error cancellation. Here, we demonstrate that these features are crucial for predicting the relative energies between low- and high-pressure polymorphs of different coordination number as, e.g., between $\alpha$-quartz and stishovite in SiO$_2$, or layered and cubic BN. Furthermore, a reliable (H$_2$O)$_2$ potential energy surface is obtained, necessary for describing the various phases of ice. The RPAx gives results comparable to other high-level methods such as coupled cluster and quantum Monte Carlo, also in cases where the RPA breaks down. Although higher computational cost than RPA we observe a faster convergence with respect to the number of eigenvalues in the response function.


title: A Passive Online Technique for Learning Hybrid Automata from Input/Output Traces

abstract: Specification synthesis is the process of deriving a model from the input-output traces of a system. It is used extensively in test design, reverse engineering, and system identification. One type of the resulting artifact of this process for cyber-physical systems is hybrid automata. They are intuitive, precise, tool independent, and at a high level of abstraction, and can model systems with both discrete and continuous variables. In this paper, we propose a new technique for synthesizing hybrid automaton from the input-output traces of a non-linear cyber-physical system. Similarity detection in non-linear behaviors is the main challenge for extracting such models. We address this problem by utilizing the Dynamic Time Warping technique. Our approach is passive, meaning that it does not need interaction with the system during automata synthesis from the logged traces; and online, which means that each input/output trace is used only once in the procedure. In other words, each new trace can be used to improve the already synthesized automaton. We evaluated our algorithm in two industrial and simulated case studies. The accuracy of the derived automata show promising results.


title: Spinning-off stringy electro-magnetic memories

abstract: We extend and generalise the string corrections to the EM memory to the Type I superstring including spin effects. Very much as in the simpler bosonic string context, the relevant corrections are non-perturbative in $\alpha'$, slowly decaying (as $1/R$) at large distances and modulated in retarded time $u=t-R$. For spin $N$ states in the first Regge trajectory they entail a sequence of $N$ derivatives wrt $u$ on the `parent' $N=0$ amplitude. We also briefly discuss how to include loop effects, that broaden and shift the string resonances, and how to modify our analysis for macroscopic semi-classical quasi-BPS coherent states, whose collisions may lead to detectable string memory signals in viable Type I models.


title: Explicit continuation methods with L-BFGS updating formulas for linearly constrained optimization problems

abstract: This paper considers an explicit continuation method with the trusty time-stepping scheme and the limited-memory BFGS (L-BFGS) updating formula (Eptctr) for the linearly constrained optimization problem. At every iteration, Eptctr only involves three pairs of the inner product of vector and one matrix-vector product, other than the traditional and representative optimization method such as the sequential quadratic programming (SQP) or the latest continuation method such as Ptctr \cite{LLS2020}, which needs to solve a quadratic programming subproblem (SQP) or a linear system of equations (Ptctr). Thus, Eptctr can save much more computational time than SQP or Ptctr. Numerical results also show that the consumed time of EPtctr is about one tenth of that of Ptctr or one fifteenth to 0.4 percent of that of SQP. Furthermore, Eptctr can save the storage space of an $(n+m) \times (n+m)$ large-scale matrix, in comparison to SQP. The required memory of Eptctr is about one fifth of that of SQP. Finally, we also give the global convergence analysis of the new method under the standard assumptions.


title: Relative Log-Symplectic structure on a semi-stable degeneration of moduli of Higgs bundles

abstract: In a recent paper \cite{3}, a semi-stable degeneration of moduli space of Higgs bundles on a curve has been constructed. In this paper, we show that there is a relative log-symplectic form on this degeneration, whose restriction to the generic fibre is the classical symplectic form discovered by Hitchin. We compute the Poisson ranks at every point and describe the symplectic foliation on the closed fibre. We also show that the closed fibre, which is a variety with normal crossing singularities, acquires a structure of an algebraically completely integrable system.


title: Stabilized finite element method for incompressible solid dynamics using an updated Lagrangian formulation

abstract: This paper proposes a novel way to solve transient linear, and non-linear solid dynamics for compressible, nearly incompressible, and incompressible material in the updated Lagrangian framework for tetrahedral unstructured finite elements. It consists of a mixed formulation in both displacement and pressure, where the momentum equation of the continuum is complemented with a pressure equation that handles incompresibility inherently. It is obtained through the deviatoric and volumetric split of the stress, that enables us to solve the problem in the incompressible limit. The Varitaional Multi-Scale method (VMS) is developed based on the orthogonal decomposition of the variables, which damps out spurious pressure fields for piece wise linear tetrahedral elements. Various numerical examples are presented to assess the robustness, accuracy and capabilities of our scheme in bending dominated problems, and for complex geometries.


title: Active instability and nonlinear dynamics of cell-cell junctions

abstract: Active cell-junction remodeling is important for tissue morphogenesis, yet its underlying physics is not understood. We study a mechanical model that describes junctions as dynamic active force dipoles. Their instability can trigger cell intercalations by a critical collapse. Nonlinearities in tissue's elastic response can stabilize the collapse either by a limit cycle or condensation of junction lengths at cusps of the energy landscape. Furthermore, active junction networks undergo collective instability to drive active in-plane ordering or develop a limit cycle of collective oscillations, which extends over regions of the energy landscape corresponding to distinct network topologies.


title: Simple closed geodesics on regular tetrahedra in spherical space

abstract: On a regular tetrahedron in spherical space there exist the finite number of simple closed geodesics. For any pair of coprime integers $(p,q)$ it was found the numbers $\alpha_1$ and $\alpha_2$ depending on $p$, $q$ and satisfying the inequalities $\pi/3< \alpha_1 < \alpha_2 < 2\pi/3$ such that on a regular tetrahedron in spherical space with the faces angle $\alpha \in \left( \pi/3, \alpha_1 \right)$ there exists unique, up to the rigid motion of the tetrahedron, simple closed geodesic of type $(p,q)$, and on a regular tetrahedron with the faces angle $\alpha \in \left( \alpha_2, 2\pi/3 \right)$ there is no simple closed geodesic of type $(p,q)$.


title: Linearization stability of reflection-asymmetric thin-shell wormholes with double shadows

abstract: Wormholes are hypothetical objects which can be black hole mimickers with strong gravitational fields. Recently, Wielgus et al. have constructed reflection-asymmetric thin-shell wormholes which are composed of the parts of a Schwarzschild spacetime and a Reissner-Nordstrom spacetime with two photon spheres with different sizes each other in [M. Wielgus, J. Horak, F. Vincent, and M. Abramowicz, Phys. Rev. D 102, 084044 (2020)]. They have discussed observational property of the shadows with two photon rings in different sizes as seen from an observer and they have named their shadows double shadows. In this paper, we study the linearization stability of the reflection-asymmetric thin-shell wormholes with the double shadows.


title: Deep Inertial Odometry with Accurate IMU Preintegration

abstract: Inertial Measurement Units (IMUs) are interceptive modalities that provide ego-motion measurements independent of the environmental factors. They are widely adopted in various autonomous systems. Motivated by the limitations in processing the noisy measurements from these sensors using their mathematical models, researchers have recently proposed various deep learning architectures to estimate inertial odometry in an end-to-end manner. Nevertheless, the high-frequency and redundant measurements from IMUs lead to long raw sequences to be processed. In this study, we aim to investigate the efficacy of accurate preintegration as a more realistic solution to the IMU motion model for deep inertial odometry (DIO) and the resultant DIO is a fusion of model-driven and data-driven approaches. The accurate IMU preintegration has the potential to outperform numerical approximation of the continuous IMU model used in the existing DIOs. Experimental results validate the proposed DIO.


title: On the Davey-Stewartson hierarchy: construction by two scalar pseudo-differential operators and compatibility for infinite many flows

abstract: The infinite many symmetries of Davey-Stewartson (DS) system are closely connected to the integrable deformations of surfaces in a four-dimensional space. In this paper, we give a direct algorithm to construct the expression of the DS hierarchy by two scalar pseudo-differential operators involving partial derivatives.


title: Persistent current noise in narrow Josephson junctions

abstract: Josephson junctions have broad applications in metrology, quantum information processing, and remote sensing. For these applications, the electronic noise is a limiting factor. In this work we study the thermal noise in narrow Josephson junctions using a tight-binding Hamiltonian. For a junction longer than the superconducting coherence length, several self-consistent gap profiles appear close to a phase difference $\pi$. They correspond to two stable solutions with an approximately constant phase gradient over the thin superconductor connected by a $2\pi$ phase slip, and a solitonic branch. The current noise power spectrum has pronounced peaks at the transition frequencies between the different states in each branch. We find that the noise is reduced in the gradient branches in comparison to the zero-length junction limit. In contrast, the solitonic branch exhibits an enhanced noise and a reduced current due to the pinning of the lowest excitation energy to close to zero energy.


title: Cartan approach to Teleparallel Equivalent to General Relativity: a review

abstract: In previous works, questioning the mathematical nature of the connection in the translations gauge theory formulation of Teleparallel Equivalent to General Relativity (TEGR) Theory led us to propose a new formulation using a Cartan connection. In this review, we summarize the presentation of that proposal and discuss it from a gauge theoretic perspective.


title: Rigorous Bounds on the Heating Rate in Thue-Morse Quasiperiodically and Randomly Driven Quantum Many-Body Systems

abstract: The nonequilibrium quantum dynamics of closed many-body systems is a rich yet challenging field. While recent progress for periodically driven (Floquet) systems has yielded a number of rigorous results, our understanding on quantum many-body systems driven by rapidly varying but a- and quasi-periodic driving is still limited. Here, we derive rigorous, non-perturbative, bounds on the heating rate in quantum many-body systems under Thue-Morse quasi-periodic driving and under random multipolar driving, the latter being a tunably randomized variant of the former. In the process, we derive a static effective Hamiltonian that describes the transient prethermal state, including the dynamics of local observables. Our bound for Thue-Morse quasi-periodic driving suggests that the heating time scales like $(\omega/g)^{-C\ln(\omega/g)}$ with a positive constant $C$ and a typical energy scale $g$ of the Hamiltonian, in agreement with our numerical simulations.


title: Reversible Computation in Petri Nets

abstract: Reversible computation is an unconventional form of computing that extends the standard forward-only mode of computation with the ability to execute a sequence of operations in reverse at any point during computation. As such, in this thesis we propose a reversible approach to Petri nets by introducing machinery and associated operational semantics to tackle the challenges of the main forms of reversibility. Our proposal concerns a variation of cyclic Petri nets, called Reversing Petri Nets (RPNs) where tokens are persistent and distinguished from each other by an identity. An immediate extension of the original model includes allowing multiple tokens of the same base/type to occur in a model. Specifically, we explore the individual token interpretation where one distinguishes different tokens residing in the same place by keeping track of where they come from. We also propose the collective token interpretation, as the opposite approach to token ambiguity, which considers all tokens of a certain type to be identical, disregarding their history during execution. Both of the proposed models of RPNs (with single or multi tokens) implement the notion of uncontrolled reversibility, meaning that it specifies how to reverse an execution and allows to do so freely, yet it places no restrictions as to when and whether to prefer backward execution over forward execution or vice versa. In this respect, a further aim is to control reversibility by extending our formal semantics where transitions are associated with conditions whose satisfaction allows the execution of transitions in the forward/reversed direction.


title: Data Obsolescence Detection in the Light of Newly Acquired Valid Observations

abstract: The information describing the conditions of a system or a person is constantly evolving and may become obsolete and contradict other information. A database, therefore, must be consistently updated upon the acquisition of new valid observations that contradict obsolete ones contained in the database. In this paper, we propose a novel approach for dealing with the information obsolescence problem. Our approach aims to detect, in real-time, contradictions between observations and then identify the obsolete ones, given a representation model. Since we work within an uncertain environment characterized by the lack of information, we choose to use a Bayesian network as our representation model and propose a new approximate concept, $\epsilon$-Contradiction. The new concept is parameterised by a confidence level of having a contradiction in a set of observations. We propose a polynomial-time algorithm for detecting obsolete information. We show that the resulting obsolete information is better represented by an AND-OR tree than a simple set of observations. Finally, we demonstrate the effectiveness of our approach on a real elderly fall-prevention database and showcase how this tree can be used to give reliable recommendations to doctors. Our experiments give systematically and substantially very good results.


title: High Order Asymptotic Expansions of a Good-Bad-Ugly Wave Equation

abstract: A heuristic method to find asymptotic solutions to a system of non-linear wave equations near null infinity is proposed. The non-linearities in this model, dubbed good-bad-ugly, are known to mimic the ones present in the Einstein field equations (EFE) and we expect to be able to exploit this method to derive an asymptotic expansion for the metric in General Relativity (GR) close to null infinity that goes beyond first order as performed by Lindblad and Rodnianski for the leading asymptotics. For the good-bad-ugly model, we derive formal expansions in which terms proportional to the logarithm of the radial coordinate appear at every order in the bad field, from the second order onward in the ugly field but never in the good field. The model is generalized to wave operators built from an asymptotically flat metric and it is shown that it admits polyhomogeneous asymptotic solutions. Finally we define stratified null forms, a generalization of standard null forms, which capture the behavior of different types of field, and demonstrate that the addition of such terms to the original system bears no qualitative influence on the type of asymptotic solutions found.


title: Emotional EEG Classification using Connectivity Features and Convolutional Neural Networks

abstract: Convolutional neural networks (CNNs) are widely used to recognize the user's state through electroencephalography (EEG) signals. In the previous studies, the EEG signals are usually fed into the CNNs in the form of high-dimensional raw data. However, this approach makes it difficult to exploit the brain connectivity information that can be effective in describing the functional brain network and estimating the perceptual state of the user. We introduce a new classification system that utilizes brain connectivity with a CNN and validate its effectiveness via the emotional video classification by using three different types of connectivity measures. Furthermore, two data-driven methods to construct the connectivity matrix are proposed to maximize classification performance. Further analysis reveals that the level of concentration of the brain connectivity related to the emotional property of the target video is correlated with classification performance.


title: Computational design of armored nanodroplets as nanocarriers for encapsulation and release under flow conditions

abstract: Nanocarriers are nanosized materials commonly used for targeted-oriented delivery of active compounds, including antimicrobials and small-molecular drugs. They equally represent fundamental and engineering challenges since sophisticated nanocarriers must show adequate structure, stability, and function in complex ambients. Here, we report on the computational design of a distinctive class of nanocarriers, built from buckled armored nanodroplets, able to selectively encapsulate or release a probe load under specific flow conditions. Mesoscopic simulations offer detailed insight into the interplay between the characteristics of laden surface coverage and evolution of the droplet morphology. First, we describe in detail the formation of \textit{pocket-like} structures in Pickering emulsion nanodroplets and their stability under external flow. Then we use that knowledge to test the capacity of these emulsion-based pockets to yield flow-assisted encapsulation or expulsion of a probe load. Finally, the rheological properties of our model carrier are put into perspective with those of delivery systems employed in pharmaceutical and cosmetic technology.


title: Thermal relaxation in metal films limited by diffuson lattice excitations of amorphous substrates

abstract: Here we examine the role of the amorphous insulating substrate in the thermal relaxation in thin NbN, InO$_x$, and Au/Ni films at temperatures above 5 K. The studied samples are made up of metal bridges on an amorphous insulating layer lying on or suspended above a crystalline substrate. Noise thermometry was used to measure the electron temperature $T_e$ of the films as a function of Joule power per unit of area $P_{2D}$. In all samples, we observe the dependence $P_{2D}\propto T_e^n$ with the exponent $n\simeq 2$, which is inconsistent with both electron-phonon coupling and Kapitza thermal resistance. In suspended samples, the functional dependence of $P_{2D}(T_e)$ on the length of the amorphous insulating layer is consistent with the linear $T$-dependence of the thermal conductivity, which is related to lattice excitations (diffusons) for the phonon mean free path smaller than the dominant phonon wavelength. Our findings are important for understanding the operation of devices embedded in amorphous dielectrics.


title: Old and New Major Mergers in the SOSIMPLE galaxy, NGC 7135

abstract: The simultaneous advancement of high resolution integral field unit spectroscopy and robust full-spectral fitting codes now make it possible to examine spatially-resolved kinematic, chemical composition, and star-formation history from nearby galaxies. We take new MUSE data from the Snapshot Optical Spectroscopic Imaging of Mergers and Pairs for Legacy Exploration (SOSIMPLE) survey to examine NGC 7135. With counter-rotation of gas, disrupted kinematics and asymmetric chemical distribution, NGC 7135 is consistent with an ongoing merger. Though well hidden by the current merger, we are able to distinguish stars originating from an older merger, occurring 6-10 Gyr ago. We further find a gradient in ex-situ material with galactocentric radius, with the accreted fraction rising from 0% in the galaxy centre, to ~7% within 0.6 effective radii.


title: Sum-Rate Maximization in Distributed Intelligent Reflecting Surfaces-Aided mmWave Communications

abstract: In this paper, we focus on the sum-rate optimization in a multi-user millimeter-wave (mmWave) system with distributed intelligent reflecting surfaces (D-IRSs), where a base station (BS) communicates with users via multiple IRSs. The BS transmit beamforming, IRS switch vector, and phase shifts of the IRS are jointly optimized to maximize the sum-rate under minimum user rate, unit-modulus, and transmit power constraints. To solve the resulting non-convex optimization problem, we develop an efficient alternating optimization (AO) algorithm. Specifically, the non-convex problem is converted into three subproblems, which are solved alternatively. The solution to transmit beamforming at the BS and the phase shifts at the IRS are derived by using the successive convex approximation (SCA)-based algorithm, and a greedy algorithm is proposed to design the IRS switch vector. The complexity of the proposed AO algorithm is analyzed theoretically. Numerical results show that the D-IRSs-aided scheme can significantly improve the sum-rate and energy efficiency performance.


title: Set Partitions and Other Bell Number Enumerated Objects

abstract: In this paper, we study classes of subexcedant functions enumerated by the Bell numbers and present bijections on set partitions. We present a set of permutations whose transposition arrays are the restricted growth functions, thus defining Bell permutations of the second kind. We describe a bijection between Bell permutations of the first kind (introduced by Ponti and Vajnovzski) and the second kind. We present two other Bell number enumerated classes of subexcedant functions. Further, we present bijections on set partitions, in particular, an involution that interchanges the set of merging blocks and the set of successions. We use the bijections to enumerate the distribution of these statistics over the set of set partitions, and also give some enumeration results.


title: Fundamental energy cost of finite-time computing

abstract: The fundamental energy cost of irreversible computing is given by the Landauer bound of $kT \ln2$~/bit. However, this limit is only achievable for infinite-time processes. We here determine the fundamental energy cost of finite-time irreversible computing \er{within the framework of nonequilibrium thermodynamics}. Comparing the lower bounds of energy required by ideal serial and parallel computers to solve a problem of a given size in a given finite time, we find that the energy cost of a serial computer fundamentally diverges with increasing problem size, whereas that of a parallel computer can stay close to the Landauer limit. We discuss the implications of this result in the context of current technology, and for different degrees of parallelization and amounts of overhead. Our findings provide a physical basis for the design of energy efficient computers.


title: A note on broken dilatation symmetry in planar noncommutative theory

abstract: A study of a riveting connection between noncommutativity and the anomalous dilatation (scale) symmetry is presented for a generalized quantum Hall system due to time dilatation transformations. On using the "Peierls substitution" scheme, it is shown that noncommutativity between spatial coordinates emerges naturally at a large magnetic field limit. Thereafter, we derive a path-integral action for the corresponding noncommutative quantum system and discuss the equivalence between the considered noncommutative system and the generalized Landau problem thus rendering an effective commmutative description. By exploiting the path-integral method due to Fujikawa, we derive an expression for the unintegrated scale or dilatation anomaly for the generalized Landau system, wherein the anomalies are identified with Jacobian factors arising from measure change under scale transformation and is subsequently renormalised. In fact, we derive exact expressions of anomalous Ward identities from which one may point out the existence of scale anomaly which is a purely quantum effect induced from the noncommutative structure between spatial coordinates.


title: Yet Another Representation of Binary Decision Trees: A Mathematical Demonstration

abstract: A decision tree looks like a simple computational graph without cycles, where only the leaf nodes specify the output values and the non-terminals specify their tests or split conditions. From the numerical perspective, we express decision trees in the language of computational graph. We explicitly parameterize the test phase, traversal phase and prediction phase of decision trees based on the bitvectors of non-terminal nodes. As shown later, the decision tree is a shallow binary network in some sense. Especially, we introduce the bitvector matrix to implement the tree traversal in numerical approach, where the core is to convert the logical `AND' operation to arithmetic operations. And we apply this numerical representation to extend and unify diverse decision trees in concept.


title: SoK: Fully Homomorphic Encryption Compilers

abstract: Fully Homomorphic Encryption (FHE) allows a third party to perform arbitrary computations on encrypted data, learning neither the inputs nor the computation results. Hence, it provides resilience in situations where computations are carried out by an untrusted or potentially compromised party. This powerful concept was first conceived by Rivest et al. in the 1970s. However, it remained unrealized until Craig Gentry presented the first feasible FHE scheme in 2009. The advent of the massive collection of sensitive data in cloud services, coupled with a plague of data breaches, moved highly regulated businesses to increasingly demand confidential and secure computing solutions. This demand, in turn, has led to a recent surge in the development of FHE tools. To understand the landscape of recent FHE tool developments, we conduct an extensive survey and experimental evaluation to explore the current state of the art and identify areas for future development. In this paper, we survey, evaluate, and systematize FHE tools and compilers. We perform experiments to evaluate these tools' performance and usability aspects on a variety of applications. We conclude with recommendations for developers intending to develop FHE-based applications and a discussion on future directions for FHE tools development.


title: Some Examples of Family Floer Mirrors

abstract: In this article, we give explicit calculations for the family Floer mirrors of some non-compact Calabi-Yau surfaces. We compare it with the mirror construction of Gross-Hacking-Keel for suitably chosen log Calabi-Yau pairs and the rank two cluster varieties of finite type. In particular, the analytifications of the later two give partial compactifications of the family Floer mirrors that we computed.


title: The Stellar Distribution Function and Local Vertical Potential from Gaia DR2

abstract: We develop a novel method to simultaneously determine the vertical potential, force and stellar $z-v_z$ phase space distribution function (DF) in our local patch of the Galaxy. We assume that the Solar Neighborhood can be treated as a one-dimensional system in dynamical equilibrium and directly fit the number density in the $z-v_z$ plane to what we call the Rational Linear DF (RLDF) model. This model can be regarded as a continuous sum of isothermal DFs though it has only one more parameter than the isothermal model. We apply our method to a sample of giant stars from Gaia Data Release 2 and show that the RLDF provides an excellent fit to the data. The well-known phase space spiral emerges in the residual map of the $z-v_z$ plane. We use the best-fit potential to plot the residuals in terms of the frequency and angle of vertical oscillations and show that the spiral maps into a straight line. From its slope, we estimate that the phase spirals were generated by a perturbation $\sim540$ Myr years ago. We also determine the differential surface density as a function of vertical velocity dispersion, a.k.a. the vertical temperature distribution. The result is qualitatively similar to what was previously found for SDSS/SEGUE G dwarfs. Finally, we address parameter degeneracies and the validity of the 1D approximation. Particularly, the mid-plane density derived from a cold subsample, where the 1D approximation is more secure, is closer to literature values than that derived from the sample as a whole.


title: Merging-Free Partitions and Run-Sorted Permutations

abstract: In this paper, we study merging-free partitions with their canonical forms and run-sorted permutations. We give a combinatorial proof of the conjecture made by Nabawanda et al. We describe the distribution of the statistics of runs and right-to-left minima over the set of run-sorted permutations and we give the exponential generating function for their joint distribution. We show the number of right-to-left minima is given by the shifted distribution of the Stirling number of the second kind. We also prove that the non-crossing merging-free partitions are enumerated by powers of 2. We use one of the constructive proofs given in the paper to implement an algorithm for the exhaustive generation of run-sorted permutations by number of runs.


title: Some paradoxes in the Tilman model, how to avoid or accept them

abstract: In making models for biological systems one expects to grasp the biology and be able to use ones intuition to predict the outcome. This paper is about the discrepancy between what is expected and what is the outcome of the analysis of the model, a paradox. The Tilman model of consumers competing for resources has some aspects that may appear counterintuitive. Within the standard Tilman model for a single consumer and a single resource, when a very efficient consumer rapidly eats all available food, the resource density becomes zero, but if there is no food, how can the consumer survive? The paradox can be resolved by realising that in the short term indeed rapid consumption may lead to starvation and a decline in the consumer population, but in the long term a small resource and consumer density remain. A single consumer living on two essential nutrients leaves the density of the non-limiting nutrient above its critical level, so what is done with the extra food? We explain that the extra food is not consumed, just because it is available, but that in fact the additional food is not consumed at all, because the consumer has no use for it. For a model with two consumer species competing for a single nutrient one of the consumers will eventually disappear, even if the food economics of the surviving one is worse and there is still much food wasted. Both are inherent aspects of the model, and the paradox can be avoided if the difference between the consumer species is small, in which case it will take very long to reach equilibrium. We argue that in general an extended stability analysis, in which not only the asymptotically stable state is considered, but also the unstable steady states and all time scales involved in the transient dynamics, gives more insight in the dynamics involved and thus can help in avoiding apparent contradictions in ecological models, or accepting them.


title: Rigidity of minimal Lagrangian diffeomorphisms between spherical cone surfaces

abstract: We prove that any minimal Lagrangian diffeomorphism between two closed spherical surfaces with cone singularities is an isometry, without any assumption on the multiangles of the two surfaces. As an application, we show that every branched immersion of a closed surface of constant positive Gaussian curvature in Euclidean three-space is a branched covering onto a round sphere, thus generalizing the classical rigidity theorem of Liebmann to branched immersions.


title: Beating the Market with Generalized Generating Portfolios

abstract: Stochastic portfolio theory aims at finding relative arbitrages, i.e. trading strategies which outperform the market with probability one. Functionally generated portfolios, which are deterministic functions of the market weights, are an invaluable tool in doing so. Driven by a practitioner point of view, where investment decisions are based upon consideration of various financial variables, we generalize functionally generated portfolios and allow them to depend on continuous-path semimartingales, in addition to the market weights. By means of examples we demonstrate how the inclusion of additional processes can reduce time horizons beyond which relative arbitrage is possible, boost performance of generated portfolios, and how investor preferences and specific investment views can be included in the context of stochastic portfolio theory. Striking is also the construction of a relative arbitrage opportunity which is generated by the volatility of the additional semimartingale. An in-depth empirical analysis of the performance of the proposed strategies confirms our theoretical findings and demonstrates that our portfolios represent profitable investment opportunities even in the presence of transaction costs.


title: Training collective variables for enhanced sampling via neural networks based discriminant analysis

abstract: A popular way to accelerate the sampling of rare events in molecular dynamics simulations is to introduce a potential that increases the fluctuations of selected collective variables. For this strategy to be successful, it is critical to choose appropriate variables. Here we review some recent developments in the data-driven design of collective variables, with a focus on the combination of Fisher's discriminant analysis and neural networks. This approach allows to compress the fluctuations of metastable states into a low-dimensional representation. We illustrate through several examples the effectiveness of this method in accelerating the sampling, while also identifying the physical descriptors that undergo the most significant changes in the process.


title: Model Compression for Domain Adaptation through Causal Effect Estimation

abstract: Recent improvements in the predictive quality of natural language processing systems are often dependent on a substantial increase in the number of model parameters. This has led to various attempts of compressing such models, but existing methods have not considered the differences in the predictive power of various model components or in the generalizability of the compressed models. To understand the connection between model compression and out-of-distribution generalization, we define the task of compressing language representation models such that they perform best in a domain adaptation setting. We choose to address this problem from a causal perspective, attempting to estimate the average treatment effect (ATE) of a model component, such as a single layer, on the model's predictions. Our proposed ATE-guided Model Compression scheme (AMoC), generates many model candidates, differing by the model components that were removed. Then, we select the best candidate through a stepwise regression model that utilizes the ATE to predict the expected performance on the target domain. AMoC outperforms strong baselines on dozens of domain pairs across three text classification and sequence tagging tasks.


title: Higher-order error estimates of the discrete-time Clark--Ocone formula

abstract: In this article, we investigate the convergence rate of the discrete-time Clark--Ocone formula provided by Akahori--Amaba--Okuma [1]. In that paper, they mainly focus on the $L_{2}$-convergence rate of the first-order error estimate related to the tracking error of the delta hedge in mathematical finance. Here, as two extensions, we estimate "the higher order error" for Wiener functionals with an integrability index $2$ and "an arbitrary differentiability index."


title: A fast spectral method for electrostatics in doubly-periodic slit channels

abstract: We develop a fast method for computing the electrostatic energy and forces for a collection of charges in doubly-periodic slabs with jumps in the dielectric permittivity at the slab boundaries. Our method achieves spectral accuracy by using Ewald splitting to replace the original Poisson equation for nearly-singular sources with a smooth far-field Poisson equation, combined with a localized near-field correction. Unlike existing spectral Ewald methods, which make use of the Fourier transform in the aperiodic direction, we recast the problem as a two-point boundary value problem in the aperiodic direction for each transverse Fourier mode, for which exact analytic boundary conditions are available. We solve each of these boundary value problems using a fast, well-conditioned Chebyshev method. In the presence of dielectric jumps, combining Ewald splitting with the classical method of images results in smoothed charge distributions which overlap the dielectric boundaries themselves. We show how to preserve spectral accuracy in this case through the use of a harmonic correction which involves solving a simple Laplace equation with smooth boundary data. We implement our method on Graphical Processing Units, and combine our doubly-periodic Poisson solver with Brownian Dynamics to study the equilibrium structure of double layers in binary electrolytes confined by dielectric boundaries. Consistent with prior studies, we find strong charge depletion near the interfaces due to repulsive interactions with image charges, which points to the need for incorporating polarization effects in understanding confined electrolytes, both theoretically and computationally.


title: Extended flexibility of Lyapunov exponents for Anosov diffeomorphisms

abstract: Bochi-Katok-Rodriguez Hertz proposed in [BKH21] a program on the flexibility of Lyapunov exponents for conservative Anosov diffeomorphisms, and obtained partial results in this direction. For conservative Anosov diffeomorphisms with strong hyperbolic properties we establish extended flexibility results for their Lyapunov exponents. We give examples of Anosov diffeomorphisms with the strong unstable exponent larger than the strong unstable exponent of the linear part. We also give examples of derived from Anosov diffeomorphisms with the metric entropy entropy larger than the entropy of the linear part. These results rely on a new type of deformation which goes beyond the previous Shub-Wilkinson and Baraviera-Bonatti techniques for conservative systems having some invariant directions. In order to estimate the Lyapunov exponents even after breaking the invariant bundles, we obtain an abstract result (Theorem 3.1) which gives bounds on exponents of some specific cocycles and which can be applied in various other settings. We also include various interesting comments in the appendices: our examples are $\mathcal{C}^2$ robust, the Lyapunov exponents are continuous (with respect to the map) even after breaking the invariant bundles, a similar construction can be obtained for the case of multiple eigenvalues.


title: Systematic electrochemical etching of various metal tips for tunneling spectroscopy and scanning probe microscopy

abstract: Hard point-contact spectroscopy and scanning probe microscopy/spectroscopy are powerful techniques for investigating materials with strong expandability. To support these studies, tips with various physical and chemical properties are required. To ensure the reproducibility of experimental results, the fabrication of tips should be standardized, and a controllable and convenient system should be set up. Here a systematic methodology to fabricate various tips is proposed, involving electrochemical etching reactions. The reaction parameters fall into four categories: solution, power supply, immersion depth, and interruption. An etching system was designed and built so that these parameters could be accurately controlled. With this system, etching parameters for copper, silver, gold, platinum/iridium alloy, tungsten, lead, niobium, iron, nickel, cobalt, and permalloy were explored and standardized. Among these tips, silver and niobium's new recipes were explored and standardized. Optical and scanning electron microscopies were performed to characterize the sharp needles. Relevant point-contact experiments were carried out with an etched silver tip to confirm the suitability of the fabricated tips.


title: Bringing SOUL on sky

abstract: The SOUL project is upgrading the 4 SCAO systems of LBT, pushing the current guide star limits of about 2 magnitudes fainter thanks to Electron Multiplied CCD detector. This improvement will open the NGS SCAO correction to a wider number of scientific cases from high contrast imaging in the visible to extra-galactic source in the NIR. The SOUL systems are today the unique case where pyramid WFS, adaptive secondary and EMCCD are used together. This makes SOUL a pathfinder for most of the ELT SCAO systems like the one of GMT, MICADO and HARMONI of E-ELT, where the same key technologies will be employed. Today we have 3 SOUL systems installed on the telescope in commissioning phase. The 4th system will be installed in a few months. We will present here the results achieved during daytime testing and commissioning nights up to the present date.


title: Low star formation efficiency due to turbulent adiabatic compression in the Taffy bridge

abstract: The Taffy system (UGC 12914/15) consists of two massive spiral galaxies which had a head-on collision about 20 Myr ago. New sensitive, high-resolution CO(1-0) observations of the Taffy system with the IRAM PdBI are presented. About 25% of the total interferometric CO luminosity stems from the bridge region. Assuming a Galactic N(H2)/ICO conversion factor for the galactic disks and a third of this value for the bridge gas, about 10% of the molecular gas mass is located in the bridge region. The giant HII region close to UGC 12915 is located at the northern edge of the high-surface brightness giant molecular cloud association (GMA), which has the highest velocity dispersion among the bridge GMAs. The bridge GMAs are clearly not virialized because of their high velocity dispersion. Three dynamical models are presented and while no single model reproduces all of the observed features, they are all present in at least one of the models. Most of the bridge gas detected in CO does not form stars. We suggest that turbulent adiabatic compression is responsible for the exceptionally high velocity dispersion of the molecular ISM and the suppression of star formation in the Taffy bridge. In this scenario the turbulent velocity dispersion of the largest eddies and turbulent substructures/clouds increase such that giant molecular clouds are no longer in global virial equilibrium. The increase of the virial parameter leads to a decrease of the star formation efficiency. Most of the low-surface density, CO-emitting gas will disperse without forming stars but some of the high-density gas will probably collapse and form dense star clusters, such as the luminous HII region close to UGC 12915. We suggest that globular clusters and super star clusters formed and still form through the gravitational collapse of gas previously compressed by turbulent adiabatic compression during galaxy interactions.


title: Manipulating the anisotropic phase separation in strained VO2 epitaxial films by nanoscale ion-implantation

abstract: Manipulating the strain induced poly-domains and phase transition in correlated oxide material are important for high performance devices fabrication. Though the electronic transport in the strained oxide film at macroscopic scales can be directly measured, the anisotropic electronic state and the controllable phase separation cross the insulator-to-metal transition within nanoscale size are still elusive. Here, we selected VO2 crystal film as a prototypical oxide and achieved the manipulation of anisotropy electronic phase separation via injecting He+ nanobeam into VO2 film at room temperature. In addition, this nanoscale phase separation was directly visualized by infrared near-field imaging measurements, showing the pronounced and unique cR-axis dependent anisotropy on VO2 surface. Our results offered new insights towards understanding the anisotropic nanoscale phase separation in strained metal oxide films.


title: TMDC-based topological nanospaser: single and double threshold behavior

abstract: We theoretically study a topological nanospaser, which consists of a silver nanospheroid and MoS$_2$ monolayer ake of a circular shape. The metal nanospheroid acts as a plasmonic nanoresonator that supports two rotating modes, which are coupled to the corresponding valleys of MoS$_2$. We apply external circularly polarized light that selectively pumps only one of the valleys of MoS$_2$. The generated spaser dynamics strongly depends on the size (radius) of the MoS$_2$ nanoflake. For small radius, the system has only one spasing regime when only chirally-matched plasmon mode is generated, while at larger size of MoS$_2$, depending on the pump intensity, there are two regimes. In one regime, only the chirally-matched plasmon mode is generated, while in the other regime both chirally-matched and chirally-mismatched modes exist. Different regimes of spaser operation have also opposite handedness of the far-field radiated of the spaser system. Such topological nanospaser has potential applications in different areas of infrared spectroscopy, sensing, probing, and biomedical treatment.


title: Byzantine Generals in the Permissionless Setting

abstract: Consensus protocols have traditionally been studied in a setting where all participants are known to each other from the start of the protocol execution. In the parlance of the 'blockchain' literature, this is referred to as the permissioned setting. What differentiates Bitcoin from these previously studied protocols is that it operates in a permissionless setting, i.e. it is a protocol for establishing consensus over an unknown network of participants that anybody can join, with as many identities as they like in any role. The arrival of this new form of protocol brings with it many questions. Beyond Bitcoin, what can we prove about permissionless protocols in a general sense? How does recent work on permissionless protocols in the blockchain literature relate to the well-developed history of research on permissioned protocols in distributed computing? To answer these questions, we describe a formal framework for the analysis of both permissioned and permissionless systems. Our framework allows for "apples-to-apples" comparisons between different categories of protocols and, in turn, the development of theory to formally discuss their relative merits. A major benefit of the framework is that it facilitates the application of a rich history of proofs and techniques in distributed computing to problems in blockchain and the study of permissionless systems. Within our framework, we then address the questions above. We consider the Byzantine Generals Problem as a formalisation of the problem of reaching consensus, and address a programme of research that asks, "Under what adversarial conditions, and for what types of permissionless protocol, is consensus possible?" We prove a number of results for this programme, our main result being that deterministic consensus is not possible for decentralised permissionless protocols. To close, we give a list of eight open questions.


title: Atomistic simulations of surface reactions in ultra-high-temperature ceramics: O2, H2O and CO adsorption and dissociation on ZrB2 (0001) surfaces

abstract: Understanding surface reactivity is crucial in many fields, going from heterogeneous catalysis to materials oxidation and corrosion. In order to decipher the surface reactions of ZrB2 exposed to the harsh environment of aerospace components, the chemical activity of both Zr- and B-surfaces is predicted and compared by using density functional theory and nudged elastic band methods. In particular the adsorption, dissociation and diffusion of O2, CO and H2O are extensively examined through the calculation of surface adsorption energies and reaction pathways. We find the dissociative adsorption of O2 dominating the reactivity of ZrB2 surfaces, while the dissociation of H2O and CO is weakly active on Zr-surfaces, and even less activated on B-terminated ones. Importantly, we discover that the reaction of O2 and H2O can trigger strong surface reconstruction at B-surfaces. Our work thus provides significant insights into the diverse adsorption and reaction mechanisms of ZrB2 surfaces.


title: The Sources of Statistical Bias Series: Simulated Demonstrations to Illustrate the Causes and Effects of Biases in Statistical Estimates

abstract: When teaching and discussing statistical assumptions, our focus is oftentimes placed on how to test and address potential violations rather than the effects of violating assumptions on the estimates produced by our statistical models. The latter represents a potential avenue to help us better understand the impact of researcher degrees of freedom on the statistical estimates we produce. The Violating Assumptions Series is an endeavor I have undertaken to demonstrate the effects of violating assumptions on the estimates produced across various statistical models. The series will review assumptions associated with estimating causal associations, as well as more complicated statistical models including, but not limited to, multilevel models, path models, structural equation models, and Bayesian models. In addition to the primary goal, the series of posts is designed to illustrate how simulations can be used to develop a comprehensive understanding of applied statistics.


title: The ExoGRAVITY project: using single mode interferometry to characterize exoplanets

abstract: Combining adaptive optics and interferometric observations results in a considerable contrast gain compared to single-telescope, extreme AO systems. Taking advantage of this, the ExoGRAVITY project is a survey of known young giant exoplanets located in the range of 0.1'' to 2'' from their stars. The observations provide astrometric data of unprecedented accuracy, being crucial for refining the orbital parameters of planets and illuminating their dynamical histories. Furthermore, GRAVITY will measure non-Keplerian perturbations due to planet-planet interactions in multi-planet systems and measure dynamical masses. Over time, repetitive observations of the exoplanets at medium resolution ($R=500$) will provide a catalogue of K-band spectra of unprecedented quality, for a number of exoplanets. The K-band has the unique properties that it contains many molecular signatures (CO, H$_2$O, CH$_4$, CO$_2$). This allows constraining precisely surface gravity, metallicity, and temperature, if used in conjunction with self-consistent models like Exo-REM. Further, we will use the parameter-retrieval algorithm petitRADTRANS to constrain the C/O ratio of the planets. Ultimately, we plan to produce the first C/O survey of exoplanets, kick-starting the difficult process of linking planetary formation with measured atomic abundances.


title: Gaussian Process Modelling for Improved Resolution in Faraday Depth Reconstruction

abstract: The incomplete sampling of data in complex polarization measurements from radio telescopes negatively affects both the rotation measure (RM) transfer function and the Faraday depth spectra derived from these data. Such gaps in polarization data are mostly caused by flagging of radio frequency interference and their effects worsen as the percentage of missing data increases. In this paper we present a novel method for inferring missing polarization data based on Gaussian processes (GPs). Gaussian processes are stochastic processes that enable us to encode prior knowledge in our models. They also provide a comprehensive way of incorporating and quantifying uncertainties in regression modelling. In addition to providing non-parametric model estimates for missing values, we also demonstrate that Gaussian process modelling can be used for recovering rotation measure values directly from complex polarization data, and that inferring missing polarization data using this probabilistic method improves the resolution of reconstructed Faraday depth spectra.


title: Online detection of failures generated by storage simulator

abstract: Modern large-scale data-farms consist of hundreds of thousands of storage devices that span distributed infrastructure. Devices used in modern data centers (such as controllers, links, SSD- and HDD-disks) can fail due to hardware as well as software problems. Such failures or anomalies can be detected by monitoring the activity of components using machine learning techniques. In order to use these techniques, researchers need plenty of historical data of devices in normal and failure mode for training algorithms. In this work, we challenge two problems: 1) lack of storage data in the methods above by creating a simulator and 2) applying existing online algorithms that can faster detect a failure occurred in one of the components. We created a Go-based (golang) package for simulating the behavior of modern storage infrastructure. The software is based on the discrete-event modeling paradigm and captures the structure and dynamics of high-level storage system building blocks. The package's flexible structure allows us to create a model of a real-world storage system with a configurable number of components. The primary area of interest is exploring the storage machine's behavior under stress testing or exploitation in the medium- or long-term for observing failures of its components. To discover failures in the time series distribution generated by the simulator, we modified a change point detection algorithm that works in online mode. The goal of the change-point detection is to discover differences in time series distribution. This work describes an approach for failure detection in time series data based on direct density ratio estimation via binary classifiers.


title: Commensurations of the outer automorphism group of a universal Coxeter group

abstract: This paper studies the rigidity properties of the abstract commensurator of the outer automorphism group of a universal Coxeter group of rank $n$, which is the free product $W_n$ of $n$ copies of $\mathbb{Z}/2\mathbb{Z}$. We prove that for $n\geq 5$ the natural map $\mathrm{Out}(W_n) \to \mathrm{Comm}(\mathrm{Out}(W_n))$ is an isomorphism and that every isomorphism between finite index subgroups of $\mathrm{Out}(W_n)$ is given by a conjugation by an element of $\mathrm{Out}(W_n)$.


title: Diffuseness parameter as a bottleneck for accurate half-life calculations

abstract: An investigation of the calculated $\alpha$ decay half-lives of super heavy nuclei (SHN) reveals that the diffuseness parameter is a great bottleneck for achieving accurate results and predictions. In particular, when universal proximity function is adopted for nuclear potential, half-life is found to vary significantly and nonlinearly as a function of diffuseness parameter. To overcome this limiting hurdle, a new semiempirical formula for diffuseness that is dependent on charge and neutron numbers is proposed in this work. With the model at hand, half-lives of 218 SHN are computed, for 68 of which there exists available experimental data and 150 of which are predicted. The calculations of half-lives for 68 SHN are compared against experimental data and the calculated data obtained by using deformed Woods-Saxon, deformed Coulomb potentials model, and six semiempirical formulas. The predictions of 150 SHN are compared against the predictions of seven of the current best semiempirical formulas. Calculations of the present study are in good agreement with the experimental half-lives outperforming all but ImSahu semiempirical formula. Moreover, the predictions of our model are consistent with predictions of the semiempirical formulas. We strongly conclude that more attention should be directed toward obtaining accurate diffuseness parameter values for using it in nuclear calculations.


title: Hubs-biased resistance distances on graphs and networks

abstract: We define and study two new kinds of "effective resistances" based on hubs-biased -- hubs-repelling and hubs-attracting -- models of navigating a graph/network. We prove that these effective resistances are squared Euclidean distances between the vertices of a graph. They can be expressed in terms of the Moore-Penrose pseudoinverse of the hubs-biased Laplacian matrices of the graph. We define the analogous of the Kirchhoff indices of the graph based of these resistance distances. We prove several results for the new resistance distances and the Kirchhoff indices based on spectral properties of the corresponding Laplacians. After an intensive computational search we conjecture that the Kirchhoff index based on the hubs-repelling resistance distance is not smaller than that based on the standard resistance distance, and that the last is not smaller than the one based on the hubs-attracting resistance distance. We also observe that in real-world brain and neural systems the efficiency of standard random walk processes is as high as that of hubs-attracting schemes. On the contrary, infrastructures and modular software networks seem to be designed to be navigated by using their hubs.


title: An efficient dynamical low-rank algorithm for the Boltzmann-BGK equation close to the compressible viscous flow regime

abstract: It has recently been demonstrated that dynamical low-rank algorithms can provide robust and efficient approximation to a range of kinetic equations. This is true especially if the solution is close to some asymptotic limit where it is known that the solution is low-rank. A particularly interesting case is the fluid dynamic limit that is commonly obtained in the limit of small Knudsen number. However, in this case the Maxwellian which describes the corresponding equilibrium distribution is not necessarily low-rank; because of this, the methods known in the literature are only applicable to the weakly compressible case. In this paper, we propose an efficient dynamical low-rank integrator that can capture the fluid limit -- the Navier-Stokes equations -- of the Boltzmann-BGK model even in the compressible regime. This is accomplished by writing the solution as $f=Mg$, where $M$ is the Maxwellian and the low-rank approximation is only applied to $g$. To efficiently implement this decomposition within a low-rank framework requires, in the isothermal case, that certain coefficients are evaluated using convolutions, for which fast algorithms are known. Using the proposed decomposition also has the advantage that the rank required to obtain accurate results is significantly reduced compared to the previous state of the art. We demonstrate this by performing a number of numerical experiments and also show that our method is able to capture sharp gradients/shock waves.


title: Chiral effective Lagrangian for excited heavy-light mesons from QCD

abstract: We derive the chiral effective Lagrangian for excited heavy-light mesons from QCD under proper approximations. We focus on the chiral partners with $j_l^P=\frac{3}{2}^+$ and $j_l^P=\frac{3}{2}^-$ which amounts to ($1^+,2^+$) and ($1^-,2^-$) states respectively. The low energy constants including the masses of the chiral partners are calculated. The calculated spectrum for the excited mesons are found roughly consistent with experimental data. In addition, our results indicate that quantum numbers of $B_J(5970)$ can be identified with $1^-$ or $2^-$.


title: Uplink Beam Management for Millimeter Wave Cellular MIMO Systems with Hybrid Beamforming

abstract: Hybrid analog and digital BeamForming (HBF) is one of the enabling transceiver technologies for millimeter Wave (mmWave) Multiple Input Multiple Output (MIMO) systems. This technology offers highly directional communication, which is able to confront the intrinsic characteristics of mmWave signal propagation. However, the small coherence time in mmWave systems, especially under mobility conditions, renders efficient Beam Management (BM) in standalone mmWave communication a very difficult task. In this paper, we consider HBF transceivers with planar antenna panels and design a multi-level beam codebook for the analog beamformer comprising flat top beams with variable widths. These beams exhibit an almost constant array gain for the whole desired angle width, thereby facilitating efficient hierarchical BM. Focusing on the uplink communication, we present a novel beam training algorithm with dynamic beam ordering, which is suitable for the stringent latency requirements of the latest mmWave standard discussions. Our simulation results showcase the latency performance improvement and received signal-to-noise ratio with different variations of the proposed scheme over the optimum beam training scheme based on exhaustive narrow beam search.


title: Deep Reinforcement Learning for Active High Frequency Trading

abstract: We introduce the first end-to-end Deep Reinforcement Learning (DRL) based framework for active high frequency trading in the stock market. We train DRL agents to trade one unit of Intel Corporation stock by employing the Proximal Policy Optimization algorithm. The training is performed on three contiguous months of high frequency Limit Order Book data, of which the last month constitutes the validation data. In order to maximise the signal to noise ratio in the training data, we compose the latter by only selecting training samples with largest price changes. The test is then carried out on the following month of data. Hyperparameters are tuned using the Sequential Model Based Optimization technique. We consider three different state characterizations, which differ in their LOB-based meta-features. Analysing the agents' performances on test data, we argue that the agents are able to create a dynamic representation of the underlying environment. They identify occasional regularities present in the data and exploit them to create long-term profitable trading strategies. Indeed, agents learn trading strategies able to produce stable positive returns in spite of the highly stochastic and non-stationary environment.


title: Conditional Independence Testing in Hilbert Spaces with Applications to Functional Data Analysis

abstract: We study the problem of testing the null hypothesis that X and Y are conditionally independent given Z, where each of X, Y and Z may be functional random variables. This generalises testing the significance of X in a regression model of scalar response Y on functional regressors X and Z. We show however that even in the idealised setting where additionally (X, Y, Z) has a Gaussian distribution, the power of any test cannot exceed its size. Further modelling assumptions are needed and we argue that a convenient way of specifying these is based on choosing methods for regressing each of X and Y on Z. We propose a test statistic involving inner products of the resulting residuals that is simple to compute and calibrate: type I error is controlled uniformly when the in-sample prediction errors are sufficiently small. We show this requirement is met by ridge regression in functional linear model settings without requiring any eigen-spacing conditions or lower bounds on the eigenvalues of the covariance of the functional regressor. We apply our test in constructing confidence intervals for truncation points in truncated functional linear models and testing for edges in a functional graphical model for EEG data.