title: Robot-Supervised Learning for Object Segmentation

abstract: To be effective in unstructured and changing environments, robots must learn to recognize new objects. Deep learning has enabled rapid progress for object detection and segmentation in computer vision; however, this progress comes at the price of human annotators labeling many training examples. This paper addresses the problem of extending learning-based segmentation methods to robotics applications where annotated training data is not available. Our method enables pixelwise segmentation of grasped objects. We factor the problem of segmenting the object from the background into two sub-problems: (1) segmenting the robot manipulator and object from the background and (2) segmenting the object from the manipulator. We propose a kinematics-based foreground segmentation technique to solve (1). To solve (2), we train a self-recognition network that segments the robot manipulator. We train this network without human supervision, leveraging our foreground segmentation technique from (1) to label a training set of images containing the robot manipulator without a grasped object. We demonstrate experimentally that our method outperforms state-of-the-art adaptable in-hand object segmentation. We also show that a training set composed of automatically labelled images of grasped objects improves segmentation performance on a test set of images of the same objects in the environment.


title: On a generalization of the Howe-Moore property

abstract: We define a Howe-Moore property relative to a set of subgroups. Namely, a group $G$ has the Howe-Moore property relative to a set $\mathcal{F}$ of subgroups if for every unitary representation $\pi$ of $G$, whenever the restriction of $\pi$ to any element of $\mathcal{F}$ has no non-trivial invariant vectors, the matrix coefficients vanish at infinity. We prove that a semisimple group has the Howe-Moore property relatively to the family of its factors.


title: Some variations on Lyndon words

abstract: In this paper we compare two finite words $u$ and $v$ by the lexicographical order of the infinite words $u^\omega$ and $v^\omega$. Informally, we say that we compare $u$ and $v$ by the infinite order. We show several properties of Lyndon words expressed using this infinite order. The innovative aspect of this approach is that it allows to take into account also non trivial conditions on the prefixes of a word, instead that only on the suffixes. In particular, we derive a result of Ufnarovskij [V. Ufnarovskij, "Combinatorial and asymptotic methods in algebra", 1995] that characterizes a Lyndon word as a word which is greater, with respect to the infinite order, than all its prefixes. Motivated by this result, we introduce the prefix standard permutation of a Lyndon word and the corresponding (left) Cartesian tree. We prove that the left Cartesian tree is equal to the left Lyndon tree, defined by the left standard factorization of Viennot [G. Viennot, "Alg\`ebres de Lie libres et mono\"ides libres", 1978]. This result is dual with respect to a theorem of Hohlweg and Reutenauer [C. Hohlweg and C. Reutenauer, "Lyndon words, permutations and trees", 2003].


title: Frobenius and Homological Dimensions of Complexes

abstract: It is proved that a module $M$ over a Noetherian local ring $R$ of prime characteristic and positive dimension has finite flat dimension if Tor$_i^R({}^e R, M)=0$ for dim $R$ consecutive positive values of $i$ and infinitely many $e$. Here ${}^e R$ denotes the ring $R$ viewed as an $R$-module via the $e$th iteration of the Frobenius endomorphism. In the case $R$ is Cohen-Macualay, it suffices that the Tor vanishing above holds for a single $e\geq \log_p e(R)$, where $e(R)$ is the multiplicity of the ring. This improves a result of D. Dailey, S. Iyengar, and the second author, as well as generalizing a theorem due to C. Miller from finitely generated modules to arbitrary modules. We also show that if $R$ is a complete intersection ring then the vanishing of Tor$_i^R({}^e R, M)$ for single positive values of $i$ and $e$ is sufficient to imply $M$ has finite flat dimension. This extends a result of L. Avramov and C. Miller.


title: Guided Meta-Policy Search

abstract: Reinforcement learning (RL) algorithms have demonstrated promising results on complex tasks, yet often require impractical numbers of samples since they learn from scratch. Meta-RL aims to address this challenge by leveraging experience from previous tasks so as to more quickly solve new tasks. However, in practice, these algorithms generally also require large amounts of on-policy experience during the meta-training process, making them impractical for use in many problems. To this end, we propose to learn a reinforcement learning procedure in a federated way, where individual off-policy learners can solve the individual meta-training tasks, and then consolidate these solutions into a single meta-learner. Since the central meta-learner learns by imitating the solutions to the individual tasks, it can accommodate either the standard meta-RL problem setting or a hybrid setting where some or all tasks are provided with example demonstrations. The former results in an approach that can leverage policies learned for previous tasks without significant amounts of on-policy data during meta-training, whereas the latter is particularly useful in cases where demonstrations are easy for a person to provide. Across a number of continuous control meta-RL problems, we demonstrate significant improvements in meta-RL sample efficiency in comparison to prior work as well as the ability to scale to domains with visual observations.


title: Recursive Green's function approach to Feenberg perturbation theory

abstract: We propose a new procedure by using the recursive Green's functions which remove all the repetition terms from the time-independent perturbation series for finite-level quantum systems. These Green's functions are introduced as a generalization of the Brillouin-Wigner perturbation theory and the calculations of their diagonal elements can naturally give the effective propagators which are equal to the ones in the Feenberg perturbation theory.


title: On the Efficiency of the Iterative Techniques for Solving Incompressible Navier-Stokes Equations

abstract: It is well known that the choice of the iterative method is crucial in determining the speed of the converged solution. This article presents a detailed comparison between several iterative techniques for solving incmopressible Navier-Stokes equations. The numerical approaches implemented in the solver include Jacobi, Gauss-Siedel, Successive Over Relaxation, Alternating Direct Implicit and Multigrid methods. The results reveal that multigrid method is the most powerful iterative method among all other methods investigated in terms of the computational time and the number of iterations.


title: Novel scheme for parametrizing the chemical freeze-out surface in Heavy Ion Collision Experiments

abstract: We introduce a new prescription for obtaining the chemical freeze-out parameters in the heavy-ion collision experiments using the Hadron Resonance Gas model. The scheme is found to reliably estimate the freeze-out parameters and predict the hadron yield ratios, which themselves were never used in the parametrization procedure.


title: A characterization of 3D steady Euler flows using commuting zero-flux homologies

abstract: We characterize, using commuting zero-flux homologies, those volume-preserving vector fields on a $3$-manifold that are steady solutions of the Euler equations for some Riemannian metric. This result extends Sullivan's homological characterization of geodesible flows in the volume-preserving case. As an application, we show that the steady Euler flows cannot be constructed using plugs (as in Wilson's or Kuperberg's constructions). Analogous results in higher dimensions are also proved.


title: Giant Casimir torque between rotated gratings and the $\theta=0$ anomaly

abstract: We study the Casimir torque between two metallic one-dimensional gratings rotated by an angle $\theta$ with respect to each other. We find that, for infinitely extended gratings, the Casimir energy is anomalously discontinuous at $\theta=0$, due to a critical zero-order geometric transition between a 2D- and a 1D-periodic system. This transition is a peculiarity of the grating geometry and does not exist for intrinsically anisotropic materials. As a remarkable practical consequence, for finite-size gratings, the torque per area can reach extremely large values, increasing without bounds with the size of the system. We show that for finite gratings with only 10 period repetitions, the maximum torque is already 60 times larger than the one predicted in the case of infinite gratings. These findings pave the way to the design of a contactless quantum vacuum torsional spring, with possible relevance to micro- and nano-mechanical devices.


title: Large Batch Optimization for Deep Learning: Training BERT in 76 minutes

abstract: Training large deep neural networks on massive datasets is computationally very challenging. There has been recent surge in interest in using large batch stochastic optimization methods to tackle this issue. The most prominent algorithm in this line of research is LARS, which by employing layerwise adaptive learning rates trains ResNet on ImageNet in a few minutes. However, LARS performs poorly for attention models like BERT, indicating that its performance gains are not consistent across tasks. In this paper, we first study a principled layerwise adaptation strategy to accelerate training of deep neural networks using large mini-batches. Using this strategy, we develop a new layerwise adaptive large batch optimization technique called LAMB; we then provide convergence analysis of LAMB as well as LARS, showing convergence to a stationary point in general nonconvex settings. Our empirical results demonstrate the superior performance of LAMB across various tasks such as BERT and ResNet-50 training with very little hyperparameter tuning. In particular, for BERT training, our optimizer enables use of very large batch sizes of 32868 without any degradation of performance. By increasing the batch size to the memory limit of a TPUv3 Pod, BERT training time can be reduced from 3 days to just 76 minutes (Table 1). The LAMB implementation is available at https://github.com/tensorflow/addons/blob/master/tensorflow_addons/optimizers/lamb.py


title: Well-posedness of cubic Horndeski theories

abstract: We study the local well-posedness of the initial value problem for cubic Horndeski theories. Three different strongly hyperbolic modifications of the ADM formulation of the Einstein equations are extended to cubic Horndeski theories in the weak field regime. In the first one, the equations of motion are rewritten as a coupled elliptic-hyperbolic system of partial differential equations. The second one is based on the BSSN formulation with a generalised Bona-Mass\'o slicing (covering the 1+log slicing) and non-dynamical shift vector. The third one is an extension of the CCZ4 formulation with a generalised Bona-Mass\'o slicing (also covering the 1+log slicing) and a gamma driver shift condition. This final formulation may be particularly suitable for applications in non-linear numerical simulations.


title: Complexity and Algorithms for Semipaired Domination in Graphs

abstract: For a graph $G=(V,E)$ with no isolated vertices, a set $D\subseteq V$ is called a semipaired dominating set of G if $(i)$ $D$ is a dominating set of $G$, and $(ii)$ $D$ can be partitioned into two element subsets such that the vertices in each two element set are at distance at most two. The minimum cardinality of a semipaired dominating set of $G$ is called the semipaired domination number of $G$, and is denoted by $\gamma_{pr2}(G)$. The \textsc{Minimum Semipaired Domination} problem is to find a semipaired dominating set of $G$ of cardinality $\gamma_{pr2}(G)$. In this paper, we initiate the algorithmic study of the \textsc{Minimum Semipaired Domination} problem. We show that the decision version of the \textsc{Minimum Semipaired Domination} problem is NP-complete for bipartite graphs and split graphs. On the positive side, we present a linear-time algorithm to compute a minimum cardinality semipaired dominating set of interval graphs and trees. We also propose a $1+\ln(2\Delta+2)$-approximation algorithm for the \textsc{Minimum Semipaired Domination} problem, where $\Delta$ denote the maximum degree of the graph and show that the \textsc{Minimum Semipaired Domination} problem cannot be approximated within $(1-\epsilon) \ln|V|$ for any $\epsilon > 0$ unless NP $\subseteq$ DTIME$(|V|^{O(\log\log|V|)})$.


