title: Two-Stage Fourth-order Gas-kinetic Scheme for Three-dimensional Euler and Navier-Stokes Solutions

abstract: For the one-stage third-order gas-kinetic scheme (GKS), success applications have been achieved for the three-dimensional compressible flow computations [33]. The high-order accuracy of the scheme is obtained directly by integrating a multidimensional time-accurate gas distribution function over the cell interface within a time step without implementing Gaussian quadrature points and Runge-Kutta time-stepping technique. However, for the further increasing the order of the scheme, such as the fourth-order one, the formulation becomes very complicated for the multidimensional flow. Recently, a two-stage fourth-order GKS with high efficiency has been constructed for two-dimensional inviscid and viscous flow computations [22,32], and the scheme uses the time accurate flux function and its time derivatives. In this paper, a fourth-order GKS is developed for the threedimensional flows under the two-stage framework. Based on the three-dimensional WENO reconstruction and flux evaluation at Gaussian quadrature points on a cell interface, the high-order accuracy in space is achieved first. Then, the two-stage time stepping method provides the high accuracy in time. In comparison with the formal third-order GKS [33], the current fourth-order method not only improves the accuracy of the scheme, but also reduces the complexity of the gas-kinetic solver greatly. More importantly, the fourth-order GKS has the same robustness as the second-order shock capturing scheme. This scheme is applied to both inviscid and viscous, and low and high speed flow computations. Numerical results validate the outstanding reliability and applicability of the scheme for three-dimensional flows, such as turbulent one.


title: Introduction to total dominator edge chromatic number

abstract: We introduce the total dominator edge chromatic number of a graph $G$. A total dominator edge coloring (briefly TDE-coloring) of $G$ is a proper edge coloring of $G$ in which each edge of the graph is adjacent to every edge of some color class. The total dominator edge chromatic number (briefly TDEC-number) $\chi'^t_d(G)$ of $G$ is the minimum number of color classes in a TDE-coloring of $G$. We obtain some properties of $\chi'^t_d(G)$ and compute this parameter for specific graphs. We examine the effects on $\chi'^t_d(G)$ when $G$ is modified by operations on vertex and edge of $G$. Finally, we consider the $k$-subdivison of $G$ and study TDEC-number of this kind of graphs.


title: Bardeen regular black hole as a quantum-corrected Schwarzschild black hole

abstract: Bardeen regular black hole is commonly considered as a solution of general relativity coupled to a nonlinear electrodynamics. In this paper, it is shown that the Bardeen solution may be interpreted as a quantum-corrected Schwarzschild black hole. This new interpretation is obtained by means of a generalized uncertainty principle applied to the Hawking temperature. Moreover, using the regular black hole of Bardeen, it is possible to evaluate the quantum gravity parameter of the generalized uncertainty principle or, assuming the recent upper bounds for such a parameter, to verify an enormous discrepancy between a cosmological constant and that measured by recent cosmological observations $(\sim 10^{120})$.


title: Mirrored and Hybrid Disk Arrays: Organization, Scheduling, Reliability, and Performance

abstract: Basic mirroring (BM) classified as RAID level 1 replicates data on two disks, thus doubling disk access bandwidth for read requests. RAID1/0 is an array of BM pairs with balanced loads due to striping. When a disk fails the read load on its pair is doubled, which results in halving the maximum attainable bandwidth. We review RAID1 organizations which attain a balanced load upon disk failure, but as shown by reliability analysis tend to be less reliable than RAID1/0. Hybrid disk arrays which store XORed instead of replicated data tend to have a higher reliability than mirrored disks, but incur a higher overhead in updating data. Read request response time can be improved by processing them at a higher priority than writes, since they have a direct effect on application response time. Shortest seek distance and affinity based routing both shorten seek time. Anticipatory arm placement places arms optimally to minimize the seek distance. The analysis of RAID1 in normal, degraded, and rebuild mode is provided to quantify RAID1/0 performance. We compare the reliability of mirrored disk organizations against each other and hybrid disks and erasure coded disk arrays.


title: A better communicator is always a better scientist, or the reason why every PhD student should engage in science outreach

abstract: The ability to communicate with all audiences is a skill that is rapidly becoming a must-have for any future scientist. As more physicists engage in communicating science to non-expert audiences, research shows that this experience helps them to get a better understanding of their own research and the impact on society, improves the perception of science by lay audiences and can also become an area of personal growth as a citizen. A recent deployment of a PhD student to the Amundsen Scott South Pole Station, as part of the IceCube Collaboration, provided a ready opportunity to spark interest. We present results of the efforts made by the Universit\'e libre de Bruxelles (ULB), the Vrije Universiteit Brussel (VUB) and the Interuniversity Institute for High Energies, IIHE (ULB-VUB), to introduce Belgian students and citizens to science and the life of a scientist. The essential parts of this program will be identified to show why the contributions of a PhD student to the organization of these activities are beneficial to the development of new skills as a scientist, but also to broaden the audiences and the impact of the local university and/or the specific research outreach program.


title: Refinement of Parameters of Six Selected Galactic Potential Models

abstract: This paper is devoted to the refinement of the parameters of the six three-component (bulge, disk, halo) axisymmetric Galactic gravitational potential models on the basis of modern data on circular velocities of Galactic objects located at distances up to 200 kpc from the Galactic center. In all models the bulge and disk are described by the Miyamoto--Nagai expressions. To describe the halo, the models of Allen-Santill\'an (I), Wilkinson-Evans (II), Navarro-Frenk-White (III), Binney (IV), Plummer (V), and Hernquist (VI) are used. The sought-for parameters of potential models are determined by fitting the model rotation curves to the measured velocities, taking into account restrictions on the local dynamical matter density $\rho_\odot=0.1 M_\odot$ pc$^{-3}$ and the vertical force $|K_{z=1.1}|/2\pi G=77 M_\odot $ pc$^{-2}$. A comparative analysis of the refined potential models is made and for each of the models the estimates of a number of the Galactic characteristics are presented.


title: On the algorithmic complexity of decomposing graphs into regular/irregular structures

abstract: A locally irregular graph is a graph whose adjacent vertices have distinct degrees, a regular graph is a graph where each vertex has the same degree and a locally regular graph is a graph where for every two adjacent vertices u, v, their degrees are equal. In this work, we study the set of all problems which are related to decomposition of graphs into regular, locally regular and/or locally irregular subgraphs and we present some polynomial time algorithms, NP-completeness results, lower bounds and upper bounds for them. Among our results, one of our lower bounds makes use of mutually orthogonal Latin squares which is relatively novel.


title: A Curie-Weiss Theory of the Continuum Widom-Rowlinson Model

abstract: A version of the continuum Widom-Rowlinson model is introduced and studied. It is a two-component gas of point particles placed in $\mathbf{R}^d$ in which like particles do not interact and unlike particles contained in a given vessel of volume $V$ repel each other with intensity $a/V$. This model is thermodynamically equivalent to a one-component gas with multi-particle interaction. For both models, a rigorous theory of a phase transition is presented and the ways of its construction in the framework of the grand canonical formalism are outlined.


title: The Structure of Martian Magnetosphere at the Dayside Terminator Region as Observed on MAVEN Spacecraft

abstract: We analyzed 44 passes of the MAVEN spacecraft through the magnetosphere, arranged by the angle between electric field vector and the projection of spacecraft position radius vector in the YZ plane in MSE coordinate system (${\theta}$ E ). All passes were divided into 3 angular sectors near 0{\deg}, 90{\deg} and 180{\deg} ${\theta}$ E angles in order to estimate the role of IMF direction in plasma and magnetic properties of dayside Martian magnetosphere. The time interval chosen was from January 17 through February 4, 2016 when MAVEN was crossing the dayside magnetosphere at SZA ~ 70{\deg}. Magnetosphere as the region with prevailing energetic planetary ions is always found between the magnetosheath and the ionosphere. 3 angular sectors of dayside interaction region in MSE coordinate system with different orientation of the solar wind electric field vector E = -1/c V x B showed that for each sector one can find specific profiles of the magnetosheath, the magnetic barrier and the magnetosphere. Plume ions originate in the northern MSE sector where motion electric field is directed from the planet. This electric field ejects magnetospheric ions leading to dilution of magnetospheric heavy ions population, and this effect is seen in some magnetospheric profiles. Magnetic barrier forms in front of the magnetosphere, and relative magnetic field magnitudes in these two domains vary. The average height of the boundary with ionosphere is ~530 km and the average height of the magnetopause is ~730 km. We discuss the implications of the observed magnetosphere structure to the planetary ions loss mechanism.


title: Purcell Effect in the Stimulated and Spontaneous Emission Rates of Nanoscale Semiconductor Lasers

abstract: Nanoscale semiconductor lasers have been developed recently using either metal, metallo-dielectric or photonic crystal nanocavities. While the technology of nanolasers is steadily being deployed, their expected performance for on-chip optical interconnects is still largely unknown due to a limited understanding of some of their key features. Specifically, as the cavity size is reduced with respect to the emission wavelength, the stimulated and the spontaneous emission rates are modified, which is known as the Purcell effect in the context of cavity quantum electrodynamics. This effect is expected to have a major impact in the 'threshold-less' behavior of nanolasers and in their modulation speed, but its role is poorly understood in practical laser structures, characterized by significant homogeneous and inhomogeneous broadening and by a complex spatial distribution of the active material and cavity field. In this work, we investigate the role of Purcell effect in the stimulated and spontaneous emission rates of semiconductor lasers taking into account the carriers' spatial distribution in the volume of the active region over a wide range of cavity dimensions and emitter/cavity linewidths, enabling the detailed modeling of the static and dynamic characteristics of either micro- or nano-scale lasers using single-mode rate-equations analysis. The ultimate limits of scaling down these nanoscale light sources in terms of Purcell enhancement and modulation speed are also discussed showing that the ultrafast modulation properties predicted in nanolasers are a direct consequence of the enhancement of the stimulated emission rate via reduction of the mode volume.


title: Marked clustering statistics in $f(R)$ gravity cosmologies

