title: The Alfv\'enic nature of chromospheric swirls

abstract: We investigate the evolution and origin of small-scale chromospheric swirls by analyzing numerical simulations of the quiet solar atmosphere, using the radiative magnetohydrodynamic code CO$5$BOLD. We are interested in finding their relation with magnetic field perturbations and in the processes driving their evolution. For the analysis, the swirling strength criterion and its evolution equation are applied in order to identify vortical motions and to study their dynamics. We introduce a new criterion, the magnetic swirling strength, which allows us to recognize torsional perturbations in the magnetic field. We find a strong correlation between swirling strength and magnetic swirling strength, in particular in intense magnetic flux concentrations, which suggests a tight relation between vortical motions and torsional magnetic field perturbations. Furthermore, we find that swirls propagate upward with the local Alfv\'en speed as unidirectional swirls, in the form of pulses, driven by magnetic tension forces alone. In the photosphere and low chromosphere, the rotation of the plasma co-occurs with a twist in the upwardly directed magnetic field that is in the opposite direction of the plasma flow. All together, these are characteristics of torsional Alfv\'en waves. We also find indications of an imbalance between the hydrodynamic and magnetohydrodynamic baroclinic effects being at the origin of the swirls. At the base of the chromosphere, we find a net upwardly directed Poynting flux, which is mostly associated with large and complex swirling structures that we interpret as the superposition of various small-scale vortices. We conclude that the ubiquitous swirling events observed in simulations are tightly correlated with perturbations of the magnetic field. At photospheric and chromospheric levels, they form Alfv\'en pulses that propagate upward and may contribute to chromospheric heating.


title: Tight Competitive Analyses of Online Car-sharing Problems

abstract: The car-sharing problem, proposed by Luo, Erlebach and Xu in 2018, mainly focuses on an online model in which there are two locations: 0 and 1, and $k$ total cars. Each request which specifies its pick-up time and pick-up location (among 0 and 1, and the other is the drop-off location) is released in each stage a fixed amount of time before its specified start (i.e. pick-up) time. The time between the booking (i.e. released) time and the start time is enough to move empty cars between 0 and 1 for relocation if they are not used in that stage. The model, called $k$S2L-F, assumes that requests in each stage arrive sequentially regardless of the same booking time and the decision (accept or reject) must be made immediately. The goal is to accept as many requests as possible. In spite of only two locations, the analysis does not seem easy and the (tight) competitive ratio (CR) is only known to be 2.0 for $k=2$ and 1.5 for a restricted value of $k$, i.e., a multiple of three. In this paper, we remove all the holes of unknown CR's; namely we prove that the CR is $\frac{2k}{k + \lfloor k/3 \rfloor}$ for all $k\geq 2$. Furthermore, if the algorithm can delay its decision until all requests have come in each stage, the CR is improved to roughly 4/3. We can take this advantage even further, precisely we can achieve a CR of $\frac{2+R}{3}$ if the number of requests in each stage is at most $Rk$, $1 \leq R \leq 2$, where we do not have to know the value of $R$ in advance. Finally we demonstrate that randomization also helps to get (slightly) better CR's.


title: Information Maximization Clustering via Multi-View Self-Labelling

abstract: Image clustering is a particularly challenging computer vision task, which aims to generate annotations without human supervision. Recent advances focus on the use of self-supervised learning strategies in image clustering, by first learning valuable semantics and then clustering the image representations. These multiple-phase algorithms, however, increase the computational time and their final performance is reliant on the first stage. By extending the self-supervised approach, we propose a novel single-phase clustering method that simultaneously learns meaningful representations and assigns the corresponding annotations. This is achieved by integrating a discrete representation into the self-supervised paradigm through a classifier net. Specifically, the proposed clustering objective employs mutual information, and maximizes the dependency between the integrated discrete representation and a discrete probability distribution. The discrete probability distribution is derived though the self-supervised process by comparing the learnt latent representation with a set of trainable prototypes. To enhance the learning performance of the classifier, we jointly apply the mutual information across multi-crop views. Our empirical results show that the proposed framework outperforms state-of-the-art techniques with the average accuracy of 89.1% and 49.0%, respectively, on CIFAR-10 and CIFAR-100/20 datasets. Finally, the proposed method also demonstrates attractive robustness to parameter settings, making it ready to be applicable to other datasets.


title: A hole spin qubit in a fin field-effect transistor above 4 kelvin

abstract: The greatest challenge in quantum computing is achieving scalability. Classical computing previously faced a scalability issue, solved with silicon chips hosting billions of fin field-effect transistors (FinFETs). These FinFET devices are small enough for quantum applications: at low temperatures, an electron or hole trapped under the gate serves as a spin qubit. Such an approach potentially allows the quantum hardware and its classical control electronics to be integrated on the same chip. However, this requires qubit operation at temperatures above 1K, where the cooling overcomes heat dissipation. Here, we show that silicon FinFETs can host spin qubits operating above 4K. We achieve fast electrical control of hole spins with driving frequencies up to 150MHz, single-qubit gate fidelities at the fault-tolerance threshold, and a Rabi oscillation quality factor greater than 87. Our devices feature both industry compatibility and quality, and are fabricated in a flexible and agile way that should accelerate further development.


title: Test case generation for agent-based models: A systematic literature review

abstract: Agent-based models play an important role in simulating complex emergent phenomena and supporting critical decisions. In this context, a software fault may result in poorly informed decisions that lead to disastrous consequences. The ability to rigorously test these models is therefore essential. In this systematic literature review, we answer five research questions related to the key aspects of test case generation in agent-based models: What are the information artifacts used to generate tests? How are these tests generated? How is a verdict assigned to a generated test? How is the adequacy of a generated test suite measured? What level of abstraction of an agent-based model is targeted by a generated test? Our results show that whilst the majority of techniques are effective for testing functional requirements at the agent and integration levels of abstraction, there are comparatively few techniques capable of testing society-level behaviour. Additionally, we identify a need for more thorough evaluation using realistic case studies that feature challenging properties associated with a typical agent-based model.


title: PatchNet -- Short-range Template Matching for Efficient Video Processing

abstract: Object recognition is a fundamental problem in many video processing tasks, accurately locating seen objects at low computation cost paves the way for on-device video recognition. We propose PatchNet, an efficient convolutional neural network to match objects in adjacent video frames. It learns the patchwise correlation features instead of pixel features. PatchNet is very compact, running at just 58MFLOPs, $5\times$ simpler than MobileNetV2. We demonstrate its application on two tasks, video object detection and visual object tracking. On ImageNet VID, PatchNet reduces the flops of R-FCN ResNet-101 by 5x and EfficientDet-D0 by 3.4x with less than 1% mAP loss. On OTB2015, PatchNet reduces SiamFC and SiamRPN by 2.5x with no accuracy loss. Experiments on Jetson Nano further demonstrate 2.8x to 4.3x speed-ups associated with flops reduction. Code is open sourced at https://github.com/RalphMao/PatchNet.


title: ACTION-Net: Multipath Excitation for Action Recognition

