title: Stochastic approach to epidemic spreading

abstract: We analyze four models of epidemic spreading using a stochastic approach in which the primary stochastic variables are the numbers of individuals in each class. The stochastic approach is described by a master equation and the transition rate for each process such as infection or recovery are set up by using the law of mass action. We perform numerical simulations as well as numerical integration of the evolution equations for the average number of each class of individuals. The onset of the epidemic spreading is obtained by a linear analysis of the disease free state, from which follows the initial exponential increase of the infected and the frequency of new cases. The order parameter and the variance in the number of individuals are also obtained characterizing the onset of epidemic spreading as a critical phase transition.


title: Tridiagonal kernels and left-invertible operators with applications to Aluthge transforms

abstract: Given scalars $a_n (\neq 0)$ and $b_n$, $n \geq 0$, the tridiagonal kernel or band kernel with bandwidth $1$ is the positive definite kernel $k$ on the open unit disc $\mathbb{D}$ defined by \[ k(z, w) = \sum_{n=0}^\infty \Big((a_n + b_n z)z^n\Big) \Big((\bar{a}_n + \bar{b}_n \bar{w}) \bar{w}^n \Big) \qquad (z, w \in \mathbb{D}). \] This defines a reproducing kernel Hilbert space $\mathcal{H}_k$ (known as tridiagonal space) of analytic functions on $\mathbb{D}$ with $\{(a_n + b_nz) z^n\}_{n=0}^\infty$ as an orthonormal basis. We consider shift operators $M_z$ on $\mathcal{H}_k$ and prove that $M_z$ is left-invertible if and only if $\{|{a_n}/{a_{n+1}}|\}_{n\geq 0}$ is bounded away from zero. We find that, unlike the case of weighted shifts, Shimorin's models for left-invertible operators fail to bring to the foreground the tridiagonal structure of shifts. In fact, the tridiagonal structure of a kernel $k$, as above, is preserved under Shimorin model if and only if $b_0=0$ or that $M_z$ is a weighted shift. We prove concrete classification results concerning invariance of tridiagonality of kernels, Shimorin models, and positive operators. We also develop a computational approach to Aluthge transforms of shifts. Curiously, in contrast to direct kernel space techniques, often Shimorin models fails to yield tridiagonal Aluthge transforms of shifts defined on tridiagonal spaces.


title: Deep Local and Global Spatiotemporal Feature Aggregation for Blind Video Quality Assessment

abstract: In recent years, deep learning has achieved promising success for multimedia quality assessment, especially for image quality assessment (IQA). However, since there exist more complex temporal characteristics in videos, very little work has been done on video quality assessment (VQA) by exploiting powerful deep convolutional neural networks (DCNNs). In this paper, we propose an efficient VQA method named Deep SpatioTemporal video Quality assessor (DeepSTQ) to predict the perceptual quality of various distorted videos in a no-reference manner. In the proposed DeepSTQ, we first extract local and global spatiotemporal features by pre-trained deep learning models without fine-tuning or training from scratch. The composited features consider distorted video frames as well as frame difference maps from both global and local views. Then, the feature aggregation is conducted by the regression model to predict the perceptual video quality. Finally, experimental results demonstrate that our proposed DeepSTQ outperforms state-of-the-art quality assessment algorithms.


title: Analysis of landscape hierarchy during coarsening and aging in Ising spin glasses

abstract: We use record dynamics (RD), a coarse-grained description of the ubiquitous relaxation phenomenology known as "aging", as a diagnostic tool to find universal features that distinguish between the energy landscapes of Ising spin models and the ferromagnet. According to RD, a non-equilibrium system after a quench relies on fluctuations that randomly generate a sequence of irreversible record-sized events (quakes or avalanches) that allow the system to escape ever-higher barriers of meta-stable states within a complex, hierarchical energy landscape. Once these record events allow the system to overcome such barriers, the system relaxes by tumbling into the following meta-stable state that is marginally more stable. Within this framework, a clear distinction can be drawn between the coarsening dynamics of an Ising ferromagnet and the aging of the spin glass, which are often put in the same category. To that end, we interpolate between the spin glass and ferromagnet by varying the admixture $p$ of ferromagnetic over anti-ferromagnetic bonds from the glassy state (at 50% each) to wherever clear ferromagnetic behavior emerges. The accumulation of record events grows logarithmic with time in the glassy regime, with a sharp transition at a specific admixture into the ferromagnetic regime where such activations saturate quickly. We show this effect both for the Edwards-Anderson model on a cubic lattice as well as the Sherrington-Kirkpatrick (mean-field) spin glass. While this transition coincides with a previously observed zero-temperature equilibrium transition in the former, that transition has not yet been described for the latter.


title: Precision measurements of A=3 nuclei in Hall B

abstract: We propose a high-statistics measurement of few body nuclear structure and short range correlations in quasi-elastic scattering at 6.6 GeV from $^2$H, $^3$He and $^3$H targets in Hall B with the CLAS12 detector. We will measure absolute cross sections for $(e,e'p)$ and $(e,e'pN)$ quasi-elastic reaction channels up to a missing momentum $p_{miss} \approx 1$ GeV/c over a wide range of $Q^2$ and $x_B$ and construct the isoscalar sum of $^3$H and $^3$He. We will compare $(e,e'p)$ cross sections to nuclear theory predictions using a wide variety of techniques and $NN$ interactions in order to constrain the $NN$ interaction at short distances. We will measure $(e,e'pN)$ quasi-elastic reaction cross sections and $(e,e'pN)/(e,e'p)$ ratios to understand short range correlated (SRC) $NN$ pairs in the simplest non-trivial system. $^3$H and $^3$He, being mirror nuclei, exploit the maximum available isospin asymmetry. They are light enough that their ground states are readily calculable, but they already exhibit complex nuclear behavior, including $NN$ SRCs. We will also measure $^2$H$(e,e'p)$ in order to help theorists constrain non-quasielastic reaction mechanisms in order to better calculate reactions on $A=3$ nuclei. Measuring all three few body nuclei together is critical, in order to understand and minimize different reaction effects, such as single charge exchange final state interactions, in order to test ground-state nuclear models. We will also measure the ratio of inclusive $(e,e')$ quasi-elastic cross sections (integrated over $x_B$) from $^3$He and $^3$H in order to extract the neutron magnetic form factor $G_M^n$ at small and moderate values of $Q^2$. We will measure this at both 6.6 GeV and 2.2 GeV.


title: Attack-resilient observer pruning for path-tracking control of Wheeled Mobile Robot

abstract: Path-tracking control of wheeled mobile robot (WMR) has gained a lot of research attention, primarily because of its wide applicability -- for example intelligent wheelchairs, exploration-assistant remote WMR. Recent increase in remote and autonomous operations\requirements for WMR has led to more and more use of IoT devices within the control loop. Consequently, providing interfaces for malicious interactions through false data injection attacks (FDIA). Moreover, optimization-based FDIAs have been shown to cause catastrophic consequences in feedback control systems while by-passing any residual-based monitoring system. Since these attacks target system measurement process, this paper focuses on the problem of improving the resiliency of dynamical observers against FDIA. Specifically, we propose an attack-resilient pruning algorithm which attempts to exclude compromised channels from being processed by the observer. The proposed pruning algorithm improves attack-localization precision to $100\%$ with high probability, which correspondingly improves the resiliency of the underlying UKF to FDIA. The improvements due to the developed resilient pruning-based observer is validated through a numerical simulation of a two-layer path-tracking control platform of differential-driven wheeled mobile robot (DDWMR) under FDIA.


title: Exponential speedup in measuring out-of-time-ordered correlators with a single bit of quantum information

abstract: Out-of-time-ordered correlators (OTOC) are a quantifier of quantum information scrambling and quantum chaos. We propose an efficient quantum algorithm to measure OTOCs that provides an exponential speed-up over the best known classical algorithm provided the OTOC operator to be estimated admits an efficient gate decomposition. We also discuss a scheme to obtain information about the eigenvalue spectrum and the spectral density of OTOCs.


title: Probabilistic analysis of algorithms for cost constrained minimum weighted combinatorial objects

abstract: We consider cost constrained versions of the minimum spanning tree problem and the assignment problem. We assume edge weights are independent copies of a continuous random variable $Z$ that satisfies $F(x)=\Pr(Z\leq x)\approx x^\alpha$ as $x\to0$, where $\alpha\geq 1$. Also, there are $r=O(1)$ budget constraints with edge costs chosen from the same distribution. We use Lagrangean duality to construct polynomial time algorithms that produce asymptotically optimal solutions. For the spanning tree problem, we allow $r>1$, but for the assignment problem we can only analyse the case $r=1$.


title: Learning Interpretable Feature Context Effects in Discrete Choice

abstract: The outcomes of elections, product sales, and the structure of social connections are all determined by the choices individuals make when presented with a set of options, so understanding the factors that contribute to choice is crucial. Of particular interest are context effects, which occur when the set of available options influences a chooser's relative preferences, as they violate traditional rationality assumptions yet are widespread in practice. However, identifying these effects from observed choices is challenging, often requiring foreknowledge of the effect to be measured. In contrast, we provide a method for the automatic discovery of a broad class of context effects from observed choice data. Our models are easier to train and more flexible than existing models and also yield intuitive, interpretable, and statistically testable context effects. Using our models, we identify new context effects in widely used choice datasets and provide the first analysis of choice set context effects in social network growth.


title: On a conjecture by Anthony Hill

abstract: In the 1950's, English painter Anthony Hill described drawings of complete graphs $K_n$ in the plane having precisely $$H(n) = \tfrac{1}{4}\lfloor \tfrac{n}{2}\rfloor \, \lfloor \tfrac{n-1}{2}\rfloor \, \lfloor \tfrac{n-2}{2}\rfloor \,\lfloor \tfrac{n-3}{2}\rfloor$$ crossings. It became a conjecture that this number is minimum possible and, despite serious efforts, the conjecture is still widely open. Another way of drawing $K_n$ with the same number of crossings was found by Bla\v{z}ek and Koman in 1963. In this note we provide, for the first time, a very general construction of drawings attaining the same bound. Surprisingly, the proof is extremely short and may as well qualify as a "book proof". In particular, it gives a very simple explanation of the phenomenon discovered by Moon in 1968 that a random set of $n$ points on the unit sphere $\SS^2$ in $\RR^3$ joined by geodesics gives rise to a drawing whose number of crossings asymptotically approaches the Hill value $H(n)$.


