title: MHD-shock structures of astrospheres: $\lambda$ Cephei-like astrospheres

abstract: The interpretation of recent observations of bow shocks around O-stars and the creation of corresponding models require a detailed understanding of the associated (magneto-)hydrodynamic structures. We base our study on three-dimensional numerical magneto-)hydrodynamical models, which are analyzed using the dynamically relevant parameters, in particular, the (magneto)sonic Mach numbers. The analytic Rankine-Hugoniot relation for HD and MHD are compared with those obtained by the numerical model. In that context we also show that the only distance which can be approximately determined is that of the termination shock, if it is a hydrodynamical shock. For MHD shocks the stagnation point does not, in general, lie on the inflow line, which is the line parallel to the inflow vector and passing through the star. Thus an estimate via the Bernoulli equation as in the HD case is, in general, not possible. We also show that in O-star astrospheres, distinct regions exist in which the fast, slow, Alfv\'enic, and sonic Mach numbers become lower than one, implying sub-slow magnetosonic as well as sub-fast and sub-sonic flows. Nevertheless, the analytic MHD Rankine Hugoniot relations can be used for further studies of turbulence and cosmic ray modulation.


title: Investigating the Compositional Structure Of Deep Neural Networks

abstract: The current understanding of deep neural networks can only partially explain how input structure, network parameters and optimization algorithms jointly contribute to achieve the strong generalization power that is typically observed in many real-world applications. In order to improve the comprehension and interpretability of deep neural networks, we here introduce a novel theoretical framework based on the compositional structure of piecewise linear activation functions. By defining a direct acyclic graph representing the composition of activation patterns through the network layers, it is possible to characterize the instances of the input data with respect to both the predicted label and the specific (linear) transformation used to perform predictions. Preliminary tests on the MNIST dataset show that our method can group input instances with regard to their similarity in the internal representation of the neural network, providing an intuitive measure of input complexity.


title: Pandora's Box Problem with Order Constraints

abstract: The Pandora's Box Problem, originally formalized by Weitzman in 1979, models selection from set of random, alternative options, when evaluation is costly. This includes, for example, the problem of hiring a skilled worker, where only one hire can be made, but the evaluation of each candidate is an expensive procedure. Weitzman showed that the Pandora's Box Problem admits an elegant, simple solution, where the options are considered in decreasing order of reservation value,i.e., the value that reduces to zero the expected marginal gain for opening the box. We study for the first time this problem when order - or precedence - constraints are imposed between the boxes. We show that, despite the difficulty of defining reservation values for the boxes which take into account both in-depth and in-breath exploration of the various options, greedy optimal strategies exist and can be efficiently computed for tree-like order constraints. We also prove that finding approximately optimal adaptive search strategies is NP-hard when certain matroid constraints are used to further restrict the set of boxes which may be opened, or when the order constraints are given as reachability constraints on a DAG. We complement the above result by giving approximate adaptive search strategies based on a connection between optimal adaptive strategies and non-adaptive strategies with bounded adaptivity gap for a carefully relaxed version of the problem.


title: CoBeam: Beamforming-based Spectrum Sharing With Zero Cross-Technology Signaling for 5G Wireless Networks

abstract: This article studies an essential yet challenging problem in 5G wireless networks: \emph{Is it possible to enable spectrally-efficient spectrum sharing for heterogeneous wireless networks with different, possibly incompatible, spectrum access technologies on the same spectrum bands; without modifying the protocol stacks of existing wireless networks?} To answer this question, this article explores the system challenges that need to be addressed to enable a new spectrum sharing paradigm based on beamforming, which we refer to as {\em CoBeam}. In CoBeam, a secondary wireless network is allowed to access a spectrum band based on {\em cognitive beamforming} without mutual temporal exclusion, i.e., without interrupting the ongoing transmissions of coexisting wireless networks on the same bands; and without cross-technology communication. We first describe the main components of CoBeam, including \emph{programmable physical layer driver}, \emph{cognitive sensing engine}, and \emph{beamforming engine}, and then we showcase the potential of the CoBeam framework by designing a practical coexistence scheme between Wi-Fi and LTE on unlicensed bands. We present a prototype of the resulting coexisting Wi-Fi/U-LTE network built on off-the-shelf software radios based on which we evaluate the performance of CoBeam through an extensive experimental campaign. Performance evaluation results indicate that CoBeam can achieve on average $169\%$ throughput gain while requiring \emph{no} signaling exchange between the coexisting wireless networks.


title: The mapping class group of the Cantor tree has only geometric normal subgroups

abstract: A normal subgroup of the (extended) mapping class group of a surface is said to be geometric if its automorphism group is the mapping class group. We prove that in the case of the Cantor tree surface, every normal subgroup is geometric. We note that there is no non-trivial finite-type mapping class group for which this statement is true. We study a generalisation of the curve graph, proving that its automorphism group is again the mapping class group. This strategy is adapted from that of Brendle-Margalit and the author for certain normal subgroups in the finite-type setting.


title: Formation of complex molecules in translucent clouds: Acetaldehyde, vinyl alcohol, ketene, and ethanol via nonenergetic processing of C2H2 ice

abstract: Complex organic molecules (COMs) have been identified toward high- and low-mass protostars as well as molecular clouds, suggesting that these interstellar species originate from the early stage(s) of starformation. The reaction pathways resulting in COMs described by the formula C$_2$H$_\text{n}$O are still under debate. In this work, we investigate the laboratory possible solid-state reactions that involve simple hydrocarbons and OH-radicals along with H$_2$O ice under translucent cloud conditions (1$\leq$A$_V$$\leq$5 and \textit{n}$_\text{H}$$\sim$10$^3$ cm$^{-3}$). We focus on the interactions of C$_2$H$_2$ with H-atoms and OH-radicals, which are produced along the H$_2$O formation sequence on grain surfaces at 10 K. Ultra-high vacuum (UHV) experiments were performed to study the surface chemistry observed during C$_2$H$_2$ + O$_2$ + H codeposition, where O$_2$ was used for the in-situ generation of OH-radicals. Reflection absorption infrared spectroscopy (RAIRS) was applied to in situ monitor the initial and newly formed species. After that, a temperature-programmed desorption experiment combined with a Quadrupole mass spectrometer (TPD-QMS) was used as a complementary analytical tool. The investigated 10 K surface chemistry of C$_2$H$_2$ with H-atoms and OH-radicals not only results in semi and fully saturated hydrocarbons, such as ethylene (C$_2$H$_4$) and ethane (C$_2$H$_6$), but it also leads to the formation of COMs, such as vinyl alcohol, acetaldehyde, ketene, ethanol, and possibly acetic acid. It is concluded that OH-radical addition reactions to C$_2$H$_2$, acting as a molecular backbone, followed by isomerization (i.e., keto-enol tautomerization) via an intermolecular pathway and successive hydrogenation provides a so far experimentally unreported solid-state route for the formation of these species without the need of energetic input.


title: Color Tuning by Spontaneous Propagation of Gap Surface Plasmons in Epsilon Near Zero Nano-Cavity

abstract: This work presents numerical and experimental results of plasmonic Metal-Insulator optical nano-cavities. The systems are composed of Silver as metal, polyvinylpyrrolidone or indium tin oxide or zinc oxide as insulator. The proposed nano-cavities exhibit extraordinary effects as colors changing in function of the incident/view angles, enhancement of the experimental real and imaginary parts of the pseudo dielectric constant, an extraordinary transmission of 50$\%$ and zero reflection for both polarizations at the resonant wavelengths. These results lead to the formation and propagation of surface plasmon polaritons and their hybridization in gap surface plasmons. The concurrence of these effects allow studying the Goos-H\"anchen shift by exploiting the very narrow $\Delta$ and $\Psi$ ellipsometer parameters. Thank to their optical properties, the proposed nano-cavities could be applied in several fields such tunable color filters, optics, photonics and sensing.


title: Tridiagonalization of systems of coupled linear differential equations with variable coefficients by a Lanczos-like method

abstract: We show constructively that, under certain regularity assumptions, any system of coupled linear differential equations with variable coefficients can be tridiagonalized by a time-dependent Lanczos-like method. The proof we present formally establishes the convergence of the Lanczos-like algorithm and yields a full characterization of algorithmic breakdowns. From there, the solution of the original differential system is available in closed form. This is a key piece in evaluating the elusive ordered exponential function both formally and numerically.


title: The inconsistency of h-index: a mathematical analysis

abstract: Citation distributions are lognormal. We use 30 lognormally distributed synthetic series of numbers that simulate real series of citations to investigate the consistency of the h index. Using the lognormal cumulative distribution function, the equation that defines the h index can be formulated; this equation shows that h has a complex dependence on the number of papers (N). We also investigate the correlation between h and the number of papers exceeding various citation thresholds, from 5 to 500 citations. The best correlation is for the 100 threshold but numerous data points deviate from the general trend. The size-independent indicator h/N shows no correlation with the probability of publishing a paper exceeding any of the citation thresholds. In contrast with the h index, the total number of citations shows a high correlation with the number of papers exceeding the thresholds of 10 and 50 citations; the mean number of citations correlates with the probability of publishing a paper that exceeds any level of citations. Thus, in synthetic series, the number of citations and the mean number of citations are much better indicators of research performance than h and h/N. We discuss that in real citation distributions there are other difficulties.


title: Cross-sectional Stock Price Prediction using Deep Learning for Actual Investment Management

abstract: Stock price prediction has been an important research theme both academically and practically. Various methods to predict stock prices have been studied until now. The feature that explains the stock price by a cross-section analysis is called a "factor" in the field of finance. Many empirical studies in finance have identified which stocks having features in the cross-section relatively increase and which decrease in terms of price. Recently, stock price prediction methods using machine learning, especially deep learning, have been proposed since the relationship between these factors and stock prices is complex and non-linear. However, there are no practical examples for actual investment management. In this paper, therefore, we present a cross-sectional daily stock price prediction framework using deep learning for actual investment management. For example, we build a portfolio with information available at the time of market closing and invest at the time of market opening the next day. We perform empirical analysis in the Japanese stock market and confirm the profitability of our framework.


