title: Simulating gravitational motion, gas dynamics, and structure in the cosmos

abstract: We provide introductory explanations and illustrations of the $N$-body hydrodynamics code Charm N-body GrAvity solver (ChaNGa). ChaNGa simulates the gravitational motion and gas dynamics of matter in space, with the goal of modeling galactic and/or cosmological structure and evolution. We discuss the algorithm for leapfrog integration and smoothed particle hydrodynamics and computer science concepts used by the program, including the binary data structure for the particle positions. Our presentation borrows from the doctoral dissertation of J.\ G.\ Stadel. Problems are provided in order to use ChaNGa to learn or solidify some cosmological concepts.


title: Kaluza-Klein Five-Point Functions from $\textrm{AdS}_5\times S_5$ Supergravity

abstract: We continue to explore the bootstrap approach to five-point correlation functions for IIB supergravity on $AdS_5\times S^5$. Building on the result of [1], we develop an improved algorithm that allows us to more efficiently compute correlators of higher Kaluza-Klein modes. The new method uses only factorization and a superconformal twist, and is entirely within Mellin space where the analytic structure of holographic correlators is simpler. Using this method, we obtain in a closed form all five-point functions of the form $\langle pp222\rangle$, extending the earlier result for $p=2$. As a byproduct of our analysis, we also obtain explicit results for spinning four-point functions of higher Kaluza-Klein modes.


title: Quantifying Theory in Politics: Identification, Interpretation and the Role of Structural Methods

abstract: The best empirical research in political science clearly defines substantive parameters of interest, presents a set of assumptions that guarantee its identification, and uses an appropriate estimator. We argue for the importance of explicitly integrating rigorous theory into this process and focus on the advantages of doing so. By integrating theoretical structure into one's empirical strategy, researchers can quantify the effects of competing mechanisms, consider the ex-ante effects of new policies, extrapolate findings to new environments, estimate model-specific theoretical parameters, evaluate the fit of a theoretical model, and test competing models that aim to explain the same phenomena. As a guide to such a methodology, we provide an overview of structural estimation, including formal definitions, implementation suggestions, examples, and comparisons to other methods.


title: Emulating the measurement postulates of quantum mechanics via non-Hermitian Hamiltonian

abstract: Ever since the formulation of quantum mechanics, there is very little understanding of the process of the collapse of a wavefunction. We have proposed a dynamical model to emulate the measurement postulates of quantum mechanics. We postulate that a non-Hermitian Hamiltonian operates during the process of measurement, which evolves any state to an attracting equilibrium state, thus, mimicking a "collapse". We demonstrate this using a 2-level system and then extend it to an N-level system. For a 2-level system, we also demonstrate that the dynamics generated by the Lindblad master equation can be replicated as an incoherent sum of the evolution by two separate non-Hermitian Hamiltonians.


title: Linear functionals and $\Delta$- coherent pairs of the second kind

abstract: We classify all the \emph{$\Delta$-}coherent pairs of measures of the second kind on the real line. We obtain $5$ cases, corresponding to all the families of discrete semiclassical orthogonal polynomials of class $s\leq1.$


title: Some Refinements of Formulae Involving Floor and Ceiling Functions

abstract: The floor and ceiling functions appear often in mathematics and manipulating sums involving floors and ceilings is a subtle game. Fortunately, the well-known textbook Concrete Mathematics provides a nice introduction with a number of techniques explained and a number of single or double sums treated as exercises. For two such double sums we provide their single-sum analogues. These closed-form identities are given in terms of a dual partition of the multiset (regarded as a partition) of all b-ary digits of a nonnegative integer. We also present the double- and single-sum analogues involving the fractional part function and the shifted fractional part function.


title: Turing instability in a diffusive predator-prey model with multiple Allee effect and herd behavior

abstract: Diffusion-driven instability and bifurcation analysis are studied in a predator-prey model with herd behavior and quadratic mortality by incorporating multiple Allee effect into prey species. The existence and stability of the equilibria of the system are studied. And bifurcation behaviors of the system without diffusion are shown. The sufficient and necessary conditions for Turing instability occurring are obtained. And the stability and the direction of Hopf and steady state bifurcations are explored by using the normal form method. Furthermore, some numerical simulations are presented to support our theoretical analysis. We found that too large diffusion rate of prey prevents Turing instability from emerging. Finally, we summarize our findings in the conclusion.


title: Power Availability of PV plus Thermal Batteries in real-world electric power grids

abstract: As variable renewable energy sources comprise a growing share of total electricity generation, energy storage technologies are becoming increasingly critical for balancing energy generation and demand. In this study, we modeled an existing thermal energy storage unit with estimated capital costs that are sufficiently low to enable large-scale deployment in the electric power system. Our analysis emphasizes the value of using such units to cost-effectively improve renewable energy dispatchability. This study modeled an existing real-world grid rather than simulating hypothetical future electric power systems. The storage unit coupled with a photovoltaic (PV) system was modeled with different storage capacities, whereas each storage unit size had various discharge capacities. The modeling was performed under a baseline case with no emission constraints and under hypothetical scenarios in which CO$_2$ emissions were reduced. The results show that power availability increases with increasing storage size and vastly increases in the hypothetical CO$_2$ reduction scenarios, as the storage unit is utilized differently. When CO$_2$ emissions are reduced, the power system must be less dependent on fossil fuel technologies that currently serve the grid, and thus rely more on the power that is served from the PV + storage unit. The proposed approach can provide increased knowledge to power system planners regarding how adding PV + storage systems to existing grids can contribute to the efficient stepwise decarbonization of electric power systems.


title: The role of electron correlations in the electronic structure of putative Chern magnet TbMn$_6$Sn$_6$

abstract: A member of the RMn$_6$Sn$_6$ rare-earth family materials, TbMn$_6$Sn$_6$, recently showed experimental signatures of the realization of a quantum-limit Chern magnet. In this work, we use quantum Monte Carlo (QMC) and density functional theory with Hubbard $U$ (DFT$+U$) calculations to examine the electronic structure of TbMn$_6$Sn$_6$. To do so, we optimize accurate, correlation-consistent pseudopotentials for Tb and Sn using coupled-cluster and configuration-interaction (CI) methods. We find that DFT$+U$ and single-reference QMC calculations suffer from the same overestimation of the magnetic moments as meta-GGA and hybrid density functional approximations. Our findings point to the need for improved orbitals/wavefunctions for this class of materials, such as natural orbitals from CI, or for the inclusion of multi-reference effects that capture the static correlations for an accurate prediction of magnetic properties. DFT$+U$ with Mn magnetic moments adjusted to experiment predict the Dirac crossing in bulk to be close to the Fermi level, within $\sim 120$ meV, in agreement with the experiments. Our non-stoichiometric slab calculations show that the Dirac crossing approaches even closer to the Fermi level, suggesting the possible realization of Chern magnetism in this limit.

categories: cond-mat.str-el cond-mat.mtrl-sci physics.chem-ph physics.comp-ph


title: Amazing behavior and transition to chaos of some sequences using Collatz like problems and Quibic duffing

abstract: In this paper we shall show amazing behavior of some discrete maps using Collatze like problems and some advanced theories in analytic number theory and dynamical system,we have investigated the driven cubic-quintic Duffing equation such that , We were able to predict the number of limit cycles around the equilibrium and to develop a theoretical approach to chaos suppression in damped driven systems using Collatze like problem sequences , some new results regarding behavior of that sequence are presented.


title: On the Maximum Atom-Bond Sum-Connectivity Index of Graphs

abstract: The atom-bond sum-connectivity (ABS) index of a graph $G$ with edges $e_1,\cdots,e_m$ is the sum of the numbers $\sqrt{1-2(d_{e_i}+2)^{-1}}$ over $1\le i \le m$, where $d_{e_i}$ is the number of edges adjacent with $e_i$. In this paper, we study the maximum values of the ABS index over graphs with given parameters. More specifically, we determine the maximum ABS index of connected graphs of a given order and with a fixed (i) minimum degree, (ii) maximum degree, (iii) chromatic number, (iv) independence number, or (v) number of pendent vertices. We also characterize the graphs attaining the maximum ABS values in all of these classes.


title: Two-population Bayesian hierarchical model of type Ia supernovae

abstract: The currently used standardisation of type Ia supernovae results in Hubble residuals whose physical origin is unaccounted for. Here, we present a complete physical interpretation of the Hubble residuals based on a novel Bayesian hierarchical model of type Ia supernovae in which latent variables describing intrinsic and extrinsic (dust related) supernova properties originate from two populations. Fitting the model to SALT2 light curve parameters of supernovae in the Hubble flow we find strong ($4\sigma$) evidence for the presence of two overlapping, but distinct, populations differentiated primarily by their mean SALT2 shape parameter (stretch) x_1. Supernovae from the population with predominantly slow decliners (higher average x_1) are found to be intrinsically bluer (mean SALT2 colour c=-0.11) and twice as reddened by dust (mean reddening E(B-V)=0.10) than those from the opposite population dominated by fast decliners (lower average x_1) with c=-0.04$ and E(B-V)=0.05. The inferred extinction coefficient R_B in both supernova populations follows a broad (scatter 0.9) distribution with a mean of 4.1, which coincides closely with the mean extinction law measured in the Milky Way. We also find that the supernova data favour a peaked (two-tailed) distribution of selective extinction E(B-V) over the commonly adopted exponential model (one-tailed with maximum fixed at 0). Our model provides a complete explanation of the distribution of supernova light curve parameters in terms of extinction properties and the above-mentioned differences between the two populations, without the need for introducing any intrinsic scatter.

categories: astro-ph.SR astro-ph.CO astro-ph.GA astro-ph.HE


title: Extended Watson-Harkins Sum

abstract: The Watson-Harkins sum involving the product of the cosine and cosecant functions is extended to derive the finite sum of generalized Hurwitz-Lerch Zeta functions is derived in terms of the Hurwitz-Lerch Zeta function. A transformation formula arises for various finite values of the parameters involved. The finite product of trigonometric functions are also derived. All the results in this work are new.


title: Linear Absorption Spectrum of the Spin-Boson Model Studied by Extended Hierarchical Equations of Motion

