title: Distributed Feedback Optimisation for Robotic Coordination

abstract: Feedback optimisation is an emerging technique aiming at steering a system to an optimal steady state for a given objective function. We show that it is possible to employ this control strategy in a distributed manner. Moreover, we prove asymptotic convergence to the set of optimal configurations. To this scope, we show that exponential stability is needed only for the portion of the state that affects the objective function. This is showcased by driving a swarm of agents towards a target location while maintaining a target formation. Finally, we provide a sufficient condition on the topological structure of the specified formation to guarantee convergence of the swarm in formation around the target location.


title: All-optical control of phase singularities using strong light-matter coupling

abstract: Strong light-matter coupling occurs when the coupling strength between a confined electromagnetic mode and a molecular resonance exceeds losses to the environment. The study of strong coupling has been motivated by applications such as lasing and the modification of chemical processes. Here we show that strong coupling can be used to create phase singularities. Many nanophotonic structures have been designed to generate phase singularities for use in sensing and optoelectronics. We utilise the concept of cavity-free strong coupling, where electromagnetic modes sustained by a material are strong enough to strongly couple to the material's own molecular resonance, to create phase singularities in a simple thin film of organic molecules. We show that the use of photochromic molecules allows for all-optical control of phase singularities. Our results suggest a new application for strong light-matter coupling and a new, simplified, more versatile pathway to singular phase optics.


title: The structure of digraphs with excess one

abstract: A digraph $G$ is \emph{$k$-geodetic} if for any (not necessarily distinct) vertices $u,v$ there is at most one directed walk from $u$ to $v$ with length not exceeding $k$. The order of a $k$-geodetic digraph with minimum out-degree $d$ is bounded below by the directed Moore bound $M(d,k) = 1+d+d^2+\dots +d^k$. The Moore bound can be met only in the trivial cases $d = 1$ and $k = 1$, so it is of interest to look for $k$-geodetic digraphs with out-degree $d$ and smallest possible order $M(d,k)+\epsilon $, where $\epsilon $ is the \emph{excess} of the digraph. Miller, Miret and Sillasen recently ruled out the existence of digraphs with excess one for $k = 3,4$ and $d \geq 2$ and for $k = 2$ and $d \geq 8$. We conjecture that there are no digraphs with excess one for $d,k \geq 2$ and in this paper we investigate the structure of minimal counterexamples to this conjecture. We severely constrain the possible structures of the outlier function and prove the non-existence of certain digraphs with degree three and excess one, as well closing the open cases $k = 2$ and $d = 3,4,5,6,7$ left by the analysis of Miller et al. We further show that there are no involutary digraphs with excess one, i.e. the outlier function of any such digraph must contain a cycle of length $\geq 3$.


title: A portable muon telescope for multidisciplinary applications

abstract: Muon tomography or muography is an emerging imaging technique that uses cosmogenic muons as the radiation source. Due to its diverse range of applications and the use of natural radiation, muography is being applied across many fields such as geology, archaeology, civil engineering, nuclear reactor monitoring, nuclear waste characterization, underground surveys, etc. Muons can be detected using various detector technologies, among which, resistive plate chambers (RPC) are a very cost effective choice. RPCs are planar detectors which use ionization in a thin gas gap to detect cosmic muons, already used since years in major particle accelerator experiments. We have developed a muon telescope (or muoscope) composed of small scale RPCs. The design goal for our muoscope is to be portable and autonomous, in order to take data in places that are not easily accessible. The whole setup is light and compact, such to be easily packed in a car trunk. Individual RPCs are hosted in gas-tight aluminium cases. There is no need for gas bottles, once the chambers are filled. The muoscope can be controlled from a reasonable distance using wireless connection. In this paper, we summarize the guiding principles of our project and present some recent developments and future prospects, including a long-term stability study of the resistivity of the semiconductive coating obtained with serigraphy.


title: Might I Get Pwned: A Second Generation Compromised Credential Checking Service

abstract: Credential stuffing attacks use stolen passwords to log into victim accounts. To defend against these attacks, recently deployed compromised credential checking (C3) services provide APIs that help users and companies check whether a username, password pair is exposed. These services however only check if the exact password is leaked, and therefore do not mitigate credential tweaking attacks - attempts to compromise a user account with variants of a user's leaked passwords. Recent work has shown credential tweaking attacks can compromise accounts quite effectively even when the credential stuffing countermeasures are in place. We initiate work on C3 services that protect users from credential tweaking attacks. The core underlying challenge is how to identify passwords that are similar to their leaked passwords while preserving honest clients' privacy and also preventing malicious clients from extracting breach data from the service. We formalize the problem and explore ways to measure password similarity that balance efficacy, performance, and security. Based on this study, we design "Might I Get Pwned" (MIGP), a new kind of breach alerting service. Our simulations show that MIGP reduces the efficacy of state-of-the-art 1000-guess credential tweaking attacks by 94%. MIGP preserves user privacy and limits potential exposure of sensitive breach entries. We show that the protocol is fast, with response time close to existing C3 services. We worked with Cloudflare to deploy MIGP in practice.


title: Hydrodynamic limit of Exclusion Processes with slow boundaries on hypercubes

abstract: We study the hydrodynamic limit of SSEP with slow boundaries on hypercubes in dimension at least two. The hydrodynamic limit equation is shown to be a heat equation with three different types of boundary conditions according to the slowness of the boundary dynamics.The proof is based on Yau's relative entropy method.


title: Variational Inference for Continuous-Time Switching Dynamical Systems

abstract: Switching dynamical systems provide a powerful, interpretable modeling framework for inference in time-series data in, e.g., the natural sciences or engineering applications. Since many areas, such as biology or discrete-event systems, are naturally described in continuous time, we present a model based on an Markov jump process modulating a subordinated diffusion process. We provide the exact evolution equations for the prior and posterior marginal densities, the direct solutions of which are however computationally intractable. Therefore, we develop a new continuous-time variational inference algorithm, combining a Gaussian process approximation on the diffusion level with posterior inference for Markov jump processes. By minimizing the path-wise Kullback-Leibler divergence we obtain (i) Bayesian latent state estimates for arbitrary points on the real axis and (ii) point estimates of unknown system parameters, utilizing variational expectation maximization. We extensively evaluate our algorithm under the model assumption and for real-world examples.


title: Automatic Discovery and Description of Human Planning Strategies

abstract: Scientific discovery concerns finding patterns in data and creating insightful hypotheses that explain these patterns. Traditionally, this process required human ingenuity, but with the galloping advances in artificial intelligence (AI) it becomes feasible to automate some parts of scientific discovery. In this work we leverage AI for strategy discovery for understanding human planning. In the state-of-the-art methods data about the process of human planning is often used to group similar behaviors together and formulate verbal descriptions of the strategies which might underlie those groups. Here, we automate these two steps. Our method utilizes a new algorithm, called Human-Interpret, that performs imitation learning to describe sequences of planning operations in terms of a procedural formula and then translates that formula to natural language. We test our method on a benchmark data set that researchers have previously scrutinized manually. We find that the descriptions of human planning strategies obtained automatically are about as understandable as human-generated descriptions. They also cover a substantial proportion of of relevant types of human planning strategies that had been discovered manually. Our method saves scientists' time and effort as all the reasoning about human planning is done automatically. This might make it feasible to more rapidly scale up the search for yet undiscovered cognitive strategies to many new decision environments, populations, tasks, and domains. Given these results, we believe that the presented work may accelerate scientific discovery in psychology, and due to its generality, extend to problems from other fields.


title: Bounding the separable rank via polynomial optimization

abstract: We investigate questions related to the set $\mathcal{SEP}_d$ consisting of the linear maps $\rho$ acting on $\mathbb{C}^d\otimes \mathbb{C}^d$ that can be written as a convex combination of rank one matrices of the form $xx^*\otimes yy^*$. Such maps are known in quantum information theory as the separable bipartite states, while nonseparable states are called entangled. In particular we introduce bounds for the separable rank $\mathrm{rank_{sep}}(\rho)$, defined as the smallest number of rank one states $xx^*\otimes yy^*$ entering the decomposition of a separable state $\rho$. Our approach relies on the moment method and yields a hierarchy of semidefinite-based lower bounds, that converges to a parameter $\tau_{\mathrm{sep}}(\rho)$, a natural convexification of the combinatorial parameter $\mathrm{rank_{sep}}(\rho)$. A distinguishing feature is exploiting the positivity constraint $\rho-xx^*\otimes yy^* \succeq 0$ to impose positivity of a polynomial matrix localizing map, the dual notion of the notion of sum-of-squares polynomial matrices. Our approach extends naturally to the multipartite setting and to the real separable rank, and it permits strengthening some known bounds for the completely positive rank. In addition, we indicate how the moment approach also applies to define hierarchies of semidefinite relaxations for the set $\mathcal{SEP}_d$ and permits to give new proofs, using only tools from moment theory, for convergence results on the DPS hierarchy from (A.C. Doherty, P.A. Parrilo and F.M. Spedalieri. Distinguishing separable and entangled states. Phys. Rev. Lett. 88(18):187904, 2002).


title: On bounded two-dimensional globally dissipative Euler flows

abstract: We examine the two-dimensional Euler equations including the local energy (in)equality as a differential inclusion and show that the associated relaxation essentially reduces to the known relaxation for the Euler equations considered without local energy (im)balance. Concerning bounded solutions we provide a sufficient criterion for a globally dissipative subsolution to induce infinitely many globally dissipative solutions having the same initial data, pressure and dissipation measure as the subsolution. The criterion can easily be verified in the case of a flat vortex sheet giving rise to the Kelvin-Helmholtz instability. As another application we show that there exists initial data, for which associated globally dissipative solutions realize every dissipation measure from an open set in $\mathcal{C}^0(\mathbb{T}^2\times[0,T])$. In fact the set of such initial data is dense in the space of solenoidal $L^2(\mathbb{T}^2;\mathbb{R}^2)$ vector fields.


title: Stationary generalizations for the Bronnikov-Ellis wormhole and for the vacuum ring wormhole