title: A simplified expression for the solution of cubic polynomial equations using function evaluation

abstract: This paper presents a simplified method of expressing the solution to cubic equations in terms of function evaluation only. The method eliminates the need to manipulate the original coefficients of the cubic polynomial and makes the solution free from such coefficients. In addition, the usual substitution needed to reduce the cubic equation is implicit in that the final solution is expressed in terms only of the function and derivative values of the given cubic polynomial at a single point. The proposed methodology simplifies the solution to cubic equations making them easy to remember and solve.


title: Convergent non complete interpolatory quadrature rules

abstract: We find a family of convergent schemes of nodes for non-complete interpolatory quadrature rules.


title: A sharp bound on the expected local time of a continuous ${\cal L}_2$-bounded Martingale

abstract: For a continuous ${\cal L}_2$-bounded Martingale with no intervals of constancy, starting at $0$ and having final variance $\sigma^2$, the expected local time at $x \in \cal{R}$ is at most $\sqrt{\sigma^2+x^2}-|x|$. This sharp bound is attained by Standard Brownian Motion stopped at the first exit time from the interval $(x-\sqrt{\sigma^2+x^2},x+\sqrt{\sigma^2+x^2})$. Sharp bounds for the expected maximum, maximal absolute value, maximal diameter and maximal number of upcrossings of intervals, have been established by Dubins and Schwarz (1988), Dubins, Gilat and Meilijson (2009) and by the authors (2017).


title: Convergence of End-to-End Training in Deep Unsupervised Contrastive Learning

abstract: Unsupervised contrastive learning has gained increasing attention in the latest research and has proven to be a powerful method for learning representations from unlabeled data. However, little theoretical analysis was known for this framework. In this paper, we study the optimization of deep unsupervised contrastive learning. We prove that, by applying end-to-end training that simultaneously updates two deep over-parameterized neural networks, one can find an approximate stationary solution for the non-convex contrastive loss. This result is inherently different from the existing over-parameterized analysis in the supervised setting because, in contrast to learning a specific target function, unsupervised contrastive learning tries to encode the unlabeled data distribution into the neural networks, which generally has no optimal solution. Our analysis provides theoretical insights into the practical success of these unsupervised pretraining methods.


title: ProSpect: Generating Spectral Energy Distributions with Complex Star Formation and Metallicity Histories

abstract: We introduce ProSpect, a generative galaxy spectral energy distribution (SED) package that encapsulates the best practices for SED methodologies in a number of astrophysical domains. ProSpect comes with two popular families of stellar population libraries (BC03 and EMILES), and a large variety of methods to construct star formation and metallicity histories. It models dust through the use of a Charlot & Fall attenuation model, with re-emission using Dale far-infrared templates. It also has the ability to model AGN through the inclusion of a simple AGN and hot torus model. Finally, it makes use of MAPPINGS-III photoionisation tables to produce line emission features. We test the generative and inversion utility of ProSpect through application to the Shark galaxy formation semi-analytic code, and informed by these results produce fits to the final ultraviolet to far-infrared photometric catalogues produces by the Galaxy and Mass Assembly Survey (GAMA). As part of the testing of ProSpect, we also produce a range of simple photometric stellar mass approximations covering a range of filters for both observed frame and rest frame photometry.


title: Residue-Torsion and the Laplacian on Riemannian Manifolds

abstract: The Ray-Singer analytic torsion is the zeta-function trace of a certain sum of logarithm operators on the de Rham complex. In this note we examine the residue analytic torsion, defined using the residue-trace instead of the spectral zeta function quasi-trace.


title: Periodic Magnetic Geodesics on Heisenberg Manifolds

abstract: We study the dynamics of magnetic flows on Heisenberg groups. Let $H$ denote the three-dimensional simply connected Heisenberg Lie group endowed with a left-invariant Riemannian metric and an exact, left-invariant magnetic field. Let $\Gamma$ be a lattice subgroup of $H,$ so that $\Gamma\backslash H$ is a closed nilmanifold. We first find an explicit description of magnetic geodesics on $H$, then determine all closed magnetic geodesics and their lengths for $\Gamma \backslash H$. We then consider two applications of these results: the density of periodic magnetic geodesics and marked magnetic length spectrum rigidity.


title: Flow-induced order-order transitions in amyloid fibril liquid crystalline tactoids

abstract: Understanding and controlling the director field configuration, shape, and orientation in nematic and cholesteric liquid crystals is of fundamental importance in several branches of science. Liquid crystalline droplets, also known as tactoids, which spontaneously form by nucleation and growth within the biphasic region of the phase diagram where isotropic and nematic phases coexist, challenge our current understanding of liquid crystals under confinement, due to the influence of anisotropic surface boundaries at vanishingly small interfacial tension and are mostly studied under quiescent, quasi-equilibrium conditions. Here, we show that different classes of amyloid fibril nematic and cholesteric tactoids undergo out-of-equilibrium order-order transitions by flow-induced deformations of their shape. The tactoids align under extensional flow and undergo extreme deformation into highly elongated oblate shapes, allowing the cholesteric pitch to decrease as an inverse power law of the tactoids aspect ratio. Energy functional theory and experimental measurements are combined to rationalize the critical elongation ratio above which the director-field configuration of tactoids transforms from bipolar and uniaxial cholesteric to homogenous and to debate on the thermodynamic nature of these transitions. Our findings suggest new opportunities in designing self-assembled liquid crystalline materials where structural and dynamical properties may be tuned by non-equilibrium phase transitions.


title: Model-independent determination of the nucleon charge radius from lattice QCD

abstract: Lattice QCD calculations of nucleon form factors are restricted to discrete values of the Euclidean four-momentum transfer. Therefore, the extraction of radii typically relies on parametrizing and fitting the lattice QCD data to obtain its slope close to zero momentum transfer. We investigate a new method, which allows to compute the nucleon radius directly from existing lattice QCD data, without assuming a functional form for the momentum dependence of the underlying form factor. The method is illustrated for the case of the isovector mean square charge radius of the nucleon $\langle r^2_\mathrm{isov} \rangle$ and the quark-connected contributions to $\langle r^2_p\rangle$ and $\langle r^2_n \rangle$ for the proton and neutron, respectively. Computations are performed using a single gauge ensemble with $N_f=2+1+1$ maximally twisted mass clover-improved fermions at physical quark mass and a lattice spacing of $a=0.08\mathrm{fm}$.


title: A simple proof for a $C^\infty$ Nekhoroshev theorem

abstract: We prove a $C^\infty$ version of the Nekhoroshev's estimate on the stability times of the actions in close to integrable Hamiltonian systems. The proof we give is a variant of the original Nekhoroshev's proof and it consists in first conjugating, globally in the phase space, and up to a small remainder, the system to a normal form. Then we perform the geometric part of the proof in the normalized variables. As a result, we obtain a proof which is simpler than the usual ones.


title: Geometrical aspects in the analysis of microcanonical phase-transitions

abstract: In the present work, we discuss how the functional form of thermodynamic observables can be deduced from the geometric properties of subsets of phase space. The geometric quantities taken into account are mainly extrinsic curvatures of the energy level sets of the Hamiltonian of a system under investigation. In particular, it turns out that peculiar behaviours of thermodynamic observables at a phase transition point are rooted in more fundamental changes of the geometry of the energy level sets in phase space. More specifically, we discuss how microcanonical and geometrical descriptions of phase-transitions are shaped in the special case of $\phi^4$ models with either nearest-neighbours and mean-field interactions.


title: DeepLight: Deep Lightweight Feature Interactions for Accelerating CTR Predictions in Ad Serving

abstract: Click-through rate (CTR) prediction is a crucial task in online display advertising. The embedding-based neural networks have been proposed to learn both explicit feature interactions through a shallow component and deep feature interactions using a deep neural network (DNN) component. These sophisticated models, however, slow down the prediction inference by at least hundreds of times. To address the issue of significantly increased serving delay and high memory usage for ad serving in production, this paper presents \emph{DeepLight}: a framework to accelerate the CTR predictions in three aspects: 1) accelerate the model inference via explicitly searching informative feature interactions in the shallow component; 2) prune redundant layers and parameters at intra-layer and inter-layer level in the DNN component; 3) promote the sparsity of the embedding layer to preserve the most discriminant signals. By combining the above efforts, the proposed approach accelerates the model inference by 46X on Criteo dataset and 27X on Avazu dataset without any loss on the prediction accuracy. This paves the way for successfully deploying complicated embedding-based neural networks in production for ad serving.


title: Anisotropic Fluid Cosmology: an Alternative to Dark Matter?

abstract: We use anisotropic fluid cosmology to describe the present, dark energy-dominated, universe. Similarly to what has been proposed for galactic dynamics, the anisotropic fluid gives an effective description of baryonic matter, dark energy and their possible interaction, without assuming the presence of dark matter. The resulting anisotropic fluid spacetime naturally generates inhomogeneities at small scales, triggered by an anisotropic stress, and could therefore be responsible for structure formation at these scales. Solving the cosmological equations, we show that the dynamics of the scale factor $a$ is described by usual FLRW cosmology and decouples completely from that describing inhomogeneities. We assume that the cosmological anisotropic fluid inherits the equation of state from that used to describe galaxy rotation curves. We show that, in the large scale regime, the fluid can be described as a generalized Chaplygin gas and fits well the distance modulus experimental data of type Ia supernovae, thus correctly modelling the observed accelerated expansion of the universe. Conversely, in the small scale regime, we use cosmological perturbation theory to derive the power spectrum $P(k)$ for mass density distribution. At short wavelengths, we find a $1/k^4$ behaviour, in good accordance with the observed correlation function for matter distribution at small scales.


title: Theoretical study of quantum emitters in two-dimensional silicon carbide monolayer