abstract: With a decomposition scheme for the bath correlation function, the hierarchical equations of motion (HEOM) are extended to the zero-temperature sub-Ohmic and Ohmic spin-boson model. We investigate the linear absorption spectrum of the sub-Ohmic and Ohmic spin-boson model at zero temperature. By applying the extended HEOM, the equilibrium spin dynamics are obtained approximately. Then the linear response function is calculated according to the Kubo formula rewritten in Liouville space. To explore the essence of phase transition, we compute the linear absorption spectrum defined by the linear response function of a dipole moment. By analyzing the peak structure of the linear response spectrum, we get the dependence of linear absorption spectrum with Kondo parameter and different bath exponents. The spin relaxation dynamics are also calculated to explore the coherent-incoherent dynamic transition (CI) and delocalizedlocalized phase transition (DL). The corresponding phase diagram of DL and CI transition are also obtained. We pose a energy level picture to understand the different mechanism of DL phase transition between the deep sub-Ohmic and Ohmic spin-boson model.


title: Symmetry groups for social preference functions

abstract: We introduce the anonymity group, the neutrality group and the symmetry group of a social preference function. Inspired by a problem posed by Kelly in 1991 and remained unsolved, we investigate the problem of recognizing which permutation groups may arise as anonymity, neutrality and symmetry group of a social preference function. A complete description is found for the neutrality groups and a sufficient condition, which largely encompasses the problem, is found for the anonymity groups. Using the concept of orbit extension of a group $U$, we formulate manageable necessary conditions for being $U$ an anonymity or a symmetry group. Our research deeply interacts with problems of representability by Boolean functions shedding light on them.


title: Low progress math in a high performing system

abstract: Math anxiety negatively relates to math performance. This negative relationship may be exacerbated in low-progress math learners. However, there are limited studies on math anxiety among low progress learners in a paradoxically high performing education system like Singapore. To fill this research gap, this research analysed the anxiety profiles of 151 students who were in the math learning support intervention program administered by the Ministry of Education, Singapore (MOE). We examined the complex relationship centred in math anxiety with relevant variables such as demographic characteristics, working memory and math performance. Limitations and future directions are discussed.


title: A generalized series expansion of the arctangent function based on the enhanced midpoint integration

abstract: In this work we derive a generalized series expansion of the acrtangent function by using the enhanced midpoint integration (EMI). Algorithmic implementation of the generalized series expansion utilizes two-step iteration without surd and complex numbers. The computational test we performed reveals that such a generalization improves accuracy in computation of the arctangent function by many orders of the magnitude with increasing integer $M$, associated with subintervals in the EMI formula. The generalized series expansion may be promising for practical applications. It may be particularly useful in practical tasks, where extensive computations with arbitrary precision floating points are needed. The algorithmic implementation of the generalized series expansion of the arctangent function shows a rapid convergence rate in the computation of digits of $\pi$ in the Machin-like formulas.


title: A note on the Lambert W function: Bernstein and Stieltjes properties for a creep model in Linear Viscoelasticity

abstract: The purpose of this note is to propose an application of the Lambert W function in linear viscoelasticity based on the Bernstein and Stieltjes properties of this function. In particular we recognize the role of the main branch W_0(t) in a peculiar model of creep with two spectral functions in frequency that completely characterize the creep model. In order to calculate these spectral functions, it turns out that the conjugate symmetry property of the Lambert W function along its branch cut on the negative real axis is essential. We supplement our analysis by computing the corresponding relaxation function and providing the plots of all computed functions.


title: Fire smoke dispersion inside and outside of a warehouse building in moderate and strong wind conditions

abstract: A moderate-sized (8 MW) fire was modeled in a warehouse located in an urban area with diverse architecture, under wind conditions. The study investigates the impact of neighboring architecture on the wind field and the resulting effect on smoke control performance and smoke dispersion in the vicinity of the building. 25 CFD simulations were performed using ANSYS Fluent CFD code, for wind velocities of 5 m/s (moderate), 10 m/s (strong) at 12 wind angles each (0 - 330$^\circ$, 30$^\circ$ increment), and at 0 m/s for reference. The results show that the performance of the building smoke venting system was impacted by the wind, reaching 74% - 114% and 78% - 158% of the reference mass flow at moderate and strong wind velocities, respectively. Interestingly, in multiple cases, the exhaust flow rate of the ventilators was increased rather than hindered by the wind. This was attributed to the arrangement of inlets and outlets on the facades and the resulting pressure difference between the wall and roof openings. The smoke plume was highly dependent on the wind angle and the type of architecture upwind and downwind of the fire. Significant urban canyon effects and large vortices forming behind tall buildings were observed, leading to smoke accumulation at a significant distance from the fire.


title: Some hyperbolicity revisited and robust transitivity