abstract: We analyze possibilities to obtain a globally regular stationary generalization for the ultrastatic wormhole with a repulsive scalar field found by Bronnikov and by Ellis in 1973. The extreme simplicity of this static solution suggests that its stationary version could be obtainable analytically and should be globally regular, but no such generalization has been found. We analyze the problem and find that the difficulty originates in the vacuum theory, since the scalar field can be eliminated within the Eris-Gurses procedure. The problem then reduces to constructing the spinning generalization for the vacuum wormhole sourced by a thin ring of negative tension. Solving the vacuum Ernst equations determines the $g_{00}$, $g_{0\varphi}$ metric components and hence the AMD mass ${\rm M}$ and angular momentum ${\rm J}$, all of these being specified by the ring source. The scalar field can be included into consideration afterwords, but this only affects $g_{rr}$ and $g_{\vartheta\vartheta}$ without changing the rest. Within this approach, we analyze a number of exact stationary generalizations for the wormhole, but none of them are satisfactory. However, the perturbative expansion around the static vacuum background contains only bounded functions and presumably converges to an exact solution. Including the scalar field screens the singularity at the ring source and renders the geometry regular. This solution describes a globally regular spinning wormhole with two asymptotically flat regions. Even though the source itself is screened and not visible, the memory of it remains in $g_{00}$ and $g_{0\varphi}$ and accounts for the ${\rm M}\propto {\rm J}^2$ relation typical for a rotating extended source. Describing stationary spacetimes with an extended source is a complicated problem, which presumably explains the difficulty in finding the solution.


title: Ruler Wrapping

abstract: In 1985 Hopcroft, Joseph and Whitesides showed it is NP-complete to decide whether a carpenter's ruler with segments of given positive lengths can be folded into a line of at most a given length, such that the folded hinges alternate between 180 degrees clockwise and 180 degrees counter-clockwise. At the open-problem session of 33rd Canadian Conference on Computational Geometry (CCCG '21), O'Rourke proposed a natural variation of this problem called {\em ruler wrapping}, in which all folded hinges must be folded the same way. In this paper we show O'Rourke's variation has an linear-time solution. We also show how, given a sequence of positive numbers, in linear time we can partition it into the maximum number of substrings whose totals are non-decreasing.


title: Density conditions with stabilizers for lattice orbits of Bergman kernels on bounded symmetric domains

abstract: Let $\pi_{\alpha}$ be a holomorphic discrete series representation of a connected semi-simple Lie group $G$ with finite center, acting on a weighted Bergman space $A^2_{\alpha} (\Omega)$ on a bounded symmetric domain $\Omega$, of formal dimension $d_{\pi_{\alpha}} > 0$. It is shown that if the Bergman kernel $k^{(\alpha)}_z$ is a cyclic vector for the restriction $\pi_{\alpha} |_{\Gamma}$ to a lattice $\Gamma \leq G$ (resp. $(\pi_{\alpha} (\gamma) k^{(\alpha)}_z)_{\gamma \in \Gamma}$ is a frame for $A^2_{\alpha}(\Omega)$), then $\mathrm{vol}(G/\Gamma) d_{\pi_{\alpha}} \leq |\Gamma_z|^{-1}$. The estimate $\mathrm{vol}(G/\Gamma) d_{\pi_{\alpha}} \geq |\Gamma_z|^{-1}$ holds for $k^{(\alpha)}_z$ being a $p_z$-separating vector (resp. $(\pi_{\alpha} (\gamma) k^{(\alpha)}_z)_{\gamma \in \Gamma / \Gamma_z}$ being a Riesz sequence in $A^2_{\alpha} (\Omega)$). These estimates improve on general density theorems for restricted discrete series through the dependence on the stabilizers, while recovering in part sharp results for $G = \mathrm{PSU}(1, 1)$.


title: $2$-distance list $(\Delta+2)$-coloring of planar graphs with girth at least 10

abstract: Given a graph $G$ and a list assignment $L(v)$ for each vertex of $v$ of $G$. A proper $L$-list-coloring of $G$ is a function that maps every vertex to a color in $L(v)$ such that no pair of adjacent vertices have the same color. We say that a graph is list $k$-colorable when every vertex $v$ has a list of colors of size at least $k$. A $2$-distance coloring is a coloring where vertices at distance at most 2 cannot share the same color. We prove the existence of a $2$-distance list ($\Delta+2$)-coloring for planar graphs with girth at least $10$ and maximum degree $\Delta\geq 4$.


title: Scaling Behavior of the Hirsch Index for Failure Avalanches, Percolation Clusters and Paper Citations

abstract: A popular measure for citation inequalities of individual scientists has been the Hirsch index ($h$). If for any scientist the number $n_c$ of citations is plotted against the serial number $n_p$ of the paper having those many citations (when the papers are ordered from highest cited to lowest) then $h$ corresponds to the nearest lower integer value of $n_p$ below the fixed point of the non-linear citation function (or given by $n_c = h = n_p$ if both $n_p$ and $n_c$ are dense set of integers near the $h$ value). The same index can be estimated (from $h=s=n_{s}$) for the avalanche or cluster of size ($s$) distributions ($n_s$) in elastic fiber bundle or percolation models. Another such inequality index, called the Kolkata index ($k$) says that $(1-k)$ fraction of papers attract $k$ fraction of citations ($k=0.80$ corresponds to the 80-20 law of Pareto). We find, for stress ($\sigma$), lattice occupation probability ($p$) or Kolkata index ($k$) near the bundle failure threshold ($\sigma_c$) or percolation threshold ($p_c$) or critical value of Kolkata index $k_c$, good fit to Widom-Stauffer like scaling $h/[\sqrt{N}/log N]$ = $f(\sqrt{N}[\sigma_c -\sigma]^\alpha)$, $h/[\sqrt{N}/log N]=f(\sqrt{N}|p_c -p|^\alpha)$ or $h/[\sqrt{N_c}/log N_c]=f(\sqrt{N_c}|k_c -k|^\alpha)$ respectively, with asymptotically defined scaling function $f$, for systems of size $N$ (total number of fibers or lattice sites) or $N_c$ (total number of citations), and $\alpha$ denoting the appropriate scaling exponent. We also show that if the number ($N_m$) of members of parliaments or national assemblies of different countries (with population $N$) is identified as their respective $h-$index, then the data fit the scaling relation $N_m \sim \sqrt N /log N$, resolving a major recent controversy.


title: Towards a theory of out-of-distribution learning

abstract: Learning is a process wherein a learning agent enhances its performance through exposure of experience or data. Throughout this journey, the agent may encounter diverse learning environments. For example, data may be presented to the leaner all at once, in multiple batches, or sequentially. Furthermore, the distribution of each data sample could be either identical and independent (iid) or non-iid. Additionally, there may exist computational and space constraints for the deployment of the learning algorithms. The complexity of a learning task can vary significantly, depending on the learning setup and the constraints imposed upon it. However, it is worth noting that the current literature lacks formal definitions for many of the in-distribution and out-of-distribution learning paradigms. Establishing proper and universally agreed-upon definitions for these learning setups is essential for thoroughly exploring the evolution of ideas across different learning scenarios and deriving generalized mathematical bounds for these learners. In this paper, we aim to address this issue by proposing a chronological approach to defining different learning tasks using the provably approximately correct (PAC) learning framework. We will start with in-distribution learning and progress to recently proposed lifelong or continual learning. We employ consistent terminology and notation to demonstrate how each of these learning frameworks represents a specific instance of a broader, more generalized concept of learnability. Our hope is that this work will inspire a universally agreed-upon approach to quantifying different types of learning, fostering greater understanding and progress in the field.


title: Untangling Braids with Multi-agent Q-Learning

abstract: We use reinforcement learning to tackle the problem of untangling braids. We experiment with braids with 2 and 3 strands. Two competing players learn to tangle and untangle a braid. We interface the braid untangling problem with the OpenAI Gym environment, a widely used way of connecting agents to reinforcement learning problems. The results provide evidence that the more we train the system, the better the untangling player gets at untangling braids. At the same time, our tangling player produces good examples of tangled braids.


title: Nonreciprocal transport in a Rashba ferromagnet, delafossite PdCoO$_2$

abstract: Rashba interfaces yield efficient spin-charge interconversion and give rise to nonreciprocal transport phenomena. Here, we report magnetotransport experiments in few-nanometer-thick films of PdCoO$_2$, a delafossite oxide known to display a large Rashba splitting and surface ferromagnetism. By analyzing the angle dependence of the first- and second-harmonic longitudinal and transverse resistivities, we identify a Rashba-driven unidirectional magnetoresistance that competes with the anomalous Nernst effect below the Curie point. We estimate a Rashba coefficient of 0.75 {\pm} 0.3 eV {\AA} and argue that our results qualify delafossites as a new family of oxides for nano-spintronics and spin-orbitronics, beyond perovskite materials.


title: Random sections of $\ell_p$-ellipsoids, optimal recovery and Gelfand numbers of diagonal operators

abstract: We study the circumradius of a random section of an $\ell_p$-ellipsoid, $0<p\le \infty$, and compare it with the minimal circumradius over all sections with subspaces of the same codimension. Our main result is an upper bound for random sections, which we prove using techniques from asymptotic geometric analysis if $1\leq p \leq \infty$ and compressed sensing if $0<p \leq 1$. This can be interpreted as a bound on the quality of random (Gaussian) information for the recovery of vectors from an $\ell_p$-ellipsoid for which the radius of optimal information is given by the Gelfand numbers of a diagonal operator. In the case where the semiaxes decay polynomially and $1\le p\le \infty$, we conjecture that, as the amount of information increases, the radius of random information either decays like the radius of optimal information or is bounded from below by a constant, depending on whether the exponent of decay is larger than the critical value $1-\frac{1}{p}$ or not. If $1\leq p\leq 2$, we prove this conjecture by providing a matching lower bound. This extends the recent work of Hinrichs et al. [Random sections of ellipsoids and the power of random information, Trans. Amer. Math. Soc., 2021+] for the case $p=2$.


