title: South American auroral reports during the Carrington storm

abstract: The importance of the investigation of magnetic superstorms is not limited to academic interest, because these superstorms can cause catastrophic impact on the modern civilisation due to our increasing dependency on technological infrastructure. In this context, the Carrington storm in September 1859 is considered as a benchmark of observational history owing to its magnetic disturbance and equatorial extent of the auroral oval. So far, several recent auroral reports at that time have been published but those reports are mainly derived from the Northern Hemisphere. In this study, we analyse datable auroral reports from South America and its vicinity, assess the auroral extent using philological and astrometric approaches, identify the auroral visibility at - 17.3{\deg} magnetic latitude and further poleward and reconstruct the equatorial boundary of the auroral oval to be 25.1{\deg} +/- 0.5{\deg} in invariant latitude. Interestingly, brighter and more colourful auroral displays were reported in the South American sector than in the Northern Hemisphere. This north-south asymmetry is presumably associated with variations of their magnetic longitude and the weaker magnetic field over South America compared to the magnetic conjugate point and the increased amount of magnetospheric electron precipitation into the upper atmosphere. These results attest that the magnitude of the Carrington storm indicates that its extent falls within the range of other superstorms, such as those that occurred in May 1921 and February 1872, in terms of the equatorial boundary of the auroral oval.


title: Projectivity in (bounded) integral residuated lattices

abstract: In this paper we study projective algebras in varieties of (bounded) commutative integral residuated lattices from an algebraic (as opposed to categorical) point of view. In particular we use a well-established construction in residuated lattices: the ordinal sum. Its interaction with divisibility makes our results have a better scope in varieties of divisibile commutative integral residuated lattices, and it allows us to show that many such varieties have the property that every finitely presented algebra is projective. In particular, we obtain results on (Stonean) Heyting algebras, certain varieties of hoops, and product algebras. Moreover, we study varieties with a Boolean retraction term, showing for instance that in a variety with a Boolean retraction term all finite Boolean algebras are projective. Finally, we connect our results with the theory of Unification.


title: Ion-pair dissociation dynamics in electron collision with carbon dioxide probed by velocity slice imaging

abstract: Ion-pair dissociation (IPD) to gas phase carbon dioxide molecule has been studied using time of flight (TOF) based mass spectroscopy in combination with the highly differential velocity slice imaging (VSI) technique. The appearance energy of the fragmented anion provides the experimental threshold energy value for ion-pair production. The kinetic energy (KE) distributions and angular distributions (AD) of the fragment anion dispense the detailed insight into the IPD dynamics. The KE distribution clearly reveals that the IPD dynamics may be due to the direct access to the ion-pair states. However, indirect mechanism can't be ruled out at higher incident electron energies. The angular distribution data unambiguously identified the involvement of the ion-pair state associated with Sigma symmetry and a minor contribution from Pi symmetric states. Computational calculations using density functional theory (DFT) strongly support the experimental observations.


title: A New Approach to Generalised Uncertainty Relations

