title: Implications of dark matter cascade decay from DAMPE, HESS, Fermi-LAT and AMS02 data

abstract: Recent high-energy cosmic $e^\pm$ measurement from the DArk Matter Particle Explorer (DAMPE) satellite confirms the deviation of total cosmic ray electron spectrum above 700-900 GeV from a simple power law. In this paper we demonstrate that the cascade decay of dark matter (DM) can account for DAMPE's TeV $e^+e^-$ spectrum. We select the least constraint DM decay channel into four muons as the benchmark scenario, and perform an analysis with propagation variance in both DM signal and the Milky Way's electron background. The best-fit of the model is obtained for joint DAMPE, Fermi-Large Area Telescope (Fermi-LAT), High Energy Stereoscopic System (HESS), high energy electron data sets, and with an $\mathcal{O}(10^{26})$ second decay lifetime, which is consistent with existing gamma ray and cosmic microwave background limits. We compare the spectral difference between the cascade decay of typical final-state channels. The least constrained $4\mu$ channels give good fits to the electron spectrum's TeV scale down-turn, yet their low energy spectrum has tension with sub-TeV positron data from AMS02. We also consider a three-step cascade decay into eight muons, and also a gamma-ray constrained $4\mu,4b$ mixed channel, to demonstrate that a further softened cascade decay signal would be required for the agreement with all the data sets.


title: Deep Neural Network Detects Quantum Phase Transition

abstract: We detect the quantum phase transition of a quantum many-body system by mapping the observed results of the quantum state onto a neural network. In the present study, we utilized the simplest case of a quantum many-body system, namely a one-dimensional chain of Ising spins with the transverse Ising model. We prepared several spin configurations, which were obtained using repeated observations of the model for a particular strength of the transverse field, as input data for the neural network. Although the proposed method can be employed using experimental observations of quantum many-body systems, we tested our technique with spin configurations generated by a quantum Monte Carlo simulation without initial relaxation. The neural network successfully classified the strength of transverse field only from the spin configurations, leading to consistent estimations of the critical point of our model $\Gamma_c =J$.


title: Bayesian analysis of the break in DAMPE lepton spectra

abstract: Recently, DAMPE has released its first results on the high-energy cosmic-ray electrons and positrons (CREs) from about $25$ GeV to $4.6$ TeV, which directly detect a break at $\sim 1$ TeV. This result gives us an excellent opportunity to study the source of the CREs excess. In this work, we used the data fo proton and helium flux (from AMS-02 and CREAM), $\bar{\mathrm{p}}/\mathrm{p}$ ratio (from AMS-02), positron flux (from AMS-02) and CREs flux (from DAMPE without the peak signal point at $\sim 1.4$ TeV) to do global fitting simultaneously, which can account the influence from the propagation model, the nuclei and electron primary source injection and the secondary lepton production precisely. For extra source to interpret the excess in lepton spectrum, we consider two separate scenarios (pulsar and dark matter annihilation via leptonic channels) to construct the bump ($\gtrsim 100$ GeV) and the break at $\sim 1$ TeV. The result shows: (i) in pulsar scenario, the spectral index of the injection should be $\nu_{\mathrm{psr}} \sim 0.65$ and the cut-off should be $R_{c} \sim 650$ GV; (ii) in dark matter scenario, the dark matter particle's mass is $m_{\chi} \sim 1208$ GeV and the cross section is $\langle \sigma v \rangle \sim 1.48 \times 10^{-23} \mathrm{cm}^{3} \mathrm{s}^{-1}$. Moreover, in the dark matter scenario, the $\tau \bar{\tau}$ annihilation channel is highly suppressed, and a DM model is built to satisfy the fitting results.


title: Joint Access-Backhaul Perspective on Mobility Management in 5G Networks

abstract: The ongoing efforts in the research development and standardization of 5G, by both industry and academia, have resulted in the identification of enablers (Software Defined Networks, Network Function Virtualization, Distributed Mobility Management, etc.) and critical areas (Mobility management, Interference management, Joint access-backhaul mechanisms, etc.) that will help achieve the 5G objectives. During these efforts, it has also been identified that the 5G networks due to their high degree of heterogeneity, high QoS demand and the inevitable density (both in terms of access points and users), will need to have efficient joint backhaul and access mechanisms as well as enhanced mobility management mechanisms in order to be effective, efficient and ubiquitous. Therefore, in this paper we first provide a discussion on the evolution of the backhaul scenario, and the necessity for joint access and backhaul optimization. Subsequently, and since mobility management mechanisms can entail the availability, reliability and heterogeneity of the future backhaul/fronthaul networks as parameters in determining the most optimal solution for a given context, a study with regards to the effect of future backhaul/fronthaul scenarios on the design and implementation of mobility management solutions in 5G networks has been performed.


title: Precision Learning: Towards Use of Known Operators in Neural Networks

abstract: In this paper, we consider the use of prior knowledge within neural networks. In particular, we investigate the effect of a known transform within the mapping from input data space to the output domain. We demonstrate that use of known transforms is able to change maximal error bounds. In order to explore the effect further, we consider the problem of X-ray material decomposition as an example to incorporate additional prior knowledge. We demonstrate that inclusion of a non-linear function known from the physical properties of the system is able to reduce prediction errors therewith improving prediction quality from SSIM values of 0.54 to 0.88. This approach is applicable to a wide set of applications in physics and signal processing that provide prior knowledge on such transforms. Also maximal error estimation and network understanding could be facilitated within the context of precision learning.


title: Maximum-width Axis-Parallel Empty Rectangular Annulus

abstract: Given a set $P$ of $n$ points on $\mathbb R^{2}$, we address the problem of computing an axis-parallel empty rectangular annulus $A$ of maximum-width such that no point of $P$ lies inside $A$ but all points of $P$ must lie inside, outside and on the boundaries of two parallel rectangles forming the annulus $A$. We propose an $O(n^3)$ time and $O(n)$ space algorithm to solve the problem. In a particular case when the inner rectangle of an axis-parallel empty rectangular annulus reduces to an input point we can solve the problem in $O(n \log n)$ time and $O(n)$ space.


title: Polar Coding for the Large Hadron Collider: Challenges in Code Concatenation

abstract: In this work, we present a concatenated repetition-polar coding scheme that is aimed at applications requiring highly unbalanced unequal bit-error protection, such as the Beam Interlock System of the Large Hadron Collider at CERN. Even though this concatenation scheme is simple, it reveals significant challenges that may be encountered when designing a concatenated scheme that uses a polar code as an inner code, such as error correlation and unusual decision log-likelihood ratio distributions. We explain and analyze these challenges and we propose two ways to overcome them.


title: Don't Just Assume; Look and Answer: Overcoming Priors for Visual Question Answering

abstract: A number of studies have found that today's Visual Question Answering (VQA) models are heavily driven by superficial correlations in the training data and lack sufficient image grounding. To encourage development of models geared towards the latter, we propose a new setting for VQA where for every question type, train and test sets have different prior distributions of answers. Specifically, we present new splits of the VQA v1 and VQA v2 datasets, which we call Visual Question Answering under Changing Priors (VQA-CP v1 and VQA-CP v2 respectively). First, we evaluate several existing VQA models under this new setting and show that their performance degrades significantly compared to the original VQA setting. Second, we propose a novel Grounded Visual Question Answering model (GVQA) that contains inductive biases and restrictions in the architecture specifically designed to prevent the model from 'cheating' by primarily relying on priors in the training data. Specifically, GVQA explicitly disentangles the recognition of visual concepts present in the image from the identification of plausible answer space for a given question, enabling the model to more robustly generalize across different distributions of answers. GVQA is built off an existing VQA model -- Stacked Attention Networks (SAN). Our experiments demonstrate that GVQA significantly outperforms SAN on both VQA-CP v1 and VQA-CP v2 datasets. Interestingly, it also outperforms more powerful VQA models such as Multimodal Compact Bilinear Pooling (MCB) in several cases. GVQA offers strengths complementary to SAN when trained and evaluated on the original VQA v1 and VQA v2 datasets. Finally, GVQA is more transparent and interpretable than existing VQA models.


title: Time Limits in Reinforcement Learning

abstract: In reinforcement learning, it is common to let an agent interact for a fixed amount of time with its environment before resetting it and repeating the process in a series of episodes. The task that the agent has to learn can either be to maximize its performance over (i) that fixed period, or (ii) an indefinite period where time limits are only used during training to diversify experience. In this paper, we provide a formal account for how time limits could effectively be handled in each of the two cases and explain why not doing so can cause state aliasing and invalidation of experience replay, leading to suboptimal policies and training instability. In case (i), we argue that the terminations due to time limits are in fact part of the environment, and thus a notion of the remaining time should be included as part of the agent's input to avoid violation of the Markov property. In case (ii), the time limits are not part of the environment and are only used to facilitate learning. We argue that this insight should be incorporated by bootstrapping from the value of the state at the end of each partial episode. For both cases, we illustrate empirically the significance of our considerations in improving the performance and stability of existing reinforcement learning algorithms, showing state-of-the-art results on several control tasks.