title: Optical spin transport theory of spin-1/2 topological Fermi superfluids

abstract: We theoretically investigate optical (frequency-dependent) bulk spin transport properties in a spin-1/2 topological Fermi superfluid. We specifically consider a one-dimensional system with an interspin {\it p}-wave interaction, which can be realized in ultracold atom experiments. Developing the BCS-Leggett theory to describe the Bardeen-Cooper-Schrieffer (BCS) to Bose-Einstein condensate (BEC) evolution and the $\mathbb{Z}_2$ topological phase transition in this system, we show how the spin transport reflects these many-body aspects. We find that the optical spin conductivity, which is a small AC response of a spin current, shows the spin gapped spectrum in the wide parameter region and the gap closes at $\mathbb{Z}_2$ topological phase transition point. Moreover, the validity of the low-energy effective model of the Majorana zero mode is discussed along the BCS-BEC evolution in connection with the scale invariance at {\it p}-wave unitarity.


title: Vibrational Kerr solitons in an optomechanical microresonator

abstract: Soliton microcombs based on Kerr nonlinearity in microresonators have been a prominent miniaturized coherent light source. Here, for the first time, we demonstrate the existence of Kerr solitons in an optomechanical microresonator, for which a nonlinear model is built by incorporating a single mechanical mode and multiple optical modes. Interestingly, an exotic vibrational Kerr soliton state is found, which is modulated by a self-sustained mechanical oscillation. Besides, the soliton provides extra mechanical gain through the optical spring effect, and results in phonon lasing with a red-detuned pump. Various nonlinear dynamics is also observed, including limit cycle, higher periodicity, and transient chaos. This work provides a guidance for not only exploring many-body nonlinear interactions, but also promoting precision measurements by featuring superiority of both frequency combs and optomechanics.


title: A Deformation Quantization for Non-Flat Spacetimes and Applications to QFT

abstract: We provide a deformation quantization, in the sense of Rieffel, for \textit{all} globally hyperbolic spacetimes with a Poisson structure. The Poisson structures have to satisfy Fedosov type requirements in order for the deformed product to be associative. We apply the novel deformation to quantum field theories and their respective states and we prove that the deformed state (i.e.\ a state in non-commutative spacetime) has a singularity structure resembling Minkowski, i.e.\ is \textit{Hadamard}, if the undeformed state is Hadamard. This proves that the Hadamard condition, and hence the quantum field theoretical implementation of the equivalence principle is a general concept that holds in spacetimes with quantum features (i.e. a non-commutative spacetime).


title: Cross-domain Semi-Supervised Audio Event Classification Using Contrastive Regularization

abstract: In this study, we proposed a novel semi-supervised training method that uses unlabeled data with a class distribution that is completely different from the target data or data without a target label. To this end, we introduce a contrastive regularization that is designed to be target task-oriented and trained simultaneously. In addition, we propose an audio mixing based simple augmentation strategy that performed in batch samples. Experimental results validate that the proposed method successfully contributed to the performance improvement, and particularly showed that it has advantages in stable training and generalization.


title: PAC-Bayes Information Bottleneck

abstract: Understanding the source of the superior generalization ability of NNs remains one of the most important problems in ML research. There have been a series of theoretical works trying to derive non-vacuous bounds for NNs. Recently, the compression of information stored in weights (IIW) is proved to play a key role in NNs generalization based on the PAC-Bayes theorem. However, no solution of IIW has ever been provided, which builds a barrier for further investigation of the IIW's property and its potential in practical deep learning. In this paper, we propose an algorithm for the efficient approximation of IIW. Then, we build an IIW-based information bottleneck on the trade-off between accuracy and information complexity of NNs, namely PIB. From PIB, we can empirically identify the fitting to compressing phase transition during NNs' training and the concrete connection between the IIW compression and the generalization. Besides, we verify that IIW is able to explain NNs in broad cases, e.g., varying batch sizes, over-parameterization, and noisy labels. Moreover, we propose an MCMC-based algorithm to sample from the optimal weight posterior characterized by PIB, which fulfills the potential of IIW in enhancing NNs in practice.


title: Random coordinate descent algorithm for open multi-agent systems with complete topology and homogeneous agents

abstract: We study the convergence in expectation of the Random Coordinate Descent algorithm (RCD) for solving optimal resource allocations problems in open multi-agent systems, i.e., multi-agent systems that are subject to arrivals and departures of agents. Assuming all local functions are strongly-convex and smooth, and their minimizers lie in a given ball, we analyse the evolution of the distance to the minimizer in expectation when the system is occasionally subject to replacements in addition to the usual iterations of the RCD algorithm. We focus on complete graphs where all agents interact with each other with the same probability, and provide conditions to guarantee convergence in open system. Finally, a discussion around the tightness of our results is provided.


title: Elliptica: a new pseudo-spectral code for the construction of initial data

abstract: Numerical studies of the dynamics of gravitational systems, e.g., black hole-neutron star systems, require physical and constraint-satisfying initial data. In this article, we present the newly developed pseudo-spectral code Elliptica, an infrastructure for construction of initial data for various binary and single gravitational systems of all kinds. The elliptic equations under consideration are solved on a single spatial hypersurface of the spacetime manifold. Using coordinate maps, the hypersurface is covered by patches whose boundaries can adapt to the surface of the compact objects. To solve elliptic equations with arbitrary boundary condition, Elliptica deploys a Schur complement domain decomposition method with a direct solver. In this version, we use cubed sphere coordinate maps and the fields are expanded using Chebyshev polynomials of the first kind. Here, we explain the building blocks of Elliptica and the initial data construction algorithm for a black hole-neutron star binary system. We perform convergence tests and evolve the data to validate our results. Within our framework, the neutron star can reach spin values close to breakup with arbitrary direction, while the black hole can have arbitrary spin with dimensionless spin magnitude $\sim 0.8$.


title: A search for the variation of speed of light using galaxy cluster gas mass fraction measurements

abstract: In this paper, we implement a new method to test the invariance of the speed of light ($c$) as a function of redshift, by combining the measurements of galaxy cluster gas mass fraction, $H(z)$ from cosmic chronometers, and Type-Ia supernovae (SNe Ia). In our analyses, we consider both a constant depletion factor (which corresponds to the ratio by which the cluster gas mass fraction is depleted with respect to the universal baryonic mean) and one varying with redshift. We also consider the influence of different $H_0$ estimates on our results. We look for a variation of $c$, given by $c(z)=c_0(1+c_1z)$. We find a degeneracy between our final results on $c$ variation and the assumptions on the gas mass fraction depletion factor. Most of our analyses indicate negligible variation of the speed of light.


title: $C_p$-Theory for Model Theorists

abstract: We present applications of $C_p$-theory, the branch of general topology concerned with spaces of real-valued continuous functions, to model theory, mostly in the context of continuous logics. We include $C_p$-theoretic results and proofs in a self-contained way for model theorists who are not familiar with the techniques of this field. We further generalize some results of Casazza and Iovino, and of the authors, involving the definability of Banach spaces including isomorphic copies of $c_0$ or $\ell^p$, after a problem posed by Odell and Gowers.


title: ASAS J071404+7004.3 -- a close, bright nova-like cataclysmic variable with gusty winds

abstract: Despite being bright ($V=12$) and nearby ($d=212$ pc) ASAS J071404+7004.3 has only recently been identified as a nova-like cataclysmic variable. We present time-resolved optical spectroscopy obtained at the Isaac Newton Telescope together with $\textit{Swift}$ X-ray and ultraviolet observations. We combined these with $\textit{TESS}$ photometry and find a period of 3.28h and a mass transfer rate of $4-9 \times 10^{-9} M_{sun}/yr$. Historical photometry shows at least one low state establishing the system as a VY Scl star. Our high-cadence spectroscopy also revealed rapidly changing winds emanating from the accretion disc. We have modelled these using the Monte Carlo PYTHON code and shown that all the emission lines could emanate from the wind - which would explain the lack of double-peaked lines in such systems. In passing,we discuss the effect of variability on the position of cataclysmic variables in the $\textit{Gaia}$ Hertzsprung-Russell diagram.


title: The effect of buoyancy driven convection on the growth and dissolution of bubbles on electrodes

abstract: Enhancing the efficiency of water electrolysis, which can be severely impacted by the nucleation and growth of bubbles, is key in the energy transition. In this combined experimental and numerical study, in-situ bubble evolution and dissolution processes are imaged and compared to numerical simulations employing the immersed boundary method. We find that it is crucial to include solutal driven natural convection in order to represent the experimentally observed bubble behaviour even though such effects have commonly been neglected in modelling efforts so far. We reveal how the convective patterns depend on current densities and bubble spacings, leading to distinctively different bubble growth and shrinkage dynamics. Bubbles are seen to promote the convective instability if their spacing is large ($\geq4$mm for the present conditions), whereas the onset of convection is delayed if the inter-bubble distance is smaller. Our approach and our results can help devise efficient mass transfer solutions for gas evolving electrodes.


title: On Assessing the Usefulness of Proxy Domains for Developing and Evaluating Embodied Agents

abstract: In many situations it is either impossible or impractical to develop and evaluate agents entirely on the target domain on which they will be deployed. This is particularly true in robotics, where doing experiments on hardware is much more arduous than in simulation. This has become arguably more so in the case of learning-based agents. To this end, considerable recent effort has been devoted to developing increasingly realistic and higher fidelity simulators. However, we lack any principled way to evaluate how good a "proxy domain" is, specifically in terms of how useful it is in helping us achieve our end objective of building an agent that performs well in the target domain. In this work, we investigate methods to address this need. We begin by clearly separating two uses of proxy domains that are often conflated: 1) their ability to be a faithful predictor of agent performance and 2) their ability to be a useful tool for learning. In this paper, we attempt to clarify the role of proxy domains and establish new proxy usefulness (PU) metrics to compare the usefulness of different proxy domains. We propose the relative predictive PU to assess the predictive ability of a proxy domain and the learning PU to quantify the usefulness of a proxy as a tool to generate learning data. Furthermore, we argue that the value of a proxy is conditioned on the task that it is being used to help solve. We demonstrate how these new metrics can be used to optimize parameters of the proxy domain for which obtaining ground truth via system identification is not trivial.