abstract: In this article we revisit the notion of \textit{Some Hyperbolicity} introduced by Pujals and Sambarino in \cite{PuSa}. We present a more general definition, that in particular can be applied to the symplectic context (something that the previous couldn't). As an application we construct $C^1$ robustly transitive derived from Anosov diffeomorphisms with mixed behaviour on center leaves.


title: Sample Complexity of Probability Divergences under Group Symmetry

abstract: We rigorously quantify the improvement in the sample complexity of variational divergence estimations for group-invariant distributions. In the cases of the Wasserstein-1 metric and the Lipschitz-regularized $\alpha$-divergences, the reduction of sample complexity is proportional to an ambient-dimension-dependent power of the group size. For the maximum mean discrepancy (MMD), the improvement of sample complexity is more nuanced, as it depends on not only the group size but also the choice of kernel. Numerical simulations verify our theories.


title: Spectral radius of graphs of given size with forbidden subgraphs

abstract: Let $\rho(G)$ be the spectral radius of a graph $G$ with $m$ edges. Let $S_{m-k+1}^{k}$ be the graph obtained from $K_{1,m-k}$ by adding $k$ disjoint edges within its independent set. Nosal's theorem states that if $\rho(G)>\sqrt{m}$, then $G$ contains a triangle. Zhai and Shu showed that any non-bipartite graph $G$ with $m\geq26$ and $\rho(G)\geq\rho(S_{m}^{1})>\sqrt{m-1}$ contains a quadrilateral unless $G\cong S_{m}^{1}$ [M.Q. Zhai, J.L. Shu, Discrete Math. 345 (2022) 112630]. Wang proved that if $\rho(G)\geq\sqrt{m-1}$ for a graph $G$ with size $m\geq27$, then $G$ contains a quadrilateral unless $G$ is one of four exceptional graphs [Z.W. Wang, Discrete Math. 345 (2022) 112973]. In this paper, we show that any non-bipartite graph $G$ with size $m\geq51$ and $\rho(G)\geq\rho(S_{m-1}^{2})>\sqrt{m-2}$ contains a quadrilateral unless $G$ is one of three exceptional graphs. Moreover, we show that if $\rho(G)\geq\rho(S_{\frac{m+4}{2},2}^{-})$ for a graph $G$ with even size $m\geq74$, then $G$ contains a $C_{5}^{+}$ unless $G\cong S_{\frac{m+4}{2},2}^{-}$, where $C_{t}^{+}$ denotes the graph obtained from $C_{t}$ and $C_{3}$ by identifying an edge, $S_{n,k}$ denotes the graph obtained by joining each vertex of $K_{k}$ to $n-k$ isolated vertices and $S_{n,k}^{-}$ denotes the graph obtained by deleting an edge incident to a vertex of degree two, respectively.


title: Topological Order from Measurements and Feed-Forward on a Trapped Ion Quantum Computer

abstract: Quantum systems evolve in time in one of two ways: through the Schr\"odinger equation or wavefunction collapse. So far, deterministic control of quantum many-body systems in the lab has focused on the former, due to the probabilistic nature of measurements. This imposes serious limitations: preparing long-range entangled states, for example, requires extensive circuit depth if restricted to unitary dynamics. In this work, we use mid-circuit measurement and feed-forward to implement deterministic non-unitary dynamics on Quantinuum's H1 programmable ion-trap quantum computer. Enabled by these capabilities, we demonstrate for the first time a constant-depth procedure for creating a toric code ground state in real-time. In addition to reaching high stabilizer fidelities, we create a non-Abelian defect whose presence is confirmed by transmuting anyons via braiding. This work clears the way towards creating complex topological orders in the lab and exploring deterministic non-unitary dynamics via measurement and feed-forward.


title: PyGlove: Efficiently Exchanging ML Ideas as Code

abstract: The increasing complexity and scale of machine learning (ML) has led to the need for more efficient collaboration among multiple teams. For example, when a research team invents a new architecture like "ResNet," it is desirable for multiple engineering teams to adopt it. However, the effort required for each team to study and understand the invention does not scale well with the number of teams or inventions. In this paper, we present an extension of our PyGlove library to easily and scalably share ML ideas. PyGlove represents ideas as symbolic rule-based patches, enabling researchers to write down the rules for models they have not seen. For example, an inventor can write rules that will "add skip-connections." This permits a network effect among teams: at once, any team can issue patches to all other teams. Such a network effect allows users to quickly surmount the cost of adopting PyGlove by writing less code quicker, providing a benefit that scales with time. We describe the new paradigm of organizing ML through symbolic patches and compare it to existing approaches. We also perform a case study of a large codebase where PyGlove led to an 80% reduction in the number of lines of code.


title: Critical numerosity in collective behavior

abstract: Natural collectives, despite comprising individuals who may not know their numerosity, can exhibit behaviors that depend sensitively on it. This paper proves that the collective behavior of number-oblivious individuals can even have a critical numerosity, above and below which it qualitatively differs. We formalize the concept of critical numerosity in terms of a family of zero--one laws and introduce a model of collective motion, called chain activation and transport (CAT), that has one. CAT describes the collective motion of $n \geq 2$ individuals as a Markov chain that rearranges $n$-element subsets of the $d$-dimensional grid, $m < n$ elements at a time. According to the individuals' dynamics, with each step, CAT removes $m$ elements from the set and then progressively adds $m$ elements to the boundary of what remains, in a way that favors the consecutive addition and removal of nearby elements. This paper proves that, if $d \geq 3$, then CAT has a critical numerosity of $n_c = 2m+2$ with respect to the behavior of its diameter. Specifically, if $n < n_c$, then the elements form one "cluster," the diameter of which has an a.s.--finite limit infimum. However, if $n \geq n_c$, then there is an a.s.--finite time at which the set consists of clusters of between $m+1$ and $2m+1$ elements, and forever after which these clusters grow apart, resulting in unchecked diameter growth. The existence of critical numerosities means that collectives can exhibit "phase transitions" that are governed purely by their numerosity and not, for example, their density or the strength of their interactions. This fact challenges prevalent beliefs about collective behavior and suggests new functionality for programmable matter. More broadly, it demonstrates an opportunity to explore the possible behaviors collectives through the study of random processes that rearrange sets.


title: Efficient Constrained Codes That Enable Page Separation in Modern Flash Memories

abstract: The pivotal storage density win achieved by solid-state devices over magnetic devices recently is a result of multiple innovations in physics, architecture, and signal processing. Constrained coding is used in Flash devices to increase reliability via mitigating inter-cell interference. Recently, capacity-achieving constrained codes were introduced to serve that purpose. While these codes result in minimal redundancy, they result in non-negligible complexity increase and access speed limitation since pages cannot be read separately. In this paper, we suggest new constrained coding schemes that have low-complexity and preserve the desirable high access speed in modern Flash devices. The idea is to eliminate error-prone patterns by coding data either only on the left-most page (binary coding) or only on the two left-most pages ($4$-ary coding) while leaving data on all the remaining pages uncoded. Our coding schemes are systematic and capacity-approaching. We refer to the proposed schemes as read-and-run (RR) constrained coding schemes. The $4$-ary RR coding scheme is introduced to limit the rate loss. We analyze the new RR coding schemes and discuss their impact on the probability of occurrence of different charge levels. We also demonstrate the performance improvement achieved via RR coding on a practical triple-level cell Flash device.


title: Transformers in Action Recognition: A Review on Temporal Modeling

abstract: In vision-based action recognition, spatio-temporal features from different modalities are used for recognizing activities. Temporal modeling is a long challenge of action recognition. However, there are limited methods such as pre-computed motion features, three-dimensional (3D) filters, and recurrent neural networks (RNN) for modeling motion information in deep-based approaches. Recently, transformers success in modeling long-range dependencies in natural language processing (NLP) tasks has gotten great attention from other domains; including speech, image, and video, to rely entirely on self-attention without using sequence-aligned RNNs or convolutions. Although the application of transformers to action recognition is relatively new, the amount of research proposed on this topic within the last few years is astounding. This paper especially reviews recent progress in deep learning methods for modeling temporal variations. It focuses on action recognition methods that use transformers for temporal modeling, discussing their main features, used modalities, and identifying opportunities and challenges for future research.


title: Variational waveguide QED simulators

abstract: Waveguide QED simulators are analogue quantum simulators made by quantum emitters interacting with one-dimensional photonic band-gap materials. One of their remarkable features is that they can be used to engineer tunable-range emitter interactions. Here, we demonstrate how these interactions can be a resource to develop more efficient variational quantum algorithms for certain problems. In particular, we illustrate their power in creating wavefunction ans\"atze that capture accurately the ground state of quantum critical spin models (XXZ and Ising) with less gates and optimization parameters than other variational ans\"atze based on nearest-neighbor or infinite-range entangling gates. Finally, we study the potential advantages of these waveguide ans\"atze in the presence of noise. Overall, these results evidence the potential of using the interaction range as a variational parameter and place waveguide QED simulators as a promising platform for variational quantum algorithms.


title: Real-Time Traffic End-of-Queue Detection and Tracking in UAV Video

abstract: Highway work zones are susceptible to undue accumulation of motorized vehicles which calls for dynamic work zone warning signs to prevent accidents. The work zone signs are placed according to the location of the end-of-queue of vehicles which usually changes rapidly. The detection of moving objects in video captured by Unmanned Aerial Vehicles (UAV) has been extensively researched so far, and is used in a wide array of applications including traffic monitoring. Unlike the fixed traffic cameras, UAVs can be used to monitor the traffic at work zones in real-time and also in a more cost-effective way. This study presents a method as a proof of concept for detecting End-of-Queue (EOQ) of traffic by processing the real-time video footage of a highway work zone captured by UAV. EOQ is detected in the video by image processing which includes background subtraction and blob detection methods. This dynamic localization of EOQ of vehicles will enable faster and more accurate relocation of work zone warning signs for drivers and thus will reduce work zone fatalities. The method can be applied to detect EOQ of vehicles and notify drivers in any other roads or intersections too where vehicles are rapidly accumulating due to special events, traffic jams, construction, or accidents.


title: CLASH: Contrastive learning through alignment shifting to extract stimulus information from EEG

abstract: Stimulus-evoked EEG data has a notoriously low signal-to-noise ratio and high inter-subject variability. We propose a novel paradigm for the self-supervised extraction of stimulus-related brain response data: a model is trained to extract similar information between two time-aligned segments of EEG in response to the same stimulus. The extracted information can subsequently be used to obtain better results in downstream tasks that utilize the response to the stimulus. We show the efficacy of our method for a downstream task of decoding the speech envelope from auditory EEG. Our method outperforms other state-of-the-art denoising techniques, improving reconstruction scores by 45\%. Additionally, we show that in contrast to the baseline denoising techniques, our method can be used with data of unseen subjects and stimuli without retraining, improving decoding performance by 19\% and 34\% over raw EEG for two holdout datasets. Finally, the last experiment reveals that the accuracies obtained in the CLASH paradigm are significantly correlated with the percentile of obtained reconstruction correlation on the null distribution. In general, we showed that the proposed paradigm is suitable to train deep learning models to extract stimulus information from EEG while being stimulus feature agnostic.


title: Learning a Fourier Transform for Linear Relative Positional Encodings in Transformers

abstract: We propose a new class of linear Transformers called FourierLearner-Transformers (FLTs), which incorporate a wide range of relative positional encoding mechanisms (RPEs). These include regular RPE techniques applied for sequential data, as well as novel RPEs operating on geometric data embedded in higher-dimensional Euclidean spaces. FLTs construct the optimal RPE mechanism implicitly by learning its spectral representation. As opposed to other architectures combining efficient low-rank linear attention with RPEs, FLTs remain practical in terms of their memory usage and do not require additional assumptions about the structure of the RPE mask. Besides, FLTs allow for applying certain structural inductive bias techniques to specify masking strategies, e.g. they provide a way to learn the so-called local RPEs introduced in this paper and give accuracy gains as compared with several other linear Transformers for language modeling. We also thoroughly test FLTs on other data modalities and tasks, such as image classification, 3D molecular modeling, and learnable optimizers. To the best of our knowledge, for 3D molecular data, FLTs are the first Transformer architectures providing linear attention and incorporating RPE masking.


title: Geometric-algebraic approach to acid-base equilibrium

abstract: Closed-form expressions for the pH of ideal weak acid and weak acid buffer solutions (and their titration with a strong base) were obtained and analyzed with the aid of computer algebra systems. These expressions are used to evaluate, without the use of numerical approximations, the precision and accuracy of different approximations commonly employed to calculate the pH. The closed-form expression for the pH of a buffer solution is used to obtain an analytical expression for the pH stability when a strong base is added. Finally, it is shown that the pH expressions for all the systems under study can be obtained from an unique closed-form expression in terms of an effective weak acid constant and an effective acid concentration.


title: Constraints on compensated isocurvature perturbations from BOSS DR12 galaxy data

abstract: We use the BOSS DR12 galaxy power spectrum to constrain compensated isocurvature perturbations (CIP), which are opposite-sign primordial baryon and dark matter perturbations that leave the total matter density unchanged. Long-wavelength CIP $\sigma(\vec{x})$ enter the galaxy density contrast as $\delta_g(\vec{x}) \supset b_\sigma\sigma(\vec{x})$, with $b_\sigma$ the linear CIP galaxy bias parameter. We parameterize the CIP spectra as $P_{\sigma\sigma} = A^2P_{\mathcal{R}\mathcal{R}}$ and $P_{\sigma\mathcal{R}} = \xi\sqrt{P_{\sigma\sigma}P_{\mathcal{R}\mathcal{R}}}$, where $A$ is the CIP amplitude and $\xi$ is the correlation with the curvature perturbations $\mathcal{R}$. We find a significance of detection of $Ab_\sigma \neq 0$ of $1.8\sigma$ for correlated ($\xi = 1$) and $3.7\sigma$ for uncorrelated ($\xi = 0$) CIP. Large-scale data systematics have a bigger impact for uncorrelated CIP, which may explain the large significance of detection. The constraints on $A$ depend on the assumed priors for the $b_\sigma$ parameter, which we estimate using separate universe simulations. Assuming $b_\sigma$ values representative of all halos we find $\sigma_A = 145$ for correlated CIP and $\sigma_{|A|} = 475$ for uncorrelated CIP. Our strongest uncorrelated CIP constraint is for $b_\sigma$ representative of the $33\%$ most concentrated halos, $\sigma_{|A|} = 197$, which is better than the current CMB bounds $|A| \lesssim 360$. We also discuss the impact of the local primordial non-Gaussianity parameter $f_{\rm NL}$ in CIP constraints. Our results demonstrate the power of galaxy data to place tight constraints on CIP, and motivate works to understand better the impact of data systematics, as well as to determine theory priors for $b_\sigma$.


title: Aligning Robot and Human Representations

abstract: To act in the world, robots rely on a representation of salient task aspects: for example, to carry a coffee mug, a robot may consider movement efficiency or mug orientation in its behavior. However, if we want robots to act for and with people, their representations must not be just functional but also reflective of what humans care about, i.e. they must be aligned. We observe that current learning approaches suffer from representation misalignment, where the robot's learned representation does not capture the human's representation. We suggest that because humans are the ultimate evaluator of robot performance, we must explicitly focus our efforts on aligning learned representations with humans, in addition to learning the downstream task. We advocate that current representation learning approaches in robotics should be studied from the perspective of how well they accomplish the objective of representation alignment. We mathematically define the problem, identify its key desiderata, and situate current methods within this formalism. We conclude by suggesting future directions for exploring open challenges.


title: Mathematical results on harmonic polynomials

abstract: Some years ago, the harmonic polynomial was introduced in order to understand better the harmonic topological index; for instance, it allows to obtain bounds of the harmonic index of the main products of graphs. Here, we obtain several properties of this polynomial, and we prove that several properties of graphs can be deduced from their harmonic polynomials. Also, we show that two graphs with the same harmonic polynomial have to be similar.


title: A phase field model for high-cycle fatigue: total-life analysis

abstract: We present a generalised phase field formulation for predicting high-cycle fatigue in metals. Different fatigue degradation functions are presented, together with new damage accumulation strategies, to account for (i) a typical S-N curve slope, (ii) the fatigue endurance limit, and (iii) the mean stress effect. The numerical implementation exploits an efficient quasi-Newton monolithic solution strategy and Virtual S-N curves are computed for both smooth and notched samples. The comparison with experiments reveals that the model can accurately predict fatigue lives and endurance limits, as well as naturally capture the influence of the stress concentration factor and the load ratio.


title: Characterization and Generation of 3D Realistic Geological Particles with Metaball Descriptor based on X-Ray Computed Tomography

abstract: The morphology of geological particles is crucial in determining its granular characteristics and assembly responses. In this paper, Metaball-function based solutions are proposed for morphological characterization and generation of three-dimensional realistic particles according to the X-ray Computed Tomography (XRCT) images. For characterization, we develop a geometric-based Metaball-Imaging algorithm. This algorithm can capture the main contour of parental particles with a series of non-overlapping spheres and refine surface-texture details through gradient search. Four types of particles, hundreds of samples, are applied for evaluations. The result shows good matches on key morphological indicators(i.e., volume, surface area, sphericity, circularity, corey-shape factor, nominal diameter and surface-equivalent-sphere diameter), confirming its characterization precision. For generation, we propose the Metaball Variational Autoencoder. Assisted by deep neural networks, this method can generate 3D particles in Metaball form, while retaining coessential morphological features with parental particles. Additionally, this method allows for control over the generated shapes through an arithmetic pattern, enabling the generation of particles with specific shapes. Two sets of XRCT images different in sample number and geometric features are chosen as parental data. On each training set, one thousand particles are generated for validations. The generation fidelity is demonstrated through comparisons of morphologies and shape-feature distributions between generated and parental particles. Examples are also provided to demonstrate controllability on the generated shapes. With Metaball-based simulations frameworks previously proposed by the authors, these methods have the potential to provide valuable insights into the properties and behavior of actual geological particles.


title: Sequential pattern mining in educational data: The application context, potential, strengths, and limitations

abstract: Increasingly, researchers have suggested the benefits of temporal analysis to improve our understanding of the learning process. Sequential pattern mining (SPM), as a pattern recognition technique, has the potential to reveal the temporal aspects of learning and can be a valuable tool in educational data science. However, its potential is not well understood and exploited. This chapter addresses this gap by reviewing work that utilizes sequential pattern mining in educational contexts. We identify that SPM is suitable for mining learning behaviors, analyzing and enriching educational theories, evaluating the efficacy of instructional interventions, generating features for prediction models, and building educational recommender systems. SPM can contribute to these purposes by discovering similarities and differences in learners' activities and revealing the temporal change in learning behaviors. As a sequential analysis method, SPM can reveal unique insights about learning processes and be powerful for self-regulated learning research. It is more flexible in capturing the relative arrangement of learning events than the other sequential analysis methods. Future research may improve its utility in educational data science by developing tools for counting pattern occurrences as well as identifying and removing unreliable patterns. Future work needs to establish a systematic guideline for data preprocessing, parameter setting, and interpreting sequential patterns.


title: Weak $(p,k)$-Dirac manifolds

abstract: In this paper, we introduce the notion of a weak $(p,k)$-Dirac structure in $TM\oplus \Lambda^pT^*M$, where $0\leq k \leq p-1$. The weak $(p,k)$-Lagrangian condition has more informations than the $(p,k)$-Lagrangian condition and contains the $(p,k)$-Lagrangian condition. The weak $(p,0)$-Dirac structures are exactly the higher Dirac structures of order p introduced by N. Martinez Alba and H. Bursztyn in [23] and [6], respectively. The regular weak $(p,p-1)$-Dirac structure together with $(p,p-1)$-Lagrangian subspace at each point $m\in M$ have the multisymplectic foliation. Finally, we introduce the notion of weak $(p,k)$-Dirac morphism. We give the condition that a weak $(p,k)$-Dirac manifold is also a weak $(p,k)$-Dirac manifold after pulling back.


title: A neural operator-based surrogate solver for free-form electromagnetic inverse design

abstract: Neural operators have emerged as a powerful tool for solving partial differential equations in the context of scientific machine learning. Here, we implement and train a modified Fourier neural operator as a surrogate solver for electromagnetic scattering problems and compare its data efficiency to existing methods. We further demonstrate its application to the gradient-based nanophotonic inverse design of free-form, fully three-dimensional electromagnetic scatterers, an area that has so far eluded the application of deep learning techniques.


title: CAB: Empathetic Dialogue Generation with Cognition, Affection and Behavior

abstract: Empathy is an important characteristic to be considered when building a more intelligent and humanized dialogue agent. However, existing methods did not fully comprehend empathy as a complex process involving three aspects: cognition, affection and behavior. In this paper, we propose CAB, a novel framework that takes a comprehensive perspective of cognition, affection and behavior to generate empathetic responses. For cognition, we build paths between critical keywords in the dialogue by leveraging external knowledge. This is because keywords in a dialogue are the core of sentences. Building the logic relationship between keywords, which is overlooked by the majority of existing works, can improve the understanding of keywords and contextual logic, thus enhance the cognitive ability. For affection, we capture the emotional dependencies with dual latent variables that contain both interlocutors' emotions. The reason is that considering both interlocutors' emotions simultaneously helps to learn the emotional dependencies. For behavior, we use appropriate dialogue acts to guide the dialogue generation to enhance the empathy expression. Extensive experiments demonstrate that our multi-perspective model outperforms the state-of-the-art models in both automatic and manual evaluation.


title: Classical Electromagnetic Interaction of a Charge with a Solenoid or Toroid

abstract: The Aharonov-Bohm phase shift in a particle interference pattern when electrons pass a long solenoid is identical in form with the optical interference pattern shift when a piece of retarding glass is introduced into one path of a two-beam optical interference pattern. The particle interference-pattern deflection is a relativistic effect of order 1/c^2, though this relativity aspect is rarely mentioned in the literature. Here we give a thorough analysis of the classical electromagnetic aspects of the interaction between a solenoid or toroid and a charged particle. We point out the magnetic Lorentz force which the solenoid or toroid experiences due to a passing charge. Although analysis in the rest frame of the solenoid or toroid will involve back Faraday fields on the charge, the analysis in the inertial frame in which the charge is initially at rest involves forces due to only electric fields where forces are equal in magnitude and opposite in direction. The classical analysis is made using the Darwin Lagrangian. We point out that the classical analysis suggests an angular deflection independent of Planck's constant where the deflection magnitude is identical with that given by the traditional quantum analysis, but where the deflection direction is unambiguous.


title: The Classical Aharonov-Bohm Interaction as a Relativity Paradox

abstract: The situation of a charged particle passing down the symmetry axis through a magnetic toroid presents a relativity paradox; different inertial frames suggest different forces on the charge and on the toroid due to the unperturbed systems. We review the charge-toroid interaction and suggest that the magnetic Aharonov-Bohm situation is misunderstood because of unfamiliarity with the acceleration fields following from the Darwin Lagrangian, which go unmentioned in recent textbooks of classical electromagnetism.


title: Large $N$ von Neumann algebras and the renormalization of Newton's constant

abstract: I derive a family of Ryu--Takayanagi formulae that are valid in the large $N$ limit of holographic quantum error-correcting codes, and parameterized by a choice of UV cutoff in the bulk. The bulk entropy terms are matched with a family of von Neumann factors nested inside the large $N$ von Neumann algebra describing the bulk effective field theory. These factors are mapped onto one another by a family of conditional expectations, which are interpreted as a renormalization group flow for the code subspace. Under this flow, I show that the renormalizations of the area term and the bulk entropy term exactly compensate each other. This result provides a concrete realization of the ER=EPR paradigm, as well as an explicit proof of a conjecture due to Susskind and Uglum.


title: Rescattering effects in nucleon-to-meson form factors and application to tau-lepton-induced proton decay

abstract: Nucleon decays put extremely stringent bounds on baryon-number-violating interactions. However, in case the corresponding operators involve only $\tau$ leptons, the direct two-body decays, e.g., $p\to\pi^0 \tau^+$, are kinematically not allowed and nucleon decay can only proceed via an off-shell $\tau$, leading to $p\to \pi^0\ell^+\nu_\ell\bar\nu_\tau$. To calculate such processes, the momentum dependence of the form factors for nucleon-to-meson transitions, which describe the hadronization of the underlying process at the quark level, is needed. In this work, we point out new isospin and Fierz relations among such proton-kaon matrix elements and calculate the momentum dependence of the nucleon-to-meson form factors from the universal final-state interactions in terms of pion-nucleon or kaon-nucleon scattering phase shifts. We use these results to derive novel limits on the Wilson coefficients of the baryon-number-violating dimension-$6$ operators involving a $\tau$ lepton, which were previously unconstrained.


title: Hamiltonian cycles on Ammann-Beenker Tilings

abstract: We provide a simple algorithm for constructing Hamiltonian graph cycles (visiting every vertex exactly once) on the set of aperiodic two-dimensional Ammann-Beenker (AB) tilings. Using this result, and the discrete scale symmetry of AB tilings, we find exact solutions to a range of other problems which lie in the complexity class NP-Complete for general graphs. These include the equal-weight travelling salesperson problem, providing for example the most efficient route a scanning tunneling microscope tip could take to image the atoms of physical quasicrystals with AB symmetries; the longest path problem, whose solution demonstrates that collections of flexible molecules of any length can adsorb onto AB quasicrystal surfaces at density one, with possible applications to catalysis; and the 3-colouring problem, giving ground states for the $q$-state Potts model ($q\ge 3$) of magnetic interactions defined on the planar dual to AB, which may provide useful models for protein folding.

categories: cond-mat.stat-mech cond-mat.str-el math-ph math.MP


title: Efficient optimization of deep neural quantum states toward machine precision

abstract: Neural quantum states (NQSs) have emerged as a novel promising numerical method to solve the quantum many-body problem. However, it has remained a central challenge to train modern large-scale deep network architectures to desired quantum state accuracy, which would be vital in utilizing the full power of NQSs and making them competitive or superior to conventional numerical approaches. Here, we propose a minimum-step stochastic reconfiguration (MinSR) method that reduces the optimization complexity by orders of magnitude while keeping similar accuracy as compared to conventional stochastic reconfiguration. MinSR allows for accurate training on unprecedentedly deep NQS with up to 64 layers and more than $10^5$ parameters in the spin-1/2 Heisenberg $J_1$-$J_2$ models on the square lattice. We find that this approach yields better variational energies as compared to existing numerical results and we further observe that the accuracy of our ground state calculations approaches different levels of machine precision on modern GPU and TPU hardware. The MinSR method opens up the potential to make NQS superior as compared to conventional computational methods with the capability to address yet inaccessible regimes for two-dimensional quantum matter in the future.


title: GLASS: Generator for Large Scale Structure

abstract: We present GLASS, the Generator for Large Scale Structure, a new code for the simulation of galaxy surveys for cosmology, which iteratively builds a light cone with matter, galaxies, and weak gravitational lensing signals as a sequence of nested shells. This allows us to create deep and realistic simulations of galaxy surveys at high angular resolution on standard computer hardware and with low resource consumption. GLASS also introduces a new technique to generate transformations of Gaussian random fields (including lognormal) to essentially arbitrary precision, an iterative line-of-sight integration over matter shells to obtain weak lensing fields, and flexible modelling of the galaxies sector. We demonstrate that GLASS readily produces simulated data sets with per cent-level accurate two-point statistics of galaxy clustering and weak lensing, thus enabling simulation-based validation and inference that is limited only by our current knowledge of the input matter and galaxy properties.


title: New determination of the production cross section for $\gamma$ rays in the Galaxy

abstract: The flux of $\gamma$ rays is measured with unprecedented accuracy by the $\textit{Fermi}$ Large Area Telescope from 100 MeV to almost 1 TeV. In the future, the Cherenkov Telescope Array will have the capability to measure photons up to 100 TeV. To accurately interpret this data, precise predictions of the production processes, specifically the cross section for the production of photons from the interaction of cosmic-ray protons and helium with atoms of the ISM, are necessary. In this study, we determine new analytical functions describing the Lorentz-invariant cross section for $\gamma$-ray production in hadronic collisions. We utilize the limited total cross section data for $\pi^0$ production channels and supplement this information by drawing on our previous analyses of charged pion production to infer missing details. In this context, we highlight the need for new data on $\pi^0$ production. Our predictions include the cross sections for all production channels that contribute down to the 0.5% level of the final cross section, namely $\eta$, $K^+$, $K^-$, $K^0_S$, and $K^0_L$ mesons as well as $\Lambda$, $\Sigma$, and $\Xi$ baryons. We determine the total differential cross section $d\sigma(p+p\rightarrow \gamma+X)/dE_{\gamma}$ from 10 MeV to 100 TeV with an uncertainty of 10% below 10 GeV of $\gamma$-ray energies, increasing to 20% at the TeV energies. We provide numerical tables and a script for the community to access our energy-differential cross sections, which are provided for incident proton (nuclei) energies from 0.1 to $10^7$ GeV (GeV/n).


title: Spin supplementary condition in quantum field theory: covariant SSC and physical state projection

abstract: The spin supplementary conditions are constraints on spin degrees of freedom in classical relativity which restricts physical degrees of freedom to rotations. It is argued that the equivalent constraints in quantum field theory are the projection conditions on polarisation tensors, which remove timelike/longitudinal polarisations from the physical spectrum. The claim is supported by three examples of massive spinning particles coupled to electromagnetism: Dirac and Proca fields in quantum field theory, and $\mathcal{N} = 1$ worldline QFT for classical worldline theory. This suggests a resolution to the apparent discrepancy between effective field theory description of massive higher-spin fields~\cite{Bern:2020buy,Bern:2022kto} and post-Newtonian effective field theory of spinning classical particles~\cite{Levi:2015msa}, where the former admits more unfixed parameters compared to the latter; the additional parameters are fixed by projection conditions and therefore are not tunable parameters.


title: The divergence theorem and nonlocal counterparts

abstract: We present a new proof of the classical divergence theorem in bounded domains. Our proof is based on a nonlocal analog of the divergence theorem and a rescaling argument. Main ingredients in the proof are nonlocal versions of the divergence and the normal derivative. We employ these to provide definitions of well-known nonlocal concepts such as the fractional perimeter.


title: Can the Universe decelerate in the future?

abstract: The possibility of an expanding decelerating Universe in the distant future is investigated in the context of a quintessence scalar field cosmology. Such a conceivable evolution is tested against SNe Ia and $H(z)$ cosmic chronometers data, and also through a model independent method based on Gaussian Processes. The scalar field model is an extension of the exponential Ratra-Peebles (RP) quintessential cosmology whose potential now depends on a pair of parameters ($\alpha, \lambda)$ and predicts a decelerated expansion in the future. Different from RP approach the $\alpha$ parameter allows for a decelerating cosmology in the future thereby frustrating the inevitable evolution for a de Sitter Cosmology as predicted by the cosmic concordance model ($\Lambda$CDM). The statistical model analysis is updated with the most recent SNe Ia and $H(z)$ data thereby obtaining $H_0 = 68.6\pm3.7$ km/s/Mpc, $\Omega_{\Phi0} = 0.735^{+0.083}_{-0.069} $, $\alpha < 6.56$ and $\lambda< 0.879 $ (at $2\sigma$ c.l.). It is also found that the extended RP model allows for a future deceleration both for $H(z)$ and SNe Ia data. In the (model-independent) Gaussian Processes analysis, however, future deceleration is allowed only in the case of $H(z)$ data.


title: Extracting the gamma-ray source-count distribution below the Fermi-LAT detection limit with deep learning

abstract: We reconstruct the extra-galactic gamma-ray source-count distribution, or $dN/dS$, of resolved and unresolved sources by adopting machine learning techniques. Specifically, we train a convolutional neural network on synthetic 2-dimensional sky-maps, which are built by varying parameters of underlying source-counts models and incorporate the Fermi-LAT instrumental response functions. The trained neural network is then applied to the Fermi-LAT data, from which we estimate the source count distribution down to flux levels a factor of 50 below the Fermi-LAT threshold. We perform our analysis using 14 years of data collected in the $(1,10)$ GeV energy range. The results we obtain show a source count distribution which, in the resolved regime, is in excellent agreement with the one derived from catalogued sources, and then extends as $dN/dS \sim S^{-2}$ in the unresolved regime, down to fluxes of $5 \cdot 10^{-12}$ cm$^{-2}$ s$^{-1}$. The neural network architecture and the devised methodology have the flexibility to enable future analyses to study the energy dependence of the source-count distribution.

categories: astro-ph.CO astro-ph.HE astro-ph.IM cs.LG


title: An empirical relation to estimate host galaxy stellar light from AGN spectra

abstract: Measurement of black hole mass for low-$z$ ($z\leq 0.8$) Active Galactic Nuclei (AGNs) is difficult due to the strong contribution from host galaxy stellar light necessitating detailed spectral decomposition to estimate the AGN luminosity. Here, we present an empirical relation to estimate host galaxy stellar luminosity from the optical spectra of AGNs at $z\leq 0.8$. The spectral data were selected from the fourteenth data release of the Sloan Digital Sky Survey (SDSS-DR14) quasar catalog having a signal-to-noise ratio at 5100 \AA (SNR$_{5100}$) $>$10 containing 11415 quasars. The median total luminosity (log ($L_\text{total}$/[erg s$^{-1}$])), stellar luminosity (log ($L_\text{star}$/[erg s$^{-1}$])), and AGN continuum luminosity (log ($L_\text{cont}$/[erg s$^{-1}$])) in our sample are 44.52, 44.06, and 44.30, respectively. We fit the AGN power-law continuum, host galaxy, and iron blend contribution, simultaneously over the entire available spectrum. We found the host galaxy fraction to anti-correlate with continuum luminosity and can be well-represented by a polynomial function, which can be used to correct the stellar light contribution from AGN spectra. We also found anti-correlation between host galaxy fraction and iron strength, Eddington ratio, and redshift. The empirical relation gives comparable results of host-fraction with the image decomposition method.


title: Modest dust settling in the IRAS04302+2247 Class I protoplanetary disk

abstract: We present new VLA observations, between 6.8mm and 66mm, of the edge-on Class~I disk IRAS04302+2247. Observations at 6.8mm and 9.2mm lead to the detection of thermal emission from the disk, while shallow observations at the other wavelengths are used to correct for emission from other processes. The disk radial brightness profile transitions from broadly extended in previous ALMA 0.9mm and 2.1mm observations to much more centrally brightened at 6.8mm and 9.2mm, which can be explained by optical depth effects. The radiative transfer modeling of the 0.9mm, 2.1mm, and 9.2mm data suggests that the grains are smaller than 1cm in the outer regions of the disk and allows us to obtain the first lower limit for the scale height of grains emitting at millimeter wavelengths in a protoplanetary disk. We find that the millimeter dust scale height is between 1au and 6au at a radius 100au from the central star, while the gas scale height is estimated to be about 7au, indicating a modest level of settling. The estimated dust height is intermediate between less evolved Class 0 sources, that are found to be vertically thick, and more evolved Class II sources, which show a significant level of settling. This suggests that we are witnessing an intermediate stage of dust settling.


title: $Q$-Deformed Rainbows: a Universal Simulator of Free Entanglement Spectra

abstract: The behaviour of correlations across a bipartition is an indispensable tool in diagnosing quantum phases of matter. Here we present a spin chain with position-dependent XX couplings and magnetic fields, that can reproduce arbitrary structure of free fermion correlations across a bipartition. In particular, by choosing appropriately the strength of the magnetic fields we can obtain any single particle energies of the entanglement spectrum with high fidelity. The resulting ground state can be elegantly formulated in terms of $q$-deformed singlets. To demonstrate the versatility of our method we consider certain examples, such as a system with homogeneous correlations and a system with correlations that follow a prime number decomposition. Hence, our entanglement simulator can be easily employed for the generation of arbitrary entanglement spectra with possible applications in quantum technologies and condensed matter physics.


title: Dynamical coherence in isotopy classes of fibered lifted partially hyperbolic diffeomorphisms

abstract: We introduce the notion of \textit{fibered lifted partially hyperbolic diffeomorphisms} and we prove that any partially hyperbolic diifeomorphism isotopic to a fibered lifted one where the isotopy take place inside partially hyperbolic systems is \textit{dynamically coherent.} Moreover we prove some global stability result: every two partially hyperbolic diffeomorphisms in the same connected component of a fibered lifted partially hyperbolic diffeomorphisms, are leaf conjugate.


title: On a continuous time model of gradient descent dynamics and instability in deep learning

abstract: The recipe behind the success of deep learning has been the combination of neural networks and gradient-based optimization. Understanding the behavior of gradient descent however, and particularly its instability, has lagged behind its empirical success. To add to the theoretical tools available to study gradient descent we propose the principal flow (PF), a continuous time flow that approximates gradient descent dynamics. To our knowledge, the PF is the only continuous flow that captures the divergent and oscillatory behaviors of gradient descent, including escaping local minima and saddle points. Through its dependence on the eigendecomposition of the Hessian the PF sheds light on the recently observed edge of stability phenomena in deep learning. Using our new understanding of instability we propose a learning rate adaptation method which enables us to control the trade-off between training stability and test set evaluation performance.


title: Trapping microparticles in a structured dark focus

abstract: We experimentally demonstrate stable trapping and controlled manipulation of silica microspheres in a structured optical beam consisting of a dark focus surrounded by light in all directions - the so-called Dark Focus Tweezer. Results from power spectrum and potential analysis demonstrate the non-harmonicity of the trapping potential landspace, which is reconstructed from experimental data in agreement to Lorentz-Mie numerical simulations. Applications of the dark tweezer in levitated optomechanics and biophysics are discussed.


title: Pauli-Villars' regularization of ghosts in path-integral string formulation

abstract: I consider Pauli-Villars' regulators for the ghosts in the path-integral string formulation and show how they preserve conformal invariance. I calculate the regulator contributions to the effective action and to the central charge and demonstrate the consistency of the mean-field quantization of the Nambu-Goto string in $2<d\leq26$. The higher-derivative corrections to the Liouville action are briefly considered for the Pauli-Villars and proper-time regularizations.


title: Fixed-kinetic Neural Hamiltonian Flows for enhanced interpretability and reduced complexity

abstract: Normalizing Flows (NF) are Generative models which are particularly robust and allow for exact sampling of the learned distribution. They however require the design of an invertible mapping, whose Jacobian determinant has to be computable. Recently introduced, Neural Hamiltonian Flows (NHF) are based on Hamiltonian dynamics-based Flows, which are continuous, volume-preserving and invertible and thus make for natural candidates for robust NF architectures. In particular, their similarity to classical Mechanics could lead to easier interpretability of the learned mapping. However, despite being Physics-inspired architectures, the originally introduced NHF architecture still poses a challenge to interpretability. For this reason, in this work, we introduce a fixed kinetic energy version of the NHF model. Inspired by physics, our approach improves interpretability and requires less parameters than previously proposed architectures. We then study the robustness of the NHF architectures to the choice of hyperparameters. We analyze the impact of the number of leapfrog steps, the integration time and the number of neurons per hidden layer, as well as the choice of prior distribution, on sampling a multimodal 2D mixture. The NHF architecture is robust to these choices, especially the fixed-kinetic energy model. Finally, we adapt NHF to the context of Bayesian inference and illustrate our method on sampling the posterior distribution of two cosmological parameters knowing type Ia supernovae observations.


title: The massless non-adjacent double off-shell scalar box integral -- branch cut structure and all-order epsilon expansion

abstract: We generalize the result of our recent paper on the massless single off-shell scalar box integral to the case of two non-adjacent end points off the light cone. An analytic result in $d=4-2\varepsilon$ dimensions is established in terms of four Gauss hypergeometric ${_2}\mathrm{F}_1$ functions respectively their single-valued counterparts. This allows for an explicit splitting of real and imaginary parts, as well as an all-order $\varepsilon$-expansion in terms of single-valued polylogarithms.


title: Accelerating exploration of Marine Cloud Brightening impacts on tipping points Using an AI Implementation of Fluctuation-Dissipation Theorem

abstract: Marine cloud brightening (MCB) is a proposed climate intervention technology to partially offset greenhouse gas warming and possibly avoid crossing climate tipping points. The impacts of MCB on regional climate are typically estimated using computationally expensive Earth System Model (ESM) simulations, preventing a thorough assessment of the large possibility space of potential MCB interventions. Here, we describe an AI model, named AiBEDO, that can be used to rapidly projects climate responses to forcings via a novel application of the Fluctuation-Dissipation Theorem (FDT). AiBEDO is a Multilayer Perceptron (MLP) model that uses maps monthly-mean radiation anomalies to surface climate anomalies at a range of time lags. By leveraging a large existing dataset of ESM simulations containing internal climate noise, we use AiBEDO to construct an FDT operator that successfully projects climate responses to MCB forcing, when evaluated against ESM simulations. We propose that AiBEDO-FDT can be used to optimize MCB forcing patterns to reduce tipping point risks while minimizing negative side effects in other parts of the climate.


title: Notes on the type classification of von Neumann algebras

abstract: These notes provide an explanation of the type classification of von Neumann algebras, which has made many appearances in recent work on entanglement in quantum field theory and quantum gravity. The goal is to bridge a gap in the literature between resources that are too technical for the non-expert reader, and resources that seek to explain the broad intuition of the theory without giving precise definitions. Reading these notes will provide you with: (i) an argument for why "factors" are the fundamental von Neumann algebras that one needs to study; (ii) an intuitive explanation of the type classification of factors in terms of renormalization schemes that turn unnormalizable positive operators into "effective density matrices;" (iii) a mathematical explanation of the different types of renormalization schemes in terms of the allowed traces on a factor; (iv) an intuitive characterization of type I and II factors in terms of their "standard forms;" and (v) a list of some interesting connections between type classification and modular theory, including the argument for why type III$_1$ factors are believed to be the relevant ones in quantum field theory. None of the material is new, but the pedagogy is different from other sources I have read; it is most similar in spirit to the recent work on gravity and the crossed product by Chandrasekaran, Longo, Penington, and Witten.


title: Surface tension of bulky colloids, capillarity under gravity and the microscopic origin of the Kardar-Parisi-Zhang equation

abstract: Experimental measurements of the surface tension of colloidal interfaces have long been in conflict with computer simulations. In this work we show that the surface tension of colloids as measured by surface fluctuations picks up a gravity dependent contribution which removes the discrepancy. The presence of this term puts a strong constraint on the structure of the interface which allows one to identify corrections to the fundamental equation of equilibrium capillarity and deduce bottom-up the microscopic origin of a growth model with close relation to the Kardar-Parisi-Zhang equation.

categories: cond-mat.soft cond-mat.mes-hall physics.app-ph physics.chem-ph physics.flu-dyn


title: Dynamic Arctic weather variability and connectivity

abstract: The rapidly shrinking Arctic sea ice is changing weather patterns and disrupting the balance of nature. Dynamics of Arctic weather variability (WV) plays a crucial role in weather forecasting and is closely related to extreme weather events. Yet, assessing and quantifying the WV for both local Arctic regions and its planetary impacts under anthropogenic climate change is still unknown. Here, we develop a complexity-based approach to systematically evaluate and analyze the dynamic behaviour of WV. We reveal that the WV within and around the Arctic is statistically correlated to the Arctic Oscillation at the intraseasonal time scale. We further find that the variability of the daily Arctic sea ice is increasing due to its dramatic decline under a warming climate. Unstable Arctic weather conditions can disturb regional weather patterns through atmospheric teleconnection pathways, resulting in higher risk to human activities and greater weather forecast uncertainty. A multivariate climate network analysis reveals the existence of such teleconnections and implies a positive feedback loop between the Arctic and global weather instabilities. This enhances the mechanistic understanding of the influence of Arctic amplification on mid-latitude severe weather. Our framework provides a fresh perspective on the linkage of complexity science, WV and the Arctic.


title: Asymmetric Certified Robustness via Feature-Convex Neural Networks

abstract: Recent works have introduced input-convex neural networks (ICNNs) as learning models with advantageous training, inference, and generalization properties linked to their convex structure. In this paper, we propose a novel feature-convex neural network architecture as the composition of an ICNN with a Lipschitz feature map in order to achieve adversarial robustness. We consider the asymmetric binary classification setting with one "sensitive" class, and for this class we prove deterministic, closed-form, and easily-computable certified robust radii for arbitrary $\ell_p$-norms. We theoretically justify the use of these models by characterizing their decision region geometry, extending the universal approximation theorem for ICNN regression to the classification setting, and proving a lower bound on the probability that such models perfectly fit even unstructured uniformly distributed data in sufficiently high dimensions. Experiments on Malimg malware classification and subsets of MNIST, Fashion-MNIST, and CIFAR-10 datasets show that feature-convex classifiers attain state-of-the-art certified $\ell_1$-radii as well as substantial $\ell_2$- and $\ell_{\infty}$-radii while being far more computationally efficient than any competitive baseline.


title: Correspondence between open bosonic systems and stochastic differential equations

abstract: Bosonic mean-field theories can approximate the dynamics of systems of $n$ bosons provided that $n \gg 1$. We show that there can also be an exact correspondence at finite $n$ when the bosonic system is generalized to include interactions with the environment and the mean-field theory is replaced by a stochastic differential equation. When the $n \to \infty$ limit is taken, the stochastic terms in this differential equation vanish, and a mean-field theory is recovered. Besides providing insight into the differences between the behavior of finite quantum systems and their classical limits given by $n \to \infty$, the developed mathematics can provide a basis for quantum algorithms that solve some stochastic nonlinear differential equations. We discuss conditions on the efficiency of these quantum algorithms, with a focus on the possibility for the complexity to be polynomial in the log of the stochastic system size. A particular system with the form of a stochastic discrete nonlinear Schr\"{o}dinger equation is analyzed in more detail.


title: Formulas for Hitting Times and Cover Times for Random Walks on Groups

abstract: Using the results of Ding, Lee, Peres [3], we develop formulas to compute the hitting times and cover times for random walks on groups. We developed an explicit formula for hitting times in terms of the irreducible representations of the group. We also have a way of computing cover times in terms of these hitting times. This computation is based on a quantity we indentified, which we call the volume growth function. And we believe that it is the right object to study in order to understand the cover time.


title: Elementary excitations of a system of one-dimensional chiral fermions with short-range interactions

abstract: We study general features of the excitation spectrum of a system of one-dimensional chiral spinless fermions with short-range interactions. We show that the nature of the elementary excitations of such a system depends strongly on the nonlinearity of the underlying dispersion of the fermions. In the case of quadratic nonlinearity, the low-momentum excitations are essentially fermionic quasiparticles and quasiholes, whereas the high-momentum ones are classical harmonic waves and solitons. In the case of cubic nonlinearity, the nature of the elementary excitations does not depend on momentum and is determined by the strength of the interactions. At a certain critical value of the interaction strength the excitation spectrum changes qualitatively, pointing to a dynamic phase transition in the system.


title: A first-principles method to calculate fourth-order elastic constants of solid materials

abstract: A first-principles method is presented to calculate elastic constants up to the fourth order of crystals with the cubic and hexagonal symmetries. The method relies on the numerical differentiation of the second Piola-Kirchhoff stress tensor and a density functional theory approach to compute the Cauchy stress tensors for a minimal list of strained configurations of a reference state. The number of strained configurations required to calculate the independent elastic constants of the second, third, and fourth order is 24 and 37 for crystals with the cubic and hexagonal symmetries, respectively. Here, this method is applied to five crystalline materials with the cubic symmetry (diamond, silicon, aluminum, silver, and gold) and two metals with the hexagonal close packing structure (beryllium and magnesium). Our results are compared to available experimental data and previous computational studies. Calculated linear and nonlinear elastic constants are also used, within a nonlinear elasticity treatment of a material, to predict values of volume and bulk modulus at zero temperature over an interval of pressures. To further validate our method, these predictions are compared to results obtained from explicit density functional theory calculations.


title: Towards an Understanding of Distributed Asymmetric Collaborative Visualization on Problem-solving

abstract: This paper provided empirical knowledge of the user experience for using collaborative visualization in a distributed asymmetrical setting through controlled user studies. With the ability to access various computing devices, such as Virtual Reality (VR) head-mounted displays, scenarios emerge when collaborators have to or prefer to use different computing environments in different places. However, we still lack an understanding of using VR in an asymmetric setting for collaborative visualization. To get an initial understanding and better inform the designs for asymmetric systems, we first conducted a formative study with 12 pairs of participants. All participants collaborated in asymmetric (PC-VR) and symmetric settings (PC-PC and VR-VR). We then improved our asymmetric design based on the key findings and observations from the first study. Another ten pairs of participants collaborated with enhanced PC-VR and PC-PC conditions in a follow-up study. We found that a well-designed asymmetric collaboration system could be as effective as a symmetric system. Surprisingly, participants using PC perceived less mental demand and effort in the asymmetric setting (PC-VR) compared to the symmetric setting (PC-PC). We provided fine-grained discussions about the trade-offs between different collaboration settings.


title: Measurement of the top quark mass using a profile likelihood approach with the lepton+jets final states in proton-proton collisions at $\sqrt{s}$ = 13 TeV

abstract: The mass of the top quark is measured in 36.3 fb$^{-1}$ of LHC proton-proton collision data collected with the CMS detector at $\sqrt{s}$ = 13 TeV. The measurement uses a sample of top quark pair candidate events containing one isolated electron or muon and at least four jets in the final state. For each event, the mass is reconstructed from a kinematic fit of the decay products to a top quark pair hypothesis. A profile likelihood method is applied using up to four observables to extract the top quark mass. The top quark mass is measured to be 171.77 $\pm$ 0.37 GeV. This approach significantly improves the precision over previous measurements.


title: A new method for interpreting well-to-well interference tests and quantifying the magnitude of production impact: Theory and applications in a multi-basin case study

abstract: Interference tests are used in shale reservoirs to evaluate the strength of connectivity between wells. The results inform engineering decisions about well spacing. In this paper, we propose a new procedure for interpreting interference tests. We fit the initial interference response with the solution to the 1D diffusivity equation at an offset observation point. It is advantageous to use the initial interference response, rather than the subsequent trend, because it is less affected by nonlinearities, time-varying boundary conditions, and uncertainties about flow geometry and regime. From the curve fit, we estimate the hydraulic diffusivity and conductivity of the fractures connecting the wells. For engineering purposes, it would be useful to quantify the impact of interference on well production. Thus, we seek a relationship between the degree of production interference (DPI) and an appropriate dimensionless quantity that can be derived from the estimate of fracture conductivity. Using simulations run under a wide range of conditions, we find that the classical definition for dimensionless fracture conductivity does not achieve a consistent prediction of DPI. This occurs because the dimensionless fracture conductivity is derived assuming radial flow geometry, but the dominant flow geometry during shale production is linear. We also find that the CPG (Chow Pressure Group) metric does not yield consistently accurate predictions of DPI. As an alternative, we derive a dimensionless quantity similar to the classical dimensionless fracture conductivity, but derived for linear, flow geometry. Using this approach, we calculate a dimensionless interference length that collapses all cases onto a single curve that predicts DPI as a function of fracture conductivity, well spacing, and formation properties. We conclude by applying the new method to field cases from the Anadarko and Delaware Basins.


title: A possible mechanism of the Kirkwood gap formations at the very beginning

abstract: The orbits of asteroids from the MPC catalogue of May 31, 2020 with absolute magnitudes H < 16m, in the 3:1, 5:2 and 2:1 mean motion resonances (MMRs) with Jupiter were selected. The number of the orbits in the 2:1 MMR is dozens times greater than in the other two resonances. There are fragments of parent bodies of neighbour asteroid families, in particular the Themis family, among bodies in the 2:1 MMR. Numerical calculations were performed to evaluate the evolution of the selected orbits over hundreds of thousand years. Perturbations from all eight major planets and the relativistic effects of orbital perihelion displacement were taken into account. For all orbits in the 3:1 and 5:2 MMRs an increase in the orbit eccentricities, which are sufficient for the bodies to approach Mars, was obtained. In the 2:1 MMR, a sufficient increase in the orbit eccentricities was not detected. An increase in orbit eccentricities in this resonance can occur due to the action of non-gravitational effects (NGEs). The action of the Yarkovsky effect can explain the exit of an asteroid with a size of 5 km from the 2:1 MMR over a period about 1 billion years or more. More than 2 billion years ago, there were dozens of bodies over 50 km in size in the 2:1 gap. To form the gap in the 2:1 resonance at the very beginning, the physical conditions in the asteroid belt had to be significantly different from the today ones. In particular, the intensity of the solar radiation in the early Solar system could be much higher as compared to the today one.


title: Efficient Gradient Approximation Method for Constrained Bilevel Optimization

abstract: Bilevel optimization has been developed for many machine learning tasks with large-scale and high-dimensional data. This paper considers a constrained bilevel optimization problem, where the lower-level optimization problem is convex with equality and inequality constraints and the upper-level optimization problem is non-convex. The overall objective function is non-convex and non-differentiable. To solve the problem, we develop a gradient-based approach, called gradient approximation method, which determines the descent direction by computing several representative gradients of the objective function inside a neighborhood of the current estimate. We show that the algorithm asymptotically converges to the set of Clarke stationary points, and demonstrate the efficacy of the algorithm by the experiments on hyperparameter optimization and meta-learning.


title: How Bad is Top-$K$ Recommendation under Competing Content Creators?

abstract: Content creators compete for exposure on recommendation platforms, and such strategic behavior leads to a dynamic shift over the content distribution. However, how the creators' competition impacts user welfare and how the relevance-driven recommendation influences the dynamics in the long run are still largely unknown. This work provides theoretical insights into these research questions. We model the creators' competition under the assumptions that: 1) the platform employs an innocuous top-$K$ recommendation policy; 2) user decisions follow the Random Utility model; 3) content creators compete for user engagement and, without knowing their utility function in hindsight, apply arbitrary no-regret learning algorithms to update their strategies. We study the user welfare guarantee through the lens of Price of Anarchy and show that the fraction of user welfare loss due to creator competition is always upper bounded by a small constant depending on $K$ and randomness in user decisions; we also prove the tightness of this bound. Our result discloses an intrinsic merit of the myopic approach to the recommendation, i.e., relevance-driven matching performs reasonably well in the long run, as long as users' decisions involve randomness and the platform provides reasonably many alternatives to its users.


