title: An Intra-Cluster Model with Diffuse Scattering for mmWave Communications: RT-ICM

abstract: In millimeter-wave (mmWave) channels, to overcome the high path loss, beamforming is required. Hence, the spatial representation of the channel is essential. Further, for accurate beam alignment and minimizing the outages, inter-beam interferences, etc., cluster-level spatial modeling is also necessary. Since, statistical channel models fail to reproduce the intra-cluster parameters due to the site-specific nature of the mmWave channel, in this paper, we propose a ray tracing intra-cluster model (RT-ICM) for mmWave channels. The model considers only the first-order reflection; thereby reducing the computation load while capturing most of the energy in a large number of important cases. The model accounts for diffuse scattering as it contributes significantly to the received power. Finally, since the clusters are spatially well-separated due to the sparsity of first-order reflectors, we generalize the intra-cluster model to the mmWave channel model via replication. Since narrow beamwidth increases the number of single-order clusters, we show that the proposed model suits well to MIMO and massive MIMO applications. We illustrate that the model gives matching results with published measurements made in a classroom at 60 GHz. For this specific implementation, while the maximum cluster angle of arrival (AoA) error is 1 degree, mean angle spread error is 9 degrees. The RMS error for the cluster peak power is found to be 2.2 dB.


title: Gauge invariance of the background average effective action

abstract: Using the background field method for the functional renormalization group approach in the case of a generic gauge theory, we study the background field symmetry and gauge dependence of the background average effective action, when the regulator action depends on external fields. The final result is that the symmetry of the average effective action can be maintained for a wide class of regulator functions, but in all cases the dependence of the gauge fixing remains on-shell. The Yang-Mills theory is considered as the main particular example.


title: Better Technical Debt Detection via SURVEYing

abstract: Software analytics can be improved by surveying; i.e. rechecking and (possibly) revising the labels offered by prior analysis. Surveying is a time-consuming task and effective surveyors must carefully manage their time. Specifically, they must balance the cost of further surveying against the additional benefits of that extra effort. This paper proposes SURVEY0, an incremental Logistic Regression estimation method that implements cost/benefit analysis. Some classifier is used to rank the as-yet-unvisited examples according to how interesting they might be. Humans then review the most interesting examples, after which their feedback is used to update an estimator for estimating how many examples are remaining. This paper evaluates SURVEY0 in the context of self-admitted technical debt. As software project mature, they can accumulate "technical debt" i.e. developer decisions which are sub-optimal and decrease the overall quality of the code. Such decisions are often commented on by programmers in the code; i.e. it is self-admitted technical debt (SATD). Recent results show that text classifiers can automatically detect such debt. We find that we can significantly outperform prior results by SURVEYing the data. Specifically, for ten open-source JAVA projects, we can find 83% of the technical debt via SURVEY0 using just 16% of the comments (and if higher levels of recall are required, SURVEY0can adjust towards that with some additional effort).


title: On the perturbed photogravitational restricted five-body problem: the analysis of fractal basins of convergence

abstract: In the framework of photogravitational version of the restricted five-body problem, the existence and stability of the in-plane equilibrium points, the possible regions for motion are explored and analysed numerically, under the combined effect of small perturbations in the Coriolis and centrifugal forces. Moreover, the multivariate version of the Newton-Raphson iterative scheme is applied in an attempt to unveil the topology of the basins of convergence linked with the libration points as function of radiation parameters, and the parameters corresponding to Coriolis and centrifugal forces.


title: A strongly irreducible affine iterated function system with two invariant measures of maximal dimension

abstract: A classical theorem of Hutchinson asserts that if an iterated function system acts on $\mathbb{R}^d$ by similitudes and satisfies the open set condition then it admits a unique self-similar measure with Hausdorff dimension equal to the dimension of the attractor. In the class of measures on the attractor which arise as the projections of shift-invariant measures on the coding space, this self-similar measure is the unique measure of maximal dimension. In the context of affine iterated function systems it is known that there may be multiple shift-invariant measures of maximal dimension if the linear parts of the affinities share a common invariant subspace, or more generally if they preserve a finite union of proper subspaces of $\mathbb{R}^d$. In this note we construct examples where multiple invariant measures of maximal dimension exist even though the linear parts of the affinities do not preserve a finite union of proper subspaces.


title: A Neural Network Architecture for Learning Word-Referent Associations in Multiple Contexts

abstract: This article proposes a biologically inspired neurocomputational architecture which learns associations between words and referents in different contexts, considering evidence collected from the literature of Psycholinguistics and Neurolinguistics. The multi-layered architecture takes as input raw images of objects (referents) and streams of word's phonemes (labels), builds an adequate representation, recognizes the current context, and associates label with referents incrementally, by employing a Self-Organizing Map which creates new association nodes (prototypes) as required, adjusts the existing prototypes to better represent the input stimuli and removes prototypes that become obsolete/unused. The model takes into account the current context to retrieve the correct meaning of words with multiple meanings. Simulations show that the model can reach up to 78% of word-referent association accuracy in ambiguous situations and approximates well the learning rates of humans as reported by three different authors in five Cross-Situational Word Learning experiments, also displaying similar learning patterns in the different learning conditions.


title: Impact of non-Hemiticity on modal strength and correlation in transmission through random open cavities

abstract: The nonorthogonality of eigenfunctions over the volume of non-Hermitian systems determines the nature of waves in complex systems. Here, we show in microwave measurements of the transmission matrix that the non-Hermiticity of open random systems leads to enhanced modal excitation and strong correlation between modes. Modal transmission coefficients reach values comparable to the dimensionless conductance which may be much larger than unity. This is accompanied by strong negative correlation between modal speckle patterns ensuring that net transmission is never larger than the incident power.


title: Inference under Information Constraints II: Communication Constraints and Shared Randomness

abstract: A central server needs to perform statistical inference based on samples that are distributed over multiple users who can each send a message of limited length to the center. We study problems of distribution learning and identity testing in this distributed inference setting and examine the role of shared randomness as a resource. We propose a general-purpose simulate-and-infer strategy that uses only private-coin communication protocols and is sample-optimal for distribution learning. This general strategy turns out to be sample-optimal even for distribution testing among private-coin protocols. Interestingly, we propose a public-coin protocol that outperforms simulate-and-infer for distribution testing and is, in fact, sample-optimal. Underlying our public-coin protocol is a random hash that when applied to the samples minimally contracts the chi-squared distance of their distribution to the uniform distribution.

categories: cs.DS cs.DM cs.IT cs.LG math.IT math.ST stat.TH


title: Interpolation of equation-of-state data

abstract: Aims. We use Hermite splines to interpolate pressure and its derivatives simultaneously, thereby preserving mathematical relations between the derivatives. The method therefore guarantees that thermodynamic identities are obeyed even between mesh points. In addition, our method enables an estimation of the precision of the interpolation by comparing the Hermite-spline results with those of frequent cubic (B-) spline interpolation. Methods. We have interpolated pressure as a function of temperature and density with quintic Hermite 2D-splines. The Hermite interpolation requires knowledge of pressure and its first and second derivatives at every mesh point. To obtain the partial derivatives at the mesh points, we used tabulated values if given or else thermodynamic equalities, or, if not available, values obtained by differentiating B-splines. Results. The results were obtained with the grid of the SAHA-S equation-of-state (EOS) tables. The maximum $lg P$ difference lies in the range from $10^{-9}$ to $10^{-4}$, and $\Gamma_1$ difference varies from $10^{-9}$ to $10^{-3}$. Specifically, for the points of a solar model, the maximum differences are one order of magnitude smaller than the aforementioned values. The poorest precision is found in the dissociation and ionization regions, occurring at $T \sim 1.5\cdot 10^3 - 10^5$ K. The best precision is achieved at higher temperatures, $T>10^5$ K. To discuss the significance of the interpolation errors we compare them with the corresponding difference between two different equation-of-state formalisms, SAHA-S and OPAL 2005. We find that the interpolation errors of the pressure are a few orders of magnitude less than the differences from between the physical formalisms, which is particularly true for the solar-model points.


title: Einstein-Maxwell-scalar black holes: classes of solutions, dyons and extremality

abstract: Spherical black hole (BH) solutions in Einstein-Maxwell-scalar (EMS) models wherein the scalar field is non-minimally coupled to the Maxwell invariant by some coupling function are discussed. We suggest a classification for these models into two classes, based on the properties of the coupling function, which, in particular, allow, or not, the Reissner-Nordstr\"om (RN) BH solution of electrovacuum to solve a given model. Then, a comparative analysis of two illustrative families of solutions, one belonging to each class is performed: $dilatonic$ versus $scalarised$ BHs. By including magnetic charge, that is considering dyons, we show that scalarised BHs can have a smooth extremal limit, unlike purely electric or magnetic solutions. In particular, we study this extremal limit using the entropy function formalism, which provides insight on why both charges are necessary for extremal solutions to exist.


title: Balanced algebraic unknotting, linking forms, and surfaces in three- and four-space

abstract: We provide three 3-dimensional characterizations of the Z-slice genus of a knot, the minimal genus of a locally-flat surface in 4-space cobounding the knot whose complement has cyclic fundamental group: in terms of balanced algebraic unknotting, in terms of Seifert surfaces, and in terms of presentation matrices of the Blanchfield pairing. This result generalizes to a knot in an integer homology 3-sphere and surfaces in certain simply connected signature zero 4-manifolds cobounding this homology sphere. Using the Blanchfield characterization, we obtain effective lower bounds for the Z-slice genus from the linking pairing of the double branched covering of the knot. In contrast, we show that for odd primes p, the linking pairing on the first homology of the p-fold branched covering is determined up to isometry by the action of the deck transformation group on said first homology. As an application of the new upper and lower bounds, we complete the calculation of the Z-slice genus for all prime knots with crossing number up to 12.


title: Tikhonov-Fenichel reduction for parameterized critical manifolds with applications to chemical reaction networks

