title: Electric streamers as a nonlinear instability: the model details

abstract: We propose a new approach to unambiguous determination of parameters of positive and negative electric streamer discharges. From hydrodynamic equations, in the assumption of a solution in the shape of a streamer, it is possible to derive several relations between streamer parameters, which form a system of algebraic equations (SAE). Because of the made approximations, the error in the solution of this system is expected to be probably up to a few tens of percent. Solving the SAE allows us to express all streamer parameters in terms of the streamer length $L$, the constant uniform external electric field $E_e$, and the streamer radius. The solutions with different radii are valid solutions of the hydrodynamic equations, and are analogous to the propagation modes of flat-front perturbations with different transverse wavelengths. We interpret the streamer as a nonlinear instability, whose behavior is determined by choosing the radius at which the velocity is maximized, because, as we show, the velocity plays the same role as the exponential growth rate in the case of linear instabilities. Thus, streamer behavior is unambiguously determined by $E_e$ and $L$, in a relatively computationally economical way. In contrast, numerical methods of solving the microscopic equations, such as hydrodynamic simulations, are more computationally demanding, and the preferred solution in them arises automatically because of numerical fluctuations. The calculations for air at sea level conditions produce reasonable values for commonly observed streamer parameters. The calculated positive streamer velocities and negative threshold fields are compatible with experimental measurements. The physical reason for the positive threshold fields is also discussed. A much simplified analytical model (Appendix B) reproduces many of the presented results, at least qualitatively.


title: Learning reduced systems via deep neural networks with memory

abstract: We present a general numerical approach for constructing governing equations for unknown dynamical systems when only data on a subset of the state variables are available. The unknown equations for these observed variables are thus a reduced system of the complete set of state variables. Reduced systems possess memory integrals, based on the well known Mori-Zwanzig (MZ) formulism. Our numerical strategy to recover the reduced system starts by formulating a discrete approximation of the memory integral in the MZ formulation. The resulting unknown approximate MZ equations are of finite dimensional, in the sense that a finite number of past history data are involved. We then present a deep neural network structure that directly incorporates the history terms to produce memory in the network. The approach is suitable for any practical systems with finite memory length. We then use a set of numerical examples to demonstrate the effectiveness of our method.


title: Mass loss in tidally stripped systems; the energy-based truncation method

abstract: The ability to accurately predict the evolution of tidally stripped haloes is important for understanding galaxy formation and testing the properties of dark matter. Most studies of substructure evolution make predictions based on empirical models of tidal mass loss that are calibrated using numerical simulations. This approach can be accurate in the cases considered, but lacks generality and does not provide a physical understanding of the processes involved. Recently, we demonstrated that truncating NFW distribution functions sharply in energy results in density profiles that resemble those of tidally stripped systems, offering a path to constructing physically motivated models of tidal mass loss. In this work, we review calculations of mass loss based on energy truncation alone, and then consider what secondary effects may modulate mass loss beyond this. We find that a combination of dependence on additional orbital parameters and variations in individual particle energies over an orbit results in a less abrupt truncation in energy space as a subhalo loses mass. Combining the energy truncation approach with a simple prediction for the mass-loss rate, we construct a full model of mass loss that can accurately predict the evolution of a subhalo in terms of a single parameter $\eta_{\rm eff}$. This parameter can be fully determined from the initial orbital and halo properties, and does not require calibration with numerical simulations.


title: Cartesian bicategories with choice

abstract: Relational structures are emerging as ubiquitous mathematical machinery in the semantics of open systems of various kinds. Cartesian bicategories are a well-known categorical algebra of relations that has proved especially useful in recent applications. The passage between a category and its bicategory of relations is an important question that has been widely studied for decades. We study an alternative construction that yields a cartesian bicategory of relations. Its behaviour is closely related to the axiom of choice, which itself can be expressed in the language of cartesian bicategories.


title: Probabilistic learning of boolean functions applied to the binary classification problem with categorical covariates

abstract: In this work we cast the problem of binary classification in terms of estimating a partition on Bernoulli data. When the explanatory variables are all categorical, the problem can be modelled using the language of boolean functions. We offer a probabilistic analysis of the problem, and propose two algorithms for learning boolean functions from binary data.


title: Aqueous Nanoclusters Govern Ion Partitioning in Dense Polymer Membranes

abstract: The uptake and sorption of charged molecules by responsive polymer membranes and hydrogels in aqueous solutions is of key importance for the development of soft functional materials. Here we investigate the partitioning of simple monoatomic (Na$^+$, K$^+$, Cs$^+$, Cl$^-$, I$^-$) and one molecular ion (4-nitrophenolate; NP$^-$) within a dense, electroneutral poly($N$-isopropylacrylamide) membrane using explicit-water molecular dynamics simulations. Inside the predominantly hydrophobic environment water distributes in a network of polydisperse water nanoclusters. The average cluster size determines the mean electrostatic self-energy of the simple ions, which preferably reside deeply inside them; we therefore find substantially larger partition ratios $K\simeq\>$10$^{-1}$ than expected from a simple Born picture using a uniform dielectric constant. Despite their irregular shapes we observe that the water clusters possess a universal negative electrostatic potential with respect to their surrounding, as is known for aqueous liquid-vapor interfaces. This potential, which we find concealed in cases of symmetric monoatomic salts, can dramatically impact the transfer free energies of larger charged molecules because of their weak hydration and increased affinity to interfaces. Consequently, and in stark contrast to the simple ions, the molecular ion NP$^-$ can have a partition ratio much larger than unity, $K\simeq\>$10-30 (depending on the cation type) or even $10^3$ in excess of monovalent salt, which explains recent observations of enhanced reaction kinetics of NP$^-$ reduction catalyzed within dense polymer networks. These results also suggest that ionizing a molecule can even enhance the partitioning in a collapsed, rather hydrophobic gel, which strongly challenges the traditional simplistic reasoning.


title: Reconstructing gravitational wave signals from binary black hole mergers with minimal assumptions

abstract: We present a systematic comparison of the binary black hole (BBH) signal waveform reconstructed by two independent and complementary approaches used in LIGO and Virgo source inference: a template-based analysis, and a morphology-independent analysis. We apply the two approaches to real events and to two sets of simulated observations made by adding simulated BBH signals to LIGO and Virgo detector noise. The first set is representative of the 10 BBH events in the first Gravitational Wave Transient Catalog (GWTC-1). The second set is constructed from a population of BBH systems with total mass and signal strength in the ranges that ground based detectors are typically sensitive. We find that the reconstruction quality of the GWTC-1 events is consistent with the results of both sets of simulated signals. We also demonstrate a simulated case where the presence of a mismodelled effect in the observed signal, namely higher order modes, can be identified through the morphology-independent analysis. This study is relevant for currently progressing and future observational runs by LIGO and Virgo.


title: Compactly supported $\mathbb{A}^{1}$-Euler characteristic and the Hochschild complex

abstract: We show the $\mathbb{A}^{1}$-Euler characteristic of a smooth, projective scheme over a characteristic $0$ field is represented by its Hochschild complex together with a canonical bilinear form, and give an exposition of the compactly supported $\mathbb{A}^{1}$-Euler characteristic $\chi^{c}_{\mathbb{A}^{1}}: K_0(\mathbf{Var}_{k}) \to \text{GW}(k)$ from the Grothendieck group of varieties to the Grothendieck--Witt group of bilinear forms. We also provide example computations.


title: Cantor-solus and Cantor-multus Distributions

abstract: The Cantor distribution is obtained from bitstrings; the Cantor-solus distribution (a new name) admits only strings without adjacent 1 bits. We review moments and order statistics associated with these. The Cantor-multus distribution is introduced -- which instead admits only strings without isolated 1 bits -- and more complicated formulas emerge.


title: Adaptive Numerical Simulation of a Phase-field Fracture Model in Mixed Form tested on an L-shaped Specimen with High Poisson Ratios

abstract: This work presents a new adaptive approach for the numerical simulation of a phase-field model for fractures in nearly incompressible solids. In order to cope with locking effects, we use a recently proposed mixed form where we have a hydro-static pressure as additional unknown besides the displacement field and the phase-field variable. To fulfill the fracture irreversibility constraint, we consider a formulation as a variational inequality in the phase-field variable. For adaptive mesh refinement, we use a recently developed residual-type a posteriori error estimator for the phase-field variational inequality which is efficient and reliable, and robust with respect to the phase-field regularization parameter. The proposed model and the adaptive error-based refinement strategy are demonstrated by means of numerical tests derived from the L-shaped panel test, originally developed for concrete. Here, the Poisson's ratio is changed from the standard setting towards the incompressible limit $\nu \to 0.5$.


title: Explained Variation under the Additive Hazards Model

abstract: We study explained variation under the additive hazards regression model for right-censored data. We consider different approaches for developing such a measure, and focus on one that estimates the proportion of variation in the failure time explained by the covariates. We study the properties of the measure both analytically, and through extensive simulations. We apply the measure to a well-known survival data set as well as the linked Surveillance, Epidemiology and End Results (SEER)-Medicare database for prediction of mortality in early-stage prostate cancer patients using high dimensional claims codes.


title: Adversarial Robustness on In- and Out-Distribution Improves Explainability

abstract: Neural networks have led to major improvements in image classification but suffer from being non-robust to adversarial changes, unreliable uncertainty estimates on out-distribution samples and their inscrutable black-box decisions. In this work we propose RATIO, a training procedure for Robustness via Adversarial Training on In- and Out-distribution, which leads to robust models with reliable and robust confidence estimates on the out-distribution. RATIO has similar generative properties to adversarial training so that visual counterfactuals produce class specific features. While adversarial training comes at the price of lower clean accuracy, RATIO achieves state-of-the-art $l_2$-adversarial robustness on CIFAR10 and maintains better clean accuracy.


title: Probing the Possibilities of Ergodicity in the 1D Spin-1/2 XY Chain with Quench Dynamics