title: Shuffle algebras for quivers and R-matrices

abstract: We define slope subalgebras in the shuffle algebra associated to a (doubled) quiver, thus yielding a factorization of the universal R-matrix of the double of the shuffle algebra in question. We conjecture that this factorization matches the one defined by [1, 18, 32, 33, 34] using Nakajima quiver varieties.


abstract: We propose Generative Probabilistic Image Colorization, a diffusion-based generative process that trains a sequence of probabilistic models to reverse each step of noise corruption. Given a line-drawing image as input, our method suggests multiple candidate colorized images. Therefore, our method accounts for the ill-posed nature of the colorization problem. We conducted comprehensive experiments investigating the colorization of line-drawing images, report the influence of a score-based MCMC approach that corrects the marginal distribution of estimated samples, and further compare different combinations of models and the similarity of their generated images. Despite using only a relatively small training dataset, we experimentally develop a method to generate multiple diverse colorization candidates which avoids mode collapse and does not require any additional constraints, losses, or re-training with alternative training conditions. Our proposed approach performed well not only on color-conditional image generation tasks using biased initial values, but also on some practical image completion and inpainting tasks.


title: Functional a posteriori error estimates for parabolic obstacle problems

abstract: The paper is concerned with functional type a posteriori estimates for the initial boundary value problem for a parabolic partial differential equation with an obstacle. We deduce a guaranteed and computable bound of the distance between the exact minimizer and any function from the admissible (energy) class of functions. Applications to the analysis of modeling errors caused by data implification are discussed. An important case of time incremental approximations is specially studied. Numerical examples presented in the last section show how the estimates work in practice.


title: Improving DRAM Performance, Security, and Reliability by Understanding and Exploiting DRAM Timing Parameter Margins

abstract: This dissertation rigorously characterizes many modern commodity DRAM devices and shows that by exploiting DRAM access timing margins within manufacturer-recommended DRAM timing specifications, we can significantly improve system performance, reduce power consumption, and improve device reliability and security. First, we characterize DRAM timing parameter margins and find that certain regions of DRAM can be accessed faster than other regions due to DRAM cell process manufacturing variation. We exploit this by enabling variable access times depending on the DRAM cells being accessed, which not only improves overall system performance, but also decreases power consumption. Second, we find that we can uniquely identify DRAM devices by the locations of failures that result when we access DRAM with timing parameters reduced below specification values. Because we induce these failures with DRAM accesses, we can generate these unique identifiers significantly more quickly than prior work. Third, we propose a random number generator that is based on our observation that timing failures in certain DRAM cells are randomly induced and can thus be repeatedly polled to very quickly generate true random values. Finally, we characterize the RowHammer security vulnerability on a wide range of modern DRAM chips while violating the DRAM refresh requirement in order to directly characterize the underlying DRAM technology without the interference of refresh commands. We demonstrate with our characterization of real chips, that existing RowHammer mitigation mechanisms either are not scalable or suffer from prohibitively large performance overheads in projected future devices and it is critical to research more effective solutions to RowHammer. Overall, our studies build a new understanding of modern DRAM devices to improve computing system performance, reliability and security all at the same time.


title: Finite volume methods for the computation of statistical solutions of the incompressible Euler equations

abstract: We present an efficient numerical scheme based on Monte Carlo integration to approximate statistical solutions of the incompressible Euler equations. The scheme is based on finite volume methods, which provide a more flexible framework than previously existing spectral methods for the computation of statistical solutions for incompressible flows. This finite volume scheme is rigorously proven to, under experimentally verifiable assumptions, converge in an appropriate topology and with increasing resolution to a statistical solution. The convergence obtained is stronger than that of measure-valued solutions, as it implies convergence of multi-point correlation marginals. We present results of numerical experiments which support the claim that the aforementioned assumptions are very natural, and appear to hold in practice.


title: Lipschitz Analysis of Generalized Phase Retrievable Matrix Frames

abstract: The classical phase retrieval problem arises in contexts ranging from speech recognition to x-ray crystallography and quantum state tomography. The generalization to matrix frames is natural in the sense that it corresponds to quantum tomography of impure states. We provide computable global stability bounds for the quasi-linear analysis map $\beta$ and a path forward for understanding related problems in terms of the differential geometry of key spaces. In particular, we manifest a Whitney stratification of the positive semidefinite matrices of low rank which allows us to ``stratify'' the computation of the global stability bound. We show that for the impure state case no such global stability bounds can be obtained for the non-linear analysis map $\alpha$ with respect to certain natural distance metrics. Finally, our computation of the global lower Lipschitz constant for the $\beta$ analysis map provides novel conditions for a frame to be generalized phase retrievable.


title: Online Robust Reinforcement Learning with Model Uncertainty

abstract: Robust reinforcement learning (RL) is to find a policy that optimizes the worst-case performance over an uncertainty set of MDPs. In this paper, we focus on model-free robust RL, where the uncertainty set is defined to be centering at a misspecified MDP that generates a single sample trajectory sequentially and is assumed to be unknown. We develop a sample-based approach to estimate the unknown uncertainty set and design a robust Q-learning algorithm (tabular case) and robust TDC algorithm (function approximation setting), which can be implemented in an online and incremental fashion. For the robust Q-learning algorithm, we prove that it converges to the optimal robust Q function, and for the robust TDC algorithm, we prove that it converges asymptotically to some stationary points. Unlike the results in [Roy et al., 2017], our algorithms do not need any additional conditions on the discount factor to guarantee the convergence. We further characterize the finite-time error bounds of the two algorithms and show that both the robust Q-learning and robust TDC algorithms converge as fast as their vanilla counterparts(within a constant factor). Our numerical experiments further demonstrate the robustness of our algorithms. Our approach can be readily extended to robustify many other algorithms, e.g., TD, SARSA, and other GTD algorithms.


title: Search for a bound di-neutron by comparing $^3$He(e,e'p)d and $^3$H(e,e'p)X measurements

abstract: We report on a search for a bound di-neutron by comparing electron-induced proton-knockout $(e,e'p)$ measurements from Helium-3 ($^3$He) and Tritium ($^3$H). The measurements were performed at Jefferson Lab Hall A with a 4.326 GeV electron beam, and kinematics of large momentum transfer $Q^2 \approx 1.9$ (GeV/$c$)$^2$ and $x_B>1$, to minimize contributions from non quasi-elastic (QE) reaction mechanisms. Analyzing the measured $^3$He missing mass ($M_{miss}$) and missing energy ($E_{miss}$) distributions, we can distinguish the two-body break-up reaction, in which the residual proton-neutron system remains bound as a deuteron. In the $^3$H mirror case, under the exact same kinematic conditions, we do not identify a signature for a bound di-neutron with similar binding energy to that of the deuteron. We calculate exclusion limits as a function of the di-neutron binding energy and find that, for binding equivalent to the deuteron, the two-body break-up cross section on $^3$H is less than 0.9% of that on $^3$He in the measured kinematics at the 95% confidence level.


title: DRAN: Detailed Region-Adaptive Normalization for Conditional Image Synthesis

abstract: In recent years, conditional image synthesis has attracted growing attention due to its controllability in the image generation process. Although recent works have achieved realistic results, most of them have difficulty handling fine-grained styles with subtle details. To address this problem, a novel normalization module, named Detailed Region-Adaptive Normalization~(DRAN), is proposed. It adaptively learns both fine-grained and coarse-grained style representations. Specifically, we first introduce a multi-level structure, Spatiality-aware Pyramid Pooling, to guide the model to learn coarse-to-fine features. Then, to adaptively fuse different levels of styles, we propose Dynamic Gating, making it possible to adaptively fuse different levels of styles according to different spatial regions. Finally, we collect a new makeup dataset (Makeup-Complex dataset) that contains a wide range of complex makeup styles with diverse poses and expressions. To evaluate the effectiveness and show the general use of our method, we conduct a set of experiments on makeup transfer and semantic image synthesis. Quantitative and qualitative experiments show that equipped with DRAN, simple baseline models are able to achieve promising improvements in complex style transfer and detailed texture synthesis. Both the code and the proposed dataset will be available at https://github.com/Yueming6568/DRAN-makeup.git.


title: On the reliability of published findings using the regression discontinuity design in political science

abstract: The regression discontinuity (RD) design offers identification of causal effects under weak assumptions, earning it a position as a standard method in modern political science research. But identification does not necessarily imply that causal effects can be estimated accurately with limited data. In this paper, we highlight that estimation under the RD design involves serious statistical challenges and investigate how these challenges manifest themselves in the empirical literature in political science. We collect all RD-based findings published in top political science journals in the period 2009-2018. The distribution of published results exhibits pathological features; estimates tend to bunch just above the conventional level of statistical significance. A reanalysis of all studies with available data suggests that researcher discretion is not a major driver of these features. However, researchers tend to use inappropriate methods for inference, rendering standard errors artificially small. A retrospective power analysis reveals that most of these studies were underpowered to detect all but large effects. The issues we uncover, combined with well-documented selection pressures in academic publishing, cause concern that many published findings using the RD design may be exaggerated.


title: Human-centric Data Dissemination in the IoP: Large-scale Modeling and Evaluation

abstract: Data management using Device-to-Device (D2D) communications and opportunistic networks (ONs) is one of the main focuses of human-centric pervasive Internet services. In the recently proposed "Internet of People" paradigm, accessing relevant data dynamically generated in the environment nearby is one of the key services. Moreover, personal mobile devices become proxies of their human users while exchanging data in the cyber world and, thus, largely use ONs and D2D communications for exchanging data directly. Recently, researchers have successfully demonstrated the viability of embedding human cognitive schemes in data dissemination algorithms for ONs. In this paper, we consider one such scheme based on the recognition heuristic, a human decision-making scheme used to efficiently assess the relevance of data. While initial evidence about its effectiveness is available, the evaluation of its behaviour in large-scale settings is still unsatisfactory. To overcome these limitations, we have developed a novel hybrid modelling methodology, which combines an analytical model of data dissemination within small-scale communities of mobile users, with detailed simulations of interactions between different communities. This methodology allows us to evaluate the algorithm in large-scale city- and country-wide scenarios. Results confirm the effectiveness of cognitive data dissemination schemes, even when content popularity is very heterogenous.


title: Sublinear Time and Space Algorithms for Correlation Clustering via Sparse-Dense Decompositions

abstract: We present a new approach for solving (minimum disagreement) correlation clustering that results in sublinear algorithms with highly efficient time and space complexity for this problem. In particular, we obtain the following algorithms for $n$-vertex $(+/-)$-labeled graphs $G$: -- A sublinear-time algorithm that with high probability returns a constant approximation clustering of $G$ in $O(n\log^2{n})$ time assuming access to the adjacency list of the $(+)$-labeled edges of $G$ (this is almost quadratically faster than even reading the input once). Previously, no sublinear-time algorithm was known for this problem with any multiplicative approximation guarantee. -- A semi-streaming algorithm that with high probability returns a constant approximation clustering of $G$ in $O(n\log{n})$ space and a single pass over the edges of the graph $G$ (this memory is almost quadratically smaller than input size). Previously, no single-pass algorithm with $o(n^2)$ space was known for this problem with any approximation guarantee. The main ingredient of our approach is a novel connection to sparse-dense graph decompositions that are used extensively in the graph coloring literature. To our knowledge, this connection is the first application of these decompositions beyond graph coloring, and in particular for the correlation clustering problem, and can be of independent interest.


title: Strong Solutions for 1D Compressible Navier-Stokes/Allen-Cahn System with Phase Variable Dependent Viscosity

abstract: This paper is concerned with a non-isentropic compressible Navier-Stokes/Allen-Cahn system with phase variable dependent viscosity $\eta(\chi)=\chi^\alpha$ and temperature dependent heat-conductivity $\kappa(\theta)=\theta^\beta$. We show the global existence of strong solutions under some assumptions on growth exponent $\alpha$ and initial data. It is worth noting that the initial data could be large if $\alpha\ge0$ is small, and the growth exponent $\beta>0$ can be arbitrary large.


title: Deep Spatio-Temporal Wind Power Forecasting

abstract: Wind power forecasting has drawn increasing attention among researchers as the consumption of renewable energy grows. In this paper, we develop a deep learning approach based on encoder-decoder structure. Our model forecasts wind power generated by a wind turbine using its spatial location relative to other turbines and historical wind speed data. In this way, we effectively integrate spatial dependency and temporal trends to make turbine-specific predictions. The advantages of our method over existing work can be summarized as 1) it directly predicts wind power based on historical wind speed, without the need for prediction of wind speed first, and then using a transformation; 2) it can effectively capture long-term dependency 3) our model is more scalable and efficient compared with other deep learning based methods. We demonstrate the efficacy of our model on the benchmark real-world datasets.