title: Optimization and Growth in First-Passage Resetting

abstract: We combine the processes of resetting and first-passage to define \emph{first-passage resetting}, where the resetting of a random walk to a fixed position is triggered by a first-passage event of the walk itself. In an infinite domain, first-passage resetting of isotropic diffusion is non-stationary, with the number of resetting events growing with time as $\sqrt{t}$. We calculate the resulting spatial probability distribution of the particle analytically, and also obtain this distribution by a geometric path decomposition. In a finite interval, we define an optimization problem that is controlled by first-passage resetting; this scenario is motivated by reliability theory. The goal is to operate a system close to its maximum capacity without experiencing too many breakdowns. However, when a breakdown occurs the system is reset to its minimal operating point. We define and optimize an objective function that maximizes the reward (being close to maximum operation) minus a penalty for each breakdown. We also investigate extensions of this basic model to include delay after each reset and to two dimensions. Finally, we study the growth dynamics of a domain in which the domain boundary recedes by a specified amount whenever the diffusing particle reaches the boundary after which a resetting event occurs. We determine the growth rate of the domain for the semi-infinite line and the finite interval and find a wide range of behaviors that depend on how much the recession occurs when the particle hits the boundary.


title: A Hybrid Neuro-Symbolic Approach for Complex Event Processing

abstract: Training a model to detect patterns of interrelated events that form situations of interest can be a complex problem: such situations tend to be uncommon, and only sparse data is available. We propose a hybrid neuro-symbolic architecture based on Event Calculus that can perform Complex Event Processing (CEP). It leverages both a neural network to interpret inputs and logical rules that express the pattern of the complex event. Our approach is capable of training with much fewer labelled data than a pure neural network approach, and to learn to classify individual events even when training in an end-to-end manner. We demonstrate this comparing our approach against a pure neural network approach on a dataset based on Urban Sounds 8K.


title: Releasing latent chirality in magnetic two-dimensional materials

abstract: Dzyaloshinskii-Moriya interaction is at heart of chiral magnetism and causes emergence of rich non-collinear and unique topological spin textures in magnetic materials, including cycloids, helices, skyrmions and other. Here we show that strong intrinsic DMI lives in recently discovered van der Waals magnetic two-dimensional materials, due to the sizeable spin-orbit coupling on the non-magnetic ions. In a perfect crystal, this intrinsic DMI remains hidden, but is released with any break of point-inversion symmetry between magnetic ions, unavoidable at the sample edges, at ever present structural defects, with any buckling of the material, or with non-uniform strain on an uneven substrate. We demonstrate such release of latent chirality on an archetypal magnetic monolayer - CrI3, and discuss the plethora of realizable DMI patterns, their control by nanoengineering and tuning by external electric field, thereby opening novel routes in 2D magnetoelectronics.


title: Electrothermal equivalent three-dimensional Finite Element Model of a single neuron

abstract: Objective: We propose a novel approach for modelling the inter-dependence of electrical and mechanical phenomena in nervous cells, by using electro-thermal equivalences in finite element (FE) analysis so that existing thermo-mechanical tools can be applied. Methods: First, the equivalence between electrical and thermal properties of the nerve materials is established, and results of a pure heat conduction analysis performed in Abaqus CAE Software 6.13-3 are validated with analytical solutions for a range of steady and transient conditions. This validation includes the definition of equivalent active membrane properties that enable prediction of the action potential. Then, as a step towards fully coupled models, electromechanical coupling is implemented through the definition of equivalent piezoelectric properties of the nerve membrane using the thermal expansion coefficient, enabling prediction of the mechanical response of the nerve to the action potential. Results: Results of the coupled electro-mechanical model are validated with previously published experimental results of deformation for the squid giant axon, crab nerve fibre and garfish olfactory nerve fibre. Conclusion: A simplified coupled electro-mechanical modelling approach is established through an electro-thermal equivalent FE model of a nervous cell for biomedical applications. Significance: One of the key findings is the mechanical characterization of the neural activity in a coupled electro-mechanical domain, which provides insights into the electromechanical behaviour of nervous cells, such as thinning of the membrane. This is a first step towards modelling 3D electromechanical alteration induced by trauma at nerve bundle, tissue and organ levels.


title: Snowmass 2021 Letter of Interest: Decays of Heavy Flavors Beauty, Charm, and Tau

abstract: The Heavy Flavor Averaging Group provides this Letter of Interest (LOI) as input to the Snowmass 2021 Particle Physics Community Planning Exercise organized by the Division of Particles and Fields of the American Physical Society. Research in heavy flavor physics is an essential component of particle physics, both within and beyond the Standard Model. To fully realize the potential of this field, we advocate strong support within the U.S. high energy physics program for ongoing and future experimental and theory research in heavy flavor physics.


title: Neutral vortex necklace in a trapped planar superfluid

abstract: We study quantum vortex states consisting of a ring of vortices with alternating sign, in a homogeneous superfluid confined to a circular domain. We find an exact stationary solution of the point vortex model for the neutral vortex necklace. We investigate the stability of the necklace state within both the point-vortex model and the Gross-Pitaevskii equation describing a trapped atomic Bose-Einstein condensate at low temperature. The point-vortex stationary states are found to also be stationary states of the Gross-Pitaevskii equation provided the finite thickness of the outer fluid boundary is accounted for. Under significant perturbation, the Gross-Pitaevskii evolution and point-vortex model exhibit instability as expected for metastable states. The perturbed vortex necklace exhibits sensitivity to the perturbation, suggesting a route to seeding vortex chaos or quantum turbulence.


title: Wrinkles in the opening angle method

abstract: We investigate the stability of the deformation modeled by the opening angle method, often used to give a measure of residual stresses in arteries and other biological soft tubular structures. Specifically, we study the influence of stiffness contrast, dimensions and inner pressure on the onset of wrinkles when an open sector of a soft tube, coated with a stiffer film, is bent into a full cylinder. The tube and its coating are made of isotropic, incompressible, hyperelastic materials. We provide a full analytical exposition of the governing equations and the associated boundary value problem for the large deformation and for the superimposed small-amplitude wrinkles. For illustration, we solve them numerically with a robust algorithm in the case of Mooney-Rivlin materials. We confront the results to experimental data that we collected for soft silicone sectors. We study the influence of axial stretch and inner pressure on the stability of closed-up coated tubes with material parameters comparable with those of soft biological tubes such as arteries and veins, although we do not account for anisotropy. We find that the large deformation described in the opening angle method does not always exist, as it can become unstable for certain combinations of dimensions and material parameters.


title: A $K$-rough path above the space-time fractional Brownian motion

abstract: We construct a $K$-rough path above either a space-time or a spatial fractional Brownian motion, in any space dimension $d$. This allows us to provide an interpretation and a unique solution for the corresponding parabolic Anderson model, understood in the renormalized sense. We also consider the case of a spatial fractional noise.


title: Complex valued analytic torsion and dynamical zeta function on locally symmetric spaces

abstract: We show that the Ruelle dynamical zeta function on a closed odd dimensional locally symmetric space twisted by an arbitrary flat vector bundle has a meromorphic extension to the whole complex plane and that its leading term in the Laurent series at the zero point is related to the regularised determinant of the flat Laplacian of Cappell-Miller. When the flat vector bundle is close to an acyclic and unitary one, we show that the dynamical zeta function is regular at the zero point and that its value is equal to the complex valued analytic torsion of Cappell-Miller. This generalises author's previous results for unitarily flat vector bundles as well as M\"uller and Spilioti's results on hyperbolic manifolds.


title: Making a Quantum Universe: Symmetry and Gravity

abstract: So far, none of attempts to quantize gravity has led to a satisfactory model that not only describe gravity in the realm of a quantum world, but also its relation to elementary particles and other fundamental forces. Here, we outline the preliminary results for a model of quantum universe, in which gravity is fundamentally and by construction quantic. The model is based on three well motivated assumptions with compelling observational and theoretical evidence: quantum mechanics is valid at all scales; quantum systems are described by their symmetries; universe has infinite independent degrees of freedom. The last assumption means that the Hilbert space of the Universe has $SU(N\rightarrow \infty) \cong \text{area preserving Diff.} (S_2)$ symmetry, which is parameterized by two angular variables. We show that, in the absence of a background spacetime, this Universe is trivial and static. Nonetheless, quantum fluctuations break the symmetry and divide the Universe to subsystems. When a subsystem is singled out as reference -- observer -- and another as clock, two more continuous parameters arise, which can be interpreted as distance and time. We identify the classical spacetime with parameter space of the Hilbert space of the Universe. Therefore, its quantization is meaningless. In this view, the Einstein equation presents the projection of quantum dynamics in the Hilbert space into its parameter space. Finite dimensional symmetries of elementary particles emerge as a consequence of symmetry breaking when the Universe is divided to subsystems/particles, without having any implication for the infinite dimensional symmetry and its associated interaction - perceived as gravity. This explains why gravity is a universal force.


title: Polariton-assisted donor-acceptor role reversal in resonant energy transfer between organic dyes strongly coupled to electromagnetic modes of a tuneable microcavity

abstract: Resonant interaction between excitonic transitions of molecules and localized electromagnetic field allows the formation of hybrid light-matter polaritonic states. This hybridization of the light and the matter states has been shown to be able to significantly alter the intrinsic properties of molecular ensembles placed inside the optical cavity. Here, we have achieved strong coupling between the excitonic transition in typical oligonucleotide-based molecular beacons labelled with a pair of organic dye molecules, demonstrating an efficient donor to acceptor resonance energy transfer, and the tuneable open-access cavity mode. The photoluminescence of this hybrid system under non-resonant laser excitation and the dependence of the relative population of light-matter hybrid states on cavity detuning have been characterized. Furthermore, by analysing the dependence of the relaxation pathways between energy states in this system, we have demonstrated that predominant strong coupling of the cavity photon to the exciton transition in the donor dye molecule can lead to such a large an energy shift that the energy transfer from the acceptor exciton reservoir to the mainly donor lower polaritonic state can be achieved, thus yielding the chromophores donor-acceptor role reversal or carnival effect. Our experimental data confirm the theoretically predicted possibility for confined electromagnetic fields to control and mediate polariton-assisted remote energy transfer thus paving the way to new approaches to remote-controlled chemistry, energy harvesting, energy transfer and sensing.