abstract: Ergodicity sits at the heart of the connection between statistical mechanics and dynamics of a physical system. By fixing the initial state of the system into the ground state of the Hamiltonian at zero temperature and tuning a control parameter, we consider the occurrence of the ergodicity with quench dynamics in the one-dimensional (1D) spin-1/2 XY model in a transverse magnetic field. The ground-state phase diagram consists of two ferromagnetic and paramagnetic phases. It is known the magnetization in this spin system is non-ergodic. We set up two different experiments as we call them single and double quenches and test the dynamics of the magnetization along the $Z$-axis and the spin-spin correlation function along the $X$-axis which are the order parameters of the zero-temperature phases . Our exact results reveal that for single quenches at zero-temperature, the ergodicity depends on the initial state and the order parameter. Interestingly on the other setup, a double quench on a cyclic path, ergodicity is completely broken for starting from the phase corresponding to the order parameter. Otherwise, it depends on the first quenched point, and the quench time $T$ when the model spent before a second quench in the way back which gives an ability to controlling the ergodicity in the system. Therefore, and contrary to expectations, in the mentioned model the ergodicity can be observed with probing quench dynamics at zero-temperature. Our results provide further insight into the zero-temperature dynamical behavior of quantum systems and their connections to the ergodicity phenomenon.


title: Addressing the Grid-size Sensitivity Issue in Large-eddy Simulations of Stable Boundary Layers

abstract: In this study, we have identified certain fundamental limitations of a mixing length parameterization used in a popular turbulent kinetic energy-based subgrid-scale model. Replacing this parameterization with a more physically realistic one significantly improves the overall quality of the large-eddy simulations (LESs) of stable boundary layers. For the range of grid sizes considered here (specifically, 1 m -- 12.5 m), the revision dramatically reduces the grid-size sensitivity of the simulations. Most importantly, the revised scheme allows us to reliably estimate the first- and second-order statistics of a well-known LES intercomparison case, even with a coarse grid-size of O(10 m).


title: Addressing via N-body simulations the distribution of the satellite tidal debris in the Milky Way environment

abstract: We study the distribution of the Milky Way satellites stellar and dark matter debris. For the first time we address the question of the tidal disruption of satellites in simulations by utilising simultaneously a) a realistic set of orbits extracted from cosmological simulations, b) a three component host galaxy with live halo, disc and bulge components, and c) satellites from hydrodynamical simulations. We analyse the statistical properties of the satellite debris of all massive galaxies reaching the inner Milky Way on a timescale of 2 Gyr. Up to 80$\%$ of the dark matter is stripped from the satellites, while this happens for up to 30$\%$ of their stars. The stellar debris ends mostly in the inner Milky Way halo, whereas the dark matter debris shows a flat mass distribution over the full main halo. The dark matter debris follows a density profile with inner power law index $\alpha_{\rm DM}=-0.66$ and outer index $\beta_{\rm DM}=2.94$, while for stars $\alpha_{*}=-0.44$ and $\beta_{*}=6.17$. In the inner 25 kpc, the distribution of the stellar debris is flatter than that of the dark matter debris and the orientations of their short axes differ significantly. Changing the orientation of the stellar disc by 90$^{\rm{o}}$ has only a minor impact on the distribution of the satellite debris. Our results indicate that the dark matter is more easily stripped than stars from the Milky Way satellites. The structure of the debris is dominated by the satellite orbital properties. The radial profiles, the flattening and the orientation of the stellar and dark matter debris are significantly different, which prevents the prediction of the dark matter distribution from the observed stellar component.


title: Weighted Meta-Learning

abstract: Meta-learning leverages related source tasks to learn an initialization that can be quickly fine-tuned to a target task with limited labeled examples. However, many popular meta-learning algorithms, such as model-agnostic meta-learning (MAML), only assume access to the target samples for fine-tuning. In this work, we provide a general framework for meta-learning based on weighting the loss of different source tasks, where the weights are allowed to depend on the target samples. In this general setting, we provide upper bounds on the distance of the weighted empirical risk of the source tasks and expected target risk in terms of an integral probability metric (IPM) and Rademacher complexity, which apply to a number of meta-learning settings including MAML and a weighted MAML variant. We then develop a learning algorithm based on minimizing the error bound with respect to an empirical IPM, including a weighted MAML algorithm, $\alpha$-MAML. Finally, we demonstrate empirically on several regression problems that our weighted meta-learning algorithm is able to find better initializations than uniformly-weighted meta-learning algorithms, such as MAML.


title: Optimal Local Explainer Aggregation for Interpretable Prediction

abstract: A key challenge for decision makers when incorporating black box machine learned models into practice is being able to understand the predictions provided by these models. One proposed set of methods is training surrogate explainer models which approximate the more complex model. Explainer methods are generally classified as either local or global, depending on what portion of the data space they are purported to explain. The improved coverage of global explainers usually comes at the expense of explainer fidelity. One way of trading off the advantages of both approaches is to aggregate several local explainers into a single explainer model with improved coverage. However, the problem of aggregating these local explainers is computationally challenging, and existing methods only use heuristics to form these aggregations. In this paper we propose a local explainer aggregation method which selects local explainers using non-convex optimization. In contrast to other heuristic methods, we use an integer optimization framework to combine local explainers into a near-global aggregate explainer. Our framework allows a decision-maker to directly tradeoff coverage and fidelity of the resulting aggregation through the parameters of the optimization problem. We also propose a novel local explainer algorithm based on information filtering. We evaluate our algorithmic framework on two healthcare datasets---the Parkinson's Progression Marker Initiative (PPMI) data set and a geriatric mobility dataset---which is motivated by the anticipated need for explainable precision medicine. Our method outperforms existing local explainer aggregation methods in terms of both fidelity and coverage of classification and improves on fidelity over existing global explainer methods, particularly in multi-class settings where state-of-the-art methods achieve 70% and ours achieves 90%.


title: Snowball sampling from graphs

abstract: We develop unbiased strategies to probabilistic T-wave snowball sampling from graphs, where the interest of estimation may concern finite-order subgraphs such as triangles, cycles or stars. Our approaches encompass also the finite-population sampling strategies to multiplicity sampling and adaptive cluster sampling, both of which can be recast as snowball sampling aimed at graph node totals. A general snowball sampling theory offers greater flexibility in terms of scope and efficiency of graph sampling, in addition to the existing random node or edge sampling methods.


title: Accelerating Auxiliary-Field Quantum Monte Carlo Simulations of Solids with Graphical Processing Unit

abstract: We outline how auxiliary-field quantum Monte Carlo (AFQMC) can leverage graphical processing units (GPUs) to accelerate the simulation of solid state sytems. By exploiting conservation of crystal momentum in the one- and two-electron integrals we show how to efficiently formulate the algorithm to best utilize current GPU architectures. We provide a detailed description of different optimization strategies and profile our implementation relative to standard approaches, demonstrating a factor of 40 speed up over a CPU implementation. With this increase in computational power we demonstrate the ability of AFQMC to systematically converge solid state calculations with respect to basis set and system size by computing the cohesive energy of Carbon in the diamond structure to within 0.02 eV of the experimental result.


title: Spectral classification of B stars: The empirical sequence using SDSS-IV/APOGEE near-IR data

abstract: We present a semi-empirical spectral classification scheme for normal B-type stars using near-infrared spectra (1.5-1.7 $\mu$m) from the SDSS APOGEE2-N DR14 database. The main motivation for working with B-type stars is their importance in the evolution of young stellar clusters, however we also take advantage of having a numerous sample (316 stars) of B-type star candidates in APOGEE2-N, for which we also have optical (3600-9100 \angstrom) counterparts from the LAMOST survey. By first obtaining an accurate spectral classification of the sources using the LAMOST DR3 spectra and the canonical spectral classification scheme Gray & Corbally 2009, we found a linear relation between optical spectral types and the equivalent widths of the hydrogen lines of the Brackett series in the APOGEE2-N NIR spectra. This relation extends smoothly from a similar relation for O and early-B stars found by Roman-Lopes et al. (2018). This way, we obtain a catalog of B-type sources with features in both the optical and NIR, and a classification scheme refined down to one spectral sub-class.


title: Bounding the k-rainbow total domination number

abstract: Recently the notion of $k$-rainbow total domination was introduced for a graph $G$, motivated by a desire to reduce the problem of computing the total domination number of the generalized prism $G \Box K_k$ to an integer labeling problem on $G$. In this paper we further demonstrate usefulness of the labeling approach, presenting bounds on the rainbow total domination number in terms of the total domination number, the rainbow domination number and the rainbow total domination number, as well as the usual domination number, where the latter presents a generalization of a result by Goddard and Henning (2018). We establish Vizing-like results for rainbow domination and rainbow total domination. By stating a Vizing-like conjecture for rainbow total domination we present a different viewpoint on Vizing's original conjecture in the case of bipartite graphs.


title: Skellam Type Processes of Order K and Beyond

abstract: In this article, we introduce Skellam process of order k and its running average. We also discuss the time-changed Skellam process of order k. In particular we discuss space-fractional Skellam process and tempered space-fractional Skellam process via time changes in Poisson process by independent stable subordinator and tempered stable subordinator, respectively. We derive the marginal probabilities, Levy measures, governing difference-differential equations of the introduced processes. Our results generalize Skellam process and running average of Poisson process in several directions.


title: A comparison of wave observations in the Arctic marginal ice zone with spectral models