title: Proton wave function in a water molecule: Breakdown of degeneration caused by interactions with the magnetic field of a Magnetic Resonance Imaging device

abstract: The concept of a Magnetic Resonance Imaging (MRI) device is based on the emission of radio waves produced by the protons of the hydrogen atoms in water molecules when placed in a constant magnetic field after they interact with a pulsed radio frequency (RF) current. When the RF field is turned on, the protons are brought to a spin excited state. When the RF field is turned off, the MRI sensors are able to detect the energy released as the protons realign their spins with the magnetic field. In this work we provide a simple model to describe the basic physical mechanism responsible for the operation of MRI devices. We model the water molecule in terms of a central force problem, where the protons move around the (unstructured) doubly negatively charged oxygen atom. First, we employ an analytical treatment to obtain the system's wave function as well as its energy levels, which we show are degenerate. Next, the energy levels from the water molecule are studied in the presence of a uniform external magnetic field. As a result, they get shifted and the degeneration is lifted. We provide numerical results for a magnetic field strength commonly used in MRI devices.


title: Multispectral CCD-in-CMOS Time Delay Integration imager for high resolution Earth observation

abstract: Many future small satellite missions are aimed to provide low-cost remote sensing data at unprecedented revisit rates, with a ground resolution of less than one meter. This requires high resolution, fast and sensitive line-scan imagers operating at low power consumption and ideally featuring spectral sensitivity. In this paper we present comprehensive characterization results of our 7 band Back-Side Illuminated (BSI) CCD-in-CMOS sensor with a pixel pitch of 5.4 um. We have extensively characterized the key performance parameters of our CCD-in-CMOS sensor, such as quantum efficiency (QE), full well capacity (FWC), read noise, conversion gain, non-linearity, dark current etc. Novelty of this device is the combination of 7 TDI bands on the same imager allowing simultaneous multispectral TDI capture. Glass-based broadband filters with a typical band-pass width of about 100 nm have been developed and glued together to form a filter assembly of 6 band-pass filters and one panchromatic channel. Multispectral capability of this sensor is particularly interesting for Low Earth Observation (LEO) applications such as environmental monitoring, precision agriculture, disaster detection and monitoring. To highlight its ad-vantages for use in vegetation observation, we demonstrated a fake leaf and a real leaf imaging using a 7 band BSI sensor with integrated filters operating in 7-band mode at 15 kHz.


title: A Mass Conserving Mixed hp-FEM Scheme for Stokes Flow. Part III: Implementation and Preconditioning

abstract: This is the third part in a series on a mass conserving, high order, mixed finite element method for Stokes flow. In this part, we study a block-diagonal preconditioner for the indefinite Schur complement system arising from the discretization of the Stokes equations using these elements. The underlying finite element method is uniformly stable in both the mesh size h and polynomial order p, and we prove bounds on the eigenvalues of the preconditioned system which are independent of h and grow modestly in p. The analysis relates the Schur complement system to an appropriate variational setting with subspaces for which exact sequence properties and inf-sup stability hold. Several numerical examples demonstrate agreement with the theoretical results.


title: A verified algebraic representation of Cairo program execution

abstract: Cryptographic interactive proof systems provide an efficient and scalable means of verifying the results of computation on blockchain. A prover constructs a proof, off-chain, that the execution of a program on a given input terminates with a certain result. The prover then publishes a certificate that can be verified efficiently and reliably modulo commonly accepted cryptographic assumptions. The method relies on an algebraic encoding of execution traces of programs. Here we report on a verification of the correctness of such an encoding of the Cairo model of computation with respect to the STARK interactive proof system, using the Lean 3 proof assistant.


title: Analyse or Transmit: Utilising Correlation at the Edge with Deep Reinforcement Learning

abstract: Millions of sensors, cameras, meters, and other edge devices are deployed in networks to collect and analyse data. In many cases, such devices are powered only by Energy Harvesting(EH) and have limited energy available to analyse acquired data. When edge infrastructure is available, a device has a choice: to perform analysis locally or offload the task to other resource-rich devices such as cloudlet servers. However, such a choice carries a price in terms of consumed energy and accuracy. On the one hand, transmitting raw data can result in a higher energy cost in comparison to the required energy to process data locally. On the other hand, performing data analytics on servers can improve the task's accuracy. Additionally, due to the correlation between information sent by multiple devices, accuracy might not be affected if some edge devices decide to neither process nor send data and preserve energy instead. For such a scenario, we propose a Deep Reinforcement Learning (DRL) based solution capable of learning and adapting the policy to the time-varying energy arrival due to EH patterns. We leverage two datasets, one to model energy an EH device can collect and the other to model the correlation between cameras. Furthermore, we compare the proposed solution performance to three baseline policies. Our results show that we can increase accuracy by 15% in comparison to conventional approaches while preventing outages.


title: PINNup: Robust neural network wavefield solutions using frequency upscaling and neuron splitting

abstract: Solving for the frequency-domain scattered wavefield via physics-informed neural network (PINN) has great potential in seismic modeling and inversion. However, when dealing with high-frequency wavefields, its accuracy and training cost limits its applications. Thus, we propose a novel implementation of PINN using frequency upscaling and neuron splitting, which allows the neural network model to grow in size as we increase the frequency while leveraging the information from the pre-trained model for lower-frequency wavefields, resulting in fast convergence to high-accuracy solutions. Numerical results show that, compared to the commonly used PINN with random initialization, the proposed PINN exhibits notable superiority in terms of convergence and accuracy and can achieve neuron based high-frequency wavefield solutions with a two-hidden-layer model.


title: EAS observation conditions in the SPHERE-2 balloon experiment

abstract: The SPHERE project studies primary cosmic rays by detection of the Cherenkov light of extensive air showers reflected from the snowy surface of the earth. Measurements with the aerial-based detector SPHERE-2 were performed in 2011-2013. The detector was lifted by the balloon at altitudes up to 900 m above snowed surface of Lake Baikal, Russia. The results of the experiment are summarized now in a series of papers that opens with this article. An overview of the SPHERE-2 detector telemetry monitoring systems is presented along with the analysis of the measurements conditions including atmosphere profile. The analysis of the detector state and environment atmosphere conditions monitoring provided various cross-checks of detector calibration, positioning and performance.


title: Unification of Gauge Symmetries... including their breaking

abstract: In this talk, we present a minimal viable scenario that unifies the gauge symmetries of the Standard Model (SM) and their breaking sector. Our Gauge-Higgs Grand Unification setup employs 5D warped space with a $SU(6)$ bulk gauge field that includes both a $SU(5)$ grand unified theory (GUT) and a Higgs sector as a scalar component of the 5D vector field, solving the hierarchy problem. By appropriately breaking the gauge symmetry on the boundaries of the extra dimension the issue of light exotic new states, appearing generically in such models, is eliminated and the SM fermion spectrum is naturally reproduced. The Higgs potential is computed at one-loop, finding straightforward solutions with a realistic $m_h = 125$ GeV. The problem of proton decay is addressed by showing that baryon number is a hidden symmetry of the model. The presence of a scalar leptoquark and a scalar singlet is highlighted, which might play a role in solving further problems of the SM, allowing for example for electroweak baryogenesis. Finally, the $X$ and $Y$ gauge bosons from $SU(5)$ GUTs are found at collider accessible masses, opening a window to the unification structure at low energies.