title: Combinatorics-Based Approaches to Controllability Characterization for Bilinear Systems

abstract: The control of bilinear systems has attracted considerable attention in the field of systems and control for decades, owing to their prevalence in diverse applications across science and engineering disciplines. Although much work has been conducted on analyzing controllability properties, the mostly used tool remains the Lie algebra rank condition. In this paper, we develop alternative approaches based on theory and techniques in combinatorics to study controllability of bilinear systems. The core idea of our methodology is to represent vector fields of a bilinear system by permutations or graphs, so that Lie brackets are represented by permutation multiplications or graph operations, respectively. Following these representations, we derive combinatorial characterization of controllability for bilinear systems, which consequently provides novel applications of symmetric group and graph theory to control theory. Moreover, the developed combinatorial approaches are compatible with Lie algebra decompositions, including the Cartan and non-intertwining decomposition. This compatibility enables the exploitation of representation theory for analyzing controllability, which allows us to characterize controllability properties of bilinear systems governed by semisimple and reductive Lie algebras.


title: Oblique wrinkles

abstract: We prove theoretically that when a soft solid is subjected to an extreme deformation, wrinkles can form on its surface at an angle that is oblique to a principal direction of stretch. These oblique wrinkles occur for a strain that is smaller than the one required to obtain wrinkles normal to the direction of greatest compression. We go on to explain why they will probably never be observed in real-world experiments.


title: Is Everything Fine, Grandma? Acoustic and Linguistic Modeling for Robust Elderly Speech Emotion Recognition

abstract: Acoustic and linguistic analysis for elderly emotion recognition is an under-studied and challenging research direction, but essential for the creation of digital assistants for the elderly, as well as unobtrusive telemonitoring of elderly in their residences for mental healthcare purposes. This paper presents our contribution to the INTERSPEECH 2020 Computational Paralinguistics Challenge (ComParE) - Elderly Emotion Sub-Challenge, which is comprised of two ternary classification tasks for arousal and valence recognition. We propose a bi-modal framework, where these tasks are modeled using state-of-the-art acoustic and linguistic features, respectively. In this study, we demonstrate that exploiting task-specific dictionaries and resources can boost the performance of linguistic models, when the amount of labeled data is small. Observing a high mismatch between development and test set performances of various models, we also propose alternative training and decision fusion strategies to better estimate and improve the generalization performance.


title: A Distributed Power Control Algorithm for Energy Efficiency Maximization in Wireless Cellular Networks

abstract: In this paper, we propose a distributed power control algorithm for addressing the global energy efficiency (GEE) maximization problem subject to satisfying a minimum target SINR for all user equipments (UEs) in wireless cellular networks. We state the problem as a multi-objective optimization problem which targets minimizing total power consumption and maximizing total throughput, simultaneously, while a minimum target SINR is guaranteed for all UEs. We propose an iterative scheme executed in the UEs to control their transmit power using individual channel state information (CSI) such that the GEE is maximized in a distributed manner. We prove the convergence of the proposed iterative algorithm to its corresponding unique fixed point also shown by our numerical results. Additionally, simulation results demonstrate that our proposed scheme outperforms other algorithms in the literature and performs like the centralized algorithm executed in the base station and maximizes the GEE using the global CSI.


title: A Fast Parametric Ellipse Algorithm

abstract: This paper describes a 2-D graphics algorithm that uses shifts and adds to precisely plot a series of points on an ellipse of any shape and orientation. The algorithm can also plot an elliptic arc that starts and ends at arbitrary angles. The ellipse algorithm described here is largely based on earlier papers by Van Aken and Simar [1,2], which extend Marvin Minsky's well-known circle algorithm [3,4,5] to ellipses, and show how to cancel out the sources of error in Minsky's original algorithm. A new flatness test is presented for automatically controlling the spacing between points plotted on an ellipse or elliptic arc. Most of the calculations performed by the ellipse algorithm and flatness test use fixed-point addition and shift operations, and thus are well-suited to run on less-powerful processors. C++ source code listings are included. Keywords: parametric ellipse algorithm, rotated ellipse, Minsky circle algorithm, flatness, elliptic arc, conjugate diameters, affine invariance


title: Time Evolution and Probability in Quantum Theory: The Central Role of Born's Rule

abstract: In this treatise I introduce the time dependent Generalized Born's Rule for the probabilities of quantum events, including conditional and consecutive probabilities, as the unique fundamental time evolution equation of quantum theory. Then these probabilities, computed from states and events, are to be compared with relative frequencies of observations. Schrodinger's equation still is valid in one model of the axioms of quantum theory, which I call the Schrodinger model. However, the role of Schrodinger's equation is auxiliary, since it serves to help compute the continuous temporal evolution of the probabilities given by the Generalized Born's Rule. In other models, such as the Heisenberg model, the auxiliary equations are quite different, but the Generalized Born's Rule is the same formula (covariance) and gives the same results (invariance). Also some aspects of the Schrodinger model are not found in the isomorphic Heisenberg model, and they therefore do not have any physical significance. One example of this is the infamous collapse of the quantum state. Other quantum phenomena, such as entanglement, are easy to analyze in terms of the Generalized Born's Rule without any reference to the unnecessary concept of collapse. Finally, this leads to the possibility of quantum theory with other sorts of auxiliary equations instead of Schrodinger's equation, and examples of this are given. Throughout this treatise the leit motif is the central importance of quantum probability and most especially of the simplifying role of the time dependent Generalized Born's Rule in quantum theory.


title: Globalization? Trade War? A Counterbalance Perspective

abstract: The embrace of globalization and protectionism among economies has ebbed and flowed over the past few decades. These fluctuations call for quantitative analytics to help countries improve their trade policies. Changing attitudes about globalization also imply that the best trade policies may vary over time and be country-specific. We argue that the imports and exports of all economies constitute a counterbalanced network where conflict and cooperation are two sides of the same coin. Quantitative competitiveness is then formulated for each country using a network counterbalance equilibrium. A country could improve its relative strength in the network by embracing globalization, protectionism, trade collaboration, or conflict. This paper presents the necessary conditions for globalization and trade wars, evaluates their side effects, derives national bargaining powers, identifies appropriate targets for conflict or collaboration, and recommends fair resolutions for trade conflicts. Data and events from the past twenty years support these conditions.


title: Measurement of low-energy cosmic-ray electron and positron spectra at 1 AU with the AESOP-Lite spectrometer

abstract: We report on a new measurement of the cosmic ray (CR) electron and positron spectra in the energy range of 20 MeV -- 1 GeV. The data were taken during the first flight of the balloon-borne spectrometer AESOP-Lite (Anti Electron Sub Orbital Payload), which was flown from Esrange, Sweden, to Ellesmere Island, Canada, in May 2018. The instrument accumulated over 130 hours of exposure at an average altitude of 3 g.cm$^{-2}$ of residual atmosphere. The experiment uses a gas Cherenkov detector and a magnetic spectrometer, consisting of a permanent dipole magnet and silicon strip detectors (SSDs), to identify particle type and measure the rigidity. Electrons and positrons were detected against a background of protons and atmospheric secondary particles. The primary cosmic ray spectra of electrons and positrons, as well as the re-entrant albedo fluxes, were extracted between 20 MeV -- 1 GeV during a positive solar magnetic polarity epoch. The positron fraction below 100 MeV appears flat, suggesting diffusion dominated solar modulation at low rigidity. The all-electron spectrum is presented and compared with models from a heliospheric numerical transport code.

categories: astro-ph.HE astro-ph.IM astro-ph.SR physics.ins-det physics.space-ph


title: On the rate of core collapse supernovae in the Milky Way

abstract: Several large neutrino telescopes, operating at various sites around the world, have as their main objective the first detection of neutrinos emitted by a gravitational collapse in the Milky Way. The success of these observation programs depends on the rate of supernova core collapse in the Milky Way, $R$. In this work, standard statistical techniques are used to combine several independent results. Their consistency is discussed and the most critical input data are identified. The inference on $R$ is further tested and refined by including direct information on the occurrence rate of gravitational collapse events in the Milky Way and in the Local Group, obtained from neutrino telescopes and electromagnetic surveys. A conservative treatment of the errors yields a combined rate $R=1.63 \pm 0.46$ (100 yr)$^{-1}$; the corresponding time between core collapse supernova events turns out to be $T=61_{-14}^{+24}$~yr. The importance to update the analysis of the stellar birthrate method is emphasized.

categories: astro-ph.HE astro-ph.GA astro-ph.SR hep-ex hep-ph


title: ZBLLAC -- A spectroscopic database of BL Lacertae objects

abstract: This paper describes the database of optical spectra of BL Lacertae (BLL) objects ZBLLAC available at the URL https://web.oapd.inaf.it/zbllac/. At the present date, it contains calibrated spectra for 295 BLL. For about 35% of them, we report a firm measure of redshift z while for 35 sources we set a lower limit on z based on the detection of intervening absorption systems, mainly ascribed to Mg II (2800 Ang). We report here on the architecture of the database and on its web front-end that permits to filter, query and interactively explore the data. We discuss some properties of the objects in the present dataset by giving the distribution of the redshifts and reporting on the detected emission lines, which turn out to be mainly forbidden and ascribed to [O II] (3737 Ang) and [O III] (5007 Ang). Finally, we discuss on intervening absorption systems detected in 35 BLLs that allow us to set lower limits to their distance.


title: Non-thermal Sunyaev-Zeldovich signal from radio galaxy lobes