abstract: We derive a reduction formula for singularly perturbed ordinary differential equations (in the sense of Tikhonov and Fenichel) with a known parameterization of the critical manifold. No a priori assumptions concerning separation of slow and fast variables are made, or necessary.We apply the theoretical results to chemical reaction networks with mass action kinetics admitting slow and fast reactions. For some relevant classes of such systems there exist canonical parameterizations of the variety of stationary points, hence the theory is applicable in a natural manner. In particular we obtain a closed form expression for the reduced system when the fast subsystem admits complex balanced steady states.


title: Out of time ordered effective dynamics of a quartic oscillator

abstract: We study the dynamics of a quantum Brownian particle weakly coupled to a thermal bath. Working in the Schwinger-Keldysh formalism, we develop an effective action of the particle up to quartic terms. We demonstrate that this quartic effective theory is dual to a stochastic dynamics governed by a non-linear Langevin equation. The Schwinger-Keldysh effective theory, or the equivalent non-linear Langevin dynamics, is insufficient to determine the out of time order correlators (OTOCs) of the particle. To overcome this limitation, we construct an extended effective action in a generalised Schwinger-Keldysh framework. We determine the additional quartic couplings in this OTO effective action and show their dependence on the bath's 4-point OTOCs. We analyse the constraints imposed on the OTO effective theory by microscopic reversibility and thermality of the bath. We show that these constraints lead to a generalised fluctuation-dissipation relation between the non-Gaussianity in the distribution of the thermal noise experienced by the particle and the thermal jitter in its damping coefficient. The quartic effective theory developed in this work provides extension of several results previously obtained for the cubic OTO dynamics of a Brownian particle.


title: Detection of similar successive groups in a model with diverging number of variable groups

abstract: In this paper, a linear model with grouped explanatory variables is considered. The idea is to perform an automatic detection of different successive groups of the unknown coefficients under the assumption that the number of groups is of the same order as the sample size. The standard least squares loss function and the quantile loss function are both used together with the fused and adaptive fused penalty to simultaneously estimate and group the unknown parameters. The proper convergence rate is given for the obtained estimators and the upper bound for the number of different successive group is derived. A simulation study is used to compare the empirical performance of the proposed fused and adaptive fused estimators and a real application on the air quality data demonstrates the practical applicability of the proposed methods.


title: He abundances in disc galaxies -- I. Predictions from cosmological chemodynamical simulations

abstract: We investigate how the stellar and gas-phase He abundances evolve as functions of time within simulated star-forming disc galaxies with different star formation histories. We make use of a cosmological chemodynamical simulation for galaxy formation and evolution, which includes star formation, as well as energy and chemical enrichment feedback from asymptotic giant branch stars, core-collapse supernovae, and Type Ia supernovae. The predicted relations between the He mass fraction, $Y$, and the metallicity, $Z$, in the interstellar medium of our simulated disc galaxies depend on the past galaxy star formation history. In particular, $dY/dZ$ is not constant and evolves as a function of time, depending on the specific chemical element that we choose to trace $Z$; in particular, $dY/dX_{\text{O}}$ and $dY/dX_{\text{C}}$ increase as functions of time, whereas $dY/dX_{\text{N}}$ decreases. In the gas-phase, we find negative radial gradients of $Y$, due to the inside-out growth of our simulated galaxy discs as a function of time; this gives rise to longer chemical enrichment time scales in the outer galaxy regions, where we find lower average values for $Y$ and $Z$. Finally, by means of chemical evolution models, in the galactic bulge and inner disc, we predict steeper $Y$ versus age relations at high $Z$ than in the outer galaxy regions. We conclude that, for calibrating the assumed $Y$-$Z$ relation in stellar models, C, N, and C+N are better proxies for the metallicity than O, because they show steeper and less scattered relations.


title: Magnetic activity of the young solar analogue V1358 Ori

abstract: Young, fast rotating single stars can show dramatically different magnetic signatures and levels of magnetic activity as compared with the Sun. While losing angular momentum due to magnetic breaking and mass loss through stellar winds, the stars gradually spin down resulting in decreasing levels of activity. Studying magnetic activity on such solar analogues plays a key role in understanding the evolution of solar-like stars and allows a glimpse into the past of the Sun as well. In order to widen our knowledge of the magnetic evolution of the Sun and solar-like stars, magnetic activity of the young solar analogue V1358 Ori is investigated. Fourier analysis of long-term photometric data is used to derive rotational period and activity cycle length, while spectral synthesis is applied on high resolution spectroscopic data in order to derive precise astrophysical parameters. Doppler imaging is performed to recover surface temperature maps for two subsequent intervals. Cross-correlation of the consecutive Doppler maps is used to derive surface differential rotation. The rotational modulation of the chromospheric activity indicators is also investigated. An activity cycle of $\approx 1600$ days is detected for V1358 Ori. Doppler imaging revealed a surface temperature distribution dominated by a large polar cap with a few weaker features around the equator. This spot configuration is similar to other maps of young solar analogues from the literature, and supports recent model predictions. We detected solar-like surface differential rotation with a surface shear parameter of $\alpha = 0.016 \pm 0.010$ which fits pretty well to our recently proposed empirical relation between rotation and differential rotation. The chromospheric activity indicators showed a rotational modulation.


title: A shuffling theorem for lozenge tilings of doubly-dented hexagons

