title: Largest small polygons: A sequential convex optimization approach

abstract: A small polygon is a polygon of unit diameter. The maximal area of a small polygon with $n=2m$ vertices is not known when $m\ge 7$. Finding the largest small $n$-gon for a given number $n\ge 3$ can be formulated as a nonconvex quadratically constrained quadratic optimization problem. We propose to solve this problem with a sequential convex optimization approach, which is an ascent algorithm guaranteeing convergence to a locally optimal solution. Numerical experiments on polygons with up to $n=128$ sides suggest that the optimal solutions obtained are near-global. Indeed, for even $6 \le n \le 12$, the algorithm proposed in this work converges to known global optimal solutions found in the literature.


title: SwarmCCO: Probabilistic Reactive Collision Avoidance for Quadrotor Swarms under Uncertainty

abstract: We present decentralized collision avoidance algorithms for quadrotor swarms operating under uncertain state estimation. Our approach exploits the differential flatness property and feedforward linearization to approximate the quadrotor dynamics and performs reciprocal collision avoidance. We account for the uncertainty in position and velocity by formulating the collision constraints as chance constraints, which describe a set of velocities that avoid collisions with a specified confidence level. We present two different methods for formulating and solving the chance constraints: our first method assumes a Gaussian noise distribution. Our second method is its extension to the non-Gaussian case by using a Gaussian Mixture Model (GMM). We reformulate the linear chance constraints into equivalent deterministic constraints, which are used with an MPC framework to compute a local collision-free trajectory for each quadrotor. We evaluate the proposed algorithm in simulations on benchmark scenarios and highlight its benefits over prior methods. We observe that both the Gaussian and non-Gaussian methods provide improved collision avoidance performance over the deterministic method. On average, the Gaussian method requires ~5ms to compute a local collision-free trajectory, while our non-Gaussian method is computationally more expensive and requires ~9ms on average in scenarios with 4 agents.


title: A categorical duality for algebras of partial functions

abstract: We prove a categorical duality between a class of abstract algebras of partial functions and a class of (small) topological categories. The algebras are the isomorphs of collections of partial functions closed under the operations of composition, antidomain, range, and preferential union (or 'override'). The topological categories are those whose space of objects is a Stone space, source map is a local homeomorphism, target map is open, and all of whose arrows are epimorphisms.


title: Captum: A unified and generic model interpretability library for PyTorch

abstract: In this paper we introduce a novel, unified, open-source model interpretability library for PyTorch [12]. The library contains generic implementations of a number of gradient and perturbation-based attribution algorithms, also known as feature, neuron and layer importance algorithms, as well as a set of evaluation metrics for these algorithms. It can be used for both classification and non-classification models including graph-structured models built on Neural Networks (NN). In this paper we give a high-level overview of supported attribution algorithms and show how to perform memory-efficient and scalable computations. We emphasize that the three main characteristics of the library are multimodality, extensibility and ease of use. Multimodality supports different modality of inputs such as image, text, audio or video. Extensibility allows adding new algorithms and features. The library is also designed for easy understanding and use. Besides, we also introduce an interactive visualization tool called Captum Insights that is built on top of Captum library and allows sample-based model debugging and visualization using feature importance metrics.


title: Critical properties of the Susceptible-Exposed-Infected model with correlated temporal disorder

abstract: In this paper we study the critical properties of the non-equilibrium phase transition of the Susceptible-Exposed-Infected model under the effects of long-range correlated time-varying environmental noise on the Bethe lattice. We show that temporal noise is perturbatively relevant changing the universality class from the (mean-field) dynamical percolation to the exotic infinite-noise universality class of the contact process model. Our analytical results are based on a mapping to the one-dimensional fractional Brownian motion with an absorbing wall and is confirmed by Monte Carlo simulations. Unlike the contact process, our theory also predicts that it is quite difficult to observe the associated active temporal Griffiths phase in the long-time limit. Finally, we also show an equivalence between the infinite-noise and the compact directed percolation universality classes by relating the SEI model in the presence of temporal disorder to the Domany-Kinzel cellular automaton in the limit of compact clusters.


title: Bayesian phase estimation with adaptive grid refinement

abstract: We introduce a novel Bayesian phase estimation technique based on adaptive grid refinement method. This method automatically chooses the number particles needed for accurate phase estimation using grid refinement and cell merging strategies such that the total number of particles needed at each step is minimal. The proposed method provides a powerful alternative to traditional sampling based sequential Monte Carlo method which tend to fail in certain instances such as when the posterior distribution is bimodal. We also combine grid based and sampling based methods as hybrid particle filter where grid based method can be used to estimate a small but dominant set of parameters and Liu-West (LW) based SMC for the remaining set of parameters. Principal kurtosis analysis can be used to decide the choice of parameters for grid refinement method and for sampling based methods. We provide numerical results comparing the performance of the proposed grid refinement method with Liu-West resampling based SMC. Numerical results suggest that the proposed method is quite promising for quantum phase estimation. It can be easily adapted to Hamiltonian learning which is a very useful technique for estimating unknown parameters of a Hamiltonian and for characterizing unknown quantum devices.


title: Comparison Lift: Bandit-based Experimentation System for Online Advertising

abstract: Comparison Lift is an experimentation-as-a-service (EaaS) application for testing online advertising audiences and creatives at JD.com. Unlike many other EaaS tools that focus primarily on fixed sample A/B testing, Comparison Lift deploys a custom bandit-based experimentation algorithm. The advantages of the bandit-based approach are two-fold. First, it aligns the randomization induced in the test with the advertiser's goals from testing. Second, by adapting experimental design to information acquired during the test, it reduces substantially the cost of experimentation to the advertiser. Since launch in May 2019, Comparison Lift has been utilized in over 1,500 experiments. We estimate that utilization of the product has helped increase click-through rates of participating advertising campaigns by 46% on average. We estimate that the adaptive design in the product has generated 27% more clicks on average during testing compared to a fixed sample A/B design. Both suggest significant value generation and cost savings to advertisers from the product.


title: Strong scattering and parallel guiding of ultracold neutrons

abstract: For ultracold neutrons with a kinetic energy below 10 neV, strong scattering, characterized by $2\pi l_{c} / \lambda\leq 1$, can be obtained in metamaterials of C and $^7$Li. Here $l_{c}$ and $\lambda$ are the coherent scattering mean free path and the neutron wavelength, respectively. UCN interferometry and high-resolution spectroscopy (nano-electronvolt to pico-electronvolt resolution) in parallel waveguide arrays of neutronic metamaterials are given as examples of new experimental possibilities.


title: Braids with the symmetries of Platonic polyhedra in the Coulomb (N+1)-body problem

abstract: We take into account the Coulomb (N + 1)-body problem with N = 12, 24, 60. One of the particles has positive charge Q > 0, and the remaining N have all the same negative charge q < 0. These particles move under the Coulomb force, and the positive charge is assumed to be at rest at the center of mass. Imposing a symmetry constraint, given by the symmetry group of the Platonic polyhedra, we were able to compute periodic orbits, using a shooting method and continuation with respect to the value Q of the positive charge. In the setting of the classical N -body problem, the existence of such orbits is proved with Calculus of Variation techniques, by minimizing the action functional. Here this approach does not seem to work, and numerical computations show that the orbits we compute are not minimizers of the action.


title: Characterization of (semi-)Eberlein compacta using retractional skeletons

abstract: We deeply study retractions associated to suitable models in compact spaces admitting a retractional skeleton and find several interesting consequences. Most importantly, we provide a new characterization of Valdivia compacta using the notion of retractional skeletons, which seems to be helpful when characterizing its subclasses. Further, we characterize Eberlein and semi-Eberlein compacta in terms of retractional skeletons and show that our new characterizations give an alternative proof of the fact that continuous image of an Eberlein compact is Eberlein as well as new stability results for the class of semi-Eberlein compacta, solving in particular an open problem posed by Kubis and Leiderman.


title: Genuine Nonlinearity and its Connection to the Modified Korteweg - de Vries Equation in Phase Dynamics

abstract: The study of hyperbolic waves involves various notions which help characterise how these structures evolve. One important facet is the notion of \emph{genuine nonlinearity}, namely the ability for shocks and rarefactions to form instead of contact discontinuities. In the context of the Whitham Modulation equations, this paper demonstrate that a loss of genuine nonlinearity leads to the appearance of a dispersive set of dynamics in the form of the modified Korteweg de-Vries equation governing the evolution of the waves instead. Its form is universal in the sense that its coefficients can be written entirely using linear properties of the underlying waves such as the conservation laws and linear dispersion relation. This insight is applied to two systems of physical interest, one an optical model and the other a stratified hydrodynamics experiment, to demonstrate how it can be used to provide insight into how waves in these systems evolve when genuine nonlinearity is lost.


title: Instability in large bounded domains -- branched versus unbranched resonances

abstract: We study transitions from convective to absolute instability near a trivial state in large bounded domains for prototypical model problems in the presence of transport and negative nonlinear feedback. We identify two generic scenarios, depending on the nature of the linear mechanism for instability, which both lead to different, universal bifurcation diagrams. In the first, classical case of a linear branched resonance the transition is hard, that is, small changes in a control parameter lead to a finite-size state. In the second, novel case of an unbranched resonance, the transition is gradual. In both cases, the bifurcation diagram is determined by interaction of the leading edge of an invasion front with upstream boundary conditions. Technically, we analyze this interaction in a heteroclinic gluing bifurcation analysis that uses geometric desingularization of the trivial state.


title: Optical spectroscopic observations of low-energy counterparts of Fermi-LAT gamma-ray sources