abstract: Energetic electrons in the lobes of radio galaxies make them potential sources for not only radio and X-rays but also Sunyaev-Zeldovich (SZ) distortions in the cosmic microwave background (CMB) radiation. Previous works have discussed the energetics of radio galaxy lobes, but assuming thermal SZ effect, coming from the non-thermal electron population. We use an improved evolutionary model for radio galaxy lobes to estimate the observed parameters such as the radio luminosity and intensity of SZ-distortions at the redshifts of observation. We, further, quantify the effects of various relevant physical parameters of the radio galaxies, such as the jet power, the time scale over which the jet is active, the evolutionary time scale for the lobe, etc on the observed parameters. For current SZ observations towards galaxy clusters, we find that the non-thermal SZ distortions from radio lobes embedded in galaxy clusters can be non-negligible compared to the amount of thermal SZ distortion from the intra-cluster medium and, hence, can not be neglected. We show that small and young (and preferably residing in a cluster environment) radio galaxies offer better prospects for the detection of the non-thermal SZ signal from these sources. We further discuss the limits on different physical parameters for some sources for which SZ effect has been either detected or upper limits are available. The evolutionary models enable us to obtain limits, previously unavailable, on the low energy cut-off of electron spectrum ($p_{min} \sim 1\hbox{--}2$) in order to explain the recent non-thermal SZ detection. Finally, we discuss how future CMB experiments, which would cover higher frequency bands ($>$400 GHz), may provide clear signatures for non-thermal SZ effect.


title: A spacetime outlook on CFD: Spacetime correlated models and spacetime coupled algorithms

abstract: A spacetime outlook on Computational Fluid Dynamics is advocated: models in fluid mechanics often have the spacetime correlation property, which should be inherited and preserved in the corresponding numerical algorithms. Starting from the fundamental formulation of fluid mechanics under continuum hypothesis, this paper defines the meaning of spacetime correlation of the models, establishes the fundamental principle of finite volume schemes, expounds the necessity of spacetime coupling of algorithms, as well as realizes the physical and mathematical unification of basic governing equations of fluid mechanics and finite volume schemes. In practice, the design methodology of spacetime coupling high order numerical algorithms is presented, and the difference from spacetime decoupling method is compared. It should be pointed out that most of the contents in this paper are suitable for computational fluid dynamics under the assumption of continuous medium, and some are only suitable for compressible flow.


title: Machine learning-based method for linearization and error compensation of an absolute rotary encoder

abstract: The main objective of this work is to develop a miniaturized, high accuracy, single-turn absolute, rotary encoder called ASTRAS360. Its measurement principle is based on capturing an image that uniquely identifies the rotation angle. To evaluate this angle, the image first has to be classified into its sector based on its color, and only then can the angle be regressed. In-spired by machine learning, we built a calibration setup, able to generate labeled training data automatically. We used these training data to test, characterize, and compare several machine learning algorithms for the classification and the regression. In an additional experiment, we also characterized the tolerance of our rotary encoder to eccentric mounting. Our findings demonstrate that various algorithms can perform these tasks with high accuracy and reliability; furthermore, providing extra-inputs (e.g. rotation direction) allows the machine learning algorithms to compensate for the mechanical imperfections of the rotary encoder.


title: Ensemble Riemannian Data Assimilation over the Wasserstein Space

abstract: In this paper, we present an ensemble data assimilation paradigm over a Riemannian manifold equipped with the Wasserstein metric. Unlike the Eulerian penalization of error in the Euclidean space, the Wasserstein metric can capture translation and difference between the shapes of square-integrable probability distributions of the background state and observations -- enabling to formally penalize geophysical biases in state-space with non-Gaussian distributions. The new approach is applied to dissipative and chaotic evolutionary dynamics and its potential advantages and limitations are highlighted compared to the classic variational and filtering data assimilation approaches under systematic and random errors.


title: The Effect of a Magnetic Field on the Dynamics of Debris Discs Around White Dwarfs

abstract: Observational estimates of the lifetimes and inferred accretion rates from debris discs around polluted white dwarfs are often inconsistent with the predictions of models of shielded Poynting-Robertson drag on the dust particles in the discs. Moreover, many cool polluted white dwarfs do not show any observational evidence of accompanying discs. This may be explained, in part, if the debris discs had shorter lifetimes and higher accretion rates than predicted by Poynting-Robertson drag alone. We consider the role of a magnetic field on tidally disrupted diamagnetic debris and its subsequent effect on the formation, evolution, and accretion rate of a debris disc. We estimate that magnetic field strengths greater than $\sim$10kG may decrease the time needed for circularisation and the disc lifetimes by several orders of magnitude and increase the associated accretion rates by a similar factor, relative to Poynting-Robertson drag. We suggest some polluted white dwarfs may host magnetic fields below the typical detectable limit and that these fields may account for a proportion of polluted white dwarfs with missing debris discs. We also suggest that diamagnetic drag may account for the higher accretion rate estimates among polluted white dwarfs that cannot be predicted solely by Poynting-Robertson drag and find a dependence on magnetic field strength, orbital pericentre distance, and particle size on predicated disc lifetimes and accretion rates.


title: The Submillimeter Active Region Excess Brightness Temperature during Solar Cycles 23 and 24

abstract: We report the temporal evolution of the excess brightness temperature above solar active regions (ARs) observed with the Solar Submillimeter Telescope (SST) at 212 ({\lambda} = 1.4 mm) and 405 GHz ({\lambda} = 0.7 mm) during Cycles 23 and 24. Comparison with the sunspot number (SSN) yields a Pearson's correlation coefficient R = 0.88 and 0.74 for 212 and 405 GHz, respectively. Moreover, when only Cycle 24 is taken into account the correlation coefficients go to 0.93 and 0.81 for each frequency. We derive the spectral index {\alpha} between SST frequencies and found a slight anti-correlation with the SSN (R = -0.25); however, since the amplitude of the variation is lower than the standard deviation we cannot draw a definite conclusion. Indeed, {\alpha} remains almost constant within the uncertainties with a median value approximate to 0 characteristic of an optically thick thermal source. Since the origin of the AR submillimeter radiation is thermal continuum produced at chromospheric heights, the strong correlation between the excess brightness temperature and the magnetic cycle evolution could be related to the available free magnetic energy to be released in reconnection events.


title: Musical tone coloring via bifurcation control of Eulerian n-tuple Hopf singularities

abstract: An intrinsic essence of sounds in music is the evolution of their qualitative types while in mathematics we interpret each qualitative change by a bifurcation. Hopf bifurcation is an important venue to generate a signal with an arbitrary frequency. Hence, the investigations of musical sounds via bifurcation control theory are long-overdue and natural contributions. In this paper, we address the tone coloring of sounds by dynamical modeling of spectral and temporal envelopes. Multiple number of leading harmonic partials of a note (modulo a hearing sound velocity threshold) are attributed into an Eulerian differential system with n-tuple Hopf singularity. The qualitative evolution of the temporal envelop is then simulated over a set of consecutive time-intervals via bifurcation control of the differential system. For an instance, our proposed approach is applied on audio C#4 files obtained from piano and violin. Fourier analysis is used to generate the amplitude spectral vectors. Then, we associate each amplitude spectral vector with an Eulerian flow-invariant leaf. Bifurcation control suffices to accurately construct the desired spectral and amplitude envelopes of musical notes. These correspond with a rich bifurcation scenarios involving Clifford toral manifolds for the Eulerian differential system. In order to reduce the technicalities, we employ several reduction techniques and use one bifurcation parameter. We show how different ordered sets of elementary bifurcations such as pitchfork and (double) saddle-node bifurcations are associated with the qualitative temporal envelop changes of a C]4 played by either a piano or a violin. A complete hysteresis type cycle is observed within the temporal envelop bifurcations of the C#4 played by violin.


title: Connections between Non-Rotating, Slowly Rotating, and Rapidly Rotating Turbulent Convection Transport Scalings

abstract: Buoyancy-driven convection is likely the dominant driver of turbulent motions in the universe, and thus, is widely studied by physicists, engineers, geophysicists and astrophysicists. Maybe unsurprisingly, these different communities discuss the gross convective behaviors in different ways, often without significant cross-talk existing between them. Here, we seek to draw connections between these communities. We do so by carrying out a set of basic scale estimations for how heat and fluid momentum transport should behave in non-rotating, slowly rotating and rapidly rotating buoyancy-driven convective environments. We find that slowly and rapidly rotating scalings can be inter-related via one parameter, the so-called convective Rossby number $\RoC$, a dissipation-free parameter measuring the importance of buoyancy driving relative to rotation. Further, we map between non-flux-based and the flux-based, buoyancy-driven scalings used by different groups. In doing so, these scalings show that there are clean connections between the different communities' approaches and that a number of the seemingly different scalings are actually synonymous with one another.


title: On Local Kirigami Mechanics I: Isometric Conical Solutions

abstract: Over the past decade, kirigami--the Japanese art of paper cutting--has been playing an increasing role in the emerging field of mechanical metamaterials and a myriad of other mechanical applications. Nonetheless, a deep understanding of the mathematics and mechanics of kirigami structures is yet to be achieved in order to unlock their full potential to pioneer more advanced applications in the field. In this work, we study the most fundamental geometric building block of kirigami: a thin sheet with a single cut. We consider a reduced two-dimensional plate model of a circular thin disk with a radial slit and investigate its deformation following the opening of the slit and the rotation of its lips. In the isometric limit--as the thickness of the disk approaches zero--the elastic energy has no stretching contribution and the thin sheet takes a conical shape known as the e-cone. We solve the post-buckling problem for the e-cone in the geometrically nonlinear setting assuming a Saint Venant-Kirchhoff constitutive plate model; we find closed-form expressions for the stress fields and show the geometry of the e-cone to be governed by the spherical elastica problem. This allows us to fully map out the space of solutions and investigate the stability of the post-buckled e-cone problem assuming mirror symmetric boundary conditions on the rotation of the lips on the open slit.


title: Survival Analysis via Ordinary Differential Equations