abstract: MacMahon's theorem on plane partitions yields a simple product formula for tiling number of a hexagon, and Cohn, Larsen and Propp's theorem provides an explicit enumeration for tilings of a dented semihexagon via semi-strict Gelfand--Tsetlin patterns. In this paper, we prove a natural hybrid of the two theorems for hexagons with an arbitrary set of unit triangles removed along a horizontal axis. In particular, we show that the `shuffling' of removed unit triangles only changes the tiling number of the region by a simple multiplicative factor. Our main result generalizes a number of known enumerations and asymptotic enumerations of tilings. We also reveal connections of the main result to the study of symmetric functions and $q$-series.


title: The many facets of exceptional field theory

abstract: Exceptional field theories are the manifestly duality covariant formulations of the target space theories of string/M-theory in the low-energy limit (supergravity) or for certain truncations. These theories feature a rich system of sub-theories, corresponding to different consistent truncations. We review the structure of exceptional field theory and elaborate on various potential applications. In particular, we discuss how exceptional field theories capture some of the the magic supergravity triangles, Hull's $*$-theories related to timelike dualities, as well as the generalized supergravities related to integrable deformations of AdS/CFT.


title: Inferring Global Dynamics of a Black-Box System Using Machine Learning

abstract: We present that, instead of establishing the equations of motion, one can model-freely reveal the dynamical properties of a black-box system using a learning machine. Trained only by a segment of time series of a state variable recorded at present parameters values, the dynamics of the learning machine at different training stages can be mapped to the dynamics of the target system along a particular path in its parameter space, following an appropriate training strategy that monotonously decreases the cost function. This path is important, because along that, the primary dynamical properties of the target system will subsequently emerge, in the simple-to-complex order, matching closely the evolution law of certain self-evolved systems in nature. Why such a path can be reproduced is attributed to our training strategy. This particular function of the learning machine opens up a novel way to probe the global dynamical properties of a black-box system without artificially establish the equations of motion, and as such it might have countless applications. As an example, this method is applied to infer what dynamical stages a variable star has experienced and how it will evolve in future, by using the light curve observed presently.


title: Longitudinal Dynamic versus Kinematic Models for Car-Following Control Using Deep Reinforcement Learning

abstract: The majority of current studies on autonomous vehicle control via deep reinforcement learning (DRL) utilize point-mass kinematic models, neglecting vehicle dynamics which includes acceleration delay and acceleration command dynamics. The acceleration delay, which results from sensing and actuation delays, results in delayed execution of the control inputs. The acceleration command dynamics dictates that the actual vehicle acceleration does not rise up to the desired command acceleration instantaneously due to dynamics. In this work, we investigate the feasibility of applying DRL controllers trained using vehicle kinematic models to more realistic driving control with vehicle dynamics. We consider a particular longitudinal car-following control, i.e., Adaptive Cruise Control (ACC), problem solved via DRL using a point-mass kinematic model. When such a controller is applied to car following with vehicle dynamics, we observe significantly degraded car-following performance. Therefore, we redesign the DRL framework to accommodate the acceleration delay and acceleration command dynamics by adding the delayed control inputs and the actual vehicle acceleration to the reinforcement learning environment state, respectively. The training results show that the redesigned DRL controller results in near-optimal control performance of car following with vehicle dynamics considered when compared with dynamic programming solutions.


title: Multitask Learning of Temporal Connectionism in Convolutional Networks using a Joint Distribution Loss Function to Simultaneously Identify Tools and Phase in Surgical Videos

abstract: Surgical workflow analysis is of importance for understanding onset and persistence of surgical phases and individual tool usage across surgery and in each phase. It is beneficial for clinical quality control and to hospital administrators for understanding surgery planning. Video acquired during surgery typically can be leveraged for this task. Currently, a combination of convolutional neural network (CNN) and recurrent neural networks (RNN) are popularly used for video analysis in general, not only being restricted to surgical videos. In this paper, we propose a multi-task learning framework using CNN followed by a bi-directional long short term memory (Bi-LSTM) to learn to encapsulate both forward and backward temporal dependencies. Further, the joint distribution indicating set of tools associated with a phase is used as an additional loss during learning to correct for their co-occurrence in any predictions. Experimental evaluation is performed using the Cholec80 dataset. We report a mean average precision (mAP) score of 0.99 and 0.86 for tool and phase identification respectively which are higher compared to prior-art in the field.


title: Asymptotic Semigroups and Two-sided Weak Orders

abstract: Various partial orders related to the structures of dual canonical monoids are investigated. It is shown that the nilpotent variety of a dual canonical monoid is equidimensional; its dimension is found. It is shown in type A that certain intervals of the Putcha poset of a dual canonical monoid are isomorphic to the Renner monoids of matrices. The notion of a two-sided weak order on a normal reductive monoid is introduced. A criterion, in terms of type maps, for the covering relations in a two-sided weak order to have degree 2 is found. It is shown that, for the unique equivariant divisor of a dual canonical monoid (the asymptotic semigroup), the covering relations of the two-sided weak order are always of degree 1. These computations provide new insights for the two-sided weak orders on Coxeter groups. In type A, some enumerative results for the covering relations are presented.


title: Corson reflections

abstract: A reflection principle for Corson compacta holds in the forcing extension obtained by Levy-collapsing a supercompact cardinal to~$\aleph_2$. In this model, a compact Hausdorff space is Corson if and only if all of its continuous images of weight~$\aleph_1$ are Corson compact. We use the Gelfand--Naimark duality, and our results are stated in terms of unital abelian \cstar-algebras.


title: DeepCABAC: Context-adaptive binary arithmetic coding for deep neural network compression

abstract: We present DeepCABAC, a novel context-adaptive binary arithmetic coder for compressing deep neural networks. It quantizes each weight parameter by minimizing a weighted rate-distortion function, which implicitly takes the impact of quantization on to the accuracy of the network into account. Subsequently, it compresses the quantized values into a bitstream representation with minimal redundancies. We show that DeepCABAC is able to reach very high compression ratios across a wide set of different network architectures and datasets. For instance, we are able to compress by x63.6 the VGG16 ImageNet model with no loss of accuracy, thus being able to represent the entire network with merely 8.7MB.


title: A Framework for Data-Based Turbulent Combustion Closure: A Priori Validation

abstract: Experimental multi-scalar measurements in laboratory flames have provided important databases for the validation of turbulent combustion closure models. In this work, we present a framework for data-based closure in turbulent combustion and establish an a priori validation of this framework. The approach is based on the construction of joint probability density functions (PDFs) and conditional statistics using experimental data based on the parameterization of the composition space using principal component analysis (PCA). The PCA on the data identifies key parameters, principal components (PCs), on which joint scalar PDFs and conditional scalar means can be constructed. To enable a generic implementation for the construction of joint scalar PDFs, we use the multi-dimensional kernel density estimation (KDE) approach. An a priori validation of the PCA-KDE approach is carried out using the Sandia piloted jet turbulent flames D, E and F. The analysis of the results suggests that a few PCs are adequate to reproduce the statistics, resulting in a low-dimensional representation of the joint scalars PDFs and the scalars' conditional means. A reconstruction of the scalars' means and RMS statistics are in agreement with the corresponding statistics extracted directly from the experimental data. We also propose one strategy to recover missing species and construct conditional means for the reaction rates based on a variation of the pairwise-mixing stirred reactor (PMSR) model. The model is demonstrated using numerical simulations based on the one-dimensional turbulence (ODT) model for the same flames.


title: Locally Differentially Private Frequency Estimation with Consistency

abstract: Local Differential Privacy (LDP) protects user privacy from the data collector. LDP protocols have been increasingly deployed in the industry. A basic building block is frequency oracle (FO) protocols, which estimate frequencies of values. While several FO protocols have been proposed, the design goal does not lead to optimal results for answering many queries. In this paper, we show that adding post-processing steps to FO protocols by exploiting the knowledge that all individual frequencies should be non-negative and they sum up to one can lead to significantly better accuracy for a wide range of tasks, including frequencies of individual values, frequencies of the most frequent values, and frequencies of subsets of values. We consider 10 different methods that exploit this knowledge differently. We establish theoretical relationships between some of them and conducted extensive experimental evaluations to understand which methods should be used for different query tasks.


title: On local analysis

abstract: We extend to Gaussian distributions a result providing smoothed analysis estimates for condition numbers given as relativized distances to illposedness. We also introduce a notion of local analysis meant to capture the behavior of these condition numbers around a point.


title: Semidefinite relaxation of multi-marginal optimal transport for strictly correlated electrons in second quantization

abstract: We consider the strictly correlated electron (SCE) limit of the fermionic quantum many-body problem in the second-quantized formalism. This limit gives rise to a multi-marginal optimal transport (MMOT) problem. Here the marginal state space for our MMOT problem is the binary set $\{0,1\}$, and the number of marginals is the number $L$ of sites in the model. The costs of storing and computing the exact solution of the MMOT problem both scale exponentially with respect to $L$. We propose an efficient convex relaxation which can be solved by semidefinite programming (SDP). In particular, the semidefinite constraint is only of size $2L\times 2L$. Moreover, the SDP-based method yields an approximation of the dual potential needed to the perform self-consistent field iteration in the so-called Kohn-Sham SCE framework, which, once converged, yields a lower bound for the total energy of the system. We demonstrate the effectiveness of our methods on spinless and spinful Hubbard-type models. Numerical results indicate that our relaxation methods yield tight lower bounds for the optimal cost, in the sense that the error due to the semidefinite relaxation is much smaller than the intrinsic modeling error of the Kohn-Sham SCE method. We also describe how our relaxation methods generalize to arbitrary MMOT problems with pairwise cost functions.


title: Unconventional light scattering from glassy photonic films and metasurfaces

abstract: The propagation of light through a random medium is an important problem in photonics. When the random fluctuations of the orientation for individual rods were introduced to the ideal woodpile photonic structure, a crossover from Laue diffraction to randomly scattered fields which is similar in appearance to speckle patterns was observed and investigated. Unexpected interplay between order and disorder was discovered from anisotropic glassy samples when orientational disorder was added only in one direction of square woodpile structure. It is found that the ordered sets of rods produced disordered patterns and vice versa the disordered sets of rods produced ordered patterns that continue to be bright and sharp with increasing disorder. To explain this effect, it was demonstrated theoretically and experimentally that the light scattering can be described purely in terms of the intersection points of the rods.


title: Rare top quark production and decays at ATLAS and CMS

abstract: The most recent studies in the top quark sector are reviewed with the focus on the rare production mechanisms and suppressed decays. The experimental results obtained with the ATLAS and CMS detectors in proton-proton collisions at the center-of-mass energy of 13 TeV include the measurements of the associated production of top quark pairs with vector bosons (${\rm t\bar{t}W}$, ${\rm t\bar{t}Z}$, ${\rm t\bar{t}\gamma}$), the first evidence for the ${\rm t(\bar{t})\gamma q}$ process, the first observation of the ${\rm t(\bar{t})Zq}$ production, the study of the ${\rm t\bar{t}+b\bar{b}}$ and ${\rm t\bar{t}+t\bar{t}}$ processes, as well as searches for lepton flavour violation in top quark decays and effective field theory interpretations. The experimental results show good agreement with the theoretical predictions.


title: Towards Neural Decompilation

abstract: We address the problem of automatic decompilation, converting a program in low-level representation back to a higher-level human-readable programming language. The problem of decompilation is extremely important for security researchers. Finding vulnerabilities and understanding how malware operates is much easier when done over source code. The importance of decompilation has motivated the construction of hand-crafted rule-based decompilers. Such decompilers have been designed by experts to detect specific control-flow structures and idioms in low-level code and lift them to source level. The cost of supporting additional languages or new language features in these models is very high. We present a novel approach to decompilation based on neural machine translation. The main idea is to automatically learn a decompiler from a given compiler. Given a compiler from a source language S to a target language T , our approach automatically trains a decompiler that can translate (decompile) T back to S . We used our framework to decompile both LLVM IR and x86 assembly to C code with high success rates. Using our LLVM and x86 instantiations, we were able to successfully decompile over 97% and 88% of our benchmarks respectively.


title: Planar pure braids on six strands

abstract: The group of planar (or flat) pure braids on $n$ strands, also known as the pure twin group, is the fundamental group of the configuration space $F_{n,3}(\mathbb{R})$ of $n$ labelled points in $\mathbb{R}$ no three of which coincide. The planar pure braid groups on 3, 4 and 5 strands are free. In this note we describe the planar pure braid group on 6 strands: it is a free product of the free group on 71 generators and 20 copies of the free abelian group of rank two.


title: Tools for analyzing R code the tidy way

abstract: With the current emphasis on reproducibility and replicability, there is an increasing need to examine how data analyses are conducted. In order to analyze the between researcher variability in data analysis choices as well as the aspects within the data analysis pipeline that contribute to the variability in results, we have created two R packages: matahari and tidycode. These packages build on methods created for natural language processing; rather than allowing for the processing of natural language, we focus on R code as the substrate of interest. The matahari package facilitates the logging of everything that is typed in the R console or in an R script in a tidy data frame. The tidycode package contains tools to allow for analyzing R calls in a tidy manner. We demonstrate the utility of these packages as well as walk through two examples.


title: Improved Description of Perovskite Oxide Crystal Structure and Electronic Properties using Self-Consistent Hubbard $U$ Corrections from ACBN0

abstract: The wide variety of complex physical behavior exhibited in transition metal oxides, particularly the perovskites A$B$O$_3$, makes them a material family of interest in many research areas, but the drastically different electronic structures possible in these oxides raises challenges in describing them accurately within density functional theory (DFT) and related methods. Here we evaluate the ability of the ACBN0, a recently developed first-principles approach to computing the Hubbard $U$ correction self-consistently, to describe the structural and electronic properties of the first-row transition metal perovskites with $\left(B=\textrm{V}-\textrm{Ni} \right)$. ACBN0 performs competitively with hybrid functional approaches such as the Heyd-Scuseria-Ernzerhof (HSE) functional even when they are optimized empirically, at a fraction of the computational cost. ACBN0 also describes both the structure and band gap of the oxides more accurately than a conventional Hubbard $U$ correction performed by using $U$ values taken from the literature.


title: Cosmicflows-3: Cosmography of the Local Void

abstract: Cosmicflows-3 distances and inferred peculiar velocities of galaxies have permitted the reconstruction of the structure of over and under densities within the volume extending to 0.05c. This study focuses on the under dense regions, particularly the Local Void that lies largely in the zone of obscuration and consequently has received limited attention. Major over dense structures that bound the Local Void are the Perseus-Pisces and Norma-Pavo-Indus filaments separated by 8,500 km/s. The void network of the universe is interconnected and void passages are found from the Local Void to the adjacent very large Hercules and Sculptor voids. Minor filaments course through voids. A particularly interesting example connects the Virgo and Perseus clusters, with several substantial galaxies found along the chain in the depths of the Local Void. The Local Void has a substantial dynamical effect, causing a deviant motion of the Local Group of 200-250 km/s. The combined perturbations due to repulsion from the Local Void and attraction toward the Virgo Cluster account for about 50% of the motion of the Local Group in the rest frame given by the cosmic microwave background.


title: Raking and Regression Calibration: Methods to Address Bias from Correlated Covariate and Time-to-Event Error

abstract: Medical studies that depend on electronic health records (EHR) data are often subject to measurement error, as the data are not collected to support research questions under study. These data errors, if not accounted for in study analyses, can obscure or cause spurious associations between patient exposures and disease risk. Methodology to address covariate measurement error has been well developed; however, time-to-event error has also been shown to cause significant bias but methods to address it are relatively underdeveloped. More generally, it is possible to observe errors in both the covariate and the time-to-event outcome that are correlated. We propose regression calibration (RC) estimators to simultaneously address correlated error in the covariates and the censored event time. Although RC can perform well in many settings with covariate measurement error, it is biased for nonlinear regression models, such as the Cox model. Thus, we additionally propose raking estimators which are consistent estimators of the parameter defined by the population estimating equation. Raking can improve upon RC in certain settings with failure-time data, require no explicit modeling of the error structure, and can be utilized under outcome-dependent sampling designs. We discuss features of the underlying estimation problem that affect the degree of improvement the raking estimator has over the RC approach. Detailed simulation studies are presented to examine the performance of the proposed estimators under varying levels of signal, error, and censoring. The methodology is illustrated on observational EHR data on HIV outcomes from the Vanderbilt Comprehensive Care Clinic.


title: ROMEO: A Plug-and-play Software Platform of Robotics-inspired Algorithms for Modeling Biomolecular Structures and Motions

abstract: Motivation: Due to the central role of protein structure in molecular recognition, great computational efforts are devoted to modeling protein structures and motions that mediate structural rearrangements. The size, dimensionality, and non-linearity of the protein structure space present outstanding challenges. Such challenges also arise in robot motion planning, and robotics-inspired treatments of protein structure and motion are increasingly showing high exploration capability. Encouraged by such findings, we debut here ROMEO, which stands for Robotics prOtein Motion ExplOration framework. ROMEO is an open-source, object-oriented platform that allows researchers access to and reproducibility of published robotics-inspired algorithms for modeling protein structures and motions, as well as facilitates novel algorithmic design via its plug-and-play architecture. Availability and implementation: ROMEO is written in C++ and is available in GitLab (https://github.com/). This software is freely available under the Creative Commons license (Attribution and Non-Commercial). Contact: amarda@gmu.edu


title: Behavior Identification and Prediction for a Probabilistic Risk Framework

abstract: Operation in a real world traffic requires autonomous vehicles to be able to plan their motion in complex environments (multiple moving participants). Planning through such environment requires the right search space to be provided for the trajectory or maneuver planners so that the safest motion for the ego vehicle can be identified. Given the current states of the environment and its participants, analyzing the risks based on the predicted trajectories of all the traffic participants provides the necessary search space for the planning of motion. This paper provides a fresh taxonomy of safety / risks that an autonomous vehicle should be able to handle while navigating through traffic. It provides a reference system architecture that needs to be implemented as well as describes a novel way of identifying and predicting the behaviors of the traffic participants using classic Multi Model Adaptive Estimation (MMAE). Preliminary simulation results of the implemented model are included.


title: Determining the Best Method of Calculating the Large Frequency Separation For Stellar Models

abstract: Asteroseismology of solar-like oscillators often relies on the comparisons between stellar models and stellar observations in order to determine the properties of stars. The values of the global seismic parameters, $\nu_\mathrm{max}$ (the frequency where the smoothed amplitude of the oscillations peak) and $\Delta \nu$ (the large frequency separation), are frequently used in grid-based modeling searches. However, the methods by which $\Delta \nu$ is calculated from observed data and how $\Delta \nu$ is calculated from stellar models are not the same. Typically for observed stars, especially for those with low signal-to-noise data, $\Delta \nu$ is calculated by taking the power spectrum of a power spectrum, or with autocorrelation techniques. However, for stellar models, the actual individual mode frequencies are calculated and the average spacing between them directly determined. In this work we try to determine the best way to combine model frequencies in order to obtain $\Delta \nu$ that can be compared with observations. For this we use stars with high signal-to-noise observations from Kepler as well as simulated TESS data of Ball et al. (2018). We find that when determining $\Delta \nu$ from individual mode frequencies the best method is to use the $\ell=0$ modes with either no weighting or with a Gaussian weighting around $\nu_\mathrm{max}$.


title: Geodesic rays, the "Lion-Man" game, and the fixed point property

abstract: This paper focuses on the relation among the existence of different types of curves (such as directional ones, quasi-geodesic or geodesic rays), the (approximate) fixed point property for nonexpansive mappings, and a discrete lion and man game. Our main result holds in the setting of CAT(0) spaces that are additionally Gromov hyperbolic.


title: Cavity-induced mirror-mirror entanglement in a single-atom Raman laser

abstract: We address an experimental scheme to analyze the optical bistability and the entanglement of two movable mirrors coupled to a two-mode laser inside a doubly resonant cavity. With this aim we investigate the master equations of the atom-cavity subsystem in conjunction with the quantum Langevin equations that describe the interaction of the mirror cavity. The parametric amplification-type coupling induced by the two-photon coherence on the optical bistability of the intracavity mean photon numbers is found and investigated. Under this condition, the optical intensities exhibit bistability for all large values of cavity laser detuning. We also provide numerical evidence for the generation of strong entanglement between the movable mirrors and show that it is robust against environmental thermalization.


title: A new non-convex model of the binary asteroid (809) Lundia obtained with the SAGE modelling technique

abstract: We present a new non-convex model of the binary asteroid (809) Lundia. A SAGE (Shaping Asteroids with Genetic Evolution) method using disc-integrated photometry only was used for deriving physical parameters of this binary system. The model of (809) Lundia improves former system's pole solution and gives the ecliptic coordinates of the orbit pole - $\lambda=122^{\circ}$, $\beta=22^{\circ}$, $\sigma=\pm5^{\circ}$ - and the orbital period of $15.41574 \pm 0.00001$ h. For scaling our results we used effective diameter of $D_{eff} = 9.6 \pm 1.1$ km obtained from Spitzer observations. The non-convex shape description of the components permitted a refined calculation of the components' volumes, leading to a density estimation of $2.5\pm0.2$ g/cm$^3$ and macroporosity of 13-23\%. The intermediate-scale features of the model may also offer new clues on the components' origin and evolution.


title: Ingesting High-Velocity Streaming Graphs from Social Media Sources

abstract: Many data science applications like social network analysis use graphs as their primary form of data. However, acquiring graph-structured data from social media presents some interesting challenges. The first challenge is the high data velocity and bursty nature of the social media data. The second challenge is that the complex nature of the data makes the ingestion process expensive. If we want to store the streaming graph data in a graph database, we face a third challenge -- the database is very often unable to sustain the ingestion of high-velocity, high-burst data. We have developed an adaptive buffering mechanism and a graph compression technique that effectively mitigates the problem. A novel aspect of our method is that the adaptive buffering algorithm uses the data rate, the data content as well as the CPU resources of the database machine to determine an optimal data ingestion mechanism. We further show that an ingestion-time graph-compression strategy improves the efficiency of the data ingestion into the database. We have verified the efficacy of our ingestion optimization strategy through extensive experiments.


title: A response to critiques of "The reproducibility of research and the misinterpretation of p-values"

abstract: I proposed (8, 1, 3) that p values should be supplemented by an estimate of the false positive risk (FPR). FPR was defined as the probability that, if you claim that there is a real effect on the basis of p value from a single unbiased experiment, that you will be mistaken and the result has occurred by chance. This is a Bayesian quantity and that means that there is an infinitude of ways to calculate it. My choice of a way to estimate FPR was, therefore, arbitrary. I maintain that it is a reasonable way, and has the advantage of being mathematically simpler than other proposals and easier to understand than other methods. This might make it more easily accepted by users. As always, not every statistician agrees. This paper is a response to a critique of my 2017 paper (1) by Arandjelovic (2)


title: Measuring the Complexity of Packet Traces

abstract: This paper studies the structure of several real-world traces (including Facebook, High-Performance Computing, Machine Learning, and simulation generated traces) and presents a systematic approach to quantify and compare the structure of packet traces based on the entropy contained in the trace file. Insights into the structure of packet traces can lead to improved network algorithms that are optimized toward specific traffic patterns. We then present a methodology to quantify the temporal and non-temporal components of entropy contained in a packet trace, called the trace complexity, using randomization and compression. We show that trace complexity provides unique insights into the characteristics of various applications and argue that there is a need for traffic generation models that preserve the intrinsic structure of empirically measured application traces. We then propose a traffic generator model that is able to produce a synthetic trace that matches the complexity level of its corresponding real-world trace.


title: Coupling Matrix Representation of Nonreciprocal Filters Based on Time Modulated Resonators

abstract: This paper addresses the analysis and design of non-reciprocal filters based on time modulated resonators. We analytically show that time modulating a resonator leads to a set of harmonic resonators composed of the unmodulated lumped elements plus a frequency invariant element that accounts for differences in the resonant frequencies. We then demonstrate that harmonic resonators of different order are coupled through non-reciprocal admittance inverters whereas harmonic resonators of the same order couple with the admittance inverter coming from the unmodulated filter network. This coupling topology provides useful insights to understand and quickly design non-reciprocal filters and permits their characterization using an asynchronously tuned coupled resonators network together with the coupling matrix formalism. Two designed filters, of orders three and four, are experimentally demonstrated using quarter wavelength resonators implemented in microstrip technology and terminated by a varactor on one side. The varactors are biased using coplanar waveguides integrated in the ground plane of the device. Measured results are found to be in good agreement with numerical results, validating the proposed theory.


title: Illusion Mechanisms with Cylindrical Metasurfaces: A General Synthesis Approach

abstract: We explore the use of cylindrical metasurfaces in providing several illusion mechanisms including scattering cancellation and creating fictitious line sources. We present the general synthesis approach that leads to such phenomena by modeling the metasurface with effective polarizability tensors and by applying boundary conditions to connect the tangential components of the desired fields to the required surface polarization current densities that generate such fields. We then use these required surface polarizations to obtain the effective polarizabilities for the synthesis of the metasurface. We demonstrate the use of this general method for the synthesis of metasurfaces that lead to scattering cancellation and illusion effects, and discuss practical scenarios by using loaded dipole antennas to realize the discretized polarization current densities. This study is the first fundamental step that may lead to interesting electromagnetic applications, like stealth technology, antenna synthesis, wireless power transfer, sensors, cylindrical absorbers, etc.


title: Macroscopic forcing method: a tool for turbulence modeling and analysis of closures

abstract: This study presents a numerical procedure, which we call the macroscopic forcing method (MFM), which reveals the differential operators acting on the mean fields of quantities transported by underlying fluctuating flows. Specifically, MFM can precisely determine the eddy diffusivity operator, or more broadly said, it can reveal differential operators associated with turbulence closure for scalar and momentum transport. We present this methodology by considering canonical problems with increasing complexity. Starting from the well-known problem of dispersion of passive scalars by parallel flows we elucidate the basic steps in quantitative determination of the eddy viscosity operators using MFM. Utilizing the operator representation in Fourier space, we obtain a stand-alone compact analytical operator that can accurately capture the nonlocal mixing effects. Furthermore, a cost-effective generalization of MFM for analysis of nonhomogeneous and wall-bounded flows is developed and is comprehensively discussed through a demonstrative example. Extension of MFM for analysis of momentum transport is theoretically constructed through the introduction of a generalized momentum transport equation. We show that closure operators obtained through MFM analysis of this equation provide the exact RANS solutions obtained through averaging of the Navier-Stokes equation. We introduce MFM as an effective tool for quantitative understanding of non-Boussinesq effects and assessment of model forms in turbulence closures, particularly, the effects associated with anisotropy and nonlocality of macroscopic mixing.


title: The $A_2$ Rogers-Ramanujan identities revisited

abstract: In this note we show how to rederive the $A_2$ Rogers-Ramanujan identities proven by Andrews, Schilling and Warnaar using cylindric partitions. This paper is dedicated to George Andrews for his $80^{th}$ birthday.


title: Regularity of the Density of SRB Measures for Solenoidal Attractors

abstract: We show that a class of higher-dimensional hyperbolic endomorphisms admit absolutely continuous invariant probabilities whose density are regular. The maps we consider are given by $T(x,y) = (E (x), C(y) + f(x) )$, where $E$ is a linear expanding map of $\mathbb{T}$, $C$ is a linear contracting map of $\mathbb{R}^d$, $f$ is in $C^r(\mathbb{T}^u,\mathbb{R}^d)$ and $r \geq 2$. We prove that if $|(\det C)(\det E)| \|C^{-1}\|^{-2s}>1$ for some $s<r-(\frac{u+d}{2}+1)$ and $T$ satisfies a certain transversality condition, then the density of the SRB measure of $T$ is contained in the Sobolev space $H^s(\mathbb{T}^u\times \mathbb{R}^d)$, in particular, if $s>\frac{u+d}{2}$ then the density is $C^k$ for every $k<s-\frac{u+d}{2}$. We also exhibit a condition involving $E$ and $C$ under which this tranversality condition is valid for almost every $f$.


title: Enumeration of a special class of irreducible polynomials in characteristic 2

abstract: A-polynomials were introduced by Meyn and play an important role in the iterative construction of high degree self-reciprocal irreducible polynomials over the field F_2, since they constitute the starting point of the iteration. The exact number of A-polynomials of each degree was given by Niederreiter. Kyuregyan extended the construction of Meyn to arbitrary even finite fields. We relate the A-polynomials in this more general setting to inert places in a certain extension of elliptic function fields and obtain an explicit counting formula for their number. In particular, we are able to show that, except for an isolated exception, there exist A-polynomials of every degree.


title: Stress relaxation in F-actin solutions by severing

abstract: Networks of filamentous actin (F-actin) are important for the mechanics of most animal cells. These cytoskeletal networks are highly dynamic, with a variety of actin-associated proteins that control cross-linking, polymerization and force generation in the cytoskeleton. Inspired by recent rheological experiments on reconstituted solutions of dynamic actin filaments, we report a theoretical model that describes stress relaxation behavior of these solutions in the presence of severing proteins. We show that depending on the kinetic rates of assembly, disassembly, and severing, one can observe both length-dependent and length-independent relaxation behavior.


title: Decrement Operators in Belief Change

abstract: While research on iterated revision is predominant in the field of iterated belief change, the class of iterated contraction operators received more attention in recent years. In this article, we examine a non-prioritized generalisation of iterated contraction. In particular, the class of weak decrement operators is introduced, which are operators that by multiple steps achieve the same as a contraction. Inspired by Darwiche and Pearl's work on iterated revision the subclass of decrement operators is defined. For both, decrement and weak decrement operators, postulates are presented and for each of them a representation theorem in the framework of total preorders is given. Furthermore, we present two sub-types of decrement operators.


title: Leaking Information Through Cache LRU States

abstract: The Least-Recently Used cache replacement policy and its variants are widely deployed in modern processors. This paper shows for the first time in detail that the LRU states of caches can be used to leak information: any access to a cache by a sender will modify the LRU state, and the receiver is able to observe this through a timing measurement. This paper presents LRU timing-based channels both when the sender and the receiver have shared memory, e.g., shared library data pages, and when they are separate processes without shared memory. In addition, the new LRU timing-based channels are demonstrated on both Intel and AMD processors in scenarios where the sender and the receiver are sharing the cache in both hyper-threaded setting and time-sliced setting. The transmission rate of the LRU channels can be up to 600Kbps per cache set in the hyper-threaded setting. Different from the majority of existing cache channels which require the sender to trigger cache misses, the new LRU channels work with the sender only having cache hits, making the channel faster and more stealthy. This paper also demonstrates that the new LRU channels can be used in transient execution attacks, e.g., Spectre. Further, this paper shows that the LRU channels pose threats to existing secure cache designs, and this work demonstrates the LRU channels affect the secure PL cache. The paper finishes by discussing and evaluating possible defenses.


title: Slow-Roll Inflation in Scalar-Tensor Models

abstract: The linear and quadratic perturbations for a scalar-tensor model with non-minimal coupling to curvature, coupling to the Gauss-Bonnet invariant and non-minimal kinetic coupling to the Einstein tensor are developed. The quadratic action for the scalar and tensor perturbations is constructed and the power spectra for the primordial scalar and tensor fluctuations are given. A consistency relation that is useful to discriminate the model from the standard inflation with canonical scalar field was found. For some power-law potentials it is shown that the Introduction of additional interactions, given by non-minimal, kinetic and Gauss-Bonnet couplings, can lower the tensor-to-scalar ratio to values that are consistent with latest observational constraints, and the problem of large fields in chaotic inflation can be avoided.


title: Correspondence between dark energy quantum cosmology and Maxwell equations

abstract: A Friedmann-Robertson-Walker cosmology with dark energy can be modelled using a quintessence field. That system is equivalent to a relativistic particle moving on a two-dimensional conformal spacetime. When the quintessence behaves as a free massless scalar field in a Universe with cosmological constant, the quantized version of that theory can lead to a Supersymmetric Majorana quantum cosmology. The purpose of this work is to show that such quantum cosmological model correspond to the Maxwell equations for electromagnetic waves propagating in a medium with specific values for its relative permittivity and relative permeability. The form of those media parameters are calculated, implying that a Majorana quantum cosmology can be studied in an analogue electromagnetic system.


title: Weak polynomial identities for a vector space with a symmetric bilinear form

abstract: Let $V_k$ be a $k$-dimensional vector space with a non-degenerate symmetric bilinear form over a field $K$ of characteristic 0 and let $C_k$ be the Clifford algebra on $V_k$. We study the weak polynomial identities of the pair $(C_k,V_k)$. We establish that all they follow from $[x_1^2,x_2]=0$ when $k=\infty$ and from $[x_1^2,x_2]=0$ and $S_{k+1}(x_1,\ldots,x_{k+1})=0$ when $k<\infty$. We also prove that the weak identity $[x_1^2,x_2]=0$ satisfies the Specht property. As a consequence we obtain a new proof of the theorem of Razmyslov that the weak Lie polynomial identities of the pair $(M_2(K),sl_2(K))$ follow from $[x_1^2,x_2]=0$.


title: Robust sound event detection in bioacoustic sensor networks

abstract: Bioacoustic sensors, sometimes known as autonomous recording units (ARUs), can record sounds of wildlife over long periods of time in scalable and minimally invasive ways. Deriving per-species abundance estimates from these sensors requires detection, classification, and quantification of animal vocalizations as individual acoustic events. Yet, variability in ambient noise, both over time and across sensors, hinders the reliability of current automated systems for sound event detection (SED), such as convolutional neural networks (CNN) in the time-frequency domain. In this article, we develop, benchmark, and combine several machine listening techniques to improve the generalizability of SED models across heterogeneous acoustic environments. As a case study, we consider the problem of detecting avian flight calls from a ten-hour recording of nocturnal bird migration, recorded by a network of six ARUs in the presence of heterogeneous background noise. Starting from a CNN yielding state-of-the-art accuracy on this task, we introduce two noise adaptation techniques, respectively integrating short-term (60 milliseconds) and long-term (30 minutes) context. First, we apply per-channel energy normalization (PCEN) in the time-frequency domain, which applies short-term automatic gain control to every subband in the mel-frequency spectrogram. Secondly, we replace the last dense layer in the network by a context-adaptive neural network (CA-NN) layer. Combining them yields state-of-the-art results that are unmatched by artificial data augmentation alone. We release a pre-trained version of our best performing system under the name of BirdVoxDetect, a ready-to-use detector of avian flight calls in field recordings.


title: Harmonic surfaces in the Cayley plane

abstract: We consider the twistor theory of nilconformal harmonic maps from a Riemann surface into the Cayley plane $\mathbf{O} P^2=F_4/\mathrm{Spin}(9)$. By exhibiting this symmetric space as a submanifold of the Grassmannian of $10$-dimensional subspaces of the fundamental representation of $F_4$, techniques and constructions similar to those used in earlier works on twistor constructions of nilconformal harmonic maps into classical Grassmannians can also be applied in this case. The originality of our approach lies on the use of the classification of nilpotent orbits in Lie algebras as described by D. Djokovi\'{c}.


title: Discrete Prolate Spheroidal Wave Functions: Further spectral analysis and some related applications

abstract: For fixed $W\in \big(0,\frac{1}{2}\big)$ and positive integer $N\geq 1,$ the discrete prolate spheroidal wave functions (DPSWFs), denoted by $U_{k,W}^N,$ $0\leq k\leq N-1$ form the set of the eigenfunctions of the positive and finite rank integral operator $\widetilde Q_{N,W},$ defined on $L^2(-1/2,1/2),$ with kernel $K_N(x,y)=\frac{\sin(N\pi(x-y))}{\sin(\pi(x-y))}\, \mathbf 1_{[-W,W]}(y).$ It is well known that the DPSWF's have a wide range of classical as well as recent signal processing applications. These applications rely heavily on the properties of the DPSWFs as well as the behaviour of their eigenvalues $\widetilde \lambda_{k,N}(W).$ In his pioneer work \cite{Slepian}, D. Slepian has given the properties of the DPSWFs, their asymptotic approximations as well as the asymptotic behaviour and asymptotic decay rate of these eigenvalues. In this work, we give further properties as well as new non-asymptotic decay rates of the spectrum of the operator $\widetilde Q_{N,W}.$ In particular, we show that each eigenvalue $\widetilde \lambda_{k,N}(W)$ is up to a small constant bounded above by the corresponding eigenvalue, associated with the classical prolate spheroidal wave functions (PSWFs). Then, based on the well established results concerning the distribution and the decay rates of the eigenvalues associated with the PSWFs, we extend these results to the eigenvalues $\widetilde \lambda_{k,N}(W)$. Also, we show that the DPSWFs can be used for the approximation of classical band-limited functions and they are well adapted for the approximation of functions from periodic Sobolev spaces. Finally, we provide the reader with some numerical examples that illustrate the different results of this work.


title: Importance of Coordination and Cultural Diversity for an Efficient and Flexible Manufacturing System

abstract: Manufacturing systems of the future need to have flexible resources and flexible routing to produce extremely personalized products, even of lot size equal to one. In this paper, we have proposed a framework, which is designed to achieve this goal. Towards this, we have integrated an established cultural evolution model to achieve desired flexibility of resources and acceptable routing time. Promising results are evidenced through a simple proof-of-concept simulation.


title: Pump frequency resonances for light-induced incipient superconductivity in YBa$_2$Cu$_3$O$_{6.5}$

abstract: Optical excitation in the cuprates has been shown to induce transient superconducting correlations above the thermodynamic transition temperature, $T_C$, as evidenced by the terahertz frequency optical properties in the non-equilibrium state. In YBa$_2$Cu$_3$O$_{6+x}$ this phenomenon has so far been associated with the nonlinear excitation of certain lattice modes and the creation of new crystal structures. In other compounds, like La$_{2-x}$Ba$_x$CuO$_4$, similar effects were reported also for excitation at near infrared frequencies, and were interpreted as a signature of the melting of competing orders. However, to date it has not been possible to systematically tune the pump frequency widely in any one compound, to comprehensively compare the frequency dependent photo-susceptibility for this phenomenon. Here, we make use of a newly developed optical parametric amplifier, which generates widely tunable high intensity femtosecond pulses, to excite YBa$_2$Cu$_3$O$_{6.5}$ throughout the entire optical spectrum (3 - 750 THz). In the far-infrared region (3 - 25 THz), signatures of non-equilibrium superconductivity are induced only for excitation of the 16.4 THz and 19.2 THz vibrational modes that drive $c$-axis apical oxygen atomic positions. For higher driving frequencies (25 - 750 THz), a second resonance is observed around the charge transfer band edge at ~350 THz. These observations highlight the importance of coupling to the electronic structure of the CuO$_2$ planes, either mediated by a phonon or by charge transfer.


title: Double ionization of a three-electron atom: Spin correlation effects

abstract: We study the effects of spin degrees of freedom and wave function symmetries on double ionization in three-electron systems. Each electron is assigned one spatial degree of freedom. The resulting three-dimensional Schr\"odinger equation is integrated numerically using grid-based Fourier transforms. We reveal three-electron effects on the double ionization yield by comparing signals for different ionization channels. We explain our findings by the existence of fundamental differences between three-electronic and truly two-electronic spin-resolved ionization schemes. We find, for instance, that double ionization from a three-electron system is dominated by electrons that have the opposite spin.


title: Persistent Threshold Dynamics with Recovery in Complex Networks

abstract: Threshold rules of spreading in binary-state networks lead to cascades. We study persistent cascade-recovery dynamics on quasi-robust networks, i.e., networks which are robust against small trigger but may collapse for larger one. It is observed that depending on the relative rate of triggering and recovery, the network falls into one of the two dynamical phases: collapsing or active phase. We devise an analytical framework which characterizes not only the critical behavior but also the temporal evolution of network activity in both phases. Agent-based simulation results show good agreement with the analytical calculations, indicating strong predicative power of our method for persistent cascade dynamics in complex networks.


title: Why Machines Cannot Learn Mathematics, Yet

abstract: Nowadays, Machine Learning (ML) is seen as the universal solution to improve the effectiveness of information retrieval (IR) methods. However, while mathematics is a precise and accurate science, it is usually expressed by less accurate and imprecise descriptions, contributing to the relative dearth of machine learning applications for IR in this domain. Generally, mathematical documents communicate their knowledge with an ambiguous, context-dependent, and non-formal language. Given recent advances in ML, it seems canonical to apply ML techniques to represent and retrieve mathematics semantically. In this work, we apply popular text embedding techniques to the arXiv collection of STEM documents and explore how these are unable to properly understand mathematics from that corpus. In addition, we also investigate the missing aspects that would allow mathematics to be learned by computers.


title: Conditionally-additive-noise Models for Structure Learning

abstract: Constraint-based structure learning algorithms infer the causal structure of multivariate systems from observational data by determining an equivalent class of causal structures compatible with the conditional independencies in the data. Methods based on additive-noise (AN) models have been proposed to further discriminate between causal structures that are equivalent in terms of conditional independencies. These methods rely on a particular form of the generative functional equations, with an additive noise structure, which allows inferring the directionality of causation by testing the independence between the residuals of a nonlinear regression and the predictors (nrr-independencies). Full causal structure identifiability has been proven for systems that contain only additive-noise equations and have no hidden variables. We extend the AN framework in several ways. We introduce alternative regression-free tests of independence based on conditional variances (cv-independencies). We consider conditionally-additive-noise (CAN) models, in which the equations may have the AN form only after conditioning. We exploit asymmetries in nrr-independencies or cv-independencies resulting from the CAN form to derive a criterion that infers the causal relation between a pair of variables in a multivariate system without any assumption about the form of the equations or the presence of hidden variables.


title: Introduction to the physics of solar eruptions and their space weather impact

abstract: The physical processes, which drive powerful solar eruptions, play an important role in our understanding of the Sun-Earth connection. In this Special Issue, we firstly discuss how magnetic fields emerge from the solar interior to the solar surface, to build up active regions, which commonly host large-scale coronal disturbances, such as coronal mass ejections (CMEs). Then, we discuss the physical processes associated with the driving and triggering of these eruptions, the propagation of the large-scale magnetic disturbances through interplanetary space and the interaction of CMEs with Earth's magnetic field. The acceleration mechanisms for the solar energetic particles related to explosive phenomena (e.g. flares and/or CMEs) in the solar corona are also discussed. The main aim of this Issue, therefore, is to encapsulate the present state-of-the-art in research related to the genesis of solar eruptions and their space-weather implications. This article is part of the theme issue 'Solar eruptions and their space weather impact'.


title: Ultrafast Quantum-path Interferometry Revealing the Generation Process of Coherent Phonons

abstract: Optical dual-pulse pumping actively creates quantum-mechanical superposition of the electronic and phononic states in a bulk solid. We here made transient reflectivity measurements in an n-GaAs using a pair of relative-phase-locked femtosecond pulses and found characteristic interference fringes. This is a result of quantum-path interference peculiar to the dual-pulse excitation as indicated by theoretical calculation. Our observation reveals that the pathway of coherent phonon generation in the n-GaAs is impulsive stimulated Raman scattering at the displaced potential due to the surface-charge field, even though the photon energy lies in the opaque region.


title: Non-Abelian basis tensor gauge theory

abstract: Basis tensor gauge theory is a vierbein analog reformulation of ordinary gauge theories in which the difference of local field degrees of freedom has the interpretation of an object similar to a Wilson line. Here we present a non-Abelian basis tensor gauge theory formalism. Unlike in the Abelian case, the map between the ordinary gauge field and the basis tensor gauge field is nonlinear. To test the formalism, we compute the beta function and the two-point function at the one-loop level in non-Abelian basis tensor gauge theory and show that it reproduces the well-known results from the usual formulation of non-Abelian gauge theory.


title: Efficient Synthesis with Probabilistic Constraints

abstract: We consider the problem of synthesizing a program given a probabilistic specification of its desired behavior. Specifically, we study the recent paradigm of distribution-guided inductive synthesis (DIGITS), which iteratively calls a synthesizer on finite sample sets from a given distribution. We make theoretical and algorithmic contributions: (i) We prove the surprising result that DIGITS only requires a polynomial number of synthesizer calls in the size of the sample set, despite its ostensibly exponential behavior. (ii) We present a property-directed version of DIGITS that further reduces the number of synthesizer calls, drastically improving synthesis performance on a range of benchmarks.


title: Propagating large open quantum systems towards their steady states: cluster implementation of the time-evolving block decimation scheme

abstract: Many-body quantum systems are subjected to the Curse of Dimensionality: The dimension of the Hilbert space $\mathcal{H}$, where these systems live in, grows exponentially with systems' 'size' (number of their components, "bodies"). It means that, in order to specify a state of a quantum system, we need a description whose length grows exponentially with the system size. However, with some systems it is possible to escape the curse by using low-rank tensor approximations known as `matrix-product state/operator (MPS/O) representation' in the quantum community and `tensor-train decomposition' among applied mathematicians. Motivated by recent advances in computational quantum physics, we consider chains of $N$ spins coupled by nearest-neighbor interactions. The spins are subjected to an action coming from the environment. Spatially disordered interaction and environment-induced decoherence drive systems into non-trivial asymptotic states. The dissipative evolution is modeled with a Markovian master equation in the Lindblad form. By implementing the MPO technique and propagating system states with the time-evolving block decimation (TEBD) scheme (which allows to keep the length of the state descriptions fixed), it is in principle possible to reach the corresponding steady states. We propose and realize a cluster implementation of this idea. The implementation on four nodes allowed us to resolve steady states of the model systems with $N = 128$ spins.