abstract: A significant fraction of all $\gamma$-ray sources detected by the Large Area Telescope aboard the \fer\ satellite is still lacking a low-energy counterpart. In addition, there is still a large population of $\gamma$-ray sources with associated low-energy counterparts that lack firm classifications. In the last 10 years we have undertaken an optical spectroscopic campaign to address the problem of unassociated/unidentified $\gamma$-ray sources (UGSs), mainly devoted to observing blazars and blazar candidates because they are the largest population of $\gamma$-ray sources associated to date. Here we describe the overall impact of our optical spectroscopic campaign on sources associated in \fer-LAT catalogs, coupled with objects found in the literature. In the literature search, we kept track of efforts by different teams that presented optical spectra of counterparts or potential counterparts of \fer-LAT catalog sources. Our summary includes an analysis of an additional 30 newly-collected optical spectra of counterparts or potential counterparts of \fer-LAT sources of previously unknown nature.New spectra were acquired at the Blanco 4-m and OAN-SPM 2.1-m telescopes, and those available in the Sloan Digital Sky Survey (data release 15) archive. All new sources with optical spectra analyzed here are classified as blazars. Thanks to our campaign, we altogether discovered and classified 394 targets with an additional 123 objects collected from a literature search. We began our optical spectroscopic campaign between the release of the second and third \fer-LAT source catalogs (2FGL and 3FGL, respectively), and classified about 25\% of the sources with uncertain nature and discovered a blazar-like potential counterpart for $\sim$10\% of UGSs listed therein. In the 4FGL catalog, about 350 \fer-LAT sources are classified to date thanks to our campaign. [incomplete abstract]


title: Temporal mode transformations by sequential time and frequency phase modulation for applications in quantum information science

abstract: Controlling the temporal mode shape of quantum light pulses has wide ranging application to quantum information science and technology. Techniques have been developed to control the bandwidth, allow shifting in the time and frequency domains, and perform mode-selective beam-splitter-like transformations. However, there is no present scheme to perform targeted multimode unitary transformations on temporal modes. Here we present a practical approach to realize general transformations for temporal modes. We show theoretically that any unitary transformation on temporal modes can be performed using a series of phase operations in the time and frequency domains. Numerical simulations show that several key transformations on temporal modes can be performed with greater than 95% fidelity using experimentally feasible specifications.


title: Matrix Profile XXII: Exact Discovery of Time Series Motifs under DTW

abstract: Over the last decade, time series motif discovery has emerged as a useful primitive for many downstream analytical tasks, including clustering, classification, rule discovery, segmentation, and summarization. In parallel, there has been an increased understanding that Dynamic Time Warping (DTW) is the best time series similarity measure in a host of settings. Surprisingly however, there has been virtually no work on using DTW to discover motifs. The most obvious explanation of this is the fact that both motif discovery and the use of DTW can be computationally challenging, and the current best mechanisms to address their lethargy are mutually incompatible. In this work, we present the first scalable exact method to discover time series motifs under DTW. Our method automatically performs the best trade-off between time-to-compute and tightness-of-lower-bounds for a novel hierarchy of lower bounds representation we introduce. We show that under realistic settings, our algorithm can admissibly prune up to 99.99% of the DTW computations.


title: Generic spectrum of the weighted Laplacian operator on Cayley graphs

abstract: In this paper, we investigate the spectrum of a class of weighted Laplacians on Cayley graphs and determine under what conditions the corresponding eigenspaces are generically irreducible. Specifically, we analyze the spectrum on left-invariant Cayley graphs endowed with an invariant metric, and we give some criteria for generically irreducible eigenspaces. Additionally, we introduce an operator that is comparable to the Laplacian and show that the same criterion holds.


title: Direct Detection Limits on Heavy Dark Matter

abstract: Multiply-interacting massive particles (MIMPs) are heavy (>10^10 GeV/c^2) dark matter particles that interact strongly with regular matter, but may have evaded detection due to the low number density required to make up the local dark matter halo. These particles could leave track-like signatures in current experiments, similar to lightly-ionizing particles. We show that previously calculated limits from the MAJORANA Demonstrator on the flux of lightly-ionizing particles can be used to exclude MIMP dark matter parameter space up to a mass of 10^15 GeV/c^2. We also calculate limits from the standard XENON1T analysis in this high-mass regime, properly taking into account flux limitations and multi-scatter effects. Finally, we show that a dedicated MIMP analysis using the XENON1T dark matter search could probe unexplored parameter space up to masses of 10^18 GeV/c^2.


title: Characteristic and Necessary Minutiae in Fingerprints

abstract: Fingerprints feature a ridge pattern with moderately varying ridge frequency (RF), following an orientation field (OF), which usually features some singularities. Additionally at some points, called minutiae, ridge lines end or fork and this point pattern is usually used for fingerprint identification and authentication. Whenever the OF features divergent ridge lines (e.g. near singularities), a nearly constant RF necessitates the generation of more ridge lines, originating at minutiae. We call these the necessary minutiae. It turns out that fingerprints feature additional minutiae which occur at rather arbitrary locations. We call these the random minutiae or, since they may convey fingerprint individuality beyond the OF, the characteristic minutiae. In consequence, the minutiae point pattern is assumed to be a realization of the superposition of two stochastic point processes: a Strauss point process (whose activity function is given by the divergence field) with an additional hard core, and a homogeneous Poisson point process, modelling the necessary and the characteristic minutiae, respectively. We perform Bayesian inference using an MCMC-based minutiae separating algorithm (MiSeal). In simulations, it provides good mixing and good estimation of underlying parameters. In application to fingerprints, we can separate the two minutiae patterns and verify by example of two different prints with similar OF that characteristic minutiae convey fingerprint individuality.


title: The multi-returning secretary problem

abstract: In this paper we consider the so-called Multi-returning secretary problem, a version of the Secretary problem in which each candidate has $m$ identical copies. The case $m=2$ has already been completely solved by several authors using different methods both the case $m>2$ had not been satisfactorily solved yet. Here, we provide and efficient algorithm to compute the optimal threshold and the probability of success for every $m$. Moreover, we give a method to determine their asymtoptic values based on the solution of a system of $m$ ODEs.


title: A Monte Carlo simulation of a protein (CoVE) in a matrix of random barriers

abstract: Monte Carlo simulations are performed to study structure and dynamics of a protein CoVE in random media generated by a random distribution of barriers at concentration c with a coarse-grained model in its native (low temperature) and denatured (high temperature) phase. The stochastic dynamics of the protein is diffusive in denature phase at low c, it slows down on increasing c and stops moving beyond a threshold (cth = 0.10). In native phase, the protein moves extremely slow at low c but speeds up on further increasing c in a characteristic range (c = 0.10 - 0.20) before getting trapped at high c (cth = 0.30). The radius of gyration (Rg) of CoVE shows different non-monotonic dependence on c (increase followed by decay) in native and denature phase with a higher and sharper rate of change in farmer. Effective dimension (D) of CoVE is estimated from the scaling of structure factor: in denatured phase, D = 2 (a random coil conformation) at low c (= 0.01 - 0.10) with appearance of some globularization i.e. D ? 2.3, 2.5 at higher c (= 0.2, 0.3). Increasing c seems to reduce the globularity (D = 3) of CoVE in native phase.


title: A structured modified Newton approach for solving systems of nonlinear equations arising in interior-point methods for quadratic programming

abstract: The focus in this work is on interior-point methods for inequality-constrained quadratic programs, and particularly on the system of nonlinear equations to be solved for each value of the barrier parameter. Newton iterations give high quality solutions, but we are interested in modified Newton systems that are computationally less expensive at the expense of lower quality solutions. We propose a structured modified Newton approach where each modified Jacobian is composed of a previous Jacobian, plus one low-rank update matrix per succeeding iteration. Each update matrix is, for a given rank, chosen such that the distance to the Jacobian at the current iterate is minimized, in both 2-norm and Frobenius norm. The approach is structured in the sense that it preserves the nonzero pattern of the Jacobian. The choice of update matrix is supported by results in an ideal theoretical setting. We also produce numerical results with a basic interior-point implementation to investigate the practical performance within and beyond the theoretical framework. In order to improve performance beyond the theoretical framework, we also motivate and construct two heuristics to be added to the method.


title: WarpCore: A Library for fast Hash Tables on GPUs

abstract: Hash tables are ubiquitous. Properties such as an amortized constant time complexity for insertion and querying as well as a compact memory layout make them versatile associative data structures with manifold applications. The rapidly growing amount of data emerging in many fields motivated the need for accelerated hash tables designed for modern parallel architectures. In this work, we exploit the fast memory interface of modern GPUs together with a parallel hashing scheme tailored to improve global memory access patterns, to design WarpCore -- a versatile library of hash table data structures. Unique device-sided operations allow for building high performance data processing pipelines entirely on the GPU. Our implementation achieves up to 1.6 billion inserts and up to 4.3 billion retrievals per second on a single GV100 GPU thereby outperforming the state-of-the-art solutions cuDPP, SlabHash, and NVIDIA RAPIDS cuDF. This performance advantage becomes even more pronounced for high load factors of over $90\%$. To overcome the memory limitation of a single GPU, we scale our approach over a dense NVLink topology which gives us close-to-optimal weak scaling on DGX servers. We further show how WarpCore can be used for accelerating a real world bioinformatics application (metagenomic classification) with speedups of over two orders-of-magnitude against state-of-the-art CPU-based solutions. WC is written in C++/CUDA-C and is openly available at https://github.com/sleeepyjack/warpcore.


title: A semi-analytical solution to the maximum likelihood fit of Poisson data to a linear model using the Cash statistic

abstract: [ABRIDGED] The Cash statistic, also known as the C stat, is commonly used for the analysis of low-count Poisson data, including data with null counts for certain values of the independent variable. The use of this statistic is especially attractive for low-count data that cannot be combined, or re-binned, without loss of resolution. This paper presents a new maximum-likelihood solution for the best-fit parameters of a linear model using the Poisson-based Cash statistic. The solution presented in this paper provides a new and simple method to measure the best-fit parameters of a linear model for any Poisson-based data, including data with null counts. In particular, the method enforces the requirement that the best-fit linear model be non-negative throughout the support of the independent variable. The method is summarized in a simple algorithm to fit Poisson counting data of any size and counting rate with a linear model, by-passing entirely the use of the traditional $\chi^2$ statistic.