abstract: The features of some potential single-photon sources in two-dimensional silicon carbide monolayers is studied via ab-initio calculation and group theory analysis. A few point defects in three charge states (negative, positive and neutral) are considered. By applying performance criteria, Stone-Wales defects without and with combination of antisite defects are studied in details. The formation energy calculations reveal that neutral and positive charge states of these defects are stable. We compute the zero-phonon-line energy, the Huang-Rhys (HR) factor and the photoluminescence spectrum for the available transitions in different charge states. The calculated HR values and the related Debye-Waller factors guarantee that the Stone-Wales defects have a high potential of performing as a promising single-photon emitter.


title: Footprints of quantum pigeons

abstract: We show that in the mathematical framework of the quantum theory the classical pigeonhole principle can be violated more directly than previously suggested, i.e., in a setting closer to the traditional statement of the principle. We describe how the counterfactual reasoning of the paradox may be operationally grounded in the analysis of the tiny footprints left in the environment by the pigeons. After identifying the drawbacks of recent experiments of the quantum pigeonhole effect, we argue that a definitive experimental violation of the pigeonhole principle is still needed and propose such an implementation using modern quantum computing hardware: a superconducting circuit with transmon qubits.


title: Learning Group Structure and Disentangled Representations of Dynamical Environments

abstract: Learning disentangled representations is a key step towards effectively discovering and modelling the underlying structure of environments. In the natural sciences, physics has found great success by describing the universe in terms of symmetry preserving transformations. Inspired by this formalism, we propose a framework, built upon the theory of group representation, for learning representations of a dynamical environment structured around the transformations that generate its evolution. Experimentally, we learn the structure of explicitly symmetric environments without supervision from observational data generated by sequential interactions. We further introduce an intuitive disentanglement regularisation to ensure the interpretability of the learnt representations. We show that our method enables accurate long-horizon predictions, and demonstrate a correlation between the quality of predictions and disentanglement in the latent space.


title: Backward stochastic Volterra integral equations with jumps in a general filtration

abstract: In this paper, we study backward stochastic Volterra integral equations introduced in [26, 45] and extend the existence, uniqueness or comparison results for general filtration as in [31] (not only Brownian-Poisson setting). We also consider Lp-data and explore the time regularity of the solution in the It{\^o} setting, which is also new in this jump setting.


title: In Search for an Optimal Authenticated Byzantine Agreement

abstract: In this paper, we challenge the conventional approach of state machine replication systems to design deterministic agreement protocols in the eventually synchronous communication model. We first prove that no such protocol can guarantee bounded communication cost before the global stabilization time and propose a different approach that hopes for the best (synchrony) but prepares for the worst (asynchrony). Accordingly, we design an optimistic byzantine agreement protocol that first tries an efficient deterministic algorithm that relies on synchrony for termination only, and then, only if an agreement was not reached due to asynchrony, the protocol uses a randomized asynchronous protocol for fallback that guarantees termination with probability 1. We formally prove that our protocol achieves optimal communication complexity under all network conditions and failure scenarios. We first prove a lower bound of $\Omega(ft+ t)$ for synchronous deterministic byzantine agreement protocols, where $t$ is the failure threshold, and $f$ is the actual number of failures. Then, we present a tight upper bound and use it for the synchronous part of the optimistic protocol. Finally, for the asynchronous fallback, we use a variant of the (optimal) VABA protocol, which we reconstruct to safely combine it with the synchronous part. We believe that our adaptive to failures synchronous byzantine agreement protocol has an independent interest since it is the first protocol we are aware of which communication complexity optimally depends on the actual number of failures.


title: Separating Bohr denseness from measurable recurrence

abstract: We prove that there is a set of integers $A$ having positive upper Banach density whose difference set $A-A:=\{a-b:a,b\in A\}$ does not contain a Bohr neighborhood of any integer, answering a question asked by Bergelson, Hegyv\'ari, Ruzsa, and the author, in various combinations. In the language of dynamical systems, this result shows that there is a set of integers $S$ which is dense in the Bohr topology of $\mathbb Z$ and which is not a set of measurable recurrence. Our proof yields the following stronger result: if $S\subseteq \mathbb Z$ is dense in the Bohr topology of $\mathbb Z$, then there is a set $S'\subseteq S$ such that $S'$ is dense in the Bohr topology of $\mathbb Z$ and for all $m\in \mathbb Z,$ the set $(S'-m)\setminus \{0\}$ is not a set of measurable recurrence.


title: New structure canditates for the experimentally synthesized heptazine-based and triazine-based two dimensional graphitic carbon nitride

abstract: The widely used crystal structures for both heptazine-based and triazine-based two-dimensional (2D) graphitic carbon nitride (g-C$_3$N$_4$) are the flat P-6m2 configurations. However, the experimentally synthesized 2D g-C$_3$N$_4$ possess thickness ranging in 0.2-0.5 nm, indicating that the theoretically used flat P-6m2 configurations are not the correct ground states. In this work, we propose three new corrugated structures P321, P3m1 and Pca21 with energies of 66 (86), 77 (87) and 78 (89) meV/atom lower than that of the corresponding heptazine-based (triazine-based) g-C$_3$N$_4$ in flat P-6m2 configuration, respectively. These corrugated structures have very similar periodic patterns to the flat P-6m2 ones and they are difficult to be distinguished from each other according to their top-views. The optimized thicknesses of the three corrugated structures ranging in 1.347-3.142 {\AA} are in good agreement with the experimental results. The first-principles results show that these corrugated structural candidates are also semiconductors with band gaps slightly larger than those of the correspondingly flat P-6m2 ones. Furthermore, they possess also suitable band edge positions for sun-light-driven water-splitting at both $pH=0$ and $pH=7$ environments. Our results show that these three new structures are more promising candidates for the experimentally synthesized g-C$_3$N$_4$.


title: On the isoperimetric inequality of Coulhon and Saloff-Coste

abstract: We improve the isoperimetric inequality of Coulhon and Saloff-Coste following a method of mass transport proposed by Gromov.


title: Hypergraph Isomorphism for Groups with Restricted Composition Factors

abstract: We consider the isomorphism problem for hypergraphs taking as input two hypergraphs over the same set of vertices $V$ and a permutation group $\Gamma$ over domain $V$, and asking whether there is a permutation $\gamma \in \Gamma$ that proves the two hypergraphs to be isomorphic. We show that for input groups, all of whose composition factors are isomorphic to a subgroup of the symmetric group on $d$ points, this problem can be solved in time $(n+m)^{O((\log d)^{c})}$ for some absolute constant $c$ where $n$ denotes the number of vertices and $m$ the number of hyperedges. In particular, this gives the currently fastest isomorphism test for hypergraphs in general. The previous best algorithm for this problem due to Schweitzer and Wiebking (STOC 2019) runs in time $n^{O(d)}m^{O(1)}$. As an application of this result, we obtain, for example, an algorithm testing isomorphism of graphs excluding $K_{3,h}$ ($h \geq 3$) as a minor in time $n^{O((\log h)^{c})}$. In particular, this gives an isomorphism test for graphs of Euler genus at most $g$ running in time $n^{O((\log g)^{c})}$.


title: RapidLayout: Fast Hard Block Placement of FPGA-optimized Systolic Arrays using Evolutionary Algorithms

abstract: Evolutionary algorithms can outperform conventional placement algorithms such as simulated annealing, analytical placement as well as manual placement on metrics such as runtime, wirelength, pipelining cost, and clock frequency when mapping FPGA hard block intensive designs such as systolic arrays on Xilinx UltraScale+ FPGAs. For certain hard-block intensive, systolic array accelerator designs, the commercial-grade Xilinx Vivado CAD tool is unable to provide a legal routing solution without tedious manual placement constraints. Instead, we formulate an automatic FPGA placement algorithm for these hard blocks as a multi-objective optimization problem that targets wirelength squared and maximum bounding box size metrics. We build an end-to-end placement and routing flow called RapidLayout using the Xilinx RapidWright framework. RapidLayout runs 5-6$\times$ faster than Vivado with manual constraints and eliminates the weeks-long effort to generate placement constraints manually for the hard blocks. We also perform automated post-placement pipelining of the long wires inside each convolution block to target 650MHz URAM-limited operation. RapidLayout outperforms (1) the simulated annealer in VPR by 33% in runtime, 1.9-2.4$\times$ in wirelength, and 3-4$\times$ in bounding box size, while also (2) beating the analytical placer UTPlaceF by 9.3$\times$ in runtime, 1.8-2.2$\times$ in wirelength, and 2-2.7$\times$ in bounding box size. We employ transfer learning from a base FPGA device to speed-up placement optimization for similar FPGA devices in the UltraScale+ family by 11-14$\times$ than learning the placements from scratch.


title: Approximation of a Cauchy-Jensen additive mapping in various normed spaces

abstract: In this paper, using the fixed point and direct methods, we prove the generalized Hyers-Ulam-Rassias stability of a Cauchy-Jensen additive functional equation in various normed spaces. The concept of Hyers-Ulam-Rassias stability originated from Th. M. Rassias stability theorem that appeared in his paper: On the stability of the linear mapping in Banach spaces, Proc. Amer. Math. Soc. 72 (1978), 297-300.


title: Verifiable type-II seesaw and dark matter in a gauged $U(1)_{B-L}$ model

abstract: We propose a gauged $U(1)_{B-L}$ extension of the standard model (SM) to explain simultaneously the light neutrino masses and dark matter (DM). The generation of neutrino masses occurs through a variant of type-II seesaw mechanism in which one of the scalar triplets lies at the TeV scale yet have a large dilepton coupling, which paves a path for probing this model at colliders. The gauging of $U(1)_{B-L}$ symmetry in a type-II seesaw framework introduces $B-L$ anomalies. Therefore we invoke three right handed neutrinos $\nu_{R_{i}}$(i=1,2,3) with $B-L$ charges -4,-4,+5 to cancel the anomalies. We further show that the lightest one among the three right handed neutrinos can be a viable DM candidate. The stability of DM can be owed to a remnant $Z_2$ symmetry under which the right handed neutrinos are odd while all other particles are even. We then discuss the constraints on the model parameters from observed DM abundance and the search at direct detection experiments.