abstract: We analyse the two-point and marked correlation functions of haloes and galaxies in three variants of the chameleon $f(R)$ gravity model using N-body simulations, and compare to a fiducial $\Lambda$CDM model based on general relativity (GR). Using a halo occupation distribution prescription (HOD) we populate dark matter haloes with galaxies, where the HOD parameters have been tuned such that the galaxy number densities and the real-space galaxy two-point correlation functions in the modified gravity models match those in GR to within $1\sim3\%$. We test the idea that since the behaviour of gravity is dependent on environment, marked correlation functions may display a measurable difference between the models. For this we test marks based on the density field and the Newtonian gravitational potential. We find that the galaxy marked correlation function shows significant differences measured in different models on scales smaller than $r\lesssim 20~h^{-1}$ Mpc. Guided by simulations to identify a suitable mark, this approach could be used as a new probe of the accelerated expansion of the Universe.


title: Correlated Components Analysis - Extracting Reliable Dimensions in Multivariate Data

abstract: How does one find dimensions in multivariate data that are reliably expressed across repetitions? For example, in a brain imaging study one may want to identify combinations of neural signals that are reliably expressed across multiple trials or subjects. For a behavioral assessment with multiple ratings, one may want to identify an aggregate score that is reliably reproduced across raters. Correlated Components Analysis (CorrCA) addresses this problem by identifying components that are maximally correlated between repetitions (e.g. trials, subjects, raters). Here we formalize this as the maximization of the ratio of between-repetition to within-repetition covariance. We show that this criterion maximizes repeat-reliability, defined as mean over variance across repeats, and that it leads to CorrCA or to multi-set Canonical Correlation Analysis, depending on the constraints. Surprisingly, we also find that CorrCA is equivalent to Linear Discriminant Analysis for zero-mean signals, which provides an unexpected link between classic concepts of multivariate analysis. We present an exact parametric test of statistical significance based on the F-statistic for normally distributed independent samples, and present and validate shuffle statistics for the case of dependent samples. Regularization and extension to non-linear mappings using kernels are also presented. The algorithms are demonstrated on a series of data analysis applications, and we provide all code and data required to reproduce the results.


title: Symmetric Polynomials in Upper-bound Semirings

abstract: The fundamental theorem of symmetric polynomials over rings is a classical result which states that every unital commutative ring is fully elementary, i.e. we can express symmetric polynomials with elementary ones in a unique way. The result does not extend directly to polynomials over semirings, but we do have analogous results for some special semirings, for example, the tropical, extended and supertropical semirings. These all fall into a larger class of upper-bound semirings. In this paper we extend the known results and give a complete characterization of fully elementary upper-bound semirings. We further improve this characterization statement in the case of linearly ordered upper-bound semirings.


title: Resonant-state expansion for open optical systems: Generalization to magnetic, chiral, and bi-anisotropic materials

abstract: The resonant-state expansion, a recently developed powerful method in electrodynamics, is generalized here for open optical systems containing magnetic, chiral, or bi-anisotropic materials. It is shown that the key matrix eigenvalue equation of the method remains the same, but the matrix elements of the perturbation now contain variations of the permittivity, permeability, and bi-anisotropy tensors. A general normalization of resonant states in terms of the electric and magnetic fields is presented.


title: First Demonstration of a Pixelated Charge Readout for Single-Phase Liquid Argon Time Projection Chambers

abstract: Liquid Argon Time Projection Chambers (LArTPCs) have been selected for the future long-baseline Deep Underground Neutrino Experiment (DUNE). To allow LArTPCs to operate in the high-multiplicity near detector environment of DUNE, a new charge readout technology is required. Traditional charge readout technologies introduce intrinsic ambiguities, combined with a slow detector response, these ambiguities have limited the performance of LArTPCs, until now. Here, we present a novel pixelated charge readout that enables the full 3D tracking capabilities of LArTPCs. We characterise the signal to noise ratio of charge readout chain, to be about 14, and demonstrate track reconstruction on 3D space points produced by the pixel readout. This pixelated charge readout makes LArTPCs a viable option for the DUNE near detector complex.


title: Fractional powers and singular perturbations of quantum differential Hamiltonians

abstract: We consider the fractional powers of singular (point-like) perturbations of the Laplacian, and the singular perturbations of fractional powers of the Laplacian, and we compare such two constructions focusing on their perturbative structure for resolvents and on the local singularity structure of their domains. In application to the linear and non-linear Schr\"{o}dinger equations for the corresponding operators we outline a programme of relevant questions that deserve being investigated.


title: First-principles study of (Ba,Ca)TiO$_3$ and Ba(Ti,Zr)O$_3$ solid solutions

abstract: (Ba,Ca)TiO$_3$ and Ba(Ti,Zr)O$_3$ solid solutions are the building blocks of lead-free piezoelectric materials that attract a renewed interest. We investigate the properties of these systems by means of first-principles calculations, with a focus on the lattice dynamics and the competition between different ferroelectric phases. We first analyze the four parent compounds in order to compare their properties and their different tendency towards ferroelectricity. The core of our study is systematic characterization of the binary systems (Ba,Ca)TiO$_3$ and Ba(Ti,Zr)O$_3$ within both the Virtual Crystal Approximation and direct supercell calculations. In the case of Ca doping, we find a gradual transformation from $A$-site to $B$-site ferroelectricity due to steric effects that largely determines the behavior of the system. In the case of Zr doping, in contrast, the behavior is eventually dominated by cooperative Zr-Ti motions and the local electrostatics. In addition, our comparative study reveals that the specific microscopic physics of these solid sets severe limits to the applicability of the Virtual Crystal Approximation for these systems.


title: Neutrino electromagnetic properties: a window to new physics - II

abstract: There is merely a short note on the selected issues of neutrino electromagnetic properties with focus on effects of new physics. The meaning of "new physics" is twofold: 1) a massive neutrino have nonzero electromagnetic properties that can be considered as manifestation of new physics beyond the Standard Model, and 2) in studies of neutrinos electromagnetic interactions new effects are predicted that can lead to new phenomena accessible for observations.


title: Antisymetric matrices and Lotka-Volterra algebras

abstract: The purpose of this paper is to study the structure of Lotka-Volterra algebras, the set of their idempotent elements and their group of automorphisms. These algebras are defined through antisymmetric matrices and they emerge in connection with biological problems and Lotka-Volterra systems for the interactions of neighboring individuals.


title: Non reflection and perfect reflection via Fano resonance in waveguides

abstract: We investigate a time-harmonic wave problem in a waveguide. By means of asymptotic analysis techniques, we justify the so-called Fano resonance phenomenon. More precisely, we show that the scattering matrix considered as a function of a geometrical parameter $\varepsilon$ and of the frequency $\lambda$ is in general not continuous at a point $(\varepsilon,\lambda)=(0,\lambda^0)$ where trapped modes exist. In particular, we prove that for a given $\varepsilon\ne0$ small, the scattering matrix exhibits a rapid change for frequencies varying in a neighbourhood of $\lambda^0$. We use this property to construct examples of waveguides such that the energy of an incident wave propagating through the structure is perfectly transmitted (non reflection) or perfectly reflected in monomode regime. We provide numerical results to illustrate our theorems.


title: Cosmic ray production in supernovae

abstract: We give a brief review of the origin and acceleration of cosmic rays (CRs), emphasizing the production of CRs at different stages of supernova evolution by the first-order Fermi shock acceleration mechanism. We suggest that supernovae with trans-relativistic outflows, despite being rather rare, may accelerate CRs to energies above 10$^{18}$ eV over the first year of their evolution. Supernovae in young compact clusters of massive stars, and interaction powered superluminous supernovae, may accelerate CRs well above the PeV regime. We discuss the acceleration of the bulk of the galactic CRs in isolated supernova remnants and re-acceleration of escaped CRs by the multiple shocks present in superbubbles produced by associations of OB stars. The effects of magnetic field amplification by CR driven instabilities, as well as superdiffusive CR transport, are discussed for nonthermal radiation produced by nonlinear shocks of all speeds including trans-relativistic ones.


title: Linear build-up of Fano resonance spectral profiles

abstract: The build-up dynamics of a continuous spectrum under the action of a weak laser field on a Fano resonance with the use of the pulses with the Lorentz spectrum and ultrashort pulses in the wavelet form is investigated. A dispersion-time excitation dependence of the Fano resonances in a He atom, in an InP impurity semiconductor, in longitudinal optical LO-phonons of a shallow donor exciton in pure ZnO crystals, and in metamaterials are calculated. The numerical simulation of the dynamics has shown time-dependent formation a Fano spectral profile in the systems of different physical natures under the action of ultrashort pulses with attosecond and femtosecond durations.


title: Comparison Between Robust and Stochastic Optimisation for Long-term Reservoir Operations Under Uncertainty

abstract: Long-term reservoir management often uses bounds on the reservoir level, between which the operator can work. However, these bounds are not always kept up-to-date with the latest knowledge about the reservoir drainage area, and thus become obsolete. The main difficulty with bounds computation is to correctly take into account the high uncertainty about the inflows to the reservoir. In this article, we propose a methodology to derive minimum bounds while providing formal guarantees about the quality of the obtained solutions. The uncertainty is embedded using either stochastic or robust programming in a model-predictive-control framework. We compare the two paradigms to the existing solution for a case study and find that the obtained solutions vary substantially. By combining the stochastic and the robust approaches, we also assign a confidence level to the solutions obtained by stochastic programming. The proposed methodology is found to be both efficient and easy to implement. It relies on sound mathematical principles, ensuring that a global optimum is reached in all cases.


title: Structural, Electronic, and Dynamical Properties of Liquid Water by ab initio Molecular Dynamics based on SCAN Functional within the Canonical Ensemble

abstract: We perform ab initio molecular dynamics (AIMD) simulation of liquid water in the canonical ensemble at ambient conditions using the SCAN meta-GGA functional approximation, and carry out systematic comparisons with the results obtained from the GGA-level PBE functional, and Tkatchenko-Scheffler van der Waals (vdW) dispersion correction inclusive PBE functional. We analyze various properties of liquid water including radial distribution functions, oxygen-oxygen-oxygen triplet angular distribution, tetrahedrality, hydrogen bonds, diffusion coefficients, ring statistics, density of states, band gaps, and dipole moments. We find that the SCAN functional is generally more accurate than the other two functionals for liquid water by not only capturing the intermediate-range vdW interactions but also mitigating the overly strong hydrogen bonds prescribed in PBE simulations. We also compare the results of SCAN-based AIMD simulations in the canonical and isothermal-isobaric ensembles. Our results suggest that SCAN provides a reliable description for most structural, electronic, and dynamical properties in liquid water.

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