title: Causal Bandits without prior knowledge using separating sets

abstract: The Causal Bandit is a variant of the classic Bandit problem where an agent must identify the best action in a sequential decision-making process, where the reward distribution of the actions displays a non-trivial dependence structure that is governed by a causal model. Methods proposed for this problem thus far in the literature rely on exact prior knowledge of the full causal graph. We formulate new causal bandit algorithms that no longer necessarily rely on prior causal knowledge. Instead, they utilize an estimator based on separating sets, which we can find using simple conditional independence tests or causal discovery methods. We show that, given a true separating set, for discrete i.i.d. data, this estimator is unbiased, and has variance which is upper bounded by that of the sample mean. We develop algorithms based on Thompson Sampling and UCB for discrete and Gaussian models respectively and show increased performance on simulation data as well as on a bandit drawing from real-world protein signaling data.


title: Classical particles in the continuum subjected to high density boundary conditions

abstract: We consider a continuous system of classical particles confined in a finite region $\Lambda$ of $\mathbb{R}^d$ interacting through a superstable and tempered pair potential in presence of non free boundary conditions. We prove that the thermodynamic limit of the pressure of the system at any fixed inverse temperature $\beta$ and any fixed fugacity $\lambda$ does not depend on boundary conditions produced by particles outside $\Lambda$ whose density may increase sub-linearly with the distance from the origin at a rate which depends on how fast the pair potential decays at large distances. In particular, if the pair potential $v(x-y)$ is of Lennard-Jones type, i.e. it decays as $C/\|x-y\|^{d+p}$ (with $p>0$) where $\|x-y\|$ is the Euclidean distance between $x$ and $y$, then the existence of the thermodynamic limit of the pressure is guaranteed in presence of boundary conditions generated by external particles which may be distributed with a density increasing with the distance $r$ from the origin as $\rho(1+ r^q)$, where $\rho$ is any positive constant (even arbitrarily larger than the density $\rho_0(\beta,\lambda)$ of the system evaluated with free boundary conditions) and $q\le {1\over 2}\min\{1, p\}$.


title: Graphical Enumeration and Stained Glass Windows, 1: Rectangular Grids

abstract: A survey of enumeration problems arising from the study of graphs formed when the edges of a polygon are marked with evenly spaced points and every pair of points is joined by a line. A few of these problems have been solved, a classical example being the the graph K_n formed when all pairs of vertices of a regular n-gon are joined by chords, which was analyzed by Poonen and Rubinstein in 1998. Most of these problems are unsolved, however, and this two-part article provides data from a number of such problems as well as colored illustrations, which are often reminiscent of stained glass windows. The polygons considered include rectangles, hollow rectangles (or frames), triangles, pentagons, pentagrams, crosses, etc., as well as figures formed by drawing semicircles joining equally-spaced points on a line. %The paper ends with a brief discussion of the problem of how to %design aesthetically pleasing colorings for these graphs. This first part discusses rectangular grids. The 1 X n grids, or equally the graphs K_{n+1,n+1}, were studied by Legendre and Griffiths, and here we investigate the number of cells with a given number of edges and the number of nodes with a given degree. We have only partial results for the m X n rectangles, including upper bounds on the numbers of nodes and cells.


title: A topologically-derived dislocation theory for twist and stretch moir\'e superlattices in bilayer graphene

abstract: We develop a continuum dislocation description of twist and stretch moire superlattices in 2D material bilayers. The continuum formulation is based on the topological constraints introduced by the periodic dislocation network associated with the moire structure. The approach is based on solving analytically for the structural distortion and displacement fields that satisfy the topological constraints, and which minimize the total energy. The total energy is described by both the strain energy of each individual distorted layer, and a Peierls-Nabarro like interfacial contribution arising from stacking disregistry. The dislocation core emerges naturally within the formalism as a result of the competition between the two contributions. The approach presented here captures the structure and energetics of twist and stretch moire superlattices of dislocations with arbitrary direction and character, without assuming an analytical solution a priori, with no adjustable parameters, while accounting naturally for dislocation-dislocation image interactions. In comparisons to atomistic simulations using classical potentials, the maximum structure deviation is 6%, while the maximum line energy deviation is 0.019 eV/Angstrom. Several applications of our model are shown, including predicting the variation of structure with twist angle, and describing dislocation line tension and junction energies.


title: Geometric series expansion of the Neumann-Poincar\'{e} operator: application to composite materials

abstract: The Neumann-Poincar\'{e} operator, a singular integral operator on the boundary of a domain, naturally appears when one solves a conductivity transmission problem via the boundary integral formulation. Recently, a series expression of the Neumann-Poincar\'{e} operator was developed in two dimensions based on geometric function theory. In this paper, we investigate geometric properties of composite materials by using this series expansion. In particular, we obtain explicit formulas for the polarization tensor and the effective conductivity for an inclusion or a periodic array of inclusions of arbitrary shape with extremal conductivity, in terms of the associated exterior conformal mapping. Also, we observe by numerical computations that the spectrum of the Neumann--Poincar\'{e} operator has a monotonic behavior with respect to the shape deformation of the inclusion. Additionally, we derive inequality relations of the coefficients of the Riemann mapping of an arbitrary Lipschitz domain by using the properties of the polarization tensor corresponding to the domain.


title: Variational properties of $\sigma_u$-curvature for closed submanifolds of arbitrary codimension in Riemannian manifolds

abstract: The objet of this paper is the study of variations of a functional whose integrant is the $\sigma_u$-curvature of closed submanifolds of arbitrary codimension in Riemannian manifolds.


title: Nutty Kaluza-Klein dyons revisited

abstract: We extend the previous analysis of (locally) asymptotically flat solutions of Kaluza-Klein (KK) theory by assuming that the dilaton charge is an independent parameter. This corresponds to a general nondegenerate matrix of charges within the geodesic sigma model approach and comes into contact with singular solutions of the four-dimensional Einstein-scalar theory. New features of the degenerate class of solutions, which includes regular KK black holes, are also revealed. Solving the constraint equation, we find three distinct branches of the dilaton charge as a function of the other asymptotic charges, one of which contains the previously known solutions, and the other two, related by electric/magnetic duality, are new and singular. We also investigate whether a super-extreme non-rotating solution in the presence of a Newman-Unti-Tamburino (NUT) charge can become a wormhole, as is the case in Einstein-Maxwell theory. It is shown that the dilaton prevents this possibility, while non-traversable five-dimensional vacuum gravitational wormholes can exist. Finally, we analyze the geodesic structure within the chronosphere around the Misner string of Nutty KK dyons, showing that there are no closed timelike geodesics.


title: Hiding in Plain Sight: A Measurement and Analysis of Kids' Exposure to Malicious URLs on YouTube

abstract: The Internet has become an essential part of children's and adolescents' daily life. Social media platforms are used as educational and entertainment resources on daily bases by young users, leading enormous efforts to ensure their safety when interacting with various social media platforms. In this paper, we investigate the exposure of those users to inappropriate and malicious content in comments posted on YouTube videos targeting this demographic. We collected a large-scale dataset of approximately four million records, and studied the presence of malicious and inappropriate URLs embedded in the comments posted on these videos. Our results show a worrisome number of malicious and inappropriate URLs embedded in comments available for children and young users. In particular, we observe an alarming number of inappropriate and malicious URLs, with a high chance of kids exposure, since the average number of views on videos containing such URLs is 48 million. When using such platforms, children are not only exposed to the material available in the platform, but also to the content of the URLs embedded within the comments. This highlights the importance of monitoring the URLs provided within the comments, limiting the children's exposure to inappropriate content.


title: A deep Chandra survey for faint X-ray sources in the Galactic globular cluster M30, and searches for optical and radio counterparts

abstract: We present a deep ($\sim 330~\mathrm{ks}$) {\it Chandra} survey of the Galactic globular cluster M30 (NGC 7099). Combining the new Cycle 18 with the previous Cycle 3 observations we report a total of 10 new X-ray point sources within the $1.03$ arcmin half-light radius, compiling an extended X-ray catalogue of a total of 23 sources. We incorporate imaging observations by the {\it Hubble Space Telescope} and the {\it Karl G. Jansky Very Large Array} from the MAVERIC survey to search for optical and radio counterparts to the new and old sources. Two X-ray sources are found to have a radio counterpart, including the known millisecond pulsar PSR J2140$-$2310A, the radio position of which also matches a previously reported faint optical counterpart which is slightly redder than the main sequence. We found optical counterparts to $18$ of the $23$ X-ray sources, identifying $2$ new cataclysmic variables (CVs), $5$ new CV candidates, $2$ new candidates of RS CVn type of active binary (AB), and $2$ new candidates of BY Dra type of AB. The remaining unclassified X-ray sources are likely background active galactic nuclei (AGN), as their number is consistent with the expected number of AGN at our X-ray sensitivity. Finally, our analysis of radial profiles of different source classes suggests that bright CVs are more centrally distributed than faint CVs in M30, consistent with other core-collapsed globular clusters.


title: Competitive Ratios for Online Multi-capacity Ridesharing