title: On admissible singular drifts of symmetric $\alpha$-stable process

abstract: We consider the problem of existence of a (unique) weak solution to the SDE describing symmetric $\alpha$-stable process with a locally unbounded drift $b:\mathbb R^d \rightarrow \mathbb R^d$, $d \geq 3$, $1<\alpha<2$. In this paper, $b$ belongs to the class of weakly form-bounded vector fields. The latter arises as the class providing the $L^2$ theory of the non-local operator behind the SDE, i.e.\,$(-\Delta)^{\frac{\alpha}{2}} + b \cdot \nabla$, and contains as proper sub-classes the other classes of singular vector fields studied in the literature in connection with this operator, such as the Kato class, weak $L^{\frac{d}{\alpha-1}}$ class and the Campanato-Morrey class (thus, $b$ can be so singular that it destroys the standard heat kernel estimates in terms of the heat kernel of the fractional Laplacian). We show that for such $b$ the operator $-(-\Delta)^{\frac{\alpha}{2}} - b \cdot \nabla$ admits a realization as a Feller generator, and that the probability measures determined by the Feller semigroup (uniquely in appropriate sense) admit description as weak solutions to the corresponding SDE. The proof is based on detailed regularity theory of $(-\Delta)^{\frac{\alpha}{2}} + b \cdot \nabla$ in $L^p$, $p>d-\alpha+1$.


title: Jittering Samples using a kd-Tree Stratification

abstract: Monte Carlo sampling techniques are used to estimate high-dimensional integrals that model the physics of light transport in virtual scenes for computer graphics applications. These methods rely on the law of large numbers to estimate expectations via simulation, typically resulting in slow convergence. Their errors usually manifest as undesirable grain in the pictures generated by image synthesis algorithms. It is well known that these errors diminish when the samples are chosen appropriately. A well known technique for reducing error operates by subdividing the integration domain, estimating integrals in each \emph{stratum} and aggregating these values into a stratified sampling estimate. Na\"{i}ve methods for stratification, based on a lattice (grid) are known to improve the convergence rate of Monte Carlo, but require samples that grow exponentially with the dimensionality of the domain. We propose a simple stratification scheme for $d$ dimensional hypercubes using the kd-tree data structure. Our scheme enables the generation of an arbitrary number of equal volume partitions of the rectangular domain, and $n$ samples can be generated in $O(n)$ time. Since we do not always need to explicitly build a kd-tree, we provide a simple procedure that allows the sample set to be drawn fully in parallel without any precomputation or storage, speeding up sampling to $O(\log n)$ time per sample when executed on $n$ cores. If the tree is implicitly precomputed ($O(n)$ storage) the parallelised run time reduces to $O(1)$ on $n$ cores. In addition to these benefits, we provide an upper bound on the worst case star-discrepancy for $n$ samples matching that of lattice-based sampling strategies, which occur as a special case of our proposed method. We use a number of quantitative and qualitative tests to compare our method against state of the art samplers for image synthesis.


title: A Newton Frank-Wolfe Method for Constrained Self-Concordant Minimization

abstract: We demonstrate how to scalably solve a class of constrained self-concordant minimization problems using linear minimization oracles (LMO) over the constraint set. We prove that the number of LMO calls of our method is nearly the same as that of the Frank-Wolfe method in the L-smooth case. Specifically, our Newton Frank-Wolfe method uses $\mathcal{O}(\epsilon^{-\nu})$ LMO's, where $\epsilon$ is the desired accuracy and $\nu:= 1 + o(1)$. In addition, we demonstrate how our algorithm can exploit the improved variants of the LMO-based schemes, including away-steps, to attain linear convergence rates. We also provide numerical evidence with portfolio design with the competitive ratio, D-optimal experimental design, and logistic regression with the elastic net where Newton Frank-Wolfe outperforms the state-of-the-art.


title: Seeds of Life in Space (SOLIS).VII. Discovery of a cold dense methanol blob toward the L1521F VeLLO system

abstract: The SOLIS (Seeds Of Life In Space) IRAM/NOEMA Large Program aims at studying a set of crucial complex organic molecules in a sample of sources, with well-known physical structure, covering the various phases of Solar-type star formation. One representative object of the transition from the prestellar core to the protostar phases has been observed toward the Very Low Luminosity Object (VeLLO) called L1521F. This type of source is important to study to make the link between prestellar cores and Class 0 sources and also to constrain the chemical evolution during the process of star formation. Two frequency windows (81.6-82.6 GHz and 96.65-97.65 GHz) were used to observe the emission from several complex organics toward the L1521F VeLLO. Only 2 transitions of methanol (A+, E2) have been detected in the narrow window centered at 96.7 GHz (with an upper limit on E1) in a very compact emission blob (~7'' corresponding to ~1000au) toward the NE of the L1521F protostar. The CS 2-1 transition is also detected within the WideX bandwidth. Consistently, with what has been found in prestellar cores, the methanol emission appears ~1000au away from the dust peak. The location of the methanol blob coincides with one of the filaments previously reported in the literature. The Tex of the gas inferred from methanol is (10$\pm$2) K, while the H2 gas density (estimated from the detected CS 2-1 emission and previous CS 5-4 ALMA obs.) is a factor >25 higher than the density in the surrounding environment (n(H2) >10$^{7}$ cm$^{-3}$). From its compactness, low excitation temperature and high gas density, we suggest that the methanol emission detected with NOEMA is either a cold and dense shock-induced blob, recently formed ($\leq$ few hundred years) by infalling gas or a cold and dense fragment that may have just been formed as a result of the intense gas dynamics found within the L1521F VeLLO system.


title: Verification benchmarks for single-phase flow in three-dimensional fractured porous media

abstract: Flow in fractured porous media occurs in the earth's subsurface, in biological tissues, and in man-made materials. Fractures have a dominating influence on flow processes, and the last decade has seen an extensive development of models and numerical methods that explicitly account for their presence. To support these developments, we present a portfolio of four benchmark cases for single-phase flow in three-dimensional fractured porous media. The cases are specifically designed to test the methods' capabilities in handling various complexities common to the geometrical structures of fracture networks. Based on an open call for participation, results obtained with 17 numerical methods were collected. This paper presents the underlying mathematical model, an overview of the features of the participating numerical methods, and their performance in solving the benchmark cases.


title: Tracer diffusion in crowded solutions of sticky polymers

abstract: Macromolecular diffusion in strongly confined geometries and crowded environments is still to a large extent an open subject in soft matter physics and biology. In this paper, we employ large-scale Langevin dynamics simulations to investigate how the diffusion of a tracer is influenced by the combined action of excluded-volume and weak attractive crowder-tracer interactions. We consider two species of tracers, standard hard-core particles described by the Weeks-Chandler-Andersen (WCA) repulsive potential and core-softened (CS) particles, which model, e.g., globular proteins, charged colloids and nanoparticles covered by polymeric brushes. These systems are characterized by the presence of two length scales in the interaction and can show water-like anomalies in their diffusion, stemming from the inherent competition between different length scales. Here we report a comprehensive study of both diffusion and structure of these two tracer species in an environment crowded by quenched configurations of polymers at increasing density. We analyze in detail how the tracer-polymer affinity and the system density affect transport as compared to the emergence of specific static spatial correlations. In particular, we find that, while hardly any differences emerge in the diffusion properties of WCA and CS particles, the propensity to develop structural order for large crowding is strongly frustrated for CS particles. Surprisingly, for large enough affinity for the crowding matrix, the diffusion coefficient of WCA tracers display a non-monotonic trend as their density is increased when compared to the zero affinity scenario. This water-like anomaly turns out to be even larger than what observed for CS particle and appears to be rooted in a similar competition between excluded-volume and affinity effects.


title: The Synthesizability of Molecules Proposed by Generative Models

abstract: The discovery of functional molecules is an expensive and time-consuming process, exemplified by the rising costs of small molecule therapeutic discovery. One class of techniques of growing interest for early-stage drug discovery is de novo molecular generation and optimization, catalyzed by the development of new deep learning approaches. These techniques can suggest novel molecular structures intended to maximize a multi-objective function, e.g., suitability as a therapeutic against a particular target, without relying on brute-force exploration of a chemical space. However, the utility of these approaches is stymied by ignorance of synthesizability. To highlight the severity of this issue, we use a data-driven computer-aided synthesis planning program to quantify how often molecules proposed by state-of-the-art generative models cannot be readily synthesized. Our analysis demonstrates that there are several tasks for which these models generate unrealistic molecular structures despite performing well on popular quantitative benchmarks. Synthetic complexity heuristics can successfully bias generation toward synthetically-tractable chemical space, although doing so necessarily detracts from the primary objective. This analysis suggests that to improve the utility of these models in real discovery workflows, new algorithm development is warranted.


title: Effects of initial state fluctuations on the mean transverse momentum

abstract: We show that in ideal hydrodynamic simulations of heavy-ion collisions, initial state fluctuations result in an increase of the mean transverse momentum of outgoing hadrons, $\langle p_t\rangle$. Specifically, $\langle p_t\rangle$ is larger by a few percent if realistic fluctuations are implemented than with smooth initial conditions, the multiplicity and impact parameter being kept fixed. We show that result can be traced back to the fact that for a given total entropy, the initial energy contained in the fluid is larger if the density profile is bumpy. We discuss the implication of these results for the extraction of the equation of state of QCD from experimental data.


title: Restructuring and Breakup of Nanowires with the Diamond Cubic Crystal Structure into Nanoparticles