title: DCA: Delayed Charging Attack on the Electric Shared Mobility System

abstract: An efficient operation of the electric shared mobility system (ESMS) relies heavily on seamless interconnections among shared electric vehicles (SEV), electric vehicle supply equipment (EVSE), and the grid. Nevertheless, this interconnectivity also makes the ESMS vulnerable to cyberattacks that may cause short-term breakdowns or long-term degradation of the ESMS. This study focuses on one such attack with long-lasting effects, the Delayed Charge Attack (DCA), that stealthily delays the charging service by exploiting the physical and communication vulnerabilities. To begin, we present the ESMS threat model by highlighting the assets, information flow, and access points. We next identify a linked sequence of vulnerabilities as a viable attack vector for launching DCA. Then, we detail the implementation of DCA, which can effectively bypass the detection in the SEV's battery management system and the cross-verification in the cloud environment. We test the DCA model against various Anomaly Detection (AD) algorithms by simulating the DCA dynamics in a Susceptible-Infectious-Removed-Susceptible process, where the EVSE can be compromised by the DCA or detected for repair. Using real-world taxi trip data and EVSE locations in New York City, the DCA model allows us to explore the long-term impacts and validate the system consequences. The results show that a 10-min delay results in 12-min longer queuing times and 8% more unfulfilled requests, leading to a 10.7% (\$311.7) weekly revenue loss per driver. With the AD algorithms, the weekly revenue loss remains at least 3.8% (\$111.8) with increased repair costs of \$36,000, suggesting the DCA's robustness against the AD.