abstract: We outline a new model in which generalised uncertainty relations are obtained without modified commutation relations. While existing models introduce modified phase space volumes for the canonical degrees of freedom, we introduce new degrees of freedom for the background geometry. The phase space is therefore enlarged but remains Euclidean. The spatial background is treated as a genuinely quantum object, with an associated state vector, and the model naturally gives rise to the extended generalised uncertainty principle (EGUP). Importantly, this approach solves (or rather, evades) well known problems associated with modified commutators, including violation of the equivalence principle, the `soccer ball' problem for multiparticle states, and the velocity dependence of the minimum length. However, it implies two radical conclusions. The first is that space must be quantised on a different scale to matter and the second is that the fundamental quanta of geometry are fermions. We explain how, in the context of the model, these do not contradict established results including the no go theorems for multiple quantisation constants, which still hold for species of material particles, and the spin-$2$ nature of gravitons.


title: Enhancement of low-mass dileptons in ultraperipheral collisions

abstract: It is shown that production of low-mass $e^+e^-$-pairs in ultraperipheral nuclear collisions is enhanced due to the Sommerfeld-Gamow-Sakharov (SGS) factor. This effect is especially strong near the threshold of creation of unbound $e^+e^-$-pairs with low masses in the two-photon fusion. Coulomb attraction of the non-relativistic components of such pairs may lead to the increased intensity of 511 keV photons. It can be recorded at the NICA collider and has some astrophysical implications. The analogous effect can be observed at LHC in dilepton production.

categories: hep-ph astro-ph.HE hep-ex nucl-ex nucl-th


title: Retarded resonance Casimir-Polder interaction of a uniformly rotating two-atom system

abstract: We consider here, a two-atom system is uniformly moving through a circular ring at an ultra-relativistic speed and weakly interacting with common external fields. The vacuum fluctuations of the quantum fields generate the entanglement between the atoms. Hence an effective energy shift is originated, which depends on the inter-atomic distance. This is commonly known as resonance Casimir-Polder interaction (RCPI). It is well known that, for a linearly accelerated system coupled with a massless scalar field, we get a thermal response when the local inertial approximation is valid. On the contrary, the non-thermality arises in the presence of the centripetal acceleration. We use the quantum master equation formalism to calculate the second-order energy shift of the entangled states in the presence of two kinds of fields. They are the massive free scalar field and the electromagnetic vector field. For both cases, we observe the non-thermal behavior. A unique retarded response is also noticed in comparison to the free massless case, which can be observed via the polarization transfer technique.


title: Nonequilibrium Theory of Photoinduced Valley Hall Effect

abstract: A recent scientific debate has arisen: Which processes underlie the actual ground of the valley Hall effect (VHE) in two-dimensional materials? The original VHE emerges in samples with ballistic transport of electrons due to the anomalous velocity terms resulting from the Berry phase effect. In disordered samples though, alternative mechanisms associated with electron scattering off impurities have been suggested: (i) asymmetric electron scattering, called skew scattering, and (ii) a shift of the electron wave packet in real space, called a side jump. It has been claimed that the side jump not only contributes to the VHE but fully offsets the anomalous terms regardless of the drag force for fundamental reasons and, thus, the side-jump together with skew scattering become the dominant mechanisms. However, this claim is based on equilibrium theories without any external valley-selective optical pumping, which makes the results fundamentally interesting but incomplete and impracticable. We develop in this paper a microscopic theory of the photoinduced VHE using the Keldysh nonequilibrium diagrammatic technique and find that the asymmetric skew scattering mechanism is dominant in the vicinity of the interband absorption edge. This allows us to explain the operation of optical transistors based on the VHE.


title: A Novel Training Protocol for Performance Predictors of Evolutionary Neural Architecture Search Algorithms

abstract: Evolutionary Neural Architecture Search (ENAS) can automatically design the architectures of Deep Neural Networks (DNNs) using evolutionary computation algorithms. However, most ENAS algorithms require intensive computational resource, which is not necessarily available to the users interested. Performance predictors are a type of regression models which can assist to accomplish the search, while without exerting much computational resource. Despite various performance predictors have been designed, they employ the same training protocol to build the regression models: 1) sampling a set of DNNs with performance as the training dataset, 2) training the model with the mean square error criterion, and 3) predicting the performance of DNNs newly generated during the ENAS. In this paper, we point out that the three steps constituting the training protocol are not well though-out through intuitive and illustrative examples. Furthermore, we propose a new training protocol to address these issues, consisting of designing a pairwise ranking indicator to construct the training target, proposing to use the logistic regression to fit the training samples, and developing a differential method to building the training instances. To verify the effectiveness of the proposed training protocol, four widely used regression models in the field of machine learning have been chosen to perform the comparisons on two benchmark datasets. The experimental results of all the comparisons demonstrate that the proposed training protocol can significantly improve the performance prediction accuracy against the traditional training protocols.


title: Effective viscosity of random suspensions without uniform separation

abstract: This work is devoted to the definition and the analysis of the effective viscosity associated with a random suspension of small rigid particles in a steady Stokes fluid. While previous works on the topic have been conveniently assuming that particles are uniformly separated, we relax this restrictive assumption in form of mild moment bounds on interparticle distances.


title: Performance Analysis of the Gaussian Quasi-Maximum Likelihood Approach for Independent Vector Analysis

abstract: Maximum Likelihood (ML) estimation requires precise knowledge of the underlying statistical model. In Quasi ML (QML), a presumed model is used as a substitute to the (unknown) true model. In the context of Independent Vector Analysis (IVA), we consider the Gaussian QML Estimate (QMLE) of the demixing matrices set and present an (approximate) analysis of its asymptotic separation performance. In Gaussian QML the sources are presumed to be Gaussian, with covariance matrices specified by some "educated guess". The resulting quasi-likelihood equations of the demixing matrices take a special form, recently termed an extended "Sequentially Drilled" Joint Congruence (SeDJoCo) transformation, which is reminiscent of (though essentially different from) classical joint diagonalization. We show that asymptotically this QMLE, i.e., the solution of the resulting extended SeDJoCo transformation, attains perfect separation (under some mild conditions) regardless of the sources' true distributions and/or covariance matrices. In addition, based on the "small-errors" assumption, we present a first-order perturbation analysis of the extended SeDJoCo solution. Using the resulting closed-form expressions for the errors in the solution matrices, we provide closed-form expressions for the resulting Interference-to-Source Ratios (ISRs) for IVA. Moreover, we prove that asymptotically the ISRs depend only on the sources' covariances, and not on their specific distributions. As an immediate consequence of this result, we provide an asymptotically attainable lower bound on the resulting ISRs. We also present empirical results, corroborating our analytical derivations, of three simulation experiments concerning two possible model errors - inaccurate covariance matrices and sources' distribution mismodeling.


title: Liouville type theorems and periodic solutions for the nonhomogeneous parabolic systems

abstract: In the present paper we derive Liouville type results and existence of periodic solutions for $\chi^{(2)}$ type systems with non-homogeneous nonlinearities. Moreover, we prove both universal bounds as well as singularity and decay estimates for this class of problems. In this study, we have to face new difficulties due to the non-homogenous nonlinearities. To overcome this issue, we carry out delicate integral estimates for this class of nonlinearities and modify the usual scaling and blow up arguments. This seems to be the first result for parabolic systems with non-homogeneous nonlinearities.


title: Caching Transient Content for IoT Sensing: Multi-Agent Soft Actor-Critic

abstract: Edge nodes (ENs) in Internet of Things commonly serve as gateways to cache sensing data while providing accessing services for data consumers. This paper considers multiple ENs that cache sensing data under the coordination of the cloud. Particularly, each EN can fetch content generated by sensors within its coverage, which can be uploaded to the cloud via fronthaul and then be delivered to other ENs beyond the communication range. However, sensing data are usually transient with time whereas frequent cache updates could lead to considerable energy consumption at sensors and fronthaul traffic loads. Therefore, we adopt age of information to evaluate data freshness and investigate intelligent caching policies to preserve data freshness while reducing cache update costs. Specifically, we model the cache update problem as a cooperative multi-agent Markov decision process with the goal of minimizing the long-term average weighted cost. To efficiently handle the exponentially large number of actions, we devise a novel reinforcement learning approach, which is a discrete multi-agent variant of soft actor-critic (SAC). Furthermore, we generalize the proposed approach into a decentralized control, where each EN can make decisions based on local observations only. Simulation results demonstrate the superior performance of the proposed SAC-based caching schemes.


title: Neural Codes With Three Maximal Codewords: Convexity and Minimal Embedding Dimension

abstract: Neural codes, represented as collections of binary strings called codewords, are used to encode neural activity. A code is called convex if its codewords are represented as an arrangement of convex open sets in Euclidean space. Previous work has focused on addressing the question: how can we tell when a neural code is convex? Giusti and Itskov identified a local obstruction and proved that convex neural codes have no local obstructions. The converse is true for codes on up to four neurons, but false in general. Nevertheless, we prove this converse holds for codes with up to three maximal codewords, and moreover the minimal embedding dimension of such codes is at most two.


title: Learn by Observation: Imitation Learning for Drone Patrolling from Videos of A Human Navigator

abstract: We present an imitation learning method for autonomous drone patrolling based only on raw videos. Different from previous methods, we propose to let the drone learn patrolling in the air by observing and imitating how a human navigator does it on the ground. The observation process enables the automatic collection and annotation of data using inter-frame geometric consistency, resulting in less manual effort and high accuracy. Then a newly designed neural network is trained based on the annotated data to predict appropriate directions and translations for the drone to patrol in a lane-keeping manner as humans. Our method allows the drone to fly at a high altitude with a broad view and low risk. It can also detect all accessible directions at crossroads and further carry out the integration of available user instructions and autonomous patrolling control commands. Extensive experiments are conducted to demonstrate the accuracy of the proposed imitating learning process as well as the reliability of the holistic system for autonomous drone navigation. The codes, datasets as well as video demonstrations are available at https://vsislab.github.io/uavpatrol


title: IBP reduction coefficients made simple

abstract: We present an efficient method to shorten the analytic integration-by-parts (IBP) reduction coefficients of multi-loop Feynman integrals. For our approach, we develop an improved version of Leinartas' multivariate partial fraction algorithm, and provide a modern implementation based on the computer algebra system Singular. Furthermore, We observe that for an integral basis with uniform transcendental (UT) weights, the denominators of IBP reduction coefficients with respect to the UT basis are either symbol letters or polynomials purely in the spacetime dimension $D$. With a UT basis, the partial fraction algorithm is more efficient both with respect to its performance and the size reduction. We show that in complicated examples with existence of a UT basis, the IBP reduction coefficients size can be reduced by a factor of as large as $\sim 100$. We observe that our algorithm also works well for settings without a UT basis.


title: Ray-tracing log-normal simulation for weak gravitational lensing: application to the cross-correlation with galaxies

abstract: We present an algorithm to self-consistently generate mock weak gravitational lensing convergence fields and galaxy distributions in redshift space. We generate three-dimensional cosmic density fields that follow a log-normal distribution, and ray-trace them to produce convergence maps. As we generate the galaxy distribution from the same density fields in a manner consistent with ray-tracing, the galaxy-convergence cross-power spectrum measured from the mock agrees with the theoretical expectation with high precision. We use this simulation to forecast the quality of galaxy-shear cross-correlation measurements from the Subaru Hyper Suprime-Cam (HSC) and Prime Focus Spectrograph (PFS) surveys. We find that the nominal HSC and PFS surveys would detect the cross power spectra with signal-to-noise ratios of 20 and 5 at the lowest ($z = 0.7$) and highest ($z = 2.2$) redshift bins, respectively.


title: Finding Action Tubes with a Sparse-to-Dense Framework

abstract: The task of spatial-temporal action detection has attracted increasing attention among researchers. Existing dominant methods solve this problem by relying on short-term information and dense serial-wise detection on each individual frames or clips. Despite their effectiveness, these methods showed inadequate use of long-term information and are prone to inefficiency. In this paper, we propose for the first time, an efficient framework that generates action tube proposals from video streams with a single forward pass in a sparse-to-dense manner. There are two key characteristics in this framework: (1) Both long-term and short-term sampled information are explicitly utilized in our spatiotemporal network, (2) A new dynamic feature sampling module (DTS) is designed to effectively approximate the tube output while keeping the system tractable. We evaluate the efficacy of our model on the UCF101-24, JHMDB-21 and UCFSports benchmark datasets, achieving promising results that are competitive to state-of-the-art methods. The proposed sparse-to-dense strategy rendered our framework about 7.6 times more efficient than the nearest competitor.


title: Searching for new physics using optically levitated sensors

abstract: We describe a variety of searches for new physics beyond the Standard Model of particle physics which may be enabled in the coming years by the use of optically levitated masses in high vacuum. Such systems are expected to reach force and acceleration sensitivities approaching (and possibly eventually exceeding) the standard quantum limit over the next decade. For new forces or phenomena that couple to mass, high precision sensing using objects with masses in the fg-ng range have significant discovery potential for new physics. Such applications include tests of fundamental force laws, searches for non-neutrality of matter, high-frequency gravitational wave detectors, dark matter searches, and tests of quantum foundations using massive objects.

categories: quant-ph hep-ex hep-ph physics.ins-det physics.optics


title: Measuring and Managing Carbon Risk in Investment Portfolios

abstract: This article studies the impact of carbon risk on stock pricing. To address this, we consider the seminal approach of G\"orgen \textsl{et al.} (2019), who proposed estimating the carbon financial risk of equities by their carbon beta. To achieve this, the primary task is to develop a brown-minus-green (or BMG) risk factor, similar to Fama and French (1992). Secondly, we must estimate the carbon beta using a multi-factor model. While G\"orgen \textsl{et al.} (2019) considered that the carbon beta is constant, we propose a time-varying estimation model to assess the dynamics of the carbon risk. Moreover, we test several specifications of the BMG factor to understand which climate change-related dimensions are priced in by the stock market. In the second part of the article, we focus on the carbon risk management of investment portfolios. First, we analyze how carbon risk impacts the construction of a minimum variance portfolio. As the goal of this portfolio is to reduce unrewarded financial risks of an investment, incorporating the carbon risk into this approach fulfils this objective. Second, we propose a new framework for building enhanced index portfolios with a lower exposure to carbon risk than capitalization-weighted stock indices. Finally, we explore how carbon sensitivities can improve the robustness of factor investing portfolios.


title: The Extended "Sequentially Drilled" Joint Congruence Transformation and its Application in Gaussian Independent Vector Analysis

abstract: Independent Vector Analysis (IVA) has emerged in recent years as an extension of Independent Component Analysis (ICA) into multiple sets of mixtures, where the source signals in each set are independent, but may depend on source signals in the other sets. In a semi-blind IVA (or ICA) framework, information regarding the probability distributions of the sources may be available, giving rise to Maximum Likelihood (ML) separation. In recent work we have shown that under the multivariate Gaussian model, with arbitrary temporal covariance matrices (stationary or non-stationary) of the source signals, ML separation requires the solution of a "Sequentially Drilled" Joint Congruence (SeDJoCo) transformation of a set of matrices, which is reminiscent of (but different from) classical joint diagonalization. In this paper we extend our results to the IVA problem, showing how the ML solution for the Gaussian model (with arbitrary covariance and cross-covariance matrices) takes the form of an extended SeDJoCo problem. We formulate the extended problem, derive a condition for the existence of a solution, and propose two iterative solution algorithms. In addition, we derive the induced Cram\'er-Rao Lower Bound (iCRLB) on the resulting Interference-to-Source Ratios (ISR) matrices, and demonstrate by simulation how the ML separation obtained by solving the extended SeDJoCo problem indeed attains the iCRLB (asymptotically), as opposed to other separation approaches, which cannot exploit prior knowledge regarding the sources' distributions.


title: Arithmetic of Some Sequences Via $2$-determinants

abstract: We extend our investigation of $2$-determinants, which we defined in a previous paper. For a linear homogenous recurrence of the second order, we consider relations between different sequences satisfying the same linear homogeneous recurrence of the second order. After we prove a generalized identity of d'Ocagne, we derive, from a single identity, a number of classical identities (and their generalizations) such as d'Ocagne's, Cassini's, Catalan's and Vajda's. Along the way, the corresponding combinatorial interpretations in terms of restricted words over a finite alphabet are stated for the sequences we investigate.


title: Bottom-up reconstruction of non-singular bounce in F(R) gravity from observational indices

abstract: We apply the bottom-up reconstruction technique in the context of bouncing cosmology in F(R) gravity, where the starting point is a suitable ansatz of observable quantity (like spectral index or tensor to scalar ratio) rather than a priori form of Hubble parameter. In inflationary scenario, the slow roll conditions are assumed to hold true, and thus the observational indices have general expressions in terms of the slow-roll parameters, as for example the tensor to scalar ratio in F(R) inflation can be expressed as $r = 48\epsilon_F^2$ with $\epsilon_F = -\frac{1}{H_F^2}\frac{dH_F}{dt_F}$ and $H_F$, $t_F$ are the Hubble parameter, cosmic time respectively. However, in the bouncing cosmology (say in F(R) gravity theory), the slow-roll conditions are not satisfied, in general, and thus the observable quantities do not have any general expressions that will hold true irrespective of the form of F(R). Thus, in order to apply the bottom-up reconstruction procedure in F(R) bouncing model, we use the conformal correspondence between F(R) and scalar-tensor model where the conformal factor in the present context is chosen in a way such that it leads to an inflationary scenario in the scalar-tensor frame. Due to the reason that the scalar and tensor perturbations remain invariant under conformal transformation, the observable viability of the scalar-tensor inflationary model confirms the viability of the conformally connected F(R) bouncing model. Motivated by these arguments, here we construct a viable non-singular bounce in F(R) gravity directly from the observable indices of the corresponding scalar-tensor inflationary model.


title: 3D-3D Correspondence from Seifert Fibering Operators

abstract: Using recently developed Seifert fibering operators for 3D $\mathcal{N} = 2$ gauge theories, we formulate the necessary ingredients for a state-integral model of the topological quantum field theory dual to a given Seifert manifold under the 3D-3D correspondence, focusing on the case of Seifert homology spheres with positive orbifold Euler characteristic. We further exhibit a set of difference operators that annihilate the wavefunctions of this TQFT on hyperbolic three-manifolds, generalizing similar constructions for lens space partition functions and holomorphic blocks. These properties offer intriguing clues as to the structure of the underlying TQFT.


title: On some $p$-transitive association schemes

abstract: In this paper, for any prime $p$, we propose the notion of a $p$-transitive association scheme. This notion aims to generalize the fact that the regular module of a group algebra of a finite group has a unique trivial submodule to the case of the regular modules of modular adjacency algebras. We completely determine the $p$-transitive quasi-thin association schemes and the $p$-transitive association schemes with thin thin residue by their structure theory properties.


title: Exact Solutions of the 2D Dunkl--Klein--Gordon Equation: The Coulomb Potential and the Klein--Gordon Oscillator

abstract: We introduce the Dunkl--Klein--Gordon (DKG) equation in 2D by changing the standard partial derivatives by the Dunkl derivatives in the standard Klein--Gordon (KG) equation. We show that the generalization with Dunkl derivative of the $z$-component of the angular momentum is what allows the separation of variables of the DKG equation. Then, we compute the energy spectrum and eigenfunctions of the DKG equations for the 2D Coulomb potential and the Klein--Gordon oscillator analytically and from an ${\rm su}(1,1)$ algebraic point of view. Finally, we show that if the parameters of the Dunkl derivative vanish, the obtained results suitably reduce to those reported in the literature for these 2D problems.


title: Large-area femtosecond laser milling of silicon employing trench analysis

abstract: A femtosecond laser is a powerful tool for micromachining of silicon. In this work, large-area laser ablation of crystalline silicon is comprehensively studied using a laser source of pulse width 300 fs at two wavelengths of 343 nm and 1030 nm. We develop a unique approach to gain insight into the laser milling process by means of detailed analysis of trenches. Laser scribed trenches and milled areas are characterized using optical profilometry to extract dimensional and roughness parameters with accuracy and repeatability. In a first step, multiple measures of the trench including the average depth, the volume of recast material, the average longitudinal profile roughness, the inner trench width and the volume removal rate are studied. This allows for delineation of ablation regimes and associated characteristics allowing to determine the impact of fluence and repetition rate on laser milling. In a second step, additional factors of debris formation and material redeposition that come into play during laser milling are further elucidated. These results are utilized for processing large-area (up to few mm2) with milling depths up to 200 {\mu}m to enable the fabrication of cavities with low surface roughness at high removal rates of up to 6.9 {\mu}m3 {\mu}s-1. Finally, laser processing in combination with XeF2 etching is applied on SOI-CMOS technology in the fabrication of radio-frequency (RF) functions standing on suspended membranes. Performance is considerably improved on different functions like RF switch (23 dB improvement in 2nd harmonic), inductors (near doubling of Q-factor) and LNA (noise figure improvement of 0.1 dB) demonstrating the applicability of milling to radio-frequency applications.


title: Free energy sources in current sheets formed in collisionless plasma turbulence

abstract: Collisionless dissipation of macroscopic energy into heat is an unsolved problem of space and astrophysical plasmas, e.g., solar wind and Earth's magnetosheath. The most viable process under consideration is the turbulent-cascade of macroscopic energy to kinetic-scales where collisionless-plasma-processes dissipate the energy. Space observations and numerical simulations show the formation of kinetic scale current sheets in turbulent plasmas. Instabilities in these CS can provide collisionless dissipation and influence the turbulence. Spatial gradients of physical quantities and non-Maxwellian velocity distribution functions provide the free-energy-sources for CS plasma instabilities. To determine the free-energy-sources provided by the spatial gradients of plasma density and electron/ion bulk velocities in CS formed in collisionless turbulent plasmas with an external magnetic field $\mathbf{B}_0$, we carried out two-dimensional PIC-hybrid simulations and interpret the results within the limitations of the simulation model. We found that ion-scale CS in a collisionless turbulent plasma are formed primarily by electron shear flows, i.e., electron bulk velocity inside CS is much larger than ion bulk velocity while the density variations through the CS are relatively small ($<$ 10\%). The electron-bulk-velocity and, thus, the current density inside the sheets are directed mainly parallel to $\mathbf{B}_0$. The shear in the perpendicular electron- and ion-bulk-velocities generates parallel electron- and ion-flow-vorticities. Inside CS, parallel electron-flow-vorticity exceeds the parallel ion-flow-vorticity, changes sign around the CS centers and peaks near the CS edges. An ion temperature anisotropy develops near CS during the CS formation. It has positive correlation with the parallel ion- and electron-flow-vorticities. Theoretical estimates support the simulation results.


title: The Personalized A-Theory of Time and Perspective

abstract: A-theorists and B-theorists debate whether the "Now" is metaphysically distinguished from other time slices. Analogously, one may ask whether the "I" is metaphysically distinguished from other perspectives. Few philosophers would answer the second question in the affirmative. An exception is Caspar Hare, who has devoted two papers and a book to arguing for such a positive answer. In this paper, I argue that those who answer the first question in the affirmative -- A-theorists -- should also answer the second question in the affirmative. This is because key arguments in favor of the A-theory are more effective as arguments in favor of the resulting combined position, and key arguments against the A-theory are ineffective against the combined position.


title: A note on finite determinacy of matrices

abstract: In this note, we give a necessary and sufficient condition for a matrix A in M to be finitely G-determined, where M is the ring of 2 x 2 matrices whose entries are formal power series over an infinite field, and G is a group acting on M by change of coordinates together with multiplication by invertible matrices from both sides.


title: Tight Bound for the Number of Distinct Palindromes in a Tree

abstract: For an undirected tree with $n$ edges labelled by single letters, we consider its substrings, which are labels of the simple paths between pairs of nodes. We prove that there are $O(n^{1.5})$ different palindromic substrings. This solves an open problem of Brlek, Lafreni\`ere, and Proven\c{c}al (DLT 2015), who gave a matching lower-bound construction. Hence, we settle the tight bound of $\Theta(n^{1.5})$ for the maximum palindromic complexity of trees. For standard strings, i.e., for paths, the palindromic complexity is $n+1$. We also propose $O(n^{1.5} \log{n})$-time algorithm for reporting all distinct palindromes in an undirected tree with $n$ edges.


title: Multiway $p$-spectral graph cuts on Grassmann manifolds

abstract: Nonlinear reformulations of the spectral clustering method have gained a lot of recent attention due to their increased numerical benefits and their solid mathematical background. We present a novel direct multiway spectral clustering algorithm in the $p$-norm, for $p \in (1, 2]$. The problem of computing multiple eigenvectors of the graph $p$-Laplacian, a nonlinear generalization of the standard graph Laplacian, is recasted as an unconstrained minimization problem on a Grassmann manifold. The value of $p$ is reduced in a pseudocontinuous manner, promoting sparser solution vectors that correspond to optimal graph cuts as $p$ approaches one. Monitoring the monotonic decrease of the balanced graph cuts guarantees that we obtain the best available solution from the $p$-levels considered. We demonstrate the effectiveness and accuracy of our algorithm in various artificial test-cases. Our numerical examples and comparative results with various state-of-the-art clustering methods indicate that the proposed method obtains high quality clusters both in terms of balanced graph cut metrics and in terms of the accuracy of the labelling assignment. Furthermore, we conduct studies for the classification of facial images and handwritten characters to demonstrate the applicability in real-world datasets.


title: The human behavioural immune system is a product of cultural evolution

abstract: Disgust is a basic emotion that serves to avoid contaminants, and is central to the behavioural immune system. While disgust-motivated avoidance occurs in bonobos and chimpanzees, humans show uniquely high levels of contamination sensitivity. Current theory postulates that human disgust is primarily a genetic adaptation, side-lining social transmission and denying parental modelling. Here, we test whether this strong view is warranted by simulating 100000 years of cultural and biological evolution in various scenarios. We modelled disgust as a trait that governed the extent to which individuals forewent potentially contaminated nutrition. This lowered their risk of gastrointestinal illness, at the cost of increasing starvation risk and birth-interval. As environmental levels and costs of contamination rose, so did relative fitness of high-disgust individuals. Evolutionary shifts in disgust as polygenic trait occurred as a consequence of germ-cell mutations, but were most prominent in populations with high initial genetic variance. Crucially, cultural transmission between generations operated at a higher rate, even if parental modelling was eliminated. This study serves not only as evidence of cultural evolution shaping the behavioural immune system, but also as an illustration of emerging theories that paint affective and cognitive mechanisms as socially transmitted rather than biologically determined functions.


title: Reinforcement Learning Based Penetration Testing of a Microgrid Control Algorithm

abstract: Microgrids (MGs) are small-scale power systems which interconnect distributed energy resources and loads within clearly defined regions. However, the digital infrastructure used in an MG to relay sensory information and perform control commands can potentially be compromised due to a cyberattack from a capable adversary. An MG operator is interested in knowing the inherent vulnerabilities in their system and should regularly perform Penetration Testing (PT) activities to prepare for such an event. PT generally involves looking for defensive coverage blindspots in software and hardware infrastructure, however the logic in control algorithms which act upon sensory information should also be considered in PT activities. This paper demonstrates a case study of PT for an MG control algorithm by using Reinforcement Learning (RL) to uncover malicious input which compromises the effectiveness of the controller. Through trial-and-error episodic interactions with a simulated MG, we train an RL agent to find malicious input which reduces the effectiveness of the MG controller.


title: Mixture of Speaker-type PLDAs for Children's Speech Diarization

abstract: In diarization, the PLDA is typically used to model an inference structure which assumes the variation in speech segments be induced by various speakers. The speaker variation is then learned from the training data. However, human perception can differentiate speakers by age, gender, among other characteristics. In this paper, we investigate a speaker-type informed model that explicitly captures the known variation of speakers. We explore a mixture of three PLDA models, where each model represents an adult female, male, or child category. The weighting of each model is decided by the prior probability of its respective class, which we study. The evaluation is performed on a subset of the BabyTrain corpus. We examine the expected performance gain using the oracle speaker type labels, which yields an 11.7% DER reduction. We introduce a novel baby vocalization augmentation technique and then compare the mixture model to the single model. Our experimental result shows an effective 0.9% DER reduction obtained by adding vocalizations. We discover empirically that a balanced dataset is important to train the mixture PLDA model, which outperforms the single PLDA by 1.3% using the same training data and achieving a 35.8% DER. The same setup improves over a standard baseline by 2.8% DER.


title: Native Skywatchers and the Ojibwe Giizhig Anung Masinaaigan-Ojibwe Sky Star Map

abstract: The Native Skywatchers programming addresses the loss of Ojibwe and Dakota star knowledge. The goal is to help preserve indigenous astronomy and pass it on to future generations. The star knowledge will serve as a cornerstone inspiration for native communities and at the same time be influential to Native students interested in science and math. Presented here is the Ojibwe Giizhig Anung Masinaaigan-Ojibwe Sky Star Map created by Annette Lee and collaborators as part of the Native Skywatchers Middle School Teachers Workshop held June 2012.


title: Astraeus III: The environment and physical properties of reionization sources

abstract: In this work, we use the {\sc astraeus} (seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework which couples galaxy formation and reionization in the first billion years. Exploring a number of models for reionization feedback and the escape fraction of ionizing radiation from the galactic environment ($f_\mathrm{esc}$), we quantify how the contribution of star-forming galaxies {(with halo masses $M_h>10^{8.2}$M$_\odot$)} to reionization depends on the radiative feedback model, $f_\mathrm{esc}$, and the environmental over-density. Our key findings are: (i) for constant $f_\mathrm{esc}$ models, intermediate-mass galaxies (with halo masses of $M_h\simeq10^{9-11}$M$_\odot$ and absolute UV magnitudes of $M_{UV} \sim -15$ to $-20$) in intermediate-density regions drive reionization; (ii) scenarios where $f_\mathrm{esc}$ increases with decreasing halo mass shift the galaxy population driving reionization to lower-mass galaxies ($M_h\lesssim10^{9.5}$M$_\odot$) with lower luminosities ($M_{UV} \gtrsim-16$) and over-densities; (iii) reionization imprints its topology on the ionizing emissivity of low-mass galaxies ($M_h\lesssim10^{9}$M$_\odot$) through radiative feedback. Low-mass galaxies experience a stronger suppression of star formation by radiative feedback and show lower ionizing emissivities in over-dense regions; (iv) a change in $f_\mathrm{esc}$ with galaxy properties has the largest impact on the sources of reionization and their detectability, with the radiative feedback strength and environmental over-density playing a sub-dominant role; (v) JWST-surveys (with a limiting magnitude of $M_{UV} = -16$) will be able to detect the galaxies providing $\sim 60-70\%$ ($\sim 10\%$) of reionization photons at $z=7$ for constant $f_\mathrm{esc}$ models (scenarios where $f_\mathrm{esc}$ increases with decreasing halo mass).


title: Frame-Based Filter-Function Formalism for Quantum Characterization and Control

abstract: We introduce a theoretical framework for resource-efficient characterization and control of non-Markovian open quantum systems, which naturally allows for the integration of given, experimentally motivated, control capabilities and constraints. This is achieved by developing a transfer filter-function formalism based on the general notion of a frame and by appropriately tying the choice of frame to the available control. While recovering the standard frequency-based filter-function formalism as a special instance, this control-adapted generalization affords intrinsic flexibility and, crucially, it permits an efficient representation of the relevant control matrix elements and dynamical integrals if an appropriate finite-size frame condition is obeyed. Our frame-based formulation overcomes important limitations of existing approaches. In particular, we show how to implement quantum noise spectroscopy in the presence of nonstationary noise sources, and how to effectively achieve control-driven model reduction for noise-optimized prediction and quantum gate design.


title: Singular function emerging from one-dimensional elementary cellular automaton Rule 150

abstract: In this paper, we give a singular function on a unit interval derived from the dynamic of the one-dimensional elementary cellular automaton Rule 150. We describe properties of the resulting function, that is strictly increasing, uniformly continuous, and differentiable almost everywhere, and we show that it is not differentiable at dyadic rational points. We also give functional equations that the function satisfies, and show that the function is the only solution of the functional ones.


title: Covering numbers of commutative rings

abstract: A cover of a unital, associative (not necessarily commutative) ring $R$ is a collection of proper subrings of $R$ whose set-theoretic union equals $R$. If such a cover exists, then the covering number $\sigma(R)$ of $R$ is the cardinality of a minimal cover, and a ring $R$ is called $\sigma$-elementary if $\sigma(R) < \sigma(R/I)$ for every nonzero two-sided ideal $I$ of $R$. In this paper, we show that if $R$ has a finite covering number, then the calculation of $\sigma(R)$ can be reduced to the case where $R$ is a finite ring of characteristic $p$ and the Jacobson radical $J$ of $R$ has nilpotency 2. Our main result is that if $R$ has a finite covering number and $R/J$ is commutative (even if $R$ itself is not), then either $\sigma(R)=\sigma(R/J)$, or $\sigma(R)=p^d+1$ for some $d \geqslant 1$. As a byproduct, we classify all commutative $\sigma$-elementary rings with a finite covering number and characterize the integers that occur as the covering number of a commutative ring.


title: Magnetic field turbulence in the solar wind at sub-ion scales: in situ observations and numerical simulations

abstract: We investigate the transition of the solar wind turbulent cascade from MHD to sub-ion range by means of a detail comparison between in situ observations and hybrid numerical simulations. In particular we focus on the properties of the magnetic field and its component anisotropy in Cluster measurements and hybrid 2D simulations. First, we address the angular distribution of wave-vectors in the kinetic range between ion and electron scales by studying the variance anisotropy of the magnetic field components. When taking into account the single-direction sampling performed by spacecraft in the solar wind, the main properties of the fluctuations observed in situ are also recovered in our numerical description. This result confirms that solar wind turbulence in the sub-ion range is characterized by a quasi-2D gyrotropic distribution of k-vectors around the mean field. We then consider the magnetic compressibility associated with the turbulent cascade and its evolution from large-MHD to sub-ion scales. The ratio of field-aligned to perpendicular fluctuations, typically low in the MHD inertial range, increases significantly when crossing ion scales and its value in the sub-ion range is a function of the total plasma beta only, as expected from theoretical predictions, with higher magnetic compressibility for higher beta. Moreover, we observe that this increase has a gradual trend from low to high beta values in the in situ data; this behaviour is well captured by the numerical simulations. The level of magnetic field compressibility that is observed in situ and in the simulations is in fairly good agreement with theoretical predictions, especially at high beta, suggesting that in the kinetic range explored the turbulence is supported by low-frequency and highly-oblique fluctuations in pressure balance, like kinetic Alfv\'en waves or other slowly evolving coherent structures.


title: Shell model structure of proxy-SU(3) pairs of orbitals

abstract: The Nilsson orbitals used in the substitutions occurring in the proxy-SU(3) scheme, which are the orbitals bearing the maximum value of total angular momentum in each shell, have an extremely simple structure in the shell model basis |N l j Omega>, with each Nilsson orbital corresponding to a single shell model eigenvector. This simple structure is valid at all deformations for these orbitals, while in other orbitals it is valid only at small deformations. Nilsson 0[110] pairs are found to correspond to |1 1 1 0> pairs in the spherical shell model basis, paving the way for using the proxy-SU(3) approximation within the shell model.


title: Human-in-the-Loop Methods for Data-Driven and Reinforcement Learning Systems

abstract: Recent successes combine reinforcement learning algorithms and deep neural networks, despite reinforcement learning not being widely applied to robotics and real world scenarios. This can be attributed to the fact that current state-of-the-art, end-to-end reinforcement learning approaches still require thousands or millions of data samples to converge to a satisfactory policy and are subject to catastrophic failures during training. Conversely, in real world scenarios and after just a few data samples, humans are able to either provide demonstrations of the task, intervene to prevent catastrophic actions, or simply evaluate if the policy is performing correctly. This research investigates how to integrate these human interaction modalities to the reinforcement learning loop, increasing sample efficiency and enabling real-time reinforcement learning in robotics and real world scenarios. This novel theoretical foundation is called Cycle-of-Learning, a reference to how different human interaction modalities, namely, task demonstration, intervention, and evaluation, are cycled and combined to reinforcement learning algorithms. Results presented in this work show that the reward signal that is learned based upon human interaction accelerates the rate of learning of reinforcement learning algorithms and that learning from a combination of human demonstrations and interventions is faster and more sample efficient when compared to traditional supervised learning algorithms. Finally, Cycle-of-Learning develops an effective transition between policies learned using human demonstrations and interventions to reinforcement learning. The theoretical foundation developed by this research opens new research paths to human-agent teaming scenarios where autonomous agents are able to learn from human teammates and adapt to mission performance metrics in real-time and in real world scenarios.


title: Improved Lite Audio-Visual Speech Enhancement

abstract: Numerous studies have investigated the effectiveness of audio-visual multimodal learning for speech enhancement (AVSE) tasks, seeking a solution that uses visual data as auxiliary and complementary input to reduce the noise of noisy speech signals. Recently, we proposed a lite audio-visual speech enhancement (LAVSE) algorithm for a car-driving scenario. Compared to conventional AVSE systems, LAVSE requires less online computation and to some extent solves the user privacy problem on facial data. In this study, we extend LAVSE to improve its ability to address three practical issues often encountered in implementing AVSE systems, namely, the additional cost of processing visual data, audio-visual asynchronization, and low-quality visual data. The proposed system is termed improved LAVSE (iLAVSE), which uses a convolutional recurrent neural network architecture as the core AVSE model. We evaluate iLAVSE on the Taiwan Mandarin speech with video dataset. Experimental results confirm that compared to conventional AVSE systems, iLAVSE can effectively overcome the aforementioned three practical issues and can improve enhancement performance. The results also confirm that iLAVSE is suitable for real-world scenarios, where high-quality audio-visual sensors may not always be available.


title: Deep Reinforcement Learning for Contact-Rich Skills Using Compliant Movement Primitives

abstract: In recent years, industrial robots have been installed in various industries to handle advanced manufacturing and high precision tasks. However, further integration of industrial robots is hampered by their limited flexibility, adaptability and decision making skills compared to human operators. Assembly tasks are especially challenging for robots since they are contact-rich and sensitive to even small uncertainties. While reinforcement learning (RL) offers a promising framework to learn contact-rich control policies from scratch, its applicability to high-dimensional continuous state-action spaces remains rather limited due to high brittleness and sample complexity. To address those issues, we propose different pruning methods that facilitate convergence and generalization. In particular, we divide the task into free and contact-rich sub-tasks, perform the control in Cartesian rather than joint space, and parameterize the control policy. Those pruning methods are naturally implemented within the framework of dynamic movement primitives (DMP). To handle contact-rich tasks, we extend the DMP framework by introducing a coupling term that acts like the human wrist and provides active compliance under contact with the environment. We demonstrate that the proposed method can learn insertion skills that are invariant to space, size, shape, and closely related scenarios, while handling large uncertainties. Finally we demonstrate that the learned policy can be easily transferred from simulations to real world and achieve similar performance on UR5e robot.


title: Oriented cycles in digraphs of large outdegree

abstract: In 1985, Mader conjectured that for every acyclic digraph $F$ there exists $K=K(F)$ such that every digraph $D$ with minimum out-degree at least $K$ contains a subdivision of $F$. This conjecture remains widely open, even for digraphs $F$ on five vertices. Recently, Aboulker, Cohen, Havet, Lochet, Moura and Thomass\'{e} studied special cases of Mader's problem and made the following conjecture: for every $\ell \geq 2$ there exists $K = K(\ell)$ such that every digraph $D$ with minimum out-degree at least $K$ contains a subdivision of every orientation of a cycle of length $\ell$. We prove this conjecture and answer further open questions raised by Aboulker et al.


title: SOLAR: Sparse Orthogonal Learned and Random Embeddings

abstract: Dense embedding models are commonly deployed in commercial search engines, wherein all the document vectors are pre-computed, and near-neighbor search (NNS) is performed with the query vector to find relevant documents. However, the bottleneck of indexing a large number of dense vectors and performing an NNS hurts the query time and accuracy of these models. In this paper, we argue that high-dimensional and ultra-sparse embedding is a significantly superior alternative to dense low-dimensional embedding for both query efficiency and accuracy. Extreme sparsity eliminates the need for NNS by replacing them with simple lookups, while its high dimensionality ensures that the embeddings are informative even when sparse. However, learning extremely high dimensional embeddings leads to blow up in the model size. To make the training feasible, we propose a partitioning algorithm that learns such high dimensional embeddings across multiple GPUs without any communication. This is facilitated by our novel asymmetric mixture of Sparse, Orthogonal, Learned and Random (SOLAR) Embeddings. The label vectors are random, sparse, and near-orthogonal by design, while the query vectors are learned and sparse. We theoretically prove that our way of one-sided learning is equivalent to learning both query and label embeddings. With these unique properties, we can successfully train 500K dimensional SOLAR embeddings for the tasks of searching through 1.6M books and multi-label classification on the three largest public datasets. We achieve superior precision and recall compared to the respective state-of-the-art baselines for each of the tasks with up to 10 times faster speed.


title: Norm inequalities related to operator monotone functions

abstract: Let $A$ be a positive definite operator on a Hilbert space $H$, and $|||.|||$ be a unitarily invariant norm on $B(H)$. We show that if $f$ is an operator monotone function on $(0,\infty)$ and $n\in \mathbb{N}$, then $|||D^n f(A)|||\leq\|f^{(n)}(A)\|$ and $\|f^{(n)}(\cdot)\|$ is a quasi-convex function on the set of all positive definite operators in $B(H)$. We establish some estimates of the right hand side of some Hermite-Hadamard type inequalities in which differentiable functions are involved, and norms of the maps induced by them on the set of self adjoint operators are convex, quasi-convex or $s$-convex. As applications, we obtain some of bounds for $|||f(B)-f(A)|||$ in term of $|||B-A|||$. For instance, Let $f,g$ be two operator monotone functions on $(0,\infty)$. Then, for every unitarily invariant norm $|||.|||$ and every positive definite operators $A,B$, \begin{align*} &\left|\left|\left|f(A)g(A)-f(B)g(B)\right|\right|\right|\notag\\ &\leq|||B-A|||\Big[\max\left\{\|f'(A)\|,\|f'(B)\|\right\}\times\max\left\{\|g(A)\|,\|g(B)\|\right\}\notag\\ &+\max\left\{\|f(A)\|,\|f(B)\|\right\}\times \max\left\{\|g'(A)\|,\|g'(B)\|\right\}\Big]. \end{align*}


title: A Compact Deep Architecture for Real-time Saliency Prediction

abstract: Saliency computation models aim to imitate the attention mechanism in the human visual system. The application of deep neural networks for saliency prediction has led to a drastic improvement over the last few years. However, deep models have a high number of parameters which makes them less suitable for real-time applications. Here we propose a compact yet fast model for real-time saliency prediction. Our proposed model consists of a modified U-net architecture, a novel fully connected layer, and central difference convolutional layers. The modified U-Net architecture promotes compactness and efficiency. The novel fully-connected layer facilitates the implicit capturing of the location-dependent information. Using the central difference convolutional layers at different scales enables capturing more robust and biologically motivated features. We compare our model with state of the art saliency models using traditional saliency scores as well as our newly devised scheme. Experimental results over four challenging saliency benchmark datasets demonstrate the effectiveness of our approach in striking a balance between accuracy and speed. Our model can be run in real-time which makes it appealing for edge devices and video processing.


title: Stochastically Realized Observables for Excitonic Molecular Aggregates

abstract: We show that a stochastic approach enables calculations of the optical properties of large 2-dimensional and nanotubular excitonic molecular aggregates. Previous studies of such systems relied on numerically diagonalizing the dense and disordered Frenkel Hamiltonian, which scales approximately as $\mathcal{O}(N^3)$ for $N$ dye molecules. Our approach scales much more efficiently as $\mathcal{O}(N\log(N))$, enabling quick study of systems with a million of coupled molecules on the micron size scale. We calculate several important experimental observable including the optical absorption spectrum and density of states, and develop a stochastic formalism for the participation ratio. Quantitative agreement with traditional matrix diagonalization methods is demonstrated for both small- and intermediate-size systems. The stochastic methodology enables the study of the effects of spatial-correlation in site energies on the optical signatures of large 2D aggregates. Our results demonstrate that stochastic methods present a path forward for screening structural parameters and validating experiments and theoretical predictions in large excitonic aggregates.


title: An evolutionary perspective on the design of neuromorphic shape filters

abstract: A substantial amount of time and energy has been invested to develop machine vision using connectionist (neural network) principles. Most of that work has been inspired by theories advanced by neuroscientists and behaviorists for how cortical systems store stimulus information. Those theories call for information flow through connections among several neuron populations, with the initial connections being random (or at least non-functional). Then the strength or location of connections are modified through training trials to achieve an effective output, such as the ability to identify an object. Those theories ignored the fact that animals that have no cortex, e.g., fish, can demonstrate visual skills that outpace the best neural network models. Neural circuits that allow for immediate effective vision and quick learning have been preprogrammed by hundreds of millions of years of evolution and the visual skills are available shortly after hatching. Cortical systems may be providing advanced image processing, but most likely are using design principles that had been proven effective in simpler systems. The present article provides a brief overview of retinal and cortical mechanisms for registering shape information, with the hope that it might contribute to the design of shape-encoding circuits that more closely match the mechanisms of biological vision.


title: A continuous-time asset market game with short-lived assets

abstract: We consider a continuous-time game-theoretic model of an investment market with short-lived assets and endogenous asset prices. The first goal of the paper is to formulate a stochastic equation which determines wealth processes of investors and to provide conditions for the existence of its solution. The second goal is to show that there exists a strategy such that the logarithm of the relative wealth of an investor who uses it is a submartingale regardless of the strategies of the other investors, and the relative wealth of any other essentially different strategy vanishes asymptotically. This strategy can be considered as an optimal growth portfolio in the model.


title: Direct measurement of the $7s\ ^2S_{1/2}\rightarrow$ $7p\ ^2P_{3/2}$ transition frequency in $^{226}$Ra$^{+}$

abstract: We report a direct measurement of the $7s\ ^2S_{1/2}\rightarrow$ $7p\ ^2P_{3/2}$ electric dipole transition frequency in $^{226}$Ra$^{+}$. With a single laser-cooled radium ion we determine the transition frequency to be 785722.11(3) GHz by directly driving the transition with frequency-doubled light and measuring the frequency of the undoubled light with an iodine reference. This measurement addresses a discrepancy of five combined standard deviations between previously reported values.


title: On a rank-unimodality conjecture of Morier-Genoud and Ovsienko

abstract: Let alpha = (a,b,...) be a composition. Consider the associated poset F(alpha), called a fence, whose covering relations are x_1 < x_2 < ... < x_{a+1} > x_{a+2} > ... > x_{a+b+1} < x_{a+b+2} < ... . We study the associated distributive lattice L(alpha) consisting of all lower order ideals of F(alpha). These lattices are important in the theory of cluster algebras and their rank generating functions can be used to define q-analogues of rational numbers. In particular, we make progress on a recent conjecture of Morier-Genoud and Ovsienko that L(alpha) is rank unimodal. We show that if one of the parts of alpha is greater than the sum of the others, then the conjecture is true. We conjecture that L(alpha) enjoys the stronger properties of having a nested chain decomposition and having a rank sequence which is either top or bottom interlacing, the latter being a recently defined property of sequences. We verify that these properties hold for compositions with at most three parts and for what we call d-divided posets, generalizing work of Claussen and simplifying a construction of Gansner.


title: Finite and divergent parts of the self-force of a point charge from its spherically averaged self-field

abstract: The electromagnetic self-force of a point charge moving arbitrarily on a rectilinear trajectory is calculated by averaging its retarded electric self-field over a sphere of infinitesimal radius centered on the charge's present position. The finite part of the self-force obtained is the well-established relativistic radiation reaction, while its divergent part implies the pre-relativistic longitudinal electromagnetic mass of Abraham.


title: An all-optical intrinsic atomic gradiometer with sub-20 fT/cm/$\sqrt{\rm Hz}$ sensitivity in a 22 $\mu$T earth-scale magnetic field

abstract: In this work we demonstrate a high sensitivity atomic gradiometer capable of operation in earth-field level environments. We apply a light-pulse sequence at four times the Larmor frequency to achieve gradiometer sensitivity <20 fT/cm/$\sqrt{\rm Hz}$ at the finite field strength of 22 $\mu$T. The experimental timing sequence can be tuned to the field magnitude of interest. Our one dimensional all-optical gradiometer performs a differential measurement between two regions of a single vapor cell on a 4 cm baseline. Our results pave the way for extensions to operation in higher dimensions, vector sensitivity, and more advanced gradiometers.


title: An Objective for Hierarchical Clustering in Euclidean Space and its Connection to Bisecting K-means

abstract: This paper explores hierarchical clustering in the case where pairs of points have dissimilarity scores (e.g. distances) as a part of the input. The recently introduced objective for points with dissimilarity scores results in every tree being a 1/2 approximation if the distances form a metric. This shows the objective does not make a significant distinction between a good and poor hierarchical clustering in metric spaces. Motivated by this, the paper develops a new global objective for hierarchical clustering in Euclidean space. The objective captures the criterion that has motivated the use of divisive clustering algorithms: that when a split happens, points in the same cluster should be more similar than points in different clusters. Moreover, this objective gives reasonable results on ground-truth inputs for hierarchical clustering. The paper builds a theoretical connection between this objective and the bisecting k-means algorithm. This paper proves that the optimal 2-means solution results in a constant approximation for the objective. This is the first paper to show the bisecting k-means algorithm optimizes a natural global objective over the entire tree.


title: Discrete Curvature and Torsion from Cross-Ratios

abstract: Motivated by a M\"obius invariant subdivision scheme for polygons, we study a curvature notion for discrete curves where the cross-ratio plays an important role in all our key definitions. Using a particular M\"obius invariant point-insertion-rule, comparable to the classical four-point-scheme, we construct circles along discrete curves. Asymptotic analysis shows that these circles defined on a sampled curve converge to the smooth curvature circles as the sampling density increases. We express our discrete torsion for space curves, which is not a M\"obius invariant notion, using the cross-ratio and show its asymptotic behavior in analogy to the curvature.


title: On Solutions to the Nonlocal $\overline{\partial}$ Problem and (2+1) Dimensional Completely Integrable Systems

abstract: In this short note we discuss a new formula for solving the nonlocal $\overline{\partial}$-problem, and discuss application to the Manakov--Zakharov dressing method. We then explicitly apply this formula to solving the complex (2+1)D Kadomtsev--Petviashvili equation and complex (2+1)D completely integrable generalization of the (2+1)D Kaup--Broer (or Kaup--Boussinesq) system. We will also discuss how real (1+1)D solutions are expressed using this formalism. It is simple to express the formalism for finite gap primitive solutions from [10], [8] using the formalism of this note. We also discuss recent results on the infinite soliton limit for the (1+1)D Korteweg--de Vries equation and the (2+1)D Kaup--Broer system. In an appendix, the classical solutions to the 3D Laplace equation (2+1)D d'Alembert wave equation by Whittaker are described. This appendix is included to elucidate an analogy between the dressing method and the Whittaker solutions.


title: Utilization of 3D segmentation for measurement of pediatric brain tumor outcomes after treatment: review of available free tools, step-by-step instructions, and applications to clinical practice

abstract: Volumetric measurements are known to provide more information when it comes to segmenting tumors, in comparison to one- and two-dimensional measurements, and thus can lead to better informed therapy. In this work, we review the free and easily accessible computer platforms available for conducting these 3D measurements, such as Horos and 3D Slicer and compare the segmentations to commercial Visage software. We compare the time for 3D segmentation of tumors and demonstrate how to use a novel plugin that we developed in 3D slicer for the efficient and accurate segmentation of the cystic component of a tumor.


title: Convex Optimization of Launch Vehicle Ascent Trajectory with Heat-Flux and Splash-Down Constraints

abstract: This paper presents a convex programming approach to the optimization of a multistage launch vehicle ascent trajectory, from the liftoff to the payload injection into the target orbit, taking into account multiple nonconvex constraints, such as the maximum heat flux after fairing jettisoning and the splash-down of the burned-out stages. Lossless and successive convexification are employed to convert the problem into a sequence of convex subproblems. Virtual controls and buffer zones are included to ensure the recursive feasibility of the process and a state-of-the-art method for updating the reference solution is implemented to filter out undesired phenomena that may hinder convergence. A hp pseudospectral discretization scheme is used to accurately capture the complex ascent and return dynamics with a limited computational effort. The convergence properties, computational efficiency, and robustness of the algorithm are discussed on the basis of numerical results. The ascent of the VEGA launch vehicle toward a polar orbit is used as case study to discuss the interaction between the heat flux and splash-down constraints. Finally, a sensitivity analysis of the launch vehicle carrying capacity to different splash-down locations is presented.


title: Employing Simulation to Facilitate the Design of Dynamic Code Generators

abstract: Dynamic Translation (DT) is a sophisticated technique that allows the implementation of high-performance emulators and high-level-language virtual machines. In this technique, the guest code is compiled dynamically at runtime. Consequently, achieving good performance depends on several design decisions, including the shape of the regions of code being translated. Researchers and engineers explore these decisions to bring the best performance possible. However, a real DT engine is a very sophisticated piece of software, and modifying one is a hard and demanding task. Hence, we propose using simulation to evaluate the impact of design decisions on dynamic translators and present RAIn, an open-source DT simulator that facilitates the test of DT's design decisions, such as Region Formation Techniques (RFTs). RAIn outputs several statistics that support the analysis of how design decisions may affect the behavior and the performance of a real DT. We validated RAIn running a set of experiments with six well known RFTs (NET, MRET2, LEI, NETPlus, NET-R, and NETPlus-e-r) and showed that it can reproduce well-known results from the literature without the effort of implementing them on a real and complex dynamic translator engine.


title: Minimum Dielectric-Resonator Mode Volumes

abstract: We show that global lower bounds to the mode volume of a dielectric resonator can be computed via Lagrangian duality. State-of-the-art designs rely on sharp tips, but such structures appear to be highly sub-optimal at nanometer-scale feature sizes, and we demonstrate that computational inverse design offers orders-of-magnitude possible improvements. Our bound can be applied for geometries that are simultaneously resonant at multiple frequencies, for high-efficiency nonlinear-optics applications, and we identify the unavoidable penalties that must accompany such multiresonant structures.


title: V1460 Her: A fast spinning white dwarf accreting from an evolved donor star

abstract: We present time-resolved optical and ultraviolet spectroscopy and photometry of V1460~Her, an eclipsing cataclysmic variable with a 4.99\,h orbital period and an overluminous K5-type donor star. The optical spectra show emission lines from an accretion disc along with absorption lines from the donor. We use these to measure radial velocities, which, together with constraints upon the orbital inclination from photometry, imply masses of $M_1=0.869\pm0.006\,\mathrm{M}_\odot$ and $M_2=0.295\pm0.004\,\mathrm{M}_\odot$ for the white dwarf and the donor. The radius of the donor, $R_2=0.43\pm0.002\,\mathrm{R}_\odot$, is $\approx 50$ per cent larger than expected given its mass, while its spectral type is much earlier than the M3.5 type that would be expected from a main sequence star with a similar mass. HST spectra show strong $\mathrm{N{\small V}}$ 1240 A emission but no $\mathrm{C{\small IV}}$ 1550 A emission, evidence for CNO-processed material. The donor is therefore a bloated, over-luminous remnant of a thermal-timescale stage of high mass transfer and has yet to re-establish thermal equilibrium. Remarkably, the HST ultraviolet data also show a strong 30 per cent peak-to-peak, $38.9\,$s pulsation that we explain as being due to the spin of the white dwarf, potentially putting V1460 Her in a similar category to the propeller system AE Aqr in terms of its spin frequency and evolutionary path. AE Aqr also features a post-thermal timescale mass donor, and V1460 Her may therefore be its weak magnetic field analogue since the accretion disc is still present, with the white dwarf spin-up a result of a recent high accretion rate.


title: Semitoric degenerations of Hibi varieties and flag varieties

abstract: We construct a family of flat semitoric degenerations for the Hibi variety of every finite distributive lattice. The irreducible components of each degeneration are the toric varieties associated with polytopes forming a regular subdivision of the order polytope of the underlying poset. These components are themselves Hibi varieties. For each degeneration in our family we also define the corresponding weight polytope and embed the degeneration into the associated toric variety as the union of orbit closures given by a set of faces. Every such weight polytope projects onto the order polytope with the chosen faces projecting into the parts of the corresponding regular subdivision. We apply these constructions to obtain a family of flat semitoric degenerations for every type A Grassmannian and complete flag variety.


title: The prestige and status of research fields within mathematics

abstract: While the ``hierarchy of science'' has been widely analysed, there is no corresponding study of the status of subfields within a given scientific field. We use bibliometric data to show that subfields of mathematics have a different ``standing'' within the mathematics community. Highly ranked departments tend to specialize in some subfields more than in others, and the same subfields are also over-represented in the most selective mathematics journals or among recipients of top prizes. Moreover this status of subfields evolves markedly over the period of observation (1984--2016), with some subfields gaining and others losing in standing. The status of subfields is related to different publishing habits, but some of those differences are opposite to those observed when considering the hierarchy of scientific fields. We examine possible explanations for the ``status'' of different subfields. Some natural explanations -- availability of funding, importance of applications -- do not appear to function, suggesting that factors internal to the discipline are at work. We propose a different type of explanation, based on a notion of ``focus'' of a subfield, that might or might not be specific to mathematics.


title: Floating-Point Multiplication Using Neuromorphic Computing

abstract: Neuromorphic computing describes the use of VLSI systems to mimic neuro-biological architectures and is also looked at as a promising alternative to the traditional von Neumann architecture. Any new computing architecture would need a system that can perform floating-point arithmetic. In this paper, we describe a neuromorphic system that performs IEEE 754-compliant floating-point multiplication. The complex process of multiplication is divided into smaller sub-tasks performed by components Exponent Adder, Bias Subtractor, Mantissa Multiplier and Sign OF/UF. We study the effect of the number of neurons per bit on accuracy and bit error rate, and estimate the optimal number of neurons needed for each component.


title: On partial geometries arising from maximal arcs

abstract: The subject of this paper are partial geometries $pg(s,t,\alpha)$ with parameters $s=d(d'-1), \ t=d'(d-1), \ \alpha=(d-1)(d'-1)$, $d, d' \ge 2$. In all known examples, $q=dd'$ is a power of 2 and the partial geometry arises from a maximal arc of degree $d$ or $d'$ in a projective plane of order $q$ via a known construction due to Thas \cite{Thas73} and Wallis \cite{W}, with a single known exception of a partial geometry $pg(4,6,3)$ found by Mathon \cite{Math} that is not associated with a maximal arc in the projective plane of order 8. A parallel class of lines is a set of pairwise disjoint lines that covers the point set. Two parallel classes are called orthogonal if they share exactly one line. An upper bound on the maximum number of pairwise orthogonal parallel classes in a partial geometry $G$ with parameters $pg(d(d'-1),d'(d-1),(d-1)(d'-1))$ is proved, and it is shown that a necessary and sufficient condition for $G$ to arise from a maximal arc of degree $d$ or $d'$ in a projective plane of order $q=dd'$ is that both $G$ and its dual geometry contain sets of pairwise orthogonal parallel classes that meet the upper bound. An alternative construction of Mathon's partial geometry is presented, and the new necessary condition is used to demonstrate why this partial geometry is not associated with any maximal arc in the projective plane of order 8. The partial geometries associated with all known maximal arcs in projective planes of order 16 are classified up to isomorphism, and their parallel classes of lines and the 2-rank of their incidence matrices are computed. Based on these results, some open problems and conjectures are formulated.


title: Hochschild lattices and shuffle lattices

abstract: In his study of a Hochschild complex arising in connection with the free loop fibration, S. Saneblidze defined the freehedron, a certain polytope constructed via a truncation process from the hypercube. It was recently conjectured by F. Chapoton and proven by C. Combe that a certain orientation of the $1$-skeleton of the freehedron carries a lattice structure. The resulting lattice was dubbed the Hochschild lattice and is congruence uniform and extremal. These properties allow for the definition of three associated structures: the Galois graph, the canonical join complex and the core label order. In this article, we study and characterize these structures. We exhibit an isomorphism from the core label order of the Hochschild lattice to a particular shuffle lattice of C. Greene. We also uncover an enumerative connection between the core label order of the Hochschild lattice, a certain order extension of its poset of irreducibles and the freehedron. These connections nicely parallel the situation surrounding the better-known Tamari lattices, noncrossing partition lattices and associahedra.


title: Lattice dynamics of topological Dirac semimetal LaAgSb$_{2}$ with charge density wave ordering

abstract: LaAgSb$_{2}$ is a rare material, which offers the opportunity to investigate the complex interplay between charge density wave (CDW) ordering and topology protected electronic band structure. As both of these phenomena are governed by the structural symmetries, a comprehensive study of the lattice dynamics is highly desirable. In this report, we present the results of temperature and pressure dependent Raman spectroscopy and x-ray diffraction in single crystalline LaAgSb$_{2}$. Our results confirm that Raman spectroscopy is a highly sensitive tool to probe CDW ordering phenomenon, particularly the low-temperature second CDW transition in LaAgSb$_{2}$, which appears as a very weak anomaly in most experiments. The crystal orientation-dependent measurements provide the evolution of Raman modes with crystallographic symmetries and can be further studied through group symmetry analysis. The low-temperature x-ray diffraction data show the emergence of structural modulations corresponding to the CDW instability. The combined high-pressure Raman spectroscopy and synchrotron x-ray diffraction reveal multiple structural phase transitions through lowering of crystalline symmetries, which are also expected to lead to electronic topological transitions.


title: Advances in modeling nuclear matrix elements of neutrinoless double beta decay

abstract: Accurate nuclear matrix elements (NMEs) for neutrinoless double beta decays of candidate nuclei are important for the design and interpretation of future experiments. Significant progress has been made in the modeling of these NMEs from first principles. The NME for 48Ca shows a good agreement among three different ab initio calculations starting from the same nuclear interaction constructed within the chiral EFT and the same decay operator. These studies open the door to ab initio calculations of the matrix elements for the decay of heavier nuclei such as 76Ge, 130Te, and 136Xe. The ultimate goal is the computation of NMEs in many-body calculations with controllable approximations, using nuclear interactions and weak transition operators derived consistently from chiral EFT. We are expecting more progress towards this goal in the near future.


title: The Arveson boundary of a Free Quadrilateral is given by a noncommutative variety

abstract: Let $SM_n(\mathbb{R})^g$ denote $g$-tuples of $n \times n$ real symmetric matrices and set $SM(\mathbb{R})^g = \cup_n SM_n(\mathbb{R})^g$. A free quadrilateral is the collection of tuples $X \in SM(\mathbb{R})^2$ which have positive semidefinite evaluation on the linear equations defining a classical quadrilateral. Such a set is closed under a rich class of convex combinations called matrix convex combination. That is, given elements $X=(X_1, \dots, X_g) \in SM_{n_1}(\mathbb{R})^g$ and $Y=(Y_1, \dots, Y_g) \in SM_{n_2}(\mathbb{R})^g$ of a free quadrilateral $\mathcal{Q}$, one has \[ V_1^T X V_1+V_2^T Y V_2 \in \mathcal{Q} \] for any contractions $V_1:\mathbb{R}^n \to \mathbb{R}^{n_1}$ and $V_2:\mathbb{R}^n \to \mathbb{R}^{n_2}$ satisfying $V_1^T V_1+V_2^T V_2=I_n$. These matrix convex combinations are a natural analogue of convex combinations in the dimension free setting. A natural class of extreme point for free quadrilaterals is free extreme points: elements of a free quadrilateral which cannot be expressed as a nontrivial matrix convex combination of elements of the free quadrilateral. These free extreme points serve as the minimal set which recovers a free quadrilateral through matrix convex combinations. In this article we show that the set of free extreme points of a free quadrilateral is determined by the zero set of a collection of noncommutative polynomials. More precisely, given a free quadrilateral $\mathcal{Q}$, we construct noncommutative polynomials $p_1,p_2,p_3,p_4$ such that a tuple $X \in SM (\mathbb{R})^2$ is a free extreme point of a $\mathcal{Q}$ if and only if $X \in \mathcal{Q}$ and $p_i(X) =0 $ for $i=1,2,3,4$ and $X$ is irreducible. In addition we establish several basic results for projective maps of free spectrahedra and for homogeneous free spectrahedra.


title: Note on Shape Moduli Stabilization, String Gas Cosmology and the Swampland Criteria

abstract: In String Gas Cosmology, the simplest shape modulus fields are naturally stabilized by taking into account the presence of string winding and momentum modes. We determine the resulting effective potential for these fields and show that it obeys the de Sitter conjecture, one of the swampland criteria for effective field theories to be consistent with superstring theory.


title: Differentiation of measures on complete Riemannian manifolds

abstract: In this note we give a new proof of a version of the Besicovitch covering theorem, given in \cite{EG1992}, \cite{Bogachev2007} and extended in \cite{Federer1969}, for locally finite Borel measures on finite dimensional complete Riemannian manifolds $(M,g)$. As a consequence, we prove a differentiation theorem for Borel measures on $(M,g)$, which gives a formula for the Radon-Nikodym density of two nonnegative locally finite Borel measures $\nu_1, \nu_2$ on $(M, g)$ such that $\nu_1 \ll \nu_2$, extending the known case when $(M, g)$ is a standard Euclidean space.


title: Statistical Tree-based Population Seeding for Rolling Horizon EAs in General Video Game Playing

abstract: Multiple Artificial Intelligence (AI) methods have been proposed over recent years to create controllers to play multiple video games of different nature and complexity without revealing the specific mechanics of each of these games to the AI methods. In recent years, Evolutionary Algorithms (EAs) employing rolling horizon mechanisms have achieved extraordinary results in these type of problems. However, some limitations are present in Rolling Horizon EAs making it a grand challenge of AI. These limitations include the wasteful mechanism of creating a population and evolving it over a fraction of a second to propose an action to be executed by the game agent. Another limitation is to use a scalar value (fitness value) to direct evolutionary search instead of accounting for a mechanism that informs us how a particular agent behaves during the rolling horizon simulation. In this work, we address both of these issues. We introduce the use of a statistical tree that tackles the latter limitation. Furthermore, we tackle the former limitation by employing a mechanism that allows us to seed part of the population using Monte Carlo Tree Search, a method that has dominated multiple General Video Game AI competitions. We show how the proposed novel mechanism, called Statistical Tree-based Population Seeding, achieves better results compared to vanilla Rolling Horizon EAs in a set of 20 games, including 10 stochastic and 10 deterministic games.


title: Deep Volumetric Universal Lesion Detection using Light-Weight Pseudo 3D Convolution and Surface Point Regression

abstract: Identifying, measuring and reporting lesions accurately and comprehensively from patient CT scans are important yet time-consuming procedures for physicians. Computer-aided lesion/significant-findings detection techniques are at the core of medical imaging, which remain very challenging due to the tremendously large variability of lesion appearance, location and size distributions in 3D imaging. In this work, we propose a novel deep anchor-free one-stage VULD framework that incorporates (1) P3DC operators to recycle the architectural configurations and pre-trained weights from the off-the-shelf 2D networks, especially ones with large capacities to cope with data variance, and (2) a new SPR method to effectively regress the 3D lesion spatial extents by pinpointing their representative key points on lesion surfaces. Experimental validations are first conducted on the public large-scale NIH DeepLesion dataset where our proposed method delivers new state-of-the-art quantitative performance. We also test VULD on our in-house dataset for liver tumor detection. VULD generalizes well in both large-scale and small-sized tumor datasets in CT imaging.


title: Formulating Module Assessment for Improved Academic Performance Predictability in Higher Education

abstract: Various studies have shown that students tend to get higher marks when assessed through coursework based assessment methods which include either modules that are fully assessed through coursework or a mixture of coursework and examinations than assessed by examination alone. There are a large number of educational data mining studies that preprocess data through conventional data mining processes including data preparation process, but they are using transcript data as they stand without looking at examination and coursework results weighting which could affect prediction accuracy. This paper proposes a different data preparation process through investigating more than 230000 student records in order to prepare students marks based on the assessment methods of enrolled modules. The data have been processed through different stages in order to extract a categorical factor through which students module marks are refined during the data preparation process. The results of this work show that students final marks should not be isolated from the nature of the enrolled modules assessment methods. They must rather be investigated thoroughly and considered during EDMs data preprocessing phases. More generally, it is concluded that educational data should not be prepared in the same way as other data types due to differences as data sources, applications, and types of errors in them. Therefore, an attribute, coursework assessment ratio, is proposed to be used in order to take the different modules assessment methods into account while preparing student transcript data. The effect of CAR on prediction process using the random forest classification technique has been investigated. It is shown that considering CAR as an attribute increases the accuracy of predicting students second year averages based on their first year results.


title: Breaking SU(3) spectral degeneracies in heavy deformed nuclei

abstract: Symmetries are manifested in nature through degeneracies in the spectra of physical systems. In the case of heavy deformed nuclei, when described in the framework of the Interacting Boson Model, within which correlated proton (neutron) pairs are approximated as bosons, the ground state band has no symmetry partner, while the degeneracy between the first excited beta and gamma bands is broken through the use of three-body and/or four-body terms. In the framework of the proxy-SU(3) model, in which an approximate SU(3) symmetry of fermions is present, the same three-body and/or four-body operators are used for breaking the degeneracy between the ground state band and the first excited gamma band. Experimentally accessible quantities being independent of any free parameters are pointed out in the latter case.


title: The exact modulus of the generalized concave Kurdyka-{\L}ojasiewicz property

abstract: We introduce a generalized version of the concave Kurdyka-\L ojasiewicz (KL) property by employing nonsmooth desingularizing functions. We also present the exact modulus of the generalized concave KL property, which provides an answer to the open question regarding the optimal concave desingularizing function. The exact modulus is designed to be the smallest among all possible concave desingularizing functions. Examples are given to illustrate this pleasant property. In turn, using the exact modulus we provide the sharpest upper bound for the total length of iterates generated by the celebrated Bolte-Sabach-Teboulle PALM algorithm.


title: Finite-size effects in the rough phase of the 3d Ising model

abstract: Using Monte Carlo simulations, finite-size effects of interfacial properties in the rough phase of the Ising on a cubic lattice with $L\times L\times R$ sites are studied. In particular, magnetization profiles perpendicular to the flat interface of size L$\times$R are studied, with $L$ being considerably larger than $R$, in the (pre)critical temperature range. The resulting $R$-dependences are compared with predictions of the standard capillary-wave theory, in the Gaussian approximation, and with a field theory based on effective string actions, for $L$=$\infty$.


title: Extraplanar gas in Edge-on Galaxies traced by SOFIA observations of [C II]

abstract: Bursts of localized star formation in galaxies can levitate material from their midplanes. Spiral galaxies that are edge-on allow clear distinction of material that is levitated off the galaxies' midplanes. We used SOFIA to measure the vertical distribution of [C II] 157.7 micron line emission for two nearby, edge-on galaxies, NGC 891 and NGC 5907. We find that for the central region and actively-star-forming regions in the northern portion of NGC 891, and for NGC 5907, a thin (0.3 kpc) disk is supplemented by a thick disk with an exponential scale height of about 2 kpc. The [C II] is far more extended than mid-infrared emission (0.1 kpc, tracing present-day massive star formation) but not as extended as the H I (100 kpc, tracing low-metallicity circum/inter-galactic matter). The extraplanar [C II] may arise in walls of chimneys that connect the disk to the halo.


title: On the non-existence of extended perfect codes and some perfect colorings

abstract: In this paper we obtain the necessary condition for the existence of perfect $k$-colorings (equitable $k$-partitions) in Hamming graphs $H(n,q)$, where $q=2,3,4$ and Doob graphs $D(m,n)$. As an application, we prove the non-existence of extended perfect codes in $H(n,q)$, where $q=3,4$, $n>q+2$, and in $D(m,n)$, where $2m+n>6$.


title: Benchmarking adversarial attacks and defenses for time-series data

abstract: The adversarial vulnerability of deep networks has spurred the interest of researchers worldwide. Unsurprisingly, like images, adversarial examples also translate to time-series data as they are an inherent weakness of the model itself rather than the modality. Several attempts have been made to defend against these adversarial attacks, particularly for the visual modality. In this paper, we perform detailed benchmarking of well-proven adversarial defense methodologies on time-series data. We restrict ourselves to the $L_{\infty}$ threat model. We also explore the trade-off between smoothness and clean accuracy for regularization-based defenses to better understand the trade-offs that they offer. Our analysis shows that the explored adversarial defenses offer robustness against both strong white-box as well as black-box attacks. This paves the way for future research in the direction of adversarial attacks and defenses, particularly for time-series data.


title: LinkRing: A Wearable Haptic Display for Delivering Multi-contact and Multi-modal Stimuli at the Finger Pads

abstract: LinkRing is a novel wearable tactile display for providing multi-contact and multi-modal stimuli at the finger. The system of two five-bar linkage mechanisms is designed to operate with two independent contact points, which combined can provide such stimulation as shear force and twist stimuli, slippage, and pressure. The proposed display has a lightweight and easy to wear structure. Two experiments were carried out in order to determine the sensitivity of the finger surface, the first one aimed to determine the location of the contact points, and the other for discrimination the slippage with varying rates. The results of the experiments showed a high level of pattern recognition.


title: Discrete Lebedev-Skalskaya transforms

abstract: Discrete analogs of the Lebedev-Skalskaya transforms are introduced and investigated. It involves series and integrals with respect to the kernels ${\rm Re} K_{\alpha+in}(x), {\rm Im} K_{\alpha+in}(x), x >0, n \in \mathbb{N}, |\alpha | < 1,\ i $ is the imaginary unit and $K_\nu(z)$ is the modified Bessel function. The corresponding inversion formulas for suitable functions and sequences in terms of these series and integrals are established when $\alpha = \pm 1/2$. The case $\alpha=0$ reduces to the Kontorovich-Lebedev transform.


title: High-Throughput Virtual Screening of 4487 flavonoids: New insights on the structural inhibition of SARS-CoV-2 Main Protease

abstract: COVID-19 presents a great threat to public health worldwide and the infectious agent SARS-CoV-2 is currently the target of much research aiming at inhibition. The virus' main protease is a dimeric enzyme that has only recently begun to be thoroughly described, opening the door for virtual screening more broadly. Here, a PAIN-filtered flavonoid database was screened against four sites of the protease: a free (normal) conformation of the Substrate Binding Site (NSBS), an induced-fit state of the SBS (ISBS), a Dimerization Site (DS) and a Cryptic Site (CS). The mean binding energies of the top five ligands from each site were -9.52, -11.512, -7.042 and -10.348 kcal/mol for the NSBS, the ISBS, the DS and the CS, respectively. For the DS and CS, these top five compounds were selected as candidates to bind their respective site. In the case of SBS, the top 30 ligands with the lowest binding energies from NSBS and ISBS were contrasted and the ones present in both lists were selected as the final candidates. The final list was: Dorsilurin E (FL3FQUNP0001), Euchrenone a11 (FL2FALNP0014), Kurziflavolactone C (FL2FA9NC0016), Licorice glycoside E (FL2F1AGSN001) and Taxifolin 3'- (6"-phenyl- acetylglucoside) (FL4DACGS0020) for the SBS; Sanggenol O (FL2FALNP0020), CHEMBL2171573, Kanzonol E (FL3F1ANP0001), CHEMBL2171584 and Abyssynoflavanone VI (FL2FACNP0014) for DS and CHEMBL2171598, CHEMBL2171577, Denticulaflavanol (FL5FAANR0001), Kurzichalcolactone (FL1CA9NC0001) and CHEMBL2171578 for CS. Virtual screening integrated several confirmation methods, including cross-docking assays and positive and negative controls. All 15 compounds are currently subjected to molecular dynamics so as to theoretically validate their binding to the protease.


title: A Meta-Learning Control Algorithm with Provable Finite-Time Guarantees

abstract: In this work we provide provable regret guarantees for an online meta-learning control algorithm in an iterative control setting, where in each iteration the system to be controlled is a linear deterministic system that is different and unknown, the cost for the controller in an iteration is a general additive cost function and the control input is required to be constrained, which if violated incurs an additional cost. We prove (i) that the algorithm achieves a regret for the controller cost and constraint violation that are $O(T^{3/4})$ for an episode of duration $T$ with respect to the best policy that satisfies the control input control constraints and (ii) that the average of the regret for the controller cost and constraint violation with respect to the same policy vary as $O((1+\log(N)/N)T^{3/4})$ with the number of iterations $N$, showing that the worst regret for the learning within an iteration continuously improves with experience of more iterations.