abstract: In multi-capacity ridesharing, multiple requests (e.g., customers, food items, parcels) with different origin and destination pairs travel in one resource. In recent years, online multi-capacity ridesharing services (i.e., where assignments are made online) like Uber-pool, foodpanda, and on-demand shuttles have become hugely popular in transportation, food delivery, logistics and other domains. This is because multi-capacity ridesharing services benefit all parties involved { the customers (due to lower costs), the drivers (due to higher revenues) and the matching platforms (due to higher revenues per vehicle/resource). Most importantly these services can also help reduce carbon emissions (due to fewer vehicles on roads). Online multi-capacity ridesharing is extremely challenging as the underlying matching graph is no longer bipartite (as in the unit-capacity case) but a tripartite graph with resources (e.g., taxis, cars), requests and request groups (combinations of requests that can travel together). The desired matching between resources and request groups is constrained by the edges between requests and request groups in this tripartite graph (i.e., a request can be part of at most one request group in the final assignment). While there have been myopic heuristic approaches employed for solving the online multi-capacity ridesharing problem, they do not provide any guarantees on the solution quality. To that end, this paper presents the first approach with bounds on the competitive ratio for online multi-capacity ridesharing (when resources rejoin the system at their initial location/depot after serving a group of requests).


title: Gravitational form-factors of the $ \pi $ and $ K $ mesons in QCD sum rules

abstract: Taking into account only the quark part of the energy-momentum tensor, the gravitational form-factors of the $ \pi $ and $ K $ mesons are calculated within the light-cone sum rules method using the distribution amplitudes of these pseudoscalar mesons. The $ Q^2 $-dependence of the relevant form-factors are obtained. Moreover, the mean square mass radii of $ \pi $ and $ K $ mesons are calculated and we find that $ \sqrt{\langle r_{\pi}^2 \rangle} = 0.31 {\rm\ fm} $ and $ \sqrt{\langle r_{K}^2 \rangle} = 0.24 {\rm\ fm} $, respectively. We compare our results for gravitational form-factors with the existing literature results. We obtain that our results on $ \sqrt{\langle r_\pi^2 \rangle} $ are in good agreement with the predictions of NJL model, AdS/QCD model, and other models as well as with the existing experimental data.


title: A Reanalysis of Public Galactic Bulge Gravitational Microlensing Events from OGLE-III and IV

abstract: Modern surveys of gravitational microlensing events have progressed to detecting thousands per year. Surveys are capable of probing Galactic structure, stellar evolution, lens populations, black hole physics, and the nature of dark matter. One of the key avenues for doing this is studying the microlensing Einstein radius crossing time distribution ($t_E$). However, systematics in individual light curves as well as over-simplistic modeling can lead to biased results. To address this, we developed a model to simultaneously handle the microlensing parallax due to Earth's motion, systematic instrumental effects, and unlensed stellar variability with a Gaussian Process model. We used light curves for nearly 10,000 OGLE-III and IV Milky Way bulge microlensing events and fit each with our model. We also developed a forward model approach to infer the timescale distribution by forward modeling from the data rather than using point estimates from individual events. We find that modeling the variability in the baseline removes a source of significant bias in individual events, and previous analyses over-estimated the number of long timescale ($t_E>100$ days) events due to their over simplistic models ignoring parallax effects and stellar variability. We use our fits to identify hundreds of events that are likely black holes.

categories: astro-ph.GA astro-ph.EP astro-ph.SR physics.data-an


title: Insight into Delay Based Reservoir Computing via Eigenvalue Analysis

abstract: In this paper we give a profound insight into the computation capability of delay-based reservoir computing via an eigenvalue analysis. We concentrate on the task-independent memory capacity to quantify the reservoir performance and compare these with the eigenvalue spectrum of the dynamical system. We show that these two quantities are deeply connected, and thus the reservoir computing performance is predictable by analyzing the small signal response of the reservoir. Our results suggest that any dynamical system used as a reservoir can be analyzed in this way. We apply our method exemplarily to a photonic laser system with feedback and compare the numerically computed recall capabilities with the eigenvalue spectrum. Optimal performance is found for a system with the eigenvalues having real parts close to zero and off-resonant imaginary parts.


title: Exploration of Fine-Grained Parallelism for Load Balancing Eager K-truss on GPU and CPU

abstract: In this work we present a performance exploration on Eager K-truss, a linear-algebraic formulation of the K-truss graph algorithm. We address performance issues related to load imbalance of parallel tasks in symmetric, triangular graphs by presenting a fine-grained parallel approach to executing the support computation. This approach also increases available parallelism, making it amenable to GPU execution. We demonstrate our fine-grained parallel approach using implementations in Kokkos and evaluate them on an Intel Skylake CPU and an Nvidia Tesla V100 GPU. Overall, we observe between a 1.261. 48x improvement on the CPU and a 9.97-16.92x improvement on the GPU due to our fine-grained parallel formulation.


title: Correlative Synchrotron X-ray Imaging and Diffraction of Directed Energy Deposition Additive Manufacturing

abstract: The governing mechanistic behaviour of Directed Energy Deposition Additive Manufacturing (DED-AM) is revealed by a combined in situ and operando synchrotron X-ray imaging and diffraction study of a nickel-base superalloy, IN718. Using a unique process replicator, real-space phase-contrast imaging enables quantification of the melt-pool boundary and flow dynamics during solidification. This imaging knowledge informed precise diffraction measurements of temporally resolved microstructural phases during transformation and stress development with a spatial resolution of 100 $\mu$m. The diffraction quantified thermal gradient enabled a dendritic solidification microstructure to be predicted and coupled to the stress orientation and magnitude. The fast cooling rate entirely suppressed the formation of secondary phases or recrystallisation in the solid-state. Upon solidification, the stresses rapidly increase to the yield strength during cooling. This insight, combined with IN718 $'$s large solidification range suggests that the accumulated plasticity exhausts the alloy$'$s ductility, causing liquation cracking. This study has revealed additional fundamental mechanisms governing the formation of highly non-equilibrium microstructures during DED-AM.


title: Free minimal actions of solvable Lie groups which are not affable

abstract: We construct an uncountable family of transversely Cantor laminations of compact spaces defined by free minimal actions of solvable groups, which are not affable and whose orbits are not quasi-isometric to Cayley graphs.


title: On Weak Flexibility in Planar Graphs

abstract: Recently, Dvo\v{r}\'ak, Norin, and Postle introduced flexibility as an extension of list coloring on graphs [JGT 19']. In this new setting, each vertex $v$ in some subset of $V(G)$ has a request for a certain color $r(v)$ in its list of colors $L(v)$. The goal is to find an $L$ coloring satisfying many, but not necessarily all, of the requests. The main studied question is whether there exists a universal constant $\epsilon >0$ such that any graph $G$ in some graph class $\mathcal{C}$ satisfies at least $\epsilon$ proportion of the requests. More formally, for $k > 0$ the goal is to prove that for any graph $G \in \mathcal{C}$ on vertex set $V$, with any list assignment $L$ of size $k$ for each vertex, and for every $R \subseteq V$ and a request vector $(r(v): v\in R, ~r(v) \in L(v))$, there exists an $L$-coloring of $G$ satisfying at least $\epsilon|R|$ requests. If this is true, then $\mathcal{C}$ is called $\epsilon$-flexible for lists of size $k$. Choi et al. [arXiv 20'] introduced the notion of weak flexibility, where $R = V$. We further develop this direction by introducing a tool to handle weak flexibility. We demonstrate this new tool by showing that for every positive integer $b$ there exists $\epsilon(b)>0$ so that the class of planar graphs without $K_4, C_5 , C_6 , C_7, B_b$ is weakly $\epsilon(b)$-flexible for lists of size $4$ (here $K_n$, $C_n$ and $B_n$ are the complete graph, a cycle, and a book on $n$ vertices, respectively). We also show that the class of planar graphs without $K_4, C_5 , C_6 , C_7, B_5$ is $\epsilon$-flexible for lists of size $4$. The results are tight as these graph classes are not even 3-colorable.


title: Stable Surfaces and Free Boundary Marginally Outer Trapped Surfaces

abstract: We explore various notions of stability for surfaces embedded and immersed in spacetimes and initial data sets. The interest in such surfaces lies in their potential to go beyond the variational techniques which often underlie the study of minimal and CMC surfaces. We prove two versions of Christodoulou-Yau estimate for $\mathbf{H}$-stable surfaces, a Cohn-Vossen type inequality for non-compact stable marginally outer trapped surface (MOTS), and a global theorem on the topology of $\mathbf{H}$-stable surfaces. Moreover, we give a definition of capillary stability for MOTS with boundary. This notion of stability leads to an area inequality and a local splitting theorem for free boundary stable MOTS. Finally, we establish an index estimate and a diameter estimate for free boundary MOTS. These are straightforward generalizations of Chen-Fraser-Pang and Carlotto-Franz results for free boundary minimal surfaces, respectively.


title: Computationally Efficient Deep Bayesian Unit-Level Modeling of Survey Data under Informative Sampling for Small Area Estimation

abstract: The topic of deep learning has seen a surge of interest in recent years both within and outside of the field of Statistics. Deep models leverage both nonlinearity and interaction effects to provide superior predictions in many cases when compared to linear or generalized linear models. However, one of the main challenges with deep modeling approaches is quantification of uncertainty. The use of random weight models, such as the popularized "Extreme Learning Machine," offer a potential solution in this regard. In addition to uncertainty quantification, these models are extremely computationally efficient as they do not require optimization through stochastic gradient descent, which is what is typically done for deep learning. We show how the use of random weights in a deep model can fit into a likelihood based framework to allow for uncertainty quantification of the model parameters and any desired estimates. Furthermore, we show how this approach can be used to account for informative sampling of survey data through the use of a pseudo-likelihood. We illustrate the effectiveness of this methodology through simulation and with a real survey data application involving American National Election Studies data.


title: Towards an Objective Metric for the Performance of Exact Triangle Count

abstract: The performance of graph algorithms is often measured in terms of the number of traversed edges per second (TEPS). However, this performance metric is inadequate for a graph operation such as exact triangle counting. In triangle counting, execution times on graphs with a similar number of edges can be distinctly different as demonstrated by results from the past Graph Challenge entries. We discuss the need for an objective performance metric for graph operations and the desired characteristics of such a metric such that it more accurately captures the interactions between the amount of work performed and the capabilities of the hardware on which the code is executed. Using exact triangle counting as an example, we derive a metric that captures how certain techniques employed in many implementations improve performance. We demonstrate that our proposed metric can be used to evaluate and compare multiple approaches for triangle counting, using a SIMD approach as a case study against a scalar baseline.


title: How to marry a star: probabilistic constraints for meaning in context

abstract: In this paper, we derive a notion of 'word meaning in context' that characterizes meaning as both intensional and conceptual. We introduce a framework for specifying local as well as global constraints on word meaning in context, together with their interactions, thus modelling the wide range of lexical shifts and ambiguities observed in utterance interpretation. We represent sentence meaning as a 'situation description system', a probabilistic model which takes utterance understanding to be the mental process of describing to oneself one or more situations that would account for an observed utterance. We show how the system can be implemented in practice, and apply it to examples containing various contextualisation phenomena.


title: Post Quantum Secure Command and Control of Mobile Agents : Inserting quantum-resistant encryption schemes in the Secure Robot Operating System

abstract: The secure command and control (C&C) of mobile agents arises in various settings including unmanned aerial vehicles, single pilot operations in commercial settings, and mobile robots to name a few. As more and more of these applications get integrated into aerospace and defense use cases, the security of the communication channel between the ground station and the mobile agent is of increasing importance. The development of quantum computing devices poses a unique threat to secure communications due to the vulnerability of asymmetric ciphers to Shor's algorithm. Given the active development of new quantum resistant encryption techniques, we report the first integration of post-quantum secure encryption schemes with the robot operating system (ROS) and C&C of mobile agents, in general. We integrate these schemes in the application and network layers, and study the performance of these methods by comparing them to present day security schemes such as the widely used RSA algorithm.


title: Type-augmented Relation Prediction in Knowledge Graphs

abstract: Knowledge graphs (KGs) are of great importance to many real world applications, but they generally suffer from incomplete information in the form of missing relations between entities. Knowledge graph completion (also known as relation prediction) is the task of inferring missing facts given existing ones. Most of the existing work is proposed by maximizing the likelihood of observed instance-level triples. Not much attention, however, is paid to the ontological information, such as type information of entities and relations. In this work, we propose a type-augmented relation prediction (TaRP) method, where we apply both the type information and instance-level information for relation prediction. In particular, type information and instance-level information are encoded as prior probabilities and likelihoods of relations respectively, and are combined by following Bayes' rule. Our proposed TaRP method achieves significantly better performance than state-of-the-art methods on four benchmark datasets: FB15K, FB15K-237, YAGO26K-906, and DB111K-174. In addition, we show that TaRP achieves significantly improved data efficiency. More importantly, the type information extracted from a specific dataset can generalize well to other datasets through the proposed TaRP model.


title: Design of the New Wideband Vivaldi Feed for the HERA Radio-Telescope Phase II

abstract: This paper presents the design of a new dual-polarised Vivaldi feed for the Hydrogen Epoch of Reionization Array (HERA) radio-telescope. This wideband feed has been developed to replace the Phase I dipole feed, and is used to illuminate a 14-m diameter dish. It aims to improve the science capabilities of HERA, by allowing it to characterise the redshifted 21-cm hydrogen signal from the Cosmic Dawn as well as from the Epoch of Reionization. This is achieved by increasing the bandwidth from 100 -- 200 MHz to 50 -- 250 MHz, optimising the time response of the antenna - receiver system, and improving its sensitivity. This new Vivaldi feed is directly fed by a differential front-end module placed inside the circular cavity and connected to the back-end via cables which pass in the middle of the tapered slot. We show that this particular configuration has minimal effects on the radiation pattern and on the system response.


title: Scattering Amplitudes and Soft Theorems in Multi-Flavor Galileon Theories

abstract: In this paper, we initiate the study of multi-flavor Galileon theories using the methods of scattering amplitudes. We explore this topic from different perspectives and extend the techniques employed so far mainly in the single-flavor case. This includes soft theorems, generalized soft theorems with non-trivial right-hand side, Galileon dualities, soft bootstrap and bonus relations. We demonstrate new properties on two examples, the multi-flavor U(N) Galileon and the three-flavor U(2)/U(1) Galileon.


title: Stabilization of Finite-Energy Gottesman-Kitaev-Preskill States

abstract: We introduce a new approach to Gottesman-Kitaev-Preskill (GKP) states that treats their finite-energy version in an exact manner. Based on this analysis, we develop new qubit-oscillator circuits that autonomously stabilize a GKP manifold, correcting errors without relying on qubit measurements. Finally, we show numerically that logical information encoded in GKP states is very robust against typical oscillator noise sources when stabilized by these new circuits.


title: Thermal Properties of Rung Disordered Two-leg Quantum Spin Ladders: Quantum Monte Carlo Study

abstract: A two-leg quenched random bond disordered antiferromagnetic spin$-1/2$ Heisenberg ladder system is investigated by means of stochastic series expansion (SSE) quantum Monte Carlo (QMC) method. Thermal properties of the uniform and staggered susceptibilities, the structure factor, the specific heat and the spin gap are calculated over a large number of random realizations in a wide range of disorder strength. According to our QMC simulation results, the considered system has a special temperature point at which the specific heat take the same value regardless of the strength of the disorder. Moreover, the uniform susceptibility is shown to display the same character except for a small difference in the location of the special point. Finally, the spin gap values are found to decrease with increasing disorder parameter and the smallest gap value found in this study is well above the weak coupling limit of the clean case.


title: An analysis of deep neural networks for predicting trends in time series data

abstract: Recently, a hybrid Deep Neural Network (DNN) algorithm, TreNet was proposed for predicting trends in time series data. While TreNet was shown to have superior performance for trend prediction to other DNN and traditional ML approaches, the validation method used did not take into account the sequential nature of time series data sets and did not deal with model update. In this research we replicated the TreNet experiments on the same data sets using a walk-forward validation method and tested our optimal model over multiple independent runs to evaluate model stability. We compared the performance of the hybrid TreNet algorithm, on four data sets to vanilla DNN algorithms that take in point data, and also to traditional ML algorithms. We found that in general TreNet still performs better than the vanilla DNN models, but not on all data sets as reported in the original TreNet study. This study highlights the importance of using an appropriate validation method and evaluating model stability for evaluating and developing machine learning models for trend prediction in time series data.


title: A Magnus extension for locally indicable groups

abstract: A group $G$ possesses the Magnus property if for every two elements $u$, $v \in G$ with the same normal closure, $u$ is conjugate to $v$ or $v^{-1}$. O. Bogopolski and J. Howie proved independently that the fundamental groups of all closed orientable surfaces possess the Magnus property. The analogous result for closed non-orientable surfaces was proved by O. Bogopolski and K. Sviridov except for one case that was later covered by the author. In this article, we generalize those results, which can be viewed as Magnus extensions for free groups, to a Magnus extension for locally indicable groups and consider the influence of adding a group as a direct factor. For this purpose, we also prove versions of the Freiheitssatz for locally indicable groups and of a result by M. Edjvet adding a group as a direct factor.


title: Search for high-redshift blazars with Fermi/LAT

abstract: High-$z$ blazars (z $\geq 2.5$) are the most powerful class of persistent $\gamma$-ray sources in the Universe. These objects possess the highest jet powers and luminosities and have black hole masses often in excess of $10^9$ solar masses. In addition, high-$z$ blazars are important cosmological probes and serve as test objects for blazar evolution models. Due to their large distance, their high-energy emission typically peaks below the GeV range, which makes them difficult to study with Fermi/LAT. Therefore, only the very brightest objects are detectable and, to date, only a small number of high-z blazars have been detected with Fermi/LAT. In this work, we studied the monthly binned long-term $\gamma$-ray emission of a sample of 176 radio and optically detected blazars that have not been reported as known $\gamma$-ray sources in the 3FGL catalog. In order to account for false-positive detections, we calculated monthly Fermi/LAT light curves for a large sample of blank sky positions and derived the number of random fluctuations that we expect at various test statistic (TS) levels. For a given blazar, a detection of TS > 9 in at least one month is expected $\sim 15\%$ of the time. Although this rate is too high to secure detection of an individual source, half of our sample shows such single-month $\gamma$-ray activity, indicating a population of high-energy blazars at distances of up to z=5.2. Multiple TS > 9 monthly detections are unlikely to happen by chance, and we have detected several individual new sources in this way, including the most distant $\gamma$-ray blazar, BZQ J1430+4204 (z = 4.72). Finally, two new $\gamma$-ray blazars at redshifts of z = 3.63 and z = 3.11 are unambiguously detected via very significant (TS > 25) flares in individual monthly time bins.


title: A new partial geometry pg(5,5,2)

abstract: We construct a new partial geometry with parameters pg(5,5,2), not isomorphic to the partial geometry of van Lint and Schrijver.


title: Using Machine Learning and Alternative Data to Predict Movements in Market Risk

abstract: Using machine learning and alternative data for the prediction of financial markets has been a popular topic in recent years. Many financial variables such as stock price, historical volatility and trade volume have already been through extensive investigation. Remarkably, we found no existing research on the prediction of an asset's market implied volatility within this context. This forward-looking measure gauges the sentiment on the future volatility of an asset, and is deemed one of the most important parameters in the world of derivatives. The ability to predict this statistic may therefore provide a competitive edge to practitioners of market making and asset management alike. Consequently, in this paper we investigate Google News statistics and Wikipedia site traffic as alternative data sources to quantitative market data and consider Logistic Regression, Support Vector Machines and AdaBoost as machine learning models. We show that movements in market implied volatility can indeed be predicted through the help of machine learning techniques. Although the employed alternative data appears to not enhance predictive accuracy, we reveal preliminary evidence of non-linear relationships between features obtained from Wikipedia page traffic and movements in market implied volatility.


title: Unconventional magnon excitation by off-resonant microwaves

abstract: It is widely recognized that a physical system can only respond to a periodic driving significantly when the driving frequency matches the normal mode frequency of the system, which leads to resonance. Off-resonant phenomena are rarely considered because of the difficulty to realize strong coupling between physical systems and off-resonant waves. Here we examine the response of a magnetic system to squeezed light and surprisingly find that the magnons are maximally excited when the effective driving frequency is several orders of magnitude larger than the resonant frequency. The generated magnons are squeezed which brings the advantage of tunable squeezing through an external magnetic field. Furthermore, we demonstrate that such off-resonant quasi-particle excitation is universal in all the hybrid systems in which the coherent and parametric interaction of bosons exists and that it is purely a quantum effect, which is rooted in the quantum fluctuations of particles in the squeezed vacuum. Our findings may provide an unconventional route to study off-resonant phenomena and may further benefit the use of hybrid matter-light systems in continuous variable quantum information.


title: Light-matter interaction in open cavities with dielectric stacks

abstract: We evaluate the exact dipole coupling strength between a single emitter and the radiation field within an optical cavity, taking into account the effects of multilayer dielectric mirrors. Our model allows one to freely vary the resonance frequency of the cavity, the frequency of light or atomic transition addressing it and the design wavelength of the dielectric mirror. The coupling strength is derived for an open system with unbound frequency modes. For very short cavities, the effective length used to determine their mode volume and the lengths defining their resonances are different, and also found to diverge appreciably from their geometric length, with the radiation field being strongest within the dielectric mirror itself. Only for cavities much longer than their resonant wavelength does the mode volume asymptotically approach that normally assumed from their geometric length.


title: Cocycle Invariants and Oriented Singular Knots

abstract: We extend the quandle cocycle invariant to oriented singular knots and links using algebraic structures called \emph{oriented singquandles} and assigning weight functions at both regular and singular crossings. This invariant coincides with the classical cocycle invariant for classical knots but provides extra information about singular knots and links. The new invariant distinguishes the singular granny knot from the singular square knot.


title: Brain-Computer Interfaces and the Dangers of Neurocapitalism

abstract: We review how existing trends are relevant to the discussion of brain-computer interfaces and the data they would generate. Then, we posit how the commerce of neural data, dubbed Neurocapitalism, could be impacted by the maturation of brain-computer interface technology. We explore how this could pose fundamental changes to our way of interacting, as well as our sense of autonomy and identity. Because of the power inherent in the technology, and its potentially ruinous abuses, action must be taken before the appearance of the technology, and not come as a reaction to it. The widespread adoption of brain-computer interface technology will certainly change our way of life. Whether it is changed for the better or worse, depends on how well we prepare for its arrival.


title: A stationary solution for the mixed turbulence shell model

abstract: The aim of this work is to prove an existence result on the mixed shell model extending the classic standard existence results from $\ell^2$ initial conditions to $\mu$-almost every initial conditions, where $\mu$ is a Gaussian measure on the infinite dimensional space of initial conditions. A similar result is also shown to hold for turbulence models where a dyadic tree structure replaces the linear one.


title: Improved decoupling for the parabola

abstract: We prove an $(l^2, l^6)$ decoupling inequality for the parabola with constant $(\log R)^c$. In the appendix, we present an application to the six-order correlation of the integer solutions to $x^2+y^2=m$.


title: Detecting and monitoring long-term landslides in urbanized areas with nighttime light data and multi-seasonal Landsat imagery across Taiwan from 1998 to 2017

abstract: Monitoring long-term landslide activity is important for risk assessment and land management. Despite the widespread use of open-access 30m Landsat imagery, their utility for landslide detection is often limited when separating landslides from other anthropogenic disturbances. Here, we produce landslide maps retrospectively from 1998 to 2017 for landslide-prone and highly populated Taiwan (35,874 km2). To improve classification accuracy of landslides, we integrate nighttime light imagery from the Defense Meteorological Satellite Program (DMSP) and the Visible Infrared Imaging Radiometer Suite (VIIRS), with multi-seasonal daytime optical Landsat time-series, and digital elevation data from the Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER). We employed a non-parametric machine-learning classifier, random forest, to classify the satellite imagery. The classifier was trained with data from three years (2005, 2010, and 2015), and was validated with an independent reference sample from twelve years. Our results demonstrated that combining nighttime light data and multi-seasonal imagery significantly improved the classification (p<0.001), compared to conventional methods based on single-season optical imagery. The results confirmed that the developed classification model enabled mapping of landslides across Taiwan over a long period with annual overall accuracy varying between 96% and 97%, user's and producer's accuracies between 73% and 86%. Spatiotemporal analysis of the landslide inventories from 1998 to 2017 revealed different temporal patterns of landslide activities, showing those areas where landslides were persistent and other areas where landslides tended to reoccur after vegetation regrowth. In sum, we provide a robust method to detect long-term landslide activities based on freely available satellite imagery, which can be applied elsewhere.