title: Entanglement sudden death and revival in quantum dark-soliton qubits

abstract: We study the finite time entanglement dynamics between two dark soliton qubits due to quantum fluctuations in a quasi one dimensional Bose-Einstein condensates. Recently, dark solitons are proved to be an appealing platform for qubits due to their appreciably long life time. We explore the entanglement decay for an entangled state of two phonon coherences and the qubits to be in the diagonal basis of so called Dicke states. We observe the collapse and revival of the entanglement, depending critically on the collective damping term but independent of the qubit-qubit interaction for both the initial states. The collective behavior of the dark soliton qubits demonstrate the dependence of entanglement evolution on the interatomic distance.


title: On possible application of spin light of neutrino in astrophysics

abstract: The $ spin $ $ light$ $ of $ $neutrino $ ($SL\nu$) is a phenomenon of electromagnetic radiation by a massive neutrino moving in external media that is originated owing to neutrino magnetic moment. In this short paper we note on the importance of this effect in the light of its connection with the neutrino magnetic moment, recap its basic properties in dense matter and give some general criteria for its best efficiency in nature. On this basis we propose a set of possible astrophysical environments where the $SL\nu$ can be manifested in principle.


title: Neutrino quantum decoherence due to entanglement with magnetic field

abstract: The phenomena of neutrino oscillations emerges due to coherent superposition of different neutrino states. The entanglement of neutrinos with its environment can lead to a suppression of neutrino oscillations. The master equation for neutrino evolution in a magnetic field is derived taking into account the entanglement with a magnetic field.


title: Thermal desorption of formamide and methylamine from graphite and amorphous water ice surfaces

abstract: Formamide (NH2CHO) and methylamine (CH3NH2) are known to be the most abundant amine-containing molecules in many astrophysical environments. The presence of these molecules in the gas phase may result from thermal desorption of interstellar ices. The aim of this work is to determine the values of the desorption energies of formamide and methylamine from analogues of interstellar dust grain surfaces and to understand their interaction with water ice. TPD experiments of formamide and methylamine ices were performed in the submonolayer and monolayer regimes on graphite (HOPG) and non-porous amorphous solid water ice surfaces at 40-240 K. The desorption energy distributions of these two molecules were calculated from TPD measurements using a set of independent Polanyi-Wigner equations. The maximum of the desorption of formamide from both graphite and ASW ice surfaces occurs at 176 K after the desorption of H2O molecules, whereas the desorption profile of methylamine depends strongly on the substrate. Solid methylamine starts to desorb below 100 K from the graphite surface. Its desorption from the water ice surface occurs after 120 K and stops during the water ice sublimation around 150 K. It continues to desorb from the graphite surface at temperatures higher than 160 K. More than 95 % of NH2CHO diffuses through the water ice surface towards the graphitic substrate and is released into the gas phase with a desorption energy distribution (7460-9380 K), which is measured with the best-fit pre-exponential factor A=10^18 s-1. However, the desorption energy distribution of methylamine from the np-ASW ice surface (3850-8420 K) is measured with the best-fit pre-exponential factor A=10^12 s-1. A fraction of solid methylamine monolayer of roughly 0.15 diffuses through the water ice surface towards the HOPG substrate. This small amount desorbs later with higher binding energies (5050-8420 K).


title: A note on "MLE in logistic regression with a diverging dimension"

abstract: This short note is to point the reader to notice that the proof of high dimensional asymptotic normality of MLE estimator for logistic regression under the regime $p_n=o(n)$ given in paper: "Maximum likelihood estimation in logistic regression models with a diverging number of covariates. Electronic Journal of Statistics, 6, 1838-1846." is wrong.


title: A new relativistic component of the accretion disk wind in PDS 456

abstract: Past X-ray observations of the nearby luminous quasar PDS 456 (at $z=0.184$) have revealed a wide angle accretion disk wind (Nardini et al. 2015), with an outflow velocity of $\sim-0.25c$. Here we unveil a new, relativistic component of the wind through hard X-ray observations with NuSTAR and XMM-Newton, obtained in March 2017 when the quasar was in a low flux state. This very fast wind component, with an outflow velocity of $-0.46\pm0.02c$, is detected in the iron K band, in addition to the $-0.25c$ wind zone. The relativistic component may arise from the innermost disk wind, launched from close to the black hole at radius of $\sim10$ gravitational radii. The opacity of the fast wind also increases during a possible obscuration event lasting for 50 ks. We suggest that the very fast wind may only be apparent during the lowest X-ray flux states of PDS 456, becoming overly ionized as the luminosity increases. Overall, the total wind power may even approach the Eddington value.


title: Topological aspect of monodromy groupoid for a group-groupoid

abstract: In this paper we develop star topological and topological group-groupoid structures of monodromy groupoid and prove that the monodromy groupoid of a topological group-groupoid is also a topological group-groupoid.


title: Detecting Changes in Fully Polarimetric SAR Imagery with Statistical Information Theory

abstract: Images obtained from coherent illumination processes are contaminated with speckle. A prominent example of such imagery systems is the polarimetric synthetic aperture radar (PolSAR). For such remote sensing tool the speckle interference pattern appears in the form of a positive definite Hermitian matrix, which requires specialized models and makes change detection a hard task. The scaled complex Wishart distribution is a widely used model for PolSAR images. Such distribution is defined by two parameters: the number of looks and the complex covariance matrix. The last parameter contains all the necessary information to characterize the backscattered data and, thus, identifying changes in a sequence of images can be formulated as a problem of verifying whether the complex covariance matrices differ at two or more takes. This paper proposes a comparison between a classical change detection method based on the likelihood ratio and three statistical methods that depend on information-theoretic measures: the Kullback-Leibler distance and two entropies. The performance of these four tests was quantified in terms of their sample test powers and sizes using simulated data. The tests are then applied to actual PolSAR data. The results provide evidence that tests based on entropies may outperform those based on the Kullback-Leibler distance and likelihood ratio statistics.


title: Product interactions and feedback in diffusion-controlled reactions

abstract: Steric or attractive interactions among reactants or between reactants and inert crowders can substantially influence the total rate of a diffusion-influenced reaction in the liquid phase. However, the role of the product species, that has typically different physical properties than the reactant species, has been disregarded so far. Here we study the effects of reactant-product and product-product interactions as well as asymmetric diffusion properties on the rate of diffusion-controlled reactions in the classical Smoluchowski-setup for chemical transformations at a perfect catalytic sphere. For this we solve the diffusion equation with appropriate boundary conditions coupled by a mean-field approach on the second virial level to account for the particle interactions. We find that all particle spatial distributions and the total rate can change significantly, depending on the diffusion and interaction properties of the accumulated products. Complex competing and self-regulating (homeostatic) or self-amplifying effects are observed for the system, leading to both decrease and increase of the rates, as the presence of interacting products feeds back to the reactant flux and thus the rate with which the products are generated.


title: Covering morphisms of internal groupoids in the models of a semi-abelian theory

abstract: In this paper, for given an algebraic theory $T$ whose category $C$ of models is semi-abelian, we consider the topological models of $T$ called topological $T$-algebras and obtain some results related to the fundamental groups of topological $T$-algebras. We also deal with the internal groupoid structure in the category of models providing that the fundamental groupoid deduces a functor from topological $T$-algebras to the internal groupoids in $C$ and prove a criterion for the lifting of such an internal groupoid structure to the covering groupoids.


title: Maximum principle and its application for the nonlinear time-fractional diffusion equations with Cauchy-Dirichlet conditions

abstract: In this paper, a maximum principle for the one-dimensional sub-diffusion equation with Atangana-Baleanu fractional derivative is formulated and proved. The proof of the maximum principle is based on an extremum principle for the Atangana-Baleanu fractional derivative that is given in the paper, too. The maximum principle is then applied to show that the initial-boundary-value problem for the linear and nonlinear time-fractional diffusion equations possesses at most one classical solution and this solution continuously depends on the initial and boundary conditions.


title: On Motzkin numbers and central trinomial coefficients

abstract: The Motzkin numbers $M_n=\sum_{k=0}^n\binom n{2k}\binom{2k}k/(k+1)$ $(n=0,1,2,\ldots)$ and the central trinomial coefficients $T_n$ ($n=0,1,2,\ldots)$ given by the constant term of $(1+x+x^{-1})^n$, have many combinatorial interpretations. In this paper we establish the following surprising arithmetic properties of them with $n$ any positive integer: $$\frac2n\sum_{k=1}^n(2k+1)M_k^2\in\mathbb Z,$$ $$\frac{n^2(n^2-1)}6\,\bigg|\,\sum_{k=0}^{n-1}k(k+1)(8k+9)T_kT_{k+1},$$ and also $$\sum_{k=0}^{n-1}(k+1)(k+2)(2k+3)M_k^23^{n-1-k}=n(n+1)(n+2)M_nM_{n-1}.$$


title: Uncovering the mechanism of the impurity-selective Mott transition in paramagnetic V$_{2}$O$_{3}$

abstract: While the phase diagrams of the one- and multi-orbital Hubbard model have been well studied, the physics of real Mott insulators is often much richer, material dependent, and poorly understood. In the prototype Mott insulator V$_{2}$O$_{3}$, chemical pressure was initially believed to explain why the paramagnetic-metal to antiferromagnetic-insulator transition temperature is lowered by Ti doping while Cr doping strengthens correlations, eventually rendering the high-temperature phase paramagnetic insulating. However, this scenario has been recently shown both experimentally and theoretically to be untenable. Based on full structural optimization, we demonstrate via the charge self-consistent combination of density functional theory and dynamical mean-field theory that changes in the V$_{2}$O$_{3}$ phase diagram are driven by defect-induced local symmetry breakings resulting from dramatically different couplings of Cr and Ti dopants to the host system. This finding emphasizes the high sensitivity of the Mott metal-insulator transition to the local environment and the importance of accurately accounting for the one-electron Hamiltonian, since correlations crucially respond to it.


title: Molecular Dynamics of Lithium Ion Transport in a Model Solid Electrolyte Interphase