title: Octagons I: Combinatorics and Non-Planar Resummations

abstract: We explain how the 't Hooft expansion of correlators of half-BPS operators can be resummed in a large-charge limit in N=4 super Yang-Mills theory. The full correlator in the limit is given by a non-trivial function of two variables: One variable is the charge of the BPS operators divided by the square root of the number Nc of colors; the other variable is the octagon that contains all the 't Hooft coupling and spacetime dependence. At each genus g in the large Nc expansion, this function is a polynomial of degree 2g+2 in the octagon. We find several dual matrix model representations of the correlators in the large-charge limit. Amusingly, the number of colors in these matrix models is formally taken to zero in the relevant limit.


title: Local-global principles for norm one tori over semi-global fields

abstract: Let K be a complete discretely valued field with residue field k and F be a function field of a curve over K. Let L/F be a Galois extension of degree n. If n is coprime to char(k), then under some assumptions on k(e.g. k is algebraically closed or a finite field) and on the geometry of the curve, we show that there is a local-global principle with respect to discrete valuations for the norms from the extension L/F.


title: Synchronization of Network-Coupled Oscillators with Uncertain Dynamics

abstract: Synchronization of network-coupled dynamical units is important to a variety of natural and engineered processes including circadian rhythms, cardiac function, neural processing, and power grids. Despite this ubiquity, it remains poorly understood how complex network structures and heterogeneous local dynamics combine to either promote or inhibit synchronization. Moreover, for most real-world applications it is impossible to obtain the exact specifications of the system, and there is a lack of theory for how uncertainty affects synchronization. We address this open problem by studying the Synchrony Alignment Function (SAF), which is an objective measure for the synchronization properties of a network of heterogeneous oscillators with given natural frequencies. We extend the SAF framework to analyze network-coupled oscillators with heterogeneous natural frequencies that are drawn as a multivariate random vector. Using probability theory for quadratic forms, we obtain expressions for the expectation and variance of the SAF for given network structures. We conclude with numerical experiments that illustrate how the incorporation of uncertainty yields a more robust theoretical framework for enhancing synchronization, and we provide new perspectives for why synchronization is generically promoted by network properties including degree-frequency correlations, link directedness, and link weight delocalization.


title: An all-electron product-basis set: application to plasmon anisotropy in simple metals

abstract: An efficient basis set for products of all-electron wave functions is proposed, which comprises plane waves defined over the entire unit cell and orbitals confined to small non-overlapping spheres. The size of the set and the basis functions are, in principle, independent of the computational parameters of the band structure method. The approach is implemented in the extended LAPW method, and its properties and accuracy are discussed. The method is applied to analyze the dielectric response of the simple metals Al, Na, Li, K, Rb, and Cs with a focus on the origin of the anisotropy of the plasmon dispersion in Al and Na. The anisotropy is traced to tiny structures of the one-particle excitation spectra of Al and Na, and relevant experimental observations are explained.


title: Free quotients of fundamental groups of smooth quasi-projective varieties

abstract: We consider the structure of classes of curves on a projective simply connected surface for which fundamental groups of the complements admit free quotients having rank greater than one with irreducible components belonging to a selected subset the effective cone of the surface. In particular we show a finiteness result for such classes if the ranks of free quotients of the fundamental groups with components in the subset of effective cone are sufficiently large.


title: The Ineludible non-Gaussianity of the Primordial Black Hole Abundance

abstract: We study the formation of primordial black holes when they are generated by the collapse of large overdensities in the early universe. Since the density contrast is related to the comoving curvature perturbation by a nonlinear relation, the overdensity statistics is unavoidably non-Gaussian. We show that the abundance of primordial black holes at formation may not be captured by a perturbative approach which retains the first few cumulants of the non-Gaussian probability distribution. We provide two techniques to calculate the non-Gaussian abundance of primordial black holes at formation, one based on peak theory and the other on threshold statistics. Our results show that the unavoidable non-Gaussian nature of the inhomogeneities in the energy density makes it harder to generate PBHs. We provide simple (semi-)analytical expressions to calculate the non-Gaussian abundances of the primordial black holes and show that for both narrow and broad power spectra the gaussian case from threshold statistics is reproduced by increasing the amplitude of the power spectrum by a factor ${\cal O}(2\div 3)$.


title: Asymptotic sign coherence conjecture

abstract: The sign coherence phenomenon is an important feature of c-vectors in cluster algebras with principal coefficients. In this note, we consider a more general version of c-vectors defined for arbitrary cluster algebras of geometric type and formulate a conjecture describing their asymptotic behavior. This conjecture, which is called the asymptotic sign coherence conjecture, states that for any infinite sequence of matrix mutations that satisfies certain natural conditions, the corresponding c-vectors eventually become sign coherent. We prove this conjecture for rank 2 cluster algebras of infinite type and for a particular sequence of mutations in a cluster algebra associated with the Markov quiver.


title: Formulation of Entropy-Stable schemes for the multicomponent compressible Euler equations

abstract: In this work, Entropy-Stable (ES) schemes are formulated for the multicomponent compressible Euler equations. Entropy-conservative (EC) and ES fluxes are derived. Particular attention is paid to the limit case of zero partial densities where the structure required by ES schemes breaks down (the entropy variables are no longer defined). It is shown that while an EC flux is well-defined in this limit, a well-defined upwind ES flux requires appropriately averaged partial densities in the dissipation matrix. A similar result holds for the high-order TecNO reconstruction. However, this does not prevent the numerical solution from developing negative partial densities or internal energy. Numerical experiments were performed on one-dimensional and two-dimensional interface and shock-interface problems. The present scheme exactly preserves stationary interfaces. On moving interfaces, it produces pressure anomalies typically observed with conservative schemes [Karni, \textit{J. Comput. Phys.}, 112 (1994) 1]. We find that these anomalies, which are not present in the single-component case, violate neither entropy stability nor a minimum principle of the specific entropy. Finally, we show that the scheme is able to reproduce the physical mechanisms of the two-dimensional shock-bubble interaction problem [Haas \& Sturtevant, J. Fluid Mech. 181 (1987) 41, Quirk \& Karni, J. Fluid Mech. 318 (1996) 129].


title: Recognising and evaluating the effectiveness of extortion in the Iterated Prisoner's Dilemma

abstract: Since the introduction of zero-determinant strategies, extortionate strategies have received considerable interest. While an interesting class of strategies, the definitions of extortionate strategies are algebraically rigid, apply only to memory-one strategies, and require complete knowledge of a strategy (memory-one cooperation probabilities). We describe a method to detect extortionate behaviour from the history of play of a strategy. When applied to a corpus of 204 strategies this method detects extortionate behaviour in well-known extortionate strategies as well others that do not fit the algebraic definition. The highest performing strategies in this corpus are able to exhibit selectively extortionate behavior, cooperating with strong strategies while exploiting weaker strategies, which no memory-one strategy can do. These strategies emerged from an evolutionary selection process and their existence contradicts widely-repeated folklore in the evolutionary game theory literature: complex strategies can be extraordinarily effective, zero-determinant strategies can be outperformed by non-zero determinant strategies, and longer memory strategies are able to outperform short memory strategies. Moreover, while resistance to extortion is critical for the evolution of cooperation, the extortion of weak opponents need not prevent cooperation between stronger opponents, and this adaptability may be crucial to maintaining cooperation in the long run.


title: Halo breakup and the Coulomb-nuclear interference problem

abstract: The Coulomb-nuclear interference is studied as a function of the projectile ground-state binding energy in the 8Li + 12C and 8Li + 208Pb breakup reactions, by considering an arbitrary range for the 8Li ground-state binding energies {\epsilon}b, varying from the experimental one 2.03 MeV down to 0.01 MeV. Regardless the target mass, we first show that the Coulomb breakup cross section is stronger dependent on {\epsilon}b than the nuclear breakup cross section, due to the long-range nature of the Coulomb forces and to the electromagnetic transition matrix elements. For example, in case of 8Li + 208Pb reaction at Elab = 60MeV, it is found that |{\sigma}_int|\simeq 4\times{\sigma}_nucl$, while ${\sigma}_Coul\simeq 35\times{\sigma_nucl}$. This shows clearly that small nuclear contribution in a Coulomb-dominated reaction does not imply insignificant Coulomb-nuclear interference. Such result can be mainly attributed to peripheral interference phenomenon, represented by a function of the binding energy, which determines the peripheral range of nuclear forces, where Coulomb and nuclear forces strongly interfere destructively.


title: Quantum fluctuations stabilize an inverted pendulum

abstract: We explore analytically the quantum dynamics of a point mass pendulum using the Heisenberg equation of motion. Choosing as variables the mean position of the pendulum, a suitably defined generalised variance and a generalised skewness, we set up a dynamical system which reproduces the correct limits of simple harmonic oscillator like and free rotor like behaviour. We then find the unexpected result that the quantum pendulum released from and near the inverted position executes oscillatory motion around the classically unstable position provided the initial wave packet has a variance much greater than the variance of the well known coherent state of the simple harmonic oscillator. The behaviour of the dynamical system for the quantum pendulum is a higher dimensional analogue of the behaviour of the Kapitza pendulum where the point of support is vibrated vertically with a frequency higher than the critical value needed to stabilize the inverted position. A somewhat similar phenomenon has recently been observed in the non equilibrium dynamics of a spin - 1 Bose-Einstein Condensate.


title: Bisimulation for Feller-Dynkin Processes