title: Discovering causal factors of drought in Ethiopia

abstract: Drought is a costly natural hazard, many aspects of which remain poorly understood. It has many contributory factors, driving its outset, duration, and severity, including land surface, anthropogenic activities, and, most importantly, meteorological anomalies. Prediction plays a crucial role in drought preparedness and risk mitigation. However, this is a challenging task at socio-economically critical lead times (1-2 years), because meteorological anomalies operate at a wide range of temporal and spatial scales. Among them, past studies have shown a correlation between the Sea Surface Temperature (SST) anomaly and the amount of precipitation in various locations in Africa. In its Eastern part, the cooling phase of El Nino-Southern Oscillation (ENSO) and SST anomaly in the Indian ocean are correlated with the lack of rainfall. Given the intrinsic shortcomings of correlation coefficients, we investigate the association among SST modes of variability and the monthly fraction of grid points in Ethiopia, which are in drought conditions in terms of causality. Using the empirical extreme quantiles of precipitation distribution as a proxy for drought, We show that the level of SST second mode of variability in the prior year influences the occurrence of drought in Ethiopia. The causal link between these two variables has a negative coefficient that verifies the conclusion of past studies that rainfall deficiency in the Horn of Africa is associated with ENSO's cooling phase.


title: H\"older continuity of tangent cones in RCD(K,N) spaces and applications to non-branching