categories: cs.CR cs.SY eess.SY math.DS math.OC physics.soc-ph


title: Measuring The Impact Of Programming Language Distribution

abstract: Current benchmarks for evaluating neural code models focus on only a small subset of programming languages, excluding many popular languages such as Go or Rust. To ameliorate this issue, we present the BabelCode framework for execution-based evaluation of any benchmark in any language. BabelCode enables new investigations into the qualitative performance of models' memory, runtime, and individual test case results. Additionally, we present a new code translation dataset called Translating Python Programming Puzzles (TP3) from the Python Programming Puzzles (Schuster et al. 2021) benchmark that involves translating expert-level python functions to any language. With both BabelCode and the TP3 benchmark, we investigate if balancing the distributions of 14 languages in a training dataset improves a large language model's performance on low-resource languages. Training a model on a balanced corpus results in, on average, 12.34% higher $pass@k$ across all tasks and languages compared to the baseline. We find that this strategy achieves 66.48% better $pass@k$ on low-resource languages at the cost of only a 12.94% decrease to high-resource languages. In our three translation tasks, this strategy yields, on average, 30.77% better low-resource $pass@k$ while having 19.58% worse high-resource $pass@k$.


title: Conic Sparsity: Estimation of Regression Parameters in Closed Convex Polyhedral Cones