abstract: Bisimulation is a concept that captures behavioural equivalence. It has been studied extensively on nonprobabilistic systems and on discrete-time Markov processes and on so-called continuous-time Markov chains. In the latter time is continuous but the evolution still proceeds in jumps. We propose two definitions of bisimulation on continuous-time stochastic processes where the evolution is a \emph{flow} through time. We show that they are equivalent and we show that when restricted to discrete-time, our concept of bisimulation encompasses the standard discrete-time concept. The concept we introduce is not a straightforward generalization of discrete-time concepts.


title: Sentiment analysis with genetically evolved Gaussian kernels

abstract: Sentiment analysis consists of evaluating opinions or statements from the analysis of text. Among the methods used to estimate the degree in which a text expresses a given sentiment, are those based on Gaussian Processes. However, traditional Gaussian Processes methods use a predefined kernel with hyperparameters that can be tuned but whose structure can not be adapted. In this paper, we propose the application of Genetic Programming for evolving Gaussian Process kernels that are more precise for sentiment analysis. We use use a very flexible representation of kernels combined with a multi-objective approach that simultaneously considers two quality metrics and the computational time spent by the kernels. Our results show that the algorithm can outperform Gaussian Processes with traditional kernels for some of the sentiment analysis tasks considered.


title: Complete Matching for Quasi-distribution Functions in Large Momentum Effective Theory

abstract: We complete the procedure of extracting parton distribution functions (PDFs) using large momentum effective theory (LaMET) at leading power accuracy in the hadron momentum. We derive a general factorization formula for the quasi PDFs in the presence of mixing, and give the corresponding hard matching kernel at $\mathcal O(\alpha_s)$, both for the unpolarized and for the polarized quark and gluon quasi-PDFs. Our calculation is performed in a regularization-independent momentum subtraction scheme. The results allow us to match the nonperturbatively renormalized quasi-PDFs to normal PDFs in the presence of mixing, and therefore can be used to extract flavor-singlet quark PDFs as well as gluon PDFs from lattice simulations.


title: Regional Homogeneity: Towards Learning Transferable Universal Adversarial Perturbations Against Defenses

abstract: This paper focuses on learning transferable adversarial examples specifically against defense models (models to defense adversarial attacks). In particular, we show that a simple universal perturbation can fool a series of state-of-the-art defenses. Adversarial examples generated by existing attacks are generally hard to transfer to defense models. We observe the property of regional homogeneity in adversarial perturbations and suggest that the defenses are less robust to regionally homogeneous perturbations. Therefore, we propose an effective transforming paradigm and a customized gradient transformer module to transform existing perturbations into regionally homogeneous ones. Without explicitly forcing the perturbations to be universal, we observe that a well-trained gradient transformer module tends to output input-independent gradients (hence universal) benefiting from the under-fitting phenomenon. Thorough experiments demonstrate that our work significantly outperforms the prior art attacking algorithms (either image-dependent or universal ones) by an average improvement of 14.0% when attacking 9 defenses in the transfer-based attack setting. In addition to the cross-model transferability, we also verify that regionally homogeneous perturbations can well transfer across different vision tasks (attacking with the semantic segmentation task and testing on the object detection task). The code is available here: https://github.com/LiYingwei/Regional-Homogeneity.


title: Earthquake Arrival Association with Backprojection and Graph Theory

abstract: The association of seismic wave arrivals with causative earthquakes becomes progressively more challenging as arrival detection methods become more sensitive, and particularly when earthquake rates are high. For instance, seismic waves arriving across a monitoring network from several sources may overlap in time, false arrivals may be detected, and some arrivals may be of unknown phase (e.g., P- or S-waves). We propose an automated method to associate arrivals with earthquake sources and obtain source locations applicable to such situations. To do so we use a pattern detection metric based on the principle of backprojection to reveal candidate sources, followed by graph-theory-based clustering and an integer linear optimization routine to associate arrivals with the minimum number of sources necessary to explain the data. This method solves for all sources and phase assignments simultaneously, rather than in a sequential greedy procedure as is common in other association routines. We demonstrate our method on both synthetic and real data from the Integrated Plate Boundary Observatory Chile (IPOC) seismic network of northern Chile. For the synthetic tests we report results for cases with varying complexity, including rates of 500 earthquakes/day and 500 false arrivals/station/day, for which we measure true positive detection accuracy of > 95%. For the real data we develop a new catalog between January 1, 2010 - December 31, 2017 containing 817,548 earthquakes, with detection rates on average 279 earthquakes/day, and a magnitude-of-completion of ~M1.8. A subset of detections are identified as sources related to quarry and industrial site activity, and we also detect thousands of foreshocks and aftershocks of the April 1, 2014 Mw 8.2 Iquique earthquake. During the highest rates of aftershock activity, > 600 earthquakes/day are detected in the vicinity of the Iquique earthquake rupture zone.


title: Toric degenerations of Grassmannians and Schubert varieties from matching field tableaux

abstract: We study the combinatorics of Gr\"obner degenerations of Grassmannians and the Schubert varieties inside them. We provide a family of binomial ideals whose combinatorics is governed by tableaux induced by matching fields in the sense of Sturmfels and Zelevinsky. We prove that these ideals are all quadratically generated and they yield a SAGBI basis of the Pl\"ucker algebra. This leads to a new family of toric degenerations of Grassmannians. Moreover, we apply our results to construct a family of Gr\"obner degenerations of Schubert varieties inside Grassmannians. We provide a complete characterization of toric ideals among these degenerations in terms of the combinatorics of matching fields, permutations, and semi-standard tableaux.


title: Automatic Nonrigid Histological Image Registration with Adaptive Multistep Algorithm

abstract: In this paper, we present a short description of the method proposed to ANHIR challenge organized jointly with the IEEE ISBI 2019 conference. We propose a method consisting of preprocessing, initial alignment, nonrigid registration algorithms and a method to automatically choose the best result. The method turned out to be robust (99.792% robustness) and accurate (0.38% average median rTRE). The main drawback of the proposed method is relatively high computation time. However, this aspect can be easily improved by cleaning the code and proposing a GPU implementation.


title: Unitary equivalence of operator-valued multishifts

abstract: We systematically study various aspects of operator-valued multishifts. Beginning with basic properties, we show that the class of multishifts on the directed Cartesian product of rooted directed trees is contained in that of operator-valued multishifts. Further, we establish circularity, analyticity and wandering subspace property of these multishifts. In the rest part of the paper, we study the function theoretic behaviour of operator-valued multishifts. We determine the bounded point evaluation, reproducing kernel structure and the unitary equivalence of operator-valued multishifts with invertible operator weights. In contrast with a result of Lubin, it appears that the set of all bounded point evaluations of an operator-valued multishift may be properly contained in the joint point spectrum of the adjoint of underlying multishift.


title: Primordial black hole formation and abundance: contribution from the non-linear relation between the density and curvature perturbation

abstract: The formation and abundance of primordial black holes (PBHs) arising from the curvature perturbation $\zeta$ is studied. The non-linear relation between $\zeta$ and the density contrast $\delta$ means that, even when $\zeta$ has an exactly Gaussian distribution, significant non-Gaussianities affecting PBH formation must be considered. Numerical simulations are used to investigate the critical value and the mass of PBHs which form, and peaks theory is used to calculate the mass fraction of the universe collapsing to form PBHs at the time of formation. A formalism to calculate the total present day PBH abundance and mass function is also derived. It is found that the abundance of PBHs is very sensitive to the non-linear effects, and that the power spectrum $\mathcal{P}_\zeta$ must be a factor of $\sim2$ larger to produce the same number of PBHs as the linear model (where the exact value depends on the critical value for a region to collapse and form a PBH). This also means that the derived constraints on the small-scale power spectrum from constraints on the abundance of PBHs are weaker by the same factor.


title: Alignment of a circumbinary disc around an eccentric binary with application to KH 15D

abstract: We analyse the evolution of a mildly inclined circumbinary disc that orbits an eccentric orbit binary by means of smoother particle hydrodynamic (SPH) simulations and linear theory. We show that the alignment process of an initially misaligned circumbinary disc around an eccentric orbit binary is significantly different than around a circular orbit binary and involves tilt oscillations. The more eccentric the binary, the larger the tilt oscillations and the longer it takes to damp these oscillations. A circumbinary disc that is only mildly inclined may increase its inclination by a factor of a few before it moves towards alignment. The results of the SPH simulations agree well with those of linear theory. We investigate the properties of the circumbinary disc/ring around KH 15D. We determine disc properties based on the observational constraints imposed by the changing binary brightness. We find that the inclination is currently at a local minimum and will increase substantially before setting to coplanarity. In addition, the nodal precession is currently near its most rapid rate. The recent observations that show a reappearance of Star B impose constraints on the thickness of the layer of obscuring material. Our results suggest that disc solids have undergone substantial inward drift and settling towards to disc midplane. For disc masses $\sim 0.001 M_\odot$, our model indicates that the level of disc turbulence is low $\alpha \ll 0.001$. Another possibility is that the disc/ring contains little gas.


title: A simple method for detecting chaos in nature

abstract: Chaos, or exponential sensitivity to small perturbations, appears everywhere in nature. Moreover, chaos is predicted to play diverse functional roles in living systems. A method for detecting chaos from empirical measurements should therefore be a key component of the biologist's toolkit. But, classic chaos-detection tools are highly sensitive to measurement noise and break down for common edge cases, making it difficult to detect chaos in domains, like biology, where measurements are noisy. However, newer tools promise to overcome these limitations. Here, we combine several such tools into an automated processing pipeline, and show that our pipeline can detect the presence (or absence) of chaos in noisy recordings, even for difficult edge cases. As a first-pass application of our pipeline, we show that heart rate variability is not chaotic as some have proposed, and instead reflects a stochastic process in both health and disease. Our tool is easy-to-use and freely available.


title: Convexity and monotonicity in nonlinear optimal control under uncertainty