abstract: In this paper we prove that a metric measure space $(X,d,m)$ satisfying the finite Riemannian curvature-dimension condition ${\sf RCD}(K,N)$ is non-branching and that tangent cones from the same sequence of rescalings are H\"older continuous along the interior of every geodesic in $X$. More precisely, we show that the geometry of balls of small radius centred in the interior of any geodesic changes in at most a H\"older continuous way along the geodesic in pointed Gromov-Hausdorff distance. This improves a result in the Ricci limit setting by Colding-Naber where the existence of at least one geodesic with such properties between any two points is shown. As in the Ricci limit case, this implies that the regular set of an ${\sf RCD}(K,N)$ space has $m$-a.e. constant dimension, a result already established by Bru\`e-Semola, and is $m$-a.e convex. It also implies that the top dimension regular set is weakly convex and, therefore, connected. In proving the main theorems, we develop in the ${\sf RCD}(K,N)$ setting the expected second order interpolation formula for the distance function along the Regular Lagrangian flow of some vector field using its covariant derivative.


title: Wave-packet scattering at a normal-superconductor interface in two-dimensional materials: a generalized theoretical approach

abstract: A wave-packet time evolution method, based on the split-operator technique, is developed to investigate the scattering of quasi-particles at a normal-superconductor interface of arbitrary profile and shape. As a practical application, we consider a system where low energy electrons can be described as Dirac particles, which is the case for most two-dimensional materials, such as graphene and transition metal dichalcogenides. However the method is easily adapted for other cases such as electrons in few layer black phosphorus, or any Schr\"odinger quasi-particles within the effective mass approximation in semiconductors. We employ the method to revisit Andreev reflection in graphene, where specular and retro reflection cases are observed for electrons scattered by a step-like superconducting region. The effect of opening a zero-gap channel across the superconducting region on the electron and hole scattering is also addressed, as an example of the versatility of the technique proposed here.


title: Toward precise simulations of the coupled ultrafast dynamics of electrons and atomic vibrations in materials

abstract: Ultrafast spectroscopies can access the dynamics of electrons and nuclei at short timescales, shedding light on nonequilibrium phenomena in materials. However, development of accurate calculations to interpret these experiments has lagged behind as widely adopted simulation schemes are limited to sub-picosecond timescales or employ simplified interactions lacking quantitative accuracy. Here we show a precise approach to obtain the time-dependent populations of nonequilibrium electrons and atomic vibrations (phonons) up to tens of picoseconds, with a femtosecond time resolution. Combining first-principles electron-phonon and phonon-phonon interactions with a parallel numerical scheme to time-step the coupled electron and phonon Boltzmann equations, our method provides unprecedented microscopic insight into scattering mechanisms in excited materials. Focusing on graphene as a case study, we demonstrate calculations of ultrafast electron and phonon dynamics, transient optical absorption, structural snapshots and diffuse X-ray scattering. Our first-principles approach paves the way for quantitative atomistic simulations of ultrafast dynamics in materials.


title: Quantum weak invariants: Dynamical evolution of fluctuations and correlations

abstract: Weak invariants are time-dependent observables with conserved expectation values. Their fluctuations, however, do not remain constant in time. On the assumption that time evolution of the state of an open quantum system is given in terms of a completely positive map, the fluctuations monotonically grow even if the map is not unital, in contrast to the fact that monotonic increases of both the von Neumann entropy and R\'enyi entropy require the map to be unital. In this way, the weak invariants describe temporal asymmetry in a manner different from the entropies. A formula is presented for time evolution of the covariance matrix associated with the weak invariants in the case when the system density matrix obeys the Gorini-Kossakowski-Lindblad-Sudarshan equation.


title: Bump attractors and waves in networks of leaky integrate-and-fire neurons

abstract: Bump attractors are wandering localised patterns observed in in vivo experiments of spatially-extended neurobiological networks. They are important for the brain's navigational system and specific memory tasks. A bump attractor is characterised by a core in which neurons fire frequently, while those away from the core do not fire. We uncover a relationship between bump attractors and travelling waves in a classical network of excitable, leaky integrate-and-fire neurons. This relationship bears strong similarities to the one between complex spatiotemporal patterns and waves at the onset of pipe turbulence. Waves in the spiking network are determined by a firing set, that is, the collection of times at which neurons reach a threshold and fire as the wave propagates. We define and study analytical properties of the voltage mapping, an operator transforming a solution's firing set into its spatiotemporal profile. This operator allows us to construct localised travelling waves with an arbitrary number of spikes at the core, and to study their linear stability. A homogeneous "laminar" state exists in the network, and it is linearly stable for all values of the principal control parameter. Sufficiently wide disturbances to the homogeneous state elicit the bump attractor. We show that one can construct waves with a seemingly arbitrary number of spikes at the core; the higher the number of spikes, the slower the wave, and the more its profile resembles a stationary bump. As in the fluid-dynamical analogy, such waves coexist with the homogeneous state, are unstable, and the solution branches to which they belong are disconnected from the laminar state; we provide evidence that the dynamics of the bump attractor displays echoes of the unstable waves, which form its building blocks.