abstract: Statistical problems often involve linear equality and inequality constraints on model parameters. Direct estimation of parameters restricted to general polyhedral cones, particularly when one is interested in estimating low dimensional features, may be challenging. We use a dual form parameterization to characterize parameter vectors restricted to lower dimensional faces of polyhedral cones and use the characterization to define a notion of 'sparsity' on such cones. We show that the proposed notion agrees with the usual notion of sparsity in the unrestricted case and prove the validity of the proposed definition as a measure of sparsity. The identifiable parameterization of the lower dimensional faces allows a generalization of popular spike-and-slab priors to a closed convex polyhedral cone. The prior measure utilizes the geometry of the cone by defining a Markov random field over the adjacency graph of the extreme rays of the cone. We describe an efficient way of computing the posterior of the parameters in the restricted case. We illustrate the usefulness of the proposed methodology for imposing linear equality and inequality constraints by using wearables data from the National Health and Nutrition Examination Survey (NHANES) actigraph study where the daily average activity profiles of participants exhibit patterns that seem to obey such constraints.


title: Interdicting attack plans with boundedly-rational players and multiple attackers: An adversarial risk analysis approach

abstract: Cybersecurity planning supports the selection of and implementation of security controls in resource-constrained settings to manage risk. Doing so requires considering adaptive adversaries with different levels of strategic sophistication in modeling efforts to support risk management. However, most models in the literature only consider rational or non-strategic adversaries. Therefore, we study how to inform defensive decision-making to mitigate the risk from boundedly rational players, with a particular focus on making integrated, interdependent planning decisions. To achieve this goal, we introduce a modeling framework for selecting a portfolio of security mitigations that interdict adversarial attack plans that uses a structured approach for risk analysis. Our approach adapts adversarial risk analysis and cognitive hierarchy theory to consider a maximum reliability path interdiction problem with a single defender and multiple attackers who have different goals and levels of strategic sophistication. Instead of enumerating all possible attacks and defenses, we introduce a solution technique based on integer programming and approximation algorithms to iteratively solve the defender's and attackers' problems. A case study illustrates the proposed models and provides insights into defensive planning.