title: New Solution based on Hodge Decomposition for Abstract Games

abstract: This paper proposes Hodge Potential Choice (HPC), a new solution for abstract games with irreflexive dominance relations. This solution is formulated by involving geometric tools like differential forms and Hodge decomposition onto abstract games. We provide a workable algorithm for the proposed solution with a new data structure of abstract games. From the view of gaming, HPC overcomes several weaknesses of conventional solutions. HPC coincides with Copeland Choice in complete cases and can be extended to slove games with marginal strengths. It will be proven that the Hodge potential choice possesses three prevalent axiomatic properties: neutrality, strong monotonicity, dominance cycle s reversing independence, and sensitivity to mutual dominance. To compare the HPC with Copeland Choice in large samples of games, we design digital experiments with randomly generated abstract games with different sizes and completeness. The experimental results present the advantage of HPC in the statistical sense.


title: On a class of non-Hermitian Hamiltonians with tridiagonal matrix representation

abstract: We show that some non-Hermitian Hamiltonian operators with tridiagonal matrix representation may be quasi Hermitian or similar to Hermitian operators. In the class of Hamiltonian operators discussed here the transformation is given by a Hermitian, positive-definite, diagonal operator. We show that there is an important difference between open boundary conditions and periodic ones. We illustrate the theoretical results by means of two simple, widely used, models.


title: The Berndt-Bhargava-Garvan Transfer Principle for signature three

abstract: We offer a new proof for the Berndt-Bhargava-Garvan Transfer Principle that connects the signature-three elliptic theory of Ramanujan to the classical elliptic theory.


title: gwverse: a template for a new generic Geographically Weighted Rpackage

abstract: GWR is a popular approach for investigating the spatial variation in relationships between response and predictor variables, and critically for investigating and understanding process spatial heterogeneity. The geographically weighted (GW) framework is increasingly used to accommodate different types of models and analyses reflecting a wider desire to explore spatial variation in model parameters or components. However the growth in the use of GWR and different GW models has only been partially supported by package development in both R and Python, the major coding environments for spatial analysis. The result is that refinements have been inconsistently included (if at all) within GWR and GW functions in any given package. This paper outlines the structure of a new `gwverse` package, that will over time replace `GWmodel`, that takes advantage of recent developments in the composition of complex, integrated packages. It conceptualises `gwverse` as having a modular structure, that separates core GW functionality and applications such as GWR. It adopts a function factory approach, in which bespoke functions are created and returned to the user based on user-defined parameters. The paper introduces two demonstrator modules that can be used to undertake GWR and identifies a number of key considerations and next steps.


title: Topology of orientational defects in confined smectic liquid crystals

abstract: We propose a general formalism to characterize orientational frustration of smectic liquid crystals in confinement by interpreting the emerging networks of grain boundaries as objects with a topological charge. In a formal idealization, this charge is distributed in point-like units of quarter-integer magnitude, which we identify with tetratic disclinations located at the end points and nodes. This coexisting nematic and tetratic order is analyzed with the help of extensive Monte Carlo simulations for a broad range of two-dimensional confining geometries as well as colloidal experiments, showing how the observed defect networks can be universally reconstructed from simple building blocks. We further find that the curvature of the confining wall determines the anchoring behavior of grain boundaries, such that the number of nodes in the emerging networks and the location of their end points can be tuned by changing the number and smoothness of corners, respectively.


title: Higher-order least squares: assessing partial goodness of fit of linear causal models

abstract: We introduce a simple diagnostic test for assessing the overall or partial goodness of fit of a linear causal model with errors being independent of the covariates. In particular, we consider situations where hidden confounding is potentially present. We develop a method and discuss its capability to distinguish between covariates that are confounded with the response by latent variables and those that are not. Thus, we provide a test and methodology for partial goodness of fit. The test is based on comparing a novel higher-order least squares principle with ordinary least squares. In spite of its simplicity, the proposed method is extremely general and is also proven to be valid for high-dimensional settings.


title: Activation Functions in Deep Learning: A Comprehensive Survey and Benchmark