title: Decomposition and Adaptive Sampling for Data-Driven Inverse Linear Optimization

abstract: This work addresses inverse linear optimization where the goal is to infer the unknown cost vector of a linear program. Specifically, we consider the data-driven setting in which the available data are noisy observations of optimal solutions that correspond to different instances of the linear program. We introduce a new formulation of the problem that, compared to other existing methods, allows the recovery of a less restrictive and generally more appropriate admissible set of cost estimates. It can be shown that this inverse optimization problem yields a finite number of solutions, and we develop an exact two-phase algorithm to determine all such solutions. Moreover, we propose an efficient decomposition algorithm to solve large instances of the problem. The algorithm extends naturally to an online learning environment where it can be used to provide quick updates of the cost estimate as new data becomes available over time. For the online setting, we further develop an effective adaptive sampling strategy that guides the selection of the next samples. The efficacy of the proposed methods is demonstrated in computational experiments involving two applications, customer preference learning and cost estimation for production planning. The results show significant reductions in computation and sampling efforts.


title: Heat treatment - microstructure - hardness relationships of new nickel-rich nickel-titanium-hafnium alloys developed for tribological applications

abstract: The effects of various heat treatments on the microstructure and hardness of new Ni56Ti41Hf3 and Ni56Ti36Hf8 (atomic %) alloys were studied to evaluate the suitability of these materials for tribological applications. A solid-solution strengthening effect due to Hf atoms was observed for the solution annealed (SA) Ni56Ti36Hf8 alloy (716 HV), resulting in a comparable hardness to the Ni56Ti41Hf3 alloy containing 54 vol.% of Ni4Ti3 precipitates (707 HV). In the Ni56Ti41Hf3 alloy, the maximum hardness (752 HV), achieved after aging at 300C for 12 h, was attributed to dense, semi-coherent precipitation of the Ni4Ti3 phase. Unlike the lenticular morphology usually observed within binary NiTi alloys, a blocky Ni4Ti3 morphology formed within Ni56Ti36Hf3 due to a smaller lattice mismatch in the direction normal to the habit plane at the precipitate/matrix interface. The maximum hardness for Ni56Ti36Hf8 (769 HV) was obtained after applying an intermediate aging step (300C for 12 h) followed by normal aging (550C for 4 h). This two-step aging treatment induces dense nanoscale precipitation of two interspersed precipitate phases, namely H-phase and a new cubic Ni-rich precipitate phase, resulting in the highest hardness exhibited yet by this family alloys. The composition of cubic Ni-rich precipitates was measured using atom probe tomography to be approximately Ni61.5Ti31Hf7.5, while HAADF-STEM revealed a 54 atom motif cubic structure (a= 8.87 Angstroms), and electron diffraction showed that the structure belongs to the pm-3m (No. 221) space group.


title: Optimal Sepsis Patient Treatment using Human-in-the-loop Artificial Intelligence

abstract: Sepsis is one of the leading causes of death in Intensive Care Units (ICU). The strategy for treating sepsis involves the infusion of intravenous (IV) fluids and administration of antibiotics. Determining the optimal quantity of IV fluids is a challenging problem due to the complexity of a patient's physiology. In this study, we develop a data-driven optimization solution that derives the optimal quantity of IV fluids for individual patients. The proposed method minimizes the probability of severe outcomes by controlling the prescribed quantity of IV fluids and utilizes human-in-the-loop artificial intelligence. We demonstrate the performance of our model on 1122 ICU patients with sepsis diagnosis extracted from the MIMIC-III dataset. The results show that, on average, our model can reduce mortality by 22%. This study has the potential to help physicians synthesize optimal, patient-specific treatment strategies.


title: Tasty Burgers, Soggy Fries: Probing Aspect Robustness in Aspect-Based Sentiment Analysis

abstract: Aspect-based sentiment analysis (ABSA) aims to predict the sentiment towards a specific aspect in the text. However, existing ABSA test sets cannot be used to probe whether a model can distinguish the sentiment of the target aspect from the non-target aspects. To solve this problem, we develop a simple but effective approach to enrich ABSA test sets. Specifically, we generate new examples to disentangle the confounding sentiments of the non-target aspects from the target aspect's sentiment. Based on the SemEval 2014 dataset, we construct the Aspect Robustness Test Set (ARTS) as a comprehensive probe of the aspect robustness of ABSA models. Over 92% data of ARTS show high fluency and desired sentiment on all aspects by human evaluation. Using ARTS, we analyze the robustness of nine ABSA models, and observe, surprisingly, that their accuracy drops by up to 69.73%. We explore several ways to improve aspect robustness, and find that adversarial training can improve models' performance on ARTS by up to 32.85%. Our code and new test set are available at https://github.com/zhijing-jin/ARTS_TestSet


title: Recursive formulation and parallel implementation of multiscale mixed methods

abstract: Multiscale methods for second order elliptic equations based on non-overlapping domain decomposition schemes have great potential to take advantage of multi-core, state-of-the-art parallel computers. These methods typically involve solving local boundary value problems followed by the solution of a global interface problem. Known iterative procedures for the solution of the interface problem have typically slow convergence, increasing the overall cost of the multiscale solver. To overcome this problem we develop a scalable recursive solution method for such interface problem that replaces the global problem by a family of small interface systems associated with adjacent subdomains, in a hierarchy of nested subdomains. Then, we propose a novel parallel algorithm to implement our recursive formulation in multi-core devices using the Multiscale Robin Coupled Method by Guiraldello et al. (2018), that can be seen as a generalization of several multiscale mixed methods. Through several numerical studies we show that the new algorithm is very fast and exhibits excellent strong and weak scalability. We consider very large problems, that can have billions of discretization cells, motivated by the numerical simulation of subsurface flows.


title: Explicit analytic solution for the plane elastostatic problem with a rigid inclusion of arbitrary shape subject to arbitrary far-field loadings

abstract: We provide an analytical solution for the elastic fields in a two-dimensional unbounded isotropic body with a rigid inclusion. Our analysis is based on the boundary integral formulation of the elastostatic problem and geometric function theory. Specifically, we use the coordinate system provided by the exterior conformal mapping of the inclusion to define a density basis functions on the boundary of the inclusion, and we use the Faber polynomials associated with the inclusion for a basis inside the inclusion. The latter, which constitutes the main novelty of our approach, allows us to obtain an explicit series solution for the plane elastostatic problem for an inclusion of arbitrary shape in terms of the given arbitrary far-field loading.


title: Parabolic and elliptic equations with singular or degenerate coefficients: the Dirichlet problem

abstract: We consider the Dirichlet problem for a class of elliptic and parabolic equations in the upper-half space $\mathbb{R}^d_+$, where the coefficients are the product of $x_d^\alpha, \alpha \in (-\infty, 1),$ and a bounded uniformly elliptic matrix of coefficients. Thus, the coefficients are singular or degenerate near the boundary $\{x_d =0\}$ and they may not locally integrable. The novelty of the work is that we find proper weights under which the existence, uniqueness, and regularity of solutions in Sobolev spaces are established. These results appear to be the first of their kind and are new even if the coefficients are constant. They are also readily extended to systems of equations.


title: A Few-Shot Semantic Parser for Wizard-of-Oz Dialogues with the Precise ThingTalk Representation

abstract: Previous attempts to build effective semantic parsers for Wizard-of-Oz (WOZ) conversations suffer from the difficulty in acquiring a high-quality, manually annotated training set. Approaches based only on dialogue synthesis are insufficient, as dialogues generated from state-machine based models are poor approximations of real-life conversations. Furthermore, previously proposed dialogue state representations are ambiguous and lack the precision necessary for building an effective agent. This paper proposes a new dialogue representation and a sample-efficient methodology that can predict precise dialogue states in WOZ conversations. We extended the ThingTalk representation to capture all information an agent needs to respond properly. Our training strategy is sample-efficient: we combine (1) fewshot data sparsely sampling the full dialogue space and (2) synthesized data covering a subset space of dialogues generated by a succinct state-based dialogue model. The completeness of the extended ThingTalk language is demonstrated with a fully operational agent, which is also used in training data synthesis. We demonstrate the effectiveness of our methodology on MultiWOZ 3.0, a reannotation of the MultiWOZ 2.1 dataset in ThingTalk. ThingTalk can represent 98% of the test turns, while the simulator can emulate 85% of the validation set. We train a contextual semantic parser using our strategy, and obtain 79% turn-by-turn exact match accuracy on the reannotated test set.


title: Rolling contact fatigue deformation mechanisms of nickel-rich nickel-titanium-hafnium alloys

abstract: The tribological performance and underlying deformation behavior of Ni55Ti45, Ni54Ti45Hf1 and Ni56Ti36Hf8 alloys were studied using rolling contact fatigue (RCF) testing and transmission electron microscopy (TEM). TEM results of the as-machined RCF rods, prepared using focus ion beam, revealed some damage very close to the surface. TEM results after initial RCF cycling showed that additional damage was mainly confined to deformation bands that propagated several microns into the sample. These bands formed via localized dislocation slip, possibly on multiple slip systems, within the B2 matrix and/or within transformed B19 prime martensite phase under repeated applied contact stress. Further cycling of Ni55Ti45 and Ni54Ti45Hf1 led to shearing and dissolution of the strengthening precipitates within the deformation bands, followed by formation of nanocrystalline grains and finally amorphization of the remaining matrix material within the bands. The Ni56Ti36Hf8 alloy exhibited a notable increase in RCF performance and a smaller damage zone (1.5 microns) compared to the Ni55Ti45 and Ni54Ti45Hf1 alloys (over 6 microns). This was attributed to the low fraction of B2 matrix phase (less than or equal to 13 %) in the Ni56Ti36Hf8 alloy, leading to formation of narrow deformation bands (less than 11 nm) that were incapable of dissolving the much larger precipitates. Instead, the deformation bands were restricted to narrow channels between the dense cubic NiTiHf and H-phase precipitates. In contrast, broad deformation bands accompanied by shearing and eventual dissolution of the Ni4Ti3 precipitates were observed in the Ni55Ti45 and Ni54Ti45Hf1 alloys due to the high area fractions of B2 matrix phase (~49 %).