abstract: Increased economic activity and research interest in the Arctic raise the need for better wave forecasts in the marginal ice zone (MIZ). Mathematical and numerical models of wave propagation in sea ice would benefit from more in situ data for validation. This study presents shipborne wave measurements from the MIZ where altimeter readings are corrected for ship motion to obtain estimated single point ocean surface elevation. From the combined measurements, we obtain significant wave height and zero up-crossing period, as well as one-dimensional wave spectra. In addition, we provide spectra and integrated parameters obtained from inertial motion units (IMU) placed on ice floes inside the MIZ. The results are compared with integrated parameters from the WAM-4 spectral wave model over a period of three days in the open ocean. We also compare our measurements outside and inside the MIZ with hindcast data from the new pan-Arctic WAM-3 model and the Wave Watch III model for the European Arctic, which both model wave attenuation in sea ice. A good agreement is found with WAM-4 and WW3 in zero up-crossing period and significant wave height outside the MIZ, where deviations are less than 23%. WAM-3 is on the other hand up to 60% higher than observations. WW3 and WAM-3 are able to estimate the trends for significant wave height and zero up-crossing period inside the MIZ, although the discrepancies with respect to the observations were larger than in the open ocean. Wave damping by sea ice is investigated by looking at the spatial attenuation coefficients. Predicted attenuation coefficients are found to be 72-83% smaller for WW3 and 3-64% larger for WAM-3 compared to observations. Hence, further model tuning is necessary to better estimate wave parameters in the ice.


title: Exploring the nature and synchronicity of early cluster formation in the Large Magellanic Cloud V: Multiple Populations in ancient Globular Clusters

abstract: We examine four ancient Large Magellanic Cloud (LMC) globular clusters (GCs) for evidence of multiple stellar populations using the Advanced Camera for Surveys and Wide Field Camera 3 on the Hubble Space Telescope Programme GO-14164. NGC 1466, NGC 1841, and NGC 2257 all show evidence for a redder, secondary population along the main-sequence. Reticulum does not show evidence for the presence of a redder population, but this GC has the least number of stars and Monte Carlo simulations indicate that the sample of main sequence stars is too small to robustly infer whether a redder population exists in this cluster. The second, redder, population of the other three clusters constitutes $\sim30-40\%$ of the total population along the main-sequence. This brings the total number of ancient LMC GCs with known split or broadened main-sequences to five. However, unlike for Hodge 11 and NGC 2210 (see arXiv:1904.01434), none of the clusters show evidence for multiple populations in the horizontal branch. We also do not find evidence of a second population along the Red Giant Branch (RGB).


title: Experimental and theoretical study of tracer diffusion in a series of (CoCrFeMn)$_{100-x}$Ni$_x$ alloys

abstract: Tracer diffusion of all constituting elements is studied at various temperatures in a series of (CoCrFeMn)$_{100-x}$Ni$_x$ alloys with compositions ranging from pure Ni to the equiatomic CoCrFeMnNi high-entropy alloy. At a given homologous temperature, the measured tracer diffusion coefficients change non-monotonically along the transition from pure Ni to the concentrated alloys and finally to the equiatomic CoCrFeMnNi alloy. This is explained by atomistic Monte-Carlo simulations based on a modified embedded-atom potentials, which reveal that local heterogeneities of the atomic configurations around a vacancy cause correlation effects and induce significant deviations from predictions of the random alloy model.


title: A Graph Theoretic Formula for the Number of Primes $\pi(n)$

abstract: Let PR$[n]$ be the graph whose vertices are $2,3,\ldots,n$ with vertex $v$ adjacent to vertex $w$ if and only if $\gcd(v,w)>1$. It is shown that $\pi(n)$, the the number of primes no more than $n$, equals the Lov\'{a}sz number of this graph. This result suggests new avenues for graph-theoretic investigations of number-theoretic problems.


title: Examples of Fano manifolds with non-pseudoeffective tangent bundle

abstract: Let $X$ be a Fano manifold. While the properties of the anticanonical divisor $-K_X$ and its multiples have been studied by many authors, the positivity of the tangent bundle $T_X$ is much more elusive. We give a complete characterisation of the pseudoeffectivity of $T_X$ for del Pezzo surfaces, hypersurfaces in the projective space and del Pezzo threefolds.


title: A simple mathematical model for the evolution of the corona virus

abstract: The goal of this note is to present a simple mathematical model with two parameters for the number of deaths due to the corona (COVID-19) virus. The model only requires basic knowledge in differential calculus, and can also be understood by pupils attending secondary school. The model can easily be implemented on a computer, and we will illustrate it on the basis of case studies for different countries.


title: Dimensionally Consistent Preconditioning for Saddle-Point Problems

abstract: The preconditioned iterative solution of large-scale saddle-point systems is of great importance in numerous application areas, many of them involving partial differential equations. Robustness with respect to certain problem parameters is often a concern, and it can be addressed by identifying proper scalings of preconditioner building blocks. In this paper, we consider a new perspective to finding effective and robust preconditioners. Our approach is based on the consideration of the natural physical units underlying the respective saddle-point problem. This point of view, which we refer to as dimensional consistency, suggests a natural combination of the parameters intrinsic to the problem. It turns out that the scaling obtained in this way leads to robustness with respect to problem parameters in many relevant cases. As a consequence, we advertise dimensional consistency based preconditioning as a new and systematic way to designing parameter robust preconditoners for saddle-point systems arising from models for physical phenomena.


title: On the Pronormality of Subgroups of Odd Index in some Direct Products of Finite Groups

abstract: A subgroup $H$ of a group $G$ is said to be {\it pronormal} in $G$ if $H$ and $H^g$ are conjugate in $\langle H, H^g \rangle$ for each $g \in G$. Some problems in Finite Group Theory, Combinatorics, and Permutation Group Theory were solved in terms of pronormality, therefore, the question of pronormality of a given subgroup in a given group is of interest. Subgroups of odd index in finite groups satisfy a native necessary condition of pronormality. In this paper we continue investigations on pronormality of subgroups of odd index and consider the pronormality question for subgroups of odd index in some direct products of finite groups. In particular, in this paper we prove that the subgroups of odd index are pronormal in the direct product $G$ of finite simple symplectic groups over fields of odd characteristics if and only if the subgroups of odd index are pronormal in each direct factor of $G$. Moreover, deciding the pronormality of a given subgroup of odd index in the direct product of simple symplectic groups over fields of odd characteristics is reducible to deciding the pronormality of some subgroup $H$ of odd index in a subgroup of $\prod_{i=1}^t \mathbb{Z}_3\wr Sym_{n_i}$, where each $Sym_{n_i}$ acts naturally on $\{1,\dots, n_i\}$, such that $H$ projects onto $\prod_{i=1}^t Sym_{n_i}$. Thus, in this paper we obtain a criterion of pronormality of a subgroup $H$ of odd index in a subgroup of $\prod_{i=1}^t \mathbb{Z}_{p_i}\wr Sym_{n_i}$, where each $p_i$ is a prime and each $Sym_{n_i}$ acts naturally on $\{1,\dots, n_i\}$, such that $H$ projects onto $\prod_{i=1}^t Sym_{n_i}$.


title: Existence and Uniqueness for Nonlinear Integro-Differential Equations in Real Locally Complete Spaces

abstract: We extend existence and uniqueness results of [4] for nonlinear integro-differential equations of Volterra type between real locally complete vector spaces


title: Efficient Oblivious Database Joins

abstract: A major algorithmic challenge in designing applications intended for secure remote execution is ensuring that they are oblivious to their inputs, in the sense that their memory access patterns do not leak sensitive information to the server. This problem is particularly relevant to cloud databases that wish to allow queries over the client's encrypted data. One of the major obstacles to such a goal is the join operator, which is non-trivial to implement obliviously without resorting to generic but inefficient solutions like Oblivious RAM (ORAM). We present an oblivious algorithm for equi-joins which (up to a logarithmic factor) matches the optimal $O(n\log n)$ complexity of the standard non-secure sort-merge join (on inputs producing $O(n)$ outputs). We do not use use expensive primitives like ORAM or rely on unrealistic hardware or security assumptions. Our approach, which is based on sorting networks and novel provably-oblivious constructions, is conceptually simple, easily verifiable, and very efficient in practice. Its data-independent algorithmic structure makes it secure in various different settings for remote computation, even in those that are known to be vulnerable to certain side-channel attacks (such as Intel SGX) or with strict requirements for low circuit complexity (like secure multiparty computation). We confirm that our approach is easily realizable through a compact implementation which matches our expectations for performance and is shown, both formally and empirically, to possess the desired security characteristics.


title: Probing the wavefunctions of correlated states in magic angle graphene

abstract: Using scanning probe microscopy and spectroscopy, we explore the spatial symmetry of the electronic wavefunctions of twisted bilayer graphene at the "magic angle" of 1.1 degrees. This small twist angle leads to a long wavelength moir\'e unit cell on the order of 13 nm and the appearance of two flat bands. As the twist angle is decreased, correlation effects increase until they are maximized at the magic angle. At this angle, the wavefunctions at the charge neutrality point show only C2 symmetry due to the emergence of a charge ordered state. As the system is doped, the symmetry of the wavefunctions change at each commensurate filling of the moir\'e unit cell pointing to the correlated nature of the spin and valley degeneracy broken states.


title: Do Public Datasets Assure Unbiased Comparisons for Registration Evaluation?

abstract: With the increasing availability of new image registration approaches, an unbiased evaluation is becoming more needed so that clinicians can choose the most suitable approaches for their applications. Current evaluations typically use landmarks in manually annotated datasets. As a result, the quality of annotations is crucial for unbiased comparisons. Even though most data providers claim to have quality control over their datasets, an objective third-party screening can be reassuring for intended users. In this study, we use the variogram to screen the manually annotated landmarks in two datasets used to benchmark registration in image-guided neurosurgeries. The variogram provides an intuitive 2D representation of the spatial characteristics of annotated landmarks. Using variograms, we identified potentially problematic cases and had them examined by experienced radiologists. We found that (1) a small number of annotations may have fiducial localization errors; (2) the landmark distribution for some cases is not ideal to offer fair comparisons. If unresolved, both findings could incur bias in registration evaluation.


title: Inversion of the Indefinite Double Covering Map