abstract: Li+ transport within a solid electrolyte interphase (SEI) in lithium ion batteries has challenged molecular dynamics (MD) studies due to limited compositional control of that layer. In recent years, experiments and ab initio simulations have identified dilithium ethylene dicarbonate (Li2EDC) as the dominant component of SEI layers. Here, we adopt a parameterized, non-polarizable MD force field for Li2EDC to study transport characteristics of Li+ in this model SEI layer at moderate temperatures. The observed correlations are consistent with recent MD results using a polarizable force field, suggesting that this non-polarizable model is effective for our purposes of investigating Li+ dynamics over long time scales. Mean-squared displacements distinguish three distinct Li+ transport regimes in EDC ballistic, trapping, and diffusive. Compared to liquid ethylene carbonate (EC), the nanosecond trapping times in EDC are significantly longer and naturally decrease at higher temperatures. New materials developed for fast-charging Li-ion batteries should have smaller trapping regions. The analyses implemented in this paper can be used for testing transport of Li+ ion in novel battery materials. Non-Gaussian features of van Hove self -correlation functions for Li+ in EDC, along with the mean-squared displacements, are consistent in describing EDC as a glassy material compared with liquid EC. Vibrational modes of Li+ ion, identified by MD, characterize the trapping and are further validated by electronic structure calculations.


title: R-matrix-valued Lax pairs and long-range spin chains

abstract: In this paper we discuss $R$-matrix-valued Lax pairs for ${\rm sl}_N$ Calogero-Moser model and their relation to integrable quantum long-range spin chains of the Haldane-Shastry-Inozemtsev type. First, we construct the $R$-matrix-valued Lax pairs for the third flow of the classical Calogero-Moser model. Then we notice that the scalar parts (in the auxiliary space) of the $M$-matrices corresponding to the second and third flows have form of special spin exchange operators. The freezing trick restricts them to quantum Hamiltonians of long-range spin chains. We show that for a special choice of the $R$-matrix these Hamiltonians reproduce those for the Inozemtsev chain. In the general case related to the Baxter's elliptic $R$-matrix we obtain a natural anisotropic extension of the Inozemtsev chain. Commutativity of the Hamiltonians is verified numerically. Trigonometric limits lead to the Haldane-Shastry chains and their anisotropic generalizations.

categories: math-ph cond-mat.str-el hep-th math.MP nlin.SI


title: Antarctic Surface Reflectivity Calculations and Measurements from the ANITA-4 and HiCal-2 Experiments

abstract: The balloon-borne HiCal radio-frequency (RF) transmitter, in concert with the ANITA radio-frequency receiver array, is designed to measure the Antarctic surface reflectivity in the RF wavelength regime. The amplitude of surface-reflected transmissions from HiCal, registered as triggered events by ANITA, can be compared with the direct transmissions preceding them by O(10) microseconds, to infer the surface power reflection coefficient $\cal{R}$. The first HiCal mission (HiCal-1, Jan. 2015) yielded a sample of 100 such pairs, resulting in estimates of $\cal{R}$ at highly-glancing angles (i.e., zenith angles approaching $90^\circ$), with measured reflectivity for those events which exceeded extant calculations. The HiCal-2 experiment, flying from Dec., 2016-Jan., 2017, provided an improvement by nearly two orders of magnitude in our event statistics, allowing a considerably more precise mapping of the reflectivity over a wider range of incidence angles. We find general agreement between the HiCal-2 reflectivity results and those obtained with the earlier HiCal-1 mission, as well as estimates from Solar reflections in the radio-frequency regime. In parallel, our calculations of expected reflectivity have matured; herein, we use a plane-wave expansion to estimate the reflectivity R from both a flat, smooth surface (and, in so doing, recover the Fresnel reflectivity equations) and also a curved surface. Multiplying our flat-smooth reflectivity by improved Earth curvature and surface roughness corrections now provides significantly better agreement between theory and the HiCal 2a/2b measurements.


title: The zero forcing polynomial of a graph

abstract: Zero forcing is an iterative graph coloring process, where given a set of initially colored vertices, a colored vertex with a single uncolored neighbor causes that neighbor to become colored. A zero forcing set is a set of initially colored vertices which causes the entire graph to eventually become colored. In this paper, we study the counting problem associated with zero forcing. We introduce the zero forcing polynomial of a graph $G$ of order $n$ as the polynomial $\mathcal{Z}(G;x)=\sum_{i=1}^n z(G;i) x^i$, where $z(G;i)$ is the number of zero forcing sets of $G$ of size $i$. We characterize the extremal coefficients of $\mathcal{Z}(G;x)$, derive closed form expressions for the zero forcing polynomials of several families of graphs, and explore various structural properties of $\mathcal{Z}(G;x)$, including multiplicativity, unimodality, and uniqueness.


title: Neutrino motion and spin oscillations in magnetic field and matter currents

abstract: A brief review of a neutrino oscillations effect in an external magnetic field and matter currents is given. An ultra-relativistic neutrino propagation in moving external media is investigated. We have found a new spin operator which commutes with the corresponding Hamiltonian, exact wave functions and energy spectrum are obtained.


title: Secure Distributed State Estimation of an LTI System over Time-Varying Networks and Analog Erasure Channels

abstract: We study the problem of collaboratively estimating the state of an LTI system monitored by a network of sensors, subject to the following important practical considerations: (i) certain sensors might be arbitrarily compromised by an adversary and (ii) the underlying communication graph governing the flow of information across sensors might be time-varying. We first analyze a scenario involving intermittent communication losses that preserve certain information flow patterns over bounded intervals of time. By equipping the sensors with adequate memory, we show that one can obtain a fully distributed, provably correct state estimation algorithm that accounts for arbitrary adversarial behavior, provided certain conditions are met by the network topology. We then argue that our approach can handle bounded communication delays as well. Next, we explore a case where each communication link stochastically drops packets based on an analog erasure channel model. For this setup, we propose state estimate update and information exchange rules, along with conditions on the network topology and packet drop probabilities, that guarantee mean-square stability despite arbitrary adversarial attacks.


title: Iterative Refinement of A Modified Lavrentiev Regularization Method for De-convolution of the Discrete Helmholtz Type Differential Filter

abstract: We propose and analyze an iterative refinement of a modified Lavrentiev regularization method for deconvolution of the discrete Helmholtz-type differential filter. The modification for the Lavrentiev regularization method exploits the properties of the Helmholtz filter, and we prove that the modification reduces the error bound between the original solution and the approximated solution. Furthermore, we derive an optimal stopping condition on the number of iterations necessary for the regularization. We provide numerical examples demonstrating the benefits of this iterative modified Lavrentiev regularization over a family of Tikhonov regularization methods.


title: Algebraic Hybridization and Static Condensation with Application to Scalable H(div) Preconditioning

abstract: We propose an unified algebraic approach for static condensation and hybridization, two popular techniques in finite element discretizations. The algebraic approach is supported by the construction of scalable solvers for problems involving H(div)-spaces discretized by conforming (Raviart-Thomas) elements of arbitrary order. We illustrate through numerical experiments the relative performance of the two (in some sense dual) techniques in comparison with a state-of-the-art parallel solver, ADS, available in the software hypre and MFEM. The superior performance of the hybridization technique is clearly demonstrated with increased benefit for higher order elements.


title: Spectral gaps of frustration-free spin systems with boundary

abstract: In quantum many-body systems, the existence of a spectral gap above the ground state has far-reaching consequences. In this paper, we discuss "finite-size" criteria for having a spectral gap in frustration-free spin systems and their applications. We extend a criterion that was originally developed for periodic systems by Knabe and Gosset-Mozgunov to systems with a boundary. Our finite-size criterion says that if the spectral gaps at linear system size $n$ exceed an explicit threshold of order $n^{-3/2}$, then the whole system is gapped. The criterion takes into account both "bulk gaps" and "edge gaps" of the finite system in a precise way. The $n^{-3/2}$ scaling is robust: it holds in 1D and 2D systems, on arbitrary lattices and with arbitrary finite-range interactions. One application of our results is to give a rigorous foundation to the folklore that 2D frustration-free models cannot host chiral edge modes (whose finite-size spectral gap would scale like $n^{-1}$).

categories: quant-ph cond-mat.stat-mech math-ph math.MP


title: Non-local Coulomb correlations in pure and electron-doped ${\mathrm{Sr}}_{2}{\mathrm{IrO}}_{4}$: spectral functions, Fermi surface and pseudogap-like spectral weight distributions from oriented cluster dynamical mean field theory

abstract: We address the role of non-local Coulomb correlations and short-range magnetic fluctuations in the high-temperature phase of Sr$_2$IrO$_4$ within state-of-the-art spectroscopic and first-principles theoretical methods. Introducing a novel cluster dynamical mean field scheme, we compute momentum-resolved spectral functions, which we find to be in excellent agreement with angle-resolved photoemission spectra. We show that while short-range antiferromagnetic fluctuations are crucial to account for the electronic properties of the material even in the high-temperature paramagnetic phase, long-range magnetic order is not a necessary ingredient of the insulating state. Upon doping, an exotic metallic state is generated, exhibiting cuprate-like pseudo-gap spectral properties, for which we propose a surprisingly simple theoretical mechanism.


title: Learning to Evade Static PE Machine Learning Malware Models via Reinforcement Learning

abstract: Machine learning is a popular approach to signatureless malware detection because it can generalize to never-before-seen malware families and polymorphic strains. This has resulted in its practical use for either primary detection engines or for supplementary heuristic detection by anti-malware vendors. Recent work in adversarial machine learning has shown that deep learning models are susceptible to gradient-based attacks, whereas non-differentiable models that report a score can be attacked by genetic algorithms that aim to systematically reduce the score. We propose a more general framework based on reinforcement learning (RL) for attacking static portable executable (PE) anti-malware engines. The general framework does not require a differentiable model nor does it require the engine to produce a score. Instead, an RL agent is equipped with a set of functionality-preserving operations that it may perform on the PE file. Through a series of games played against the anti-malware engine, it learns which sequences of operations are likely to result in evading the detector for any given malware sample. This enables completely black-box attacks against static PE anti-malware, and produces functional evasive malware samples as a direct result. We show in experiments that our method can attack a gradient-boosted machine learning model with evasion rates that are substantial and appear to be strongly dependent on the dataset. We demonstrate that attacks against this model appear to also evade components of publicly hosted antivirus engines. Adversarial training results are also presented: by retraining the model on evasive ransomware samples, a subsequent attack is 33% less effective. However, there are overfitting dangers when adversarial training, which we note. We release code to allow researchers to reproduce and improve this approach.