title: Skeletonization and Reconstruction based on Graph Morphological Transformations

abstract: Multiscale shape skeletonization on pixel adjacency graphs is an advanced intriguing research subject in the field of image processing, computer vision and data mining. The previous works in this area almost focused on the graph vertices. We proposed novel structured based graph morphological transformations based on edges opposite to the current node based transformations and used them for deploying skeletonization and reconstruction of infrared thermal images represented by graphs. The advantage of this method is that many widely used path based approaches become available within this definition of morphological operations. For instance, we use distance maps and image foresting transform (IFT) as two main path based methods are utilized for computing the skeleton of an image. Moreover, In addition, the open question proposed by Maragos et al (2013) about connectivity of graph skeletonization method are discussed and shown to be quite difficult to decide in general case.


title: A Network-Based High-Level Data Classification Algorithm Using Betweenness Centrality

abstract: Data classification is a major machine learning paradigm, which has been widely applied to solve a large number of real-world problems. Traditional data classification techniques consider only physical features (e.g., distance, similarity, or distribution) of the input data. For this reason, those are called \textit{low-level} classification. On the other hand, the human (animal) brain performs both low and high orders of learning and it has a facility in identifying patterns according to the semantic meaning of the input data. Data classification that considers not only physical attributes but also the pattern formation is referred to as \textit{high-level} classification. Several high-level classification techniques have been developed, which make use of complex networks to characterize data patterns and have obtained promising results. In this paper, we propose a pure network-based high-level classification technique that uses the betweenness centrality measure. We test this model in nine different real datasets and compare it with other nine traditional and well-known classification models. The results show us a competent classification performance.


title: Spatial structure of unstable normal modes in a glass-forming liquid

abstract: The phenomenology of glass-forming liquids is often described in terms of their underlying, high-dimensional potential energy surface. In particular, the statistics of stationary points sampled as a function of temperature provides useful insight into the thermodynamics and dynamics of the system. To make contact with the real space physics, however, analysis of the spatial structure of the normal modes is required. In this work, we numerically study the potential energy surface of a glass-forming ternary mixture. Starting from liquid configurations equilibrated over a broad range of temperatures using a swap Monte Carlo method, we locate the nearby stationary points and investigate the spatial architecture and the energetics of the associated unstable modes. Through this spatially-resolved analysis, originally developed to study local minima, we corroborate recent evidence that the nature of the unstable modes changes from delocalized to localized around the mode-coupling temperature. We find that the displacement amplitudes of the delocalized modes have a slowly decaying far field, whereas the localized modes consist of a core with large displacements and a rapidly decaying far field. The fractal dimension of unstable modes around the mobility edge is equal to 1, consistent with the scaling of the participation ratio. Finally, we find that around and below the mode-coupling temperature the unstable modes are localized around structural defects, characterized by a disordered local structure markedly different from the liquid's locally favored structure. These defects are similar to those associated to quasi-localized vibrations in local minima and are good candidates to predict the emergence of localized excitations at low temperature.


title: Hidden energy flows in strongly coupled nonequilibrium systems

abstract: Quantifying the flow of energy within and through fluctuating nanoscale systems poses a significant challenge to understanding microscopic biological machines. A common approach involves coarse-graining, which allows a simplified description of such systems. This has the side effect of inducing so-called hidden contributions (due to sub-resolution dynamics) that complicate the resulting thermodynamics. Here we develop a thermodynamically consistent theory describing the nonequilibrium excess power internal to autonomous systems, and introduce a phenomenological framework to quantify the hidden excess power associated with their operation. We confirm our theoretical predictions in numerical simulations of a minimal model for both a molecular transport motor and a rotary motor.


title: Analysis of Generalizability of Deep Neural Networks Based on the Complexity of Decision Boundary

abstract: For supervised learning models, the analysis of generalization ability (generalizability) is vital because the generalizability expresses how well a model will perform on unseen data. Traditional generalization methods, such as the VC dimension, do not apply to deep neural network (DNN) models. Thus, new theories to explain the generalizability of DNNs are required. In this study, we hypothesize that the DNN with a simpler decision boundary has better generalizability by the law of parsimony (Occam's Razor). We create the decision boundary complexity (DBC) score to define and measure the complexity of decision boundary of DNNs. The idea of the DBC score is to generate data points (called adversarial examples) on or near the decision boundary. Our new approach then measures the complexity of the boundary using the entropy of eigenvalues of these data. The method works equally well for high-dimensional data. We use training data and the trained model to compute the DBC score. And, the ground truth for model's generalizability is its test accuracy. Experiments based on the DBC score have verified our hypothesis. The DBC is shown to provide an effective method to measure the complexity of a decision boundary and gives a quantitative measure of the generalizability of DNNs.


title: Noise-Aware Merging of High Dynamic Range Image Stacks without Camera Calibration

abstract: A near-optimal reconstruction of the radiance of a High Dynamic Range scene from an exposure stack can be obtained by modeling the camera noise distribution. The latent radiance is then estimated using Maximum Likelihood Estimation. But this requires a well-calibrated noise model of the camera, which is difficult to obtain in practice. We show that an unbiased estimation of comparable variance can be obtained with a simpler Poisson noise estimator, which does not require the knowledge of camera-specific noise parameters. We demonstrate this empirically for four different cameras, ranging from a smartphone camera to a full-frame mirrorless camera. Our experimental results are consistent for simulated as well as real images, and across different camera settings.


title: Betti and Hodge numbers of configuration spaces of a punctured elliptic curve from its zeta functions

abstract: Given an elliptic curve $E$ defined over $\mathbb{C}$, let $E^{\times}$ be an open subset of $E$ obtained by removing a point. In this paper, we show that the $i$-th Betti number of the unordered configuration space $\mathrm{Conf}^{n}(E^{\times})$ of $n$ points on $E^{\times}$ appears as a coefficient of an explicit rational function in two variables. We also compute its Hodge numbers as coefficients of another explicit rational function in four variables. Our result is interesting because these rational functions resemble the generating function of the $\mathbb{F}_{q}$-point counts of $\mathrm{Conf}^{n}(E^{\times})$, which can be obtained from the zeta function of $E$ over a finite field $\mathbb{F}_{q}$. We show that the mixed Hodge structure of the $i$-th singular cohomology group $H^{i}(\mathrm{Conf}^{n}(E^{\times}))$ with complex coefficients is pure of weight $w(i)$, an explicit integer we provide in this paper. This purity statement implies our main result about the Betti numbers and the Hodge numbers. Our proof uses Totaro's spectral sequence computation that describes the weight filtration of the mixed Hodge structure on $H^{i}(\mathrm{Conf}^{n}(E^{\times}))$.


title: Classical dynamical density functional theory: from fundamentals to applications

abstract: Classical dynamical density functional theory (DDFT) is one of the cornerstones of modern statistical mechanics. It is an extension of the highly successful method of classical density functional theory (DFT) to nonequilibrium systems. Originally developed for the treatment of simple and complex fluids, DDFT is now applied in fields as diverse as hydrodynamics, materials science, chemistry, biology, and plasma physics. In this review, we give a broad overview over classical DDFT. We explain its theoretical foundations and the ways in which it can be derived. The relations between the different forms of deterministic and stochastic DDFT as well as between DDFT and related theories, such as quantum-mechanical time-dependent DFT, mode coupling theory, and phase field crystal models, are clarified. Moreover, we discuss the wide spectrum of extensions of DDFT, which covers methods with additional order parameters (like extended DDFT), exact approaches (like power functional theory), and systems with more complex dynamics (like active matter). Finally, the large variety of applications, ranging from fluid mechanics and polymer physics to solidification, pattern formation, biophysics, and electrochemistry, is presented.

categories: cond-mat.soft cond-mat.mtrl-sci cond-mat.stat-mech nlin.PS physics.flu-dyn


title: Semidistributivity and Whitman Property in Implication Zroupoids

abstract: In 2012, the second author introduced and studied the variety $\mathcal{I}$ of implication zroupoids that generalize De Morgan algebras and $\lor$-semilattices with $0$. An algebra $\mathbf A = \langle A, \to, 0 \rangle$, where $\to$ is binary and $0$ is a constant, is called an \emph{implication zroupoid} ($\mathcal{I}$-zroupoid, for short) if $\mathbf A$ satisfies: $(x \to y) \to z \approx [(z' \to x) \to (y \to z)']'$, where $x' : = x \to 0$, and $ 0'' \approx 0$. Let $\mathcal{I}$ denote the variety of implication zroupoids and $\mathbf A \in \mathcal{I}$. For $x,y \in \mathbf A$, let $x \land y := (x \to y')'$ and $x \lor y := (x' \land y')'$. In an earlier paper we had proved that if $\mathbf A \in \mathcal{I}$, then the algebra $\mathbf A_{mj} = \langle A, \lor, \land \rangle$ is a bisemigroup. In this paper we generalize the notion of semi-distributivity from lattices to bisemigroups and prove that, for every $\mathbf A \in \mathcal{I}$, the bisemigroup $\mathbf A_{mj}$ is semidistributive. Secondly, we generalize the Whitman Property from lattices to bisemigroups and prove that the subvariety $\mathcal{MEJ}$ of $\mathcal I$, defined by the identity: $x \land y \approx x \lor y$, satisfies the Whitman Property.