title: Robust Computation in 2D Absolute EIT (a-EIT) Using D-bar Methods with the `exp' Approximation

abstract: Objective: Absolute images have important applications in medical Electrical Impedance Tomography (EIT) imaging, but the traditional minimization and statistical based computations are very sensitive to modeling errors and noise. In this paper, it is demonstrated that D-bar reconstruction methods for absolute EIT are robust to such errors. Approach: The effects of errors in domain shape and electrode placement on absolute images computed with 2D D-bar reconstruction algorithms are studied on experimental data. Main Results: It is demonstrated with tank data from several EIT systems that these methods are quite robust to such modeling errors, and furthermore the artefacts arising from such modeling errors are similar to those occurring in classic time-difference EIT imaging. Significance: This study is promising for clinical applications where absolute EIT images are desirable, but previously thought impossible.


title: An adaptive fast Gauss transform in two dimensions

abstract: A variety of problems in computational physics and engineering require the convolution of the heat kernel (a Gaussian) with either discrete sources, densities supported on boundaries, or continuous volume distributions. We present a unified fast Gauss transform for this purpose in two dimensions, making use of an adaptive quad-tree discretization on a unit square which is assumed to contain all sources. Our implementation permits either free-space or periodic boundary conditions to be imposed, and is efficient for any choice of variance in the Gaussian.


title: On Path-Complete Lyapunov Functions: Geometry and Comparison

abstract: We study optimization-based criteria for the stability of switching systems, known as Path-Complete Lyapunov Functions, and ask the question "can we decide algorithmically when a criterion is less conservative than another". Our contribution is twofold. First, we show that a Path-Complete Lyapunov Function, which is a multiple Lyapunov function by nature, can always be expressed as a common Lyapunov function taking the form of a combination of minima and maxima of the elementary functions that compose it. Geometrically, our results provide for each Path-Complete criterion an implied invariant set. Second, we provide a linear programming criterion allowing to compare the conservativeness of two arbitrary given Path-Complete Lyapunov functions.


title: Estimating Target-Object Shape Using Location-Unknown Mobile Fixed-Direction Distance Sensors

abstract: This paper proposes a method of estimating a target-object shape, the location of which is unknown, through the use of location-unknown mobile distance sensors. The direction of the sensor is fixed from the moving direction. Typically, mobile sensors are mounted on vehicles. Each sensor continuously measures the distance from it to the target object. The estimation method does not require any positioning function, anchor-location information, or additional mechanisms to obtain side information such as angle of arrival of signal. Under the assumption of a polygon target object, each edge length and vertex angle and their combinations are estimated to completely estimate the shape of the target object.


title: Real Seifert forms and polarizing forms of Steenbrink mixed Hodge structures

abstract: An isolated hypersurface singularity comes equipped with many different pairings on different spaces, the intersection form and the Seifert form on the Milnor lattice, a polarizing form for a mixed Hodge structure on a dual space, and a flat pairing on the cohomology bundle. This paper describes them and their relations systematically in an abstract setting. We expect applications also in other areas than singularity theory. A good part of the paper is elementary, but not well known: the classification of irreducible Seifert form pairs, the polarizing form on the generalized eigenspace with eigenvalue 1, an automorphism from a Fourier-Laplace transformation which involves the Gamma function and which relates Seifert form and polarizing form and a flat pairing on the cohomology bundle. New is a correction of a Thom-Sebastiani formula for Steenbrink's Hodge filtration in the case of singularities. It uses the Fourier-Laplace transformation. A special case is a square root of a Tate twist for Steenbrink mixed Hodge structures.


title: General Erased-Word Processes: Product-Type Filtrations, Ergodic Laws and Martin Boundaries

abstract: We study the dynamics of erasing randomly chosen letters from words by introducing a certain class of discrete-time stochastic processes, general erased-word processes(GEWPs), and investigating three closely related topics: Representation, Martin boundary and filtration theory. We use de Finetti's theorem and the random exchangeable linear order to obtain a de Finetti-type representation of GEWPs involving induced order statistics. Our studies expose connections between exchangeability theory and certain poly-adic filtrations that can be found in other exchangeable random objects as well. We show that ergodic GEWPs generate backward filtrations of product-type and by that generalize a result by S.Laurent.


title: Explicit formulas for heat kernels on diamond fractals

abstract: This paper provides explicit pointwise formulas for the heat kernel on compact metric measure spaces that belong to a $(\mathbb{N}\times\mathbb{N})$-parameter family of fractals which are regarded as projective limits of metric measure graphs and do not satisfy the volume doubling property. The formulas are applied to obtain uniform continuity estimates of the heat kernel and to derive an expression of the fundamental solution of the free Schr\"odinger equation. The results also open up the possibility to approach infinite dimensional spaces based on this model.


title: Probabilistic Adaptive Computation Time

abstract: We present a probabilistic model with discrete latent variables that control the computation time in deep learning models such as ResNets and LSTMs. A prior on the latent variables expresses the preference for faster computation. The amount of computation for an input is determined via amortized maximum a posteriori (MAP) inference. MAP inference is performed using a novel stochastic variational optimization method. The recently proposed Adaptive Computation Time mechanism can be seen as an ad-hoc relaxation of this model. We demonstrate training using the general-purpose Concrete relaxation of discrete variables. Evaluation on ResNet shows that our method matches the speed-accuracy trade-off of Adaptive Computation Time, while allowing for evaluation with a simple deterministic procedure that has a lower memory footprint.


title: Footprint and minimum distance functions

abstract: Let $S$ be a polynomial ring over a field $K$, with a monomial order $\prec$, and let $I$ be an unmixed graded ideal of $S$. In this paper we study two functions associated to $I$: the minimum distance function $\delta_I$ and the footprint function ${\rm fp}_I$. It is shown that $\delta_I$ is positive and that ${\rm fp}_I$ is positive if the initial ideal of $I$ is unmixed. Then we show that if $I$ is radical and its associated primes are generated by linear forms, then $\delta_I$ is strictly decreasing until it reaches the asymptotic value $1$. If $I$ is the edge ideal of a Cohen--Macaulay bipartite graph, we show that $\delta_I(d)=1$ for $d$ greater than or equal to the regularity of $S/I$. For a graded ideal of dimension $\geq 1$, whose initial ideal is a complete intersection, we give an exact sharp lower bound for the corresponding minimum distance function.


title: Conjectures on spectral numbers for upper triangular matrices and for singularities

abstract: Cecotti and Vafa proposed in 1993 a beautiful idea how to associate spectral numbers $\alpha_1,...,\alpha_n\in{\mathbb R}$ to real upper triangular $n\times n$ matrices $S$ with 1's on the diagonal and eigenvalues of $S^{-1}S^t$ in the unit sphere. Especially, $\exp(-2\pi i\alpha_j)$ shall be the eigenvalues of $S^{-1}S^t$. We tried to make their idea rigorous, but we succeeded only partially. This paper fixes our results and our conjectures. For certain subfamilies of matrices their idea works marvellously, and there the spectral numbers fit well to natural (split) polarized mixed Hodge structures. We formulate precise conjectures saying how this should extend to all matrices $S$ as above. The idea might become relevant in the context of semiorthogonal decompositions in derived algebraic geometry. Our main interest are the cases of Stokes like matrices which are associated to holomorphic functions with isolated singularities (Landau-Ginzburg models). Also there we formulate precise conjectures (which overlap with expectations of Cecotti and Vafa). In the case of the chain type singularities, we have positive results. We hope that this paper will be useful for further studies of the idea of Cecotti and Vafa.


title: Statistical properties of Fermi GBM GRBs' spectra

abstract: Statistical studies of gamma-ray burst (GRB) spectra may result in important information on the physics of GRBs. The Fermi GBM catalog contains GRB parameters (peak energy, spectral indices, intensity) estimated fitting the gamma-ray SED of the total emission (fluence, flnc), and during the time of the peak flux pflx. Using contingency tables we studied the relationship of the models best fitting pflx and flnc time intervals. Our analysis revealed an ordering of the spectra into a power law - Comptonized - smoothly broken power law - Band series. This result was further supported by a correspondence analysis (CA) of the pflx and flnc spectra categorical variables. We performed a linear discriminant analysis (LDA) to find a relationship between categorical (spectral) and model independent physical data. LDA resulted in highly significant physical differences among the spectral types, that is more pronounced in the case of the pflx spectra, than for the flnc spectra. We interpreted this difference as caused by the temporal variation of the spectrum during the outburst. This spectral variability is confirmed by the differences in the low energy spectral index and peak energy, between the pflx and flnc spectra. We found that the synchrotron radiation is significant in GBM spectra. The mean low energy spectral index is close to the canonical value of {\alpha} = -2/3 during the peak flux. However, $\alpha$ is ~ -0.9 for the spectra of the fluences. We interpret this difference as showing that the effect of cooling is important only for the fluence spectra.


title: Gain Scheduling LPV Control Scheme for the Autonomous Guidance Problem using a Dynamic Modelling Approach

abstract: This work proposes a solution for the longitudinal and lateral control problem of urban autonomous vehicles using a gain scheduling LPV control approach. Using the kinematic and dynamic vehicle models, a linear parameter varying (LPV) representation is adopted and a cascade control methodology is proposed for controlling both vehicle behaviours. In particular, for the control design, the use of both models separately lead to solve two LPV LMI-LQR problems. Furthermore, to achieve the desired levels of performance, an approach based on cascade design of the the kinematic and dynamic controllers has been proposed. This cascade control scheme is based on the idea that the dynamic closed loop behaviour is designed to be faster than the kinematic closed loop one. The obtained gain scheduling LPV control approach, jointly with a trajectory generation module, has presented suitable results in a simulated city driving scenario.


title: The Tightness of the Kesten-Stigum Reconstruction Bound of Symmetric Model with Multiple Mutations

abstract: It is well known that reconstruction problems, as the interdisciplinary subject, have been studied in numerous contexts including statistical physics, information theory and computational biology, to name a few. We consider a $2q$-state symmetric model, with two categories of $q$ states in each category, and 3 transition probabilities: the probability to remain in the same state, the probability to change states but remain in the same category, and the probability to change categories. We construct a nonlinear second order dynamical system based on this model and show that the Kesten-Stigum reconstruction bound is not tight when $q \geq 4$.


title: Electron localization properties in high pressure hydrogen at the liquid-liquid phase transition by Coupled Electron-Ion Monte Carlo

abstract: We analyze in detail the electronic properties of high pressure hydrogen around the liquid-liquid phase transition based on Coupled Electron-Ion Monte Carlo calculations. Computing the off-diagonal single particle density matrix and the momentum distribution we discuss localization properties of the electrons. The abrupt changes of these distributions indicate a metal to insulator transition occurring together with the structural transition from the atomic to molecular fluid. We further discuss the electron-proton and electron-electron pair correlation functions, which also change abruptly at the transition.


title: The Contribution of Faint, Failed and Defunct Stars to the "Stellar" Masses of Galaxies

abstract: A substantial fraction the stellar mass attributed to galaxies is invisible: stars close to the hydrogen burning limit, brown dwarfs, white dwarfs, neutron stars and black holes. These constituents do, however, gravitationally micro-lens background quasars, thereby permitting measurement of the total stellar contribution to the mass surface density along the line of sight. We report the results of such a measurement using a sample of ten quadruply lensed quasars. We discuss the prospects for improving upon this measurement with a larger sample and describe efforts to find new quadruple lenses. If we invert our argument and take the stellar mass to be known, we derive a value for the fraction of the dark halo in MaCHOs (including 20 solar mass primordial black holes) of somthing less than 10%, confirming the widely ignored result of Mediavilla et al (2009).


title: Chiral magnetic effect search in p+Au, d+Au and Au+Au collisions at RHIC

abstract: Metastable domains of fluctuating topological charges can change the chirality of quarks and induce local parity violation in quantum chromodynamics. This can lead to observable charge separation along the direction of the strong magnetic field produced by spectator protons in relativistic heavy-ion collisions, a phenomenon called the chiral magnetic effect (CME). A major background source for CME measurements using the charge-dependent azimuthal correlator ($\Delta\gamma$) is the intrinsic particle correlations (such as resonance decays) coupled with the azimuthal elliptical anisotropy ($v_{2}$). In heavy-ion collisions, the magnetic field direction and event plane angle are correlated, thus the CME and the $v_{2}$-induced background are entangled. In this report, we present two studies from STAR to shed further lights on the background issue. (1) The $\Delta\gamma$ should be all background in small system p+Au and d+Au collisions, because the event plane angles are dominated by geometry fluctuations uncorrelated to the magnetic field direction. However, significant $\Delta\gamma$ is observed, comparable to the peripheral Au+Au data, suggesting a background dominance in the latter, and likely also in the mid-central Au+Au collisions where the multiplicity and $v_{2}$ scaled correlator is similar. (2) A new approach is devised to study $\Delta\gamma$ as a function of the particle pair invariant mass ($m_{inv}$) to identify the resonance backgrounds and hence to extract the possible CME signal. Signal is consistent with zero within uncertainties at high $m_{inv}$. Signal at low $m_{inv}$, extracted from a two-component model assuming smooth mass dependence, is consistent with zero within uncertainties.


title: Reconstruction for the Asymmetric Ising Model on Regular Trees

abstract: It is known that the Kesten-Stigum reconstruction bound is tight for roughly symmetric binary channels. In this paper, we will adopt a refined analysis of moment recursion on a weighted version of the magnetization, which is engaged in Sly [2011] to handle the symmetric Potts model, and establish the critical condition of the asymmetric Ising model to make Kesten-Stigum bound the reconstruction threshold on regular $d$-ary trees.


title: Simulating the Coronal Evolution of AR 11437 using SDO/HMI Magnetograms

abstract: The coronal magnetic field evolution of AR 11437 is simulated by applying the magnetofrictional relaxation technique of Mackay et al. (2011). A sequence of photospheric line-of-sight magnetograms produced by SDO/HMI are used to drive the simulation and continuously evolve the coronal magnetic field of the active region through a series of non-linear force-free equilibria. The simulation is started during the first stages of the active region emergence so that its full evolution from emergence to decay can be simulated. A comparison of the simulation results with SDO/AIA observations show that many aspects of the active region's observed coronal evolution are reproduced. In particular, it shows the presence of a flux rope, which forms at the same location as sheared coronal loops in the observations. The observations show that eruptions occur on 2012 March 17 at 05:09 UT and 10:45 UT and on 2012 March 20 at 14:31 UT. The simulation reproduces the first and third eruption, with the simulated flux rope erupting roughly 1 and 10 hours before the observed ejections, respectively. A parameter study is conducted where the boundary and initial conditions are varied along with the physical effects of Ohmic diffusion, hyperdiffusion and an additional injection of helicity. When comparing the simulations, the evolution of the magnetic field, free magnetic energy, relative helicity and flux rope eruption timings do not change significantly. This indicates that the key element in reproducing the coronal evolution of AR 11437 is the use of line-of-sight magnetograms to drive the evolution of the coronal magnetic field.


title: Comparing experiments on quantum traversal time with the predictions of a space-time-symmetric formalism

abstract: The question of how long a particle takes to pass through a potential barrier is still a controversial topic in quantum mechanics. Arguably, the main theoretical problem in obtaining estimates for measurable times is the fact that previously defined time operators, that remained within the borders of standard quantum mechanics, present some kind of pathology. Recently, a time operator acting on an additional Hilbert space has been shown to support both Hermiticity and canonical relation with an energy observable. The theory is built in a framework which treats space and time as symmetrically as possible, in the nonrelativistic regime. In this work, we use this formalism to derive a closed analytic expression for the traversal time of a quantum particle impinging on a constant-potential barrier. We test our theory in the specific experimental scenario of a realization by Rafagni et al [Appl. Phys. Lett. {\bf 58}, 774 (1991)]. The proposed approach displays a much better performance in comparison with the B\"uttiker-Landauer and the phase-time approximations.


title: Inverse problem for the Helmholtz equation with Cauchy data: reconstruction with conditional well-posedness driven iterative regularization

abstract: In this paper, we study the performance of Full Waveform Inversion (FWI) from time-harmonic Cauchy data via conditional well-posedness driven iterative regularization. The Cauchy data can be obtained with dual sensors measuring the pressure and the normal velocity. We define a novel misfit functional which, adapted to the Cauchy data, allows the independent location of experimental and computational sources. The conditional well-posedness is obtained for a hierarchy of subspaces in which the inverse problem with partial data is Lipschitz stable. Here, these subspaces yield piecewise linear representations of the wave speed on given domain partitions. Domain partitions can be adaptively obtained through segmentation of the gradient. The domain partitions can be taken as a coarsening of an unstructured tetrahedral mesh associated with a finite element discretization of the Helmholtz equation. We illustrate the effectiveness of the iterative regularization through computational experiments with data in dimension three. In comparison with earlier work, the Cauchy data do not suffer from eigenfrequencies in the configurations.


title: Research assessment by percentile-based double rank analysis

abstract: In the double rank analysis of research publications, the local rank position of a country or institution publication is expressed as a function of the world rank position. Excluding some highly or lowly cited publications, the double rank plot fits well with a power law, which can be explained because citations for local and world publications follow lognormal distributions. We report here that the distribution of the number of country or institution publications in world percentiles is a double rank distribution that can be fitted to a power law. Only the data points in high percentiles deviate from it when the local and world $\mu$ parameters of the lognormal distributions are very different. The likelihood of publishing very highly cited papers can be calculated from the power law that can be fitted either to the upper tail of the citation distribution or to the percentile-based double rank distribution. The great advantage of the latter method is that it has universal application, because it is based on all publications and not just on highly cited publications. Furthermore, this method extends the application of the well-established percentile approach to very low percentiles where breakthroughs are reported but paper counts cannot be performed.


title: Near critical preferential attachment networks have small giant components

abstract: Preferential attachment networks with power law exponent $\tau>3$ are known to exhibit a phase transition. There is a value $\rho_{\rm c}>0$ such that, for small edge densities $\rho\leq \rho_c$ every component of the graph comprises an asymptotically vanishing proportion of vertices, while for large edge densities $\rho>\rho_c$ there is a unique giant component comprising an asymptotically positive proportion of vertices. In this paper we study the decay in the size of the giant component as the critical edge density is approached from above. We show that the size decays very rapidly, like $\exp(-c/ \sqrt{\rho-\rho_c})$ for an explicit constant $c>0$ depending on the model implementation. This result is in contrast to the behaviour of the class of rank-one models of scale-free networks, including the configuration model, where the decay is polynomial. Our proofs rely on the local neighbourhood approximations of [Dereich, Morters, 2013] and recent progress in the theory of branching random walks [Gantert, Hu, Shi, 2011].


title: On the martingale decompositions of Gundy, Meyer, and Yoeurp in infinite dimensions

abstract: We show that the canonical decomposition (comprising both the Meyer-Yoeurp and the Yoeurp decompositions) of a general $X$-valued local martingale is possible if and only if $X$ has the UMD property. More precisely, $X$ is a UMD Banach space if and only if for any $X$-valued local martingale $M$ there exist a continuous local martingale $M^c$, a purely discontinuous quasi-left continuous local martingale $M^q$, and a purely discontinuous local martingale $M^a$ with accessible jumps such that $M = M^c + M^q + M^a$. The corresponding weak $L^1$-estimates are provided. Important tools used in the proof are a new version of Gundy's decomposition of continuous-time martingales and weak $L^1$-bounds for a certain class of vector-valued continuous-time martingale transforms.


title: Super Periodic Potential

abstract: In this paper we introduce the concept of super periodic potential (SPP) of arbitrary order $n$, $n \in I^{+}$, in one dimension. General theory of wave propagation through SPP of order $n$ is presented and the reflection and transmission coefficients are derived in their closed analytical form by transfer matrix formulation. We present scattering features of super periodic rectangular potential and super periodic delta potential as special cases of SPP. It is found that the symmetric self-similarity is the special case of super periodicity. Thus by identifying a symmetric fractal potential as special cases of SPP, one can obtain the tunnelling amplitude for a particle from such fractal potential. By using the formalism of SPP we obtain the close form expression of tunnelling amplitude of a particle for general Cantor and Smith-Volterra-Cantor potentials.


title: Robust isogeometric preconditioners for the Stokes system based on the Fast Diagonalization method

abstract: In this paper we propose a new class of preconditioners for the isogeometric discretization of the Stokes system. Their application involves the solution of a Sylvester-like equation, which can be done efficiently thanks to the Fast Diagonalization method. These preconditioners are robust with respect to both the spline degree and mesh size. By incorporating information on the geometry parametrization and equation coefficients, we maintain efficiency on non-trivial computational domains and for variable kinematic viscosity. In our numerical tests we compare to a standard approach, showing that the overall iterative solver based on our preconditioners is significantly faster.


title: Observing and Tracking Bandlimited Graph Processes

abstract: One of the most crucial challenges in graph signal processing is the sampling of bandlimited graph signals, i.e., signals that are sparse in a well-defined graph Fourier domain. So far, the prior art is mostly focused on (sub)sampling single snapshots of graph signals ignoring their evolution over time. However, time can bring forth new insights, since many real signals like sensor measurements, biological, financial, and network signals in general, have intrinsic correlations in both domains. In this work, {we fill this lacuna} by jointly considering the graph-time nature of graph signals, named \emph{graph processes} for two main tasks: \emph{i)} observability of graph processes; and \emph{ii)} tracking of graph processes via Kalman filtering; both from a (possibly time-varying) subset of nodes. A detailed mathematical analysis ratifies the proposed methods and provides insights into the role played by the different actors, such as the graph topology, the process bandwidth, and the sampling strategy. Moreover, (sub)optimal sampling strategies that jointly exploit the nature of the graph structure and graph process are proposed. Several numerical tests on both synthetic and real data validate our theoretical findings and illustrate the performance of the proposed methods in coping with time-varying graph signals.


title: The Dirichlet Problem for the Complex Homogeneous Monge-Amp\`ere Equation

abstract: We survey the Dirichlet problem for the complex Homogeneous Monge-Amp\`ere Equation, both in the case of domains in $\mathbb C^n$ and the case of compact K\"ahler manifolds parametrized by a Riemann surface with boundary. We then give a self-contained account of previous work of the authors that connects this with the Hele-Shaw flow, and give several concrete examples illustrating various phenomena that solutions to this problem can display.


title: On the accuracy of Monte Carlo based beam dynamics models for the degrader in proton therapy facilities

abstract: In a cyclotron-based proton therapy facility, the energy changes are performed by means of a degrader of variable thickness. The interaction of the proton beam with the degrader creates energy tails and increases the beam emittance. A precise model of the degraded beam properties is important not only to better understand the performance of a facility already in operation, but also to support the development of new proton therapy concepts. The exact knowledge of the degraded beam properties, in terms of energy spectrum and transverse phase space, depends on the model used to describe the proton interaction with the degrader material. In this work the model of a graphite degrader has been developed with four Monte Carlo codes: three conventional Monte Carlo codes (FLUKA, GEANT4 and MCNPX) and the multi-purpose particle tracking code OPAL equipped with a simplified Monte Carlo routine. From the comparison between the different codes, we can deduce how the accuracy of the degrader model influences the precision of the beam dynamics model of a possible transport line downstream of the degrader.


title: Introduction to Protein Structure Prediction

abstract: This chapter gives a graceful introduction to problem of protein three- dimensional structure prediction, and focuses on how to make structural sense out of a single input sequence with unknown structure, the 'query' or 'target' sequence. We give an overview of the different classes of modelling techniques, notably template-based and template free. We also discuss the way in which structural predictions are validated within the global com- munity, and elaborate on the extent to which predicted structures may be trusted and used in practice. Finally we discuss whether the concept of a sin- gle fold pertaining to a protein structure is sustainable given recent insights. In short, we conclude that the general protein three-dimensional structure prediction problem remains unsolved, especially if we desire quantitative predictions. However, if a homologous structural template is available in the PDB model or reasonable to high accuracy may be generated.


title: Binarized octree generation for Cartesian adaptive mesh refinement around immersed geometries

abstract: We revisit the generation of balanced octrees for adaptive mesh refinement (AMR) of Cartesian domains with immersed complex geometries. In a recent short note [Hasbestan and Senocak, J. Comput. Phys. vol. 351:473-477 (2017)], we showed that the data-locality of the Z-order curve in hashed linear octree generation methods may not be perfect because of potential collisions in the hash table. Building on that observation, we propose a binarized octree generation method that complies with the Z-order curve exactly. Similar to a hashed linear octree generation method, we use Morton encoding to index the nodes of an octree, but use a red-black tree in place of the hash table. Red-black tree is a special kind of a binary tree, which we use for insertion and deletion of elements during mesh adaptation. By strictly working with the bitwise representation of the octree, we remove computer hardware limitations on the depth of adaptation on a single processor. Additionally, we introduce a geometry encoding technique for rapidly tagging the solid geometry for refinement. Our results for several geometries with different levels of adaptations show that the binarized octree generation outperforms the linear octree generation in terms of runtime performance at the expense of only a slight increase in memory usage. We provide the current AMR capability as open-source software.


title: Deep Learning Scaling is Predictable, Empirically

abstract: Deep learning (DL) creates impactful advances following a virtuous recipe: model architecture search, creating large training data sets, and scaling computation. It is widely believed that growing training sets and models should improve accuracy and result in better products. As DL application domains grow, we would like a deeper understanding of the relationships between training set size, computational scale, and model accuracy improvements to advance the state-of-the-art. This paper presents a large scale empirical characterization of generalization error and model size growth as training sets grow. We introduce a methodology for this measurement and test four machine learning domains: machine translation, language modeling, image processing, and speech recognition. Our empirical results show power-law generalization error scaling across a breadth of factors, resulting in power-law exponents---the "steepness" of the learning curve---yet to be explained by theoretical work. Further, model improvements only shift the error but do not appear to affect the power-law exponent. We also show that model size scales sublinearly with data size. These scaling relationships have significant implications on deep learning research, practice, and systems. They can assist model debugging, setting accuracy targets, and decisions about data set growth. They can also guide computing system design and underscore the importance of continued computational scaling.


title: On the few products, many sums problem

abstract: We prove new results on additive properties of finite sets $A$ with small multiplicative doubling $|AA|\leq M|A|$ in the category of real/complex sets as well as multiplicative subgroups in the prime residue field. The improvements are based on new combinatorial lemmata, which may be of independent interest. Our main results are the inequality $$ |A-A|^3|AA|^5 \gtrsim |A|^{10}, $$ over the reals, "redistributing" the exponents in the textbook Elekes sum-product inequality and the new best known additive energy bound $\mathsf E(A)\lesssim_M |A|^{49/20}$, which aligns, in a sense to be discussed, with the best known sum set bound $|A+A|\gtrsim_M |A|^{8/5}$. These bounds, with $M=1$, also apply to multiplicative subgroups of $\mathbb F^\times_p$, whose order is $O(\sqrt{p})$. We adapt the above energy bound to larger subgroups and obtain new bounds on gaps between elements in cosets of subgroups of order $\Omega(\sqrt{p})$.


title: Spectroscopy of $^{54}$Ti and the systematic behavior of low energy octupole states in Ca and Ti isotopes

abstract: Excited states of the $N=32$ nucleus $^{54}$Ti have been studied, via both inverse-kinematics proton scattering and one-neutron knockout from $^{55}$Ti by a liquid hydrogen target, using the GRETINA $\gamma$-ray tracking array. Inelastic proton-scattering cross sections and deformation lengths have been determined. A low-lying octupole state has been tentatively identified in $^{54}$Ti for the first time. A comparison of $(p,p')$ results on low-energy octupole states in the neutron-rich Ca and Ti isotopes with the results of Random Phase Approximation calculations demonstrates that the observed systematic behavior of these states is unexpected.


title: ABJM quantum spectral curve and Mellin transform

abstract: The present techniques for the perturbative solution of quantum spectral curve problems in N=4 SYM and ABJM models are limited to the situation when the states quantum numbers are given explicitly as some integer numbers. These techniques are sufficient to recover full analytical structure of the conserved charges provided that we know a finite basis of functions in terms of which they could be written explicitly. It is known that in the case of N=4 SYM both the contributions of asymptotic Bethe ansatz and wrapping or finite size corrections are expressed in terms of the harmonic sums. However, in the case of ABJM model only the asymptotic contribution can still be written in the harmonic sums basis, while the wrapping corrections part can not. Moreover, the generalization of harmonic sums basis for this problem is not known. In this paper we present a Mellin space technique for the solution of multiloop Baxter equations, which is the main ingredient for the solution of corresponding quantum spectral problems, and provide explicit results for the solution of ABJM quantum spectral curve in the case of twist 1 operators in sl(2) sector for arbitrary spin values up to four loop order with explicit account for wrapping corrections. It is shown that the result for anomalous dimensions could be expressed in terms of harmonic sums decorated by the fourth root of unity factors, so that maximum transcendentality principle holds.


title: Robust and Efficient Modular Grad-Div Stabilization

abstract: This paper presents two modular grad-div algorithms for calculating solutions to the Navier-Stokes equations (NSE). These algorithms add to an NSE code a minimally intrusive module that implements grad-div stabilization. The algorithms do not suffer from either breakdown (locking) or debilitating slow down for large values of grad-div parameters. Stability and optimal-order convergence of the methods are proven. Numerical tests confirm the theory and illustrate the benefits of these algorithms over a fully coupled grad-div stabilization.


title: Slime mould: the fundamental mechanisms of cognition

abstract: The slime mould Physarum polycephalum has been used in developing unconventional computing devices for in which the slime mould played a role of a sensing, actuating, and computing device. These devices treated the slime mould rather as an active living substrate yet the slime mould is a self-consistent living creature which evolved for millions of years and occupied most part of the world, but in any case, that living entity did not own true cognition, just automated biochemical mechanisms. To "rehabilitate" the slime mould from the rank of a purely living electronics element to a "creature of thoughts" we are analyzing the cognitive potential of P. polycephalum. We base our theory of minimal cognition of the slime mould on a bottom-up approach, from the biological and biophysical nature of the slime mould and its regulatory systems using frameworks suh as Lyon's biogenic cognition, Muller, di Primio-Lengeler\'s modifiable pathways, Bateson's "patterns that connect" framework, Maturana's autopoetic network, or proto-consciousness and Morgan's Canon.


title: Possible detection of a bimodal cloud distribution in the atmosphere of HAT-P-32\,A\,b from multi-band photometry

abstract: We present high-precision photometry of eight separate transit events in the HAT-P-32 planetary system. One transit event was observed simultaneously by two telescopes of which one obtained a simultaneous multi-band light curve in three optical bands, giving a total of 11 transit light curves. Due to the filter selection and in conjunction with using the defocussed photometry technique we were able to obtain an extremely high precision, ground-based transit in the \textit{u}-band (350\,nm), with an rms scatter of $\approx 1$\,mmag. All 11 transits were modelled using \textsc{prism} and \textsc{gemc}, and the physical properties of the system calculated. We find the mass and radius of the host star to be $1.182\pm 0.041\Msun$ and $1.225\pm0.015\Rsun$, respectively. For the planet we find a mass of $0.80\pm 0.14\Mjup$, a radius of $1.807\pm0.022\Rjup$ and a density of $0.126\pm0.023\pjup$. These values are consistent with those found in the literature. We also obtain a new orbital ephemeris for the system $ T_0 = {\rm BJD/TDB} \,\, 2\,454\,420.447187 (96) \, + \,2.15000800 (10) \times E $. We measured the transmission spectrum of HAT-P-32\,A\,b and compared it to theoretical transmission spectra. Our results indicate a bimodal cloud particle distribution consisting of Rayleigh--like haze and grey absorbing cloud particles within the atmosphere of HAT-P-32\,A\,b.


title: Reactive control of the dynamics of a fully turbulent wake using real-time PIV

abstract: In this study we focus on the control of the dynamics of 3D turbulent wake downstream a square-back Ahmed body ($Re_H=3.9\times10^5$). The peculiar dynamics of such a wake are first characterized through the trajectories of the pressure barycenter over the rear part of the model as well as the recirculation barycenter in the wake. In particular it is shown that these dynamics allow the definition of three different states: the two so-called reflectional symmetry-breaking (RSB) modes and the transient symmetric (TS) mode. It was shown recently that the time-fluctuations of the pressure barycenter could be characterized as a weak chaotic system with a well-defined attractor. We show that the dynamics of the bimodal wake can then be forced into a stable asymmetric or symmetric state in open loop control, using tangential continuous or pulsed blowing in three different regions along the upper edge of the rear part of the model. Finally, a simple closed-loop opposition control, based on real-time identification of the wake barycenter in the PIV fields, is used to force the chaotic dynamics of the wake into a regular oscillatory motion at a well-controlled frequency. Depending on the actuation parameters, the wake dynamics can also be switched from bimodal to a new multimodal behavior. We show that this new mode also exhibits a peculiar dynamics with an up-down instead of left-right chaotic oscillations. Interestingly, the recirculation area (size of the recirculation bubble) is much more reduced for the closed-loop experiments when the jets are pulsed rather than continuous. For the pulsed jets, the reduction is also increased when the proper frequency is chosen.


title: Stability and Dynamics of Microring Combs: Elliptic function solutions of the Lugiato-Lefever equation

abstract: We consider a new class of periodic solutions to the Lugiato-Lefever equations (LLE) that govern the electromagnetic field in a microresonator cavity. Specifically, we rigorously characterize the stability and dynamics of the Jacobi elliptic function solutions of LLE and show that the dn solution is stabilized by the pumping of the microresonator. In analogy with soliton perturbation theory, we also derive a microcomb perturbation theory that allows one to consider the effects of physically realizable perturbations on the comb line stability, including effects of Raman scattering and stimulated emission. Our results are verified through full numerical simulations of the LLE cavity dynamics. The perturbation theory gives a simple analytic platform for potentially engineering new resonator designs.


title: Transient Survey Rates for Orphan Afterglows from Compact Merger Jets

abstract: Orphan afterglows from short $\gamma$-ray bursts (GRB) are potential candidates for electromagnetic (EM) counterpart searches to gravitational wave (GW) detected neutron star or neutron star black hole mergers. Various jet dynamical and structure models have been proposed that can be tested by the detection of a large sample of GW-EM counterparts. We make predictions for the expected rate of optical transients from these jet models for future survey telescopes, without a GW or GRB trigger. A sample of merger jets is generated in the redshift limits $0\leq z\leq 3.0$, and the expected peak $r$-band flux and timescale above the LSST or ZTF detection threshold, $m_r=24.5$ and $20.4$ respectively, is calculated. General all-sky rates are shown for $m_r\leq26.0$ and $m_r\leq21.0$. The detected orphan and GRB afterglow rate depends on jet model, typically $16\lesssim R\lesssim 76$ yr$^{-1}$ for the LSST, and $2\lesssim R \lesssim 8$ yr$^{-1}$ for ZTF. An excess in the rate of orphan afterglows for a survey to a depth of $m_r\leq26$ would indicate that merger jets have a dominant low-Lorentz factor population, or the jets exhibit intrinsic jet structure. Careful filtering of transients is required to successfully identify orphan afterglows from either short or long GRB progenitors.


title: Significant improvement in estimation of the Young's modulus of the metal foam based on the Timoshenko's bend theory

abstract: A computational and experimental approach based on a natural vibration of a free prismatic thick beam with square cross-section is suggested. Three variants of the beam sample were used (one with skin and two without skin). From 13 to 16 lowest resonant frequencies of longitudinal, torsional and flexural vibration of each beam were analyzed. A rule for dependence of Young's modulus on the average foam density was derived from the sample without skin. It is shown that the skin presence causes known anisotropy that the stiffness in transverse direction is about 50% greater than that in the longitudinal direction. It has also been shown that the seeming frequency dependence of the Young's module can be explained by non-uniform distribution of mass density in the sample. Agreement among experimental and numerical data is excellent in most cases. The rule is also verified on solution of two lowest resonant frequencies of the free foam plates with various mass densities (from 400 to 1350 kg/m3). Agreement among experimental and numerical frequencies is acceptable.


title: Leptogenesis with TeV Scale $W_R$

abstract: Successful leptogenesis within the conventional TeV-scale left-right implementation of type-I seesaw has been shown to require that the mass of the right-handed $W_R^\pm$ boson should have a lower bound much above the reach of the Large Hadron Collider. This bound arises from the necessity to suppress the washout of lepton asymmetry due to $W_R^\pm$-mediated $\Delta L\neq 0$ processes. We show that in an alternative quark seesaw realization of left-right symmetry, the above bound can be avoided. Lepton asymmetry in this model is generated not via the usual right-handed neutrino decay but rather via the decay of new heavy scalars producing an asymmetry in the $B-L$ carrying Higgs triplets responsible for type-II seesaw, whose subsequent decay leads to the lepton asymmetry. This result implies that any evidence for $W_R$ at the LHC 14 will point towards this alternative realization of left-right symmetry, which is also known to solve the strong CP problem.


title: Hadronic light-by-light scattering amplitudes from lattice QCD versus dispersive sum rules

abstract: The hadronic contribution to the eight forward amplitudes of light-by-light scattering ($\gamma^*\gamma^*\to \gamma^*\gamma^*$) is computed in lattice QCD. Via dispersive sum rules, the amplitudes are compared to a model of the $\gamma^*\gamma^*\to {\rm hadrons}$ cross sections in which the fusion process is described by hadronic resonances. Our results thus provide an important test for the model estimates of hadronic light-by-light scattering in the anomalous magnetic moment of the muon, $a_\mu^{\rm HLbL}$. Using simple parametrizations of the resonance $M\to \gamma^*\gamma^*$ transition form factors, we determine the corresponding monopole and dipole masses by performing a global fit to all eight amplitudes. Together with a previous dedicated calculation of the $\pi^0\to \gamma^*\gamma^*$ transition form factor, our calculation provides valuable information for phenomenological estimates of $a_\mu^{\rm HLbL}$. The presented calculations are performed in two-flavor QCD with pion masses extending down to 190\,MeV at two different lattice spacings. In addition to the fully connected Wick contractions, on two lattice ensembles we also compute the (2+2) disconnected class of diagrams, and find that their overall size is compatible with a parameter-free, large-$N$ inspired prediction, where $N$ is the number of colors. Motivated by this observation, we estimate in the same way the disconnected contribution to $a_\mu^{\rm HLbL}$.


title: The AFLOW Fleet for Materials Discovery

abstract: The traditional paradigm for materials discovery has been recently expanded to incorporate substantial data driven research. With the intent to accelerate the development and the deployment of new technologies, the AFLOW Fleet for computational materials design automates high-throughput first principles calculations, and provides tools for data verification and dissemination for a broad community of users. AFLOW incorporates different computational modules to robustly determine thermodynamic stability, electronic band structures, vibrational dispersions, thermo-mechanical properties and more. The AFLOW data repository is publicly accessible online at aflow.org, with more than 1.7 million materials entries and a panoply of queryable computed properties. Tools to programmatically search and process the data, as well as to perform online machine learning predictions, are also available.


title: DAOS for Extreme-scale Systems in Scientific Applications

abstract: Exascale I/O initiatives will require new and fully integrated I/O models which are capable of providing straightforward functionality, fault tolerance and efficiency. One solution is the Distributed Asynchronous Object Storage (DAOS) technology, which is primarily designed to handle the next generation NVRAM and NVMe technologies envisioned for providing a high bandwidth/IOPS storage tier close to the compute nodes in an HPC system. In conjunction with DAOS, the HDF5 library, an I/O library for scientific applications, will support end-to-end data integrity, fault tolerance, object mapping, index building and querying. This paper details the implementation and performance of the HDF5 library built over DAOS by using three representative scientific application codes.


title: The reparameterization trick for acquisition functions

abstract: Bayesian optimization is a sample-efficient approach to solving global optimization problems. Along with a surrogate model, this approach relies on theoretically motivated value heuristics (acquisition functions) to guide the search process. Maximizing acquisition functions yields the best performance; unfortunately, this ideal is difficult to achieve since optimizing acquisition functions per se is frequently non-trivial. This statement is especially true in the parallel setting, where acquisition functions are routinely non-convex, high-dimensional, and intractable. Here, we demonstrate how many popular acquisition functions can be formulated as Gaussian integrals amenable to the reparameterization trick and, ensuingly, gradient-based optimization. Further, we use this reparameterized representation to derive an efficient Monte Carlo estimator for the upper confidence bound acquisition function in the context of parallel selection.


title: Strategies for protein structure model generation

abstract: This chapter deals with approaches for protein three-dimensional structure prediction, starting out from a single input sequence with unknown struc- ture, the 'query' or 'target' sequence. Both template based and template free modelling techniques are treated, and how resulting structural models may be selected and refined. We give a concrete flowchart for how to de- cide which modelling strategy is best suited in particular circumstances, and which steps need to be taken in each strategy. Notably, the ability to locate a suitable structural template by homology or fold recognition is crucial; without this models will be of low quality at best. With a template avail- able, the quality of the query-template alignment crucially determines the model quality. We also discuss how other, courser, experimental data may be incorporated in the modelling process to alleviate the problem of missing template structures. Finally, we discuss measures to predict the quality of models generated.


title: Extreme Scattering Events Towards Two Young Pulsars

abstract: We have measured the scintillation properties of 151 young, energetic pulsars with the Parkes radio telescope and have identified two extreme scattering events (ESEs). Towards PSR J1057-5226 we discovered a three-year span of strengthened scattering during which the variability in flux density and the scintillation bandwidth decreased markedly. The transverse size of the scattering region is ~23 au, and strong flux density enhancement before and after the ESE may arise from refractive focusing. Long observations reveal scintillation arcs characteristic of interference between rays scattered at large angles, and the clearest arcs appear during the ESE. The arcs suggest scattering by a screen 100-200 pc from the earth, perhaps ionized filamentary structure associated with the boundary of the local bubble(s). Towards PSR J1740-3015 we observed a "double dip" in the measured flux density similar to ESEs observed towards compact extragalactic radio sources. The observed shape is consistent with that produced by a many-au scale diverging plasma lens with electron density ~500cm$^{-3}$ . The continuing ESE is at least 1500 d long, making it the longest detected event to date. These detections, with materially different observational signatures, indicate that well-calibrated pulsar monitoring is a keen tool for ESE detection and ISM diagnostics. They illustrate the strong r\^ole au-scale non-Kolmogorov density fluctuations and the local ISM structure play in such events and are key to understanding both their intrinsic physics and their impact on other phenomena, particularly fast radio bursts.


title: Unsupervised Classification of PolSAR Data Using a Scattering Similarity Measure Derived from a Geodesic Distance

abstract: In this letter, we propose a novel technique for obtaining scattering components from Polarimetric Synthetic Aperture Radar (PolSAR) data using the geodesic distance on the unit sphere. This geodesic distance is obtained between an elementary target and the observed Kennaugh matrix, and it is further utilized to compute a similarity measure between scattering mechanisms. The normalized similarity measure for each elementary target is then modulated with the total scattering power (Span). This measure is used to categorize pixels into three categories i.e. odd-bounce, double-bounce and volume, depending on which of the above scattering mechanisms dominate. Then the maximum likelihood classifier of [J.-S. Lee, M. R. Grunes, E. Pottier, and L. Ferro-Famil, Unsupervised terrain classification preserving polarimetric scattering characteristics, IEEE Trans. Geos. Rem. Sens., vol. 42, no. 4, pp. 722731, April 2004.] based on the complex Wishart distribution is iteratively used for each category. Dominant scattering mechanisms are thus preserved in this classification scheme. We show results for L-band AIRSAR and ALOS-2 datasets acquired over San Francisco and Mumbai, respectively. The scattering mechanisms are better preserved using the proposed methodology than the unsupervised classification results using the Freeman-Durden scattering powers on an orientation angle (OA) corrected PolSAR image. Furthermore, (1) the scattering similarity is a completely non-negative quantity unlike the negative powers that might occur in double- bounce and odd-bounce scattering component under Freeman Durden decomposition (FDD), and (2) the methodology can be extended to more canonical targets as well as for bistatic scattering.


title: Novel Exploration Techniques (NETs) for Malaria Policy Interventions

abstract: The task of decision-making under uncertainty is daunting, especially for problems which have significant complexity. Healthcare policy makers across the globe are facing problems under challenging constraints, with limited tools to help them make data driven decisions. In this work we frame the process of finding an optimal malaria policy as a stochastic multi-armed bandit problem, and implement three agent based strategies to explore the policy space. We apply a Gaussian Process regression to the findings of each agent, both for comparison and to account for stochastic results from simulating the spread of malaria in a fixed population. The generated policy spaces are compared with published results to give a direct reference with human expert decisions for the same simulated population. Our novel approach provides a powerful resource for policy makers, and a platform which can be readily extended to capture future more nuanced policy spaces.


title: Event-Triggered Communication and Control of Networked Systems for Multi-Agent Consensus

abstract: This article provides an introduction to event-triggered coordination for multi-agent average consensus. We provide a comprehensive account of the motivations behind the use of event-triggered strategies for consensus, the methods for algorithm synthesis, the technical challenges involved in establishing desirable properties of the resulting implementations, and their applications in distributed control. We pay special attention to the assumptions on the capabilities of the network agents and the resulting features of the algorithm execution, including the interconnection topology, the evaluation of triggers, and the role of imperfect information. The issues raised in our discussion transcend the specific consensus problem and are indeed characteristic of cooperative algorithms for networked systems that solve other coordination tasks. As our discussion progresses, we make these connections clear, highlighting general challenges and tools to address them widespread in the event-triggered control of networked systems.


title: Building the peanut: simulations and observations of peanut-shaped structures and ansae in face-on disk galaxies

abstract: (X/peanut)-shaped features observed in a significant fraction of disk galaxies are thought to have formed from vertically buckled bars. Despite being three dimensional structures, they are preferentially detected in near edge-on projection. Only a few galaxies are found to have displayed such structures when their disks are relatively face-on - suggesting that either they are generally weak in face-on projection or many may be hidden by the light of their galaxy's face-on disk. Here we report on three (collisionless) simulated galaxies displaying peanut-shaped structures when their disks are seen both face-on and edge-on - resembling a three-dimensional peanut or dumbbell. Furthermore, these structures are accompanied by ansae and an outer ring at the end of the bar --- as seen in real galaxies such as IC~5240. The same set of quantitative parameters used to measure peanut structures in real galaxies have been determined for the simulated galaxies, and a broad agreement is found. In addition, the peanut length grows in tandem with the bar, and is a maximum at half the length of the bar. Beyond the cutoff of these peanut structures, towards the end of the bar, we discover a new positive/negative feature in the $B_6$ radial profile associated with the isophotes of the ansae/ring. Our simulated, self-gravitating, three-dimensional peanut structures display cylindrical rotation even in the near-face-on disk projection. In addition, we report on a kinematic pinch in the velocity map along the bar minor-axis, matching that seen in the surface density map.


title: Integrable Trotterization: Local Conservation Laws and Boundary Driving

abstract: We discuss a general procedure to construct an integrable real-time trotterization of interacting lattice models. As an illustrative example we consider a spin-$1/2$ chain, with continuous time dynamics described by the isotropic ($XXX$) Heisenberg Hamiltonian. For periodic boundary conditions local conservation laws are derived from an inhomogeneous transfer matrix and a boost operator is constructed. In the continuous time limit these local charges reduce to the known integrals of motion of the Heisenberg chain. In a simple Kraus representation we also examine the nonequilibrium setting, where our integrable cellular automaton is driven by stochastic processes at the boundaries. We show explicitly, how an exact nonequilibrium steady state density matrix can be written in terms of a staggered matrix product ansatz. This simple trotterization scheme, in particular in the open system framework, could prove to be a useful tool for experimental simulations of the lattice models in terms of trapped ion and atom optics setups.

categories: cond-mat.stat-mech math-ph math.MP nlin.CG nlin.SI


title: Mapping Jet-ISM Interactions in X-ray Binaries with ALMA: A GRS 1915$+$105 Case Study

abstract: We present Atacama Large Millimetre/Sub-Millimetre Array (ALMA) observations of IRAS 19132+1035, a candidate jet-ISM interaction zone near the black hole X-ray binary (BHXB) GRS 1915+105. With these ALMA observations (combining data from the 12 m array and the Atacama Compact Array), we map the molecular line emission across the IRAS 19132+1035 region. We detect emission from the $^{12}$CO [$J=2-1$], $^{13}$CO [$\nu=0$, $J=2-1$], C$^{18}$O [$J=2-1$], ${\rm H}_{2}{\rm CO}$ [$J=3_{0,3}-2_{0,2}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,2}-2_{2,1}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,1}-2_{2,0}$], SiO [$\nu=0$, $J=5-4$], CH$_3$OH [$J=4_{2,2}-3_{1,2}$], and CS [$\nu=0$, $J=5-4$] transitions. Given the morphological, spectral, and kinematic properties of this molecular emission, we present several lines of evidence that support the presence of a jet-ISM interaction at this site, including a jet-blown cavity in the molecular gas. This compelling new evidence identifies this site as a jet-ISM interaction zone, making GRS 1915$+$105 the third Galactic BHXB with at least one conclusive jet-ISM interaction zone. However, we find that this interaction occurs on much smaller scales than was postulated by previous work, where the BHXB jet does not appear to be dominantly powering the entire IRAS 19132+1035 region. Using estimates of the ISM conditions in the region, we utilize the detected cavity as a calorimeter to estimate the time-averaged power carried in the GRS 1915+105 jets of $(8.4^{+7.7}_{-8.1})\times10^{32}\,{\rm erg\,s}^{-1}$. Overall, our analysis demonstrates that molecular lines are excellent diagnostic tools to identify and probe jet-ISM interaction zones near Galactic BHXBs.


title: Single-Shot Object Detection with Enriched Semantics

abstract: We propose a novel single shot object detection network named Detection with Enriched Semantics (DES). Our motivation is to enrich the semantics of object detection features within a typical deep detector, by a semantic segmentation branch and a global activation module. The segmentation branch is supervised by weak segmentation ground-truth, i.e., no extra annotation is required. In conjunction with that, we employ a global activation module which learns relationship between channels and object classes in a self-supervised manner. Comprehensive experimental results on both PASCAL VOC and MS COCO detection datasets demonstrate the effectiveness of the proposed method. In particular, with a VGG16 based DES, we achieve an mAP of 81.7 on VOC2007 test and an mAP of 32.8 on COCO test-dev with an inference speed of 31.5 milliseconds per image on a Titan Xp GPU. With a lower resolution version, we achieve an mAP of 79.7 on VOC2007 with an inference speed of 13.0 milliseconds per image.


title: On the treewidth of triangulated 3-manifolds

abstract: In graph theory, as well as in 3-manifold topology, there exist several width-type parameters to describe how "simple" or "thin" a given graph or 3-manifold is. These parameters, such as pathwidth or treewidth for graphs, or the concept of thin position for 3-manifolds, play an important role when studying algorithmic problems; in particular, there is a variety of problems in computational 3-manifold topology - some of them known to be computationally hard in general - that become solvable in polynomial time as soon as the dual graph of the input triangulation has bounded treewidth. In view of these algorithmic results, it is natural to ask whether every 3-manifold admits a triangulation of bounded treewidth. We show that this is not the case, i.e., that there exists an infinite family of closed 3-manifolds not admitting triangulations of bounded pathwidth or treewidth (the latter implies the former, but we present two separate proofs). We derive these results from work of Agol, of Scharlemann and Thompson, and of Scharlemann, Schultens and Saito by exhibiting explicit connections between the topology of a 3-manifold M on the one hand and width-type parameters of the dual graphs of triangulations of M on the other hand, answering a question that had been raised repeatedly by researchers in computational 3-manifold topology. In particular, we show that if a closed, orientable, irreducible, non-Haken 3-manifold M has a triangulation of treewidth (resp. pathwidth) k then the Heegaard genus of M is at most 18(k+1) (resp. 4(3k+1)).


title: Coherent spin dynamics of ytterbium ions in yttrium orthosilicate

abstract: We investigate the electron and nuclear spin coherence properties of ytterbium ($\mathrm{Yb}^{3+}$) ions with non-zero nuclear spin, within an yttrium orthosilicate (Y$_2$SiO$_5$) crystal, with a view to their potential application in quantum memories or repeaters. We find electron spin-lattice relaxation times are maximised at low magnetic field ($<100$ mT) where $g~\sim6$, reaching 5 s at 2.5 K, while coherence times are maximised when addressing ESR transitions at higher fields where $g\sim0.7$ where a Hahn echo measurement yields $T_2$ up to 73 $\mu$s. Dynamical decoupling (XY16) can be used to suppress spectral diffusion and extend the coherence lifetime to over 0.5 ms, close to the limit of instantaneous diffusion. Using Davies electron-nuclear-double-resonance (ENDOR), we performed coherent control of the $^{173}\mathrm{Yb}^{3+}$ nuclear spin and studied its relaxation dynamics. At around 4.5 K we measure a nuclear spin $T_1$ and $T_2$ of 4 and 0.35 ms, respectively, about 4 and 14 times longer than the corresponding times for the electron spin.


title: Unsupervised Learning for Color Constancy

abstract: Most digital camera pipelines use color constancy methods to reduce the influence of illumination and camera sensor on the colors of scene objects. The highest accuracy of color correction is obtained with learning-based color constancy methods, but they require a significant amount of calibrated training images with known ground-truth illumination. Such calibration is time consuming, preferably done for each sensor individually, and therefore a major bottleneck in acquiring high color constancy accuracy. Statistics-based methods do not require calibrated training images, but they are less accurate. In this paper an unsupervised learning-based method is proposed that learns its parameter values after approximating the unknown ground-truth illumination of the training images, thus avoiding calibration. In terms of accuracy the proposed method outperforms all statistics-based and many learning-based methods. An extension of the method is also proposed, which learns the needed parameters from non-calibrated images taken with one sensor and which can then be successfully applied to images taken with another sensor. This effectively enables inter-camera unsupervised learning for color constancy. Additionally, a new high quality color constancy benchmark dataset with 1707 calibrated images is created, used for testing, and made publicly available. The results are presented and discussed. The source code and the dataset are available at http://www.fer.unizg.hr/ipg/resources/color_constancy/.


title: Tumbling motion of 1I/'Oumuamua reveals body's violent past

abstract: Models of the Solar System evolution show that almost all the primitive material leftover from the formation of the planets was ejected to the interstellar space as a result of dynamical instabilities. Accordingly, minor bodies should also be ejected from other planetary systems and should be abundant in the interstellar space. The number density of such objects, and prospects for their detection as they penetrate through the Solar System, were speculated about for decades, recently rising high hopes with the Pan-STARRS and LSST surveys. These expectations materialized on 18 October 2017 with the Pan-STARRS's discovery of 1I/'Oumuamua. Here we report homogeneous photometric observations of this body from Gemini North, which densely cover a total of 8 hr over two nights. A combined ultra-deep image of 1I/'Oumuamua shows no signs of cometary activity, confirming the results from earlier, less sensitive searches. Our data also show an enormous range of brightness variations > 2.5 mag, larger than ever observed in the population of Solar System objects, suggesting a very elongated shape of the body. But most significantly, the light curve does not repeat exactly from one rotation cycle to another and its double-peaked periodicity of 7.5483 $\pm$ 0.0073 hr from our data is inconsistent with earlier determinations. These are clear signs of a tumbling motion, a remarkable characteristic of 1I/'Oumuamua's rotation, consistent with a catastrophic collision in the distant past. This first example of an impacted minor body of exosolar origin indicates that collisional evolution of minor body populations in other planetary systems is not uncommon.


title: Double-layer force suppression between charged microspheres

abstract: In this paper we propose a protocol to suppress double-layer forces between two microspheres immersed in a dielectric medium, being one microsphere metallic at a controlled potential {\psi}M and the other a charged one either metallic or dielectric. The approach is valid for a wide range of distances between them. We show that, for a given distance between the two microspheres, the double-layer force can be totally suppressed by simply tuning {\psi}M up to values dictated by the linearized Poisson-Boltzmann equation. Our key finding is that such values can be substantially different from the ones predicted by the commonly used proximity force approximation (PFA), also known as Derjaguin approximation, even in situations where the latter is expected to be accurate. The proposed procedure can be used to suppress the double-layer interaction in force spectroscopy experiments, thus paving the way for measurements of other surface interactions, such as Casimir dispersion forces.


title: Capsize of polarization in dilute photonic crystals

abstract: We investigate, experimentally and theoretically, polarization rotation effects in dilute photonic crystals with transverse permittivity inhomogeneity perpendicular to the traveling direction of waves. A capsize, namely a drastic change of polarization to the perpendicular direction is observed in a one-dimensional photonic crystal in the frequency range $10\div 140$ GHz. To gain more insights into the rotational mechanism, we have developed a theoretical model of dilute photonic crystal, based on Maxwell's equations with a spatially dependent two dimensional inhomogeneous dielectric permittivity. We show that the polarization's rotation can be explained by an optical splitting parameter appearing naturally in Maxwell's equations for magnetic or electric fields components. This parameter is an optical analogous of Rashba like spin-orbit interaction parameter present in quantum waves, introduces a correction to the band structure of the two-dimensional Bloch states, creates the dynamical phase shift between the waves propagating in the orthogonal directions and finally leads to capsizing of the initial polarization. Excellent agreement between theory and experiment is found.


title: A compilation of LEGO Technic parts to support learning experiments on linkages

abstract: We present a compilation of LEGO Technic parts to provide easy-to-build constructions of basic planar linkages. Some technical issues and their possible solutions are discussed. To solve questions on fine details---like deciding whether the motion is an exactly straight line or not---we refer to the dynamic mathematics software tool GeoGebra.


title: Role of Skin Friction Drag during Flow-Induced Reconfiguration of a Flexible Thin Plate

abstract: We investigate drag reduction due to the flow-induced reconfiguration of a flexible thin plate in presence of skin friction drag at low Reynolds Number. The plate is subjected to a uniform free stream and is tethered at one end. We extend existing models in the literature to account for the skin friction drag. The total drag on the plate with respect to a rigid upright plate decreases due to flow-induced reconfiguration and further reconfiguration increases the total drag due to increase in skin friction drag. A critical value of Cauchy number ($Ca$) exists at which the total drag on the plate with respect to a rigid upright plate is minimum at a given Reynolds number. The reconfigured shape of the plate for this condition is unique, beyond which the total drag increases on the plate even with reconfiguration. The ratio of the form drag coefficient for an upright rigid plate and skin drag coefficient for a horizontal rigid plate ($\lambda$) determines the critical Cauchy number ($Ca_{cr}$). We propose modification in the drag scaling with free stream velocity ($F_{x}$ ${\propto}$ $U^{n}$) in presence of the skin friction drag. The following expressions of $n$ are found for $0.01 \leq Re \leq 1$, $n = 4/5 + {\lambda}/5$ for 1 $\leq$ $Ca$ $<$ $Ca_{cr}$ and $n = 1 + {\lambda}/5$ for $Ca_{cr} \leq Ca \leq 300$, where $Re$ is Reynolds number. We briefly discuss the combined effect of the skin friction drag and buoyancy on the drag reduction. An assessment of the feasibility of experiments is presented in order to translate the present model to physical systems.


title: Estimates for a DM $\rightarrow a \rightarrow \gamma$ 3.55 keV line in the radio lobes of Centaurus A

abstract: The origin of the 3.55 keV excess reported in various galaxies and galaxy clusters remains unknown. Non-observations of the line in dwarf spheroidal galaxies (dSphs) create tension with galaxy cluster observations when considering a direct Dark Matter (DM) $\rightarrow \gamma$ decay. This discrepancy could be explained by a DM $\rightarrow a \rightarrow \gamma$ decay, where the DM particles decay to axion-like particles (ALPs), which then convert to photons in a magnetic field. The 3.55 keV line from a decay to ALPs therefore has very different morphology from a direct decay to photons, as it would only appear in environments with substantial magnetic fields between the source and us. To test this hypothesis we consider targets where the line strength would be enhanced compared to the DM $\rightarrow \gamma$ model. The giant lobes of radio galaxies could represent good environments, having magnetic fields with strength $\mathcal{O}(\mu\rm{G})$ that extend for hundreds of kiloparsecs, and minimal X-ray background. In the case of Centaurus A, a large population of dSphs are associated with the galaxy, some of which are potentially in or behind its lobes. In this paper we estimate the flux of a 3.55 keV DM $\rightarrow a \rightarrow \gamma$ line from these sources, and compare it to the flux from the DM halo of Centaurus~A. We comment on the potential to detect such a signal with current satellites (such as XMM-Newton and Chandra) and future satellites (such as Athena).


title: Deep Neural Network Architectures for Modulation Classification

abstract: In this work, we investigate the value of employing deep learning for the task of wireless signal modulation recognition. Recently in [1], a framework has been introduced by generating a dataset using GNU radio that mimics the imperfections in a real wireless channel, and uses 10 different modulation types. Further, a convolutional neural network (CNN) architecture was developed and shown to deliver performance that exceeds that of expert-based approaches. Here, we follow the framework of [1] and find deep neural network architectures that deliver higher accuracy than the state of the art. We tested the architecture of [1] and found it to achieve an accuracy of approximately 75% of correctly recognizing the modulation type. We first tune the CNN architecture of [1] and find a design with four convolutional layers and two dense layers that gives an accuracy of approximately 83.8% at high SNR. We then develop architectures based on the recently introduced ideas of Residual Networks (ResNet [2]) and Densely Connected Networks (DenseNet [3]) to achieve high SNR accuracies of approximately 83.5% and 86.6%, respectively. Finally, we introduce a Convolutional Long Short-term Deep Neural Network (CLDNN [4]) to achieve an accuracy of approximately 88.5% at high SNR.


title: Large-scale structure phenomenology of viable Horndeski theories

abstract: Phenomenological functions $\Sigma$ and $\mu$, also known as $G_{\rm light}/G$ and $G_{\rm matter}/G$, are commonly used to parameterize modifications of the growth of large-scale structure in alternative theories of gravity. We study the values these functions can take in Horndeski theories, i.e. the class of scalar-tensor theories with second order equations of motion. We restrict our attention to models that are in a broad agreement with tests of gravity and the observed cosmic expansion history. In particular, we require the speed of gravity to be equal to the speed of light today, as required by the recent detection of gravitational waves and electromagnetic emission from a binary neutron star merger. We examine the correlations between the values of $\Sigma$ and $\mu$ analytically within the quasi-static approximation, and numerically, by sampling the space of allowed solutions. We confirm that the conjecture made in [Pogosian:2016pwr], that $(\Sigma-1)(\mu -1) \ge 0$ in viable Horndeski theories, holds very well. Along with that, we check the validity of the quasi-static approximation within different corners of Horndeski theory. Our results show that, even with the tight bound on the present day speed of gravitational waves, there is room within Horndeski theories for non-trivial signatures of modified gravity at the level of linear perturbations.


title: Relativistic Corrections to Nonrelativistic Effective Field Theories

abstract: In this paper we develop a formalism for studying the nonrelativistic limit of relativistic field theories in a systematic way. By introducing a simple, nonlocal field redefinition, we transform a given relativistic theory, describing a real, self-interacting scalar field, into an equivalent theory, describing a complex scalar field that encodes at each time both the original field and its conjugate momentum. Our low-energy effective theory incorporates relativistic corrections to the kinetic energy as well as the backreaction of fast-oscillating terms on the behavior of the dominant, slowly varying component of the field. Possible applications of our new approach include axion dark matter, though the methods developed here should be applicable to the low-energy limits of other field theories as well.


title: Bravyi-Kitaev Superfast simulation of fermions on a quantum computer

abstract: Present quantum computers often work with distinguishable qubits as their computational units. In order to simulate indistinguishable fermionic particles, it is first required to map the fermionic state to the state of the qubits. The Bravyi-Kitaev Superfast (BKSF) algorithm can be used to accomplish this mapping. The BKSF mapping has connections to quantum error correction and opens the door to new ways of understanding fermionic simulation in a topological context. Here, we present the first detailed exposition of BKSF algorithm for molecular simulation. We provide the BKSF transformed qubit operators and report on our implementation of the BKSF fermion-to-qubits transform in OpenFermion. In this initial study of the hydrogen molecule, we have compared BKSF, Jordan-Wigner and Bravyi-Kitaev transforms under the Trotter approximation. We considered different orderings of the exponentiated terms and found lower Trotter errors than previously reported for Jordan-Wigner and Bravyi-Kitaev algorithms. These results open the door to further study of the BKSF algorithm for quantum simulation.


title: Newton-Okounkov bodies, cluster duality, and mirror symmetry for Grassmannians

abstract: We use cluster structures and mirror symmetry to explicitly describe a natural class of Newton-Okounkov bodies for Grassmannians. We consider the Grassmannian $X=Gr_{n-k}(\mathbb C^n)$, as well as the mirror dual Landau-Ginzburg model $(\check{X}^\circ, W_q:\check{X}^\circ \to \mathbb C)$, where $\check{X}^\circ$ is the complement of a particular anti-canonical divisor in a Langlands dual Grassmannian $\check{X} = Gr_k((\mathbb C^n)^*)$, and the superpotential W_q has a simple expression in terms of Pl\"ucker coordinates. Grassmannians simultaneously have the structure of an $\mathcal{A}$-cluster variety and an $\mathcal{X}$-cluster variety. Given a cluster seed G, we consider two associated coordinate systems: a $\mathcal X$-cluster chart $\Phi_G:(\mathbb C^*)^{k(n-k)}\to X^{\circ}$ and a $\mathcal A$-cluster chart $\Phi_G^{\vee}:(\mathbb C^*)^{k(n-k)}\to \check{X}^\circ$. To each $\mathcal X$-cluster chart $\Phi_G$ and ample `boundary divisor' $D$ in $X\setminus X^{\circ}$, we associate a Newton-Okounkov body $\Delta_G(D)$ in $\mathbb R^{k(n-k)}$, which is defined as the convex hull of rational points. On the other hand using the $\mathcal A$-cluster chart $\Phi_G^{\vee}$ on the mirror side, we obtain a set of rational polytopes, described by inequalities, by writing the superpotential $W_q$ in the $\mathcal A$-cluster coordinates, and then "tropicalising". Our main result is that the Newton-Okounkov bodies $\Delta_G(D)$ and the polytopes obtained by tropicalisation coincide. As an application, we construct degenerations of the Grassmannian to toric varieties corresponding to these Newton-Okounkov bodies. Additionally, when $G$ corresponds to a plabic graph, we give a formula for the lattice points of the Newton-Okounkov bodies, which has an interpretation in terms of quantum Schubert calculus.