title: Central limit theorems for combinatorial optimization problems on sparse Erd\H{o}s-R\'enyi graphs

abstract: For random combinatorial optimization problems, there has been much progress in establishing laws of large numbers and computing limiting constants for the optimal value of various problems. However, there has not been as much success in proving central limit theorems. This paper introduces a method for establishing central limit theorems in the sparse graph setting. It works for problems which display a key property which has been variously called "endogeny", "long-range independence", and "replica symmetry" in the literature. Examples of such problems are maximum weight matching, $\lambda$-diluted minimum matching, and optimal edge cover.


title: SmartNight: Turning Off the Lights on Android

abstract: Smartphone users benefit from content with dark color schemes: increasingly common OLED displays are more power efficient the darker the display, and many users prefer a dark display for night time use. Despite these benefits, many applications and the majority of web content are drawn with white backgrounds. There are many partial solutions to darken the displayed content, but none work in all situations. Enter SmartNight, a content-aware solution to dynamically darken content on Android. By trading off content fidelity, Android with SmartNight displays content with nearly 90% lower average picture level. It is implemented in the Android framework, and requires no external support. It seamlessly incorporates existing solutions, making it a bridge between the state-of-the-art and future solutions.


title: Desfuncionalizar para Provar

abstract: This paper explores the idea of using defunctionalization as a proof technique for higher-order programs. Defunctionalization builds on substituting functional values by a first-order representation. Thus, its interest is that one can use an existing program verification tool, without further extensions in order to support higher-order. This papers illustrates and discusses this approach by means of several running examples, built and verified using the Why3 verification framework.