abstract: This paper introduces an Ordinary Differential Equation (ODE) notion for survival analysis. The ODE notion not only provides a unified modeling framework, but more importantly, also enables the development of a widely applicable, scalable, and easy-to-implement procedure for estimation and inference. Specifically, the ODE modeling framework unifies many existing survival models, such as the proportional hazards model, the linear transformation model, the accelerated failure time model, and the time-varying coefficient model as special cases. The generality of the proposed framework serves as the foundation of a widely applicable estimation procedure. As an illustrative example, we develop a sieve maximum likelihood estimator for a general semi-parametric class of ODE models. In comparison to existing estimation methods, the proposed procedure has advantages in terms of computational scalability and numerical stability. Moreover, to address unique theoretical challenges induced by the ODE notion, we establish a new general sieve M-theorem for bundled parameters and show that the proposed sieve estimator is consistent and asymptotically normal, and achieves the semi-parametric efficiency bound. The finite sample performance of the proposed estimator is examined in simulation studies and a real-world data example.


title: Covid-19 Belgium: Extended SEIR-QD model with nursing homes and long-term scenarios-based forecasts

abstract: Following the spread of the COVID-19 pandemic and pending the establishment of vaccination campaigns, several non pharmaceutical interventions such as partial and full lockdown, quarantine and measures of physical distancing have been imposed in order to reduce the spread of the disease and to lift the pressure on healthcare system. Mathematical models are important tools for estimating the impact of these interventions, for monitoring the current evolution of the epidemic at a national level and for estimating the potential long-term consequences of relaxation of measures. In this paper, we model the evolution of the COVID-19 epidemic in Belgium with a deterministic age-structured extended compartmental model. Our model takes special consideration for nursing homes which are modelled as separate entities from the general population in order to capture the specific delay and dynamics within these entities. The model integrates social contact data and is fitted on hospitalisations data (admission and discharge), on the daily number of COVID-19 deaths (with a distinction between general population and nursing home related deaths) and results from serological studies, with a sensitivity analysis based on a Bayesian approach. We present the situation as in November 2020 with the estimation of some characteristics of the COVID-19 deduced from the model. We also present several mid-term and long-term projections based on scenarios of reinforcement or relaxation of social contacts for different general sectors, with a lot of uncertainties remaining.


title: Attachment and antibiotic response of early-stage biofilms studied using resonant hyperspectral imaging

abstract: Many bacterial species readily develop biofilms that act as a protective matrix against external challenge, e.g. from antimicrobial treatment. Therefore, biofilms are often responsible for persistent and recurring infections. Established methods for studying biofilms are either destructive or they focus on the biofilm surface. A non-destructive method that is sensitive to the underside of the biofilm is highly desirable, as it allows studying the penetration of antibiotics through the film. Here, we demonstrate that the high surface sensitivity of resonant hyperspectral imaging provides this capability. The method allows us to monitor the early stages of Escherichia coli biofilm formation, cell attachment and microcolony formation, in-situ and in real time. We study the response of the biofilm to a number of different antibiotics and verify our observations using confocal microscopy. Based on this ability to closely monitor the surface-bound cells, resonant hyperspectral imaging gives new insights into the antimicrobial resistance of biofilms.


title: Stable high-order cubature formulas for experimental data

abstract: In many applications, it is impractical -- if not even impossible -- to obtain data to fit a known cubature formula (CF). Instead, experimental data is often acquired at equidistant or even scattered locations. In this work, stable (in the sense of nonnegative only cubature weights) high-order CFs are developed for this purpose. These are based on the approach to allow the number of data points N to be larger than the number of basis functions K which are integrated exactly by the CF. This yields an (N-K)-dimensional affine linear subspace from which cubature weights are selected that minimize certain norms corresponding to stability of the CF. In the process, two novel classes of stable high-order CFs are proposed and carefully investigated.


title: Method to compute the stress-energy tensor for a quantized scalar field when a black hole forms from the collapse of a null shell

abstract: A method is given to compute the stress-energy tensor for a massless minimally coupled scalar field in a spacetime where a black hole forms from the collapse of a spherically symmetric null shell in four dimensions. Part of the method involves matching the modes for the in vacuum state to a complete set of modes in Schwarzschild spacetime. The other part involves subtracting from the unrenormalized expression for the stress-energy tensor when the field is in the in vacuum state, the corresponding expression when the field is in the Unruh state and adding to this the renormalized stress-energy tensor for the field in the Unruh state. The method is shown to work in the two-dimensional case where the results are known.


title: Extraction and Visualization of Poincar\'e Map Topology for Spacecraft Trajectory Design

abstract: Mission designers must study many dynamical models to plan a low-cost spacecraft trajectory that satisfies mission constraints. They routinely use Poincar\'e maps to search for a suitable path through the interconnected web of periodic orbits and invariant manifolds found in multi-body gravitational systems. This paper is concerned with the extraction and interactive visual exploration of this structural landscape to assist spacecraft trajectory planning. We propose algorithmic solutions that address the specific challenges posed by the characterization of the topology in astrodynamics problems and allow for an effective visual analysis of the resulting information. This visualization framework is applied to the circular restricted three-body problem (CR3BP), where it reveals novel periodic orbits with their relevant invariant manifolds in a suitable format for interactive transfer selection. Representative design problems illustrate how spacecraft path planners can leverage our topology visualization to fully exploit the natural dynamics pathways for energy-efficient trajectory designs.


title: Addressing Cold Start in Recommender Systems with Hierarchical Graph Neural Networks

abstract: Recommender systems have become an essential instrument in a wide range of industries to personalize the user experience. A significant issue that has captured both researchers' and industry experts' attention is the cold start problem for new items. In this work, we present a graph neural network recommender system using item hierarchy graphs and a bespoke architecture to handle the cold start case for items. The experimental study on multiple datasets and millions of users and interactions indicates that our method achieves better forecasting quality than the state-of-the-art with a comparable computational time.


title: ePointDA: An End-to-End Simulation-to-Real Domain Adaptation Framework for LiDAR Point Cloud Segmentation

abstract: Due to its robust and precise distance measurements, LiDAR plays an important role in scene understanding for autonomous driving. Training deep neural networks (DNNs) on LiDAR data requires large-scale point-wise annotations, which are time-consuming and expensive to obtain. Instead, simulation-to-real domain adaptation (SRDA) trains a DNN using unlimited synthetic data with automatically generated labels and transfers the learned model to real scenarios. Existing SRDA methods for LiDAR point cloud segmentation mainly employ a multi-stage pipeline and focus on feature-level alignment. They require prior knowledge of real-world statistics and ignore the pixel-level dropout noise gap and the spatial feature gap between different domains. In this paper, we propose a novel end-to-end framework, named ePointDA, to address the above issues. Specifically, ePointDA consists of three modules: self-supervised dropout noise rendering, statistics-invariant and spatially-adaptive feature alignment, and transferable segmentation learning. The joint optimization enables ePointDA to bridge the domain shift at the pixel-level by explicitly rendering dropout noise for synthetic LiDAR and at the feature-level by spatially aligning the features between different domains, without requiring the real-world statistics. Extensive experiments adapting from synthetic GTA-LiDAR to real KITTI and SemanticKITTI demonstrate the superiority of ePointDA for LiDAR point cloud segmentation.


title: Robust Conversational AI with Grounded Text Generation

abstract: This article presents a hybrid approach based on a Grounded Text Generation (GTG) model to building robust task bots at scale. GTG is a hybrid model which uses a large-scale Transformer neural network as its backbone, combined with symbol-manipulation modules for knowledge base inference and prior knowledge encoding, to generate responses grounded in dialog belief state and real-world knowledge for task completion. GTG is pre-trained on large amounts of raw text and human conversational data, and can be fine-tuned to complete a wide range of tasks. The hybrid approach and its variants are being developed simultaneously by multiple research teams. The primary results reported on task-oriented dialog benchmarks are very promising, demonstrating the big potential of this approach. This article provides an overview of this progress and discusses related methods and technologies that can be incorporated for building robust conversational AI systems.


title: Horus: Using Sensor Fusion to Combine Infrastructure and On-board Sensing to Improve Autonomous Vehicle Safety

abstract: Studies predict that demand for autonomous vehicles will increase tenfold between 2019 and 2026. However, recent high-profile accidents have significantly impacted consumer confidence in this technology. The cause for many of these accidents can be traced back to the inability of these vehicles to correctly sense the impending danger. In response, manufacturers have been improving the already extensive on-vehicle sensor packages to ensure that the system always has access to the data necessary to ensure safe navigation. However, these sensor packages only provide a view from the vehicle's perspective and, as a result, autonomous vehicles still require frequent human intervention to ensure safety. To address this issue, I developed a system, called Horus, that combines on-vehicle and infrastructure-based sensors to provide a more complete view of the environment, including areas not visible from the vehicle. I built a small-scale experimental testbed as a proof of concept. My measurements of the impact of sensor failures showed that even short outages (1 second) at slow speeds (25 km/hr scaled velocity) prevents vehicles that rely on on-vehicle sensors from navigating properly. My experiments also showed that Horus dramatically improves driving safety and that the sensor fusion algorithm selected plays a significant role in the quality of the navigation. With just a pair of infrastructure sensors, Horus could tolerate sensors that fail 40% of the time and still navigate safely. These results are a promising first step towards safer autonomous vehicles.


title: Waypoint Following Dynamics of a Quaternion Error Feedback Attitude Control System

abstract: Closed-loop attitude steering can be used to implement a non-standard attitude maneuver by using a conventional attitude control system to track a non-standard attitude profile. The idea has been employed to perform zero-propellant maneuvers on the International Space Station and minimum time maneuvers on NASA's TRACE space telescope. A challenge for operational implementation of the idea is the finite capacity of a space vehicle's command storage buffer. One approach to mitigate the problem is to downsample-and-hold the attitude commands as a set of waypoints for the attitude control system to follow. In this paper, we explore the waypoint following dynamics of a quaternion error feedback control law for downsample-and-hold. It is shown that downsample-and-hold induces a ripple between downsamples that causes the satellite angular rate to significantly overshoot the desired limit. Analysis in the $z$-domain is carried out in order to understand the phenomenon. An interpolating Chebyshev-type filter is proposed that allows attitude commands to be encoded in terms of a set of filter coefficients. Using the interpolating filter, commands can be issued at the ACS rate but with significantly reduced memory requirements. The attitude control system of NASA's Lunar Reconnaissance Orbiter is used as an example to illustrate the behavior of a practical attitude control system.


title: On the distribution of modular square roots of primes