abstract: Algorithmic methods for the explicit inversion of the indefinite double covering maps are proposed. These are based on either the Givens decomposition or the polar decomposition of the given matrix in the proper, indefinite orthogonal group with signature (p,q). As a by-product we establish that the preimage in the covering group, of a positive matrix in the proper, indefinite orthogonal group, can itself be chosen to be positive definite. Inversion amounts to solving a polynomial system in several variables. Our methods solve this system by either inspection, Groebner bases or by inverting the associated Lie algebra isomorphism and computing certain exponentials explicitly. The last method extends fully constructively for all (p,q) when used together with Givens decompositions. The remaining methods require details of the matrix form of the covering map, but then provide more information about the preimage. The techniques are illustrated for (p, q)in {(2,1), (2,2), (3,2), (4,1)}.


title: A generic ontology and recovery protocols for Human-Robot Collaboration (HRC) systems

abstract: Humans are considered as integral components of Human-Robot Collaboration (HRC) systems, not only as object (e.g. in health care), but also as operators and service providers in manufacturing. Sophisticated and complex tasks are to be collaboratively executed by devices (robots) and humans. We introduce a generic ontology for HRC systems. Description of humans is a part of the ontology. Critical and hazardous (for humans) situations, as well as corresponding safeguards are defined on the basis of the ontology. The ontology is an extension of the ontology introduced in Skarzynski et al. (2018) arXiv:1709.03300. The architecture of SO-MRS (see arXiv:1709.03300), a software platform for automatic task accomplishment, is extended to HRC systems. Ongoing experiments, carried out in a simulated HRC system, are to verify the ontology and the architecture.


title: Towards Metric-Affine Quantum Gravity

abstract: I review here some motivations to consider a theory of gravity based on independent metric and connection, and its status as a quantum theory.


title: A Robotic 3D Perception System for Operating Room Environment Awareness

abstract: Purpose: We describe a 3D multi-view perception system for the da Vinci surgical system to enable Operating room (OR) scene understanding and context awareness. Methods: Our proposed system is comprised of four Time-of-Flight (ToF) cameras rigidly attached to strategic locations on the daVinci Xi patient side cart (PSC). The cameras are registered to the robot's kinematic chain by performing a one-time calibration routine and therefore, information from all cameras can be fused and represented in one common coordinate frame. Based on this architecture, a multi-view 3D scene semantic segmentation algorithm is created to enable recognition of common and salient objects/equipment and surgical activities in a da Vinci OR. Our proposed 3D semantic segmentation method has been trained and validated on a novel densely annotated dataset that has been captured from clinical scenarios. Results: The results show that our proposed architecture has acceptable registration error ($3.3\%\pm1.4\%$ of object-camera distance) and can robustly improve scene segmentation performance (mean Intersection Over Union - mIOU) for less frequently appearing classes ($\ge 0.013$) compared to a single-view method. Conclusion: We present the first dynamic multi-view perception system with a novel segmentation architecture, which can be used as a building block technology for applications such as surgical workflow analysis, automation of surgical sub-tasks and advanced guidance systems.


title: Safe Reinforcement Learning of Control-Affine Systems with Vertex Networks

abstract: This paper focuses on finding reinforcement learning policies for control systems with hard state and action constraints. Despite its success in many domains, reinforcement learning is challenging to apply to problems with hard constraints, especially if both the state variables and actions are constrained. Previous works seeking to ensure constraint satisfaction, or safety, have focused on adding a projection step to a learned policy. Yet, this approach requires solving an optimization problem at every policy execution step, which can lead to significant computational costs. To tackle this problem, this paper proposes a new approach, termed Vertex Networks (VNs), with guarantees on safety during exploration and on learned control policies by incorporating the safety constraints into the policy network architecture. Leveraging the geometric property that all points within a convex set can be represented as the convex combination of its vertices, the proposed algorithm first learns the convex combination weights and then uses these weights along with the pre-calculated vertices to output an action. The output action is guaranteed to be safe by construction. Numerical examples illustrate that the proposed VN algorithm outperforms vanilla reinforcement learning in a variety of benchmark control tasks.


title: T$_{d}$ to 1T$^{\prime}$ structural phase transition in WTe$_{2}$ Weyl semimetal

abstract: Elastic neutron scattering on a single crystal and powder X-ray diffraction measurements were carried out to investigate how the crystal structure evolves as a function of temperature in the Weyl semimetal WTe$_{2}$. A sharp transition from the low-temperature orthorhombic phase (T$_{d}$) to the high-temperature monoclinic phase (1T$^{\prime}$) was observed at ambient pressure in the single crystal near $\sim$565 K. Unlike in MoTe$_{2}$, the solid-solid transition from T$_{d}$ to 1T$^{\prime}$ occurs without the cell doubling of the intermediate T$_{d}^{*}$ phase with AABB (or ABBA) layer stacking. In powders however, the thermal transition from the T$_{d}$ to the 1T$^{\prime}$ phase is broadened and a two phase coexistence was observed until 700K, well above the structural transition.


title: Ergodicity and central limit theorem for random interval homeomorphisms

abstract: The central limit theorem for Markov chains generated by iterated function systems consisting of orientation preserving homeomorphisms of the interval is proved. We study also ergodicity of such systems.


title: Intelligent multiscale simulation based on process-guided composite database

abstract: In the paper, we present an integrated data-driven modeling framework based on process modeling, material homogenization, mechanistic machine learning, and concurrent multiscale simulation. We are interested in the injection-molded short fiber reinforced composites, which have been identified as key material systems in automotive, aerospace, and electronics industries. The molding process induces spatially varying microstructures across various length scales, while the resulting strongly anisotropic and nonlinear material properties are still challenging to be captured by conventional modeling approaches. To prepare the linear elastic training data for our machine learning tasks, Representative Volume Elements (RVE) with different fiber orientations and volume fractions are generated through stochastic reconstruction. More importantly, we utilize the recently proposed Deep Material Network (DMN) to learn the hidden microscale morphologies from data. With essential physics embedded in its building blocks, this data-driven material model can be extrapolated to predict nonlinear material behaviors efficiently and accurately. Through the transfer learning of DMN, we create a unified process-guided material database that covers a full range of geometric descriptors for short fiber reinforced composites. Finally, this unified DMN database is implemented and coupled with macroscale finite element model to enable concurrent multiscale simulations. From our perspective, the proposed framework is also promising in many other emergent multiscale engineering systems, such as additive manufacturing and compressive molding.


title: First spectroscopic study of ionized gas emission lines in the extreme low surface brightness galaxy Malin 1

abstract: Malin 1 is the largest known low surface brightness (LSB) galaxy, the archetype of so-called giant LSBs. The structure and the origin of such galaxies are still poorly understood, especially due to the lack of high-resolution kinematics and spectroscopic data. We use emission lines from spectroscopic observations of Malin 1 aiming to bring new constraints on the internal dynamics and star formation history of Malin 1. We have extracted a total of 16 spectra from different regions of Malin 1 and calculated the rotational velocities of these regions from the wavelength shifts and star formation rates from the observed H$\alpha$ emission line fluxes. We compare our data with existing data and models for Malin 1. For the first time we present the inner rotation curve of Malin 1, characterized in the radial range r < 10 kpc by a steep rise in the rotational velocity up to at least 350 km/s (with a large dispersion), which had not been observed previously. We use these data to study a suite of new mass models for Malin 1. We show that in the inner regions dynamics may be dominated by the stars (although none of our models can explain the highest velocities measured) but that at large radii a massive dark matter halo remains necessary. The H$\alpha$ fluxes derived star formation rates are consistent with an early-type disk for the inner region, and with the level found in extended UV galaxies for the outer parts of the giant disk of Malin 1. We also find signs of high metallicity but low dust content for the inner regions.


title: Design admissibility and de la Garza phenomenon in multi-factor experiments

abstract: The determination of an optimal design for a given regression problem is an intricate optimization problem, especially for models with multivariate predictors. Design admissibility and invariance are main tools to reduce the complexity of the optimization problem and have been successfully applied for models with univariate predictors. In particular several authors have developed sufficient conditions for the existence of saturated designs in univariate models, where the number of support points of the optimal design equals the number of parameters. These results generalize the celebrated de la Garza phenomenon (de la Garza, 1954) which states that for a polynomial regression model of degree $p-1$ any optimal design can be based on at most $p$ points. This paper provides - for the first time - extensions of these results for models with a multivariate predictor. In particular we study a geometric characterization of the support points of an optimal design to provide sufficient conditions for the occurrence of the de la Garza phenomenon in models with multivariate predictors and characterize properties of admissible designs in terms of admissibility of designs in conditional univariate regression models.


title: Hazard recognition in an immersive virtual environment: Framework for the simultaneous analysis of visual search and EEG patterns

abstract: Unmanaged hazards in dangerous construction environments proved to be one of the main sources of injuries and accidents. Hazard recognition is crucial to achieve effective safety management and reduce injuries and fatalities in hazardous job sites. Still, there has been lack of effort that can efficiently assist workers in improving their hazard recognition skills. This study presents virtual safety training in an Immersive Virtual Environment (IVE) to enhance worker's hazard recognition skills. A worker wearing a Virtual Reality (VR) device, that is equipped with an eye-tracker, virtually recognizes hazards on simulated construction sites while a brainwave-sensing device records brain activities. This platform can analyze the overall performance of the workers in a visual hazard recognition task and identify hazards that need additional intervention for each worker. This study provides novel insights on how a worker's brain and eye act simultaneously during a visual hazard recognition process. The presented method can take current safety training programs into another level by providing personalized feedback to the workers.


title: The $\mathcal{S}$-cone and a primal-dual view on second-order representability

abstract: The $\mathcal{S}$-cone provides a common framework for cones of polynomials or exponential sums which establish non-negativity upon the arithmetic-geometric inequality, in particular for sums of non-negative circuit polynomials (SONC) or sums of arithmetic-geometric exponentials (SAGE). In this paper, we study the $\mathcal{S}$-cone and its dual from the viewpoint of second-order representability. Extending results of Averkov and of Wang and Magron on the primal SONC cone, we provide explicit generalized second-order descriptions for rational $\mathcal{S}$-cones and theirs duals.


title: Effects of excitation light polarization on fluorescence emission in two-photon light-sheet microscopy