abstract: A kinetic Monte Carlo approach is applied to study physical mechanisms responsible for the breakup of nanowires with the diamond cubic crystal structure into a chain of nanoparticles discovered in preceding experiments on Silicon nanowires. We show that this process is based on the well-known mechanism of roughening transition, which specifically manifests itself in quasi-one-dimensional systems/nanowires with a pronounced anisotropy of the surface energy density. Depending on the temperature and orientation of the nanowire relative to its internal crystal structure, the wavelengths of substantial cross-sectional modulations exceed its initial radius by 4 to 18 times. For certain orientations, a straight nanowire at the initial stage of evolution forms a serpentine/helical structure. The scenarios of the stage of nanowire ruptures into single nanoclusters are also diverse: either each spindle-shaped region of the nanowire transforms into a separate drop (by long-wave surface perturbations), or the adjacent short-scale beads absorb each other due to the Ostwald ripening effect, which can be accompanied by the formation of long-lived many-body dumbbells. The discovered features of the dynamics of quasi-onedimensional systems expand our conceptions of the physical mechanisms involved in the breakup of nanowires (presented by Nichols and Mullins as a classical model for such instabilities) which could be useful in applications based on chains of ordered nanoparticles.


title: On cohesive almost zero-dimensional spaces

abstract: We investigate C-sets in almost zero-dimensional spaces, showing that closed $\sigma$C-sets are C-sets. As corollaries, we prove that every rim-$\sigma$-compact almost zero-dimensional space is zero-dimensional and that each cohesive almost zero-dimensional space is nowhere rational. To show these results are sharp, we construct a rim-discrete connected set with an explosion point. We also show every cohesive almost zero-dimensional subspace of $($Cantor set$)$$\times\mathbb R$ is nowhere dense.


title: New Regularity Criteria for Weak Solutions to the MHD Equations in Terms of an Associated Pressure

abstract: We prove that if $0<T_0<T\leq\infty$, $(\mathbf{u},\mathbf{b},p)$ is a suitable weak solution of the MHD equations in $\mathbb{R}^3\times(0,T)$ and either $\mathcal{F}_{\gamma}(p_-)\in L^{\infty}(0,T_0;\, L^{3/2}(\mathbb{R}^3))$ or $\mathcal{F}_{\gamma}((|\mathbf{u}|^2+ |\mathbf{b}|^2+2p)_+)\in L^{\infty}(0,T_0;\, L^{3/2}(\mathbb{R}^3))$ for some $\gamma>0$, where $\mathcal{F}_{\gamma}(s)=s\, [\ln{}(1+ s)]^{1+\gamma}$ and the subscripts "$-$" and "$+$" denote the negative and the nonnegative part, respectively, then the solution $(\mathbf{u},\mathbf{b},p)$ has no singular points in $\mathbb{R}^3\times(0,T_0]$.


title: An Experimental Study of ILP Formulations for the Longest Induced Path Problem

abstract: Given a graph $G=(V,E)$, the longest induced path problem asks for a maximum cardinality node subset $W\subseteq V$ such that the graph induced by $W$ is a path. It is a long established problem with applications, e.g., in network analysis. We propose novel integer linear programming (ILP) formulations for the problem and discuss efficient implementations thereof. Comparing them with known formulations from literature, we prove that they are beneficial in theory, yielding stronger relaxations. Moreover, our experiments show their practical superiority.


title: Ductile and brittle crack-tip response in equimolar refractory high-entropy alloys

abstract: Understanding the strengthening and deformation mechanisms in refractory high-entropy alloys (HEAs), proposed as new high-temperature material, is required for improving their typically insufficient room-temperature ductility. Here, density-functional theory simulations and a continuum mechanics analysis were conducted to systematically investigate the competition between cleavage decohesion and dislocation emission from a crack tip in the body-centered cubic refractory HEAs HfNbTiZr, MoNbTaVW, MoNbTaW, MoNbTiV, and NbTiVZr. This crack-tip competition is evaluated for tensile loading and a totality of 15 crack configurations and slip systems. Our results predict that dislocation plasticity at the crack tip is generally unfavorable -- although the competition is close for some crack orientations, suggesting intrinsic brittleness and low crack-tip fracture toughness in these five HEAs at zero temperature. Fluctuations in local alloy composition, investigated for HfNbTiZr, can locally reduce the resistance to dislocation emission for a slip system relative to the configuration average of that slip system, but do not change the dominant crack-tip response. In the case of single-crystal MoNbTaW, where an experimental, room-temperature fracture-toughness value is available for a crack on a \{100\} plane, theoretical and experimental results agree favorably. Factors that may limit the agreement are discussed. We survey the effect of material anisotropy on preferred crack tip orientations, which are found to be alloy specific. Mixed-mode loadings are found to shift the competition in favor of cleavage or dislocation nucleation, depending on crack configuration and amplified by the effect of material anisotropy on crack tip stresses.


title: Two new identities involving the Bernoulli numbers, the Euler numbers, and the Stirling numbers of the first kind

abstract: We derive two new identities involving the Bernoulli numbers, the Euler numbers, and the Stirling numbers of the first kind using analytic continuation of a well known identity for the Stirling numbers of the first kind.


title: Eigenvalue gaps for hyperbolic groups and semigroups

abstract: Given a locally constant linear cocycle over a subshift of finite type, we show that the existence of a uniform gap between the i-th and (i+1)-th Lyapunov exponents for all invariant measures implies the existence of a dominated splitting of index i. We establish a similar result for sofic subshifts coming from word hyperbolic groups, in relation with Anosov representations of such groups. We discuss the case of finitely generated semigroups, and propose a notion of Anosov representation in this setting.


title: Meta-learning Extractors for Music Source Separation

abstract: We propose a hierarchical meta-learning-inspired model for music source separation (Meta-TasNet) in which a generator model is used to predict the weights of individual extractor models. This enables efficient parameter-sharing, while still allowing for instrument-specific parameterization. Meta-TasNet is shown to be more effective than the models trained independently or in a multi-task setting, and achieve performance comparable with state-of-the-art methods. In comparison to the latter, our extractors contain fewer parameters and have faster run-time performance. We discuss important architectural considerations, and explore the costs and benefits of this approach.


title: Learning Robust Representations via Multi-View Information Bottleneck

abstract: The information bottleneck principle provides an information-theoretic method for representation learning, by training an encoder to retain all information which is relevant for predicting the label while minimizing the amount of other, excess information in the representation. The original formulation, however, requires labeled data to identify the superfluous information. In this work, we extend this ability to the multi-view unsupervised setting, where two views of the same underlying entity are provided but the label is unknown. This enables us to identify superfluous information as that not shared by both views. A theoretical analysis leads to the definition of a new multi-view model that produces state-of-the-art results on the Sketchy dataset and label-limited versions of the MIR-Flickr dataset. We also extend our theory to the single-view setting by taking advantage of standard data augmentation techniques, empirically showing better generalization capabilities when compared to common unsupervised approaches for representation learning.


title: Dimension reduction through Gamma convergence for general prestrained thin elastic sheets

abstract: We study thin films with residual strain by analyzing the $\Gamma-$limit of non-Euclidean elastic energy functionals as the material's thickness tends to $0.$ We begin by extending prior results \cite{bhattacharya2016plates} \cite{agostiniani2018heterogeneous} \cite{lewicka2018dimension} \cite{schmidt2007plate}, to a wider class of films, whose prestrain depends on both the midplate and the transversal variables. The ansatz for our $\Gamma-$convergence result uses a specific type of wrinkling, which is built on exotic solutions to the Monge-Ampere equation, constructed via convex integration \cite{lewicka2017convex}. We show that the expression for our $\Gamma-$limit has a natural interpretation in terms of the orthogonal projection of the residual strain onto a suitable subspace. We also show that some type of wrinkling phenomenon is necessary to match the lower bound of the $\Gamma-$limit in certain circumstances. These results all assume a prestrain of the same order as the thickness; we also discuss why it is natural to focus on that regime by considering what can happen when the prestrain is larger.


title: Learning to Prove Theorems by Learning to Generate Theorems

abstract: We consider the task of automated theorem proving, a key AI task. Deep learning has shown promise for training theorem provers, but there are limited human-written theorems and proofs available for supervised learning. To address this limitation, we propose to learn a neural generator that automatically synthesizes theorems and proofs for the purpose of training a theorem prover. Experiments on real-world tasks demonstrate that synthetic data from our approach improves the theorem prover and advances the state of the art of automated theorem proving in Metamath. Code is available at https://github.com/princeton-vl/MetaGen.


title: Fractional Types: Expressive and Safe Space Management for Ancilla Bits

abstract: In reversible computing, the management of space is subject to two broad classes of constraints. First, as with general-purpose computation, every allocation must be paired with a matching de-allocation. Second, space can only be safely de-allocated if its contents are restored to their initial value from allocation time. Generally speaking, the state of the art provides limited partial solutions that address the first constraint by imposing a stack discipline and by leaving the second constraint to programmers' assertions. We propose a novel approach based on the idea of fractional types. As a simple intuitive example, allocation of a new boolean value initialized to $\texttt{false}$ also creates a value $1/{\texttt{false}}$ that can be thought of as a garbage collection (GC) process specialized to reclaim, and only reclaim, storage containing the value $\texttt{false}$. This GC process is a first-class entity that can be manipulated, decomposed into smaller processes and combined with other GC processes. We formalize this idea in the context of a reversible language founded on type isomorphisms, prove its fundamental correctness properties, and illustrate its expressiveness using a wide variety of examples. The development is backed by a fully-formalized Agda implementation.


title: An Efficient Robust Approach to the Day-ahead Operation of an Aggregator of Electric Vehicles

abstract: The growing use of electric vehicles (EVs) may hinder their integration into the electricity system as well as their efficient operation due to the intrinsic stochasticity associated with their driving patterns. In this work, we assume a profit-maximizer EV-aggregator who participates in the day-ahead electricity market. The aggregator accounts for the technical aspects of each individual EV and the uncertainty in its driving patterns. We propose a hierarchical optimization approach to represent the decision-making of this aggregator. The upper level models the profit-maximizer aggregator's decisions on the EV-fleet operation, while a series of lower-level problems computes the worst-case EV availability profiles in terms of battery draining and energy exchange with the market. Then, this problem can be equivalently transformed into a mixed-integer linear single-level equivalent given the totally unimodular character of the constraint matrices of the lower-level problems and their convexity. Finally, we thoroughly analyze the benefits of the hierarchical model compared to the results from stochastic and deterministic models.