abstract: We use recent bounds on bilinear sums with modular square roots to study the distribution of solutions to congruences $x^2 \equiv p \pmod q$ with primes $p\le P$ and integers $q \le Q$. This can be considered as a combined scenario of Duke, Friedlander and Iwaniec with averaging only over the modulus $q$ and of Dunn, Kerr, Shparlinski and Zaharescu with averaging only over $p$.


title: Wave Composition, Propagation, and Polarization of MHD Turbulence within 0.3AU as Observed by PSP

abstract: Turbulence, a ubiquitous phenomenon in interplanetary space, is crucial for the energy conversion of space plasma at multiple scales. This work focuses on the propagation, polarization and wave composition properties of the solar wind turbulence within 0.3AU, and its variation with heliocentric distances at MHD scales (from 10s to 1000s in the spacecraft frame). We present the probability density function of propagation wavevectors (${\rm{PDF}}(k_\parallel,k_\perp)$) for solar wind turbulence winthin 0.3 AU for the first time: (1) wavevectors cluster quasi-(anti-)parallel to the local background magnetic field for $kd_{\rm i}<0.02$, where $d_{\rm i}$ is the ion inertial length; (2) wavevectors shift to quasi-perpendicular directions for $kd_{\rm i}>0.02$. Based on our wave composition diagnosis, we find that: the outward/anti-sunward Alfv\'en mode dominates over the whole range of scales and distances, the spectral energy density fraction of the inward/sunward fast mode decreases with distance, and the fractional energy densities of the inward and outward slow mode increase with distance. The outward fast mode and inward Alfv\'en mode represent minority populations throughout the explored range of distances and scales. On average, the degree of anisotropy of the magnetic fluctuations defined with respect to the minimum variation direction decreases with increasing scale, with no trend in distance at all scales. Our results provide comprehensive insight into the scenario of transport and transfer of the solar wind fluctuations/turbulence in the inner heliosphere.

categories: physics.space-ph astro-ph.EP astro-ph.SR physics.plasm-ph


title: Qualitative properties of systems of 2 complex homogeneous ODE's: a connection to polygonal billiards

abstract: A correspondence between the orbits of a system of 2 complex, homogeneous, polynomial ordinary differential equations with real coefficients and those of a polygonal billiard is displayed. This correspondence is general, in the sense that it applies to an open set of systems of ordinary differential equations of the specified kind. This allows to transfer results well-known from the theory of polygonal billiards, such as ergodicity, the existence of periodic orbits, the absence of exponential divergence, the existence of additional conservation laws, and the presence of discontinuities in the dynamics, to the corresponding systems of ordinary differential equations. It also shows that the considerable intricacy known to exist for polygonal billiards, also attends these apparently simpler systems of ordinary differential equations.


title: Quantitative and Three-Dimensional Assessment of Holdup Material

abstract: Nuclear material deposited in equipment, transfer lines, and ventilation systems of a processing facility is usually referred to as holdup. In this work, we propose to use an array of detectors co-axial to the inspected pipe to measure the holdup material. This method is implementable into an automated system capable of crawling on surfaces and pipes of various curvatures, which would enable faster, easier, and more accurate holdup safeguards measurements. We first demonstrated that the current holdup assay procedure could lead to a non-negligible bias in the estimate of special nuclear material mass, due to the simplified assumption of deposited geometry introduced by the Generalized Geometry Holdup (GGH) model. The new approach consists of imaging the inner holdup material by characterizing the detector array's response and unfolding it from the measured light output. Our experimental proof of principle consists of three NaI(Tl) detectors surrounding an aluminum pipe containing two cesium-137( 137Cs) sources. We derived the source distribution inside the pipe by first calculating the detector response matrix using a method adaptive to the surface geometry of the object containing the measured holdup material. Creating a matrix of the detector array's measured counts, we then proceed to solve an inverse problem, resulting in an accurately located source position and activity distribution within the response matrix's spatial resolution. We then developed a simulated model of the envisioned experimental setup, which accurately described both the activity and position of the source in 2D. Finally, we extended our model onto a discretized three-dimensional model of the system, encompassing 36 detectors. For the 3D simulation of four different source geometries, the model accurately localized the source position in 3D, while the activity retained a maximum relative error of +-5.32%.


title: Strong coupling constants of doubly heavy baryons with vector mesons in QCD

abstract: Using the most general form of the interpolating current for baryons, the strong electric and magnetic coupling constants of light vector mesons $ \rho $ and $ K^* $ with doubly heavy baryons are computed within the light-cone sum rules. We consider 2- and 3-particle distribution amplitudes of the aforementioned vector mesons. The obtained results can be useful in the analysis of experimental data on the properties of doubly heavy baryons conducted at LHC.


title: LaSOT: A High-quality Large-scale Single Object Tracking Benchmark

abstract: Despite great recent advances in visual tracking, its further development, including both algorithm design and evaluation, is limited due to lack of dedicated large-scale benchmarks. To address this problem, we present LaSOT, a high-quality Large-scale Single Object Tracking benchmark. LaSOT contains a diverse selection of 85 object classes, and offers 1,550 totaling more than 3.87 million frames. Each video frame is carefully and manually annotated with a bounding box. This makes LaSOT, to our knowledge, the largest densely annotated tracking benchmark. Our goal in releasing LaSOT is to provide a dedicated high quality platform for both training and evaluation of trackers. The average video length of LaSOT is around 2,500 frames, where each video contains various challenge factors that exist in real world video footage,such as the targets disappearing and re-appearing. These longer video lengths allow for the assessment of long-term trackers. To take advantage of the close connection between visual appearance and natural language, we provide language specification for each video in LaSOT. We believe such additions will allow for future research to use linguistic features to improve tracking. Two protocols, full-overlap and one-shot, are designated for flexible assessment of trackers. We extensively evaluate 48 baseline trackers on LaSOT with in-depth analysis, and results reveal that there still exists significant room for improvement. The complete benchmark, tracking results as well as analysis are available at http://vision.cs.stonybrook.edu/~lasot/.


title: Structure Constants in equivariant oriented cohomology of flag varieties

abstract: We obtain a formula for structure constants of certain variant form of Bott-Samelson classes for equivariant oriented cohomology of flag varieties. Specializing to singular cohomology/K-theory, we recover formulas of structure constants of Schubert classes of Goldin-Knutson, and that of structure constants of Segre-Schwartz-MacPherson classes of Su. We also obtain a formula for K-theoretic stable basis. Our method comes from the study of formal affine Demazure algebra, so is purely algebraic, while the above mentioned results are geometric.


title: Inverse problems for nonlinear Maxwell's equations with second harmonic generation

abstract: In the current paper we consider an inverse boundary value problem of electromagnetism with nonlinear Second Harmonic Generation (SHG) process. We show the unique determination of the electromagnetic material parameters and the SHG susceptibility parameter of the medium by making electromagnetic measurements on the boundary. We are interested in the case when a frequency is fixed.


title: Quad-Core RSA Processor with Countermeasure Against Power Analysis Attacks

abstract: Rivest-Shamir-Adleman (RSA) cryptosystem uses modular multiplication for encryption and decryption. So, performance of RSA can be drastically improved by optimizing modular multiplication. This paper proposes a new parallel, high-radix Montgomery multiplier for 1024 bits multi-core RSA processor. Each computation step operates in radix 4. The computation speed is increased by more than 4 times. We also implement a True Random Number Generator based resilience block to protect the coprocessor against power attacks.


title: Feedback control of transitional shear flows: Sensor selection for performance recovery

abstract: The choice and placement of sensors and actuators is an essential factor determining the performance that can be realized using feedback control. This determination is especially important, but difficult, in the context of controlling transitional flows. The highly non-normal nature of the linearized Navier-Stokes equations makes the flow sensitive to small perturbations, with potentially drastic performance consequences on closed-loop flow control performance. Full-information controllers, such as the linear quadratic regulator (LQR), have demonstrated some success in reducing transient energy growth and suppressing transition; however, sensor-based output feedback controllers with comparable performance have been difficult to realize. In this study, we propose two methods for sensor selection that enable sensor-based output feedback controllers to recover full-information control performance: one based on a sparse controller synthesis approach, and one based on a balanced truncation procedure for model reduction. Both approaches are investigated within linear and nonlinear simulations of a sub-critical channel flow with blowing and suction actuation at the walls. We find that sensor configurations identified by both approaches allow sensor-based static output feedback LQR controllers to recover full-information LQR control performance, both in reducing transient energy growth and suppressing transition. Further, our results indicate that both the sensor selection methods and the resulting controllers exhibit robustness to Reynolds number variations.


title: $l_p$ regularization for ensemble Kalman inversion

abstract: Ensemble Kalman inversion (EKI) is a derivative-free optimization method that lies between the deterministic and the probabilistic approaches for inverse problems. EKI iterates the Kalman update of ensemble-based Kalman filters, whose ensemble converges to a minimizer of an objective function. EKI regularizes ill-posed problems by restricting the ensemble to the linear span of the initial ensemble, or by iterating regularization with early stopping. Another regularization approach for EKI, Tikhonov EKI, penalizes the objective function using the $l_2$ penalty term, preventing overfitting in the standard EKI. This paper proposes a strategy to implement $l_p, 0<p\leq 1,$ regularization for EKI to recover sparse structures in the solution. The strategy transforms a $l_p$ problem into a $l_2$ problem, which is then solved by Tikhonov EKI. The transformation is explicit, and thus the proposed approach has a computational cost comparable to Tikhonov EKI. We validate the proposed approach's effectiveness and robustness through a suite of numerical experiments, including compressive sensing and subsurface flow inverse problems.


title: Borel sets without perfectly many overlapping translations, III

abstract: We expand the results of Roslanowski and Shelah arXive:1806.06283 , arXive:1909.00937 to all perfect Abelian Polish groups $(H,+)$. In particular, we show that if $\alpha<\omega_1$ and $4\leq k<\omega$, then there is a ccc forcing notion adding a $\Sigma^0_2$ set $B\subseteq H$ which has $\aleph_\alpha$ many pairwise $k$--overlapping translations but not a perfect set of such translations.


title: Shannon entropy estimation for linear processes