title: Adaptive finite element methods for sparse PDE-constrained optimization

abstract: We propose and analyze reliable and efficient a posteriori error estimators for an optimal control problem that involves a nondifferentiable cost functional, the Poisson problem as state equation and control constraints. To approximate the solution to the state and adjoint equations we consider a piecewise linear finite element method whereas three different strategies are used to approximate the control variable: piecewise constant discretization, piecewise linear discretization and the so-called variational discretization approach. For the first two aforementioned solution techniques we devise an error estimator that can be decomposed as the sum of four contributions: two contributions that account for the discretization of the control variable and the associated subgradient, and two contributions related to the discretization of the state and adjoint equations. The error estimator for the variational discretization approach is decomposed only in two contributions that are related to the discretization of the state and adjoint equations. On the basis of the devised a posteriori error estimators, we design simple adaptive strategies that yield optimal rates of convergence for the numerical examples that we perform.


title: Expected neutrino fluence from short Gamma-Ray Burst 170817A and off-axis angle constraints

abstract: We compute the expected neutrino fluence from SGRB 170817A, associated with the gravitational wave event GW 170817, directly based on Fermi observations in two scenarios: structured jet and off-axis (observed) top-hat jet. While the expected neutrino fluence for the structured jet case is very small, large off-axis angles imply high radiation densities in the jet, which can enhance the neutrino production efficiency. In the most optimistic allowed scenario, the neutrino fluence can reach only $10^{-4}$ of the sensitivity of the neutrino telescopes. We furthermore demonstrate that the fact that gamma-rays can escape limits the baryonic loading (energy in protons versus photons) and the off-axis angle for the internal shock scenario. In particular, for a baryonic loading of ten, the off-axis angle is more strongly constrained by the baryonic loading than by the time delay between the gravitational wave event and the onset of the gamma-ray emission.