title: SPARLING: Learning Latent Representations with Extremely Sparse Activations

abstract: Real-world processes often contain intermediate state that can be modeled as an extremely sparse tensor. We introduce Sparling, a technique that allows you to learn models with intermediate layers that match this state from only end-to-end labeled examples (i.e., no supervision on the intermediate state). Sparling uses a new kind of informational bottleneck that enforces levels of activation sparsity unachievable using other techniques. We find that extreme sparsity is necessary to achieve good intermediate state modeling. On our synthetic DigitCircle domain as well as the LaTeX-OCR and Audio-MNIST-Sequence domains, we are able to precisely localize the intermediate states up to feature permutation with > 90% accuracy, even though we only train end-to-end.


title: Construction of Hierarchically Semi-Separable matrix Representation using Adaptive Johnson-Lindenstrauss Sketching

abstract: We extend an adaptive partially matrix-free Hierarchically Semi-Separable (HSS) matrix construction algorithm by Gorman et al. [SIAM J. Sci. Comput. 41(5), 2019] which uses Gaussian sketching operators to a broader class of Johnson--Lindenstrauss (JL) sketching operators. We present theoretical work which justifies this extension. In particular, we extend the earlier concentration bounds to all JL sketching operators and examine this bound for specific classes of such operators including the original Gaussian sketching operators, subsampled randomized Hadamard transform (SRHT) and the sparse Johnson--Lindenstrauss transform (SJLT). We discuss the implementation details of applying SJLT efficiently and demonstrate experimentally that using SJLT instead of Gaussian sketching operators leads to 1.5--2.5x speedups of the HSS construction implementation in the STRUMPACK C++ library. The generalized algorithm allows users to select their own JL sketching operators with theoretical lower bounds on the size of the operators which may lead to faster run time with similar HSS construction accuracy.