title: A Bi-Directional Co-Design Approach to Enable Deep Learning on IoT Devices

abstract: Developing deep learning models for resource-constrained Internet-of-Things (IoT) devices is challenging, as it is difficult to achieve both good quality of results (QoR), such as DNN model inference accuracy, and quality of service (QoS), such as inference latency, throughput, and power consumption. Existing approaches typically separate the DNN model development step from its deployment on IoT devices, resulting in suboptimal solutions. In this paper, we first introduce a few interesting but counterintuitive observations about such a separate design approach, and empirically show why it may lead to suboptimal designs. Motivated by these observations, we then propose a novel and practical bi-directional co-design approach: a bottom-up DNN model design strategy together with a top-down flow for DNN accelerator design. It enables a joint optimization of both DNN models and their deployment configurations on IoT devices as represented as FPGAs. We demonstrate the effectiveness of the proposed co-design approach on a real-life object detection application using Pynq-Z1 embedded FPGA. Our method obtains the state-of-the-art results on both QoR with high accuracy (IoU) and QoS with high throughput (FPS) and high energy efficiency.


title: Stabilization of two strongly coupled hyperbolic equations in exterior domains

abstract: In this paper we study the behavior of the total energy and the $L^2$-norm of solutions of two coupled hyperbolic equations by velocities in exterior domains. Only one of the two equations is directly damped by a localized damping term. We show that, when the damping set contains the coupling one and the coupling term is effective at infinity and on captive region, then the total energy decays uniformly and the $L^2$-norm of smooth solutions is bounded. In the case of two Klein-Gordon equations with equal speeds we deduce an exponential decay of the energy.