abstract: Spatial-temporal, channel-wise, and motion patterns are three complementary and crucial types of information for video action recognition. Conventional 2D CNNs are computationally cheap but cannot catch temporal relationships; 3D CNNs can achieve good performance but are computationally intensive. In this work, we tackle this dilemma by designing a generic and effective module that can be embedded into 2D CNNs. To this end, we propose a spAtio-temporal, Channel and moTion excitatION (ACTION) module consisting of three paths: Spatio-Temporal Excitation (STE) path, Channel Excitation (CE) path, and Motion Excitation (ME) path. The STE path employs one channel 3D convolution to characterize spatio-temporal representation. The CE path adaptively recalibrates channel-wise feature responses by explicitly modeling interdependencies between channels in terms of the temporal aspect. The ME path calculates feature-level temporal differences, which is then utilized to excite motion-sensitive channels. We equip 2D CNNs with the proposed ACTION module to form a simple yet effective ACTION-Net with very limited extra computational cost. ACTION-Net is demonstrated by consistently outperforming 2D CNN counterparts on three backbones (i.e., ResNet-50, MobileNet V2 and BNInception) employing three datasets (i.e., Something-Something V2, Jester, and EgoGesture). Codes are available at \url{https://github.com/V-Sense/ACTION-Net}.


title: Compression of Solar Spectroscopic Observations: a Case Study of Mg II k Spectral Line Profiles Observed by NASA's IRIS Satellite

abstract: In this study we extract the deep features and investigate the compression of the Mg II k spectral line profiles observed in quiet Sun regions by NASA's IRIS satellite. The data set of line profiles used for the analysis was obtained on April 20th, 2020, at the center of the solar disc, and contains almost 300,000 individual Mg II k line profiles after data cleaning. The data are separated into train and test subsets. The train subset was used to train the autoencoder of the varying embedding layer size. The early stopping criterion was implemented on the test subset to prevent the model from overfitting. Our results indicate that it is possible to compress the spectral line profiles more than 27 times (which corresponds to the reduction of the data dimensionality from 110 to 4) while having a 4 DN average reconstruction error, which is comparable to the variations in the line continuum. The mean squared error and the reconstruction error of even statistical moments sharply decrease when the dimensionality of the embedding layer increases from 1 to 4 and almost stop decreasing for higher numbers. The observed occasional improvements in training for values higher than 4 indicate that a better compact embedding may potentially be obtained if other training strategies and longer training times are used. The features learned for the critical four-dimensional case can be interpreted. In particular, three of these four features mainly control the line width, line asymmetry, and line dip formation respectively. The presented results are the first attempt to obtain a compact embedding for spectroscopic line profiles and confirm the value of this approach, in particular for feature extraction, data compression, and denoising.


title: Quantum communication through a partially reflecting accelerating mirror

abstract: Motivated by the fact that the null-shell of a collapsing black hole can be described by a perfectly reflecting accelerating mirror, we investigate an extension of this model to mirror semi-transparency and derive a general implicit expression for the corresponding Bogoliubov coefficients. Then, we turn this into an explicit analytical form by focusing on mirrors that are accelerated via an impulsive force. From the so-obtained Bogoliubov coefficients we derive the particle production. Finally, we realize the field coming from left-past spacetime region, passing through the semitransparent moving mirror and ending up to right-future spacetime region as undergoing the action of a Gaussian quantum channel. We study the transmission and noisy generation properties of this channel, relating them to the Bogoliubov coefficients of the mirror's motion, through which we evaluate capacities in transmitting classical and quantum information.


title: On the edge dimension and fractional edge dimension of graphs

abstract: Let $G$ be a graph with vertex set $V(G)$ and edge set $E(G)$, and let $d(u,w)$ denote the length of a $u-w$ geodesic in $G$. For any $v\in V(G)$ and $e=xy\in E(G)$, let $d(e,v)=\min\{d(x,v),d(y,v)\}$. For distinct $e_1, e_2\in E(G)$, let $R\{e_1,e_2\}=\{z\in V(G):d(z,e_1)\neq d(z,e_2)\}$. Kelenc et al. [Discrete Appl. Math. 251 (2018) 204-220] introduced the edge dimension of a graph: A vertex subset $S\subseteq V(G)$ is an edge resolving set of $G$ if $|S\cap R\{e_1,e_2\}|\ge 1$ for any distinct $e_1, e_2\in E(G)$, and the edge dimension $edim(G)$ of $G$ is the minimum cardinality among all edge resolving sets of $G$. For a real-valued function $g$ defined on $V(G)$ and for $U\subseteq V(G)$, let $g(U)=\sum_{s\in U}g(s)$. Then $g:V(G)\rightarrow[0,1]$ is an edge resolving function of $G$ if $g(R\{e_1,e_2\})\ge1$ for any distinct $e_1,e_2\in E(G)$. The fractional edge dimension $edim_f(G)$ of $G$ is $\min\{g(V(G)):g\mbox{ is an edge resolving function of }G\}$. Note that $edim_f(G)$ reduces to $edim(G)$ if the codomain of edge resolving functions is restricted to $\{0,1\}$. We introduce and study fractional edge dimension and obtain some general results on the edge dimension of graphs. We show that there exist two non-isomorphic graphs on the same vertex set with the same edge metric coordinates. We construct two graphs $G$ and $H$ such that $H \subset G$ and both $edim(H)-edim(G)$ and $edim_f(H)-edim_f(G)$ can be arbitrarily large. We show that a graph $G$ with $edim(G)=2$ cannot have $K_5$ or $K_{3,3}$ as a subgraph, and we construct a non-planar graph $H$ satisfying $edim(H)=2$. It is easy to see that, for any connected graph $G$ of order $n\ge3$, $1\le edim_f(G) \le \frac{n}{2}$; we characterize graphs $G$ satisfying $edim_f(G)=1$ and examine some graph classes satisfying $edim_f(G)=\frac{n}{2}$. We also determine the fractional edge dimension for some classes of graphs.


title: A measure differential inclusion involving time-dependent maximal monotone operators

abstract: We establish in this paper, an existence and uniqueness result of right continuous with bounded variation solution for a perturbed differential inclusion governed by time-dependent maximal monotone operators.


title: Interface spaces based on physics for multiscale mixed methods applied to flows in fractured-like porous media

abstract: It is well known that domain-decomposition-based multiscale mixed methods rely on interface spaces, defined on the skeleton of the decomposition, to connect the solution among the non-overlapping subdomains. Usual spaces, such as polynomial-based ones, cannot properly represent high-contrast channelized features such as fractures (high permeability) and barriers (low permeability) for flows in heterogeneous porous media. We propose here new interface spaces, which are based on physics, to deal with permeability fields in the simultaneous presence of fractures and barriers, accommodated respectively, by the pressure and flux spaces. Existing multiscale methods based on mixed formulations can take advantage of the proposed interface spaces, however, in order to present and test our results, we use the newly developed Multiscale Robin Coupled Method (MRCM) [Guiraldello, et al., J. Comput. Phys., 355 (2018) pp. 1-21], which generalizes most well-known multiscale mixed methods, and allows for the independent choice of the pressure and flux interface spaces. An adaptive version of the MRCM [Rocha, et al., J. Comput. Phys., 409 (2020), 109316] is considered that automatically selects the physics-based pressure space for fractured structures and the physics-based flux space for regions with barriers, resulting in a procedure with unprecedented accuracy. The features of the proposed approach are investigated through several numerical simulations of single-phase and two-phase flows, in different heterogeneous porous media. The adaptive MRCM combined with the interface spaces based on physics provides promising results for challenging problems with the simultaneous presence of fractures and barriers.


title: On eigenvalues of the Brownian sheet matrix

abstract: We derive a system of stochastic partial differential equations satisfied by the eigenvalues of the symmetric matrix whose entries are the Brownian sheets. We prove that the sequence $\left\{L_{d}(s,t), (s,t)\in[0,S]\times [0,T]\right\}_{d\in\mathbb N}$ of empirical spectral measures of the rescaled matrices is tight on $C([0,S]\times [0,T], \mathcal P(\mathbb R))$ and hence is convergent as $d$ goes to infinity by Wigner's semicircle law. We also obtain PDEs which are satisfied by the high-dimensional limiting measure.


title: Offset-free Model Predictive Control: A Ball Catching Application with a Spherical Soft Robotic Arm

abstract: This paper presents an offset-free model predictive controller for fast and accurate control of a spherical soft robotic arm. In this control scheme, a linear model is combined with an online disturbance estimation technique to systematically compensate model deviations. Dynamic effects such as material relaxation resulting from the use of soft materials can be addressed to achieve offset-free tracking. The tracking error can be reduced by 35% when compared to a standard model predictive controller without a disturbance compensation scheme. The improved tracking performance enables the realization of a ball catching application, where the spherical soft robotic arm can catch a ball thrown by a human.


title: Differentiating densities on smooth manifolds

abstract: Lebesgue integration of derivatives of strongly-oscillatory functions is a recurring challenge in computational science and engineering. Integration by parts is an effective remedy for huge computational costs associated with Monte Carlo integration schemes. In case of Lebesgue integrals over a smooth manifold, however, integration by parts gives rise to a derivative of the density implied by charts describing the domain manifold. This paper focuses on the computation of that derivative, which we call the density gradient function, on general smooth manifolds. We analytically derive formulas for the density gradient and present examples of manifolds determined by popular differential equation-driven systems. We highlight the significance of the density gradient by demonstrating a numerical example of Monte Carlo integration involving oscillatory integrands.


title: An elementary proof for dynamical scaling for certain fractional non-homogeneous Poisson processes

abstract: Dynamical scaling is an asymptotic property typical for the dynamics of first-order phase transitions in physical systems and related to self-similarity. Based on the integral-representation for the marginal probabilities of a fractional non-homogeneous Poisson process introduced by Leonenko et al. (2017) and generalising the standard fractional Poisson process, we prove the dynamical scaling under fairly mild conditions. Our result also includes the special case of the standard fractional Poisson process.


title: Value of information from vibration-based structural health monitoring extracted via Bayesian model updating

abstract: Quantifying the value of the information extracted from a structural health monitoring (SHM) system is an important step towards convincing decision makers to implement these systems. We quantify this value by adaptation of the Bayesian decision analysis framework. In contrast to previous works, we model in detail the entire process of data generation to processing, model updating and reliability calculation, and investigate it on a deteriorating bridge system. The framework assumes that dynamic response data are obtained in a sequential fashion from deployed accelerometers, subsequently processed by an output-only operational modal analysis scheme for identifying the system's modal characteristics. We employ a classical Bayesian model updating methodology to sequentially learn the deterioration and estimate the structural damage evolution over time. This leads to sequential updating of the structural reliability, which constitutes the basis for a preposterior Bayesian decision analysis. Alternative actions are defined and a heuristic-based approach is employed for the life-cycle optimization. By solving the preposterior Bayesian decision analysis, one is able to quantify the benefit of the availability of long-term SHM vibrational data. Numerical investigations show that this framework can provide quantitative measures on the optimality of an SHM system in a specific decision context.


title: On the Analyticity of Critical Points of the Generalized Integral Menger Curvature in the Hilbert Case

abstract: We prove the analyticity of smooth critical points for generalized integral Menger curvature energies $\mathrm{intM}^{(p,2)}$, with $p \in (\tfrac 73, \tfrac 83)$, subject to a fixed length constraint. This implies, together with already well-known regularity results, that finite-energy, critical $C^1$-curves $\gamma: \mathbb{R}/\mathbb{Z} \to \mathbb{R}^n$ of generalized integral Menger curvature $\mathrm{intM}^{(p,2)}$ subject to a fixed length constraint are not only $C^\infty$ but also analytic. Our approach is inspired by analyticity results on critical points for O'Hara's knot energies based on Cauchy's method of majorants and a decomposition of the first variation. The main new idea is an additional iteration in the recursive estimate of the derivatives to obtain a sufficient difference in the order of regularity.


title: When scale is surplus

abstract: We study a long-recognised but under-appreciated symmetry called "dynamical similarity" and illustrate its relevance to many important conceptual problems in fundamental physics. Dynamical similarities are general transformations of a system where the unit of Hamilton's principal function is rescaled, and therefore represent a kind of dynamical scaling symmetry with formal properties that differ from many standard symmetries. To study this symmetry, we develop a general framework for symmetries that distinguishes the observable and surplus structures of a theory by using the minimal freely specifiable initial data for the theory that is necessary to achieve empirical adequacy. This framework is then applied to well-studied examples including Galilean invariance and the symmetries of the Kepler problem. We find that our framework gives a precise dynamical criterion for identifying the observables of those systems, and that those observables agree with epistemic expectations. We then apply our framework to dynamical similarity. First we give a general definition of dynamical similarity. Then we show, with the help of some previous results, how the dynamics of our observables leads to singularity resolution and the emergence of an arrow of time in cosmology.


title: Flow structure and loads over inclined cylindrical rodlike particles and fibers

abstract: The flow past a fixed finite-length circular cylinder, the axis of which makes a nonzero angle with the incoming stream, is studied through fully-resolved simulations, from creeping-flow conditions to strongly inertial regimes. The investigation focuses on the way the body aspect ratio $\chi$ (defined as as the length-to-diameter ratio), the inclination angle $\theta$ with respect to the incoming flow and the Reynolds number $\text{Re}$ (based on the cylinder diameter) affect the flow structure past the body and therefore the hydrodynamic loads acting on it. The configuration $\theta=0^\circ$ (where the cylinder is aligned with the flow) is considered first, from creeping-flow conditions up to $\text{Re}=400$, with aspect ratios up to $20$ ($10$) for $\text{Re}\leq10$ ($\text{Re}\geq10$). In the low-to-moderate Reynolds number regime ($\text{Re}\lesssim5$), influence or the aspect ratio, inclination (from $0^{\circ}$ to $30^{\circ}$), and inertial effects is examined by comparing numerical results for the axial and transverse force components and the spanwise torque with theoretical predictions based on the slender-body approximation, possibly incorporating finite-Reynolds-number corrections. Semiempirical models based on these predictions and incorporating finite-length and inertial corrections extracted from the numerical data are derived. For large enough Reynolds numbers ($\text{Re}\gtrsim10^2$), separation takes place along the upstream part of the lateral surface of the cylinder, deeply influencing the surface stress distribution. Numerical results are used to build empirical models for the force components and the torque, valid for moderately inclined cylinders ($|\theta|\lesssim30^\circ$) of arbitrary aspect ratio up to $\text{Re}\approx300$ and matching those obtained at low-to-moderate Reynolds number.


title: LayerPCM: An implicit scheme for dielectric screening from layered substrates

abstract: We present LayerPCM, an extension of the polarizable-continuum model coupled to real-time time-dependent density-functional theory for an efficient and accurate description of the electrostatic interactions between molecules and multilayered dielectric substrates on which they are physisorbed. The former are modelled quantum-mechanically, while the latter are treated as polarizable continua characterized by their dielectric constants. The proposed approach is purposely designed to simulate complex hybrid heterostructures, with nano-engineered substrates including a stack of anisotropic layers. LayerPCM is suitable to describe the polarization-induced renormalization of frontier energy levels of the adsorbates in the static regime. Moreover, it can be reliably applied to simulating laser-induced ultrafast dynamics of the molecules through the inclusion of electric fields generated by Fresnel-reflection at the substrate. Depending on the complexity of the underlying layer structure, such reflected fields can assume non-trivial shapes and profoundly affect the dynamics of the photo-excited charge carriers in the molecule. In particular, the interaction with the substrate can give rise to strong delayed fields which lead to interference effects resembling those of multi-pulse-based spectroscopy. The robustness of the implementation and the above-mentioned features are demonstrated with a number of examples, ranging from intuitive models to realistic systems.


title: From coarse wall measurements to turbulent velocity fields through deep learning

abstract: This work evaluates the applicability of super-resolution generative adversarial networks (SRGANs) as a methodology for the reconstruction of turbulent-flow quantities from coarse wall measurements. The method is applied both for the resolution enhancement of wall fields and the estimation of wall-parallel velocity fields from coarse wall measurements of shear stress and pressure. The analysis has been carried out with a database of a turbulent open-channel flow with friction Reynolds number $Re_{\tau}=180$ generated through direct numerical simulation. Coarse wall measurements have been generated with three different downsampling factors $f_d=[4,8,16]$ from the high-resolution fields, and wall-parallel velocity fields have been reconstructed at four inner-scaled wall-normal distances $y^+=[15,30,50,100]$. We first show that SRGAN can be used to enhance the resolution of coarse wall measurements. If compared with direct reconstruction from the sole coarse wall measurements, SRGAN provides better instantaneous reconstructions, both in terms of mean-squared error and spectral-fractional error. Even though lower resolutions in the input wall data make it more challenging to achieve highly accurate predictions, the proposed SRGAN-based network yields very good reconstruction results. Furthermore, it is shown that even for the most challenging cases the SRGAN is capable of capturing the large-scale structures that populate the flow. The proposed novel methodology has great potential for closed-loop control applications relying on non-intrusive sensing.


title: A Neural Network for Semigroups

abstract: Tasks like image reconstruction in computer vision, matrix completion in recommender systems and link prediction in graph theory, are well studied in machine learning literature. In this work, we apply a denoising autoencoder-based neural network architecture to the task of completing partial multiplication (Cayley) tables of finite semigroups. We suggest a novel loss function for that task based on the algebraic nature of the semigroup data. We also provide a software package for conducting experiments similar to those carried out in this work. Our experiments showed that with only about 10% of the available data, it is possible to build a model capable of reconstructing a full Cayley from only half of it in about 80% of cases.


title: Toward a global model of the interactions in low-lying states of methyl cyanide: rotational and rovibrational spectroscopy of the $v_4 = 1$ state and tentative interstellar detection of the $v_4 = v_8 = 1$ state in Sgr B2(N)

abstract: New and existing rotational spectra of methyl cyanide were analyzed to extend the global model of low-lying vibrational states and their interactions to $v_4=1$ at 920 cm$^{-1}$. The rotational spectra cover large portions of the 36$-$1439 GHz region and reach quantum numbers $J$ and $K$ of 79 and 16, respectively. Information on the $K$ level structure of CH$_3$CN is obtained from IR spectra. A spectrum of $2\nu_8$ around 717 cm$^{-1}$, analyzed in our previous study, covered also the $\nu_4$ band. The assignments in this band cover 880$-$952 cm$^{-1}$, attaining quantum numbers $J$ and $K$ of 61 and 13, respectively. The most important interaction of $v_4=1$ appears to be with $v_8=3$, $\Delta K=0$, $\Delta l=+3$, a previously characterized anharmonic resonance. We report new analyses of interactions with $\Delta K=-2$ and $\Delta l=+1$, with $\Delta K=-4$ and $\Delta l=-1$, and with $\Delta K=-6$ and $\Delta l=-3$; these four types of interactions connect all $l$ substates of $v_8=3$ in energy to $v_4=1$. A known $\Delta K=-2$, $\Delta l=+1$ interaction with $v_7=1$ was also analyzed, and investigations of the $\Delta K=+1$, $\Delta l=-2$ and $\Delta K=+3$, $\Delta l=0$ resonances with $v_8=2$ were improved, as were interactions between successive states with $v_8\le 3$, mainly through new $v_8\le 2$ rotational data. A preliminary single state analysis of the $v_4=v_8=1$ state was carried out based on rotational transition frequencies and on $\nu_4+\nu_8-\nu_8$ hot band data. A considerable fraction of the $K$ levels was reproduced within uncertainties in its entirety or in part, despite obvious widespread perturbations in $v_4=v_8=1$. We detect rotational transitions of methyl cyanide from within all vibrational states up to $v_4=1$ and $v_4=v_8=1$ tentatively toward the hot molecular core of Sagittarius B2(N) employing the Atacama Large Millimeter Array.

categories: astro-ph.GA astro-ph.IM astro-ph.SR physics.chem-ph


title: Signal Representations for Synthesizing Audio Textures with Generative Adversarial Networks

abstract: Generative Adversarial Networks (GANs) currently achieve the state-of-the-art sound synthesis quality for pitched musical instruments using a 2-channel spectrogram representation consisting of log magnitude and instantaneous frequency (the "IFSpectrogram"). Many other synthesis systems use representations derived from the magnitude spectra, and then depend on a backend component to invert the output magnitude spectrograms that generally result in audible artefacts associated with the inversion process. However, for signals that have closely-spaced frequency components such as non-pitched and other noisy sounds, training the GAN on the 2-channel IFSpectrogram representation offers no advantage over the magnitude spectra based representations. In this paper, we propose that training GANs on single-channel magnitude spectra, and using the Phase Gradient Heap Integration (PGHI) inversion algorithm is a better comprehensive approach for audio synthesis modeling of diverse signals that include pitched, non-pitched, and dynamically complex sounds. We show that this method produces higher-quality output for wideband and noisy sounds, such as pops and chirps, compared to using the IFSpectrogram. Furthermore, the sound quality for pitched sounds is comparable to using the IFSpectrogram, even while using a simpler representation with half the memory requirements.


title: Gibbs measures of Potts model on Cayley trees: a survey and applications

abstract: In this paper we give a systematic review of the theory of Gibbs measures of Potts model on Cayley trees (developed since 2013) and discuss many applications of the Potts model to real world situations: mainly biology, physics, and some examples of alloy behavior, cell sorting, financial engineering, flocking birds, flowing foams, image segmentation, medicine, sociology etc.


title: Temporal Logic for Social Networks

abstract: This paper introduces a logic with a class of social network models that is based on standard Linear Temporal Logic (LTL), leveraging the power of existing model checkers for the analysis of social networks. We provide a short literature overview, and then define our logic and its axiomatization, present some simple motivational examples of both models and formulas, and show its soundness and completeness via a translation into propositional formulas. Lastly, we briefly discuss model checking and time complexity analysis.


title: Short minimal codes and covering codes via strong blocking sets in projective spaces

abstract: Minimal linear codes are in one-to-one correspondence with special types of blocking sets of projective spaces over a finite field, which are called strong or cutting blocking sets. In this paper we prove an upper bound on the minimal length of minimal codes of dimension $k$ over the $q$-element Galois field which is linear in both $q$ and $k$, hence improve the previous superlinear bounds. This result determines the minimal length up to a small constant factor. We also improve the lower and upper bounds on the size of so called higgledy-piggledy line sets in projective spaces and apply these results to present improved bounds on the size of covering codes and saturating sets in projective spaces as well. The contributions rely on geometric and probabilistic arguments.


title: Despeckling Polarimetric SAR Data Using a Multi-Stream Complex-Valued Fully Convolutional Network

abstract: A Polarimetric Synthetic Aperture Radar (PolSAR) sensor is able to collect images in different polarization states, making it a rich source of information for target characterization. PolSAR images are inherently affected by speckle. Therefore, before deriving ad hoc products from the data, the polarimetric covariance matrix needs to be estimated by reducing speckle. In recent years, deep learning based despeckling methods have started to evolve from single channel SAR images to PolSAR images. To this aim, deep learning based approaches separate the real and imaginary components of the complex-valued covariance matrix and use them as independent channels in a standard convolutional neural networks. However, this approach neglects the mathematical relationship that exists between the real and imaginary components, resulting in sub-optimal output. Here, we propose a multi-stream complex-valued fully convolutional network to reduce speckle and effectively estimate the PolSAR covariance matrix. To evaluate the performance of CV-deSpeckNet, we used Sentinel-1 dual polarimetric SAR images to compare against its real-valued counterpart, that separates the real and imaginary parts of the complex covariance matrix. CV-deSpeckNet was also compared against the state of the art PolSAR despeckling methods. The results show CV-deSpeckNet was able to be trained with a fewer number of samples, has a higher generalization capability and resulted in a higher accuracy than its real-valued counterpart and state-of-the-art PolSAR despeckling methods. These results showcase the potential of complex-valued deep learning for PolSAR despeckling.


title: Empowering Visual Internet-of-Things Mashups with Self-Healing Capabilities

abstract: Internet-of-Things (IoT) systems have spread among different application domains, from home automation to industrial manufacturing processes. The rushed development by competing vendors to meet the market demand of IoT solutions, the lack of interoperability standards, and the overall lack of a defined set of best practices have resulted in a highly complex, heterogeneous, and frangible ecosystem. Several works have been pushing towards visual programming solutions to abstract the underlying complexity and help humans reason about it. As these solutions begin to meet widespread adoption, their building blocks usually do not consider reliability issues. Node-RED, being one of the most popular tools, also lacks such mechanisms, either built-in or via extensions. In this work we present SHEN (Self-Healing Extensions for Node-RED) which provides 17 nodes that collectively enable the implementation of self-healing strategies within this visual framework. We proceed to demonstrate the feasibility and effectiveness of the approach using real devices and fault injection techniques.


title: Inclusive production of a heavy-light dijet system in hybrid high-energy and collinear factorization

abstract: We propose the study of the inclusive hadroproduction of a heavy-flavored jet in association with a light jet, as a probe channel of strong interactions at high energies. We build up a hybrid factorization that encodes genuine high-energy effects, provided by a partial next-to-leading BFKL resummation, inside the standard collinear structure of the cross section. We present a detailed analysis of different distributions, shaped on kinematic ranges typical of experimental analyses at the Large Hadron Collider, and differential in rapidity, azimuthal angle and transverse momentum. The fair stability that these distributions exhibit under higher-order corrections motivates our interest toward future studies. Here, the hybrid factorization could help to deepen our understanding of heavy-flavor physics in wider kinematic ranges, like the ones accessible at the Electron-Ion Collider.


title: An extensible equality checking algorithm for dependent type theories

abstract: We present a general and user-extensible equality checking algorithm that is applicable to a large class of type theories. The algorithm has a type-directed phase for applying extensionality rules and a normalization phase based on computation rules, where both kinds of rules are defined using the type-theoretic concept of object-invertible rules. We also give sufficient syntactic criteria for recognizing such rules, as well as a simple pattern-matching algorithm for applying them. A third component of the algorithm is a suitable notion of principal arguments, which determines a notion of normal form. By varying these, we obtain known notions, such as weak head-normal and strong normal forms. We prove that our algorithm is sound. We implemented it in the Andromeda 2 proof assistant, which supports user-definable type theories. The user need only provide the equality rules they wish to use, which the algorithm automatically classifies as computation or extensionality rules, and select appropriate principal arguments.


title: Geometry-induced enhancement factor improvement in covered-gold-nanorod-dimer antennas

abstract: Illuminated gapped-gold-nanorod dimers hold surface plasmon polaritons (SPPs) that can be engineered, by an appropriate choice of geometrical parameters, to enhance the electromagnetic field at the gap, allowing applications in molecular detection via surface-enhanced Raman spectroscopy (SERS). Envisioning hybrid devices in which the SERS spectroscopy of molecules in the gap is complemented by electrical measurements, it arises the question of designing efficient geometries to contact the nanorods without decreasing the enhancement factor (EF) of the nanoantenna, i.e., the figure of merit for SERS spectroscopy. Within this framework we theoretically study the feasibility to fabricate designs based on covering with gold the far-from-the-gap areas of the dimer. We show that by tuning the geometrical parameters of the designs these systems can reach enhancement factors larger than the best achieved in the uncovered dimer: this supremacy survives even in the presence of dimer asymmetries and vacancies at the interfaces between the nanorods and the covering layers. Our results show that geometrical modifications away from the gap can improve the optical response at the gap, thus enabling the use of these devices both for hybrid and optical applications.

categories: physics.optics physics.chem-ph physics.comp-ph physics.ins-det


title: Boolean Hierarchical Tucker Networks on Quantum Annealers

abstract: Quantum annealing is an emerging technology with the potential to solve some of the computational challenges that remain unresolved as we approach an era beyond Moore's Law. In this work, we investigate the capabilities of the quantum annealers of D-Wave Systems, Inc., for computing a certain type of Boolean tensor decomposition called Boolean Hierarchical Tucker Network (BHTN). Boolean tensor decomposition problems ask for finding a decomposition of a high-dimensional tensor with categorical, [true, false], values, as a product of smaller Boolean core tensors. As the BHTN decompositions are usually not exact, we aim to approximate an input high-dimensional tensor by a product of lower-dimensional tensors such that the difference between both is minimized in some norm. We show that BHTN can be calculated as a sequence of optimization problems suitable for the D-Wave 2000Q quantum annealer. Although current technology is still fairly restricted in the problems they can address, we show that a complex problem such as BHTN can be solved efficiently and accurately.


title: Super-Macdonald polynomials: Orthogonality and Hilbert space interpretation

abstract: The super-Macdonald polynomials, introduced by Sergeev and Veselov, generalise the Macdonald polynomials to (arbitrary numbers of) two kinds of variables, and they are eigenfunctions of the deformed Macdonald-Ruijsenaars operators introduced by the same authors. We introduce a Hermitian form on the algebra spanned by the super-Macdonald polynomials, prove their orthogonality, compute their (quadratic) norms explicitly, and establish a corresponding Hilbert space interpretation of the super-Macdonald polynomials and deformed Macdonald-Ruijsenaars operators. This allows for a quantum mechanical interpretation of the models defined by the deformed Macdonald-Ruijsenaars operators. Motivated by recent results in the nonrelativistic ($q\to 1$) case, we propose that these models describe the particles and anti-particles of an underlying relativistic quantum field theory, thus providing a natural generalisation of the trigonometric Ruijsenaars model.


title: Bears with Hats and Independence Polynomials

abstract: Consider the following hat guessing game. A bear sits on each vertex of a graph $G$, and a demon puts on each bear a hat colored by one of $h$ colors. Each bear sees only the hat colors of his neighbors. Based on this information only, each bear has to guess $g$ colors and he guesses correctly if his hat color is included in his guesses. The bears win if at least one bear guesses correctly for any hat arrangement. We introduce a new parameter - fractional hat chromatic number $\hat{\mu}$, arising from the hat guessing game. The parameter $\hat{\mu}$ is related to the hat chromatic number which has been studied before. We present a surprising connection between the hat guessing game and the independence polynomial of graphs. This connection allows us to compute the fractional hat chromatic number of chordal graphs in polynomial time, to bound fractional hat chromatic number by a function of maximum degree of $G$, and to compute the exact value of $\hat{\mu}$ of cliques, paths, and cycles.


title: Subradiant-to-Subradiant Phase Transition in the Bad Cavity Laser

abstract: We show that the onset of steady-state superradiance in a bad cavity laser is preceded by a dissipative phase transition between two distinct phases of steady-state subradiance. The transition is marked by a non-analytic behavior of the cavity output power and the mean atomic inversion, as well as a discontinuity in the variance of the collective atomic inversion. In particular, for repump rates below a critical value, the cavity output power is strongly suppressed and does not increase with the atom number, while it scales linearly with atom number above this value. Remarkably, we find that the atoms are in a macroscopic entangled steady state near the critical region with a vanishing fraction of unentangled atoms in the large atom number limit.


title: Generating and Characterizing Scenarios for Safety Testing of Autonomous Vehicles

abstract: Extracting interesting scenarios from real-world data as well as generating failure cases is important for the development and testing of autonomous systems. We propose efficient mechanisms to both characterize and generate testing scenarios using a state-of-the-art driving simulator. For any scenario, our method generates a set of possible driving paths and identifies all the possible safe driving trajectories that can be taken starting at different times, to compute metrics that quantify the complexity of the scenario. We use our method to characterize real driving data from the Next Generation Simulation (NGSIM) project, as well as adversarial scenarios generated in simulation. We rank the scenarios by defining metrics based on the complexity of avoiding accidents and provide insights into how the AV could have minimized the probability of incurring an accident. We demonstrate a strong correlation between the proposed metrics and human intuition.


title: Scaling limits of branching random walks and branching stable processes

abstract: Branching-stable processes have recently appeared as counterparts of stable subordinators, when addition of real variables is replaced by branching mechanism for point processes. Here, we are interested in their domains of attraction and describe explicit conditions for a branching random walk to converge after a proper magnification to a branching-stable process. This contrasts with deep results that have been obtained during the last decade on the asymptotic behavior of branching random walks and which involve either shifting without rescaling, or demagnification.


title: Optimal sequential decision making with probabilistic digital twins

abstract: Digital twins are emerging in many industries, typically consisting of simulation models and data associated with a specific physical system. One of the main reasons for developing a digital twin, is to enable the simulation of possible consequences of a given action, without the need to interfere with the physical system itself. Physical systems of interest, and the environments they operate in, do not always behave deterministically. Moreover, information about the system and its environment is typically incomplete or imperfect. Probabilistic representations of systems and environments may therefore be called for, especially to support decisions in application areas where actions may have severe consequences. In this paper we introduce the probabilistic digital twin (PDT). We will start by discussing how epistemic uncertainty can be treated using measure theory, by modelling epistemic information via $\sigma$-algebras. Based on this, we give a formal definition of how epistemic uncertainty can be updated in a PDT. We then study the problem of optimal sequential decision making. That is, we consider the case where the outcome of each decision may inform the next. Within the PDT framework, we formulate this optimization problem. We discuss how this problem may be solved (at least in theory) via the maximum principle method or the dynamic programming principle. However, due to the curse of dimensionality, these methods are often not tractable in practice. To mend this, we propose a generic approximate solution using deep reinforcement learning together with neural networks defined on sets. We illustrate the method on a practical problem, considering optimal information gathering for the estimation of a failure probability.


title: Photonic Computing to Accelerate Data Processing in Wireless Communications

abstract: Massive multiple-input multiple-output (MIMO) systems are considered as one of the leading technologies employed in the next generations of wireless communication networks (5G), which promise to provide higher spectral efficiency, lower latency, and more reliability. Due to the massive number of devices served by the base stations (BS) equipped with large antenna arrays, massive-MIMO systems need to perform high-dimensional signal processing in a considerably short amount of time. The computational complexity of such data processing, while satisfying the energy and latency requirements, is beyond the capabilities of the conventional widely-used digital electronics-based computing, i.e., Field-Programmable Gate Arrays (FPGAs) and Application-Specific Integrated Circuits (ASICs). In this paper, the speed and lossless propagation of light is exploited to introduce a photonic computing approach that addresses the high computational complexity required by massive-MIMO systems. The proposed computing approach is based on photonic implementation of multiply and accumulate (MAC) operation achieved by broadcast-and-weight (B&W) architecture. The B&W protocol is limited to real and positive values to perform MAC operations. In this work, preprocessing steps are developed to enable the proposed photonic computing architecture to accept any arbitrary values as the input. This is a requirement for wireless communication systems that typically deal with complex values. Numerical analysis shows that the performance of the wireless communication system is not degraded by the proposed photonic computing architecture, while it provides significant improvements in time and energy efficiency for massive-MIMO systems as compared to the most powerful Graphics Processing Units (GPUs).


title: Electricity intraday price modeling with marked Hawkes processes

abstract: We consider a 2-dimensional marked Hawkes process with increasing baseline intensity in order to model prices on electricity intraday markets. This model allows to represent different empirical facts such as increasing market activity, random jump sizes but above all microstructure noise through the signature plot. This last feature is of particular importance for practitioners and has not yet been modeled on those particular markets. We provide analytic formulas for first and second moments and for the signature plot, extending the classic results of Bacry et al. (2013) in the context of Hawkes processes with random jump sizes and time dependent baseline intensity. The tractable model we propose is estimated on German data and seems to fit the data well. We also provide a result about the convergence of the price process to a Brownian motion with increasing volatility at macroscopic scales, highlighting the Samuelson effect.


title: Extension of Kodama vector and quasilocal quantities in three-dimensional axisymmetric spacetimes

abstract: Spherically symmetric spacetimes admit the so-called Kodama vector, which provides a locally conserved current and a preferred time even for dynamical spacetime without any time translation symmetry. A charge associated with this conserved current leads to a quasilocal mass which agrees with the Misner-Sharp mass. In three dimensions, spherically symmetric spacetimes correspond to axisymmetric ones, while axisymmetry allows spacetimes to be rotating with angular momentum. We extend the notion of the Kodama vector to axisymmetric rotating spacetimes in three dimensions. We also define a quasilocal mass taking into account angular momentum in three-dimensional axisymmetric spacetimes.


title: Free-stream coherent structures in parallel compressible boundary-layer flows at subsonic Mach numbers

abstract: As a first step towards the asymptotic description of coherent structures in compressible shear flows, we present a description of nonlinear equilibrium solutions of the Navier--Stokes equations in the compressible asymptotic suction boundary layer (ASBL). The free-stream Mach number is assumed to be $< 0.8$ so that the flow is in the subsonic regime and we assume a perfect gas. We extend the large-Reynolds number free-stream coherent structure theory of \cite{deguchi_hall_2014a} for incompressible ASBL flow to describe a nonlinear interaction in a thin layer situated just below the free-stream which produces streaky disturbances to both the velocity and temperature fields, which can grow exponentially towards the wall. We complete the description of the growth of the velocity and thermal streaks throughout the flow by solving the compressible boundary-region equations numerically. We show that the velocity and thermal streaks obtain their maximum amplitude in the unperturbed boundary layer. Increasing the free-stream Mach number enhances the thermal streaks, whereas varying the Prandtl number changes the location of the maximum amplitude of the thermal streak relative to the velocity streak. Such nonlinear equilibrium states have been implicated in shear transition in incompressible flows; therefore, our results indicate that a similar mechanism may also be present in compressible flows.


title: A resampling approach for causal inference on novel two-point time-series with application to identify risk factors for type-2 diabetes and cardiovascular disease

abstract: Two-point time-series data, characterized by baseline and follow-up observations, are frequently encountered in health research. We study a novel two-point time series structure without a control group, which is driven by an observational routine clinical dataset collected to monitor key risk markers of type-$2$ diabetes (T2D) and cardiovascular disease (CVD). We propose a resampling approach called 'I-Rand' for independently sampling one of the two time points for each individual and making inference on the estimated causal effects based on matching methods. The proposed method is illustrated with data from a service-based dietary intervention to promote a low-carbohydrate diet (LCD), designed to impact risk of T2D and CVD. Baseline data contain a pre-intervention health record of study participants, and health data after LCD intervention are recorded at the follow-up visit, providing a two-point time-series pattern without a parallel control group. Using this approach we find that obesity is a significant risk factor of T2D and CVD, and an LCD approach can significantly mitigate the risks of T2D and CVD. We provide code that implements our method.


title: A Normal Form Algorithm for Tensor Rank Decomposition

abstract: We propose a new numerical algorithm for computing the tensor rank decomposition or canonical polyadic decomposition of higher-order tensors subject to a rank and genericity constraint. Reformulating this computational problem as a system of polynomial equations allows us to leverage recent numerical linear algebra tools from computational algebraic geometry. We characterize the complexity of our algorithm in terms of an algebraic property of this polynomial system -- the multigraded regularity. We prove effective bounds for many tensor formats and ranks, which are of independent interest for overconstrained polynomial system solving. Moreover, we conjecture a general formula for the multigraded regularity, yielding a (parameterized) polynomial time complexity for the tensor rank decomposition problem in the considered setting. Our numerical experiments show that our algorithm can outperform state-of-the-art numerical algorithms by an order of magnitude in terms of accuracy, computation time, and memory consumption.


title: Superconductivity provides a giant enhancement to the spin battery effect

abstract: We develop a theory of the spin battery effect in superconductor/ferromagnetic insulator (SC/FI) systems taking into account the magnetic proximity effect. We demonstrate that the spin-energy mixing enabled by the superconductivity leads to the enhancement of spin accumulation by several orders of magnitude relative to the normal state. This finding can explain the recently observed giant inverse spin Hall effect generated by thermal magnons in the SC/FI system. We suggest a non-local electrical detection scheme which can directly probe the spin accumulation driven by the magnetization dynamics. We predict a giant Seebeck effect converting the magnon temperature bias into the non-local voltage signal. We also show how this can be used to enhance the sensitivity of magnon detection even up to the single-magnon level.


title: Understanding and mitigating noise in trained deep neural networks

abstract: Deep neural networks unlocked a vast range of new applications by solving tasks of which many were previously deemed as reserved to higher human intelligence. One of the developments enabling this success was a boost in computing power provided by special purpose hardware, such as graphic or tensor processing units. However, these do not leverage fundamental features of neural networks like parallelism and analog state variables. Instead, they emulate neural networks relying on binary computing, which results in unsustainable energy consumption and comparatively low speed. Fully parallel and analogue hardware promises to overcome these challenges, yet the impact of analogue neuron noise and its propagation, i.e. accumulation, threatens rendering such approaches inept. Here, we determine for the first time the propagation of noise in deep neural networks comprising noisy nonlinear neurons in trained fully connected layers. We study additive and multiplicative as well as correlated and uncorrelated noise, and develop analytical methods that predict the noise level in any layer of symmetric deep neural networks or deep neural networks trained with back propagation. We find that noise accumulation is generally bound, and adding additional network layers does not worsen the signal to noise ratio beyond a limit. Most importantly, noise accumulation can be suppressed entirely when neuron activation functions have a slope smaller than unity. We therefore developed the framework for noise in fully connected deep neural networks implemented in analog systems, and identify criteria allowing engineers to design noise-resilient novel neural network hardware.


title: Real-time Nonrigid Mosaicking of Laparoscopy Images

abstract: The ability to extend the field of view of laparoscopy images can help the surgeons to obtain a better understanding of the anatomical context. However, due to tissue deformation, complex camera motion and significant three-dimensional (3D) anatomical surface, image pixels may have non-rigid deformation and traditional mosaicking methods cannot work robustly for laparoscopy images in real-time. To solve this problem, a novel two-dimensional (2D) non-rigid simultaneous localization and mapping (SLAM) system is proposed in this paper, which is able to compensate for the deformation of pixels and perform image mosaicking in real-time. The key algorithm of this 2D non-rigid SLAM system is the expectation maximization and dual quaternion (EMDQ) algorithm, which can generate smooth and dense deformation field from sparse and noisy image feature matches in real-time. An uncertainty-based loop closing method has been proposed to reduce the accumulative errors. To achieve real-time performance, both CPU and GPU parallel computation technologies are used for dense mosaicking of all pixels. Experimental results on \textit{in vivo} and synthetic data demonstrate the feasibility and accuracy of our non-rigid mosaicking method.


title: From Nano-Communications to Body Area Networks: A Perspective on Truly Personal Communications

abstract: This article presents an overview of future truly personal communications, ranging from networking inside the human body to the exchange of data with external wireless devices in the surrounding environment. At the nano- and micro-scales, communications can be realized with the aid of molecular mechanisms, Forster resonance energy transfer phenomenon, electromagnetic or ultrasound waves. At a larger scale, in the domain of Body Area Networks, a wide range of communication mechanisms is available, including smart-textiles, inductive- and body-couplings, ultrasounds, optical and wireless radio transmissions, a number of mature technologies existing already. The main goal of this article is to identify the potential mechanisms that can be exploited to provide interfaces in between nano- and micro-scale systems and Body Area Networks. These interfaces have to bridge the existing gap between the two worlds, in order to allow for truly personal communication systems to become a reality. The extraordinary applications of such systems are also discussed, as they are strong drivers of the research in this area.


title: Performance study and calibration strategy of the HADES scintillator TOF Wall with fast digital readout

abstract: We present in this work the calibration procedure and a performance study of long scintillator bars used for the time-of-flight (TOF) measurement in the HADES experiment. The digital front-end electronics installed at the TOF detector required to develop novel calibration methods. The exceptional performance of the spectrometer for particle identification and pointing accuracy allows one to determine in great detail the response of scintillators to minimum ionizing particles. A substantial position sensitivity of the calibration parameters has been found, in particular for the signal time walk. After including the position dependence, the timing accuracy for minimum ionizing particles was improved from 190~ps to 135~ps for the shortest rods (1475 mm) and to 165~ps for the longest (2356 mm). These results are in accordance with the time degradation length of the scintillator bars, as determined from previous measurements.


title: The Role of Adaptive Ray Tracing in Analyzing Black Hole Structure

abstract: The recent advent of the Event Horizon Telescope (EHT) has made direct imaging of supermassive black holes a reality. Simulated images of black holes produced via general relativistic ray tracing and radiative transfer provide a key counterpart to these observational efforts. Black hole images have a wide range of physically interesting image structures, ranging from extremely fine scales in their lensed "photon rings" to the very large scales in their relativistic jets. The multi-scale nature of the black hole system is therefore suitable for a multi-scale approach to generating simulated images that capture all key elements of the system. Here, we present a prescription for adaptive ray tracing, which enables efficient computation of extremely high resolution images of black holes. Using the polarized ray-tracing code ipole, we image a combination of semi-analytic and GRMHD models, and we show that images can be reproduced with mean squared error of less than 0.1% even after tracing 12x fewer rays. We then use adaptive ray tracing to explore properties of the photon ring. We illustrate the behavior of individual subrings in GRMHD simulations, and we explore their signatures in interferometric visibilities.


title: Vapor flux on bumpy surfaces: condensation and transpiration on leaves

abstract: Drop condensation and evaportation as a result of the gradient in vapor concentration are important in both engineering and natural systems. One of the interesting natural examples is transpiration on plant leaves. Most of water in the inner space of the leaves escapes through stomata, whose rate depends on the surface topography and a difference in vapor concentrations inside and just outside of the leaves. Previous research on the vapor flux on various surfaces has focused on numerically solving the vapor diffusion equation or using scaling arguments based on a simple solution with a flat surface. In this present work, we present and discuss simple analytical solutions on various 2D surface shapes (e.g., semicylinder, semi-ellipse, hair). The method of solving the diffusion equation is to use the complex potential theory, which provides analytical solutions for vapor concentration and flux. We find that a high mass flux of vapor is formed near the top of the microstructures while a low mass flux is developed near the stomata at the leaf surface. Such a low vapor flux near the stomata may affect transpiration in two ways. First, condensed droplets on the stomata will not grow due to a low mass flux of vapor, which will not inhibit the gas exchange through the stomatal opening. Second, the low mass flux from the atmosphere will facilitate the release of high concentrated vapor from the substomatal space.


title: Non-universal power-law dynamics of SIR models on hierarchical modular networks

abstract: Power-law (PL) time dependent infection growth has been reported in many COVID-19 statistics. In simple SIR models the number of infections grows at the outbreak as $I(t) \propto t^{d-1}$ on $d$-dimensional Euclidean lattices in the endemic phase or follow a slower universal PL at the critical point, until finite sizes cause immunity and a crossover to an exponential decay. Heterogeneity may alter the dynamics of spreading models, spatially inhomogeneous infection rates can cause slower decays, posing a threat of a long recovery from a pandemic. COVID-19 statistics have also provided epidemic size distributions with PL tails in several countries. Here I investigate SIR like models on hierarchical modular networks, embedded in 2d lattices with the addition of long-range links. I show that if the topological dimension of the network is finite, average degree dependent PL growth of prevalence emerges. Supercritically the same exponents as of regular graphs occurs, but the topological disorder alters the critical behavior. This is also true for the epidemic size distributions. Mobility of individuals does not affect the form of the scaling behavior, except for the $d=2$ lattice, but increases the magnitude of the epidemic. The addition of a super-spreader hot-spot also does not change the growth exponent and the exponential decay in the herd immunity regime.


title: Mapping quantum chemical dynamics problems onto spin-lattice simulators

abstract: The accurate computational determination of chemical, materials, biological, and atmospheric properties has critical impact on a wide range of health and environmental problems, but is deeply limited by the computational scaling of quantum-mechanical methods. The complexity of quantum-chemical studies arises from the steep algebraic scaling of electron correlation methods, and the exponential scaling in studying nuclear dynamics and molecular flexibility. To date, efforts to apply quantum hardware to such quantum chemistry problems have focused primarily on electron correlation. Here, we provide a framework which allows for the solution of quantum chemical nuclear dynamics by mapping these to quantum spin-lattice simulators. Using the example case of a short-strong hydrogen bonded system, we construct the Hamiltonian for the nuclear degrees of freedom on a single Born-Oppenheimer surface and show how it can be transformed to a generalized Ising model Hamiltonian. We then demonstrate a method to determine the local fields and spin-spin couplings needed to identically match the molecular and spin-lattice Hamiltonians. We describe a protocol to determine the on-site and inter-site coupling parameters of this Ising Hamiltonian from the Born-Oppenheimer potential and nuclear kinetic energy operator. Our approach represents a paradigm shift in the methods used to study quantum nuclear dynamics, opening the possibility to solve both electronic structure and nuclear dynamics problems using quantum computing systems.


title: A Minkowski inequality for Horowitz-Myers geon

abstract: We prove a sharp inequality for toroidal hypersurfaces in three and four dimensional Horowitz-Myers geon. This extend previous results on Minkowski inequality in the static spacetime to toroidal surfaces in asymptotically hyperbolic manifold with flat toroidal conformal infinity.


title: Distinguished categories and the Zilber-Pink conjecture

abstract: We propose an axiomatic approach towards studying unlikely intersections by introducing the framework of distinguished categories. This includes commutative algebraic groups and mixed Shimura varieties. It allows us to define all basic concepts of the field and prove some fundamental facts about them, e.g. the defect condition. In some categories that we call very distinguished, we are able to show some implications between Zilber-Pink statements with respect to base change. This yields unconditional results, i.e. the Zilber-Pink conjecture for a complex curve in $\mathcal{A}_2$ that cannot be defined over $\bar{\mathbb{Q}}$, a complex curve in the $g$-th fibered power of the Legendre family, and a complex curve in the base change of a semiabelian variety over $\bar{\mathbb{Q}}$.


title: Radiomic Deformation and Textural Heterogeneity (R-DepTH) Descriptor to characterize Tumor Field Effect: Application to Survival Prediction in Glioblastoma

abstract: The concept of tumor field effect implies that cancer is a systemic disease with its impact way beyond the visible tumor confines. For instance, in Glioblastoma (GBM), an aggressive brain tumor, the increase in intracranial pressure due to tumor burden often leads to brain herniation and poor outcomes. Our work is based on the rationale that highly aggressive tumors tend to grow uncontrollably, leading to pronounced biomechanical tissue deformations in the normal parenchyma, which when combined with local morphological differences in the tumor confines on MRI scans, will comprehensively capture tumor field effect. Specifically, we present an integrated MRI-based descriptor, radiomic-Deformation and Textural Heterogeneity (r-DepTH). This descriptor comprises measurements of the subtle perturbations in tissue deformations throughout the surrounding normal parenchyma due to mass effect. This involves non-rigidly aligning the patients MRI scans to a healthy atlas via diffeomorphic registration. The resulting inverse mapping is used to obtain the deformation field magnitudes in the normal parenchyma. These measurements are then combined with a 3D texture descriptor, Co-occurrence of Local Anisotropic Gradient Orientations (COLLAGE), which captures the morphological heterogeneity within the tumor confines, on MRI scans. R-DepTH, on N = 207 GBM cases (training set (St) = 128, testing set (Sv) = 79), demonstrated improved prognosis of overall survival by categorizing patients into low- (prolonged survival) and high-risk (poor survival) groups (on St, p-value = 0.0000035, and on Sv, p-value = 0.0024). R-DepTH descriptor may serve as a comprehensive MRI-based prognostic marker of disease aggressiveness and survival in solid tumors.


title: Measuring Reciprocity in a Directed Preferential Attachment Network

abstract: Empirical studies show that online social networks have not only in- and out-degree distributions with Pareto-like tails but also a high proportion of reciprocal edges. A classical directed preferential attachment (PA) model generates in- and out-degree distribution with power-law tails, but theoretical properties of the reciprocity feature in this model have not yet been studied. We derive the asymptotic results on the number of reciprocal edges between two fixed nodes, as well as the proportion of reciprocal edges in the entire PA network. We see that with certain choices of parameters, the proportion of reciprocal edges in a directed PA network is close to 0, which differs from the empirical observation. This points out one potential problem of fitting a classical PA model to a given network dataset with high reciprocity and indicates alternative models need to be considered.


title: Fast, Scalable Approximations to Posterior Distributions in Extended Latent Gaussian Models

abstract: We define a novel class of additive models, called Extended Latent Gaussian Models, that allow for a wide range of response distributions and flexible relationships between the additive predictor and mean response. The new class covers a broad range of interesting models including multi-resolution spatial processes, partial likelihood-based survival models, and multivariate measurement error models. Because computation of the exact posterior distribution is infeasible, we develop a fast, scalable approximate Bayesian inference methodology for this class based on nested Gaussian, Laplace, and adaptive quadrature approximations. We prove that the error in these approximate posteriors is op(1) under standard conditions, and provide numerical evidence suggesting that our method runs faster and scales to larger datasets than methods based on Integrated Nested Laplace Approximations and Markov Chain Monte Carlo, with comparable accuracy. We apply the new method to the mapping of malaria incidence rates in continuous space using aggregated data, mapping leukaemia survival hazards using a Cox Proportional-Hazards model with a continuously-varying spatial process, and estimating the mass of the Milky Way Galaxy using noisy multivariate measurements of the positions and velocities of star clusters in its orbit.


title: Magnetoelectric torque and edge currents in spin-orbit coupled graphene nanoribbons

abstract: For graphene nanoribbons with Rashba spin-orbit coupling, the peculiar magnetic response due to the presence of a magnetization and geometric confinement are analyzed within a tight-binding model. We observe a sizable transverse susceptibility that can be considered as a gate voltage-induced magnetoelectric torque without the need of a bias voltage, with different directions for zigzag and armchair ribbons. The local torque generates non-collinear spin polarization between the two edges and/or along the ribbon, and the net torque averages to zero if the magnetization is homogeneous. Nevertheless, a nonzero net torque can appear in partially magnetized nanoribbons or in nanoflakes of irregular shapes. The equilibrium spin current produced by the spin-orbit coupling also appears in nanoribbons, but the component flowing in the direction of confinement is strongly suppressed. Even without the magnetization, an out-of-plane polarized chiral edge spin current is produced, resembling that in the quantum spin Hall effect. Moreover, a magnetization pointing perpendicular to the edge produces a laminar flow of edge charge currents, whose flow direction is symmetric (non chiral) or antisymmetric (chiral) between the two edges depends on whether the magnetization points in-plane or out-of-plane.


title: Cold quark matter at N3LO: soft contributions

abstract: High-order perturbative calculations for thermodynamic quantities in QCD are complicated by the physics of dynamical screening that affects the soft, long-wavelength modes of the system. Here, we provide details for the evaluation of this soft contribution to the next-to-next-to-next-to-leading order (N3LO) pressure of high-density, zero-temperature quark matter (QM), complementing our accompanying paper in arXiv:2103.05658. Our calculation requires the determination of the pressure of the hard-thermal-loop (HTL) effective theory to full two-loop order at zero temperature, which we go through in considerable detail. In addition to this, we comprehensively discuss the structure of the weak-coupling expansion of the QM pressure, and lay out a roadmap towards the evaluation of the contributions missing from a full N3LO result for this quantity.


title: Genetic Improvement of Routing Protocols for Delay Tolerant Networks

abstract: Routing plays a fundamental role in network applications, but it is especially challenging in Delay Tolerant Networks (DTNs). These are a kind of mobile ad hoc networks made of e.g. (possibly, unmanned) vehicles and humans where, despite a lack of continuous connectivity, data must be transmitted while the network conditions change due to the nodes' mobility. In these contexts, routing is NP-hard and is usually solved by heuristic "store and forward" replication-based approaches, where multiple copies of the same message are moved and stored across nodes in the hope that at least one will reach its destination. Still, the existing routing protocols produce relatively low delivery probabilities. Here, we genetically improve two routing protocols widely adopted in DTNs, namely Epidemic and PRoPHET, in the attempt to optimize their delivery probability. First, we dissect them into their fundamental components, i.e., functionalities such as checking if a node can transfer data, or sending messages to all connections. Then, we apply Genetic Improvement (GI) to manipulate these components as terminal nodes of evolving trees. We apply this methodology, in silico, to six test cases of urban networks made of hundreds of nodes, and find that GI produces consistent gains in delivery probability in four cases. We then verify if this improvement entails a worsening of other relevant network metrics, such as latency and buffer time. Finally, we compare the logics of the best evolved protocols with those of the baseline protocols, and we discuss the generalizability of the results across test cases.


title: Constant-Depth Circuits for Dynamic Simulations of Materials on Quantum Computers

abstract: Dynamic simulation of materials is a promising application for near-term quantum computers. Current algorithms for Hamiltonian simulation, however, produce circuits that grow in depth with increasing simulation time, limiting feasible simulations to short-time dynamics. Here, we present a method for generating circuits that are constant in depth with increasing simulation time for a subset of one-dimensional materials Hamiltonians, thereby enabling simulations out to arbitrarily long times. Furthermore, by removing the effective limit on the number of feasibly simulatable time-steps, the constant-depth circuits enable Trotter error to be made negligibly small by allowing simulations to be broken into arbitrarily many time-steps. Composed of two-qubit matchgates on nearest-neighbor qubits, these constant-depth circuits are constructed based on a set of multi-matchgate identity relationships. For an $N$-spin system, the constant-depth circuit contains only $\mathcal{O}(N^2)$ CNOT gates. When compared to standard Hamiltonian simulation algorithms, our method generates circuits with order-of-magnitude fewer gates, which allows us to successfully simulate the long-time dynamics of systems with up to 5 spins on available quantum hardware. This paves the way for simulations of long-time dynamics for scientifically and technologically relevant quantum materials, enabling the observation of interesting and important atomic-level physics.


title: Flexoelectricity and the Entropic Force between Fluctuating Fluid Membranes

abstract: Biological membranes undergo noticeable thermal fluctuations at physiological temperatures. When two membranes approach each other, they hinder the out of plane fluctuations of the other. This hindrance leads to an entropic repulsive force between membranes which, in an interplay with attractive and repulsive forces due to other sources, impacts a range of biological functions: cell adhesion, membrane fusion, self-assembly, binding-unbinding transition among others. In this work, we take cognizance of the fact that biological membranes are not purely mechanical entities and, due to the phenomenon of flexoelectricity, exhibit a coupling between deformation and electric polarization. The ensuing coupled mechanics-electrostatics-statistical mechanics problem is analytically intractable. We use a variational perturbation method to analyze, in closed-form, the contribution of flexoelectricity to the entropic force between two fluctuating membranes and discuss its possible physical implications. We find that flexoelectricity leads to a correction that switches from an enhanced attraction at close membrane separations and an enhanced repulsion when the membranes are further apart.


title: Hypothesis-based acceptance sampling for modules F and F1 of the European Measuring Instruments Directive

abstract: Millions of measuring instruments are verified each year before being placed on the markets worldwide. In the EU, such initial conformity assessments are regulated by the Measuring Instruments Directive (MID). The MID modules F and F1 on product verification allow for statistical acceptance sampling, whereby only random subsets of instruments need to be inspected. This paper re-interprets the acceptance sampling conditions formulated by the MID. The new interpretation is contrasted with the one advanced in WELMEC guide 8.10, and three advantages have become apparent. Firstly, an economic advantage of the new interpretation is a producers' risk bounded from above, such that measuring instruments with sufficient quality are accepted with a guaranteed probability of no less than 95 %. Secondly, a conceptual advantage is that the new MID interpretation fits into the well-known, formal framework of statistical hypothesis testing. Thirdly, the new interpretation applies unambiguously to finite-sized lots, even very small ones. We conclude that the new interpretation is to be preferred and suggest re-formulating the statistical sampling conditions in the MID. Re-interpreting the MID conditions implies that currently available sampling plans are either not admissible or not optimal. We derive a new acceptance sampling scheme and recommend its application.


title: Characterizing and correcting the proper motion bias of the bright Gaia EDR3 sources

abstract: In this paper we characterize magnitude-dependent systematics in the proper motions of the Gaia EDR3 catalog and provide a prescription for their removal. The reference frame of bright stars (G<13) in EDR3 is known to rotate with respect to extragalactic objects, but this rotation has proven difficult to characterize and correct. We employ a sample of binary stars and a sample of open cluster members to characterize this proper motion bias as a magnitude-dependent spin of the reference frame. We show that the bias varies with G magnitude, reaching up to 80 {\mu}as/yr for sources in the range G = 11 - 13, several times the formal EDR3 proper motion uncertainties. We also show evidence for an additional dependence on the color of the source, with a magnitude up to 10 {\mu}as/yr. However, a color correction proportional to the effective wavenumber is unsatisfactory for very red or very blue stars and we do not recommend its use. We provide a recipe for a magnitude-dependent correction to align the proper motion of the Gaia EDR3 sources brighter than G=13 with the International Celestial Reference Frame.


title: Quantum Computing: Towards Industry Reference Problems

abstract: The complexity is increasing rapidly in many areas of the automotive industry. The design of an automobile involves many different engineering disciplines, e. g., mechanical, electrical, and software engineering. The software of a vehicle comprises millions of lines of code. Further, the manufacturing, logistics, distribution, and sales of a vehicle are highly complex. There is an immense need for solving simulation problems, e. g., in battery chemistry, an essential enabler for technological advancements for electric vehicles. In all these domains, myriads of optimization, simulation, and machine learning problems arise. Quantum computing-based approaches promise to overcome some of the inherent scalability limitations of classical approaches. This article investigates quantum computing applications across the automotive value chain and identifies several high-value problems that will benefit from quantum-enhanced solutions.


title: Evolution of behaviors in heterogeneous traffic models as driven annealed disorders and its relation to the n-vector model

abstract: In one-dimensional, heterogeneous systems, the whole traffic dynamics depend strongly on the behavior of the leading vehicle. This result holds for a class of vehicular traffic models satisfying the following properties. The interactions are unidirectional. The dynamics of the particles maximize the velocity or reduces the gap between particles. The particles are hard. We use this result to show a link between traffic models and graphs theory with the consequence that as driving styles spread through social contagion and appear randomly, the inhomogeneities of the system becomes dynamical, or \textit{annealed}, toward specific regions in the space of parameters. Interpreting parameters as entries of vectors defined in the parameters space appear analogies between the evolutionary dynamics of these systems and asymptotic behaviors of the \textit {n-vector model}. When the time-scale ratio of imitation to the mutation processes, $ \tau_i / \tau_m $, is small an organized state where "orientation" corresponds to the set of parameters of the slowest strategies is favored, and if this ratio is big an unorganized state without a preferential orientation is favored.


title: Markov intertwining operators, joinings, and asymptotic properties of dynamical systems

abstract: This text is written based on the author's publications during the period from 1991 to 2001. The work is devoted to the theory of Markov intertwining operators and joinings of measure-preserving group actions, as well as to their applications to study asymptotic properties of dynamical systems. Special attention is paid to Rokhlin's problems on multiple mixing and multiple spectrum. The development of these topics over the past twenty years has not been discussed. In fact many results on joinings have frozen in time, many questions have remained open without losing their relevance, but probably have ceased to excite interest due to difficulties. For example, it is not known whether the minimal self-joinings of order 2 imply all orders? Is there a non-trivial pairwise independent joining for a weakly mixing system of zero entropy? What can be said about such joinings for transformations with small local rank? These questions are ripe for a long time, and the author reminds the reader about them, combining his story with numerous partial and related results.


title: The Devinatz-Hopkins Theorem via Algebraic Geometry

abstract: In this note, we show how a continuous action of the Morava stabilizer group $\mathbb G_n$ on the Lubin-Tate spectrum $E_n$, satisfying the conclusion $E_n^{h\mathbb G_n}\simeq L_{K(n)} S$ of the Devinatz-Hopkins Theorem, may be obtained by monodromy on the stack of oriented deformations of formal groups in the context of formal spectral algebraic geometry.


title: Hyperspectral Image Denoising and Anomaly Detection Based on Low-rank and Sparse Representations

abstract: Hyperspectral imaging measures the amount of electromagnetic energy across the instantaneous field of view at a very high resolution in hundreds or thousands of spectral channels. This enables objects to be detected and the identification of materials that have subtle differences between them. However, the increase in spectral resolution often means that there is a decrease in the number of photons received in each channel, which means that the noise linked to the image formation process is greater. This degradation limits the quality of the extracted information and its potential applications. Thus, denoising is a fundamental problem in hyperspectral image (HSI) processing. As images of natural scenes with highly correlated spectral channels, HSIs are characterized by a high level of self-similarity and can be well approximated by low-rank representations. These characteristics underlie the state-of-the-art methods used in HSI denoising. However, where there are rarely occurring pixel types, the denoising performance of these methods is not optimal, and the subsequent detection of these pixels may be compromised. To address these hurdles, in this article, we introduce RhyDe (Robust hyperspectral Denoising), a powerful HSI denoiser, which implements explicit low-rank representation, promotes self-similarity, and, by using a form of collaborative sparsity, preserves rare pixels. The denoising and detection effectiveness of the proposed robust HSI denoiser is illustrated using semireal and real data.


title: Four-Spin Terms and the Origin of the Chiral Spin Liquid in Mott Insulators on the Triangular Lattice

abstract: At strong repulsion, the triangular-lattice Hubbard model is described by $s=1/2$ spins with nearest-neighbor antiferromagnetic Heisenberg interactions and exhibits conventional 120$^\circ$ order. Using the infinite density matrix renormalization group and exact diagonalization, we study the effect of the additional four-spin interactions naturally generated from the underlying Mott-insulator physics of electrons as the repulsion decreases. Although these interactions have historically been connected with a gapless ground state with emergent spinon Fermi surface, we find that at physically relevant parameters, they stabilize a chiral spin-liquid (CSL) of Kalmeyer-Laughlin (KL) type, clarifying observations in recent studies of the Hubbard model. We then present a self-consistent solution based on mean-field rewriting of the interaction to obtain a Hamiltonian with similarities to the parent Hamiltonian of the KL state, providing a physical understanding for the origin of the CSL.


title: A Lyapunov-based ISS small-gain theorem for infinite networks of nonlinear systems

abstract: In this paper, we show that an infinite network of input-to-state stable (ISS) subsystems, admitting ISS Lyapunov functions, itself admits an ISS Lyapunov function, provided that the couplings between the subsystems are sufficiently weak. The strength of the couplings is described in terms of the properties of an infinite-dimensional nonlinear positive operator, built from the interconnection gains. If this operator induces a uniformly globally asymptotically stable (UGAS) system, a Lyapunov function for the infinite network can be constructed. We analyze necessary and sufficient conditions for UGAS and relate them to small-gain conditions used in the stability analysis of finite networks.


title: A closed-form approximation for pricing geometric Istanbul options

abstract: The Istanbul options were first introduced by Michel Jacques in 1997. These derivatives are considered as an extension of the Asian options. In this paper, we propose an analytical approximation formula for a geometric Istanbul call option (GIC) under the Black-Scholes model. Our approximate pricing formula is obtained in closed-form using a second-order Taylor expansion. We compare our theoretical results with those of Monte-Carlo simulations using the control variates method. Finally, we study the effects of changes in the price of the underlying asset on the value of GIC.


title: Hard Lefschetz Property for Isometric Flows

abstract: The Hard Lefschetz Property (HLP) is an important property which has been studied in several categories of the symplectic world. For Sasakian manifolds, this duality is satisfied by the basic cohomology (so, it is a transverse property), but a new version of the HLP has been recently given in terms of duality of the cohomology of the manifold itself in arXiv:1306.2896. Both properties were proved to be equivalent (see arXiv:1311.1431) in the case of K-contact flows. In this paper we extend both versions of the HLP (transverse and not) to the more general category of isometric flows, and show that they are equivalent. We also give some explicit examples which illustrate the categories where the HLP could be considered.


title: Quark Density in Lattice QC$_2$D at Imaginary and Real Chemical Potential

abstract: We study lattice two-color QCD (QC$_2$D) with two flavors of staggered fermions at imaginary and real quark chemical potential $\mu_q$ and $T>T_c$. We employ various methods of analytic continuation of the quark number density from imaginary to real quark chemical potentials $\mu_q$, including series expansions as well as based on phenomenological models, and study their accuracy by comparing the results to the lattice data. Below the Roberge-Weiss temperature, $T<T_{RW}$, we find that the cluster expansion model provides an accurate analytic continuation of the baryon number density in the studied range of chemical potentials. On the other hand, the behavior of the reconstructed canonical partition functions indicate that the available models may require corrections at high quark densities. At $T>T_{RW}$ we show that the analytic continuation to the real values of $\mu_q$ based on trigonometric functions works equally well with the conventional method based on the Taylor expansion in powers of $\mu_q$.


title: Symmetry-resolved entanglement detection using partial transpose moments

abstract: We propose an ordered set of experimentally accessible conditions for detecting entanglement in mixed states. The $k$-th condition involves comparing moments of the partially transposed density operator up to order $k$. Remarkably, the union of all moment inequalities reproduces the Peres-Horodecki criterion for detecting entanglement. Our empirical studies highlight that the first four conditions already detect mixed state entanglement reliably in a variety of quantum architectures. Exploiting symmetries can help to further improve their detection capabilities. We also show how to estimate moment inequalities based on local random measurements of single state copies (classical shadows) and derive statistically sound confidence intervals as a function of the number of performed measurements. Our analysis includes the experimentally relevant situation of drifting sources, i.e. non-identical, but independent, state copies.


title: Phases and homogeneous ordered states in alignment-based self-propelled particle models

abstract: We study a set of models of self-propelled particles that achieve collective motion through similar alignment-based dynamics, considering versions with and without repulsive interactions that do not affect the heading directions. We explore their phase space within a broad range of values of two nondimensional parameters (coupling strength and Peclet number), characterizing their polarization and degree of clustering. The resulting phase diagrams display equivalent, similarly distributed regions for all models with repulsion. The diagrams without repulsion exhibit differences, in particular for high coupling strengths. We compare the boundaries and representative states of all regions, identifying various regimes that had not been previously characterized. We analyze in detail three types of homogeneous polarized states, comparing them to existing theoretical and numerical results by computing their velocity and density correlations, giant number fluctuations, and local order-density coupling. We find that they all deviate in one way or another from the theoretical predictions, attributing these differences either to the remaining inhomogeneities or to finite-size effects. We discuss our results in terms of the universal or specific features of each model, their thermodynamic limit, and the high mixing and low mixing regimes. Our study provides a broad, overarching perspective on the multiple phases and states found in alignment-based self-propelled particle models.


title: Efficient reconstruction of depth three circuits with top fan-in two

abstract: We develop efficient randomized algorithms to solve the black-box reconstruction problem for polynomials over finite fields, computable by depth three arithmetic circuits with alternating addition/multiplication gates, such that output gate is an addition gate with in-degree two. These circuits compute polynomials of form $G\times(T_1 + T_2)$, where $G,T_1,T_2$ are product of affine forms, and polynomials $T_1,T_2$ have no common factors. Rank of such a circuit is defined as dimension of vector space spanned by all affine factors of $T_1$ and $T_2$. For any polynomial $f$ computable by such a circuit, $rank(f)$ is defined to be the minimum rank of any such circuit computing it. Our work develops randomized reconstruction algorithms which take as input black-box access to a polynomial $f$ (over finite field $\mathbb{F}$), computable by such a circuit. Here are the results. 1 [Low rank]: When $5\leq rank(f) = O(\log^3 d)$, it runs in time $(nd^{\log^3d}\log |\mathbb{F}|)^{O(1)}$, and, with high probability, outputs a depth three circuit computing $f$, with top addition gate having in-degree $\leq d^{rank(f)}$. 2 [High rank]: When $rank(f) = \Omega(\log^3 d)$, it runs in time $(nd\log |\mathbb{F}|)^{O(1)}$, and, with high probability, outputs a depth three circuit computing $f$, with top addition gate having in-degree two. Ours is the first blackbox reconstruction algorithm for this circuit class, that runs in time polynomial in $\log |\mathbb{F}|$. This problem has been mentioned as an open problem in [GKL12] (STOC 2012)