abstract: Light-sheet microscopy (LSM) is a powerful imaging technique that uses a planar illumination oriented orthogonally to the detection axis. Two-photon (2P) LSM is a variant of LSM that exploits the 2P absorption effect for sample excitation. The light polarization state plays a significant, and often overlooked, role in 2P absorption processes. The scope of this work is to test whether using different polarization states for excitation light can affect the detected signal levels in 2P LSM imaging of typical biological samples with a spatially unordered dye population. Supported by a theoretical model, we compared the fluorescence signals obtained using different polarization states with various fluorophores (fluorescein, EGFP and GCaMP6s) and different samples (liquid solution and fixed or living zebrafish larvae). In all conditions, in agreement with our theoretical expectations, linear polarization oriented parallel to the detection plane provided the largest signal levels, while perpendicularly-oriented polarization gave low fluorescence signal with the biological samples, but a large signal for the fluorescein solution. Finally, circular polarization generally provided lower signal levels. These results highlight the importance of controlling the light polarization state in 2P LSM of biological samples. Furthermore, this characterization represents a useful guide to choose the best light polarization state when maximization of signal levels is needed, e.g. in high-speed 2P LSM.


title: Search for low-mass dark matter with the DAMIC experiment

abstract: The DAMIC (Dark Matter in CCDs) experiment searches for the interactions of dark matter particles with the nuclei and the electrons in the silicon bulk of thick fully depleted charge-coupled devices (CCDs). Because of the low noise and low dark current, DAMIC CCDs are sensitive to the ionization signals expected from low-mass dark matter particles ($< 10$ GeV). A 40-gram target detector has collected data at the SNOLAB underground laboratory since 2017. Recent results from the searches for DM-electron scattering and hidden-photon absorption will be summarized and the status of WIMPs-nucleon search reported. A new detector -- DAMIC-M (DAMIC at Modane) -- with a mass-size of 1 kg and improved CCD readout is under design and will be installed at the underground laboratory of Modane, in France. The current status of DAMIC-M and the near future plans will be presented.


title: Impurity bound states as detectors of topological band structures revisited

abstract: Band structures of topological insulators are characterized by non-local topological invariants. Consequently, proposals for the experimental detection using local probes are rare. A recent paper [Slager et al., Phys. Rev. B 92, 085126 (2015)] has argued, based on theoretical results for a particular class of models, that insulators with topologically trivial and non-trivial band structures in two space dimensions display a qualitatively different response to point-like impurities. Here we present a comprehensive investigation of the impurity response of a large set of models of non-interacting electrons on the honeycomb lattice, driven insulating by either broken inversion, broken time reversal, broken $C_3$, or broken translation symmetry. These cases include Hofstadter bands, strain-induced pseudo-Landau levels and higher-order topological insulators. Our results confirm that for hopping models respecting the lattice symmetries, the response to a single impurity can indeed distinguish between trivial and non-trivial band topology. However, for modulated or inhomogeneous host systems we find that trivial states of matter can display an impurity response akin to that of topologically non-trivial states, and thus the diagnostic fails.


title: Pulsed Neutron Spectroscopy of Low Dimensional Magnets: Past, Present, and Future

abstract: The early 1990s saw the first useful application of pulsed neutron spectroscopy to the study of excitations in low dimensional magnetic systems, with Roger Cowley as a key participant in important early experiments. Since that time the technique has blossomed as a powerful tool utilizing vastly improved neutron instrumentation coupled with more powerful pulsed sources. Here we review representative experiments illustrating some of the fascinating physics that has been revealed in quasi-one and two dimensional systems


title: Scan-less full-field fluorescence-lifetime dual-comb microscopy using two-dimensional spectral mapping and frequency multiplexing of dual-optical-comb beats

abstract: Fluorescence lifetime imaging microscopy (FLIM) is a powerful tool for quantitative fluorescence imaging because fluorescence lifetime is independent of concentration of fluorescent molecules or excitation/detection efficiency and is robust to photobleaching. However, since FLIM is based on point-to-point measurements, mechanical scanning of a focal spot is needed for forming an image, which hampers rapid imaging. In this article, we demonstrate scan-less full-field FLIM based on a one-to-one correspondence between two-dimensional (2D) image pixels and frequency-multiplexed RF signals. A vast number of dual-optical-comb beats between dual optical frequency combs is effectively adopted for 2D spectral mapping and high-density frequency multiplexing in radio-frequency region. Bimodal images of fluorescence amplitude and lifetime are obtained with high quantitativeness from amplitude and phase spectra of fluorescence RF comb modes without the need for mechanical scanning. The proposed method will be useful for rapid quantitative fluorescence imaging in life science.


title: Novel Coronavirus COVID-19 Strike on Arab Countries and Territories: A Situation Report I

abstract: The novel Coronavirus (COVID-19) is an infectious disease caused by a new virus called COVID-19 or 2019-nCoV that first identified in Wuhan, China. The disease causes respiratory illness (such as the flu) with other symptoms such as a cough, fever, and in more severe cases, difficulty breathing. This new Coronavirus seems to be very infectious and has spread quickly and globally. In this work, information about COVID-19 is provided and the situation in Arab countries and territories regarding the COVID-19 strike is presented. The next few weeks main expectations is also given.


title: Connectivity Properties of McKay Quivers

abstract: We present several results regarding the connectivity of McKay quivers of finite-dimensional complex representations of finite groups, with no restriction on the faithfulness or self-duality of the representations. We give examples of McKay quivers, as well as quivers that cannot arise as McKay quivers, and discuss a necessary and sufficient condition for two finite groups to share a connected McKay quiver.


title: Event-Based Control for Online Training of Neural Networks

abstract: Convolutional Neural Network (CNN) has become the most used method for image classification tasks. During its training the learning rate and the gradient are two key factors to tune for influencing the convergence speed of the model. Usual learning rate strategies are time-based i.e. monotonous decay over time. Recent state-of-the-art techniques focus on adaptive gradient algorithms i.e. Adam and its versions. In this paper we consider an online learning scenario and we propose two Event-Based control loops to adjust the learning rate of a classical algorithm E (Exponential)/PD (Proportional Derivative)-Control. The first Event-Based control loop will be implemented to prevent sudden drop of the learning rate when the model is approaching the optimum. The second Event-Based control loop will decide, based on the learning speed, when to switch to the next data batch. Experimental evaluationis provided using two state-of-the-art machine learning image datasets (CIFAR-10 and CIFAR-100). Results show the Event-Based E/PD is better than the original algorithm (higher final accuracy, lower final loss value), and the Double-Event-BasedE/PD can accelerate the training process, save up to 67% training time compared to state-of-the-art algorithms and even result in better performance.


title: Ellipsoidal Subspace Support Vector Data Description

abstract: In this paper, we propose a novel method for transforming data into a low-dimensional space optimized for one-class classification. The proposed method iteratively transforms data into a new subspace optimized for ellipsoidal encapsulation of target class data. We provide both linear and non-linear formulations for the proposed method. The method takes into account the covariance of the data in the subspace; hence, it yields a more generalized solution as compared to Subspace Support Vector Data Description for a hypersphere. We propose different regularization terms expressing the class variance in the projected space. We compare the results with classic and recently proposed one-class classification methods and achieve better results in the majority of cases. The proposed method is also noticed to converge much faster than recently proposed Subspace Support Vector Data Description.


title: A Reliability-aware Multi-armed Bandit Approach to Learn and Select Users in Demand Response

abstract: One challenge in the optimization and control of societal systems is to handle the unknown and uncertain user behavior. This paper focuses on residential demand response (DR) and proposes a closed-loop learning scheme to address these issues. In particular, we consider DR programs where an aggregator calls upon residential users to change their demand so that the total load adjustment is close to a target value. To learn and select the right users, we formulate the DR problem as a combinatorial multi-armed bandit (CMAB) problem with a reliability objective. We propose a learning algorithm: CUCB-Avg (Combinatorial Upper Confidence Bound-Average), which utilizes both upper confidence bounds and sample averages to balance the tradeoff between exploration (learning) and exploitation (selecting). We consider both a fixed time-invariant target and time-varying targets, and show that CUCB-Avg achieves $O(\log T)$ and $O(\sqrt{T \log(T)})$ regrets respectively. Finally, we numerically test our algorithms using synthetic and real data, and demonstrate that our CUCB-Avg performs significantly better than the classic CUCB and also better than Thompson Sampling.


title: On-Axis Tidal Forces in Kerr Spacetime

abstract: Tidal forces are an important feature of General Relativity, which are related to the curvature tensor. We analyze the tidal tensor in Kerr spacetime, with emphasis on the case along the symmetry axis of the Kerr black hole, noting that tidal forces may vanish at a certain point, unlike in the Schwarzschild spacetime, using Boyer-Lindquist coordinates. We study in detail the effects of vanishing tidal forces in a body constituted of dust infalling along the symmetry axis. We find the geodesic deviation equations and solve them, numerically and analytically, for motion along the symmetry axis of the Kerr geometry. We also point out that the intrinsic Gaussian curvature of the event horizon at the symmetry axis is equal to the component along this axis of the tidal tensor there.


title: Space Filling Split Plot Design using Fast Flexible Filling

abstract: In this article, an adaption of an algorithm for the creation of experimental designs by Lekivetz and Jones (2015) is suggested, dealing with constraints around randomization. Split-plot design of experiments is used, when the levels of some factors cannot be modified as easily as others. While most split-plot designs deal in the context of I-optimal or D-optimal designs for continuous response outputs, a space filling design strategy is suggested in here. The proposed designs are evaluated based on different design criteria, as well as an analytical example.


title: Temporal Conjunctive Query Answering in the Extended DL-Lite Family