title: Insight into the microphysics of antigorite deformation from spherical nanoindentation

abstract: The mechanical behavior of antigorite strongly influences the strength and deformation of the subduction interface. Although there is microstructural evidence elucidating the nature of brittle deformation at low pressures, there is often conflicting evidence regarding the potential for plastic deformation in the ductile regime at higher pressures. Here, we present a series of spherical nanoindentation experiments on aggregates of natural antigorite. These experiments effectively investigate the single-crystal mechanical behavior because the volume of deformed material is significantly smaller than the grain size. Individual indents reveal elastic loading followed by yield and strain hardening. The magnitude of the yield stress is a function of crystal orientation, with lower values associated with indents parallel to the basal plane. Unloading paths reveal more strain recovery than expected for purely elastic unloading. The magnitude of inelastic strain recovery is highest for indents parallel to the basal plane. We also imposed indents with cyclical loading paths, and observed strain energy dissipation during unloading-loading cycles conducted up to a fixed maximum indentation load and depth. The magnitude of this dissipated strain energy was highest for indents parallel to the basal plane. Subsequent scanning electron microscopy revealed surface impressions accommodated by shear cracks and a general lack of lattice misorientation around indents, indicating the absence of dislocations. Based on these observations, we suggest that antigorite deformation at high pressures is dominated by sliding on shear cracks. We develop a microphysical model that is able to quantitatively explain the Young's modulus and dissipated strain energy data during cyclic loading experiments, based on either frictional or cohesive sliding of an array of cracks contained in the basal plane.