title: Thermal cycling effects on static and dynamic properties of a phase separated manganite

abstract: In this work we address the interplay between two phenomena which are signatures of the out-of-equilibrium state in phase separated manganites: irreversibility against thermal cycling and aging/rejuvenation process. The sample investigated is La0.5Ca0.5MnO3, a prototypical manganite exhibiting phase separation. Two regimes for isothermal relaxation were observed according to the temperature range: for T > 100 K, aging/rejuvenation effects are observed, while for T < 100 K an irreversible aging was found. Our results show that thermal cycles act as a tool to unveil the dynamical behavior of the phase separated state in manganites, revealing the close interplay between static and dynamic properties of phase separated manganites.


title: Intervals in the Hales-Jewett theorem

abstract: The Hales-Jewett theorem states that for any $m$ and $r$ there exists an $n$ such that any $r$-colouring of the elements of $[m]^n$ contains a monochromatic combinatorial line. We study the structure of the wildcard set $S \subseteq [n]$ which determines this monochromatic line, showing that when $r$ is odd there are $r$-colourings of $[3]^n$ where the wildcard set of a monochromatic line cannot be the union of fewer than $r$ intervals. This is tight, as for $n$ sufficiently large there are always monochromatic lines whose wildcard set is the union of at most $r$ intervals.


title: Artificial Kerr effect on the self-focusing of laser in a dissipative suspension of metallic nanoparticles

abstract: Self-focusing of laser beam propagating through a dissipative suspension of metallic nanoparticles is studied. Impact of imaginary part of nanoparticle polarizability on the optical force and consequently on the particles rearrangement in the presence of laser fields with an initial Gaussian profile is considered. It is shown that considering absorption of wave, leads to the creation of optical force along the wave travelling direction which can cause longitudinal distribution of nanoparticles density. Considering fifth order nonlinearity of wave amplitude that comes from simultaneous considering of normal Kerr effect produced by the inhomogeneity of refractive index resulted from ponderomotive force acting on electrons and artificial Kerr nonlinearity caused by the polarization optical force acting on particles, set of differential equations describing nonlinear steady-state evolution of laser beam is derived by using a non-paraxial method. Evolution of laser spot size for different frequencies is investigated and optimum frequency range for improving focusing property is determined. It is shown that the artificial Kerr effect causes concentration of particles near the propagation axis that can reduce threshold power for occurring self-focusing.


title: A Redesigned Benders Decomposition Approach for Large-Scale In-Transit Freight Consolidation Operations

abstract: The growth in online shopping and third party logistics has caused a revival of interest in finding optimal solutions to the large scale in-transit freight consolidation problem. Given the shipment date, size, origin, destination, and due dates of multiple shipments distributed over space and time, the problem requires determining when to consolidate some of these shipments into one shipment at an intermediate consolidation point so as to minimize shipping costs while satisfying the due date constraints. In this paper, we develop a mixed-integer programming formulation for a multi-period freight consolidation problem that involves multiple products, suppliers, and potential consolidation points. Benders decomposition is then used to replace a large number of integer freight-consolidation variables by a small number of continuous variables that reduces the size of the problem without impacting optimality. Our results show that Benders decomposition provides a significant scale-up in the performance of the solver. We demonstrate our approach using a large-scale case with more than 27.5 million variables and 9.2 million constraints.


title: The mechanism of propulsion of a model microswimmer in a viscoelastic fluid next to a solid boundary

abstract: In this paper we study swimming of a model organism, the so-called Taylor's swimming sheet, in a viscoelastic fluid close to a solid boundary. This situation comprises natural habitats of many swimming microorganisms, and while previous investigations have considered the effects of both swimming next to a boundary and swimming in a viscoelastic fluid, seldom have both effects been considered simultaneously. We re-visit the small wave amplitude result obtained by Elfring and Lauga (Gwynn J. Elfring and Eric Lauga, in Saverio E. Spagnolie, editor, Complex Fluids in Biological Systems, Springer New York, New York, NY, 2015), and give a mechanistic explanation to the decoupling of the effects of viscoelasticity, which tend to slow the sheet, and the presence of the boundary, which tends to speed up the sheet. We also develop a numerical spectral method capable of finding the swimming speed of a waving sheet with an arbitrary amplitude and waveform. We use it to show that the decoupling mentioned above does not hold at finite wave amplitudes and that for some parameters the presence of a boundary can cause the viscoelastic effects to increase the swimming speed of microorganisms.


title: Asymptotics of quantum weighted Hurwitz numbers

abstract: This work concerns both the semiclassical and zero temperature asymptotics of quantum weighted double Hurwitz numbers. The partition function for quantum weighted double Hurwitz numbers can be interpreted in terms of the energy distri- bution of a quantum Bose gas with vanishing fugacity. We compute the leading semi- classical term of the partition function for three versions of the quantum weighted Hurwitz numbers, as well as lower order semiclassical corrections. The classical limit $\hbar \ra 0$ is shown to reproduce the simple single and double Hurwitz numbers studied by Pandharipande and Okounkov [20,22]. The KP-Toda $\tau$-function that serves as generating function for the quantum Hurwitz numbers is shown to have the $\tau$-function of [20,22] as its leading term in the classical limit, and, with suitable scaling, the same holds for the partition function, the weights and expectations of Hurwitz numbers. We also compute the zero temperature limit $T \ra 0$ of the partition function and quantum weighted Hurwitz numbers. The KP or Toda $\tau$-function serving as generating function for the quantum Hurwitz numbers are shown to give the one for Belyi curves in the zero temperature limit and, with suitable scaling, the same holds true for the partition function, the weights and the expectations of Hurwitz numbers.


title: Power spectrum of rare events in two-dimensional BTW model, violation of $1/f$ noise hypothesis

abstract: One of the primitive aims of the two-dimensional BTW model had been to explain the $1/f^{\alpha}$ noise which is widely seen in the natural systems. In this paper we study some time signals, namely the activity inside an avalanche ($x(t)$), the avalanches sizes ($s(T)$) and the rare events waiting time (REWT) ($\tau(n)$ as a new type of noise). The latter is expected to be important in predicting the period and also the vastness of the upcoming large scale events in a sequence of self-organized natural events. Especially we report some exponential anti-correlation behaviors for $s(T)$ and $\tau(n)$ which are finite-size effects. Two characteristic time scales $\delta T_{s}$ and $\delta T_{\tau}$ emerge in our analysis. It is proposed that the power spectrum of $s(T)$ and $\tau(n)$ behave like $\left( b_{s,\tau}(L)^2+\omega^2\right)^{-\frac{1}{2}}$, in which $b_{s}$ and $b_{\tau}$ are some $L$-dependent parameters and $\omega$ is the angular frequency. The $1/f$ noise is therefore obtained in the limit $\omega\gg b_{s,\tau}$. $b_{s}$ and $b_{\tau}$ decrease also in a power-law fashion with the system size $L$, which signals the fact that in the thermodynamic limit the power spectrum tends to the Dirac delta function.


title: Supersaliency: A Novel Pipeline for Predicting Smooth Pursuit-Based Attention Improves Generalizability of Video Saliency

abstract: Predicting attention is a popular topic at the intersection of human and computer vision. However, even though most of the available video saliency data sets and models claim to target human observers' fixations, they fail to differentiate them from smooth pursuit (SP), a major eye movement type that is unique to perception of dynamic scenes. In this work, we highlight the importance of SP and its prediction (which we call supersaliency, due to greater selectivity compared to fixations), and aim to make its distinction from fixations explicit for computational models. To this end, we (i) use algorithmic and manual annotations of SP and fixations for two well-established video saliency data sets, (ii) train Slicing Convolutional Neural Networks for saliency prediction on either fixation- or SP-salient locations, and (iii) evaluate our and 26 publicly available dynamic saliency models on three data sets against traditional saliency and supersaliency ground truth. Overall, our models outperform the state of the art in both the new supersaliency and the traditional saliency problem settings, for which literature models are optimized. Importantly, on two independent data sets, our supersaliency model shows greater generalization ability and outperforms all other models, even for fixation prediction.


title: Deflecting Adversarial Attacks with Pixel Deflection

abstract: CNNs are poised to become integral parts of many critical systems. Despite their robustness to natural variations, image pixel values can be manipulated, via small, carefully crafted, imperceptible perturbations, to cause a model to misclassify images. We present an algorithm to process an image so that classification accuracy is significantly preserved in the presence of such adversarial manipulations. Image classifiers tend to be robust to natural noise, and adversarial attacks tend to be agnostic to object location. These observations motivate our strategy, which leverages model robustness to defend against adversarial perturbations by forcing the image to match natural image statistics. Our algorithm locally corrupts the image by redistributing pixel values via a process we term pixel deflection. A subsequent wavelet-based denoising operation softens this corruption, as well as some of the adversarial changes. We demonstrate experimentally that the combination of these techniques enables the effective recovery of the true class, against a variety of robust attacks. Our results compare favorably with current state-of-the-art defenses, without requiring retraining or modifying the CNN.


title: Majorana bound states and subgap states in three-terminal topological superconducting nanowire-quantum dot hybrid devices

abstract: Three-terminal topological superconducting nanowire (TSNW)-quantum dot (QD) hybrid junction devices are studied. The energy spectra and the wave functions of the subgap states are calculated as a function of the superconducting phase differences between TSNWs and as a function of the QD level energy based on the Bogoliubov-de Gennes tight-binding Hamiltonians. It is shown that when the QD level is located near or inside the superconducting gap, there can exist eight subgap states. Among them, four low energy (two positive and two negative) subgap states are essentially formed by linear combinations of the six Majorana bound states (MBSs) located at the ends of the three TSNWs. The remaining four high energy subgap states are mainly built from linear combinations of the QD state and the three MBSs of the TSNWs adjacent to the QD. When there is no QD level near or inside the superconducting gap, only six subgap states built from linear combinations of the six MBSs of the three TSNWs can be present in the system. It is also shown that there exists a unique point in the parameter space of the superconducting phase differences between TSNWs, at which the energies of the six low energy subgap states move close to each other towards nearly zero energies. Simple but general effective model Hamiltonians for the three-terminal TSNW-QD hybrid devices have also been developed. Based on the effective model Hamiltonians, the subgap states of the three-terminal TSNW-QD hybrid devices in the limit of the three infinitely long TSNWs are studied. The results of the calculations and the effective model Hamiltonians could be used as a starting point to construct and investigate the braiding schemes of MBSs in TSNW junction devices.