abstract: In this paper, we estimate the Shannon entropy $S(f) = -\E[ \log (f(x))]$ of a one-sided linear process with probability density function $f(x)$. We employ the integral estimator $S_n(f)$, which utilizes the standard kernel density estimator $f_n(x)$ of $f(x)$. We show that $S_n (f)$ converges to $S(f)$ almost surely and in $\L^2$ under reasonable conditions.


title: On the Self-Repair Role of Astrocytes in STDP Enabled Unsupervised SNNs

abstract: Neuromorphic computing is emerging to be a disruptive computational paradigm that attempts to emulate various facets of the underlying structure and functionalities of the brain in the algorithm and hardware design of next-generation machine learning platforms. This work goes beyond the focus of current neuromorphic computing architectures on computational models for neuron and synapse to examine other computational units of the biological brain that might contribute to cognition and especially self-repair. We draw inspiration and insights from computational neuroscience regarding functionalities of glial cells and explore their role in the fault-tolerant capacity of Spiking Neural Networks (SNNs) trained in an unsupervised fashion using Spike-Timing Dependent Plasticity (STDP). We characterize the degree of self-repair that can be enabled in such networks with varying degree of faults ranging from 50% - 90% and evaluate our proposal on the MNIST and Fashion-MNIST datasets.


title: Graph Neural Networks for Model Recommendation using Time Series Data

abstract: Time series prediction aims to predict future values to help stakeholders make proper strategic decisions. This problem is relevant in all industries and areas, ranging from financial data to demand to forecast. However, it remains challenging for practitioners to select the appropriate model to use for forecasting tasks. With this in mind, we present a model architecture based on Graph Neural Networks to provide model recommendations for time series forecasting. We validate our approach on three relevant datasets and compare it against more than sixteen techniques. Our study shows that the proposed method performs better than target baselines and state of the art, including meta-learning. The results show the relevancy and suitability of GNN as methods for model recommendations in time series forecasting.


title: Maximal sets of mutually orthogonal frequency squares

abstract: A frequency square is a square matrix in which each row and column is a permutation of the same multiset of symbols. A frequency square is of type $(n;\lambda)$ if it contains $n/\lambda$ symbols, each of which occurs $\lambda$ times per row and $\lambda$ times per column. In the case when $\lambda=n/2$ we refer to the frequency square as binary. A set of $k$-MOFS$(n;\lambda)$ is a set of $k$ frequency squares of type $(n;\lambda)$ such that when any two of the frequency squares are superimposed, each possible ordered pair occurs equally often. A set of $k$-maxMOFS$(n;\lambda)$ is a set of $k$-MOFS$(n;\lambda)$ that is not contained in any set of $(k+1)$-MOFS$(n;\lambda)$. For even $n$, let $\mu(n)$ be the smallest $k$ such that there exists a set of $k$-maxMOFS$(n;n/2)$. It was shown in [Electron. J. Combin. 27(3) (2020), P3.7] that $\mu(n)=1$ if $n/2$ is odd and $\mu(n)>1$ if $n/2$ is even. Extending this result, we show that if $n/2$ is even, then $\mu(n)>2$. Also, we show that whenever $n$ is divisible by a particular function of $k$, there does not exist a set of $k'$-maxMOFS$(n;n/2)$ for any $k'\le k$. In particular, this means that $\limsup \mu(n)$ is unbounded. Nevertheless we can construct infinite families of maximal binary MOFS of fixed cardinality. More generally, let $q=p^u$ be a prime power and let $p^v$ be the highest power of $p$ that divides $n$. If $0\le v-uh<u/2$ for $h\ge1$ then we show that there exists a set of $(q^h-1)^2/(q-1)$-maxMOFS$(n;n/q)$.


title: Sky localization of space-based gravitational wave detectors

abstract: Localizing the sky position of the gravitational wave source is a key scientific goal for gravitational wave observations. Employing the Fisher information matrix approximation, we compute the angular resolutions of LISA and TianQin, two planned space-based gravitational wave detectors and examine how detectors' configuration properties, such as the orientation change of the detector plane, heliocentric or geocentric motion and the arm length etc. affect the accuracy of source localization. We find that the amplitude modulation due to the annual changing orientation of the detector plane helps LISA get better accuracy in the sky localization and better sky coverage at frequencies below several mHz, and its effect on TianQin is negligible although the orientation of TianQin's detector plane is fixed. At frequencies above roughly 30mHz, TianQin's ability in the sky localization is better than LISA. Further we explore potential space detector networks for fast and accurate localization of the gravitational wave sources. The LISA-TianQin network has better ability in sky localization for sources with frequencies in the range 1-100 mHz and the network has larger sky coverage for the angular resolution than the individual detector.


title: A Residual Solver and Its Unfolding Neural Network for Total Variation Regularized Models

abstract: This paper proposes to solve the Total Variation regularized models by finding the residual between the input and the unknown optimal solution. After analyzing a previous method, we developed a new iterative algorithm, named as Residual Solver, which implicitly solves the model in gradient domain. We theoretically prove the uniqueness of the gradient field in our algorithm. We further numerically confirm that the residual solver can reach the same global optimal solutions as the classical method on 500 natural images. Moreover, we unfold our iterative algorithm into a convolution neural network (named as Residual Solver Network). This network is unsupervised and can be considered as an "enhanced version" of our iterative algorithm. Finally, both the proposed algorithm and neural network are successfully applied on several problems to demonstrate their effectiveness and efficiency, including image smoothing, denoising, and biomedical image reconstruction. The proposed network is general and can be applied to solve other total variation regularized models.


title: A new route toward orthogonality

abstract: We revisit the problem of determining conditions under which a pure state, that evolves under an arbitrary unitary transformation, reaches an orthogonal state in a finite amount of the transformation parameter. Simple geometric considerations disclose the existence of a fundamental limit for the minimal amount required, providing, in particular, an intuitive hint of the Mandelstam-Tamm bound. The geometric considerations leads us to focus on a particular, yet relevant, family of states that evolve towards orthogonality. Several dynamical features are discussed, which include the (relative) entropy production during transformation, and special attention is paid to multipartite systems of $N$ bosons that are allowed to tunnel between two sites. The effects of the tunneling in the amount of transformation required for the system to attain an orthogonal state are revealed, and the relation between the latter, the tunneling intensity and the mode-entanglement is explored.


title: Higher Order Apostol-Type Poly-Genocchi Polynomials with Parameters a, b and c

abstract: In this paper, a new form of poly-Genocchi polynomials is defined by means of poly-logarithm, namely, the Apostol-type poly-Genocchi polynomials of higher order with parameters a, b and c. Several properties of these polynomials are established including some recurrence relations and explicit formulas, which express these higher order Apostol-type poly-Genocchi polynomials in terms of Stirling numbers of the second kind, Apostol-type Bernoulli and Frobenius polynomials of higher order. Moreover, certain differential identity is obtained that leads this new form of poly-Genocchi polynomials to be classified as Appell polynomials and, consequently, draw more properties using some theorems on Appell polynomials. Furthermore, a symmetrized generalization of this new form of poly-Genocchi polynomials is introduced that possesses a double generating function. Finally, the type 2 Apostol-poly-Genocchi polynomials with parameters a, b and c are defined using the concept of polyexponential function and several identities are derived, two of which show the connections of these polynomials with Stirling numbers of the first kind and the type 2 Apostol-type poly-Bernoulli polynomials.


title: Revisiting radiative leptonic $B$ decay

abstract: In this paper, we summarize the existing methods of solving the evolution equation of the leading-twist $B$-meson LCDA. Then, in the Mellin space, we derive a factorization formula with next-to-leading-logarithmic (NLL) resummation for the form factors $F_{A,V}$ in the $B \to \gamma \ell\nu$ decay at leading power in $\Lambda/m_b$. Furthermore, we investigate the power suppressed local contributions, factorizable non-local contributions (which are suppressed by $1/E_\gamma$ and $1/m_b$), and soft contributions to the form factors. In the numerical analysis, which employs the two-loop-level hard function and the jet function, we find that both the resummation effect and the power corrections can sizably decrease the form factors. Finally, the integrated branching ratios are also calculated for comparison with future experimental data.


title: Some Formulae of Genocchi Polynomials of Higher Order

abstract: In this paper, some formulae for Genoochi polynomials of higher order are derived using the fact that sets of Bernoulli and Euler polynomials of higher order form basis for the polynomial space.


abstract: Quantization of the parameters of machine learning models, such as deep neural networks, requires solving constrained optimization problems, where the constraint set is formed by the Cartesian product of many simple discrete sets. For such optimization problems, we study the performance of the Alternating Direction Method of Multipliers for Quantization ($\texttt{ADMM-Q}$) algorithm, which is a variant of the widely-used ADMM method applied to our discrete optimization problem. We establish the convergence of the iterates of $\texttt{ADMM-Q}$ to certain $\textit{stationary points}$. To the best of our knowledge, this is the first analysis of an ADMM-type method for problems with discrete variables/constraints. Based on our theoretical insights, we develop a few variants of $\texttt{ADMM-Q}$ that can handle inexact update rules, and have improved performance via the use of "soft projection" and "injecting randomness to the algorithm". We empirically evaluate the efficacy of our proposed approaches.


title: The inverse problem for a spectral asymmetry function of the Schr\"odinger operator on a finite interval

abstract: For the Schr\"odinger equation $-d^2 u/dx^2 + q(x)u = \lambda u$ on a finite $x$-interval, there is defined an "asymmetry function" $a(\lambda;q)$, which is entire of order $1/2$ and type $1$ in $\lambda$. Our main result identifies the classes of square-integrable potentials $q(x)$ that possess a common asymmetry function. For any given $a(\lambda)$, there is one potential for each Dirichlet spectral sequence.


title: Fiber cones of rational normal scrolls are Cohen-Macaulay

abstract: In this short paper, we show that the fiber cones of rational normal scrolls are Cohen-Macaulay. As an application, we compute their Castelnuovo-Mumford regularities and $\mathbf{a}$-invariants, as well as the reduction number of the defining ideals of the rational normal scrolls. We also characterize the Gorensteinness of the fiber cone.


title: Prenex normal form theorems in semi-classical arithmetic