title: Tidally induced warps of spiral galaxies in IllustrisTNG

abstract: Warps are common features in both stellar and gaseous disks of nearby spiral galaxies with the latter usually easier to detect. Several theories have been proposed in the literature to explain their formation and prevalence, including tidal interactions with external galaxies. Observational correlations also suggest the importance of tides for warp formation. Here, we use the TNG100 run from the magnetohydrodynamical cosmological simulation suite IllustrisTNG to investigate the connection between interactions and the formation of gas warps. We find that in the sample of well-resolved gas-rich spiral galaxies ($10^{10}\lesssim\mathrm{M_{*}/M_{\odot}}\lesssim10^{11}$ at $z=0$) from the simulation TNG100-1, about $16\%$ possess the characteristic S-shaped warp. Around one third of these objects have their vertical morphology induced by interactions with other galaxies. Half of these interactions end with the perturber absorbed by the host. Warps formed in interactions are more asymmetrical than the remaining sample, however after the interaction the asymmetry decreases with time. We find that warps induced by interactions survive on average for $<1$ Gyr. The angle between the orbital angular momentum of the perturber and the angular momentum of the host's disk that most likely leads to warp formation is around 45 degrees. While our main goal is to investigate tidally induced warps, we find that during interactions in addition to tides, new gas that is accreted from infalling satellites also can contribute to warp formation.


title: Provenance of classical Hamiltonian time crystals

abstract: Classical Hamiltonian systems with conserved charges and those with constraints often describe dynamics on a pre-symplectic manifold. Here we show that a pre-symplectic manifold is also the proper stage to describe autonomous energy conserving Hamiltonian time crystals. We explain how the occurrence of a time crystal relates to the wider concept of spontaneously broken symmetries; in the case of a time crystal, the symmetry breaking takes place in a dynamical context. We then analyze in detail two examples of time crystalline Hamiltonian dynamics. The first example is a piecewise linear closed string, with dynamics determined by a Lie-Poisson bracket and Hamiltonian that relates to membrane stability. We explain how the Lie-Poisson brackets descents to a time crystalline pre-symplectic bracket, and we show that the Hamiltonian dynamics supports two phases; in one phase we have a time crystal and in the other phase time crystals are absent. The second example is a discrete Hamiltonian variant of the Q-ball Lagrangian of time dependent non-topological solitons. We explain how a Q-ball becomes a time crystal, and we construct examples of time crystalline Q-balls.

categories: hep-th cond-mat.other math-ph math.MP nlin.PS physics.class-ph


title: Gaming Helps! Learning from Strategic Interactions in Natural Dynamics

abstract: We consider an online regression setting in which individuals adapt to the regression model: arriving individuals are aware of the current model, and invest strategically in modifying their own features so as to improve the predicted score that the current model assigns to them. Such feature manipulation has been observed in various scenarios -- from credit assessment to school admissions -- posing a challenge for the learner. Surprisingly, we find that such strategic manipulations may in fact help the learner recover the meaningful variables -- that is, the features that, when changed, affect the true label (as opposed to non-meaningful features that have no effect). We show that even simple behavior on the learner's part allows her to simultaneously i) accurately recover the meaningful features, and ii) incentivize agents to invest in these meaningful features, providing incentives for improvement.


title: Secure-by-synthesis network with active deception and temporal logic specifications

abstract: This paper is concerned with the synthesis of strategies in network systems with active cyber deception. Active deception in a network employs decoy systems and other defenses to conduct defensive planning against the intrusion of malicious attackers who have been confirmed by sensing systems. In this setting, the defender's objective is to ensure the satisfaction of security properties specified in temporal logic formulas. We formulate the problem of deceptive planning with decoy systems and other defenses as a two-player games with asymmetrical information and Boolean payoffs in temporal logic. We use level-2 hypergame with temporal logic objectives to capture the incomplete/incorrect knowledge of the attacker about the network system as a payoff misperception. The true payoff function is private information of the defender. Then, we extend the solution concepts of $omega$-regular games to analyze the attacker's rational strategy given her incomplete information. By generalizing the solution of level-2 hypergame in the normal form to extensive form, we extend the solutions of games with safe temporal logic objectives to decide whether the defender can ensure security properties to be satisfied with probability one, given any possible strategy that is perceived to be rational by the attacker. Further, we use the solution of games with co-safe (reachability) temporal logic objectives to determine whether the defender can engage the attacker, by directing the attacker to a high-fidelity honeypot. The effectiveness of the proposed synthesis methods is illustrated with synthetic network systems with honeypots.


title: On some families of Gushel-Mukai fourfolds

abstract: We give explicit descriptions of some Noether-Lefschetz divisors in the moduli space of Gushel-Mukai fourfolds. As a consequence we obtain that their Kodaira dimension is negative.


title: A Case for Data Centre Traffic Management on Software Programmable Ethernet Switches

abstract: Virtualisation first and cloud computing later has led to a consolidation of workload in data centres that also comprises latency-sensitive application domains such as High Performance Computing and telecommunication. These types of applications require strict latency guarantees to maintain their Quality of Service. In virtualised environments with their churn, this demands for adaptability and flexibility to satisfy. At the same time, the mere scale of the infrastructures favours commodity (Ethernet) over specialised (Infiniband) hardware. For that purpose, this paper introduces a novel traffic management algorithm that combines Rate-limited Strict Priority and Deficit round-robin for latency-aware and fair scheduling respectively. In addition, we present an implementation of this algorithm on the bmv2 P4 software switch by evaluating it against standard priority-based and best-effort scheduling.


title: Low-Rank Bottleneck in Multi-head Attention Models

abstract: Attention based Transformer architecture has enabled significant advances in the field of natural language processing. In addition to new pre-training techniques, recent improvements crucially rely on working with a relatively larger embedding dimension for tokens. Unfortunately, this leads to models that are prohibitively large to be employed in the downstream tasks. In this paper we identify one of the important factors contributing to the large embedding size requirement. In particular, our analysis highlights that the scaling between the number of heads and the size of each head in the current architecture gives rise to a low-rank bottleneck in attention heads, causing this limitation. We further validate this in our experiments. As a solution we propose to set the head size of an attention unit to input sequence length, and independent of the number of heads, resulting in multi-head attention layers with provably more expressive power. We empirically show that this allows us to train models with a relatively smaller embedding dimension and with better performance scaling.


title: Models of communication and control for brain networks: distinctions, convergence, and future outlook

abstract: Recent advances in computational models of signal propagation and routing in the human brain have underscored the critical role of white matter structure. A complementary approach has utilized the framework of network control theory to better understand how white matter constrains the manner in which a region or set of regions can direct or control the activity of other regions. Despite the potential for both of these approaches to enhance our understanding of the role of network structure in brain function, little work has sought to understand the relations between them. Here, we seek to explicitly bridge computational models of communication and principles of network control in a conceptual review of the current literature. By drawing comparisons between communication and control models in terms of the level of abstraction, the dynamical complexity, the dependence on network attributes, and the interplay of multiple spatiotemporal scales, we highlight the convergence of and distinctions between the two frameworks. Based on the understanding of the intertwined nature of communication and control in human brain networks, this work provides an integrative perspective for the field and outlines exciting directions for future work.


title: Long-lived entanglement generation of nuclear spins using coherent light

abstract: Nuclear spins of noble-gas atoms are exceptionally isolated from the environment and can maintain their quantum properties for hours at room temperature. Here we develop a mechanism for entangling two such distant macroscopic ensembles by using coherent light input. The interaction between the light and the noble-gas spins in each ensemble is mediated by spin-exchange collisions with alkali-metal spins, which are only virtually excited. The relevant conditions for experimental realizations with ^{3}\text{He} or ^{129}\text{Xe} are outlined.


title: Regularizing Semi-supervised Graph Convolutional Networks with a Manifold Smoothness Loss

abstract: Existing graph convolutional networks focus on the neighborhood aggregation scheme. When applied to semi-supervised learning, they often suffer from the overfitting problem as the networks are trained with the cross-entropy loss on a small potion of labeled data. In this paper, we propose an unsupervised manifold smoothness loss defined with respect to the graph structure, which can be added to the loss function as a regularization. We draw connections between the proposed loss with an iterative diffusion process, and show that minimizing the loss is equivalent to aggregate neighbor predictions with infinite layers. We conduct experiments on multi-layer perceptron and existing graph networks, and demonstrate that adding the proposed loss can improve the performance consistently.


title: Fully convolutional networks for structural health monitoring through multivariate time series classification

abstract: We propose a novel approach to Structural Health Monitoring (SHM), aiming at the automatic identification of damage-sensitive features from data acquired through pervasive sensor systems. Damage detection and localization are formulated as classification problems, and tackled through Fully Convolutional Networks (FCNs). A supervised training of the proposed network architecture is performed on data extracted from numerical simulations of a physics-based model (playing the role of digital twin of the structure to be monitored) accounting for different damage scenarios. By relying on this simplified model of the structure, several load conditions are considered during the training phase of the FCN, whose architecture has been designed to deal with time series of different length. The training of the neural network is done before the monitoring system starts operating, thus enabling a real time damage classification. The numerical performances of the proposed strategy are assessed on a numerical benchmark case consisting of an eight-story shear building subjected to two load types, one of which modeling random vibrations due to low-energy seismicity. Measurement noise has been added to the responses of the structure to mimic the outputs of a real monitoring system. Extremely good classification capacities are shown: among the nine possible alternatives (represented by the healthy state and by a damage at any floor), damage is correctly classified in up to 95% of cases, thus showing the strong potential of the proposed approach in view of the application to real-life cases.


title: Towards an Appropriate Query, Key, and Value Computation for Knowledge Tracing