title: A New Paradigm of Reservoir Computing Exploiting Hydrodynamics

abstract: Nonlinear waves have played a historical role in laying the foundations of the science of complexity. Recently, they have also allowed the development of a new reservoir computing paradigm: neuromorphic computing by waves. In these systems, the information transmission acts as the excitation of wave dynamics, whose evolution processes the information to perform complex tasks at low energy consumption. To enable nonlinear hydrodynamic waves to do computing, we designed the Aqua-Photonic-Advantaged Computing Machine by Artificial Neural Networks (Aqua-PACMANN), a system where wave propagation in shallow water is the leading physical phenomenon, and the presence of electronics can be reduced to a CCD camera in detection. We show the realization of an XNOR logic gate as proof of concept of the Aqua-PACMANN's architecture and pave the way to a new class of fluid dynamic neuromorphic computing.


title: Hybrid Multiparty Session Types -- Full Version

abstract: Multiparty session types (MPST) are a specification and verification framework for distributed message-passing systems. The communication protocol of the system is specified as a global type, from which a collection of local types (local process implementations) is obtained by endpoint projection. A global type is a single disciplining entity for the whole system, specified by one designer that has full knowledge of the communication protocol. On the other hand, distributed systems are often described in terms of their components: a different designer is in charge of providing a subprotocol for each component. The problem of modular specification of global protocols has been addressed in the literature, but the state of the art focuses only on dual input/output compatibility. Our work overcomes this limitation. We propose the first MPST theory of multiparty compositionality for distributed protocol specification that is semantics-preserving, allows the composition of two or more components, and retains full MPST expressiveness. We introduce hybrid types for describing subprotocols interacting with each other, define a novel compatibility relation, explicitly describe an algorithm for composing multiple subprotocols into a well-formed global type, and prove that compositionality preserves projection, thus retaining semantic guarantees, such as liveness and deadlock freedom. Finally, we test our work against real-world case studies and we smoothly extend our novel compatibility to MPST with delegation and explicit connections.


title: Sub-Bergman Hilbert spaces on the unit disk III

abstract: For a bounded analytic function $\varphi$ on the unit disk $\D$ with $\|\varphi\|_\infty\le1$ we consider the defect operators $D_\varphi$ and $D_{\overline\varphi}$ of the Toeplitz operators $T_\varphi$ and $T_{\overline\varphi}$, respectively, on the weighted Bergman space $A^2_\alpha$. The ranges of $D_\varphi$ and $D_{\overline\varphi}$, written as $H(\varphi)$ and $H(\overline\varphi)$ and equipped with appropriate inner products, are called sub-Bergman spaces. We prove the following three results in the paper: for $-1<\alpha\le0$ the space $H(\varphi)$ has a complete Nevanlinna-Pick kernel if and only if $\varphi$ is a M\"{o}bius map; for $\alpha>-1$ we have $H(\varphi)=H(\overline\varphi)=A^2_{\alpha-1}$ if and only if the defect operators $D_\varphi$ and $D_{\overline\varphi}$ are compact; and for $\alpha>-1$ we have $D^2_\varphi(A^2_\alpha)= D^2_{\overline\varphi}(A^2_\alpha)=A^2_{\alpha-2}$ if and only if $\varphi$ is a finite Blaschke product. In some sense our restrictions on $\alpha$ here are best possible.


title: Origin of Wake-Up Effect in Hafnia

abstract: We performed density functional theory (DFT) calculations on epitaxially strained hafnia. We demonstrate the stabilization of the ferroelectric ($Pca2_{1}$) phase from the antiferroelectric ($Pbcn$) in bulk hafnia in the presence of electric field. We found that the polar ($Pca2_{1}$) phase can be efficiently stabilized with an adequate choice of film orientation. We show that for a (010)-oriented Pbcn, the ferroelectric $Pca2_{1}$ phase can be reached with a relatively small electric field ($\mid\varepsilon\mid\geq 150$ KV/m). We thus provide a simple explanation to the experimental observation of polarization enhancement through electric field cycling, or wake-up effect, as a ferroelectric phase transition driven by electric field. We find, in contrast, that stress free pure hafnia does not become ferroelectric for any reasonable electric field. So we explain the wake up effect and stabilization of ferroelectric pure hafnia as coming from a combination of epitaxial strain under applied electric field perpendicular to the film. We find that strain (or doping) primarily destabilizes the baddeleyite structure, so that the antiferroelectric Pbcn and ferroelectric phases can form.