title: Learning phase transitions from dynamics

abstract: We propose the use of recurrent neural networks for classifying phases of matter based on the dynamics of experimentally accessible observables. We demonstrate this approach by training recurrent networks on the magnetization traces of two distinct models of one-dimensional disordered and interacting spin chains. The obtained phase diagram for a well-studied model of the many-body localization transition shows excellent agreement with previously known results obtained from time-independent entanglement spectra. For a periodically-driven model featuring an inherently dynamical time-crystalline phase, the phase diagram that our network traces in a previously-unexplored regime coincides with an order parameter for its expected phases.


title: Implications from GW170817 and I-Love-Q relations for relativistic hybrid stars

abstract: Gravitational wave observations of GW170817 placed bounds on the tidal deformabilities of compact stars allowing one to probe equations of state for matter at supranuclear densities. Here we design new parametrizations for hybrid hadron-quark equations of state, that give rise to low-mass twin stars, and test them against GW170817. We find that GW170817 is consistent with the coalescence of a binary hybrid star--neutron star. We also test and find that the I-Love-Q relations for hybrid stars in the third family agree with those for purely hadronic and quark stars within $\sim 3\%$ for both slowly and rapidly rotating configurations, implying that these relations can be used to perform equation-of-state independent tests of general relativity and to break degeneracies in gravitational waveforms for hybrid stars in the third family as well.