abstract: Knowledge tracing, the act of modeling a student's knowledge through learning activities, is an extensively studied problem in the field of computer-aided education. Although models with attention mechanism have outperformed traditional approaches such as Bayesian knowledge tracing and collaborative filtering, they share two limitations. Firstly, the models rely on shallow attention layers and fail to capture complex relations among exercises and responses over time. Secondly, different combinations of queries, keys and values for the self-attention layer for knowledge tracing were not extensively explored. Usual practice of using exercises and interactions (exercise-response pairs) as queries and keys/values respectively lacks empirical support. In this paper, we propose a novel Transformer based model for knowledge tracing, SAINT: Separated Self-AttentIve Neural Knowledge Tracing. SAINT has an encoder-decoder structure where exercise and response embedding sequence separately enter the encoder and the decoder respectively, which allows to stack attention layers multiple times. To the best of our knowledge, this is the first work to suggest an encoder-decoder model for knowledge tracing that applies deep self-attentive layers to exercises and responses separately. The empirical evaluations on a large-scale knowledge tracing dataset show that SAINT achieves the state-of-the-art performance in knowledge tracing with the improvement of AUC by 1.8% compared to the current state-of-the-art models.


title: Strategic advantages in mean field games with a major player

abstract: This note is concerned with a modeling question arising from the mean field games theory. We show how to model mean field games involving a major player which has a strategic advantage, while only allowing closed loop markovian strategies for all the players. We illustrate this property through three examples.


title: Unified approach to spectral properties of multipliers

abstract: Let $\mathbb B_n$ be the open unit ball in $\mathbb C^n$. We characterize the spectra of pointwise multipliers $M_u$ acting on Banach spaces of analytic functions on $\mathbb B_n$ satisfying some general conditions. These spaces include Bergman-Sobolev spaces $A^p_{\alpha,\beta}$, Bloch-type spaces $\mathcal B_\alpha$, weighted Hardy spaces $H^p_w$ with Muckenhoupt weights and Hardy-Sobolev Hilbert spaces $H^2_\beta$. Moreover, we describe the essential spectra of multipliers in most of the aforementioned spaces, in particular, in those spaces for which the set of multipliers is a subset of the ball algebra.


title: Back-and-Forth prediction for deep tensor compression

abstract: Recent AI applications such as Collaborative Intelligence with neural networks involve transferring deep feature tensors between various computing devices. This necessitates tensor compression in order to optimize the usage of bandwidth-constrained channels between devices. In this paper we present a prediction scheme called Back-and-Forth (BaF) prediction, developed for deep feature tensors, which allows us to dramatically reduce tensor size and improve its compressibility. Our experiments with a state-of-the-art object detector demonstrate that the proposed method allows us to significantly reduce the number of bits needed for compressing feature tensors extracted from deep within the model, with negligible degradation of the detection performance and without requiring any retraining of the network weights. We achieve a 62% and 75% reduction in tensor size while keeping the loss in accuracy of the network to less than 1% and 2%, respectively.


title: Dark Energy Gravitational Wave Observations and Ice Age Periodicity

abstract: Dark Energy is the dominant component of the energy density of the Universe. However, it is also very elusive since its interaction with the rest of the Universe is primarily gravitational.Since Dark Energy is a low energy phenomenon from the perspective of particle physics and field theory, a fundamental approach based on fields in curved space is sufficient to understand the current dynamics of Dark Energy. The key issue is to understand the gravitational dynamics of Dark Energy and its observational consequences. However, finding the observational consequences of Dark Energy dynamics has been a very challenging task. For something which is the dominant component of the energy density of the Universe, Dark Energy appears to be very distant and reclusive. Here we show that the Dark Energy dynamics results in the production of gravitational waves which produce the ellipticity variation in the orbit of earth that results in the periodicity of the Ice Ages observed and documented by geologists and climatologists. Previously, no observational signature of gravitational waves produced by Dark Energy dynamics has been reported. Further, no interpretation of the ellipticity variation of the orbit of earth due to gravitational waves or the linking of such gravitational waves to the Ice Age periodicity has been reported previously. We hope that the current work will lead to some fresh insights and some more interesting work.


title: Below-threshold CP-odd Higgs boson search via A -> Z* h at the LHC

abstract: We study the process q q, g g -> A -> Z* h in a 2-Higgs Doublet Model Type-II where the mass of the CP-odd Higgs state A is lower than the rest mass of the Z and h particles (the latter being the Standard Model-like Higgs state), i.e., mA < mZ + mh = 215 GeV. This is a mass range which is not being currently tested by ATLAS and CMS in this channel, yet we show that there can be sensitivity to it already during Runs 2 and 3, assuming leptonic decays of the gauge boson and bottom-antibottom quark ones for the Higgs boson.


title: Climate Change and Grain Production Fluctuations

abstract: 50 year-long time series from a Long Term Agronomic Experiment have been used to to investigate the effects of climate change on yields of Wheat and Maize. Trends and fluctuations, useful to estimate production forecasts and related risks are compared to national ones, a classical regional climatic index as Western Mediterranean Oscillation Index, and a global one given by Sun Spot Number. Data, denoised by EMD and SSA, show how SSN oscillations slowing down in the last decades, affects regional scale dynamics, where in the last two decades a range of fluctuations (7-16 years) are also evident. Both signals reflects on yield fluctuations of Wheat and Maize both a national and local level.


title: Lake Ice Detection from Sentinel-1 SAR with Deep Learning

abstract: Lake ice, as part of the Essential Climate Variable (ECV) lakes, is an important indicator to monitor climate change and global warming. The spatio-temporal extent of lake ice cover, along with the timings of key phenological events such as freeze-up and break-up, provide important cues about the local and global climate. We present a lake ice monitoring system based on the automatic analysis of Sentinel-1 Synthetic Aperture Radar (SAR) data with a deep neural network. In previous studies that used optical satellite imagery for lake ice monitoring, frequent cloud cover was a main limiting factor, which we overcome thanks to the ability of microwave sensors to penetrate clouds and observe the lakes regardless of the weather and illumination conditions. We cast ice detection as a two class (frozen, non-frozen) semantic segmentation problem and solve it using a state-of-the-art deep convolutional network (CNN). We report results on two winters ( 2016 - 17 and 2017 - 18 ) and three alpine lakes in Switzerland. The proposed model reaches mean Intersection-over-Union (mIoU) scores >90% on average, and >84% even for the most difficult lake. Additionally, we perform cross-validation tests and show that our algorithm generalises well across unseen lakes and winters.


title: A Lightweight ISA Extension for AES and SM4

abstract: We describe a lightweight RISC-V ISA extension for AES and SM4 block ciphers. Sixteen instructions (and a subkey load) is required to implement an AES round with the extension, instead of 80 without. An SM4 step (quarter-round) has 6.5 arithmetic instructions, a similar reduction. Perhaps even more importantly the ISA extension helps to eliminate slow, secret-dependent table lookups and to protect against cache timing side-channel attacks. Having only one S-box, the extension has a minimal hardware size and is well suited for ultra-low power applications. AES and SM4 implementations using the ISA extension also have a much-reduced software footprint. The AES and SM4 instances can share the same data paths but are independent in the sense that a chip designer can implement SM4 without AES and vice versa. Full AES and SM4 assembler listings, HDL source code for instruction's combinatorial logic, and C code for emulation is provided to the community under a permissive open source license. The implementation contains depth- and size-optimized joint AES and SM4 S-Box logic based on the Boyar-Peralta construction with a shared non-linear middle layer, demonstrating additional avenues for logic optimization. The instruction logic has been experimentally integrated into the single-cycle execution path of the "Pluto" RV32 core and has been tested on an FPGA system.


title: Quantum initial conditions for inflation and canonical invariance

abstract: We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since the initial and destination systems use different times. Actions related to each other via such transformations yield identical equations of motion and preserve the commutator structure. They further preserve the time-evolution of expectation values of quantum operators unless the vacuum state also changes under the transformation. These properties suggest that it is of interest to investigate vacuum prescriptions that also remain unchanged under canonical transformations. We find that initial conditions derived via minimising the vacuum expectation value of the Hamiltonian and those obtained using the Danielsson vacuum prescription are not invariant under these transformations, whereas those obtained by minimising the local energy density are. We derive the range of physically distinct initial conditions obtainable by Hamiltonian diagonalisation, and illustrate their effect on the scalar primordial power spectrum and the Cosmic Microwave Background under the just enough inflation model. We also generalise the analogy between the dynamics of a quantum scalar field on a curved, time-dependent spacetime and the gauge-invariant curvature perturbation. We argue that the invariance of the vacuum prescription obtained by minimising the renormalised stress--energy tensor should make it the preferred procedure for setting initial conditions for primordial perturbations. All other procedures reviewed in this work yield ambiguous initial conditions, which is problematic both in theory and practice.


title: Necessary conditions for binomial collisions

abstract: We shall give some necessary conditions for the equation $\binom{x}{a}=\binom{y}{b}$ to hold: if $\binom{2n+\delta}{n-m}=\binom{2n+l}{n-k}$ with $\delta=0$ or $1$, $0<m\leq 0.735k$, $k<n$ and $n$ sufficiently large, then $l>(cn/\log n)^{40/21}$ for some constant $c$.


title: Dynamic Graph Learning based on Graph Laplacian

abstract: The purpose of this paper is to infer a global (collective) model of time-varying responses of a set of nodes as a dynamic graph, where the individual time series are respectively observed at each of the nodes. The motivation of this work lies in the search for a connectome model which properly captures brain functionality upon observing activities in different regions of the brain and possibly of individual neurons. We formulate the problem as a quadratic objective functional of observed node signals over short time intervals, subjected to the proper regularization reflecting the graph smoothness and other dynamics involving the underlying graph's Laplacian, as well as the time evolution smoothness of the underlying graph. The resulting joint optimization is solved by a continuous relaxation and an introduced novel gradient-projection scheme. We apply our algorithm to a real-world dataset comprising recorded activities of individual brain cells. The resulting model is shown to not only be viable but also efficiently computable.


title: Synthesis of Deceptive Strategies in Reachability Games with Action Misperception