abstract: Ontology-based query answering (OBQA) augments classical query answering in databases by domain knowledge encoded in an ontology. Systems for OBQA use the ontological knowledge to infer new information that is not explicitly given in the data. Moreover, they usually employ the open-world assumption, which means that knowledge that is not stated explicitly in the data and that is not inferred is not assumed to be true or false. Classical OBQA however considers only a snapshot of the data, which means that information about the temporal evolution of the data is not used for reasoning and hence lost. We investigate temporal conjunctive queries (TCQs) that allow to access temporal data through classical ontologies. In particular, we study combined and data complexity of TCQ entailment for ontologies written in description logics from the extended DL-Lite family. Many of these logics allow for efficient reasoning in the atemporal setting and are successfully applied in practice. We show comprehensive complexity results for temporal reasoning with these logics.


title: Partial Hessian vibrational analysis: vapour pressure isotope effects and their relation with non-covalent interactions

abstract: In this study equilibrium isotope effects (EIEs) on evaporation to several organic solvents (bromobenzene, dibromomethane, ethanol, methanol, and trichloromethane) in the pure phase are estimated employing Kohn-Sham Density functional theory (KS-DFT) along with full Hessian vibrational analysis (FHVA) and partial Hessian vibrational analysis (PHVA). Both FHVA and PHVA qualitatively predict compounds EIEs. Weakly interacting systems (bromobenzene, dibromomethane, and trichloromethane) display identical EIE values for both, FHVA and PHVA. Whereas for strongly interacting systems (ethanol and methanol), PHVA estimates slightly smaller EIE compared to FHVA. By employing the symmetry-adapted perturbation theory (SAPT) along with independent gradient model (IGM), it was possible to establish how minimal changes in compounds interaction energy affected carbon and bromine EIEs estimation in dibromomethane. Moreover, SAPT and IGM revealed how in ethanol, apparently insignificant intermolecular interactions possessed an enormous impact in carbon position specific isotope effects, which highlight the importance of solvent effects for strongly interacting compounds. The presented results have important implications for correctly comprehending and characterising the role of non-covalent interactions in the prediction of isotope effects, as well as providing a robust test to PHVA for the prediction of EIEs.


title: Co-channel Coexistence: Let ITS-G5 and Sidelink C-V2X Make Peace

abstract: In the last few years, two technologies have been developed to enable direct exchange of information between vehicles. These technologies, currently seen as alternatives, are ITS-G5, as commonly referred in Europe, and sidelink LTE-vehicle-to-everything (LTE-V2X) (one of the solutions of the so-called cellular-V2X, C-V2X). For this reason, the attention has been mostly concentrated on comparing them and remarking their strengths and weaknesses to motivate a choice. Differently, in this work we focus on a scenario where both are used in the same area and using the same frequency channels, without the assistance from any infrastructure. Our results show that under co-channel coexistence the range of ITS-G5 is severely degraded, while impact on LTE-V2X is marginal. Additionally, a mitigation method where the CAM data generation is constrained to periodical intervals is shown to reduce the impact of co-channel coexistence, with less degradation on ITS-G5 performance and even improvement for LTE-V2X.


title: Accretion in low-mass members of the Orion Nebula Cluster with young transition disks

abstract: Although the Orion Nebula Cluster is one of the most studied clusters in the solar neighborhood, the evolution of the very low-mass members ($M_* < 0.25 \, M_\odot$) has not been fully addressed due to their faintness. Our goal is to verify if some young and very low-mass objects in the Orion Nebula Cluster show evidence of ongoing accretion using broadband VLT/X-Shooter spectra. For each target, we determined the corresponding stellar parameters, veiling, observed Balmer jump, and accretion rates. Additionally, we searched for the existence of circumstellar disks through available on-line photometry. We detected accretion activity in three young stellar objects in the Orion Nebula Cluster, two of them being in the very low-mass range. We also detected the presence of young transition disks with ages between 1 and 3.5 Myr.


title: Design and optimal control of a tiltrotor micro aerial vehicle for efficient omnidirectional flight

abstract: Omnidirectional micro aerial vehicles are a growing field of research, with demonstrated advantages for aerial interaction and uninhibited observation. While systems with complete pose omnidirectionality and high hover efficiency have been developed independently, a robust system that combines the two has not been demonstrated to date. This paper presents the design and optimal control of a novel omnidirectional vehicle that can exert a wrench in any orientation while maintaining efficient flight configurations. The system design is motivated by the result of a morphology design optimization. A six degrees of freedom optimal controller is derived, with an actuator allocation approach that implements task prioritization, and is robust to singularities. Flight experiments demonstrate and verify the system's capabilities.


title: Chiral SO(4) spin-valley density wave and degenerate topological superconductivity in magic-angle-twisted bilayer-graphene

abstract: Starting from a realistic extended Hubbard model for a $p_{x,y}$-orbital tight-binding model on the Honeycomb lattice, we perform a thorough investigation on the possible electron instabilities in the MA-TBG near the van Hove (VH) dopings. Here we focus on the interplay between the approximate SU(2)$\times$SU(2) symmetry and the $D_3$ symmetry, which leads to intriguing quantum states relevant to recent experiments, as revealed by our systematic RPA based calculations followed by a succeeding mean-field energy minimization for the ground state energy. At the SU(2)$\times$SU(2) symmetric point, the degenerate inter-valley SDW and VDW are mixed into a new state of matter dubbed as the chiral SO(4) spin-valley DW. This state simultaneously hosts three 4-component vectorial spin-valley DW orders with each adopting one wave vector, and the polarization directions of the three DW orders are mutually perpendicular to one another. %in the $\mathbb{R}^4$ space. In the presence of a tiny inter-valley exchange interaction with coefficient $J_H\to 0^{-}$ which breaks the SU(2)$\times$SU(2) symmetry, a pure chiral SDW state is obtained. In the case of $J_H\to 0^{+}$, a nematic VDW+SDW state emerges which possesses a stripy distribution of the charge density, consistent with the recent STM observations. On the aspect of SC, while the triplet $p+ip$ and singlet $d+id$ topological SCs are degenerate at $J_H=0$ near the VH dopings, the former (latter) is favored for $J_H\to 0^{-}$ ($J_H\to 0^{+}$). In addition, the two asymmetric doping-dependent behaviors of the obtained pairing phase diagram are well consistent with experiments.


title: Fast Symmetric Diffeomorphic Image Registration with Convolutional Neural Networks

abstract: Diffeomorphic deformable image registration is crucial in many medical image studies, as it offers unique, special properties including topology preservation and invertibility of the transformation. Recent deep learning-based deformable image registration methods achieve fast image registration by leveraging a convolutional neural network (CNN) to learn the spatial transformation from the synthetic ground truth or the similarity metric. However, these approaches often ignore the topology preservation of the transformation and the smoothness of the transformation which is enforced by a global smoothing energy function alone. Moreover, deep learning-based approaches often estimate the displacement field directly, which cannot guarantee the existence of the inverse transformation. In this paper, we present a novel, efficient unsupervised symmetric image registration method which maximizes the similarity between images within the space of diffeomorphic maps and estimates both forward and inverse transformations simultaneously. We evaluate our method on 3D image registration with a large scale brain image dataset. Our method achieves state-of-the-art registration accuracy and running time while maintaining desirable diffeomorphic properties.


title: A note on Riemann-Liouville fractional Sobolev spaces

abstract: Taking inspiration from a recent paper by Bergounioux, Leaci, Nardi and Tomarelli we study the Riemann-Liouville fractional Sobolev space $W^{s, p}_{RL, a+}(I)$, for $I = (a, b)$ for some $a, b \in \mathbb{R}, a < b$, $s \in (0, 1)$ and $p \in [1, \infty]$; that is, the space of functions $u \in L^{p}(I)$ such that the left Riemann-Liouville $(1 - s)$-fractional integral $I_{a+}^{1 - s}[u]$ belongs to $W^{1, p}(I)$. We prove that the space of functions of bounded variation and the fractional Sobolev space, $BV(I)$ and $W^{s, 1}(I)$, continuously embed into $W^{s, 1}_{RL, a+}(I)$. In addition, we define the space of functions with left Riemann-Liouville $s$-fractional bounded variation, $BV^{s}_{RL,a+}(I)$, as the set of functions $u \in L^{1}(I)$ such that $I^{1 - s}_{a+}[u] \in BV(I)$, and we analyze some fine properties of these functions. Finally, we prove some fractional Sobolev-type embedding results and we analyze the case of higher order Riemann-Liouville fractional derivatives.


title: Active Interaction Force Control for Contact-Based Inspection with a Fully Actuated Aerial Vehicle

abstract: This paper presents and validates active interaction force control and planning for fully actuated and omnidirectional aerial manipulation platforms, with the goal of aerial contact inspection in unstructured environments. We present a variable axis-selective impedance control which integrates direct force control for intentional interaction, using feedback from an on-board force sensor. The control approach aims to reject disturbances in free flight, while handling unintentional interaction, and actively controlling desired interaction forces. A fully actuated and omnidirectional tilt-rotor aerial system is used to show capabilities of the control and planning methods. Experiments demonstrate disturbance rejection, push-and-slide interaction, and force controlled interaction in different flight orientations. The system is validated as a tool for non-destructive testing of concrete infrastructure, and statistical results of


title: Strange Jet Tagging

abstract: Tagging jets of strongly interacting particles initiated by energetic strange quarks is one of the few largely unexplored Standard Model object classification problems remaining in high energy collider physics. In this paper we investigate the purest version of this classification problem in the form of distinguishing strange-quark jets from down-quark jets. Our strategy relies on the fact that a strange-quark jet contains on average a higher ratio of neutral kaon energy to neutral pion energy than does a down-quark jet. Long-lived neutral kaons deposit energy mainly in the hadronic calorimeter of a high energy detector, while neutral pions decay promptly to photons that deposit energy mainly in the electromagnetic calorimeter. In addition, short-lived neutral kaons that decay in flight to charged pion pairs can be identified as a secondary vertex in the inner tracking system. Using these handles we study different approaches to distinguishing strange-quark from down-quark jets, including single variable cut-based methods, a boosted decision tree (BDT) with a small number of simple variables, and a deep learning convolutional neural network (CNN) architecture with jet images. We show that modest gains are possible from the CNN compared with the BDT or a single variable. Starting from jet samples with only strange-quark and down-quark jets, the CNN algorithm can improve the strange to down ratio by a factor of roughly 2 for strange tagging efficiencies below 0.2, and by a factor of 2.5 for strange tagging efficiencies near 0.02.