abstract: We consider the problem of finite-horizon optimal control design under uncertainty for imperfectly observed discrete-time systems with convex costs and constraints. It is known that this problem can be cast as an infinite-dimensional convex program when the dynamics and measurements are linear, uncertainty is additive, and the risks associated with constraint violations and excessive costs are measured in expectation or in the worst case. In this paper, we extend this result to systems with convex or concave dynamics, nonlinear measurements, more general uncertainty structures and other coherent risk measures. In this setting, the optimal control problem can be cast as an infinite-dimensional convex program if (1) the costs, constraints and dynamics satisfy certain monotonicity properties, and (2) the measured outputs can be reversibly `purified' of the influence of the control inputs through Q- or Youla-parameterization. The practical value of this result is that the finite-dimensional subproblems arising in a variety of suboptimal control methods, notably including model predictive control and the Q-design procedure, are also convex for this class of nonlinear systems. Subproblems can therefore be solved to global optimality using convenient modeling software and efficient, reliable solvers. We illustrate these ideas in a numerical example.


title: High-strength cellulose-polyacrylamide hydrogels: mechanical behavior and structure depending on the type of cellulose

abstract: Two types of high-strength composite hydrogels possessing the structure of interpenetrating polymer networks were synthesized via free-radical polymerization of acrylamide carried out straight within the matrix of plant or bacterial cellulose swollen in the reactive solution. The mechanical behavior of synthesized hydrogels subjected to the action of compressive deformations with different amplitude values was studied. The analysis of the stress-strain curves of compression tests of the hydrogels of both types obtained in different test conditions demonstrates the substantial difference in their mechanical behavior. Both the plant cellulose-based and bacterial cellulose-based hydrogels withstand successfully the compression with the amplitude up to 80 %, but the bacterial cellulose-based compositions demonstrate in the multiple compression tests the substantial decrease of the mechanical stiffness caused by the action of compressions. This effect takes place straightly during the multiple cyclic compression tests or in the course of relaxation of previously compressed samples in water in unloaded state during 2-3 days after the first compressive tests. Being subjected to the action of the extremely sever compression (up to 80 %) the bacterial cellulose-based hydrogel varies greatly its mechanical behavior: during the subsequent compressions the stiffness of the material keeps extremely small (as compared to that registered in the first compression) up to the deformations as high as 50-60 % with the dramatic increase of the stiffness during the further loading. Submicron- and micron-scale specific features of structures of composite hydrogels of both types were studied by cryo-scanning electron microscopy to explain the peculiarities of observed mechanical effects.


title: Counterfactual Sensitivity and Robustness

abstract: We propose a framework for analyzing the sensitivity of counterfactuals to parametric assumptions about the distribution of latent variables in structural models. In particular, we derive bounds on counterfactuals as the distribution of latent variables spans nonparametric neighborhoods of a given parametric specification while other "structural" features of the model are maintained. Our approach recasts the infinite-dimensional problem of optimizing the counterfactual with respect to the distribution of latent variables (subject to model constraints) as a finite-dimensional convex program. We also develop an MPEC version of our method to further simplify computation in models with endogenous parameters (e.g., value functions) defined by equilibrium constraints. We propose plug-in estimators of the bounds and two methods for inference. We also show that our bounds converge to the sharp nonparametric bounds on counterfactuals as the neighborhood size becomes large. To illustrate the broad applicability of our procedure, we present empirical applications to matching models with transferable utility and dynamic discrete choice models.


title: Majorana corner flat bands in two-dimensional second-order topological superconductors

abstract: In this paper we find that confining a second-order topological superconductor with a harmonic potential leads to a proliferation of Majorana corner modes. As a consequence, this results in the formation of Majorana corner flat bands which have a fundamentally different origin from that of the conventional mechanism. This is due to the fact that they arise solely from the one-dimensional gapped boundary states of the hybrid system that become gapless without the bulk gap closing under the increase of the trapping potential magnitude. The Majorana corner states are found to be robust against the strength of the harmonic trap and the transition from Majorana corner states to Majorana flat bands is merely a smooth crossover. As a harmonic trap can potentially be realized in heterostructures, this proposal paves a way to observe these Majorana corner flat bands in an experimental context.


title: Reconstructing the EFT of Inflation from Cosmological Data

abstract: Reconstructions of the primordial power spectrum (PPS) of curvature perturbations from cosmic microwave background anisotropies and large-scale structure data suggest that the usually assumed power-law PPS has localised features (up to $\sim 10\%$ in amplitude), although of only marginal significance in the framework of $\Lambda$CDM cosmology. On the other hand if the underlying cosmology is assumed to be Einstein-de Sitter, larger features in the PPS (up to $\sim 20\%$) are required to accurately fit the observed acoustic peaks. Within the context of single clock inflation, we show that any given reconstruction of the PPS can be mapped on to functional parameters of the underlying effective theory of the adiabatic mode within a 2nd-order formalism, provided the best fit fractional change of the PPS, $\Delta\mathcal{P}_\mathcal{R}/\mathcal{P}_\mathcal{R}$ is such that $(\Delta\mathcal{P}_\mathcal{R}/\mathcal{P}_\mathcal{R})^3$ falls within the $1\,\sigma$ confidence interval of the reconstruction for features induced by variations of either the sound speed $c_\mathrm{s}$ or the slow-roll parameter $\epsilon$. Although there is a degeneracy amongst these functional parameters (and the models that project onto them), we can identify simple representative inflationary models that yield such features in the PPS. Thus we provide a dictionary (more accurately, a thesaurus) to go from observational data, via the reconstructed PPS, to models that reproduce them to per cent level precision.


title: Long-Time Behavior of a Point Mass in a One-Dimensional Viscous Compressible Fluid and Pointwise Estimates of Solutions

abstract: We consider the motion of a point mass in a one-dimensional viscous compressible barotropic fluid. The fluid--point mass system is governed by the barotropic compressible Navier--Stokes equations and Newton's equation of motion. Our main result concerns the long-time behavior of the fluid and the point mass, and it gives pointwise convergence estimates of the volume ratio and the velocity of the fluid to their equilibrium values. As a corollary, it is shown that the velocity $V(t)$ of the point mass satisfies a decay estimate $|V(t)|=O(t^{-3/2})$ --- a faster decay compared to $t^{-1/2}$ known for the motion of a point mass in the viscous Burgers fluid~[J.~L.~V{\'{a}}zquez and E.~Zuazua, Comm. Partial Differential Equations \textbf{28} (2003), 1705--1738]. The rate $-3/2$ is essentially related to the compressibility and the nonlinearity. As a consequence, it follows that the point mass is convected only a finite distance as opposed to the viscous Burgers case. The main tool used in the proof is the pointwise estimates of Green's function. It turns out that the understanding of the time-decay properties of the transmitted and reflected waves at the point mass is essential for the proof.


title: Equivariant Multi-View Networks

abstract: Several popular approaches to 3D vision tasks process multiple views of the input independently with deep neural networks pre-trained on natural images, achieving view permutation invariance through a single round of pooling over all views. We argue that this operation discards important information and leads to subpar global descriptors. In this paper, we propose a group convolutional approach to multiple view aggregation where convolutions are performed over a discrete subgroup of the rotation group, enabling, thus, joint reasoning over all views in an equivariant (instead of invariant) fashion, up to the very last layer. We further develop this idea to operate on smaller discrete homogeneous spaces of the rotation group, where a polar view representation is used to maintain equivariance with only a fraction of the number of input views. We set the new state of the art in several large scale 3D shape retrieval tasks, and show additional applications to panoramic scene classification.


title: Quantum 4d Yang-Mills Theory and Time-Reversal Symmetric 5d Higher-Gauge Topological Field Theory

abstract: We explore 4d Yang-Mills gauge theories (YM) living as boundary conditions of 5d gapped short/long-range entangled (SRE/LRE) topological states. Specifically, we explore 4d time-reversal symmetric pure YM of an SU(2) gauge group with a second-Chern-class topological term at $\theta=\pi$ (SU(2)$_{\theta=\pi}$ YM), by turning on background fields for both the time-reversal (i.e., on unorientable manifolds) and 1-form center global symmetry. We find four classes of time-reversal and Lorentz symmetry-enriched SU(2)$_{\theta=\pi}$ YM, labeled by $(K_1, K_2)$: $K_1=0,1$ specifies Kramers singlet/doublet Wilson line and new mixed higher 't Hooft anomalies; $K_2=0,1$ specifies boson/fermionic Wilson line and a new Wess-Zumino-Witten-like counterterm. Higher anomalies indicate that to realize all higher $n$-global symmetries locally on $n$-simplices, the 4d theory becomes a boundary of a 5d higher-symmetry-protected topological state (SPTs, as an invertible topological quantum field theory (TQFT) or a cobordism invariant in math, or as a 5d higher-symmetric interacting topological superconductor in condensed matter). By dynamically gauging the 1-form symmetry, we transform a 5d bulk SRE SPTs into an LRE symmetry-enriched topologically ordered state (SETs); thus we obtain the 4d SO(3)$_{\theta=\pi}$ YM-5d LRE-higher-SETs coupled system with higher-form gauge fields. We further derive new exotic anyonic statistics of extended objects such as 2-worldsheet of strings and 3-worldvolume of branes, physically characterizing the 5d SETs. We discover triple and quadruple link invariants associated with the 5d higher-gauge TQFTs, hinting at a relation between non-supersymmetric 4d pure YM and topological links in 5d. We provide 4d-5d lattice simplicial complex regularizations and bridge to 4d quantum spin liquids. We constrain gauge dynamics by higher anomalies and a higher symmetry-extension method.

categories: hep-th cond-mat.str-el hep-ph math.AT math.GT


title: Closed maximal ideals ideals in some Fr\'echet algebras of holomorphic functions

abstract: The space $F^p$ ($1<p<\infty$) consists of all holomorphic functions $f$ on the open unit disk $\Bbb D$ such that $ \lim_{r\to 1}(1-r)^{1/q}\log^+M_{\infty}(r,f)=0,$ where $M_{\infty}(r,f)=\max_{\vert z\vert\le r}\vert f(z)\vert$ with $0<r<1$. Stoll [5, Theorem 3.2] proved that the space $F^p$ with the topology given by the family of seminorms $\left\{\Vert \cdot\Vert_{q,c}\right\}_{c>0}$ defined for $f\in F^q$ as $\Vert f\Vert_{q,c}:=\sum_{n=0}^{\infty}\vert a_n\vert\exp\left(-cn^{1/(q+1)} \right)<\infty$, becomes a countably normed Fr\'{e}chet algebra. It is known that for every $p>1$, $F^p$ is the Fr\'{e}chet envelope of the Privalov space $N^p$. In this paper, we extend our study of [32] on the structure of maximal ideals in the algebras $F^p$ ($1<p<\infty$). Namely, the obtained characterization of closed maximal ideals in $F^p$ from [32] is extended here in terms of topology of uniform convergence on compact subsets of $\Bbb D$.


title: Distributed Ledger Technology for IoT: Parasite Chain Attacks

abstract: Directed Acyclic Graph (DAG) based Distributed Ledgers can be useful in a number of applications in the IoT domain. A distributed ledger should serve as an immutable and irreversible record of transactions, however, a DAG structure is a more complicated mathematical object than its blockchain counterparts, and as a result, providing guarantees of immutability and irreversibility is more involved. In this paper, we analyse a commonly discussed attack scenario known as a parasite chain attack for the IOTA Foundation DAG based ledger. We analyse the efficacy of IOTA core MCMC algorithm using a matrix model and present an extension which improves the ledger resistance to these attacks.


title: Quantitative phase microscopy spatial signatures of cancer cells

abstract: We present cytometric classification of live healthy and cancer cells by using the spatial morphological and textural information found in the label-free quantitative phase images of the cells. We compare both healthy cells to primary tumor cell and primary tumor cells to metastatic cancer cells, where tumor biopsies and normal tissues were isolated from the same individuals. To mimic analysis of liquid biopsies by flow cytometry, the cells were imaged while unattached to the substrate. We used low-coherence off-axis interferometric phase microscopy setup, which allows a single-exposure acquisition mode, and thus is suitable for quantitative imaging of dynamic cells during flow. After acquisition, the optical path delay maps of the cells were extracted, and used to calculate 15 parameters derived from cellular 3-D morphology and texture. Upon analyzing tens of cells in each group, we found high statistical significance in the difference between the groups in most of the parameters calculated, with the same trends for all statistically significant parameters. Furthermore, a specially designed machine learning algorithm, implemented on the phase map extracted features, classified the correct cell type (healthy/cancer/metastatic) with 81%-93% sensitivity and 81%-99% specificity. The quantitative phase imaging approach for liquid biopsies presented in this paper could be the basis for advanced techniques of staging freshly isolated live cancer cells in imaging flow cytometers.


title: Microwave photonic filters via radio frequency bandwidth scaling based on soliton crystal optical micro-combs

abstract: We demonstrate high-resolution photonic RF filters using an RF bandwidth scaling approach based on integrated Kerr optical micro-combs. By employing both an active nonlinear micro-ring resonator (MRR) as a high-quality micro-comb source and a passive high-Q MRR to slice the shaped comb, a large RF instantaneous bandwidth of 4.64 GHz and a high resolution of 117 MHz are achieved, together with a broad RF operation band covering 3.28 to 19.4 GHz (L to Ku bands) using thermal tuning. We achieve programmable RF transfer functions including binary-coded notch filters and RF equalizing filters with reconfigurable slopes. Our approach is an attractive solution for high performance RF spectral shaping with high performance and flexibility.


title: Reconstruction of unknown cavity by single measurement

abstract: In this paper we propose a domain sampling type reconstruction scheme for an inverse boundary value problem to identify an unknown cavity by single measurement on the accessible boundary of a known electric or heat conductive medium. Here the single measurement is to give single current or heat flux which can have a small support over the boundary, and we measure the corresponding voltage or temperature over the whole boundary. For this inverse boundary value problem, we adapted the single NRT (no response test) introduced by Luke and Potthast for inverse scattering problem and show that it can provide such a domain sampling type reconstruction scheme.


title: An Analysis Framework for Hardware and Software Implementations with Applications from Cryptography

abstract: With the richness of present-day hardware architectures, tightening the synergy between hardware and software has attracted a great attention. The interest in unified approaches paved the way for newborn frameworks that target hardware and software co-design. This paper confirms that a unified statistical framework can successfully classify algorithms based on a combination of the heterogeneous characteristics of their hardware and software implementations. The proposed framework produces customizable indicators for any hybridization of processing systems and can be contextualized for any area of application. The framework is used to develop the Lightness Indicator System (LIS) as a case-study that targets a set of cryptographic algorithms that are known in the literature to be tiny and light. The LIS targets state-of-the-art multi-core processors and high-end Field Programmable Gate Arrays (FPGAs). The presented work includes a generic benchmark model that aids the clear presentation of the framework and extensive performance analysis and evaluation.


title: Estimation and Prediction of technical debt: a proposal

abstract: Technical debt is a metaphor used to convey the idea that doing things in a "quick and dirty" way when designing and constructing a software leads to a situation where one incurs more and more deferred future expenses. Similarly to financial debt, technical debt requires payment of interest in the form of the additional development effort that could have been avoided if the quick and dirty design choices have not been made. Technical debt applies to all the aspects of software development, spanning from initial requirements analysis to deployment, and software evolution. Technical debt is becoming very popular from scientific and industrial perspectives. In particular, there is an increase in the number of related papers over the years. There is also an increase in the number of related tools and of their adoption in the industry, especially since technical debt is very pricey and therefore needs to be managed. However, techniques to estimate technical debt are inadequate, insufficient since they mostly focus on requirements, code, and test, disregarding key artifacts such as the software architecture and the technologies used by the software at hand. Besides, despite its high relevance, technical debt prediction is one of the least explored aspects of technical debt. To address these shortcomings, it is mandatory that I undertake research to: 1) improve existing techniques to properly estimate technical debt; 2) to determine the extent to which the use of prediction techniques to foresee and therefore avoid technical debt could help companies save money and avoid a potential bankruptcy in the subsequent years. The proposed research can have an important economic impact by helping companies save several millions. It can have a major scientific impact by leading to key findings that will be disseminated through patents, well-established scientific journals and conferences.


title: On the Vulnerability of CNN Classifiers in EEG-Based BCIs

abstract: Deep learning has been successfully used in numerous applications because of its outstanding performance and the ability to avoid manual feature engineering. One such application is electroencephalogram (EEG) based brain-computer interface (BCI), where multiple convolutional neural network (CNN) models have been proposed for EEG classification. However, it has been found that deep learning models can be easily fooled with adversarial examples, which are normal examples with small deliberate perturbations. This paper proposes an unsupervised fast gradient sign method (UFGSM) to attack three popular CNN classifiers in BCIs, and demonstrates its effectiveness. We also verify the transferability of adversarial examples in BCIs, which means we can perform attacks even without knowing the architecture and parameters of the target models, or the datasets they were trained on. To our knowledge, this is the first study on the vulnerability of CNN classifiers in EEG-based BCIs, and hopefully will trigger more attention on the security of BCI systems.


title: General framework for projection structures

abstract: In the first part, we develop a general framework for projection structures and study several inference problems within this framework. We propose procedures based on data dependent measures (DDM) and make connections with empirical Bayes and penalization methods. The main inference problem is the uncertainty quantification (UQ), but on the way we solve the estimation, DDM-contraction problems, and a weak version of the structure recovery problem. The approach is local in that the quality of the inference procedures is measured by the local quantity, the oracle rate, which is the best trade-off between the approximation error by a projection structure and the complexity of that approximating projection structure. Like in statistical learning settings, we develop distribution-free theory as no particular model is imposed, we only assume certain mild condition on the stochastic part of the projection predictor. We introduce the excessive bias restriction (EBR) under which we establish the local confidence optimality of the constructed confidence ball. The proposed general framework unifies a very broad class of high-dimensional models and structures, interesting and important on their own right. In the second part, we apply the developed theory and demonstrate how the general results deliver a whole avenue of local and global minimax results (many new ones, some known results from the literature are improved) for particular models and structures as consequences, including white noise model and density estimation with smoothness structure, linear regression and dictionary learning with sparsity structures, biclustering and stochastic block models with clustering structure, covariance matrix estimation with banding and sparsity structures, and many others. Various adaptive minimax results over various scales follow also from our local results.


title: Workflow Management on the Blockchain --- Implications and Recommendations

abstract: Blockchain technology, originally popularized by cryptocurrencies, has been proposed as an infrastructure technology with applications in many areas of business management. Blockchains provide an immutable record of transactions, which makes them useful in situations where business actors may not fully trust each other. The distributed nature of blockchains makes them particularly suitable for inter-organizational e-Business applications. In this paper we examine the use of blockchains for executing inter-organizational workflows. We discuss architectural options and describe prototype implementations of blockchain-based workflow management systems (WfMS), highlighting differences to traditional WfMS. Our main contribution is the identification of potential problems raised by blockchain infrastructure and recommendations to address them.


title: The infinite dimensional Unital 3-Lie Poisson algebra

abstract: From a commutative associative algebra $A$, the infinite dimensional unital 3-Lie Poisson algebra~$\mathfrak{L}$~is constructed, which is also a canonical Nambu 3-Lie algebra, and the structure of $\mathfrak{L}$ is discussed. It is proved that: (1) there is a minimal set of generators $S$ consisting of six vectors; (2) the quotient algebra $\mathfrak{L}/\mathbb{F}L_{0, 0}^0$ is a simple 3-Lie Poisson algebra; (3) four important infinite dimensional 3-Lie algebras: 3-Virasoro-Witt algebra $\mathcal{W}_3$, $A_\omega^\delta$, $A_{\omega}$ and the 3-$W_{\infty}$ algebra can be embedded in $\mathfrak{L}$.


title: Towards Intuitive Reasoning in Axiomatic Geometry

abstract: Proving lemmas in synthetic geometry is often a time-consuming endeavour since many intermediate lemmas need to be proven before interesting results can be obtained. Improvements in automated theorem provers (ATP) in recent years now mean they can prove many of these intermediate lemmas. The interactive theorem prover Elfe accepts mathematical texts written in fair English and verifies them with the help of ATP. Geometrical texts can thereby easily be formalized in Elfe, leaving only the cornerstones of a proof to be derived by the user. This allows for teaching axiomatic geometry to students without prior experience in formalized mathematics.


title: The structure of CaSO4 nanorods -- the precursor of gypsum

abstract: Understanding the gypsum (CaSO4.2H2O) formation pathway from aqueous solution has been the subject of intensive research in the past years. This interest stems from the fact that gypsum appears to fall into a broader category of crystalline materials whose formation does not follow classical nucleation and growth theories. The pathways involve transitory precursor cluster species, yet the actual structural properties of such clusters are not very well understood. Here, we show how in situ high-energy X-ray diffraction experiments and molecular dynamics (MD) simulations can be combined to derive the structure of small CaSO4 clusters, which are precursors of crystalline gypsum. We fitted several plausible structures to the derived pair distribution functions and explored their dynamic properties using unbiased MD simulations based on both rigid ion and polarizable force fields. Determination of the structure and (meta)stability of the primary species is important from both a fundamental and applied perspective; for example, this will allow for an improved design of additives for greater control of the nucleation pathway.


title: Algorithms Clearly Beat Gamers at Quantum Moves. A Verification

abstract: The paper [S{\o}rensen et al., Nature 532] considers how human players compare to algorithms for solving the Quantum Moves game BringHomeWater and design new algorithms based on the intuition extracted from players. The claim by [S{\o}rensen et al., Nature 532] is that players outperform widely used algorithms, in particular the KASS algorithm, based on the Krotov algorithm, and that player intuition is crucial to develop improved methods. However, as initially discussed by D. Sels [D. Sels, Phys. Rev. A 97], a standard Coordinate Ascent algorithm outperforms all players by a large margin. Albeit D. Sels only compare to player solutions, the simple algorithm outperforms all algorithms based on player solutions and Krotov, and it does so using much less time and iterations. In this paper we elaborate on the methods discussed by D. Sels and verify that the presented algorithm, solves the problem better than all players and algorithms derived from player solutions in [S{\o}rensen et al., Nature 532]. We also verify the theoretical analysis presented by D. Sels, that gives a theoretically derived protocol that outperforms all players. We add a comparison with gradient ascent or GRAPE. Starting from uniform random values, GRAPE outperforms all players by a large margin. GRAPE works at least as well as the methods from [S{\o}rensen et al., Nature 532] initialized with player solutions. A standard analysis of the results from GRAPE provides a starting point for GRAPE, that outperform all algorithms from [S{\o}rensen et al., Nature 532]. We compare with a basic Krotov algorithm, and get results similar to GRAPE, clearly outperforming players and the KASS algorithm. These experiments verify and underline the result in [D. Sels, Phys. Rev. A 97] that the conclusions from [S{\o}rensen et al., Nature 532] regarding algorithms are untenable. In fact the opposite conclusions are true.


title: A benchmark for C program verification

abstract: We present twenty-five C programs, as a benchmark for C program verification using formal methods. This benchmark can be used for system demonstration, for comparison of verification effort between systems, and as a friendly competition. For this last purpose, we give a scoring formula that allows a verification system to score up to a hundred points.


abstract: Ghost imaging LiDAR via sparsity constraints using push-broom scanning is proposed. It can image the stationary target scene continuously along the scanning direction by taking advantage of the relative movement between the platform and the target scene. Compared to conventional ghost imaging LiDAR that requires multiple speckle patterns staring the target, ghost imaging LiDAR via sparsity constraints using push-broom scanning not only simplifies the imaging system, but also reduces the sampling number. Numerical simulations and experiments have demonstrated its efficiency.


title: Oriented suspension mechanics with applications to flow linear dichrosim spectroscopy and pathogen detection

abstract: Flow linear dichroism is a biophysical spectroscopic technique that exploits the shear-induced alignment of elongated particles in suspension. Motivated by the broad aim of optimising the sensitivity of this technique, and more specifically by a handheld synthetic biotechnology prototype for waterborne-pathogen detection, a model of steady and oscillating pressure-driven channel flow and orientation dynamics of a suspension of slender microscopic fibres is developed. The model couples the Fokker-Planck equation for Brownian suspensions with the narrow channel flow equations, the latter modified to incorporate mechanical anisotropy induced by the particles. The linear dichroism signal is estimated through integrating the perpendicular components of the distribution function via an appropriate formula which takes the bi-axial nature of the orientation into account. For the specific application of pathogen detection via binding of M13 bacteriophage, it is found that increases in the channel depth are more significant in improving the linear dichroism signal than increases in the channel width. Increasing channel depth to $2~$mm and pressure gradient to $5 \times 10^4~$Pa/m essentially maximises the alignment. Oscillating flow can produce nearly equal alignment to steady flow at appropriate frequencies, which has significant potential practical value in the analysis of small sample volumes.


title: Synergic Cherenkov-Compton Radiation

abstract: In vacuum with an induced by strong electromagnetic wave refraction index larger than unity, an ultra-relativistic electron emits Cherenkov radiation. During the interaction with this wave the electron also radiates photons via the Compton scattering. Synergic Cherenkov-Compton process can be observed by colliding laser accelerated electrons with a high intensity electromagnetic pulse. Extremely high energy photons cannot be emitted via the Cherenkov radiation because the vacuum refraction index tends to unity at these energies. Experiments on studying these phenomena will reveal the properties of vacuum predicted by nonlinear quantum electrodynamics.


title: Coupling Poynting-Robertson Effect in Mass Accretion Flow Physics

abstract: In my doctoral thesis, I have focussed my attention on radiation processes in high-energy astrophysics connected with the accretion flow physics around compact objects. Generally, a radiation field beside to exert an outward radiation pressure, there is also the presence of a radiation drag force, which both can drastically change or even halt the motion of the surrounding matter. The radiation drag force, known as Poynting-Robertson effect, acts as a dissipative force against the matter's orbital motion, removing very efficiently angular momentum and energy from it. The thesis is organised in three parts: (1) for ray-tracing purposes, I have developed a mathematical method for deriving a set of high-accurate approximate polynomial formulae to easily integrate photon geodesics in a Schwarzschild spacetime; (2) I gave two fundamental contributions in the field of the general relativistic treatment of the Poynting-Robertson effect (Lagrangian formulations and extension of the model in three dimensions); (3) I reduced the data of three accreting millisecond X-ray pulsars: IGR J00291+5934, IGR J18245-2452, and SAX J1748.9-2021. This thesis offers innovative ideas in the field of radiation processes involving the Poynting-Robertson effect in high-energy astrophysics, opening thus up future interesting perspectives both in theoretical and observational physics.


title: Comparison of Possibilistic Fuzzy Local Information C-Means and Possibilistic K-Nearest Neighbors for Synthetic Aperture Sonar Image Segmentation

abstract: Synthetic aperture sonar (SAS) imagery can generate high resolution images of the seafloor. Thus, segmentation algorithms can be used to partition the images into different seafloor environments. In this paper, we compare two possibilistic segmentation approaches. Possibilistic approaches allow for the ability to detect novel or outlier environments as well as well known classes. The Possibilistic Fuzzy Local Information C-Means (PFLICM) algorithm has been previously applied to segment SAS imagery. Additionally, the Possibilistic K-Nearest Neighbors (PKNN) algorithm has been used in other domains such as landmine detection and hyperspectral imagery. In this paper, we compare the segmentation performance of a semi-supervised approach using PFLICM and a supervised method using Possibilistic K-NN. We include final segmentation results on multiple SAS images and a quantitative assessment of each algorithm.


title: Spin-charge coupled transport in van der Waals systems with random tunneling

abstract: We study the electron and spin transport in a van der Waals system formed by one layer with strong spin-orbit coupling and a second layer without spin-orbit coupling, in the regime when the interlayer tunneling is random. We find that in the layer without intrinsic spin-orbit coupling spin-charge coupled transport can be induced by two distinct mechanisms. First, the gapless diffusion modes of the two isolated layers hybridize in the presence of tunneling, which constitutes a source of spin-charge coupled transport in the second layer. Second, the random tunneling introduces spin-orbit coupling in the effective disorder-averaged single-particle Hamiltonian of the second layer. This results in non-trivial spin transport and, for sufficiently strong tunneling, in spin-charge coupling. As an example, we consider a van der Waals system formed by a two-dimensional electron gas (2DEG)--such as graphene--and the surface of a topological insulator (TI) and show that the proximity of the TI induces a coupling of the spin and charge transport in the 2DEG. In addition, we show that such coupling can be tuned by varying the doping of the TI's surface. We then obtain, for a simple geometry, the current-induced non-equilibrium spin accumulation (Edelstein effect) caused in the 2DEG by the coupling of charge and spin transport.


title: Rock 'n' Roll Solutions to the Hubble Tension

abstract: Local measurements of the Hubble parameter are increasingly in tension with the value inferred from a $\Lambda$CDM fit to the cosmic microwave background (CMB) data. In this paper, we construct scenarios in which evolving scalar fields significantly ease this tension by adding energy to the Universe around recombination in a narrow redshift window. We identify solutions with scalar field potential $V \propto \phi^{2n}$ that have simple asymptotic behavior, both oscillatory (rocking) and rolling. These solutions consistently describe both the field evolution and its fluctuations without approximation. Our findings differ qualitatively from some of the existing literature, which rely upon a coarse-grained fluid description. Combining CMB data with low-redshift measurements, the best fit model has $n=2$ with a significantly higher value of the Hubble constant as compared to a $\Lambda$CDM fit to the same data. Future measurements of the late-time amplitude of matter fluctuations and of the reionization history could help distinguish these models from competing solutions.


title: On Anomaly-Free Dark Matter Models

abstract: We investigate the predictions of anomaly-free dark matter models for direct and indirect detection experiments. We focus on gauge theories where the existence of a fermionic dark matter candidate is predicted by anomaly cancellation, its mass is defined by the new symmetry breaking scale, and its stability is guaranteed by a remnant symmetry after the breaking of the gauge symmetry. We find an upper bound on the symmetry breaking scale by applying the relic density and perturbative constraints. The anomaly-free property of the theories allows us to perform a full study of the gamma lines from dark matter annihilation. We investigate the correlation between predictions for final radiation processes and gamma lines. Furthermore, we demonstrate that the latter can be distinguished from the continuum gamma ray spectrum.


title: On the convergence of the gradient expansion in hydrodynamics

abstract: Hydrodynamic excitations corresponding to sound and shear modes in fluids are characterised by gapless dispersion relations. In the hydrodynamic gradient expansion, their frequencies are represented by power series in spatial momenta. We investigate the analytic structure and convergence properties of the hydrodynamic series by studying the associated spectral curve in the space of complexified frequency and complexified spatial momentum. For the strongly coupled ${\cal N}=4$ supersymmetric Yang-Mills plasma, we use the holographic duality methods to demonstrate that the derivative expansions have finite non-zero radii of convergence. Obstruction to the convergence of hydrodynamic series arises from level-crossings in the quasinormal spectrum at complex momenta.


title: Theoretical Reflectance Spectra of Earth-Like Planets through Their Evolutions: Impact of Clouds on the Detectability of Oxygen, Water, and Methane with Future Direct Imaging Missions

abstract: In the near-future, atmospheric characterization of Earth-like planets in the habitable zone will become possible via reflectance spectroscopy with future telescopes such as the proposed LUVOIR and HabEx missions. While previous studies have considered the effect of clouds on the reflectance spectra of Earth-like planets, the molecular detectability considering a wide range of cloud properties has not been previously explored in detail. In this study, we explore the effect of cloud altitude and coverage on the reflectance spectra of Earth-like planets at different geological epochs and examine the detectability of $\mathrm{O_2}$, $\mathrm{H_2O}$, and $\mathrm{CH_4}$ with test parameters for the future mission concept, LUVOIR, using a coronagraph noise simulator previously designed for WFIRST-AFTA. Considering an Earth-like planet located at 5 pc away, we have found that for the proposed LUVOIR telescope, the detection of the $\mathrm{O_2}$ A-band feature (0.76 $\mathrm{\mu}$m) will take approximately 100, 30, and 10 hours for the majority of the cloud parameter space modeled for the atmospheres with 10\%, 50\%, and 100\% of modern Earth O$_2$ abundances, respectively. Especially, for {the case of $\gtrsim 50$\%} of modern Earth O$_2$ abundance, the feature will be detectable with integration time $\lesssim 10$ hours as long as there are lower altitude ($\lesssim 8$ km) clouds with a global coverage of $\gtrsim 20\%$. For the 1\% of modern Earth $\mathrm{O_2}$ abundance case, however, it will take more than 100 hours for all the cloud parameters we modeled.


title: Seeking for sterile neutrinos with displaced leptons at the LHC

abstract: We study the signal of long-lived sterile neutrino at the LHC produced through the decay of the $W$ boson. It decays into charged lepton and jets. The characteristic signature is a hard prompt lepton and a lepton from the displaced decay of the sterile neutrino, which leads to a bundle of displaced tracks with large transverse impact parameter. Different from other studies, we neither reconstruct the displaced vertex nor place requirement on its invariant mass to maintain sensitivity for low sterile neutrino masses. Instead, we focus on the displaced track from the lepton. A difficulty for low mass sterile neutrino study is that the displaced lepton is usually \textit{non-isolated}. Therefore, leptons from heavy flavor quark is the major source of background. We closely follow a search for displaced electron plus muon search at CMS and study their control regions, which is related to our signal regions, in great detail to develop a robust estimation of the background for our signals. After further optimization on the signal limiting the number of jets, low $H_T$ and large lepton displacement $d_0$ to suppress SM background, we reach an exclusion sensitivity of about $10^{-8}$ ($10^{-5}$) for the mixing angle square at 10 (2) GeV sterile neutrino mass respectively. The strategy we propose can cover the light sterile masses complimentary to beam dump and forward detector experiments.


title: Gaia DR2 orbital properties for field stars with globular cluster-like CN band strengths

abstract: Large spectroscopic surveys of the Milky Way have revealed that a small population of stars in the halo have light element abundances comparable to those found in globular clusters. The favoured explanation for the peculiar abundances of these stars is that they originated inside a globular cluster and were subsequently lost. Using orbit calculations we assess the likelihood that an existing sample of 57 field stars with globular cluster-like CN band strength originated in any of the currently known Milky Way globular clusters. Using Sloan Digital Sky Survey and Gaia data, we determine orbits and integrals of motion of our sample of field stars, and use these values and metallicity to identify likely matches to globular clusters. The pivot hypothesis is that had these stars been stripped from such objects, they would have remained on very similar orbits. We find that ~ 70% of the sample of field stars have orbital properties consistent with the halo of the Milky Way; however, only 20 stars have likely orbital associations with an existing globular cluster. The remaining ~ 30% of the sample have orbits that place them in the outer Galactic disc. No cluster of similar metallicity is known on analogous disc orbits. The orbital properties of the halo stars seem to be compatible with the globular cluster escapee scenario. The stars in the outer disc are particularly surprising and deserve further investigation to establish their nature.


title: Cosmological exploitation of the size function of cosmic voids identified in the distribution of biased tracers

abstract: Cosmic voids are large underdense regions that, together with galaxy clusters, filaments and walls, build up the large-scale structure of the Universe. The void size function provides a powerful probe to test the cosmological framework. However, to fully exploit this statistics, the void sample has to be properly cleaned from spurious objects. Furthermore, the bias of the mass tracers used to detect these regions has to be taken into account in the size function model. In our work we test a cleaning algorithm and a new void size function model on a set of simulated dark matter halo catalogues, with different mass and redshift selections, to investigate the statistics of voids identified in a biased mass density field. We then investigate how the density field tracers' bias affects the detected size of voids. The main result of this analysis is a new model of the size function, parameterised in terms of the linear effective bias of the tracers used, which is straightforwardly inferred from the large-scale two-point correlation function. This represents a crucial step to exploit the method on real data catalogues. The proposed size function model has been accurately calibrated on mock catalogues, and used to validate the possibility to provide forecasts on the cosmological constraints, namely on the matter density contrast, $\Omega_{\rm M}$, and on the normalisation of the linear matter power spectrum, $\sigma_8$, at different redshifts.


title: Fractional Electromagnetism in Quantum Matter and High-Energy Physics

abstract: We present here a theory of fractional electro-magnetism which is capable of describing phenomenon as disparate as the non-locality of the Pippard kernel in superconductivity and anomalous dimensions for conserved currents in holographic dilatonic models. The starting point for our analysis is the observation that the standard current conservation equations remain unchanged if any differential operator that commutes with the total exterior derivative multiplies the current. Such an operator, effectively changing the dimension of the current, increases the allowable gauge transformations in electromagnetism and is at the heart of N\"other's second theorem. Here we develop a consistent theory of electromagnetism that exploits this hidden redundancy in which the standard gauge symmetry in electromagnetism is modified by the rotationally invariant operator, the fractional Laplacian. We show that the resultant theories all allow for anomalous (non-traditional) scaling dimensions of the gauge field and the associated current. Using well known extension theorems and the membrane paradigm, we show that either the boundary (UV) or horizon (IR) theory of holographic dilatonic models are both described by such fractional electromagnetism. We also show that the non-local Pippard kernel introduced to solve the problem of the Meissner effect in elemental superconductors can also be formulated as a special case of fractional electromagnetism. We show that the standard charge quantization rules fail when the gauge field acquires an anomalous dimension. The breakdown of charge quantization is discussed extensively in terms of the experimentally measurable modified Aharonov-Bohm effect in the strange metal phase of the cuprate superconductors.


title: Ghost imaging with engineered quantum states by Hong-Ou-Mandel interference

abstract: Traditional ghost imaging experiments exploit position correlations between correlated states of light. These correlations occur directly in spontaneous parametric down-conversion (SPDC), and in such a scenario, the two-photon state used for ghost imaging is symmetric. Here we perform ghost imaging using an anti-symmetric state, engineering the two-photon state symmetry by means of Hong-Ou-Mandel interference. We use both symmetric and anti-symmetric states and show that the ghost imaging setup configuration results in object-image rotations depending on the state selected. Further, the object and imaging arms employ spatial light modulators for the all-digital control of the projections, being able to dynamically change the measuring technique and the spatial properties of the states under study. Finally, we provide a detailed theory that explains the reported observations.


title: Kitaev honeycomb models in magnetic fields: Dynamical response and hidden symmetries

abstract: Motivated by recent reports of a field-induced intermediate phase (IP) in the antiferromagnetic honeycomb Kitaev model that may be a spin liquid whose nature is distinct from the Kitaev Z2 phase, we present a detailed numerical study on the nature and dynamical response (such as dynamical spin-structure factors and resonant inelastic x-ray scattering intensities) of this field-induced IP and neighboring phases in a family of Kitaev-based models related by hidden symmetries and duality transformations. We further show that the same field-induced IP can appear in models relevant for $\alpha$-RuCl${_3}$, which exhibit a ferromagnetic Kitaev coupling and additional interactions. In $\alpha$-RuCl${_3}$, the IP represents a new phase, that is likely independent from the putative field-induced (spin-liquid) phase recently reported from thermal Hall conductivity measurements.


title: Dicke time crystals in driven-dissipative quantum many-body systems

abstract: The Dicke model -- a paradigmatic example of superradiance in quantum optics -- describes an ensemble of atoms which are collectively coupled to a leaky cavity mode. As a result of the cooperative nature of these interactions, the system's dynamics are captured by the behavior of a single mean-field, collective spin. In this mean-field limit, it has recently been shown that the interplay between photon losses and periodic driving of light-matter coupling can lead to time-crystalline-like behavior of the collective spin. In this work, we investigate whether such a Dicke time crystal is stable to perturbations that explicitly break the mean-field solvability of the conventional Dicke model. In particular, we consider the addition of short-range interactions between atoms, which breaks the collective coupling and leads to complex many-body dynamics. In this context, the interplay between periodic driving, dissipation and interactions yields a rich set of dynamical responses including long-lived and metastable Dicke time crystals, where losses can cool down the many-body heating resulting from the continuous pump of energy from the periodic drive. Specifically, when the additional short-range interactions are ferromagnetic, we observe time crystalline behavior at non-perturbative values of the coupling strength, suggesting the possible existence of stable dynamical order in a driven-dissipative quantum many-body system. These findings illustrate the rich nature of novel dynamical responses with many-body character in quantum optics platforms.


title: Superpotential of Three Dimensional ${\cal N}=1$ Heterotic Supergravity

abstract: We dimensionally reduce the ten dimensional heterotic action on spacetimes of the form ${\cal M}_{(2,1)}\times Y$, where ${\cal M}_{(2,1)}$ is three dimensional maximally symmetric Anti de Sitter or Minkowski space, and $Y$ is a compact seven dimensional manifold with $G_2$ structure. In doing so, we derive the real superpotential functional of the corresponding three dimensional ${\cal N}=1$ theory. We confirm that extrema of this functional precisely correspond to supersymmetric heterotic compactifications on manifolds of $G_2$ structure. We make some comments on the role of the superpotential functional with respect to the coupled moduli problem of instanton bundles over $G_2$ manifolds.


title: Satellites of Satellites: The Case for Carina and Fornax

abstract: We use the Auriga cosmological simulations of Milky Way (MW)-mass galaxies and their surroundings to study the satellite populations of dwarf galaxies in $\Lambda$CDM. As expected from prior work, the number of satellites above a fixed stellar mass is a strong function of the mass of the primary dwarf. For galaxies as luminous as the Large Magellanic Cloud (LMC), and for halos as massive as expected for the LMC (determined by its rotation speed), the simulations predict about 3 satellites with stellar masses exceeding $M_*>10^5\, M_\odot$. If the LMC is on its first pericentric passage, then these satellites should be near the LMC and should have orbital angular momenta roughly coincident with that of the LMC. We use 3D positions and velocities from the 2nd data release of the Gaia mission to revisit which of the "classical" MW dwarf spheroidals could plausibly be LMC satellites. The new proper motions of the Fornax and Carina dwarf spheroidals place them on orbits closely aligned with the orbital plane of the Magellanic Clouds, hinting at a potential Magellanic association. Together with the Small Magellanic Cloud (SMC), this result raises to $3$ the number of LMC satellites with $M_*>10^5\, M_\odot$, as expected from simulations. This also fills the 12-mag luminosity gap between the SMC and the ultra-faints Hyi1, Car2, Hor1, and Car3, the few ultra-faint satellites confirmed to have orbits consistent with a Magellanic origin.


title: Probing Primordial Chirality with Galaxy Spins

abstract: Chiral symmetry is maximally violated in weak interactions, and such microscopic asymmetries in the early Universe might leave observable imprints on astrophysical scales without violating the cosmological principle. In this Letter, we propose a helicity measurement to detect primordial chiral violation. We point out that observations of halo-galaxy angular momentum directions (spins), which are frozen in during the galaxy formation process, provide a fossil chiral observable. From the clustering mode of large scale structure of the Universe, we construct a spin mode in Lagrangian space and show in simulations that it is a good probe of halo-galaxy spins. In standard model, a strong symmetric correlation between the left and right helical components of this spin mode and galaxy spins is expected. Measurements of these correlations will be sensitive to chiral breaking, providing a direct test of chiral symmetry breaking in the early Universe.


title: Beyond BAO: improving cosmological constraints from BOSS with measurement of the void-galaxy cross-correlation

abstract: We present a measurement of the anisotropic void-galaxy cross-correlation function in the CMASS galaxy sample of the BOSS DR12 data release. We perform a joint fit to the data for redshift space distortions (RSD) due to galaxy peculiar velocities and anisotropies due to the Alcock-Paczynski (AP) effect, for the first time using a velocity field reconstruction technique to remove the complicating effects of RSD in the void centre positions themselves. Fits to the void-galaxy function give a 1% measurement of the AP parameter combination $D_A(z)H(z)/c = 0.4367\pm 0.0045$ at redshift $z=0.57$, where $D_A$ is the angular diameter distance and $H$ the Hubble parameter, exceeding the precision obtainable from baryon acoustic oscillations (BAO) by a factor of ~3.5 and free of systematic errors. From voids alone we also obtain a 10% measure of the growth rate, $f\sigma_8(z=0.57)=0.501\pm0.051$. The parameter degeneracies are orthogonal to those obtained from galaxy clustering. Combining void information with that from BAO and galaxy RSD in the same CMASS sample, we measure $D_A(0.57)/r_s=9.383\pm 0.077$ (at 0.8% precision), $H(0.57)r_s=(14.05\pm 0.14)\;10^3$ kms$^{-1}$Mpc$^{-1}$ (1%) and $f\sigma_8=0.453\pm0.022$ (4.9%), consistent with cosmic microwave background (CMB) measurements from Planck. These represent a factor \sim2 improvement in precision over previous results through the inclusion of void information. Fitting a flat cosmological constant $\Lambda$CDM model to these results in combination with Planck CMB data, we find up to an 11% reduction in uncertainties on $H_0$ and $\Omega_m$ compared to use of the corresponding BOSS consensus values. Constraints on extended models with non-flat geometry and a dark energy of state that differs from $w=-1$ show an even greater improvement.


title: Modular Synthesis of Divide-and-Conquer Parallelism for Nested Loops (Extended Version)

abstract: We propose a methodology for automatic generation of divide-and-conquer parallel implementations of sequential nested loops. We focus on a class of loops that traverse read-only multidimensional collections (lists or arrays) and compute a function over these collections. Our approach is modular, in that, the inner loop nest is abstracted away to produce a simpler loop nest for parallelization. Then, the summarized version of the loop nest is parallelized. The main challenge addressed by this paper is that to perform the code transformations necessary in each step, the loop nest may have to be augmented (automatically) with extra computation to make possible the abstraction and/or the parallelization tasks. We present theoretical results to justify the correctness of our modular approach, and algorithmic solutions for automation. Experimental results demonstrate that our approach can parallelize highly non-trivial loop nests efficiently.


title: Learning to Stop in Structured Prediction for Neural Machine Translation

abstract: Beam search optimization resolves many issues in neural machine translation. However, this method lacks principled stopping criteria and does not learn how to stop during training, and the model naturally prefers the longer hypotheses during the testing time in practice since they use the raw score instead of the probability-based score. We propose a novel ranking method which enables an optimal beam search stopping criteria. We further introduce a structured prediction loss function which penalizes suboptimal finished candidates produced by beam search during training. Experiments of neural machine translation on both synthetic data and real languages (German-to-English and Chinese-to-English) demonstrate our proposed methods lead to better length and BLEU score.


title: Multitask Soft Option Learning

abstract: We present Multitask Soft Option Learning(MSOL), a hierarchical multitask framework based on Planning as Inference. MSOL extends the concept of options, using separate variational posteriors for each task, regularized by a shared prior. This ''soft'' version of options avoids several instabilities during training in a multitask setting, and provides a natural way to learn both intra-option policies and their terminations. Furthermore, it allows fine-tuning of options for new tasks without forgetting their learned policies, leading to faster training without reducing the expressiveness of the hierarchical policy. We demonstrate empirically that MSOL significantly outperforms both hierarchical and flat transfer-learning baselines.


title: Geometry of small causal diamonds

abstract: The geometry of small causal diamonds is systematically studied, based on three distinct constructions that are common in the literature, namely the geodesic ball, the Alexandrov interval and the lightcone cut. The causal diamond geometry is calculated perturbatively using Riemann normal coordinate expansion up to the leading order in both vacuum and non-vacuum. We provide a collection of results including the area of the codimension-two edge, the maximal hypersurface volume and their isoperimetric ratio for each construction, which will be useful for any applications involving the quantitative properties of causal diamonds. In particular, by solving the dynamical equations of the expansion and the shear on the lightcone, we find that intriguingly only the light-cone cut construction yields an area deficit proportional to the Bel-Robinson superenergy density W in four dimensional spacetime, but such a direct connection fails to hold in any other dimension. We also compute the volume of the Alexandrov interval causal diamond in vacuum, which we believe is important but missing from the literatures. Our work complements and extends the earlier works on the causal diamond geometry by Gibbons and Solodukhin, Jacobson, Senovilla and Speranza and others. Some potential applications of our results in mathematical general relativity and quantum gravity are discussed.


title: Entanglement Content of Quantum Particle Excitations II. Disconnected Regions and Logarithmic Negativity

abstract: In this paper we study the increment of the entanglement entropy and of the (replica) logarithmic negativity in a zero-density excited state of a free massive bosonic theory, compared to the ground state. This extends the work of two previous publications by the same authors. We consider the case of two disconnected regions and find that the change in the entanglement entropy depends only on the combined size of the regions and is independent of their connectivity. We subsequently generalize this result to any number of disconnected regions. For the replica negativity we find that its increment is a polynomial with integer coefficients depending only on the sizes of the two regions. The logarithmic negativity turns out to have a more complicated functional structure than its replica version, typically involving roots of polynomials on the sizes of the regions. We obtain our results by two methods already employed in previous work: from a qubit picture and by computing four-point functions of branch point twist fields in finite volume. We test our results against numerical simulations on a harmonic chain and find excellent agreement.

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


title: Post-Selection and Counterfactual Communication

abstract: A recent quantum protocol for counterfactual communication [Y. Aharonov and L. Vaidman, Phys. Rev. A 99, 010103(R), 2019] relies on post-selection to eliminate the weak trace in the transmission channel. We show that the post-selection in this protocol also eliminates the flow of Fisher information from transmitter to receiver. However, we also show that a classical communication protocol with post-selection can be counterfactual. Hence, we argue that post-selection should not be allowed in genuine counterfactual communication. In the quantum counterfactual communication protocol, the probability of discarding an event by post-selection tends to zero with an increasing number of ideal optical components. But the counterfactual violation strength tends to infinity at a faster rate. Consequently, the quantum protocol is not counterfactual proper.


title: An effective Lie--Kolchin theorem for quasi-unipotent matrices

abstract: We establish an effective version of the classical Lie--Kolchin Theorem. Namely, let $A,B\in\mathrm{GL}_m(\mathbb{C})$ be quasi--unipotent matrices such that the Jordan Canonical Form of $B$ consists of a single block, and suppose that for all $k\geq0$ the matrix $AB^k$ is also quasi--unipotent. Then $A$ and $B$ have a common eigenvector. In particular, $\langle A,B\rangle<\mathrm{GL}_m(\mathbb{C})$ is a solvable subgroup. We give applications of this result to the representation theory of mapping class groups of orientable surfaces.