abstract: Akama et al. systematically studied an arithmetical hierarchy of the law of excluded middle and related principles in the context of first-order arithmetic. In that paper, they first provide a prenex normal form theorem as a justification of their semi-classical principles restricted to prenex formulas. However, there are some errors in their proof. In this paper, we provide a simple counterexample of their prenex normal form theorem, then modify it in an appropriate way. In addition, we characterize several prenex normal form theorems with respect to semi-classical arithmetic.


title: On principal types and well-foundedness of the cummulativity relation in ECC

abstract: When we investigate a type system, it is helpful if we can establish the well-foundedness of types or terms with respect to a certain hierarchy, and the Extended Calculus of Constructions (called $ECC$, defined and studied comprehensively in [Luo,1994]) is no exception. However, under a very natural hierarchy relation (called the cumulativity relation in [Luo,1994]), the well-foundedness of the hierarchy does not hold generally. In this article,we show that the cumulativity relation is well-founded if it is restricted to one of the following two natural families of terms: \begin{enumerate} \item types in a valid context \item terms having normal forms \end{enumerate} Also, we give an independent proof of the existence of principal types in $ECC$ since it is used in the proof of well-foundedness of cumulativity relation in a valid context although it is often proved by utilizing the well-foundedness of the hierarchy, which would make our argument circular if adopted.


title: Null Lagrangians in linear theories of micropolar type and few other generalizations of elasticity

abstract: In the context of linear theories of generalized elasticity including those for homogeneous micropolar media, quasicrystals, piezoelectric and piezomagnetic media, we explore the concept of null Lagrangians. For obtaining the family of null Lagrangians we employ the sufficient conditions of H. Rund. In some cases a non-zero null Lagrangian is found and the stored energy admits a split into a null Lagrangian and a remainder. However, the null Lagrangian vanishes whenever the relevant elasticity tensor obeys certain symmetry conditions which can be construed as an analogue of the Cauchy relations.


title: Adversarial Attack on Large Scale Graph

abstract: Recent studies have shown that graph neural networks (GNNs) are vulnerable against perturbations due to lack of robustness and can therefore be easily fooled. Currently, most works on attacking GNNs are mainly using gradient information to guide the attack and achieve outstanding performance. However, the high complexity of time and space makes them unmanageable for large scale graphs and becomes the major bottleneck that prevents the practical usage. We argue that the main reason is that they have to use the whole graph for attacks, resulting in the increasing time and space complexity as the data scale grows. In this work, we propose an efficient Simplified Gradient-based Attack (SGA) method to bridge this gap. SGA can cause the GNNs to misclassify specific target nodes through a multi-stage attack framework, which needs only a much smaller subgraph. In addition, we present a practical metric named Degree Assortativity Change (DAC) to measure the impacts of adversarial attacks on graph data. We evaluate our attack method on four real-world graph networks by attacking several commonly used GNNs. The experimental results demonstrate that SGA can achieve significant time and memory efficiency improvements while maintaining competitive attack performance compared to state-of-art attack techniques. Codes are available via: https://github.com/EdisonLeeeee/SGAttack.


title: 2-Group Symmetries of 6d Little String Theories and T-duality

abstract: We determine the 2-group structure constants for all the six-dimensional little string theories (LSTs) geometrically engineered in F-theory without frozen singularities. We use this result as a consistency check for T-duality: the 2- groups of a pair of T-dual LSTs have to match. When the T-duality involves a discrete symmetry twist the 2-group used in the matching is modified. We demonstrate the matching of the 2-groups in several examples.


title: Null Lagrangians in Cosserat elasticity

abstract: In the framework of nonlinear theory of Cosserat elasticity, also called micropolar elasticity, we provide the complete characterization of null Lagrangians for three dimensional bodies as well as for shells. Using the Gibb's rotation vector for description of the microrotation, this task is possible by an application of a theorem stated by Olver and Sivaloganathan in `{the structure of null Lagrangians}' (Nonlinearity, {1}, 1988, pp. 389-398). A set of necessary and sufficient conditions is also provided for the elasticity tensors to correspond to a null Lagrangian in linearized micropolar theory.

categories: math-ph math.DG math.MP physics.class-ph


title: Zone-folded longnitude acoustic phonons driving self-trapped state emission in colloidal CdSe nanoplate superlattice

abstract: Colloidal cadmium chalcogenide nanoplates are two-dimensional semiconductors that have shown great application prospect for light-emitting technologies. Self-trapped state (STS), a special localized state originated from strong electron-phonon coupling (EPC), has great potential in one-step white light luminance owing to its broadband emission linewidth. However, achieving STS in cadmium chalcogenide nanocrystals is extremely challenging due to their intrinic weak EPC nature. By building hybrid superlattice (SL) structures via self-assembly of colloidal CdSe nanoplates (NPLs), we demonstrated an emergence of zone-folded longnitude acoustic phonons (ZFLAP) differ from monodispersed NPLs, and observed a broadband STS emission in spectra range of 450-600 nm. Through femtosecond transient absorption and impulsive vibrational spectroscopy, we revealed that STS is generated in time scale of ~500 fs and is driven by strong coupling of excitons and ZFLAPs with Huang-Rhys parameter as large as ~22.7. Our findings provide a new avenue for generating and manipulating STS emission by artificially designing and building hybrid periodic structures superior to single material optimization.


title: Integration and geometrization of Rota-Baxter Lie algebras

abstract: This paper first introduces the notion of a Rota-Baxter operator (of weight $1$) on a Lie group so that its differentiation gives a Rota-Baxter operator on the corresponding Lie algebra. Direct products of Lie groups, including the decompositions of Iwasawa and Langlands, carry natural Rota-Baxter operators. Formal inverse of the Rota-Baxter operator on a Lie group is precisely the crossed homomorphism on the Lie group, whose tangent map is the differential operator of weight $1$ on a Lie algebra. A factorization theorem of Rota-Baxter Lie groups is proved, deriving directly on the Lie group level, the well-known global factorization theorems of Semenov-Tian-Shansky in his study of integrable systems. As geometrization, the notions of Rota-Baxter Lie algebroids and Rota-Baxter Lie groupoids are introduced, with the former a differentiation of the latter. Further, a Rota-Baxter Lie algebroid naturally gives rise to a post-Lie algebroid, generalizing the well-known fact for Rota-Baxter Lie algebras and post-Lie algebras. It is shown that the geometrization of a Rota-Baxter Lie algebra or a Rota-Baxter Lie group can be realized by its action on a manifold. Examples and applications are provided for these new notions.


title: Quasiperiodic patterns of the complex dimensions of nonlattice self-similar strings, via the LLL algorithm

abstract: The Lattice String Approximation algorithm (or LSA algorithm) of M. L. Lapidus and M. van Frankenhuijsen is a procedure that approximates the complex dimensions of a nonlattice self-similar fractal string by the complex dimensions of a lattice self-similar fractal string. The implication of this procedure is that the set of complex dimensions of a nonlattice string has a quasiperiodic pattern. Using the LSA algorithm, together with the multiprecision polynomial solver MPSolve which is due to D. A. Bini, G. Fiorentino and L. Robol, we give a new and significantly more powerful presentation of the quasiperiodic patterns of the sets of complex dimensions of nonlattice self-similar fractal strings. The implementation of this algorithm requires a practical method for generating simultaneous Diophantine approximations, which in some cases we can accomplish by the continued fraction process. Otherwise, as was suggested by Lapidus and van Frankenhuijsen, we use the LLL algorithm of A. K. Lenstra, H. W. Lenstra, and L. Lov\'asz.


title: High order finite difference Hermite WENO fast sweeping methods for static Hamilton-Jacobi equations

abstract: In this paper, we propose a novel Hermite weighted essentially non-oscillatory (HWENO) fast sweeping method to solve the static Hamilton-Jacobi equations efficiently. During the HWENO reconstruction procedure, the proposed method is built upon a new finite difference fifth order HWENO scheme involving one big stencil and two small stencils. However, one major novelty and difference from the traditional HWENO framework lies in the fact that, we do not need to introduce and solve any additional equations to update the derivatives of the unknown function $\phi$. Instead, we use the current $\phi$ and the old spatial derivative of $\phi$ to update them. The traditional HWENO fast sweeping method is also introduced in this paper for comparison, where additional equations governing the spatial derivatives of $\phi$ are introduced. The novel HWENO fast sweeping methods are shown to yield great savings in both computational time and storage, which improves the computational efficiency of the traditional HWENO scheme. In addition, a hybrid strategy is also introduced to further reduce computational costs. Extensive numerical experiments are provided to validate the accuracy and efficiency of the proposed approaches.


title: PT-symmetric topological near-zero interface state

abstract: Photonic systems with parity-time (PT) symmetry and topology are attracting considerable attentions. In this work, topological near-zero edge states are studied in PT-symmetric photonic lattice and the results indicate that the near-zero edge states can be broken spontaneously in spite of the unbroken PT symmetry. To achieve the stable topological near-zero mode, a binary lattice with carefully designed PT-symmetric is proposed. Further study shows such a structure supports a stable topological interface state experiences phase transition similar to the bulk states in infinite lattice and thus possess real-eigenvalues even with unbroken PT phase. Our study enriches the content of non-Hermitian topological physics and might have potential applications in the fields of topological lasing and quantum computation.


title: First encounters on Bethe Lattices and Cayley Trees

abstract: In this work we consider the first encounter problems between a fixed and/or mobile target A and a moving trap B on Bethe Lattices and Cayley trees. The survival probability (SP) of the target A on the both kinds of structures are analyzed analytically and compared. On Bethe Lattices, the results show that the fixed target will still prolong its survival time, whereas, on Cayley trees, there are some initial positions where the target should move to prolong its survival time. The mean first encounter time (MFET) for mobile target A is evaluated numerically and compared with the mean first passage time (MFPT) for the fixed target A. Different initial settings are addressed and clear boundaries are obtained. These findings are helpful for optimizing the strategy to prolong the survival time of the target or to speed up the search process on Cayley trees, in relation to the target's movement and the initial position configuration of the two walkers. We also present a new method, which uses a small amount of memory, for simulating random walks on Cayley trees.