title: Deep Learning Training in Facebook Data Centers: Design of Scale-up and Scale-out Systems

abstract: Large-scale training is important to ensure high performance and accuracy of machine-learning models. At Facebook we use many different models, including computer vision, video and language models. However, in this paper we focus on the deep learning recommendation models (DLRMs), which are responsible for more than 50% of the training demand in our data centers. Recommendation models present unique challenges in training because they exercise not only compute but also memory capacity as well as memory and network bandwidth. As model size and complexity increase, efficiently scaling training becomes a challenge. To address it we design Zion - Facebook's next-generation large-memory training platform that consists of both CPUs and accelerators. Also, we discuss the design requirements of future scale-out training systems.


title: Higher-dimensional generalizations of the Thouless charge pump

abstract: We define and study analogs of the Thouless charge pump for many-body gapped systems in dimension $D$. We show how to attach a topological invariant to a $D$-dimensional family of such systems, provided all of them have an on-site $U(1)$ symmetry. For a large class of families we argue that this topological invariant is an integer. In the case of gapped systems of free fermions in two dimensions, the invariant can be expressed in terms of the curvature of the Bloch-Berry connection. We also obtain a new formula for the Thouless charge pump in 1d which involves only static linear response and is analogous to the Streda formula for Hall conductivity.


title: TArC: Incrementally and Semi-Automatically Collecting a Tunisian Arabish Corpus

abstract: This article describes the constitution process of the first morpho-syntactically annotated Tunisian Arabish Corpus (TArC). Arabish, also known as Arabizi, is a spontaneous coding of Arabic dialects in Latin characters and arithmographs (numbers used as letters). This code-system was developed by Arabic-speaking users of social media in order to facilitate the writing in the Computer-Mediated Communication (CMC) and text messaging informal frameworks. There is variety in the realization of Arabish amongst dialects, and each Arabish code-system is under-resourced, in the same way as most of the Arabic dialects. In the last few years, the focus on Arabic dialects in the NLP field has considerably increased. Taking this into consideration, TArC will be a useful support for different types of analyses, computational and linguistic, as well as for NLP tools training. In this article we will describe preliminary work on the TArC semi-automatic construction process and some of the first analyses we developed on TArC. In addition, in order to provide a complete overview of the challenges faced during the building process, we will present the main Tunisian dialect characteristics and their encoding in Tunisian Arabish.


title: A deep learning approach for lower back-pain risk prediction during manual lifting

abstract: Occupationally-induced back pain is a leading cause of reduced productivity in industry. Detecting when a worker is lifting incorrectly and at increased risk of back injury presents significant possible benefits. These include increased quality of life for the worker due to lower rates of back injury and fewer workers' compensation claims and missed time for the employer. However, recognizing lifting risk provides a challenge due to typically small datasets and subtle underlying features in accelerometer and gyroscope data. A novel method to classify a lifting dataset using a 2D convolutional neural network (CNN) and no manual feature extraction is proposed in this paper; the dataset consisted of 10 subjects lifting at various relative distances from the body with 720 total trials. The proposed deep CNN displayed greater accuracy (90.6%) compared to an alternative CNN and multilayer perceptron (MLP). A deep CNN could be adapted to classify many other activities that traditionally pose greater challenges in industrial environments due to their size and complexity.


title: Logics of Formal Inconsistency enriched with replacement: an algebraic and modal account

abstract: It is customary to expect from a logical system that it can be algebraizable, in the sense that an algebraic companion of the deductive machinery can always be found. Since the inception of da Costa's paraconsistent calculi $C_n$, algebraic equivalents for such systems have been sought. It is known, however, that these systems are not self-extensional (i.e., they do not satisfy the replacement property). More than this, they are not algebraizable in the sense of Blok-Pigozzi. The same negative results hold for several systems of the hierarchy of paraconsistent logics known as Logics of Formal Inconsistency (LFIs). Because of this, several systems belonging to this class of logics are only characterizable by semantics of a non-deterministic nature. This paper offers a solution for two open problems in the domain of paraconsistency, in particular connected to algebraization of LFIs, by extending with rules several LFIs weaker than $C_1$ , thus obtaining the replacement property (that is, such LFIs turn out to be self-extensional). Moreover, these logics become algebraizable in the standard Lindenbaum-Tarski's sense by a suitable variety of Boolean algebras extended with additional operations. The weakest LFI satisfying replacement presented here is called RmbC, which is obtained from the basic LFI called mbC. Some axiomatic extensions of RmbC are also studied. In addition, a neighborhood semantics is defined for such systems. It is shown that RmbC can be defined within the minimal bimodal non-normal logic E+E defined by the fusion of the non-normal modal logic E with itself. Finally, the framework is extended to first-order languages. RQmbC, the quantified extension of RmbC, is shown to be sound and complete w.r.t. the proposed algebraic semantics.


title: py4DSTEM: a software package for multimodal analysis of four-dimensional scanning transmission electron microscopy datasets

abstract: Scanning transmission electron microscopy (STEM) allows for imaging, diffraction, and spectroscopy of materials on length scales ranging from microns to atoms. By using a high-speed, direct electron detector, it is now possible to record a full 2D image of the diffracted electron beam at each probe position, typically a 2D grid of probe positions. These 4D-STEM datasets are rich in information, including signatures of the local structure, orientation, deformation, electromagnetic fields and other sample-dependent properties. However, extracting this information requires complex analysis pipelines, from data wrangling to calibration to analysis to visualization, all while maintaining robustness against imaging distortions and artifacts. In this paper, we present py4DSTEM, an analysis toolkit for measuring material properties from 4D-STEM datasets, written in the Python language and released with an open source license. We describe the algorithmic steps for dataset calibration and various 4D-STEM property measurements in detail, and present results from several experimental datasets. We have also implemented a simple and universal file format appropriate for electron microscopy data in py4DSTEM, which uses the open source HDF5 standard. We hope this tool will benefit the research community, helps to move the developing standards for data and computational methods in electron microscopy, and invite the community to contribute to this ongoing, fully open-source project.


title: Quasi-compact vortices

abstract: We deal with planar vortex structures in Maxwell-Higgs models in the presence of a generalized magnetic permeability. The model under investigation engenders a real parameter that controls the behavior of the tail of the solutions and of the quantities associated to them. As the parameter gets larger, the solutions attain their boundary values faster, unveiling the existence of a peculiar feature, the presence of double exponential tails. However, the solutions are not compact so we call them quasi-compact vortices.


title: Hardware-Software Codesign for Software Defined Radio: IEEE 802.11p receiver case study

abstract: Software Defined Radio (SDR) platforms are useful tools to design new wireless technologies or to improve specifications of existing ones. The IEEE 802.11p is the de-facto standard for Wireless Vehicular Ad-hoc NETworks (VANETs). It has been implemented on GNU Radio SDR [1], which experiences frames decoding/encoding latency. In this paper an FPGA based SDR is proposed as a solution to accelerate frame decoding of IEEE 802.11p waveforms. We share our experience in designing and validating an SDR on the top of real FPGA based embedded hardware architecture. On the top of an FPGA Zynq platform, we port an SDR GNU Radio following hardware/software (HW/SW) codesign approach. A real time profiling of the SDR system shows that the OFDM equalizer with FFT software function is the most time consuming task of the whole SDR receiver functions (34.74%). We suggest hardware acceleration of the FFT processing function using the free logic gates of the FPGA Zynq as a Hardware (HW) accelerator. The HW OFDM Equalizer with FFT sub-function could accelerate the SDR chain processing in comparison with a pure software implementation. Our FPGA based SDR design is a key step toward proposing an embedded SDR. It would be an alternative for not only the IEEE 802.11p standard but for any OFDM based SDR.


title: An Energy-Aware Online Learning Framework for Resource Management in Heterogeneous Platforms

abstract: Mobile platforms must satisfy the contradictory requirements of fast response time and minimum energy consumption as a function of dynamically changing applications. To address this need, system-on-chips (SoC) that are at the heart of these devices provide a variety of control knobs, such as the number of active cores and their voltage/frequency levels. Controlling these knobs optimally at runtime is challenging for two reasons. First, the large configuration space prohibits exhaustive solutions. Second, control policies designed offline are at best sub-optimal since many potential new applications are unknown at design-time. We address these challenges by proposing an online imitation learning approach. Our key idea is to construct an offline policy and adapt it online to new applications to optimize a given metric (e.g., energy). The proposed methodology leverages the supervision enabled by power-performance models learned at runtime. We demonstrate its effectiveness on a commercial mobile platform with 16 diverse benchmarks. Our approach successfully adapts the control policy to an unknown application after executing less than 25% of its instructions.


title: Predicting Real-Time Locational Marginal Prices: A GAN-Based Video Prediction Approach

abstract: In this paper, we propose an unsupervised data-driven approach to predict real-time locational marginal prices (RTLMPs). The proposed approach is built upon a general data structure for organizing system-wide heterogeneous market data streams into the format of market data images and videos. Leveraging this general data structure, the system-wide RTLMP prediction problem is formulated as a video prediction problem. A video prediction model based on generative adversarial networks (GAN) is proposed to learn the spatio-temporal correlations among historical RTLMPs and predict system-wide RTLMPs for the next hour. An autoregressive moving average (ARMA) calibration method is adopted to improve the prediction accuracy. The proposed RTLMP prediction method takes public market data as inputs, without requiring any confidential information on system topology, model parameters, or market operating details. Case studies using public market data from ISO New England (ISO-NE) and Southwest Power Pool (SPP) demonstrate that the proposed method is able to learn spatio-temporal correlations among RTLMPs and perform accurate RTLMP prediction.


title: The Endomorphisms Algebra of Translations Group and Associative Unitary Ring of Trace-Preserving Endomorphisms in Affine Plane