title: On classical solutions of the KdV equation

abstract: \begin{abstract} We show that if the initial profile $q\left( x\right) $ for the Korteweg-de Vries (KdV) equation is essentially semibounded from below and $\int^{\infty }x^{5/2}\left\vert q\left( x\right) \right\vert dx<\infty,$ (no decay at $-\infty$ is required) then the KdV has a unique global classical solution given by a determinant formula. This result is best known to date. \end{abstract}


title: The effect of a positive bound state on the KdV solution. A case study

abstract: We consider a slowly decaying oscillatory potential such that the corresponding 1D Schr\"odinger operator has a positive eigenvalue embedded into the absolutely continuous spectrum. This potential does not fall into a known class of initial data for which the Cauchy problem for the Korteweg-de Vries (KdV) equation can be solved by the inverse scattering transform. We nevertheless show that the KdV equation with our potential does admit a closed form classical solution in terms of Hankel operators. Comparing with rapidly decaying initial data our solution gains a new term responsible for the positive eigenvalue. To some extend this term resembles a positon (singular) solution but remains bounded. Our approach is based upon certain limiting arguments and techniques of Hankel operators.


title: Gaussian Process Learning via Fisher Scoring of Vecchia's Approximation

abstract: We derive a single pass algorithm for computing the gradient and Fisher information of Vecchia's Gaussian process loglikelihood approximation, which provides a computationally efficient means for applying the Fisher scoring algorithm for maximizing the loglikelihood. The advantages of the optimization techniques are demonstrated in numerical examples and in an application to Argo ocean temperature data. The new methods are more accurate and much faster than an optimization method that uses only function evaluations, especially when the covariance function has many parameters. This allows practitioners to fit nonstationary models to large spatial and spatial-temporal datasets.


title: Fast algorithm for computing nonlocal operators with finite interaction distance

abstract: Developments of nonlocal operators for modeling processes that traditionally have been described by local differential operators have been increasingly active during the last few years. One example is peridynamics for brittle materials and another is nonstandard diffusion including the use of fractional derivatives. A major obstacle for application of these methods is the high computational cost from the numerical implementation of the nonlocal operators. It is natural to consider fast methods of fast multipole or hierarchical matrix type to overcome this challenge. Unfortunately the relevant kernels do not satisfy the standard necessary conditions. In this work a new class of fast algorithms is developed and analyzed, which is some cases reduces the computational complexity of applying nonlocal operators to essentially the same order of magnitude as the complexity of standard local numerical methods.


title: Mass operator of the M2-brane on a background with constant three-form

abstract: The formulation of supermembrane theory on nontrivial backgrounds is discussed. In particular, we obtain the Hamiltonian of the supermembrane on a background with constant bosonic three form on a target space $M_9 \times T^2$.


title: Word Usage Similarity Estimation with Sentence Representations and Automatic Substitutes

abstract: Usage similarity estimation addresses the semantic proximity of word instances in different contexts. We apply contextualized (ELMo and BERT) word and sentence embeddings to this task, and propose supervised models that leverage these representations for prediction. Our models are further assisted by lexical substitute annotations automatically assigned to word instances by context2vec, a neural model that relies on a bidirectional LSTM. We perform an extensive comparison of existing word and sentence representations on benchmark datasets addressing both graded and binary similarity. The best performing models outperform previous methods in both settings.


title: Non-crossing Brownian paths and Dyson Brownian motion under a moving boundary