title: Intermediate mass black holes and dark matter at the Galactic center

abstract: Could there be a large population of intermediate-mass black holes (IMBHs) formed in the early universe? Whether primordial or formed in Population III, these are likely to be very subdominant compared to the dark matter density, but could seed early dwarf galaxy/globular cluster and supermassive black hole formation. Via survival of dark matter density spikes, we show here that a centrally concentrated relic population of IMBHs, along with ambient dark matter, could account for the Fermi gamma-ray "excess" in the Galactic center because of dark matter particle annihilations.

categories: astro-ph.GA astro-ph.CO astro-ph.HE hep-ph


title: Braids and phase gates through high-frequency virtual tunneling of Majorana Zero Modes

abstract: Braiding of non-Abelian Majorana anyons is a first step towards using them in quantum computing. We propose a protocol for braiding Majorana zero modes formed at the edges of nanowires with strong spin orbit coupling and proximity induced superconductivity. Our protocol uses high frequency virtual tunneling between the ends of the nanowires in a tri-junction, which leads to an effective low frequency coarse grained dynamics for the system, to perform the braid. The braiding operation is immune to amplitude noise in the drives, and depends only on relative phase between the drives, which can be controlled by usual phase locking techniques. We also show how a phase gate, which is necessary for universal quantum computation, can be implemented with our protocol.


title: $SU(4)$-Ward identities for QCD with restored chiral symmetry

abstract: Lattice studies of QCD at temperatures above the chiral restoration and QCD with truncated low modes of the Dirac operator indicate approximate and explicit $SU(4)$ degeneracies in hadron spectra, respectively. At the same time, the QCD classical action and the path integral measure are not invariant under $SU(4)$. Here we investigate $SU(4)$ transformations in the continuum limit by deriving corresponding Ward identities. We show that, if there is a gap in low-lying modes of the Dirac operator the obtained $SU(4)$ Ward identities are simplified and look like they would be if $SU(4)$ symmetry is preserved. Then, we discuss possible consequences for the quark matter at high temperatures.