abstract: This paper introduces a description of Endomorphisms of the translation group in an affine plane, will define the addition and composition of the set of endomorphisms and specify the neutral elements associated with these two actions and present the Endomorphism algebra thereof will distinguish the Trace-preserving endomorphism algebra in affine plane, and prove that the set of Trace-preserving endomorphism associated with the 'addition' action forms a commutative group. We also try to prove that the set of trace-preserving endomorphism, together with the two actions, in it, 'addition' and 'composition' forms an associative and unitary ring.


title: Imagination-Augmented Deep Learning for Goal Recognition

abstract: Being able to infer the goal of people we observe, interact with, or read stories about is one of the hallmarks of human intelligence. A prominent idea in current goal-recognition research is to infer the likelihood of an agent's goal from the estimations of the costs of plans to the different goals the agent might have. Different approaches implement this idea by relying only on handcrafted symbolic representations. Their application to real-world settings is, however, quite limited, mainly because extracting rules for the factors that influence goal-oriented behaviors remains a complicated task. In this paper, we introduce a novel idea of using a symbolic planner to compute plan-cost insights, which augment a deep neural network with an imagination capability, leading to improved goal recognition accuracy in real and synthetic domains compared to a symbolic recognizer or a deep-learning goal recognizer alone.


title: A Framework for Generating Explanations from Temporal Personal Health Data

abstract: Whereas it has become easier for individuals to track their personal health data (e.g., heart rate, step count, food log), there is still a wide chasm between the collection of data and the generation of meaningful explanations to help users better understand what their data means to them. With an increased comprehension of their data, users will be able to act upon the newfound information and work towards striving closer to their health goals. We aim to bridge the gap between data collection and explanation generation by mining the data for interesting behavioral findings that may provide hints about a user's tendencies. Our focus is on improving the explainability of temporal personal health data via a set of informative summary templates, or "protoforms." These protoforms span both evaluation-based summaries that help users evaluate their health goals and pattern-based summaries that explain their implicit behaviors. In addition to individual users, the protoforms we use are also designed for population-level summaries. We apply our approach to generate summaries (both univariate and multivariate) from real user data and show that our system can generate interesting and useful explanations.


title: Faster IVA: Update Rules for Independent Vector Analysis based on Negentropy and the Majorize-Minimize Principle

abstract: Algorithms for Blind Source Separation (BSS) of acoustic signals require efficient and fast converging optimization strategies to adapt to nonstationary signal statistics and time-varying acoustic scenarios. In this paper, we derive fast converging update rules from a negentropy perspective, which are based on the Majorize-Minimize (MM) principle and eigenvalue decomposition. The presented update rules are shown to outperform competing state-of-the-art methods in terms of convergence speed at a comparable runtime due to the restriction to unitary demixing matrices. This is demonstrated by experiments with recorded real-world data.


title: Message complexity of population protocols

abstract: The standard population protocol model assumes that when two agents interact, each observes the entire state of the other agent. We initiate the study of $\textit{message complexity}$ for population protocols, where the state of an agent is divided into an externally-visible $\textit{message}$ and an internal component, where only the message can be observed by the other agent in an interaction. We consider the case of $O(1)$ message complexity. When time is unrestricted, we obtain an exact characterization of the stably computable predicates based on the number of internal states $s(n)$: If $s(n) = o(n)$ then the protocol computes semilinear predicates (unlike the original model, which can compute non-semilinear predicates with $s(n) = O(\log n)$), and otherwise it computes a predicate decidable by a nondeterministic $O(n \log s(n))$-space-bounded Turing machine. We then introduce novel $O(\mathrm{polylog}(n))$ expected time protocols for junta/leader election and general purpose broadcast correct with high probability, and approximate and exact population size counting correct with probability 1. Finally, we show that the main constraint on the power of bounded-message-size protocols is the size of the internal states: with unbounded internal states, any computable function can be computed with probability 1 in the limit by a protocol that uses only $\textit{1-bit}$ messages.


title: Testing Polynomials for Vanishing on Cartesian Products of Planar Point Sets: Collinearity Testing and Related Problems

abstract: We present subquadratic algorithms, in the algebraic decision-tree model of computation, for detecting whether there exists a triple of points, belonging to three respective sets $A$, $B$, and $C$ of points in the plane, that satisfy a certain polynomial equation or two equations. The best known instance of such a problem is testing for the existence of a collinear triple of points in $A\times B\times C$, a classical 3SUM-hard problem that has so far defied any attempt to obtain a subquadratic solution, whether in the (uniform) real RAM model, or in the algebraic decision-tree model. While we are still unable to solve this problem, in full generality, in subquadratic time, we obtain such a solution, in the algebraic decision-tree model, that uses only roughly $O(n^{28/15})$ constant-degree polynomial sign tests, for the special case where two of the sets lie on two respective one-dimensional curves and the third is placed arbitrarily in the plane. Our technique is fairly general, and applies to many other problems where we seek a triple that satisfies a single polynomial equation, e.g., determining whether $A\times B\times C$ contains a triple spanning a unit-area triangle. This result extends recent work by Barba \etal~(2017) and by Chan (2018), where all three sets $A$,~$B$, and~$C$ are assumed to be one-dimensional. As a second application of our technique, we again have three $n$-point sets $A$, $B$, and $C$ in the plane, and we want to determine whether there exists a triple $(a,b,c) \in A\times B\times C$ that simultaneously satisfies two independent real polynomial equations. For example, this is the setup when testing for collinearity in the complex plane, when each of the sets $A$, $B$, $C$ lies on some constant-degree algebraic curve. We show that problems of this kind can be solved with roughly $O(n^{24/13})$ constant-degree polynomial sign tests.


title: Deep Reinforcement Learning with Robust and Smooth Policy

abstract: Deep reinforcement learning (RL) has achieved great empirical successes in various domains. However, the large search space of neural networks requires a large amount of data, which makes the current RL algorithms not sample efficient. Motivated by the fact that many environments with continuous state space have smooth transitions, we propose to learn a smooth policy that behaves smoothly with respect to states. We develop a new framework -- \textbf{S}mooth \textbf{R}egularized \textbf{R}einforcement \textbf{L}earning ($\textbf{SR}^2\textbf{L}$), where the policy is trained with smoothness-inducing regularization. Such regularization effectively constrains the search space, and enforces smoothness in the learned policy. Moreover, our proposed framework can also improve the robustness of policy against measurement error in the state space, and can be naturally extended to distribubutionally robust setting. We apply the proposed framework to both on-policy (TRPO) and off-policy algorithm (DDPG). Through extensive experiments, we demonstrate that our method achieves improved sample efficiency and robustness.


title: Generalized Curie-Weiss-Potts model and quadratic pressure in ergodic theory

abstract: We extend results on quadratic pressure and convergence of Gibbs mesures from previous joined work of the authors to the Curie-Weiss-Potts model. We define the notion of equilibrium state for the quadratic pressure and show that under some conditions on the maxima for some auxiliary function, the Gibbs measure converges to a convex combination of eigen-measures for the Transfer Operator. This extension works for dynamical systems defined by {infinite-to-one} maps. As an example, we compute the equilibrium for {the mean-field} $XY$ model as the number of particles goes to $+\infty$.


title: Mode-space-compatible inelastic scattering in atomistic nonequilibrium Green's function implementations

abstract: The nonequilibrium Green's function (NEGF) method is often used to predict transport in atomistically resolved nanodevices and yields an immense numerical load when inelastic scattering on phonons is included. To ease this load, this work extends the atomistic mode space approach of Ref. [1] to include inelastic scattering on optical and acoustic phonons in silicon nanowires. This work also includes the exact calculation of the real part of retarded scattering self-energies in the reduced basis representation using the Kramers-Kronig relations. The inclusion of the Kramers-Kronig relation for the real part of the retarded scattering self-energy increases the impact of scattering. Virtually perfect agreement with results of the original representation is achieved with matrix rank reductions of more than 97%. Time-to-solution improvements of more than 200$\times$ and peak memory reductions of more than 7$\times$ are shown. This allows for the solution of electron transport scattered on phonons in atomically resolved nanowires with cross-sections larger than 5 nm $\times$ 5 nm.


title: Covering the Relational Join

abstract: In this paper, we initiate a theoretical study of what we call the join covering problem. We are given a natural join query instance $Q$ on $n$ attributes and $m$ relations $(R_i)_{i \in [m]}$. Let $J_{Q} = \ \Join_{i=1}^m R_i$ denote the join output of $Q$. In addition to $Q$, we are given a parameter $\Delta: 1\le \Delta\le n$ and our goal is to compute the smallest subset $\mathcal{T}_{Q, \Delta} \subseteq J_{Q}$ such that every tuple in $J_{Q}$ is within Hamming distance $\Delta - 1$ from some tuple in $\mathcal{T}_{Q, \Delta}$. The join covering problem captures both computing the natural join from database theory and constructing a covering code with covering radius $\Delta - 1$ from coding theory, as special cases. We consider the combinatorial version of the join covering problem, where our goal is to determine the worst-case $|\mathcal{T}_{Q, \Delta}|$ in terms of the structure of $Q$ and value of $\Delta$. One obvious approach to upper bound $|\mathcal{T}_{Q, \Delta}|$ is to exploit a distance property (of Hamming distance) from coding theory and combine it with the worst-case bounds on output size of natural joins (AGM bound hereon) due to Atserias, Grohe and Marx [SIAM J. of Computing'13]. Somewhat surprisingly, this approach is not tight even for the case when the input relations have arity at most two. Instead, we show that using the polymatroid degree-based bound of Abo Khamis, Ngo and Suciu [PODS'17] in place of the AGM bound gives us a tight bound (up to constant factors) on the $|\mathcal{T}_{Q, \Delta}|$ for the arity two case. We prove lower bounds for $|\mathcal{T}_{Q, \Delta}|$ using well-known classes of error-correcting codes e.g, Reed-Solomon codes. We can extend our results for the arity two case to general arity with a polynomial gap between our upper and lower bounds.