abstract: We compute analytically the probability $S(t)$ that a set of $N$ Brownian paths do not cross each other and stay below a moving boundary $g(\tau)= W \sqrt{\tau}$ up to time $t$. We show that for large $t$ it decays as a power law $S(t) \sim t^{- \beta(N,W)}$. The decay exponent $\beta(N,W)$ is obtained as the ground state energy of a quantum system of $N$ non-interacting fermions in a harmonic well in the presence of an infinite hard wall at position $W$. Explicit expressions for $\beta(N,W)$ are obtained in various limits of $N$ and $W$, in particular for large $N$ and large $W$. We obtain the joint distribution of the positions of the walkers in the presence of the moving barrier $g(\tau) =W \sqrt{\tau}$ at large time. We extend our results to the case of $N$ Dyson Brownian motions (corresponding to the Gaussian Unitary Ensemble) in the presence of the same moving boundary $g(\tau)=W\sqrt{\tau}$. For $W=0$ we show that the system provides a realization of a Laguerre biorthogonal ensemble in random matrix theory. We obtain explicitly the average density near the barrier, as well as in the bulk far away from the barrier. Finally we apply our results to $N$ non-crossing Brownian bridges on the interval $[0,T]$ under a time-dependent barrier $g_B(\tau)= W \sqrt{\tau(1- \frac{\tau}{T})}$.

categories: cond-mat.stat-mech cond-mat.dis-nn math-ph math.MP math.PR nlin.SI


title: Eddy diffusivity in homogeneous isotropic turbulence

abstract: We use the recently developed Macroscopic Forcing Method [Mani and Park, Physical Review Fluids, 6:054607, 2021] to compute the scale-dependent eddy diffusivity characterizing ensemble-averaged scalar and momentum transport in incompressible homogeneous isotropic turbulence. For scales larger than the energy containing eddies, eddy diffusivity is found to be constant and consistent with the Boussinesq approximation. However for small scales eddy diffusivity is found to vanish inversely proportional to the wavenumber. Behavior at all scales is reasonably captured by a non-local eddy diffusivity operator modeled as $D/\sqrt{\mathcal{I}-l^2\nabla^2}$, where $D$ is the eddy diffusivity in the Boussinesq limit, and $l$ is a constant on the order of the large-eddy length. These results present the first direct measurement of eddy diffusivity in turbulence with implications in turbulence modeling.


title: Method for solving hyperbolic systems with initial data on non-characteristic manifolds with applications to the shallow water wave equations

abstract: We are concerned with hyperbolic systems of order-one linear PDEs originated on non-characteristic manifolds. We put forward a simple but effective method of transforming such initial conditions to standard initial conditions (i.e. when the solution is specified at an initial moment of time). We then show how our method applies in fluid mechanics. More specifically, we present a complete solution to the problem of long waves run-up in inclined bays of arbitrary shape with nonzero initial velocity.


title: Statistical methods research done as science rather than mathematics

abstract: This paper is about how we study statistical methods. As an example, it uses the random regressions model, in which the intercept and slope of cluster-specific regression lines are modeled as a bivariate random effect. Maximizing this model's restricted likelihood often gives a boundary value for the random effect correlation or variances. We argue that this is a problem; that it is a problem because our discipline has little understanding of how contemporary models and methods map data to inferential summaries; that we lack such understanding, even for models as simple as this, because of a near-exclusive reliance on mathematics as a means of understanding; and that math alone is no longer sufficient. We then argue that as a discipline, we can and should break open our black-box methods by mimicking the five steps that molecular biologists commonly use to break open Nature's black boxes: design a simple model system, formulate hypotheses using that system, test them in experiments on that system, iterate as needed to reformulate and test hypotheses, and finally test the results in an "in vivo" system. We demonstrate this by identifying conditions under which the random-regressions restricted likelihood is likely to be maximized at a boundary value. Resistance to this approach seems to arise from a view that it lacks the certainty or intellectual heft of mathematics, perhaps because simulation experiments in our literature rarely do more than measure a new method's operating characteristics in a small range of situations. We argue that such work can make useful contributions including, as in molecular biology, the findings themselves and sometimes the designs used in the five steps; that these contributions have as much practical value as mathematical results; and that therefore they merit publication as much as the mathematical results our discipline esteems so highly.


title: The Proximity Effect in a Superconductor-Quasicrystal Hybrid Ring

abstract: We compute the real-space profile of the superconducting order parameter (OP) in a hybrid ring that consists of a 1D superconductor connected to a Fibonacci chain using a self-consistent approach. In our study, the strength of the penetration, as measured by the order parameter at the center of the quasicrystal, depends on the structural parameter $\phi$, or phason angle, that characterizes different realizations of the Fibonacci chains of a given length. We show that the penetration strength dependence on $\phi$ reflects properties of the topological edge states of the Fibonacci chain. We show that the induced superconducting order parameter averaged over all chains has a power law decay as a function of distance from the S-N interface. More interestingly, we show that there are large OP fluctuations for individual chains and that the penetration strength in a finite Fibonacci chain can be significantly $larger$ than in a normal periodic conductor for special values of $\phi$.


title: Short-depth circuits for efficient expectation value estimation

abstract: The evaluation of expectation values $Tr\left[\rho O\right]$ for some pure state $\rho$ and Hermitian operator $O$ is of central importance in a variety of quantum algorithms. Near optimal techniques developed in the past require a number of measurements $N$ approaching the Heisenberg limit $N=\mathcal{O}\left(1/\epsilon\right)$ as a function of target accuracy $\epsilon$. The use of Quantum Phase Estimation requires however long circuit depths $C=\mathcal{O}\left(1/\epsilon\right)$ making their implementation difficult on near term noisy devices. The more direct strategy of Operator Averaging is usually preferred as it can be performed using $N=\mathcal{O}\left(1/\epsilon^2\right)$ measurements and no additional gates besides those needed for the state preparation. In this work we use a simple but realistic model to describe the bound state of a neutron and a proton (the deuteron) and show that the latter strategy can require an overly large number of measurement in order to achieve a reasonably small relative target accuracy $\epsilon_r$. We propose to overcome this problem using a single step of QPE and classical post-processing. This approach leads to a circuit depth $C=\mathcal{O}\left(\epsilon^\mu\right)$ (with $\mu\geq0$) and to a number of measurements $N=\mathcal{O}\left(1/\epsilon^{2+\nu}\right)$ for $0<\nu\leq1$. We provide detailed descriptions of two implementations of our strategy for $\nu=1$ and $\nu\approx0.5$ and derive appropriate conditions that a particular problem instance has to satisfy in order for our method to provide an advantage.


title: Some ground-state expectation values for the free parafermion Z(N) spin chain

abstract: We consider the calculation of ground-state expectation values for the non-Hermitian Z(N) spin chain described by free parafermions. For N=2 the model reduces to the quantum Ising chain in a transverse field with open boundary conditions. Use is made of the Hellmann-Feynman theorem to obtain exact results for particular single site and nearest-neighbour ground-state expectation values for general N which are valid for sites deep inside the chain. These results are tested numerically for N=3, along with how they change as a function of distance from the boundary.

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


title: Pseudoscalar flavor-singlet mesons from lattice QCD

abstract: We investigate the masses and mixing of $\eta$, $\eta'$ mesons in the framework of twisted mass lattice QCD with $N_f=2+1+1$ dynamical quark flavors. For the first time we perform a controlled chiral and continuum extrapolation to obtain physical results. To this end, we have analyzed 17 gauge ensembles with pion masses ranging from $220\,\mathrm{MeV}$ to $500\,\mathrm{MeV}$ and at three values of the lattice spacing. This calculation became feasible through the application of a powerful variance reduction technique for the computation of quark-disconnected diagrams together with several improvements in the final analysis, including a correction of the relevant correlation functions for residual topological finite volume effects. We obtain physical results for the masses $M_{\eta}=557(11)_\mathrm{stat}(03)_{\chi\mathrm{PT}}\,\mathrm{MeV}$ and $M_{\eta'}=911(64)_\mathrm{stat}(03)_{\chi\mathrm{PT}}\,\mathrm{MeV}$, as well as the mixing angle in the quark flavor basis $\phi=38.8(2.2)_\mathrm{stat}(2.4)_{\chi\mathrm{PT}}^\circ$, in excellent agreement with other results from phenomenology. Similarly, for the physical values of the decay constant parameters we find $f_l=125(5)_\mathrm{stat}(6)_{\chi\mathrm{PT}}\,\mathrm{MeV}$ and $f_s=178(4)_\mathrm{stat}(1)_{\chi\mathrm{PT}}\,\mathrm{MeV}$. Finally, we present a test of the Veneziano Witten formula using our lattice data.


title: Scylla: A Mesos Framework for Container Based MPI Jobs

abstract: Open source cloud technologies provide a wide range of support for creating customized compute node clusters to schedule tasks and managing resources. In cloud infrastructures such as Jetstream and Chameleon, which are used for scientific research, users receive complete control of the Virtual Machines (VM) that are allocated to them. Importantly, users get root access to the VMs. This provides an opportunity for HPC users to experiment with new resource management technologies such as Apache Mesos that have proven scalability, flexibility, and fault tolerance. To ease the development and deployment of HPC tools on the cloud, the containerization technology has matured and is gaining interest in the scientific community. In particular, several well known scientific code bases now have publicly available Docker containers. While Mesos provides support for Docker containers to execute individually, it does not provide support for container inter-communication or orchestration of the containers for a parallel or distributed application. In this paper, we present the design, implementation, and performance analysis of a Mesos framework, Scylla, which integrates Mesos with Docker Swarm to enable orchestration of MPI jobs on a cluster of VMs acquired from the Chameleon cloud [1]. Scylla uses Docker Swarm for communication between containerized tasks (MPI processes) and Apache Mesos for resource pooling and allocation. Scylla allows a policy-driven approach to determine how the containers should be distributed across the nodes depending on the CPU, memory, and network throughput requirement for each application.