title: Automatically Extracting Web API Specifications from HTML Documentation

abstract: Web API specifications are machine-readable descriptions of APIs. These specifications, in combination with related tooling, simplify and support the consumption of APIs. However, despite the increased distribution of web APIs, specifications are rare and their creation and maintenance heavily relies on manual efforts by third parties. In this paper, we propose an automatic approach and an associated tool called D2Spec for extracting specifications from web API documentation pages. Given a seed online documentation page on an API, D2Spec first crawls all documentation pages on the API, and then uses a set of machine learning techniques to extract the base URL, path templates, and HTTP methods, which collectively describe the endpoints of an API. We evaluated whether D2Spec can accurately extract endpoints from documentation on 120 web APIs. The results showed that D2Spec achieved a precision of 87.5% in identifying base URLs, a precision of 81.3% and a recall of 80.6% in generating path templates, and a precision of 84.4% and a recall of 76.2% in extracting HTTP methods. In addition, we found that D2Spec was useful when applied to APIs with pre-existing API specifications: D2Spec revealed many inconsistencies between web API documentation and their corresponding publicly available specifications. Thus, D2Spec can be used by web API providers to keep documentation and specifications in synchronization.


title: Methodological variations in lagged regression for detecting physiologic drug effects in EHR data

abstract: We studied how lagged linear regression can be used to detect the physiologic effects of drugs from data in the electronic health record (EHR). We systematically examined the effect of methodological variations ((i) time series construction, (ii) temporal parameterization, (iii) intra-subject normalization, (iv) differencing (lagged rates of change achieved by taking differences between consecutive measurements), (v) explanatory variables, and (vi) regression models) on performance of lagged linear methods in this context. We generated two gold standards (one knowledge-base derived, one expert-curated) for expected pairwise relationships between 7 drugs and 4 labs, and evaluated how the 64 unique combinations of methodological perturbations reproduce gold standards. Our 28 cohorts included patients in Columbia University Medical Center/NewYork-Presbyterian Hospital clinical database. The most accurate methods achieved AUROC of 0.794 for knowledge-base derived gold standard (95%CI [0.741, 0.847]) and 0.705 for expert-curated gold standard (95% CI [0.629, 0.781]). We observed a 0.633 mean AUROC (95%CI [0.610, 0.657], expert-curated gold standard) across all methods that re-parameterize time according to sequence and use either a joint autoregressive model with differencing or an independent lag model without differencing. The complement of this set of methods achieved a mean AUROC close to 0.5, indicating the importance of these choices. We conclude that time- series analysis of EHR data will likely rely on some of the beneficial pre-processing and modeling methodologies identified, and will certainly benefit from continued careful analysis of methodological perturbations. This study found that methodological variations, such as pre-processing and representations, significantly affect results, exposing the importance of evaluating these components when comparing machine-learning methods.


title: Recasting Gradient-Based Meta-Learning as Hierarchical Bayes

abstract: Meta-learning allows an intelligent agent to leverage prior learning episodes as a basis for quickly improving performance on a novel task. Bayesian hierarchical modeling provides a theoretical framework for formalizing meta-learning as inference for a set of parameters that are shared across tasks. Here, we reformulate the model-agnostic meta-learning algorithm (MAML) of Finn et al. (2017) as a method for probabilistic inference in a hierarchical Bayesian model. In contrast to prior methods for meta-learning via hierarchical Bayes, MAML is naturally applicable to complex function approximators through its use of a scalable gradient descent procedure for posterior inference. Furthermore, the identification of MAML as hierarchical Bayes provides a way to understand the algorithm's operation as a meta-learning procedure, as well as an opportunity to make use of computational strategies for efficient inference. We use this opportunity to propose an improvement to the MAML algorithm that makes use of techniques from approximate inference and curvature estimation.


title: Talagrand inequality at second order and application to Boolean analysis

abstract: This note is concerned with an extension, at second order, of an inequality on the discrete cube $C_n=\{-1,1\}$ (equipped with the uniform measure) due to Talagrand (\cite{TalL1L2}). As an application, the main result of this note is a Theorem in the spirit of a famous result from Kahn, Kalai and Linial (cf. \cite{KKL}) concerning the influence of Boolean functions. The notion of the influence of a couple of coordinates $(i,j)\in\{1,\ldots,n\}^2$ is introduced in section 2 and the following alternative is obtained : for any Boolean function $f\,:\, C_n\to \{0,1\}$, either there exists a coordinate with influence at least of order $(1/n)^{1/(1+\eta)}$, with $\, 0<\eta<1$ (independent of $f$ and $n$) or there exists a couple of coordinates $(i,j)\in\{1,\ldots,n\}^2$ with $i\neq j$, with influence at least of order $(\log n/n)^2$. In section 4, it is shown that this extension of Talagrand inequality can also be obtained, with minor modifications, for the standard Gaussian measure $\gamma_n$ on $\mathbb{R}^n$ ; the obtained inequality can be of independent interest. The arguments rely on interpolation methods by semigroup together with hypercontractive estimates. At the end of the article, some related open questions are presented.


title: Two-dimensional Mott variable-range hopping transport in a disordered MoS$_2$ nanoflake

abstract: The transport characteristics of a disordered MoS$_2$ nanoflake in the insulator regime are studied by electrical and magnetotransport measurements. The layered MoS$_2$ nanoflake is exfoliated from a bulk MoS$_2$ crystal and the conductance $G$ and magnetoresistance are measured in a four-probe setup over a wide range of temperatures. At high temperatures, we observe that $\log_{10}G$ exhibits a $-T^{-1}$ temperature dependence and the transport in the nanoflake dominantly arises from thermal activation. At low temperatures, where the transport in the nanoflake dominantly takes place via variable-range hopping (VRH) processes, we observe that $\log_{10}G$ exhibits a $-T^{-1/3}$ temperature dependence, an evidence for the two-dimensional (2D) Mott VRH transport. The measured low-field magnetoresistance of the nanoflake in the insulator regime exhibits a quadratic magnetic field dependence $\sim \alpha B^2$ with $\alpha\sim T^{-1}$, fully consistent with the 2D Mott VRH transport in the nanoflake.


title: Transfer of A-infinity structures to projective resolutions

abstract: We show that an A-infinity algebra structure can be transferred to a projective resolution of the complex underlying any A-infinity algebra. Under certain connectedness assumptions, this transferred structure is unique up to homotopy. In contrast to the classical results on transfer of A-infinity structures along homotopy equivalences, our result is of interest when the ground ring is not a field. We prove an analog for A-infinity module structures, and both transfer results preserve strict units.


title: Limit theorems for the least common multiple of a random set of integers

abstract: Let $L_{n}$ be the least common multiple of a random set of integers obtained from $\{1,\ldots,n\}$ by retaining each element with probability $\theta\in (0,1)$ independently of the others. We prove that the process $(\log L_{\lfloor nt\rfloor})_{t\in [0,1]}$, after centering and normalization, converges weakly to a certain Gaussian process that is not Brownian motion. Further results include a strong law of large numbers for $\log L_{n}$ as well as Poisson limit theorems in regimes when $\theta$ depends on $n$ in an appropriate way.


title: Solvable Leibniz algebras with quasi-filiform Lie algebras of maximum length nilradicals

abstract: In this paper solvable Leibniz algebras whose nilradical is quasi-filiform Lie algebra of maximum length, are classified. The rigidity of such Leibniz algebras with two-dimensional complemented space to nilradical is proved.


title: Searching for Primordial Magnetic Fields with CMB B-modes

abstract: Was the primordial universe magnetized? The answer to this question would help explain the origin of micro-Gauss strength magnetic fields observed in galaxies. It is also of fundamental importance in developing a complete theory of the early universe. While there can be other signatures of cosmological magnetic fields, a signature in the cosmic microwave background (CMB) would prove their primordial origin. The B-mode polarization of CMB is particularly promising in this regard because there are relatively few other sources of B-modes, and because the vortical modes sourced by the primordial magnetic field (PMF) survive diffusion damping up to a small fraction of the Silk length. At present, the Planck temperature and polarization spectra combined with the B-mode spectrum measured by the South Pole Telescope (SPT) constrain the PMF strength to be no more than $\sim 1$ nano-Gauss (nG). Because of the quartic scaling of the CMB anisotropy spectra with the PMF strength, this bound will not change by much even with the significantly better measurements of the B-mode spectrum by the Stage III and Stage IV CMB experiments. On the other hand, tightening the bound well below the $1$ nG threshold would rule out the purely primordial origin (requiring no dynamo action) of galactic fields. Considering Faraday rotation, which converts some of the E-modes into B-modes and scales linearly with the field strength, will help to achieve this goal. As we demonstrate, the upcoming experiments, such as SPT-3G and the Simons Observatory, will be sensitive to fields of $\sim 0.5$ nG strength thanks to the mode-coupling signature induced by Faraday rotation. A future Stage IV ground based experiment or a Space Probe will be capable of probing fields below $0.1$ nG, and would detect a scale-invariant PMF of $0.2$ nG strength without de-lensing or subtracting the galactic rotation measure.


title: Iterative map-making with two-level preconditioning for polarized Cosmic Microwave Background data sets