abstract: Neural networks have shown tremendous growth in recent years to solve numerous problems. Various types of neural networks have been introduced to deal with different types of problems. However, the main goal of any neural network is to transform the non-linearly separable input data into more linearly separable abstract features using a hierarchy of layers. These layers are combinations of linear and nonlinear functions. The most popular and common non-linearity layers are activation functions (AFs), such as Logistic Sigmoid, Tanh, ReLU, ELU, Swish and Mish. In this paper, a comprehensive overview and survey is presented for AFs in neural networks for deep learning. Different classes of AFs such as Logistic Sigmoid and Tanh based, ReLU based, ELU based, and Learning based are covered. Several characteristics of AFs such as output range, monotonicity, and smoothness are also pointed out. A performance comparison is also performed among 18 state-of-the-art AFs with different networks on different types of data. The insights of AFs are presented to benefit the researchers for doing further research and practitioners to select among different choices. The code used for experimental comparison is released at: \url{https://github.com/shivram1987/ActivationFunctions}.


title: An Energy Efficient Health Monitoring Approach with Wireless Body Area Networks

abstract: Wireless Body Area Networks (WBANs) comprise a network of sensors subcutaneously implanted or placed near the body surface and facilitate continuous monitoring of health parameters of a patient. Research endeavours involving WBAN are directed towards effective transmission of detected parameters to a Local Processing Unit (LPU, usually a mobile device) and analysis of the parameters at the LPU or a back-end cloud. An important concern in WBAN is the lightweight nature of WBAN nodes and the need to conserve their energy. This is especially true for subcutaneously implanted nodes that cannot be recharged or regularly replaced. Work in energy conservation is mostly aimed at optimising the routing of signals to minimise energy expended. In this paper, a simple yet innovative approach to energy conservation and detection of alarming health status is proposed. Energy conservation is ensured through a two-tier approach wherein the first tier eliminates `uninteresting' health parameter readings at the site of a sensing node and prevents these from being transmitted across the WBAN to the LPU. A reading is categorised as uninteresting if it deviates very slightly from its immediately preceding reading and does not provide new insight on the patient's well being. In addition to this, readings that are faulty and emanate from possible sensor malfunctions are also eliminated. These eliminations are done at the site of the sensor using algorithms that are light enough to effectively function in the extremely resource-constrained environments of the sensor nodes. We notice, through experiments, that this eliminates and thus reduces around 90% of the readings that need to be transmitted to the LPU leading to significant energy savings. Furthermore, the proper functioning of these algorithms in such constrained environments is confirmed and validated over a hardware simulation set up. The second tier of assessment includes a proposed anomaly detection model at the LPU that is capable of identifying anomalies from streaming health parameter readings and indicates an adverse medical condition. In addition to being able to handle streaming data, the model works within the resource-constrained environments of an LPU and eliminates the need of transmitting the data to a back-end cloud, ensuring further energy savings. The anomaly detection capability of the model is validated using data available from the critical care units of hospitals and is shown to be superior to other anomaly detection techniques.


title: Diamond-based Detection Systems for Accurate Pulsed X-rays Diagnostics in Radiotherapy

abstract: The widespread diffusion of precision radiotherapy techniques, geared toward the release of larger dose gradients in shorter time frames, is leading to new challenges in dosimetry. Accurate dose measurements are essential to check for beam anomalies and inaccuracies to ensure treatment efficacy and patient safety during radiotherapy. This work describes the main features of a diamond dosimeter coupled to an extremely compact front-end electronics. The detection system was tested under the X-ray pulses generated by a medical LINAC for both the 6 MV and the 18 MV accelerating voltages. Located in the LINAC's bunker, it eliminates the need for a long cable connection between the detector and the electronics, detrimental for the system response speed. Signal acquisition was performed synchronously with the impinging X-ray pulses with a sampling period as low as 20 us, allowing for a real-time beam monitoring. The dosimeter demonstrated a very good stability despite the high value of the absorbed dose during the performed experiments (about 100 Gy). The measured dose-per-pulse values of 278 uGy and 556 uGy at 6 MV and 18 MV, respectively, are in excellent agreement with the nominal values expected for the LINAC apparatus used for the tests. In addition to single-pulse measurements, fundamental for dynamic radiotherapy, the proposed system also allows for the calculation of both the total collected charge and the photocurrent generated by the detector. In this regards, despite the compactness, it demonstrates its effectiveness as a tool for source diagnostics in terms of both beam intensity and emission timing.

categories: physics.med-ph physics.app-ph physics.class-ph physics.ins-det


title: Analytic energy gradients for the driven similarity renormalization group multireference second-order perturbation theory

abstract: We derive analytic energy gradients of the driven similarity renormalization group (DSRG) multireference second-order perturbation theory (MRPT2) using the method of Lagrange multipliers. In the Lagrangian, we impose constraints for a complete-active-space self-consistent-field reference wave function and the semicanonical orthonormal molecular orbitals. Solving the associated Lagrange multipliers is found to share the same asymptotic scaling of a single DSRG-MRPT2 energy computation. A pilot implementation of the DSRG-MRPT2 analytic gradients is used to optimize the geometry of the singlet and triplet states of p-benzyne. The equilibrium bond lengths and angles are similar to those computed via other MRPT2s and Mukherjee's multireference coupled cluster theory. An approximate DSRG-MRPT2 method that neglects the contributions of three-body density cumulant is found to introduce negligible errors in the geometry of p-benzyne, lending itself to a promising low-cost approach for molecular geometry optimizations using large active spaces.


title: Vision-Guided Quadrupedal Locomotion in the Wild with Multi-Modal Delay Randomization

abstract: Developing robust vision-guided controllers for quadrupedal robots in complex environments, with various obstacles, dynamical surroundings and uneven terrains, is very challenging. While Reinforcement Learning (RL) provides a promising paradigm for agile locomotion skills with vision inputs in simulation, it is still very challenging to deploy the RL policy in the real world. Our key insight is that aside from the discrepancy in the domain gap, in visual appearance between the simulation and the real world, the latency from the control pipeline is also a major cause of difficulty. In this paper, we propose Multi-Modal Delay Randomization (MMDR) to address this issue when training RL agents. Specifically, we simulate the latency of real hardware by using past observations, sampled with randomized periods, for both proprioception and vision. We train the RL policy for end-to-end control in a physical simulator without any predefined controller or reference motion, and directly deploy it on the real A1 quadruped robot running in the wild. We evaluate our method in different outdoor environments with complex terrains and obstacles. We demonstrate the robot can smoothly maneuver at a high speed, avoid the obstacles, and show significant improvement over the baselines. Our project page with videos is at https://mehooz.github.io/mmdr-wild/.


title: A Riemannian Approach to Multivariate Geostatistical Modeling

abstract: In geosciences, the use of classical Euclidean methods is unsuitable for treating and analyzing some types of data, as this may not belong to a vector space. This is the case for correlation matrices, belonging to a subfamily of symmetric positive definite matrices, which in turn form a cone shape Riemannian manifold. We propose two novel applications for dealing with the problem of accounting with the non-linear behavior usually presented on multivariate geological data by exploiting the manifold features of correlations matrices. First, we employ an extension for the linear model of coregionalization (LMC) that alters the linear mixture, which is assumed fixed on the domain, and making it locally varying according to the local strength in the dependency of the coregionalized variables. The main challenge, once this relaxation on the LMC is assumed, is to solve appropriately the interpolation of the different known correlation matrices throughout the domain, in a reliable and coherent fashion. The present work adopts the non-euclidean framework to achieve our objective by locally averaging and interpolating the correlations between the variables, retaining the intrinsic geometry of correlation matrices. A second application deals with the problem of clustering of multivariate data.


title: DORA: Distributed Online Risk-Aware Explorer

abstract: Exploration of unknown environments is an important challenge in the field of robotics. While a single robot can achieve this task alone, evidence suggests it could be accomplished more efficiently by groups of robots, with advantages in terms of terrain coverage as well as robustness to failures. Exploration can be guided through belief maps, which provide probabilistic information about which part of the terrain is interesting to explore (either based on risk management or reward). This process can be centrally coordinated by building a collective belief map on a common server. However, relying on a central processing station creates a communication bottleneck and single point of failure for the system. In this paper, we present Distributed Online Risk-Aware (DORA) Explorer, an exploration system that leverages decentralized information sharing to update a common risk belief map. DORA Explorer allows a group of robots to explore an unknown environment discretized as a 2D grid with obstacles, with high coverage while minimizing exposure to risk, effectively reducing robot failures


title: Silence is bliss in a Platonic relationship

abstract: We describe an effective active cloaking strategy for the scalar Helmholtz equation in three dimensions where multipole active sources are located at the vertices of the Platonic solids. A "silent zone" is created interior to the imaginary Platonic solid and only the incident field remains in a defined region exterior to the silent zone and active source configuration. This distribution of the sources ensures that the implementation of the cloaking strategy is extremely efficient. In particular, once the multipole source amplitudes required at a single source location are determined, the other source amplitudes can be calculated by simple post-processing involving multiplication of the multipole source vector by a rotation matrix. The general nature of the problem means that the technique is relevant to any scalar wave field, including both acoustics and electromagnetism.


title: Ab Initio Thermodynamic Study of PbI$_2$ and CH$_3$NH$_3$PbI$_3$ Surfaces in Reaction with CH$_3$NH$_2$ Gas for Perovskite Solar Cells

abstract: Hybrid organic-inorganic halide perovskites for photovoltaics have attracted research interest due to their unique material properties, but suffered from poor material stability. In this work, we investigated the surface phase diagrams of PbI$_2$ and cubic CH$_3$NH$_3$PbI$_3$ (MAPbI$_3$) for a better understanding of precursor effect on perovskite synthesis via solid-gas reaction, by density functional theory calculations combined with thermodynamics. Using the devised slab models of PbI$_2$(001) and MAPbI$_3$(100), (110) and (111) surfaces with perfect and various vacant defect terminations, we calculated their formation energies and adsorption energies of CH$_3$NH$_2$ molecule on the PbI$_2$(001) surfaces under different synthesis conditions of temperature, pressure and pH via chemical potentials of species. Our calculations revealed that the adsorption can be facilitated by including HI or NH$_4$I molecule by dissociation of this additive and formation of CH$_3$NH$_3^+$ cation or NH$_3$-CH$_3$NH$_3^+$ complex, which is beneficial for conversion of PbI$_2$ to MAPbI$_3$ via solid-gas reaction. Furthermore, we found that among different perovskite MAPbI$_3$ surfaces, the MA-terminated (110) and MAI-terminated (100) surfaces are placed on the thermodynamically stable region of chemical potentials at pH values of 1 and 6, being agreed well with the experimental findings. We believe this work gives a fundamental understanding of solid-gas reaction for high-crystallinity perovskite synthesis towards perovskite solar cells with improved stability.


title: Coulomb-like Model for International Trade Flow and Derivation of Distribution Function for Trade Flow Strength

abstract: To describe international trade flows, we propose the coulomb force formulation, in which the magnitude of the charge represents gross domestic product (GDP) and the distance between countries is the bilateral distance, the product of spatial distance and "dielectric constant," rather than the spatial distance as used in the gravitation model, allowing it to be time dependent. The "dielectric constant" is influenced by factors such as warfare, transportation disruptions, trade agreements, social, geography, politics, culture, and others. The GDP and distance power parameters were estimated using data from high-GDP countries' export-import transactions. We also developed a trade strength distribution equation that fits World Bank data reasonably well over a decade.


title: The Impact of Network Design Interventions on CPS Security

abstract: We study a game-theoretic model of the interactions between a Cyber-Physical System's (CPS) operator (the defender) against an attacker who launches stepping-stone attacks to reach critical assets within the CPS. We consider that, in addition to optimally allocating its security budget to protect the assets, the defender may choose to modify the CPS through network design interventions. In particular, we propose and motivate four ways in which the defender can introduce additional nodes in the CPS: these nodes may be intended as additional safeguards, be added for functional or structural redundancies, or introduce additional functionalities in the system. We analyze the security implications of each of these design interventions, and evaluate their impacts on the security of an automotive network as our case study. We motivate the choice of the attack graph for this case study and elaborate how the parameters in the resulting security game are selected using the CVSS metrics and the ISO-26262 ASIL ratings as guidance. We then use numerical experiments to verify and evaluate how our proposed network interventions may be used to guide improvements in automotive security.


title: Graphene-based sponges for electrochemical degradation of persistent organic contaminants

abstract: Graphene-based sponges doped with atomic nitrogen and boron were applied for the electrochemical degradation of persistent organic contaminants in one-pass, flow-through mode, and in a low-conductivity supporting electrolyte. The B-doped anode and N-doped cathode was capable of >90% contaminant removal at the geometric anodic current density of 173 A m2. The electrochemical degradation of contaminants was achieved via the direct electron transfer, the anodically formed O3, and by the OH radicals formed by the decomposition of H2O2 produced at the cathode. The identified transformation products of iopromide show that the anodic cleavage of all three C-I bonds at the aromatic ring was preferential over scissions at the alkyl side chains, suggesting a determining role of the pi-pi interactions with the graphene surface. In the presence of 20 mM sodium chloride (NaCl), the current efficiency for chlorine production was <0.04%, and there was no chlorate and perchlorate formation, demonstrating a very low electrocatalytic activity of the graphene-based sponge anode towards chloride. Graphene-based sponges were produced using a low-cost, bottom-up method that allows easy introduction of dopants and functionalization of the reduced graphene oxide coating, and thus tailoring of the material for the removal of specific contaminants.


title: Modeling the microlensing of circumbinary systems

abstract: Gravitational microlensing is one of the methods to detect exoplanets; planets outside our solar system. Here we focus on theoretical modeling of three lens systems and in particular circumbinary systems. Circumbinary systems include two stars and a planet and are estimated to make up a sizable portion of all exoplanets. Extending a method developed for binary lenses to the three lens case, we explore the parameter space of circumbinary systems producing exact magnification maps and light curves.

categories: astro-ph.EP astro-ph.GA astro-ph.IM astro-ph.SR


title: Postural Orthostatic Tachycardia Syndrome explained using abaroreflex response model

abstract: Recent studies have shown that Postural Orthostatic Tachycardia Syndrome (POTS) patients have abnormal low frequency $\approx 0.1$ Hz blood pressure and heart rate dynamics. These dynamics are attributed to the baroreflex and can give insight into the mechanistic causes which are the basis for proposed subgroups of POTS. In this study we develop a baroreflex model replicating the low-frequency dynamics observed in POTS patient data as well as represent subgroups of POTS in our model. We utilize signal processing to quantify the effects that model parameters have on low-frequency oscillations. Results show that key physiological parameters that represent the hypothesized causes are central to our model's ability to reproduce observed dynamics from patient data.


title: Asymptotic profiles of zero points of solutions to the heat equation

abstract: In this paper, we consider the asymptotic profiles of zero points for the spatial variable of the solutions to the heat equation. By giving suitable conditions for the initial data, we prove the existence of zero points by extending the high-order asymptotic expansion theory for the heat equation. This reveals a previously unknown asymptotic profile of zero points diverging at $O(t)$. In a one-dimensional spatial case, we show the zero point's second and third-order asymptotic profiles in a general situation. We also analyze a zero-level set in high-dimensional spaces and obtain results that extend the results for the one-dimensional spatial case.


title: Unstable diffusion in social networks

abstract: How and to what extent will new activities spread through social ties? Here, we develop a more sophisticated framework than the standard mean-field approach to describe the diffusion dynamics of multiple activities on complex networks. We show that the diffusion of multiple activities follows a saddle path and can be highly unstable. In particular, when the two activities are sufficiently substitutable, either of them would dominate the other by chance even if they are equally attractive ex ante. When such symmetry-breaking occurs, any average-based approach cannot correctly calculate the Nash equilibrium - the steady state of an actual diffusion process.


title: Chiral criticality and the role of repulsive interactions in hot hadronic matter

abstract: We explore fluctuations of the net-baryon number density in hadronic matter at vanishing chemical potential. We discuss the interplay between chiral dynamics and repulsive interactions and their influence on the properties of these fluctuations near the chiral crossover in the parity doublet model. We find that the characteristic properties of cumulants are entirely linked to the critical chiral dynamics and cannot be reproduced in phenomenological models, which account only for repulsive interactions. Consequently, a description of the higher-order cumulants of the net-baryon density requires a self-consistent treatment of the chiral in-medium effects and repulsive interactions.


title: Stout-smearing, gradient flow and $c_{\text{SW}}$ at one loop order

abstract: The one-loop determination of the coefficient $c_\text{SW}$ of the Wilson quark action has been useful to push the leading cut-off effects for on-shell quantities to $\mathcal{O}(\alpha^2 a)$ and, in conjunction with non-perturbative determinations of $c_\text{SW}$, to $\mathcal{O}(a^2)$, as long as no link-smearing is employed. These days it is common practice to include some overall link-smearing into the definition of the fermion action. Unfortunately, in this situation only the tree-level value $c_\text{SW}^{(0)}=1$ is known, and cut-off effects start at $\mathcal{O}(\alpha a)$. We present some general techniques for calculating one loop quantities in lattice perturbation theory which continue to be useful for smeared-link fermion actions. Specifically, we discuss the application to the 1-loop improvement coefficient $c_\text{SW}^{(1)}$ for overall stout-smeared Wilson fermions.


title: Robust Temporal Ensembling for Learning with Noisy Labels

abstract: Successful training of deep neural networks with noisy labels is an essential capability as most real-world datasets contain some amount of mislabeled data. Left unmitigated, label noise can sharply degrade typical supervised learning approaches. In this paper, we present robust temporal ensembling (RTE), which combines robust loss with semi-supervised regularization methods to achieve noise-robust learning. We demonstrate that RTE achieves state-of-the-art performance across the CIFAR-10, CIFAR-100, ImageNet, WebVision, and Food-101N datasets, while forgoing the recent trend of label filtering and/or fixing. Finally, we show that RTE also retains competitive corruption robustness to unforeseen input noise using CIFAR-10-C, obtaining a mean corruption error (mCE) of 13.50% even in the presence of an 80% noise ratio, versus 26.9% mCE with standard methods on clean data.


title: Discrepancy and rectifiability of almost linearly repetitive Delone sets

abstract: We extend a discrepancy bound of Lagarias and Pleasants for local weight distributions on linearly repetitive Delone sets and show that a similar bound holds also for the more general case of Delone sets without finite local complexity if linear repetitivity is replaced by $\varepsilon$-linear repetitivity. As a result we establish that Delone sets that are $\varepsilon$-linear repetitive for some sufficiently small $\varepsilon$ are rectifiable, and that incommensurable multiscale substitution tilings are never almost linearly repetitive.


title: Fitting strategies of accretion column models and application to the broadband spectrum of Cen X-3

abstract: Due to the complexity of modeling the radiative transfer inside the accretion columns of neutron star binaries, their X-ray spectra are still commonly described with phenomenological models, for example, a cutoff power law. While the behavior of these models is well understood and they allow for a comparison of different sources and studying source behavior, the extent to which the underlying physics can be derived from the model parameters is very limited. During recent years, several physically motivated spectral models have been developed to overcome these limitations. Their application, however, is generally computationally much more expensive and they require a high number of parameters which are difficult to constrain. Previous works have presented an analytical solution to the radiative transfer equation inside the accretion column assuming a velocity profile that is linear in the optical depth. An implementation of this solution that is both fast and accurate enough to be fitted to observed spectra is available as a model in XSPEC. The main difficulty of this implementation is that some solutions violate energy conservation and therefore have to be rejected by the user. We propose a novel fitting strategy that ensures energy conservation during the $\chi^2$-minimization which simplifies the application of the model considerably. We demonstrate this approach as well a study of possible parameter degeneracies with a comprehensive Markov-chain Monte Carlo analysis of the complete parameter space for a combined NuSTAR and Swift/XRT dataset of Cen X-3. The derived accretion-flow structure features a small column radius of $\sim$63 m and a spectrum dominated by bulk-Comptonization of bremsstrahlung seed photons, in agreement with previous studies.


title: Generation of topologically complex three-dimensional electron beams in a plasma photocathode

abstract: Laser-triggered ionization injection is a promising way of generating controllable high-quality electrons in plasma-based acceleration. We show that ionization injection of electrons into a fully nonlinear plasma wave wake using a laser pulse comprising of one or more Laguerre-Gaussian modes with combinations of spin and orbital angular momentum can generate exotic three-dimensional (3D) spatial distributions of high-quality relativistic electrons. The phase dependent residual momenta and initial positions of the ionized electrons are encoded into their final phase space distributions, leading to complex spatiotemporal structures. The structures are formed as a result of the transverse (betatron) and longitudinal (phase slippage and energy gain) dynamics of the electrons in the wake immediately after the electrons are injected. Theoretical analysis and 3D simulations verify this mapping process leads to the generation of these complex topological beams. These beams may trigger novel beam-plasma interactions as well as produce coherent radiation with orbital angular momentum when sent through a resonant undulator.


title: Implicit Generative Copulas

abstract: Copulas are a powerful tool for modeling multivariate distributions as they allow to separately estimate the univariate marginal distributions and the joint dependency structure. However, known parametric copulas offer limited flexibility especially in high dimensions, while commonly used non-parametric methods suffer from the curse of dimensionality. A popular remedy is to construct a tree-based hierarchy of conditional bivariate copulas. In this paper, we propose a flexible, yet conceptually simple alternative based on implicit generative neural networks. The key challenge is to ensure marginal uniformity of the estimated copula distribution. We achieve this by learning a multivariate latent distribution with unspecified marginals but the desired dependency structure. By applying the probability integral transform, we can then obtain samples from the high-dimensional copula distribution without relying on parametric assumptions or the need to find a suitable tree structure. Experiments on synthetic and real data from finance, physics, and image generation demonstrate the performance of this approach.


title: Stochastic Primitive Equations with Horizontal Viscosity and Diffusivity

abstract: We establish the existence and uniqueness of pathwise strong solutions to the stochastic 3D primitive equations with only horizontal viscosity and diffusivity driven by transport noise on a cylindrical domain $M=(-h,0) \times G$, $G\subset \mathbb{R}^2$ bounded and smooth, with the physical Dirichlet boundary conditions on the lateral part of the boundary. Compared to the deterministic case where the uniqueness of $z$-weak solutions holds in $L^2$, more regular initial data are necessary to establish uniqueness in the anisotropic space $H^1_z L^2_{xy}$ so that the existence of local pathwise solutions can be deduced from the Gy\"{o}ngy-Krylov theorem. Global existence is established using the logarithmic Sobolev embedding, the stochastic Gronwall lemma and an iterated stopping time argument.


title: An Expert System for Redesigning Software for Cloud Applications

abstract: Cloud-based software has many advantages. When services are divided into many independent components, they are easier to update. Also, during peak demand, it is easier to scale cloud services (just hire more CPUs). Hence, many organizations are partitioning their monolithic enterprise applications into cloud-based microservices. Recently there has been much work using machine learning to simplify this partitioning task. Despite much research, no single partitioning method can be recommended as generally useful. More specifically, those prior solutions are "brittle"; i.e. if they work well for one kind of goal in one dataset, then they can be sub-optimal if applied to many datasets and multiple goals. In order to find a generally useful partitioning method, we propose DEEPLY. This new algorithm extends the CO-GCN deep learning partition generator with (a) a novel loss function and (b) some hyper-parameter optimization. As shown by our experiments, DEEPLY generally outperforms prior work (including CO-GCN, and others) across multiple datasets and goals. To the best of our knowledge, this is the first report in SE of such stable hyper-parameter optimization. To aid reuse of this work, DEEPLY is available on-line at https://bit.ly/2WhfFlB.


title: Hyperbolic 3-manifolds of low cusp volume

abstract: We classify the complete hyperbolic 3-manifolds admitting a maximal cusp of volume at most 2.62. We use this to show that the figure-8 knot complement is the unique 1-cusped hyperbolic 3-manifold with nine or more non-hyperbolic fillings; to show that the figure-8 knot complement and its sister are the unique hyperbolic 3-manifolds with minimal volume maximal cusps; and to extend results on determining low volume closed and cusped hyperbolic 3-manifolds.


title: Prospects for measuring off-axis spins of binary black holes with Plus-era gravitational-wave detectors

abstract: The mass and spin properties of binary black holes (BBHs) inferred from their gravitational-wave signatures reveal important clues about how these systems form. BBHs originating from isolated binary evolution are expected to have spins preferentially aligned with their orbital angular momentum, whereas there is no such preference in binaries formed via dynamical assembly. The fidelity with which near-future gravitational-wave detectors can measure off-axis spins will have implications for the study of BBH formation channels. In this work, we examine the degree to which the Advanced LIGO Plus (A+) and Advanced Virgo Plus (AdV+) interferometric detectors can measure both aligned and misaligned spins. We compare spin resolution between the LIGO-Virgo network operating at either A+/AdV+ ("Plus") sensitivity or Advanced-era design ("Design") sensitivity using simulated BBH gravitational-wave signals injected into synthetic detector noise. The signals are distributed over the mass-spin parameter space of likely BBH systems, accounting for the effects of precession and higher-order modes. We find that the Plus upgrades yield significant improvements in spin estimation for systems with unequal masses and moderate or large spins. Using simulated signals modelled after different types of hierarchical BBH mergers, we also conclude that the Plus detector network will yield substantially improved spin estimates for 1G+2G binaries compared to the Design network.