abstract: We consider a class of two-player turn-based zero-sum games on graphs with reachability objectives, known as reachability games, where the objective of Player 1 (P1) is to reach a set of goal states, and that of Player 2 (P2) is to prevent this. In particular, we consider the case where the players have asymmetric information about each other's action capabilities: P2 starts with an incomplete information (misperception) about P1's action set, and updates the misperception when P1 uses an action previously unknown to P2. When P1 is made aware of P2's misperception, the key question is whether P1 can control P2's perception so as to deceive P2 into selecting actions to P1's advantage? We show that there might exist a deceptive winning strategy for P1 that ensures P1's objective is achieved with probability one from a state otherwise losing for P1, had the information being symmetric and complete. We present three key results: First, we introduce a dynamic hypergame model to capture the reachability game with evolving misperception of P2. Second, we present a fixed-point algorithm to compute the Deceptive Almost-Sure Winning (DASW) region and DASW strategy. Finally, we show that DASW strategy is at least as powerful as Almost-Sure Winning (ASW) strategy in the game in which P1 does not account for P2's misperception. We illustrate our algorithm using a robot motion planning in an adversarial environment.


title: Object-based Metamorphic Testing through Image Structuring

abstract: Testing software is often costly due to the need of mass-producing test cases and providing a test oracle for it. This is often referred to as the oracle problem. One method that has been proposed in order to alleviate the oracle problem is metamorphic testing. Metamorphic testing produces new test cases by altering an existing test case, and uses the metamorphic relation between the inputs and the outputs of the System Under Test (SUT) to predict the expected outputs of the produced test cases. Metamorphic testing has often been used for image processing software, where changes are applied to the image's attributes to create new test cases with annotations that are the same as the original image. We refer to this existing method as the image-based metamorphic testing. In this research, we propose an object-based metamorphic testing and a composite metamorphic testing which combines different metamorphic testing approaches to relatively increase test coverage.


title: Temperature distribution inside a double-cladding optical fiber laser or amplifier

abstract: The temperature distribution inside a double-cladding optical fiber laser or amplifier is examined in detail. Traditionally, the quantum defect in the core is taken to be the main source of heating in an active optical fiber. However, contributions from the parasitic absorption of the signal and the pump may also play an important role, especially for low quantum defect or radiation-balanced lasers and amplifiers. The contributions to the heating in both the core and the inner-cladding are considered and analyzed in general terms in this paper. In particular, it is shown that if the maximum tolerable surface temperature of the fiber relative to the ambient is taken to be 300 degrees Celsius to avoid damaging the fiber's outer polymer cladding, the core temperature rises only in the range of 0-5 degrees Celsius relative to the inner-cladding for an air-cooled fiber. However, for a water-cooled fiber, the core temperature can be higher than the inner-cladding by as much as 50 degrees Celsius, potentially changing a single-mode core to multimode due to the thermo-optic effect.


title: Bit Allocation for Multi-Task Collaborative Intelligence

abstract: Recent studies have shown that collaborative intelligence (CI) is a promising framework for deployment of Artificial Intelligence (AI)-based services on mobile devices. In CI, a deep neural network is split between the mobile device and the cloud. Deep features obtained at the mobile are compressed and transferred to the cloud to complete the inference. So far, the methods in the literature focused on transferring a single deep feature tensor from the mobile to the cloud. Such methods are not applicable to some recent, high-performance networks with multiple branches and skip connections. In this paper, we propose the first bit allocation method for multi-stream, multi-task CI. We first establish a model for the joint distortion of the multiple tasks as a function of the bit rates assigned to different deep feature tensors. Then, using the proposed model, we solve the rate-distortion optimization problem under a total rate constraint to obtain the best rate allocation among the tensors to be transferred. Experimental results illustrate the efficacy of the proposed scheme compared to several alternative bit allocation methods.


title: The monitoring problem for timed automata

abstract: We study a variant of the classical membership problem in automata theory, which consists of deciding whether a given input word is accepted by a given automaton. We do so under a different perspective, that is, we consider a dynamic version of the problem, called monitoring problem, where the automaton is fixed and the input is revealed as in a stream, one symbol at a time following the natural order on positions. The goal here is to design a dynamic data structure that can be queried about whether the word consisting of symbols revealed so far is accepted by the automaton, and that can be efficiently updated when the next symbol is revealed. We provide complexity bounds for this monitoring problem, by considering timed automata that process symbols interleaved with timestamps. The main contribution is that monitoring of a one-clock timed automaton, with all its components but the clock constants fixed, can be done in amortised constant time per input symbol.


title: Ultrasensitive THz Biosensor for PCR-free cDNA detection based on frequency selective surfaces

abstract: THz technologies are a powerful tool for label-free detection of biomolecules. However, significant reduction of the lower detection limit is required to apply THz-sensors in biomedical diagnosis. This paper reports an ultrasensitive THz-biosensor based on asymmetric double split ring resonators (aDSRR) for the direct label- and PCR-free detection of DNA at physiologically relevant concentrations. We introduce selective functionalization and localized electric field concentration to enhance aDSRR sensitivity and specificity. The sensor characteristics are demonstrated using the human tumor marker MIA in cDNA samples produced from total RNA without PCR-amplification. Measurements of DNA samples with concentrations as low as $1.55\,\times\,10^{-12}\,$mol/l are presented.


title: Retrain or not retrain? -- efficient pruning methods of deep CNN networks

abstract: Convolutional neural networks (CNN) play a major role in image processing tasks like image classification, object detection, semantic segmentation. Very often CNN networks have from several to hundred stacked layers with several megabytes of weights. One of the possible methods to reduce complexity and memory footprint is pruning. Pruning is a process of removing weights which connect neurons from two adjacent layers in the network. The process of finding near optimal solution with specified drop in accuracy can be more sophisticated when DL model has higher number of convolutional layers. In the paper few approaches based on retraining and no retraining are described and compared together.


title: Nearest $\Omega$-stable matrix via Riemannian optimization

abstract: We study the problem of finding the nearest $\Omega$-stable matrix to a certain matrix $A$, i.e., the nearest matrix with all its eigenvalues in a prescribed closed set $\Omega$. Distances are measured in the Frobenius norm. An important special case is finding the nearest Hurwitz or Schur stable matrix, which has applications in systems theory. We describe a reformulation of the task as an optimization problem on the Riemannian manifold of orthogonal (or unitary) matrices. The problem can then be solved using standard methods from the theory of Riemannian optimization. The resulting algorithm is remarkably fast on small-scale and medium-scale matrices, and returns directly a Schur factorization of the minimizer, sidestepping the numerical difficulties associated with eigenvalues with high multiplicity.


title: Concurrent Reference Counting and Resource Management in Wait-free Constant Time

abstract: A common problem when implementing concurrent programs is efficiently protecting against unsafe races between processes reading and then using a resource (e.g., memory blocks, file descriptors, or network connections) and other processes that are concurrently overwriting and then destructing the same resource. Such read-destruct races can be protected with locks, or with lock-free solutions such as hazard-pointers or read-copy-update (RCU). In this paper we describe a method for protecting read-destruct races with expected constant time overhead, $O(P^2)$ space and $O(P^2)$ delayed destructs, and with just single word atomic memory operations (reads, writes, and CAS). It is based on an interface with four primitives, an acquire-release pair to protect accesses, and a retire-eject pair to delay the destruct until it is safe. We refer to this as the acquire-retire interface. Using the acquire-retire interface, we develop simple implementations for three common use cases: (1) memory reclamation with applications to stacks and queues, (2) reference counted objects, and (3) objects manage by ownership with moves, copies, and destructs. The first two results significantly improve on previous results, and the third application is original. Importantly, all operations have expected constant time overhead.


title: When symmetries are not enough: a hierarchy of hard Constraint Satisfaction Problems

abstract: We produce a class of $\omega$-categorical structures with finite signature by applying a model-theoretic construction -- a refinement of the Hrushosvki-encoding -- to $\omega$-categorical structures in a possibly infinite signature. We show that the encoded structures retain desirable algebraic properties of the original structures, but that the constraint satisfaction problems (CSPs) associated with these structures can be badly behaved in terms of computational complexity. This method allows us to systematically generate $\omega$-categorical templates whose CSPs are complete for a variety of complexity classes of arbitrarily high complexity, and $\omega$-categorical templates that show that membership in any given complexity class cannot be expressed by a set of identities on the polymorphisms. It moreover enables us to prove that recent results about the relevance of topology on polymorphism clones of $\omega$-categorical structures also apply for CSP templates, i.e., structures in a finite language. Finally, we obtain a concrete algebraic criterion which could constitute a description of the delineation between tractability and NP-hardness in the dichotomy conjecture for first-order reducts of finitely bounded homogeneous structures.


title: Optical lattice experiments at unobserved conditions and scales through generative adversarial deep learning

abstract: Machine learning provides a novel avenue for the study of experimental realizations of many-body systems, and has recently been proven successful in analyzing properties of experimental data of ultracold quantum gases. We here show that deep learning succeeds in the more challenging task of modelling such an experimental data distribution. Our generative model (RUGAN) is able to produce snapshots of a doped two-dimensional Fermi-Hubbard model that are indistinguishable from previously reported experimental realizations. Importantly, it is capable of accurately generating snapshots at conditions for which it did not observe any experimental data, such as at higher doping values. On top of that, our generative model extracts relevant patterns from small-scale examples and can use these to construct new configurations at a larger size that serve as a precursor to observations at scales that are currently experimentally inaccessible. The snapshots created by our model---which come at effectively no cost---are extremely useful as they can be employed to quantitatively test new theoretical developments under conditions that have not been explored experimentally, parameterize phenomenological models, or train other, more data-intensive, machine learning methods. We provide predictions for experimental observables at unobserved conditions and benchmark these against modern theoretical frameworks. The deep learning method we develop here is broadly applicable and can be used for the efficient large-scale simulation of equilibrium and nonequilibrium physical systems.

categories: physics.comp-ph cond-mat.dis-nn cond-mat.quant-gas cond-mat.str-el