abstract: An estimation of the sky signal from streams of Time Ordered Data (TOD) acquired by Cosmic Microwave Background (\cmb) experiments is one of the most important steps in the context of \cmb data analysis referred to as the map-making problem. The continuously growing \cmb data sets render the \cmb map-making problem more challenging in terms of computational cost and memory in particular in the context of ground based experiments. In this context, we study a novel class of the Preconditioned Conjugate Gradient (PCG) solvers which invoke two-level preconditioners. We compare them against PCG solvers commonly used in the map-making context considering their precision and time-to-solution. We compare these new methods on realistic, simulated data sets reflecting the characteristics of current and forthcoming \cmb ground-based experiment. We develop an embarrassingly parallel implementation of the approach where each processor performs a sequential map-making for a subset of the TOD. We find that considering the map level residuals the new class of solvers permits achieving tolerance of up to 3 orders of magnitude better than the standard approach, where the residual level often saturates before convergence is reached. This corresponds to an important improvement in the precision of recovered power spectra in particular on the largest angular scales. The new method also typically requires fewer iterations to reach a required precision and thus shorter runtimes for a single map-making solution. However, the construction of an appropriate two-level preconditioner can be as costly as a single standard map-making run. Nevertheless, if the same problem needs to be solved multiple times, e.g., as in Monte Carlo simulations, this cost has to be incurred only once, and the method should be competitive not only as far as its precision but also its performance is concerned.


title: Simulation for L3 Volumetric Attack Detection

abstract: The detection of a volumetric attack involves collecting statistics on the network traffic, and identifying suspicious activities. We assume that available statistical information includes the number of packets and the number of bytes passed per flow. We apply methods of machine learning to detect malicious traffic. A prototype project is implemented as a module for the Floodlight controller. The prototype was tested on the Mininet simulation platform. The simulated topology includes a number of edge switches, a connected graph of core switches, and a number of server and user hosts. The server hosts run simple web servers. The user hosts simulate web clients. The controller employs Dijkstra's algorithm to find the best flow in the graph. The controller periodically polls the edge switches and provides current and historical statistics on each active flow. The streaming analytics evaluates the traffic volume and detects volumetric attacks.


title: Calder\'on's reproducing formulas for the Weinstein $L^2$- multiplier operators

abstract: The aim of this work is the study of the Weinstein $L^2$- multiplier operators on $\mathbb{R}^{d+1}_+$ and we give for them Calder\'on's reproducing formulas and best approximation using the theory of Weinstein transform and reproducing kernels.


title: Precision mass measurements on neutron-rich rare-earth isotopes at JYFLTRAP - reduced neutron pairing and implications for the $r$-process calculations

abstract: The rare-earth peak in the $r$-process abundance pattern depends sensitively on both the astrophysical conditions and subtle changes in nuclear structure in the region. This work takes an important step elucidating the nuclear structure and reducing the uncertainties in $r$-process calculations via precise atomic mass measurements at the JYFLTRAP double Penning trap. $^{158}$Nd, $^{160}$Pm, $^{162}$Sm, and $^{164-166}$Gd have been measured for the first time and the precisions for $^{156}$Nd, $^{158}$Pm, $^{162,163}$Eu, $^{163}$Gd, and $^{164}$Tb have been improved considerably. Nuclear structure has been probed via two-neutron separation energies $S_{2n}$ and neutron pairing energy metrics $D_n$. The data do not support the existence of a subshell closure at $N=100$. Neutron pairing has been found to be weaker than predicted by theoretical mass models. The impact on the calculated $r$-process abundances has been studied. Substantial changes resulting in a smoother abundance distribution and a better agreement with the solar $r$-process abundances are observed.


title: Enhanced moments of Eu in single crystals of the metallic helical antiferromagnet EuCo{2-y}As2

abstract: The compound EuCo{2-y}As2 with the tetragonal ThCr2Si2 structure is known to contain Eu{+2} ions with spin S = 7/2 that order below a temperature TN = 47 K into an antiferromagnetic (AFM) proper helical structure with the ordered moments aligned in the tetragonal ab plane, perpendicular to the helix axis along the c axis, with no contribution from the Co atoms. Here we carry out a detailed investigation of the properties of single crystals. Enhanced ordered and effective moments of the Eu spins are found in most of our crystals. Electronic structure calculations indicate that the enhanced moments arise from polarization of the d bands, as occurs in ferromagnetic Gd metal. Electrical resistivity measurements indicate metallic behavior. The low-field in-plane magnetic susceptibilities chi{ab}(T < TN) for several crystals are reported that are fitted well by unified molecular field theory (MFT), and the Eu-Eu exchange interactions Jij are extracted from the fits. High-field magnetization M data for magnetic fields H||ab reveal what appears to be a first-order spin-flop transition followed at higher field by a second-order metamagnetic transition of unknown origin, and then by another second-order transition to the paramagnetic (PM) state. For H||c, the magnetization shows only a second-order transition from the canted AFM to the PM state, as expected. The critical fields for the AFM to PM transition are in approximate agreement with the predictions of MFT. Heat capacity Cp measurements in zero and high H are reported. Phase diagrams for H||c and H||ab versus T are constructed from the high-field M(H,T) and Cp(H,T) measurements. The magnetic part Cmag(T, H = 0) of Cp(T, H = 0) is extracted and is fitted rather well below TN by MFT, although dynamic short-range AFM order is apparent in Cmag(T) up to about 70 K, where the molar entropy attains its high-T limit of R ln8.


title: Exact edge, bulk and bound states of finite topological systems

abstract: Finite topologically non-trivial systems are often characterised by the presence of bound states at their physical edges. These topological edge modes can be distinguished from usual Shockley waves energetically, as their energies remain finite and in-gap. On a clean 1D or reducible 2D model, in either the commensurate or semi-infinite case, the edge modes can be obtained analytically, as shown in [PRL 71, 3697 (1993)] and [PRA 89, 023619 (2014)]. We put forward a method for obtaining the spectrum and wave functions of topological edge modes for arbitrary finite lattices, including the incommensurate case. A small number of parameters are easily determined numerically, with the form of the eigenstates remaining fully analytical. We also obtain the bulk modes in the finite system analytically and their eigenenergies, which lie within the infinite-size limit continuum. Our method is general and can be easily applied to obtain the properties of non-topological models and/or extended to include impurities. As an example, we consider the case of an impurity located next to one edge of a 1D system, equivalent to a softened boundary in a separable 2D model. We show that a localised impurity can have a drastic effect on the edge modes of the system. Using the periodic Harper and Hofstadter models to illustrate our method, we find that, on increasing the impurity strength, edge states can enter or exit the continuum, and a trivial Shockley state bound to the impurity may appear. The fate of the topological edge modes in the presence of impurities can be addressed by quenching the impurity strength. We find that at certain critical impurity strengths, the transition probability for a particle initially prepared in an edge mode to decay into the bulk exhibits discontinuities that mark the entry and exit points of edge modes from and into the bulk spectrum.


title: $r$-Process Nucleosynthesis from Three-dimensional Jet-driven Core-Collapse Supernovae with Magnetic Misalignments

abstract: We investigate $r$-process nucleosynthesis in three-dimensional general relativistic magnetohydrodynamic simulations of jet-driven supernovae resulting from rapidly rotating, strongly magnetized core-collapse. We explore the effect of misaligning the pre-collapse magnetic field with respect to the rotation axis by performing four simulations: one aligned model and models with 15, 30, and 45 degree misalignments. The simulations we present employ a microphysical finite-temperature equation of state and a leakage scheme that captures the overall energetics and lepton number exchange due to post-bounce neutrino emission and absorption. We track the thermodynamic properties of the ejected material with Lagrangian tracer particles and analyse its composition with the nuclear reaction network SkyNet. By using different neutrino luminosities in post-processing the tracer data with SkyNet, we constrain the impact of uncertainties in neutrino luminosities. We find that, for the aligned model considered here, the use of an approximate leakage scheme results in neutrino luminosity uncertainties corresponding to a factor of 100-1000 uncertainty in the abundance of third peak $r$-process elements. Our results show that for misalignments of 30 degrees or less, $r$-process elements are robustly produced as long as neutrino luminosities are reasonably low ($\lesssim 5 \times 10^{52}$ erg s$^{-1}$). For a more extreme misalignment of 45 degrees, we find the production of $r$-process elements beyond the second peak significantly reduced. We conclude that robust $r$-process nucleosynthesis in magnetorotational supernovae requires a progenitor stellar core with a large poloidal magnetic field component that is at least moderately (within $\sim 30$ degrees) aligned with the rotation axis.


title: Decay of solutions of the wave equation in expanding cosmological spacetimes

abstract: We study the decay of solutions of the wave equation in some expanding cosmological spacetimes, namely flat Friedmann-Lema\^itre-Robertson-Walker (FLRW) models and the cosmological region of the Reissner-Nordstr\"om-de Sitter (RNdS) solution. By introducing a partial energy and using an iteration scheme, we find that, for initial data with finite higher order energies, the decay rate of the time derivative is faster than previously existing estimates. For models undergoing accelerated expansion, our decay rate appears to be (almost) sharp.


title: Improving Weak Lensing Mass Map Reconstructions using Gaussian and Sparsity Priors: Application to DES SV

abstract: Mapping the underlying density field, including non-visible dark matter, using weak gravitational lensing measurements is now a standard tool in cosmology. Due to its importance to the science results of current and upcoming surveys, the quality of the convergence reconstruction methods should be well understood. We compare three methods: Kaiser-Squires (KS), Wiener filter, and GLIMPSE. KS is a direct inversion, not accounting for survey masks or noise. The Wiener filter is well-motivated for Gaussian density fields in a Bayesian framework. GLIMPSE uses sparsity, aiming to reconstruct non-linearities in the density field. We compare these methods with several tests using public Dark Energy Survey (DES) Science Verification (SV) data and realistic DES simulations. The Wiener filter and GLIMPSE offer substantial improvements over smoothed KS with a range of metrics. Both the Wiener filter and GLIMPSE convergence reconstructions show a 12 per cent improvement in Pearson correlation with the underlying truth from simulations. To compare the mapping methods' abilities to find mass peaks, we measure the difference between peak counts from simulated {\Lambda}CDM shear catalogues and catalogues with no mass fluctuations (a standard data vector when inferring cosmology from peak statistics); the maximum signal-to-noise of these peak statistics is increased by a factor of 3.5 for the Wiener filter and 9 for GLIMPSE. With simulations we measure the reconstruction of the harmonic phases; the phase residuals' concentration is improved 17 per cent by GLIMPSE and 18 per cent by the Wiener filter. The correlation between reconstructions from data and foreground redMaPPer clusters is increased 18 per cent by the Wiener filter and 32 per cent by GLIMPSE.


title: Expected Precision of Europa Clipper Gravity Measurements

abstract: The primary gravity science objective of NASA's Clipper mission to Europa is to confirm the presence or absence of a global subsurface ocean beneath Europa's Icy crust. Gravity field measurements obtained with a radio science investigation can reveal much about Europa's interior structure. Here, we conduct extensive simulations of the radio science measurements with the anticipated spacecraft trajectory and attitude (17F12v2) and assets on the spacecraft and the ground, including antenna orientations and beam patterns, transmitter characteristics, and receiver noise figures. In addition to two-way Doppler measurements, we also include radar altimeter crossover range measurements. We concentrate on +/-2 hour intervals centered on the closest approach of each of the 46 flybys. Our covariance analyses reveal the precision with which the tidal Love number k2, second-degree gravity coefficients C20 and C22, and higher-order gravity coefficients can be determined. The results depend on the Deep Space Network (DSN) assets that are deployed to track the spacecraft. We find that some DSN allocations are sufficient to conclusively confirm the presence or absence of a global ocean. Given adequate crossover range performance, it is also possible to evaluate whether the ice shell is hydrostatic.

categories: astro-ph.EP physics.data-an physics.geo-ph physics.space-ph


title: Heavy Neutral Leptons at FASER

abstract: We study the prospects for discovering heavy neutral leptons at ForwArd Search ExpeRiment (FASER), the newly proposed detector at the LHC. Previous studies showed that a relatively small detector with ~10 m length and ~1 m cross sectional area can probe large unconstrained parts of parameter space for dark photons and dark Higgs bosons. In this work we show that FASER will also be sensitive to heavy neutral leptons that have mixing angles with the active neutrinos that are up to an order of magnitude lower than current bounds. In particular, this is true for heavy neutral leptons produced dominantly in $B$-meson decays, in which case FASER's discovery potential is comparable to the proposed SHiP detector. We also illustrate how the search for heavy neutral leptons at FASER will be complementary to ongoing searches in high-$p_T$ experiments at the LHC and can shed light on the nature of dark matter and the process of baryogenesis in the early Universe.


title: A deeper look at the GD1 stream: density variations and wiggles

abstract: Using deep photometric data from CFHT/Megacam, we study the morphology and density of the GD-1 stream, one of the longest and coldest stellar streams in the Milky Way. Our deep data recovers the lower main sequence of the stream with unprecedented quality, clearly separating it from Milky Way foreground and background stars. An analysis of the distance to different parts of the stream shows that GD-1 lies at a heliocentric distance between 8 and 10 kpc, with only a shallow gradient across 45 deg on the sky. Matched filter maps of the stream density show clear density variations, such as deviations from a single orbital track and tentative evidence for stream fanning. We also detect a clear under-density in the middle of the stream track at $\varphi_{1}$=-45 deg surrounded by overdense stream segments on either side. This location is a promising candidate for the elusive missing progenitor of the GD-1 stream. We conclude that the GD-1 stream has clearly been disturbed by interactions with the Milky Way disk or other sub-halos.


title: Out-of-time-ordered measurements as a probe of quantum dynamics

abstract: Probing the out-of-equilibrium dynamics of quantum matter has gained renewed interest owing to immense experimental progress in artifcial quantum systems. Dynamical quantum measures such as the growth of entanglement entropy (EE) and out-of-time ordered correlators (OTOCs) have been shown, theoretically, to provide great insight by exposing subtle quantum features invisible to traditional measures such as mass transport. However, measuring them in experiments requires either identical copies of the system, an ancilla qubit coupled to the whole system, or many measurements on a single copy, thereby making scalability extremely complex and hence, severely limiting their potential. Here, we introduce an alternate quantity $-$ the out-of-time-ordered measurement (OTOM) $-$ which involves measuring a single observable on a single copy of the system, while retaining the distinctive features of the OTOCs. We show, theoretically, that OTOMs are closely related to OTOCs in a doubled system with the same quantum statistical properties as the original system. Using exact diagonalization, we numerically simulate classical mass transport, as well as quantum dynamics through computations of the OTOC, the OTOM, and the EE in quantum spin chain models in various interesting regimes (including chaotic and many-body localized systems). Our results demonstrate that an OTOM can successfully reveal subtle aspects of quantum dynamics hidden to classical measures, and crucially, provide experimental access to them.


title: Ictal and Post Ictal Impaired Consciousness due to Enhanced Mutual Information in Temporal Lobe Epilepsy

abstract: Seizure and synchronization are related to each other in complex manner. Altered synchrony has been implicated in loss of consciousness during partial seizures. However, the mechanism of altered consciousness following termination of seizures has not been studied well. In this work we used bivariate mutual information as a measure of synchronization to understand the neural correlate of altered consciousness during and after termination of mesial temporal lobe onset seizures. First, we have compared discrete bivariate mutual information (MI) measure with amplitude correlation (AC), phase synchronization (PS), nonlinear correlation and coherence, and established MI as a robust measure of synchronization. Next, we have extended MI to more than two signals by principal component method. The extended MI was applied on intracranial electroencephalogram (iEEG) before, during and after 23 temporal lobe seizures recorded from 11 patients. The analyses were carried out in delta, theta, alpha, beta and gamma bands. In 77% of the complex partial seizures MI was higher towards the seizure offset than in the first half of the seizure in the seizure onset zone (SOZ) channels in beta and gamma bands, whereas MI remained higher in the beginning or in the middle of the seizure than towards the offset across the least involved channels in the same bands. Synchronization seems built up outside the SOZ, gradually spread and culminated in SOZ and remained high beyond offset leading to impaired consciousness in 82% of the complex partial temporal lobe seizures. Consciousness impairment was scored according to a method previously applied to assess the same in patients with temporal lobe epilepsy during seizure.


title: Graph-Theoretic Framework for Unified Analysis of Observability and Data Injection Attacks in the Smart Grid

abstract: In this paper, a novel graph-theoretic framework is proposed to generalize the analysis of a broad set of security attacks, including observability and data injection attacks, that target the state estimator of a smart grid. First, the notion of observability attacks is defined based on a proposed graph-theoretic construct. In this respect, a structured approach is proposed to characterize critical sets, whose removal renders the system unobservable. It is then shown that, for the system to be observable, these critical sets must be part of a maximum matching over a proposed bipartite graph. In addition, it is shown that stealthy data injection attacks (SDIAs) constitute a special case of these observability attacks. Then, various attack strategies and defense policies, for observability and data injection attacks, are shown to be amenable to analysis using the introduced graph-theoretic framework. The proposed framework is then shown to provide a unified basis for analysis of four key security problems (among others), pertaining to the characterization of: 1) The sparsest SDIA; 2) the sparsest SDIA including a certain measurement; 3) a set of measurements which must be defended to thwart all potential SDIAs; and 4) the set of measurements, which when protected, can thwart any SDIA whose cardinality is below a certain threshold. A case study using the IEEE 14-bus system with a set of 17 measurements is used to support the theoretical findings.


title: Exact Hydrodynamic Description of Active Lattice Gases

abstract: We introduce a class of lattice gas models of active matter systems whose hydrodynamic description can be derived exactly. We illustrate our approach by considering two systems exhibiting two of the most studied collective behaviours in active matter: the motility-induced phase separation and the transition to collective motion. In both cases, we derive coupled partial differential equation describing the dynamics of the local density and polarization fields and show how they quantitatively predict the emerging properties of the macroscopic lattice gases.


title: The totally nonnegative part of G/P is a ball

abstract: We show that the totally nonnegative part of a partial flag variety (in the sense of Lusztig) is homeomorphic to a closed ball.


title: Common Occurrence of Explosive Hydrogen Burning in Type II Supernovae

abstract: We report Mo isotopic data for 16 15N-rich presolar SiC grains of type AB (14N/15N<solar, AB1) and their correlated Sr and Ba isotope ratios when available. Eight of the 16 AB1 grains show s-process Mo isotopic compositions, together with s-process Ba and/or Sr isotopic compositions. We found that a higher percentage of AB1 grains show anomalous isotopic compositions than that of AB2 grains (14N/15N>solar), thus providing further support to the division of the two AB subgroups recently proposed by Liu et al. (2017a), who showed that AB1 grains most likely originated from Type II supernovae (SNe) with explosive H burning. Comparison of the Sr, Mo, and Ba isotopic compositions of the AB1 grains with SN model predictions indicates that the s-process isotopic compositions of AB1 grains resulted from neutron-capture processes occurring during the progenitor massive stars' pre-SN evolution rather than from an explosive neutron-capture process. In addition, the observations of (1) explosive H burning occurring in the C-rich regions of the progenitor SNe of SN X grains as suggested by the isotopic compositions of X grains, and (2) explosive H burning occurring both at the bottom of the He/C zone and at the top of the He/N zone as suggested by model simulations, imply that explosive H burning is a common phenomenon in outer SN zones.


title: Stability under scalar perturbations and quasinormal modes of 4D Einstein-Born-Infeld dilaton spacetime: Exact spectrum

abstract: We study the stability under scalar perturbations, and we compute the quasinormal modes of the Einstein-Born-Infeld dilaton spacetime in 1+3 dimensions. Solving the full radial equation in terms of hypergeometric functions, we provide an exact analytical expression for the spectrum. We find that the frequencies are purely imaginary, and we confirm our results by computing them numerically. Although the scalar field that perturbs the black hole is electrically neutral, an instability similar to that seen in charged scalar perturbations of the Reissner-Nordstr\"om black hole is observed.


title: Canonical diffusions on the pattern spaces of aperiodic Delone sets

abstract: We consider pattern spaces of aperiodic and repetitive Delone sets of finite local complexity. These spaces are compact metric spaces and constitute a special class of foliated spaces. We define new Sobolev spaces with respect to the unique invariant measure and prove the existence of the unitary Schr\"odinger semigroup, which in physics terms describe the evolution of phasons. We define and study natural leafwise diffusion processes on these pattern spaces. These processes have Feller, but lack strong Feller and hypercontractivity properties, and heat kernels do not exist. The associated Dirichlet forms are regular, strongly local, irreducible and recurrent, but not strictly local. For harmonic functions we prove new Liouville and Helmholtz-Hodge type theorems.