title: High precision determination of $\alpha_s$ from a global fit of jet rates

abstract: We present state-of-the-art extractions of the strong coupling based on N$^3$LO+NNLL accurate predictions for the two-jet rate in the Durham clustering algorithm at $e^+e^-$ collisions, as well as a simultaneous fit of the two- and three-jet rates taking into account correlations between the two observables. The fits are performed on a large range of data sets collected at LEP and PETRA colliders, with energies spanning from $35$ GeV to $207$ GeV. Owing to the high accuracy of the predictions used, the perturbative uncertainty is considerably smaller than that due to hadronization. Our best determination at the $Z$ mass is $\alpha_s(M_Z) = 0.11881 \pm 0.00063(exp.) \pm 0.00101(hadr.) \pm 0.00045(ren.) \pm 0.00034(res.)$, which is in agreement with the latest world average and has a comparable total uncertainty.


title: Sculpting out quantum correlations with bosonic subtraction

abstract: We present a method to extract M-partie bosonic correlations from an N-partite maximally symmetric state (M < N) with the help of successive applications of single-boson subtractions. We also propose an experimental photonic setup to implement it that can be done with the present technologies.


title: Topology of Learning in Artificial Neural Networks

abstract: Understanding how neural networks learn remains one of the central challenges in machine learning research. From random at the start of training, the weights of a neural network evolve in such a way as to be able to perform a variety of tasks, like classifying images. Here we study the emergence of structure in the weights by applying methods from topological data analysis. We train simple feedforward neural networks on the MNIST dataset and monitor the evolution of the weights. When initialized to zero, the weights follow trajectories that branch off recurrently, thus generating trees that describe the growth of the effective capacity of each layer. When initialized to tiny random values, the weights evolve smoothly along two-dimensional surfaces. We show that natural coordinates on these learning surfaces correspond to important factors of variation.


title: Galaxy shape measurement with convolutional neural networks

abstract: We present our results from training and evaluating a convolutional neural network (CNN) to predict galaxy shapes from wide-field survey images of the first data release of the Dark Energy Survey (DES DR1). We use conventional shape measurements as ground truth from an overlapping, deeper survey with less sky coverage, the Canada-France Hawaii Telescope Lensing Survey (CFHTLenS). We demonstrate that CNN predictions from single band DES images reproduce the results of CFHTLenS at bright magnitudes and show higher correlation with CFHTLenS at fainter magnitudes than maximum likelihood model fitting estimates in the DES Y1 im3shape catalogue. Prediction of shape parameters with a CNN is also extremely fast, it takes only 0.2 milliseconds per galaxy, improving more than 4 orders of magnitudes over forward model fitting. The CNN can also accurately predict shapes when using multiple images of the same galaxy, even in different color bands, with no additional computational overhead. The CNN is again more precise for faint objects, and the advantage of the CNN is more pronounced for blue galaxies than red ones when compared to the DES Y1 metacalibration catalogue, which fits a single Gaussian profile using riz band images. We demonstrate that CNN shape predictions within the metacalibration self-calibrating framework yield shear estimates with negligible multiplicative bias, $ m < 10^{-3}$, and no significant PSF leakage. Our proposed setup is applicable to current and next generation weak lensing surveys where higher quality ground truth shapes can be measured in dedicated deep fields.


title: Asymptotics of Hankel determinants with a Laguerre-type or Jacobi-type potential and Fisher-Hartwig singularities

abstract: We obtain large $n$ asymptotics of $n \times n$ Hankel determinants whose weight has a one-cut regular potential and Fisher-Hartwig singularities. We restrict our attention to the case where the associated equilibrium measure possesses either one soft edge and one hard edge (Laguerre-type) or two hard edges (Jacobi-type). We also present some applications in the theory of random matrices. In particular, we can deduce from our results asymptotics for partition functions with singularities, central limit theorems, correlations of the characteristic polynomials, and gap probabilities for (piecewise constant) thinned Laguerre and Jacobi-type ensembles. Finally, we mention some links with the topics of rigidity and Gaussian multiplicative chaos.


title: Optimal Power Flow for DC Networks with Robust Feasibility and Stability Guarantees

abstract: With high penetrations of renewable generation and variable loads, there is significant uncertainty associated with power flows in DC networks such that stability and operational constraint satisfaction are of concern. Most existing DC network optimal power flow (DN-OPF) formulations assume exact knowledge of loading conditions and do not provide stability guarantees. In contrast, this paper studies a DN-OPF formulation which considers both stability and operational constraint satisfaction under uncertainty. The need to account for a range of uncertainty realizations in this paper's robust optimization formulation results in a challenging semi-infinite program (SIP). The proposed solution algorithm reformulates this SIP into a computationally tractable problem by constructing a tight convex inner approximation of the stability set using sufficient conditions for the existence of a feasible and stable power flow solution. Optimal generator setpoints are obtained by optimizing over the proposed convex stability set. The validity and value of the proposed algorithm are demonstrated through various DC networks adapted from IEEE test cases.


title: Learning-Based Proxy Collision Detection for Robot Motion Planning Applications

abstract: This paper demonstrates that collision detection-intensive applications such as robotic motion planning may be accelerated by performing collision checks with a machine learning model. We propose Fastron, a learning-based algorithm to model a robot's configuration space to be used as a proxy collision detector in place of standard geometric collision checkers. We demonstrate that leveraging the proxy collision detector results in up to an order of magnitude faster performance in robot simulation and planning than state-of-the-art collision detection libraries. Our results show that Fastron learns a model more than 100 times faster than a competing C-space modeling approach, while also providing theoretical guarantees of learning convergence. Using the OMPL motion planning libraries, we were able to generate initial motion plans across all experiments with varying robot and environment complexities. With Fastron, we can repeatedly perform planning from scratch at a 56 Hz rate, showing its application toward autonomous surgical assistance task in shared environments with human-controlled manipulators. All performance gains were achieved despite using only CPU-based calculations, suggesting further computational gains with a GPU approach that can parallelize tensor algebra. Code is available online.


title: The harmonicity of slice regular functions

abstract: In this article we investigate harmonicity, Laplacians, mean value theorems and related topics in the context of quaternionic analysis. We observe that a Mean Value Formula for slice regular functions holds true and it is a consequence of the well known Representation Formula for slice regular functions over $\mathbb{H}$. Motivated by this observation, we have constructed three order-two differential operators in the kernel of which slice regular functions are, answering positively to the question: is a slice regular function over $\mathbb{H}$ (analogous to an holomorphic function over $\mathbb{C}$) "harmonic" in some sense, i.e. is it in the kernel of some order-two differential operator over $\mathbb{H}$ ? Finally, some applications are deduced, such as a Poisson Formula for slice regular functions over $\mathbb{H}$ and a Jensen's Formula for semi-regular ones.


title: Marginal deformations and quasi-Hopf algebras

abstract: We establish the existence of a quasi-Hopf algebraic structure underlying the Leigh-Strassler N=1 superconformal marginal deformations of the N=4 Super-Yang-Mills theory. The scalar-sector R-matrix of these theories, which is related to their one-loop spin chain Hamiltonian, does not generically satisfy the Quantum Yang-Baxter Equation. By constructing a Drinfeld twist which relates this R-matrix to that of the N=4 SYM theory, but also produces a non-trivial co-associator, we show that the generic Leigh-Strassler R-matrix satisfies the quasi-Hopf version of the QYBE. We also use the twist to define a suitable star product which directly relates the N=4 SYM superpotential to that of the marginally deformed gauge theories. We expect our results to be relevant to studies of integrability (and its breaking) in these theories, as well as to provide useful input for supergravity solution-generating techniques.


title: A Stacked Autoencoder Application for Residential Load Curve Forecast and Peak Shaving

abstract: For the last ten years, utilities have observed on-going transitions on consumers' load curves. The previously flat load curves have more frequently turned into duck-shape. This is jointly caused by the increasing household loads as well as the popularity of rooftop solar photovoltaic. Such load curve transition challenges the operational flexibility of the existing systems and greatly increases the per-MWh energy costs. Peak shaving, in this context, becomes a critical task to demand-side management. Owing to the development of Battery Energy Storage Systems (BESS), numerous peak shaving strategies have been developed and implemented. In this paper, by applying a stacked autoencoder (SAE)-based residential peak load curve forecasting technology, we further lift the peaking shaving capabilities of BESSs to a new level. The proposed strategy takes into account the welfares of both generation-side and demand-side and reaches an optimal balance. A comprehensive case study using smart meter data demonstrates the effectiveness of the proposed method in peak shaving application.


title: On the anticipative nonlinear filtering problem and its stability

abstract: In this paper, we consider an anticipative nonlinear filtering problem, in which the observation noise is correlated with the past of the signal. This new signal-observation model has its applications in both finance models with insider trading and in engineering. We derive a new equation for the filter in this context, analyzing both the nonlinear and the linear cases. We also handle the case of a finite filter with Volterra type observation. The performance of our algorithm is presented through numerical experiments.


title: On modules $M$ with $\tau(M) \cong \nu \Omega^{d+2}(M)$ for isolated singularities of Krull dimension $d$

abstract: A classical formula for the Auslander-Reiten translate $\tau$ says that $\tau(M)\cong \nu \Omega^2(M)$ for every indecomposable module $M$ of a selfinjective Artin algebra. We generalise this by showing that for a $2d$-periodic isolated singularity $A$ of Krull dimension $d$, we have for the Auslander-Reiten translate of an indecomposable non-projective Cohen-Macaulay $A$-module $M$, $\tau(M)\cong \nu \Omega^{d+2}(M)$ if and only if $Ext_A^{d+1}(M,A)=Ext_A^{d+2}(M,A)=0$. We give several applications for Artin algebras.


title: On the choice of a conformal Gauss gauge near the cylinder representing spatial infinity

abstract: A convenient approach to analyze spatial infinity is to use a cylinder representation $I$ and impose a gauge based on a congruence of conformal geodesics. This so-called conformal Gauss gauge comes along with the freedom to specify initial data for the conformal geodesics. Such a gauge has been constructed from an ordinary Cauchy surface and from past null infinity $\mathcal{J}^-$, respectively. The purpose of this note is to compare these gauges near the critical set $I^-$, where $I$ "touches" $\mathcal{J}^-$, as it turns out that they are related in a somewhat unexpected intricate way.


title: A Dictionary Based Generalization of Robust PCA

abstract: We analyze the decomposition of a data matrix, assumed to be a superposition of a low-rank component and a component which is sparse in a known dictionary, using a convex demixing method. We provide a unified analysis, encompassing both undercomplete and overcomplete dictionary cases, and show that the constituent components can be successfully recovered under some relatively mild assumptions up to a certain $\textit{global}$ sparsity level. Further, we corroborate our theoretical results by presenting empirical evaluations in terms of phase transitions in rank and sparsity for various dictionary sizes.


title: Scattering amplitudes versus potentials in nuclear effective field theory: is there a potential compromise?

abstract: In effective field theory physical quantities, in particular observables, are expressed as a power series in terms of a small expansion parameter. For non-perturbative systems, for instance nuclear physics, this requires the non-perturbative treatment of at least a part of the interaction (or the potential, if we are dealing with a non-relativistic system). This is not entirely trivial and as a consequence different interpretations on how to treat these systems have appeared. A practical approach is to expand the effective potential, where this potential is later fully iterated in the Schr\"odinger equation for obtaining amplitudes and observables. The expectation is that this will lead to observables that will have an implicit power counting expansion. Here we explicitly check whether the amplitudes are actually following the same power counting as the potential. It happens that reality does not necessarily conform to expectations and the amplitudes will often violate the power counting with which the potential has been expanded. A more formal approach is to formulate the expansion directly in terms of amplitudes and observables, which is the original aim of the effective field theory idea. Yet this second approach is technically complicated. We explore here the possibility of constructing potentials that when fully iterated will make sure that amplitudes indeed are expansible in terms of a small expansion parameter.


title: Eventually entanglement breaking Markovian dynamics: structure and characteristic times

abstract: We investigate entanglement breaking times of Markovian evolutions in discrete and continuous time. In continuous time, we characterize which Markovian evolutions are eventually entanglement breaking, that is, evolutions for which there is a finite time after which any entanglement initially present has been destroyed by the noisy evolution. In the discrete time framework, we consider the entanglement breaking index, that is, the number of times a quantum channel has to be composed with itself before it becomes entanglement breaking. The PPT-square conjecture is that every PPT quantum channel has an entanglement breaking index of at most 2; we prove that every faithful PPT quantum channel has a finite entanglement breaking index, and more generally, any faithful PPT CP map whose Hilbert-Schmidt adjoint is also faithful is eventually entanglement breaking. We also provide a method to obtain concrete bounds on this index for any faithful quantum channel. To obtain these estimates, we use a notion of robustness of separability to obtain bounds on the radius of the largest separable ball around faithful product states. We also extend the framework of Poincar\'e inequalities for nonprimitive semigroups to the discrete setting to quantify the convergence of quantum semigroups in discrete time, which is of independent interest.


title: Stochastic thermodynamics of self-oscillations: the electron shuttle

abstract: Self-oscillation is a phenomenon studied across many scientific disciplines, including the engineering of efficient heat engines and electric generators. We investigate the single electron shuttle, a model nano-scale system that exhibits a spontaneous transition towards self-oscillation, from a thermodynamic perspective. We analyze the model at three different levels of description: The fully stochastic level based on Fokker-Planck and Langevin equations, the mean-field level, and a perturbative solution to the Fokker-Planck equation that works particularly well for small oscillation amplitudes. We provide consistent derivations of the laws of thermodynamics for this model system at each of these levels. At the mean-field level, an abrupt transition to self-oscillation arises from a Hopf bifurcation of the deterministic equations of motion. At the stochastic level, this transition is smeared out by noise, but vestiges of the bifurcation remain visible in the stationary probability density. At all levels of description, the transition towards self-oscillation is reflected in thermodynamic quantities such as heat flow, work and entropy production rate. Our analysis provides a comprehensive picture of a nano-scale self-oscillating system, with stochastic and deterministic models linked by a unifying thermodynamic perspective.


title: Mathematical analysis and numerical solution of models with dynamic Preisach hysteresis

abstract: Hysteresis is a phenomenon that is observed in a great variety of physical systems, which leads to a nonlinear and multivalued behavior, making their modeling and control difficult. Even though the analysis and mathematical properties of classical or rate-independent hysteresis models are known, this is not the case for dynamic models where current approaches lack a proper functional analytic framework which is essential to formulate optimization problems and develop stable numerics, both being crucial in practice. This paper deals with the description and mathematical analysis of the dynamic Preisach hysteresis model. Toward that end, we complete a widely accepted definition of the dynamic model commonly used to describe the constitutive relation between the magnetic field H and the magnetic induction B, in which, the values of B not only depends on the present values of H but also on the past history and its velocity. We first analyze mathematically some important properties of the model and compare them with known results for the static Preisach model. Then, we consider a parabolic problem with dynamic hysteresis motivated by electromagnetic field equations. Under suitable assumptions, we show the well posedness of a weak formulation of the problem and solve it numerically. Finally, we report a numerical test in order to assess the order of convergence and to illustrate the behavior of the numerical solution for different configurations of the dynamic Preisach model


title: Ricci $\rho$-solitons on 3-dimensional $\eta$-Einstein almost Kenmotsu mandifolds

abstract: In this paper the notion of Ricci $\rho$-soliton as a generalization of Ricci soliton is defined. We are motivated by the Ricci-Bourguignon flow to define this concept. We show that if a 3- dimensional almost Kenmotsu Einstein manifold $M$ be a $\rho$-soliton, then $M$ is a Kenmotsu manifold of constant sectional curvature $-1$ and the $\rho$-soliton is expanding, with $\lambda =2$.


title: On the growth rate of chromatic numbers of finite subgraphs

abstract: We prove that, for every function $f:\mathbb{N} \rightarrow \mathbb{N}$, there is a graph $G$ with uncountable chromatic number such that, for every $k \in \mathbb{N}$ with $k \geq 3$, every subgraph of $G$ with fewer than $f(k)$ vertices has chromatic number less than $k$. This answers a question of Erd\H{o}s, Hajnal, and Szemeredi.


title: Variational Operators, Symplectic Operators, and the Cohomology of Scalar Evolution Equations

abstract: For a scalar evolution equation $u_t=K(t,x,u,u_x,\ldots, u_n), n\geq 2$ the cohomology spaces $H^{1,s}({\mathcal R}^\infty)$ vanishes for $s\geq 3$ while the space $H^{1,2}({\mathcal R}^\infty)$ is isomorphic to the space of variational operators. The cohomology space $H^{1,2}({\mathcal R}^\infty)$ is also shown to be isomorphic to the space of symplectic operators for $u_t=K$ for which the equation is Hamiltonian. Third order scalar evolution equations admitting a first order symplectic (or variational) operator are characterized. The symplectic nature of the potential form of a bi-Hamiltonian evolution equation is also presented.


title: Online Sampling from Log-Concave Distributions

abstract: Given a sequence of convex functions $f_0, f_1, \ldots, f_T$, we study the problem of sampling from the Gibbs distribution $\pi_t \propto e^{-\sum_{k=0}^tf_k}$ for each epoch $t$ in an online manner. Interest in this problem derives from applications in machine learning, Bayesian statistics, and optimization where, rather than obtaining all the observations at once, one constantly acquires new data, and must continuously update the distribution. Our main result is an algorithm that generates roughly independent samples from $\pi_t$ for every epoch $t$ and, under mild assumptions, makes $\mathrm{polylog}(T)$ gradient evaluations per epoch. All previous results imply a bound on the number of gradient or function evaluations which is at least linear in $T$. Motivated by real-world applications, we assume that functions are smooth, their associated distributions have a bounded second moment, and their minimizer drifts in a bounded manner, but do not assume they are strongly convex. In particular, our assumptions hold for online Bayesian logistic regression, when the data satisfy natural regularity properties, giving a sampling algorithm with updates that are poly-logarithmic in $T$. In simulations, our algorithm achieves accuracy comparable to an algorithm specialized to logistic regression. Key to our algorithm is a novel stochastic gradient Langevin dynamics Markov chain with a carefully designed variance reduction step and constant batch size. Technically, lack of strong convexity is a significant barrier to analysis and, here, our main contribution is a martingale exit time argument that shows our Markov chain remains in a ball of radius roughly poly-logarithmic in $T$ for enough time to reach within $\varepsilon$ of $\pi_t$.


title: Instability of exotic compact objects and its implications for gravitational-wave echoes

abstract: Exotic compact objects (ECOs) have recently become an exciting research subject, since they are speculated to have a special response to the incident gravitational waves (GWs) that leads to GW echoes. We show that energy carried by GWs can easily cause the event horizon to form out of a static ECO --- leaving no echo signals towards spatial infinity. To show this, we use the ingoing Vaidya spacetime and take into account the back reaction due to incoming GWs. Demanding that an ECO does not collapse into a black hole puts an upper bound on the compactness of the ECO, at the cost of less distinct echo signals for smaller compactness. The trade-off between echoes' detectability and distinguishability leads to a fine tuning of ECO parameters for LIGO to find distinct echoes. We also show that an extremely compact ECO that can survive the gravitational collapse and give rise to GW echoes might have to expand its surface in a non-causal way.


title: Density distributions and depth in flocks

abstract: Recent experimental evidence suggests that interactions in flocks of birds do not involve a characteristic length scale. Bird flocks have also been revealed to have an inhomogeneous density distribution, with the density of birds near the border greater than near the centre. We introduce a strictly metric-free model for collective behaviour that incorporates a distributed motional bias, providing control of the density distribution. A simple version of this model is then able to provide a good fit to published data for the density variation across flocks of Starlings. We find that it is necessary for individuals on the edge of the flock to have an inward motional bias but that birds in the interior of the flock instead must have an outward bias. We discuss the ability of individuals to determine their depth within a flock and show how this might be achieved by relatively simple analysis of their visual environment.


title: Finding the Needle in a Haystack: Detrending Photometric Timeseries Data of Strictly Periodic Astrophysical Objects

abstract: Light curves of astrophysical objects frequently contain strictly periodic signals. In those cases we can use that property to aid the detrending algorithm to fully disentangle an unknown periodic signal and an unknown baseline signal with no power at that period. The periodic signal is modeled as a discrete probability distribution function (pdf), while the baseline signal is modeled as a residual timeseries. Those two components are disentangled by minimizing the length of the residual timeseries w.r.t. the per-bin pdf fluxes. We demonstrate the use of the algorithm on a synthetic case, on the eclipsing binary KIC 3953981 and on the eccentric ellipsoidal variable KIC 3547874. We further discuss the parameters and the limitations of the algorithm and speculate on the two most common use cases: detrending the periodic signal of interest and measuring the dependence of instrumental response on controlled instrumental variables. A more sophisticated version of the algorithm is released as open source on github and available via pip.


title: Policies for allocation of information in task-oriented groups: elitism and egalitarianism outperform welfarism

abstract: Communication or influence networks are probably the most controllable of all factors that are known to impact on the problem-solving capability of task-forces. In the case connections are costly, it is necessary to implement a policy to allocate them to the individuals. Here we use an agent-based model to study how distinct allocation policies affect the performance of a group of agents whose task is to find the global maxima of NK fitness landscapes. Agents cooperate by broadcasting messages informing on their fitness and use this information to imitate the fittest agent in their influence neighborhoods. The larger the influence neighborhood of an agent, the more links, and hence information, the agent receives. We find that the elitist policy in which agents with above-average fitness have their influence neighborhoods amplified, whereas agents with below-average fitness have theirs deflated, is optimal for smooth landscapes, provided the group size is not too small. For rugged landscapes, however, the elitist policy can perform very poorly for certain group sizes. In addition, we find that the egalitarian policy, in which the size of the influence neighborhood is the same for all agents, is optimal for both smooth and rugged landscapes in the case of small groups. The welfarist policy, in which the actions of the elitist policy are reversed, is always suboptimal, i.e., depending on the group size it is outperformed by either the elitist or the egalitarian policies.


title: Stability conditions in families

abstract: We develop a theory of Bridgeland stability conditions and moduli spaces of semistable objects for a family of varieties. Our approach is based on and generalizes previous work by Abramovich-Polishchuk, Kuznetsov, Lieblich, and Piyaratne-Toda. Our notion includes openness of stability, semistable reduction, a support property uniformly across the family, and boundedness of semistable objects. We show that such a structure exists whenever stability conditions are known to exist on the fibers. Our main application is the generalization of Mukai's theory for moduli spaces of semistable sheaves on K3 surfaces to moduli spaces of Bridgeland semistable objects in the Kuznetsov component associated to a cubic fourfold. This leads to the extension of theorems by Addington-Thomas and Huybrechts on the derived category of special cubic fourfolds, to a new proof of the integral Hodge conjecture, and to the construction of an infinite series of unirational locally complete families of polarized hyperk\"ahler manifolds of K3 type. Other applications include the deformation-invariance of Donaldson-Thomas invariants counting Bridgeland stable objects on Calabi-Yau threefolds, and a method for constructing stability conditions on threefolds via degeneration.


title: How Analytic Choices Can Affect the Extraction of Electromagnetic Form Factors from Elastic Electron Scattering Cross Section Data

abstract: Scientists often try to incorporate prior knowledge into their regression algorithms, such as a particular analytic behavior or a known value at a kinematic endpoint. Unfortunately, there is often no unique way to make use of this prior knowledge, and thus, different analytic choices can lead to very different regression results from the same set of data. To illustrate this point in the context of the proton electromagnetic form factors, we use the Mainz elastic data with its 1422 cross section points and 31 normalization parameters. Starting with a complex unbound non-linear regression, we will show how the addition of a single theory-motivated constraint removes an oscillation from the magnetic form factor and shifts the extracted proton charge radius. We then repeat both regressions using the same algorithm, but with a rebinned version of the Mainz dataset. These examples illustrate how analytic choices, such as the function that is being used or even the binning of the data, can dramatically affect the results of a complex regression. These results also demonstrate why it is critical when using regression algorithms to have either a physical model in mind or a firm mathematical basis


title: Magneto-optical Signatures of Volkov-Pankratov States in Topological Insulators

abstract: In addition to the usual chiral surface states, massive surface states can arise at a smooth interface between a topological and a trivial bulk insulator. While not subject to topological protection as the chiral states, these massive states, theorized by Volkov and Pankratov in the 1980s, reflect nevertheless emergent Dirac physics at the interface. We study theoretically the magneto-optical response of these surface states, which is strikingly different from that of the bulk states. Most saliently, we show that these states can be identified clearly in the presence of a magnetic field and its orientation with respect to the interface.


title: When Primordial Black Holes from Sound Speed Resonance Meet a Stochastic Background of Gravitational Waves

abstract: As potential candidates of dark matter, primordial black holes (PBHs) are within the core scopes of various astronomical observations. In light of the explosive development of gravitational wave (GW) and radio astronomy, we thoroughly analyze a stochastic background of cosmological GWs, induced by over large primordial density perturbations, with several spikes that was inspired by the sound speed resonance effect and can predict a particular pattern on the mass spectrum of PBHs. With a specific mechanicsm for PBHs formation, we for the first time perform the study of such induced GWs that originate from both the inflationary era and the radiation-dominated phase. We report that, besides the traditional process of generating GWs during the radiation-dominated phase, the contribution of the induced GWs in the sub-Hubble regime during inflation can become significant at critical frequency band because of a narrow resonance effect. All contributions sum together to yield a specific profile of the energy spectrum of GWs that can be of observable interest in forthcoming astronomical experiments. Our study shed light on the possible joint probe of PBHs via various observational windows of multi-messenger astronomy, including the search for electromagnetic effects with astronomical telescopes and the stochastic background of relic GWs with GW instruments.


title: Investigating dynamical properties of evolved Galactic open clusters

abstract: The stellar content of Galactic open clusters is gradually depleted during their evolution as a result of internal relaxation and external interactions. The final residues of the evolution of open clusters are called open cluster remnants, barely distinguishable from the field. We aimed to characterise and compare the dynamical states of a set of 16 such objects. The sample also includes 7 objects that are catalogued as dynamically evolved open clusters. We used photometric data from the 2MASS, astrometric data from the GAIA DR2 and a decontamination algorithm that was applied to the three-dimensional astrometric space of proper motions and parallaxes for stars in the objects' areas. The luminosity and mass functions and total masses for most open cluster remnants are derived here for the first time. Our analysis used predictions of N-body simulations to estimate the initial number of stars of the remnants from their dissolution timescales. The investigated open cluster remnants present masses and velocity dispersions within well-defined ranges: M between ~10-40M_Sun and sigma_v between ~1-7km/s. Some objects in the remnant sample have a limiting radius R_lim<~2pc, which means that they are more compact than the investigated open clusters; other remnants have R_lim between ~2-7pc, which is comparable to the open clusters. We suggest that the open cluster NGC2180 is entering a remnant evolutionary stage. In general, our clusters show signals of depletion of low-mass stars. This confirms their dynamically evolved states. We conclude that the open cluster remnants we studied are in fact remnants of initially very populous open clusters (No~10^3-10^4 stars). The outcome of the long-term evolution is to bring the final residues of the open clusters to dynamical states that are similar to each other, thus masking out the memory of the initial formation conditions of star clusters.


title: Moduli Spaces of Germs of Semiquasihomogeneous Legendrian Curves

abstract: We construct a moduli space for Legendrian curves singularities which are contactomorphic-equivalent and equisingular through a contact analogue of the Kodaira-Spencer map for curve singularities. We focus on the specific case of Legendrian curves which are the conormal of a plane curve with one Puiseux pair.


title: A count-based imaging model for the Spectrometer/Telescope for Imaging X-rays (STIX) in Solar Orbiter

abstract: The Spectrometer/Telescope for Imaging X-rays (STIX) will look at solar flares across the hard X-ray window provided by the Solar Orbiter cluster. Similarly to the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI), STIX is a visibility-based imaging instrument, which will ask for Fourier-based image reconstruction methods. However, in this paper we show that, as for RHESSI, also for STIX count-based imaging is possible. Specifically, here we introduce and illustrate a mathematical model that mimics the STIX data formation process as a projection from the incoming photon flux into a vector made of 120 count components. Then we test the reliability of Expectation Maximization for image reconstruction in the case of several simulated configurations typical of flare morphology.


title: FLAG Review 2019

abstract: We review lattice results related to pion, kaon, $D$-meson, $B$-meson, and nucleon physics with the aim of making them easily accessible to the nuclear and particle physics communities. More specifically, we report on the determination of the light-quark masses, the form factor $f_+(0)$ arising in the semileptonic $K \to \pi$ transition at zero momentum transfer, as well as the decay constant ratio $f_K/f_\pi$ and its consequences for the CKM matrix elements $V_{us}$ and $V_{ud}$. Furthermore, we describe the results obtained on the lattice for some of the low-energy constants of $SU(2)_L\times SU(2)_R$ and $SU(3)_L\times SU(3)_R$ Chiral Perturbation Theory. We review the determination of the $B_K$ parameter of neutral kaon mixing as well as the additional four $B$ parameters that arise in theories of physics beyond the Standard Model. For the heavy-quark sector, we provide results for $m_c$ and $m_b$ as well as those for $D$- and $B$-meson decay constants, form factors, and mixing parameters. These are the heavy-quark quantities most relevant for the determination of CKM matrix elements and the global CKM unitarity-triangle fit. We review the status of lattice determinations of the strong coupling constant $\alpha_s$. Finally, in this review we have added a new section reviewing results for nucleon matrix elements of the axial, scalar and tensor bilinears, both isovector and flavor diagonal.


title: Jointly Sparse Convolutional Neural Networks in Dual Spatial-Winograd Domains

abstract: We consider the optimization of deep convolutional neural networks (CNNs) such that they provide good performance while having reduced complexity if deployed on either conventional systems with spatial-domain convolution or lower-complexity systems designed for Winograd convolution. The proposed framework produces one compressed model whose convolutional filters can be made sparse either in the spatial domain or in the Winograd domain. Hence, the compressed model can be deployed universally on any platform, without need for re-training on the deployed platform. To get a better compression ratio, the sparse model is compressed in the spatial domain that has a fewer number of parameters. From our experiments, we obtain $24.2\times$ and $47.7\times$ compressed models for ResNet-18 and AlexNet trained on the ImageNet dataset, while their computational cost is also reduced by $4.5\times$ and $5.1\times$, respectively.


title: Understanding quality control of hard metals in industry -- A quantum mechanics approach

abstract: For many decades, the magnetic saturation of, e.g. hard metals (HM) such as WC-Co-based cemented carbides, has been used as process and quality control in industry to ensure consistency of product properties. In an urge of replacing cobalt as a binder phase, a demand on understanding the magnetic response as a function of composition on the atomic scale is growing. In this paper, a theoretical description of the measured weight specific magnetic saturation of hard metals as a function of the tungsten weight fraction present in the cobalt binder phase, based on first-principle calculations, has been established for standard WC-Co. The predicted magnetic saturation agrees well with the experimental one. Furthermore, it is proposed that the theoretical description can be extended to alternative and more complex binder phases which allows to transfer the production control to those hard metals.


title: Max-plus Linear Inverse Problems: 2-norm regression and system identification of max-plus linear dynamical systems with Gaussian noise

abstract: In this paper we present new theory and algorithms for 2-norm regression over the max-plus semiring. As an application we also show how max-plus 2-norm regression can be used in system identification of max-plus linear dynamical systems with Gaussian noise. We also introduce and provide methods for solving a max-plus linear inverse problem with regularization, which can be used when the the original problem is not well posed.


title: Neutrino eigenstates and flavour, spin and spin-flavour oscillations in a constant magnetic field

abstract: We further develop a recently proposed new approach to the description of the relativistic neutrino flavour $\nu_e^L \leftrightarrow \nu_{\mu}^L$, spin $\nu_e^L \leftrightarrow \nu_{e}^R$ and spin-flavour $\nu_e^L \leftrightarrow \nu_{\mu}^R$ oscillations in a constant magnetic field that is based on the use of the exact neutrino stationary states in the magnetic field. The neutrino flavour, spin and spin-flavour oscillations probabilities are calculated accounting for the whole set of possible conversions between four neutrino states. In general, the obtained expressions for the neutrino oscillations probabilities exhibit new inherent features in the oscillation patterns. It is shown, in particular, that: 1) in the presence of the transversal magnetic field for a given choice of parameters (the energy and magnetic moments of neutrinos and the strength of the magnetic field) the amplitude of the flavour oscillations $\nu_e^L \leftrightarrow \nu_{\mu}^L$ at the vacuum frequency is modulated by the magnetic field frequency, 2) the neutrino spin oscillation probability (without change of the neutrino flavour) exhibits the dependence on the mass square difference $\Delta m^2$. It is shown that the discussed interplay of neutrino oscillations in magnetic fields on different frequencies can have important consequences in astrophysical environments, in particular in those peculiar for magnetars.


title: Optimization with delay-induced bifurcations

abstract: Optimization is finding the best solution, which mathematically amounts to locating the global minimum of some cost function. Optimization is traditionally automated with digital or quantum computers, each having their limitations and none guaranteeing an optimal solution. Here, we conceive a principle behind optimization based on delay-induced bifurcations, which is potentially implementable in non-quantum analog devices. Often, optimization techniques are interpreted via a particle moving in multi-well energy landscape and to prevent confinement to a non-global minima they should incorporate mechanisms to overcome barriers between the minima. Particularly, simulated annealing digitally emulates pushing a fictitious particle over a barrier by random noise, whereas quantum computers utilize tunneling through barriers. In our principle, the barriers are effectively destroyed by delay-induced bifurcations. Although bifurcation scenarios in nonlinear delay-differential equations can be very complex and are notoriously difficult to predict, we hypothesize, verify, and utilize the finding that they become considerably more predictable in dynamical systems, where the right-hand side depends only on the delayed variable and represents a gradient of some potential energy function. By tuning the delay introduced into the gradient descent setting, thanks to global bifurcations destroying local attractors, one could force the system to spontaneously wander around all minima. This would be similar to noise-induced behavior in simulated annealing but achieved deterministically. Ideally, a slow increase and then decrease of the delay should automatically push the system toward the global minimum. We explore the possibility of this scenario and formulate some prerequisites.


title: More on the Structure of Extreme Level Sets in Branching Brownian Motion

abstract: This work is a continuation of the manuscript "the structure of extreme level sets in branching Brownian motion", in which the same authors studied the fine structure of the extreme level sets of branching Brownian motion, namely the sets of particles whose height is within a finite distance from the global maximum. It is well known that such particles congregate at large times in clusters of order-one genealogical diameter around local maxima which form a Cox process in the limit. Our main finding here is that most of the particles in an extreme level set come from only a small fraction of the clusters, which are atypically large.


title: Kinetic modeling of alcohol consumption

abstract: In most countries, alcohol consumption distributions have been shown to possess universal features. Their unimodal right-skewed shape is usually modeled in terms of the Lognormal distribution, which is easy to fit, test, and modify. However, empirical distributions often deviate considerably from the Lognormal model, and both Gamma and Weibull distributions appear to better describe the survey data. In this paper we explain the appearance of these distributions by means of classical methods of kinetic theory of multi-agent systems. The microscopic variation of alcohol consumption of agents around a universal \emph{social} accepted value of consumption, is built up introducing as main criterion for consumption a suitable value function in the spirit of the prospect theory of Kahneman and Twersky. The mathematical properties of the value function then determine the unique macroscopic equilibrium which results to be a generalized Gamma distribution. The modeling of the microscopic kinetic interaction allows to clarify the meaning of the various parameters characterizing the generalized Gamma equilibrium.


title: In Vivo Wireless Channel Modeling

abstract: In this chapter, the state of the art of in vivo wireless channel characterization has been presented. Various studies described in the literature are dedicated to the in vivo communication channel, and they consider different parameters in studying various anatomical regions. Furthermore, the location-dependent characteristics of in vivo wireless communication at 915 MHz are analyzed in detail via numerical and experimental investigations. A complete model for the in vivo channel is not available and remains an open research problem. However, considering the expected future growth of implanted technologies and their potential use for the detection and diagnosis of various health-related issues in the human body, the channel modeling studies should be further extended to develop better and more efficient communications systems for future in vivo systems.


title: Conservation of asymptotic charges from past to future null infinity: Supermomentum in general relativity

abstract: We show that the BMS-supertranslations and their associated supermomenta on past null infinity can be related to those on future null infinity, proving the conjecture of Strominger for a class of spacetimes which are asymptotically-flat in the sense of Ashtekar and Hansen. Using a cylindrical 3-manifold of both null and spatial directions of approach towards spatial infinity, we impose appropriate regularity conditions on the Weyl tensor near spatial infinity along null directions. The asymptotic Einstein equations on this 3-manifold and the regularity conditions imply that the relevant Weyl tensor components on past null infinity are antipodally matched to those on future null infinity. The subalgebra of totally fluxless supertranslations near spatial infinity provides a natural isomorphism between the BMS-supertranslations on past and future null infinity. This proves that the flux of the supermomenta is conserved from past to future null infinity in a classical gravitational scattering process provided additional suitable conditions are satisfied at the timelike infinities.


title: Blind chemical tagging with DBSCAN: prospects for spectroscopic surveys

abstract: Chemical tagging has great promise as a technique to unveil our Galaxy's history. Grouping stars based on their similar chemistry can establish details of the star formation and merger history of the Milky Way. With precise measurements of stellar chemistry, chemical tagging may be able to group together stars born from the same gas cloud, regardless of their current positions and kinematics. Successfully tagging these birth clusters requires high quality chemical space information and a good cluster-finding algorithm. To test the feasibility of chemical tagging on data from current and upcoming spectroscopic surveys, we construct a realistic set of synthetic clusters, creating both observed spectra and derived chemical abundances for each star. We use Density-Based Spatial Clustering of Applications with Noise (DBSCAN) to group stars based on their spectra or abundances; these groups are matched to input clusters and are found to be highly homogeneous and complete. The percentage of clusters with more than 10 members recovered is 40% when tagging on abundances with uncertainties achievable with current techniques. Based on our fiducial model for the Milky Way, we predict recovering over 600 clusters with at least 10 observed members and 70% membership homogeneity in a sample similar to the APOGEE survey. Tagging larger surveys like the GALAH survey and the future Milky Way Mapper in SDSS V could recover tens of thousands of clusters at high homogeneity. Access to so many unique co-eval clusters will transform how we understand the star formation history and chemical evolution of our Galaxy.


title: Hydrodynamical moving-mesh simulations of the tidal disruption of stars by supermassive black holes

abstract: When a star approaches a black hole closely, it may be pulled apart by gravitational forces in a tidal disruption event (TDE). The flares produced by TDEs are unique tracers of otherwise quiescent supermassive black holes (SMBHs) located at the centre of most galaxies. In particular, the appearance of such flares and the subsequent decay of the light curve are both sensitive to whether the star is partially or totally destroyed by the tidal field. However, the physics of the disruption and the fall-back of the debris are still poorly understood. We are here modelling the hydrodynamical evolution of realistic stars as they approach a SMBH on parabolic orbits, using for the first time the moving-mesh code AREPO, which is particularly well adapted to the problem through its combination of quasi-Lagrangian behaviour, low advection errors, and high accuracy typical of mesh-based techniques. We examine a suite of simulations with different impact parameters, allowing us to determine the critical distance at which the star is totally disrupted, the energy distribution and the fallback rate of the debris, as well as the hydrodynamical evolution of the stellar remnant in the case of a partial disruption. Interestingly, we find that the internal evolution of the remnant's core is strongly influenced by persistent vortices excited in the tidal interaction. These should be sites of strong magnetic field amplification, and the associated mixing may profoundly alter the subsequent evolution of the tidally pruned star.


title: TeV Halos are Everywhere: Prospects for New Discoveries

abstract: Milagro and HAWC have detected extended TeV gamma-ray emission around nearby pulsar wind nebulae (PWNe). Building on these discoveries, Linden et al. [1] identified a new source class -- TeV halos -- powered by the interactions of high-energy electrons and positrons that have escaped from the PWN, but which remain trapped in a larger region where diffusion is inhibited compared to the interstellar medium. Many theoretical properties of TeV halos remain mysterious, but empirical arguments suggest that they are ubiquitous. The key to progress is finding more halos. We outline prospects for new discoveries and calculate their expectations and uncertainties. We predict, using models normalized to current data, that future HAWC and CTA observations will detect in total $\sim$50--240 TeV halos, though we note that multiple systematic uncertainties still exist. Further, the existing HESS source catalog could contain $\sim$10--50 TeV halos that are presently classified as unidentified sources or PWN candidates. We quantify the importance of these detections for new probes of the evolution of TeV halos, pulsar properties, and the sources of high-energy gamma rays and cosmic rays.


title: Neutral Gas Properties and Ly$\alpha$ Escape in Extreme Green Pea Galaxies

abstract: Mechanisms regulating the escape of Ly$\alpha$ photons and ionizing radiation remain poorly understood. To study these processes we analyze VLA 21cm observations of one Green Pea (GP), J160810+352809 (hereafter J1608), and HST COS spectra of 17 GP galaxies at $z<0.2$. All are highly ionized: J1608 has the highest [O III] $\lambda5007$/[O II] $\lambda3727$ for star-forming galaxies in SDSS, and the 17 GPs have [O III]/[O II] $\geq6.6$. We set an upper limit on J1608's HI mass of $\log M_{HI}/M_\odot=8.4$, near or below average compared to similar mass dwarf galaxies. In the COS sample, eight GPs show Ly$\alpha$ absorption components, six of which also have Ly$\alpha$ emission. The HI column densities derived from Ly$\alpha$ absorption are high, $\log N_{HI}/$cm$^{-2}=19-21$, well above the LyC optically thick limit. Using low-ionization absorption lines, we measure covering fractions ($f_{\mbox{cov}}$) of $0.1-1$, and find that $f_{\mbox{cov}}$ strongly anti-correlates with Ly$\alpha$ escape fraction. Low covering fractions may facilitate Ly$\alpha$ and LyC escape through dense neutral regions. GPs with $f_{\mbox{cov}}\sim1$ all have low neutral gas velocities, while GPs with lower $f_{\mbox{cov}}=0.2-0.6$ have a larger range of velocities. Conventional mechanical feedback may help establish low $f_{\mbox{cov}}$ in some cases, whereas other processes may be important for GPs with low velocities. Finally, we compare $f_{\mbox{cov}}$ with proposed indicators of LyC escape. Ionizing photon escape likely depends on a combination of neutral gas geometry and kinematics, complicating the use of emission-line diagnostics for identifying LyC emitters.


title: Particle acceleration and synchrotron self-Compton emission in blazar jets I : an application to the quiescent emission

abstract: There are still some important unanswered questions about the detailed particle acceleration and escape occurring during the quiescent epoches. As a result, the particle distribution that is adopted in the blazar quiescent spectral model have numerous unconstrained shapes. To help remedy this problem, we introduce a analytical particle transport model to reproduce quiescent broadband spectral energy distribution of blazar. In this model, the exact electron distribution is solved from a generalized transport equation that contains the terms describing first-order and secondary-order \emph{Fermi} acceleration, escape of particle due to both the advection and spatial diffusion, energy losses due to synchrotron emission and inverse-Compton scattering of an assumed soft photon field. We suggest that the advection is a significant escape mechanism in blazar jet. We find that in our model the advection process tends to harden the particle distribution, which enhances the high energy components of resulting synchrotron and synchrotron self-Comptom spectrum from jet. Our model is able to roughly reproduce the observed spectra of extreme BL Lac object 1ES 0414+009 with reasonable assumptions about the physical parameters.


title: Unseen companions of V Hya inferred from periodic ejections

abstract: A recent study using $Hubble$ $Space$ $Telescope$ observations found periodic, high-speed, collimated ejections (or "bullets") from the star V Hya. The authors of that study proposed a model associating these bullets with the periastron passage of an unseen, substellar companion in an eccentric orbit and with an orbital period of $\sim8$ yrs. Here we propose that V Hya is part of a triple system, with a substellar companion having an orbital period of $\sim8$ yrs, and a tertiary object on a much wider orbit. In this model, the more distant object causes high-eccentricity excitations on the substellar companion's orbit via the Eccentric Kozai-Lidov mechanism. These eccentricities can reach such high values that they lead to Roche-lobe crossing, producing the observed bullet ejections via a strongly enhanced accretion episode. For example, we find that a ballistic bullet ejection mechanism can be produced by a brown-dwarf-mass companion, while magnetically driven outflows are consistent with a Jovian-mass companion. Finally, we suggest that the distant companion may reside at few a hundred AU on an eccentric orbit.


title: Observational Signatures of Quantum Gravity in Interferometers

abstract: We consider the uncertainty in the arm length of an interferometer due to metric fluctuations from the quantum nature of gravity, proposing a concrete microscopic model of energy fluctuations in holographic degrees of freedom on the surface bounding a causally connected region of spacetime. In our model, fluctuations longitudinal to the beam direction accumulate in the infrared and feature strong long distance correlation in the transverse direction. This leads to a signal that could be observed in a gravitational wave interferometer. We connect the positional uncertainty principle arising from our calculations to the 't Hooft gravitational S-matrix.


title: Quark and gluon form factors in four loop QCD: the $N_f^2$ and $N_{q\gamma} N_f$ contributions

abstract: We calculate the four-loop massless QCD corrections with two closed quark lines to quark and gluon form factors. The results for the gluon form factor and the singlet part of the quark form factor are given for the first time. From our analytic expressions for the form factors, we determine the corresponding cusp anomalous dimensions. The relevant Feynman integrals are obtained with novel integral reduction techniques and direct integration methods.


title: No classical (anti-)de Sitter solutions with O8-planes

abstract: It was recently proposed that type IIA string theory may allow classical de Sitter solutions with O8-planes as the only localized sources. We show that such solutions are incompatible with the integrated supergravity equations of motion, analogously to the no-go theorem due to Maldacena and Nu\~{n}ez. We also discuss in detail divergences and discontinuities at the O8-plane positions and argue that they do not invalidate such an argument. We furthermore show that a recently proposed class of non-supersymmetric AdS solutions with O8-planes is in contrast with our results as well.


title: EVN observations of 6.7 GHz methanol maser polarization in massive star-forming regions IV. Magnetic field strength limits and structure for 7 additional sources

abstract: MHD simulations show that the magnetic field can drive molecular outflows during the formation of massive protostars. The best probe to observationally measure both the morphology and the strength of this magnetic field at scales of 10-100 au is maser polarization. We measure the direction of magnetic fields at mas resolution around a sample of massive star-forming regions to determine whether there is a relation between the orientation of the magnetic field and of the outflows. In addition, by estimating the magnetic field strength via the Zeeman splitting measurements, the role of magnetic field in the dynamics of the massive star-forming region is investigated. We selected a flux-limited sample of 31 massive star-forming regions to perform a statistical analysis of the magnetic field properties with respect to the molecular outflows characteristics. We report the linearly and circularly polarized emission of 6.7 GHz CH3OH masers towards seven massive star-forming regions of the total sample with the EVN. The sources are: G23.44-0.18, G25.83-0.18, G25.71-0.04, G28.31-0.39, G28.83-0.25, G29.96-0.02, and G43.80-0.13. We identified a total of 219 CH3OH maser features, 47 and 2 of which showed linearly and circularly polarized emission, respectively. We measured well-ordered linear polarization vectors around all the massive YSOs and Zeeman splitting towards G25.71-0.04 and G28.83-0.25. Thanks to recent theoretical results, we were able to provide lower limits to the magnetic field strength from our Zeeman splitting measurements. We further confirm (based on ~80% of the total flux-limited sample) that the magnetic field on scales of 10-100 au is preferentially oriented along the outflow axes. The estimated magnetic field strength of |B_{||}|>61 mG and >21 mG towards G25.71-0.04 and G28.83-0.2, respectively, indicates that it dominates the dynamics of the gas in both regions.


title: Composite Twin Dark Matter

abstract: We consider Fraternal Twin Higgs models where the twin bottom quark, $b'$, is much heavier than the twin confinement scale. In this limit aspects of quark bound states, like the mass and binding energy, can be accurately calculated. We show that in this regime, dark matter can be primarily made of twin baryons containing $b' b' b'$ or, when twin hypercharge is gauged, twin atoms, composed of a baryon bound to a twin $\tau'$ lepton. We find that there are significant regions of parameter space which are allowed by current constraints but within the realm of detection in the near future. The case with twin atoms can alleviate the tension between dark matter properties inferred from dwarf galaxies and clusters.


title: Relaxion Stars and their detection via Atomic Physics

abstract: The cosmological relaxion can address the hierarchy problem, while its coherent oscillations can constitute dark matter in the present universe. We consider the possibility that the relaxion forms gravitationally bound objects that we denote as relaxion stars. The density of these stars would be higher than that of the local dark matter density, resulting in enhanced signals in table-top detectors, among others. Furthermore, we raise the possibility that these objects may be trapped by an external gravitational potential, such as that of the Earth or the Sun. This leads to formation of relaxion halos of even greater density. We discuss several interesting implications of relaxion halos, as well as detection strategies to probe them.


title: Towards Visually Grounded Sub-Word Speech Unit Discovery

abstract: In this paper, we investigate the manner in which interpretable sub-word speech units emerge within a convolutional neural network model trained to associate raw speech waveforms with semantically related natural image scenes. We show how diphone boundaries can be superficially extracted from the activation patterns of intermediate layers of the model, suggesting that the model may be leveraging these events for the purpose of word recognition. We present a series of experiments investigating the information encoded by these events.


title: Statistics of Square-tiled Surfaces: Symmetry and Short Loops

abstract: Square-tiled surfaces are a class of translation surfaces that are of particular interest in geometry and dynamics because, as covers of the square torus, they share some of its simplicity and structure. In this paper, we study counting problems that result from focusing on properties of the square torus one by one. After drawing insights from experimental evidence, we consider the implications between these properties and their frequency in each stratum of translation surfaces.


title: Modularity from Monodromy

abstract: In this note we describe a method to calculate the action of a particular Fourier-Mukai transformation on a basis of brane charges on elliptically fibered Calabi-Yau threefolds with and without a section. The Fourier-Mukai kernel is the ideal sheaf of the relative diagonal and for fibrations that admit a section this is essentially the Poincar\'e sheaf. We find that in this case it induces an action of the modular group on the charges of 2-branes.


title: Environmental Effect on the Interstellar Medium in Galaxies across the Cosmic Web at z=0.73

abstract: We present new ALMA dust continuum observations of 101 $\log(\mathrm{M}_* / \mathrm{M}_\odot)$ > 9.5 galaxies in the COSMOS field to study the effect of environment on the interstellar medium at z ~ 0.7. At this redshift, our targets span a wide range of environments allowing for a diverse sample of galaxies with densities, $\Sigma$ = 0.16-10.5 Mpc$^{-2}$ (per $\Delta$ z = 0.024). Using the ALMA observations, we calculate the total ISM mass and look for depletion as a function of galaxy density in order to understand the quenching or triggering of star formation in galaxies in different environments. ISM mass is found to have a small dependence on environment, while the depletion timescale remains constant (~200 Myrs) across all environments. We find elevated ISM mass values at intermediate densities and lower values at high densities compared to low (field) densities. Our observed evolution in gas fraction with density in this single redshift slice is equivalent to the observed evolution with cosmic time over 2-3 Gyr. To explain the change in gas mass fraction seen in galaxies in intermediate and high densities, these results suggest environmental processes such as mergers and ram pressure stripping are likely playing a role in dense filamentary-cluster environments.


title: Ultra-low loss integrated visible photonics using thin-film lithium niobate

abstract: Integrated photonics is a powerful platform that can improve the performance and stability of optical systems, while providing low-cost, small-footprint and scalable alternatives to implementations based on free-space optics. While great progress has been made on the development of low-loss integrated photonics platforms at telecom wavelengths, visible wavelength range has received less attention. Yet, many applications utilize visible or near-visible light, including those in optical imaging, optogenetics, and quantum science and technology. Here we demonstrate an ultra-low loss integrated visible photonics platform based on thin film lithium niobate on insulator. Our waveguides feature ultra-low propagation loss of 6 dB/m, while our microring resonators have an intrinsic quality factor of 11 million, both measured at 637 nm wavelength. Additionally, we demonstrate an on-chip visible intensity modulator with an electro-optic bandwidth of 10 GHz, limited by the detector used. The ultra-low loss devices demonstrated in this work, together with the strong second- and third-order nonlinearities in lithium niobate, open up new opportunities for creating novel passive, and active devices for frequency metrology and quantum information processing in the visible spectrum range.


title: Covering a tree with rooted subtrees

abstract: We consider the multiple traveling salesman problem on a weighted tree. In this problem there are $m$ salesmen located at the root initially. Each of them will visit a subset of vertices and return to the root. The goal is to assign a tour to every salesman such that every vertex is visited and the longest tour among all salesmen is minimized. The problem is equivalent to the subtree cover problem, in which we cover a tree with rooted subtrees such that the weight of the maximum weighted subtree is minimized. The classical machine scheduling problem can be viewed as a special case of our problem when the given tree is a star. We observe that, the problem remains NP-hard even if tree height and edge weight are constant, and present an FPT algorithm for this problem parameterized by the largest tour length. To achieve the FPT algorithm, we show a more general result. We prove that, integer linear programming that has a tree-fold structure is in FPT, which extends the FPT result for the $n$-fold integer programming by Hemmecke, Onn and Romanchuk.


title: An Approach to Stability Analyses in General Relativity via Symplectic Geometry

abstract: We begin with a review of the statements of non-linear, linear and mode stability of autonomous dynamical systems in classical mechanics, using symplectic geometry. We then discuss what the phase space and the Hamiltonian of general relativity are, what constitutes a dynamical system, and subsequently draw a formal analogy between the notions of stability in these two theories. Our approach here is pedagogical and geometric, and considerably simplifies a formal understanding of the statements regarding the stability of stationary solutions of general relativity. In particular, the governing equations of motion of a Hamiltonian dynamical system are simply the flow equations of the associated symplectic Hamiltonian vector field, defined on phase space, and the non-linear stability analysis of its critical points have simply to do with the divergence of its flow there. Further, the linear stability of a critical point is related to the properties of the tangent flow of the Hamiltonian vector field. Further, we posit that a study of the genericity of a particular black hole or naked singularity spacetime forming as an endstate of gravitational collapse is equivalent to an inquiry of how sensitive the orbits of the symplectic Hamiltonian vector field of general relativity are to changes in initial data. We demonstrate this by conducting a restricted non-linear stability analysis of the formation of a Schwarzschild black hole, working in the usual initial value formulation of general relativity.


title: Existence of Solutions to Nonlinear Legendre Boundary Value Problems

abstract: In this paper, we consider nonlinearly perturbed Legendre differential equations subject to the usual boundary conditions. For such problems we establish sufficient conditions for the existence of solutions and in some cases we provide a qualitative description of solutions depending on a parameter. The results presented depend on the size and limiting behavior of the nonlinearities.


title: Compact minimal hypersurfaces of index one and the width of real projective spaces

abstract: We characterize the first min-max width of real projective spaces of any dimension. The width is the minimum area over the Clifford hypersurfaces. We also compute the Morse index of the Clifford hypersurfaces in the complex and quaternionic projective spaces.


title: Learning requirements for stealth attacks

abstract: The learning data requirements are analyzed for the construction of stealth attacks in state estimation. In particular, the training data set is used to compute a sample covariance matrix that results in a random matrix with a Wishart distribution. The ergodic attack performance is defined as the average attack performance obtained by taking the expectation with respect to the distribution of the training data set. The impact of the training data size on the ergodic attack performance is characterized by proposing an upper bound for the performance. Simulations on the IEEE 30-Bus test system show that the proposed bound is tight in practical settings.


title: A Note on Covering Young Diagrams with Applications to Local Dimension of Posets

abstract: We prove that in every cover of a Young diagram with $\binom{2k}{k}$ steps with generalized rectangles there is a row or a column in the diagram that is used by at least $k+1$ rectangles. We show that this is best-possible by partitioning any Young diagram with $\binom{2k}{k}-1$ steps into actual rectangles, each row and each column used by at most $k$ rectangles. This answers two questions by Kim et al. (2018). Our results can be rephrased in terms of local covering numbers of difference graphs with complete bipartite graphs, which has applications in the recent notion of local dimension of partially ordered sets.


title: Blind Hyperspectral-Multispectral Image Fusion via Graph Laplacian Regularization

abstract: Fusing a low-resolution hyperspectral image (HSI) and a high-resolution multispectral image (MSI) of the same scene leads to a super-resolution image (SRI), which is information rich spatially and spectrally. In this paper, we super-resolve the HSI using the graph Laplacian defined on the MSI. Unlike many existing works, we don't assume prior knowledge about the spatial degradation from SRI to HSI, nor a perfectly aligned HSI and MSI pair. Our algorithm progressively alternates between finding the blur kernel and fusing HSI with MSI, generating accurate estimations of the blur kernel and the SRI at convergence. Experiments on various datasets demonstrate the advantages of the proposed algorithm in the quality of fusion and its capability in dealing with unknown spatial degradation.


title: Flexible Conductive Composites with Programmed Electrical Anisotropy Using Acoustophoresis

abstract: 3D printing mechanically flexible composite materials with high electrical conductivity is currently hindered by the need to use high loading of conductive filler, which severely limits flexibility. Here, microstructural patterning of composite materials via acoustophoresis imparts these materials with high conductivity and flexibility simultaneously, filling a technology gap in the field. Acoustophoresis patterns filler particles into highly efficient percolated networks which utilize up to 97\% of the particles in the composite, whereas the inefficient stochastic networks of conventional dispersed-fiber composites utilize $<5$\%. These patterned materials have conductivity an order of magnitude higher than conventional composites made with the same ink, reaching 48\% the conductivity of bulk silver within the assembled silver-particle networks (at 2.6v\% loading). They also have low particle loading so that they're flexible, withstanding $>$500 bending cycles without losses in conductivity and changing conductivity only $5$\% within cycles on average (for 2.6v\% composites). In contrast, conventional unpatterned composites with the same conductivity require such high loading that they're prohibitively brittle. Finally, modulating the shape of the applied acoustic fields allows control over the anisotropy of the conductive networks and produces materials which are either 2-D conductive, 1-D conductive, or insulating, all using the same nozzle and ink.


title: Graph Adversarial Training: Dynamically Regularizing Based on Graph Structure

abstract: Recent efforts show that neural networks are vulnerable to small but intentional perturbations on input features in visual classification tasks. Due to the additional consideration of connections between examples (\eg articles with citation link tend to be in the same class), graph neural networks could be more sensitive to the perturbations, since the perturbations from connected examples exacerbate the impact on a target example. Adversarial Training (AT), a dynamic regularization technique, can resist the worst-case perturbations on input features and is a promising choice to improve model robustness and generalization. However, existing AT methods focus on standard classification, being less effective when training models on graph since it does not model the impact from connected examples. In this work, we explore adversarial training on graph, aiming to improve the robustness and generalization of models learned on graph. We propose Graph Adversarial Training (GraphAT), which takes the impact from connected examples into account when learning to construct and resist perturbations. We give a general formulation of GraphAT, which can be seen as a dynamic regularization scheme based on the graph structure. To demonstrate the utility of GraphAT, we employ it on a state-of-the-art graph neural network model --- Graph Convolutional Network (GCN). We conduct experiments on two citation graphs (Citeseer and Cora) and a knowledge graph (NELL), verifying the effectiveness of GraphAT which outperforms normal training on GCN by 4.51% in node classification accuracy. Codes are available via: https://github.com/fulifeng/GraphAT.


title: Slowing the Spins of Stellar Cores

abstract: The angular momentum (AM) evolution of stellar interiors, along with the resulting rotation rates of stellar remnants, remains poorly understood. Asteroseismic measurements of red giant stars reveal that their cores rotate much faster than their surfaces, but much slower than theoretically predicted, indicating an unidentified source of AM transport operates in their radiative cores. Motivated by this, we investigate the magnetic Tayler instability and argue that it saturates when turbulent dissipation of the perturbed magnetic field energy is equal to magnetic energy generation via winding. This leads to larger magnetic field amplitudes, more efficient AM transport, and smaller shears than predicted by the classic Tayler-Spruit dynamo. We provide prescriptions for the effective AM diffusivity and incorporate them into numerical stellar models, finding they largely reproduce (1) the nearly rigid rotation of the Sun and main sequence stars, (2) the core rotation rates of low-mass red giants during hydrogen shell and helium burning, and (3) the rotation rates of white dwarfs. We discuss implications for stellar rotational evolution, internal rotation profiles, rotational mixing, and the spins of compact objects.


title: A Comprehensive Theory and Variational Framework for Anti-aliasing Sampling Patterns

abstract: In this paper, we provide a comprehensive theory of anti-aliasing sampling patterns that explains and revises known results, and show how patterns as predicted by the theory can be generated via a variational optimization framework. We start by deriving the exact spectral expression for expected error in reconstructing an image in terms of power spectra of sampling patterns, and analyzing how the shape of power spectra is related to anti-aliasing properties. Based on this analysis, we then formulate the problem of generating anti-aliasing sampling patterns as constrained variational optimization on power spectra. This allows us to not rely on any parametric form, and thus explore the whole space of realizable spectra. We show that the resulting optimized sampling patterns lead to reconstructions with less visible aliasing artifacts, while keeping low frequencies as clean as possible.


title: Towards More Flexible False Positive Control in Phase III Randomized Clinical Trials

abstract: Phase III randomized clinical trials play a monumentally critical role in the evaluation of new medical products. Because of the intrinsic nature of uncertainty embedded in our capability in assessing the efficacy of a medical product, interpretation of trial results relies on statistical principles to control the error of false positives below desirable level. The well-established statistical hypothesis testing procedure suffers from two major limitations, namely, the lack of flexibility in the thresholds to claim success and the lack of capability of controlling the total number of false positives that could be yielded by the large volume of trials. We propose two general theoretical frameworks based on the conventional frequentist paradigm and Bayesian perspectives, which offer realistic, flexible and effective solutions to these limitations. Our methods are based on the distribution of the effect sizes of the population of trials of interest. The estimation of this distribution is practically feasible as clinicaltrials.gov provides a centralized data repository with unbiased coverage of clinical trials. We provide a detailed development of the two frameworks with numerical results obtained for industry sponsored Phase III randomized clinical trials.


title: Meridional Structure and Future Changes of Tropopause Height and Temperature

abstract: We use a simple, semi-analytic, column model to better understand the meridional structure of the tropopause height and the future changes in its height and temperature associated with global warming. The model allows us to separate the effects of tropospheric lapse rate (TLR), optical depth, outgoing longwave radiation (OLR) and stratospheric cooling on the tropopause height. When applied locally at each latitudinal band the model predicts the overall meridional structure of the tropopause height, with a tropical tropopause substantially higher than that in higher latitudes and a sharp transition at the edge of the extratropics. The large optical depth of the tropics, due mainly to the large water vapour path (WVP), is the dominant tropospheric effect producing the higher tropical tropopause, whereas the larger tropical lapse rate acts to lower the tropopause height. The dynamical cooling induced by the stratospheric circulation further lifts the thermal tropopause in the tropics resulting in it being significantly cooler and higher than in mid- and high latitudes. The model quantifies the causes of the tropopause height increase with global warming that is robustly found in climate integrations from the fifth Coupled Model Intercomparison Project (CMIP5). The large spread in the increase rate of tropopause height in the CMIP5 models is captured by the simple model, which attributes the dominant contributions to changes in WVP and TLR, with changes in CO2 concentration and OLR having much smaller effects. The CMIP5 models also show a small but robust increase in the tropopause temperature in low latitudes, with a much smaller increase in higher latitudes. We suggest that the tropical increase may in part be caused by non-grey effects in the radiative transfer associated with the higher levels of water vapour in the tropics, with near constant tropopause temperatures predicted otherwise.


title: Online Multi-Object Tracking with Instance-Aware Tracker and Dynamic Model Refreshment

abstract: Recent progresses in model-free single object tracking (SOT) algorithms have largely inspired applying SOT to \emph{multi-object tracking} (MOT) to improve the robustness as well as relieving dependency on external detector. However, SOT algorithms are generally designed for distinguishing a target from its environment, and hence meet problems when a target is spatially mixed with similar objects as observed frequently in MOT. To address this issue, in this paper we propose an instance-aware tracker to integrate SOT techniques for MOT by encoding awareness both within and between target models. In particular, we construct each target model by fusing information for distinguishing target both from background and other instances (tracking targets). To conserve uniqueness of all target models, our instance-aware tracker considers response maps from all target models and assigns spatial locations exclusively to optimize the overall accuracy. Another contribution we make is a dynamic model refreshing strategy learned by a convolutional neural network. This strategy helps to eliminate initialization noise as well as to adapt to the variation of target size and appearance. To show the effectiveness of the proposed approach, it is evaluated on the popular MOT15 and MOT16 challenge benchmarks. On both benchmarks, our approach achieves the best overall performances in comparison with published results.


title: Overcoming Multi-Model Forgetting

abstract: We identify a phenomenon, which we refer to as multi-model forgetting, that occurs when sequentially training multiple deep networks with partially-shared parameters; the performance of previously-trained models degrades as one optimizes a subsequent one, due to the overwriting of shared parameters. To overcome this, we introduce a statistically-justified weight plasticity loss that regularizes the learning of a model's shared parameters according to their importance for the previous models, and demonstrate its effectiveness when training two models sequentially and for neural architecture search. Adding weight plasticity in neural architecture search preserves the best models to the end of the search and yields improved results in both natural language processing and computer vision tasks.


title: A pore-scale study of transport of inertial particles by water in porous media

abstract: We study the transport of inertial particles in water flow in porous media. Our interest lies in understanding the accumulation of particles including the possibility of clogging. We propose that accumulation can be a result of hydrodynamic effects: the tortuous paths of the porous medium generate regions of dominating strain/vorticity, which favour the accumulation/dispersion of the inertial particles. Numerical simulations show that essentially two accumulation regimes are identified: for low and for high flow velocities. When particles accumulate in high-velocity regions, at the entrance of a pore throat, a clog is formed. The formation of a clog significantly modifies the flow, as the partial blockage of the pore causes a local redistribution of pressure. This redistribution can divert the upstream water flow into neighbouring pores. Moreover, we show that accumulation in high velocity regions occurs in heterogeneous media, but not in homogeneous media, where we refer to homogeneity with respect to the distribution of the pore throat diameters.


title: An Empirical Study of Large-Batch Stochastic Gradient Descent with Structured Covariance Noise

abstract: The choice of batch-size in a stochastic optimization algorithm plays a substantial role for both optimization and generalization. Increasing the batch-size used typically improves optimization but degrades generalization. To address the problem of improving generalization while maintaining optimal convergence in large-batch training, we propose to add covariance noise to the gradients. We demonstrate that the learning performance of our method is more accurately captured by the structure of the covariance matrix of the noise rather than by the variance of gradients. Moreover, over the convex-quadratic, we prove in theory that it can be characterized by the Frobenius norm of the noise matrix. Our empirical studies with standard deep learning model-architectures and datasets shows that our method not only improves generalization performance in large-batch training, but furthermore, does so in a way where the optimization performance remains desirable and the training duration is not elongated.


title: Natural quark mixing and inverse seesaw in a left-right model with an axion

abstract: We consider a minimal left-right model with a Peccei-Quinn symmetry, where generalised charge conjugation plays the role of the left-right symmetry. We show how the spontaneous breaking of the Peccei-Quinn symmetry by a scalar singlet can provide us with solutions not only to the strong CP and dark matter problems but can also help to generate naturally suppressed off-diagonal CKM elements and small neutrino masses via the inverse seesaw mechanism. For this, we make use of an economical scalar sector composed of a bi-doublet, two doublets and a singlet only. As a result of the new gauge bosons and neutrinos, the neutrinoless double beta decay, as well as lepton flavour violating processes, receives new contributions which can, in principle, become relevant due to the low-scale nature of the inverse seesaw mechanism. The model can easily accommodate all the current data on fermion masses and mixing even if the left-right scale is only high enough to evade the current experimental constraints.


title: Lung Cancer Detection using Co-learning from Chest CT Images and Clinical Demographics

abstract: Early detection of lung cancer is essential in reducing mortality. Recent studies have demonstrated the clinical utility of low-dose computed tomography (CT) to detect lung cancer among individuals selected based on very limited clinical information. However, this strategy yields high false positive rates, which can lead to unnecessary and potentially harmful procedures. To address such challenges, we established a pipeline that co-learns from detailed clinical demographics and 3D CT images. Toward this end, we leveraged data from the Consortium for Molecular and Cellular Characterization of Screen-Detected Lesions (MCL), which focuses on early detection of lung cancer. A 3D attention-based deep convolutional neural net (DCNN) is proposed to identify lung cancer from the chest CT scan without prior anatomical location of the suspicious nodule. To improve upon the non-invasive discrimination between benign and malignant, we applied a random forest classifier to a dataset integrating clinical information to imaging data. The results show that the AUC obtained from clinical demographics alone was 0.635 while the attention network alone reached an accuracy of 0.687. In contrast when applying our proposed pipeline integrating clinical and imaging variables, we reached an AUC of 0.787 on the testing dataset. The proposed network both efficiently captures anatomical information for classification and also generates attention maps that explain the features that drive performance.


title: Global Hypoellipticity for Strongly Invariant Operators

abstract: In this note, by analyzing the behavior at infinity of the matrix symbol of an invariant operator $P$ with respect to a fixed elliptic operator, we obtain a necessary and sufficient condition to guarantee that $P$ is globally hypoelliptic. We also investigate relations between the global hypoellipticity of $P$ and global subelliptic estimates.


title: Open source software for simulations and inversions of airborne electromagnetic data

abstract: Inversions of airborne EM data are often an iterative process, not only requiring that the researcher be able to explore the impact of changing components such as the choice of regularization functional or model parameterization, but also often requiring that forward simulations be run and fields and fluxes visualized in order to build an understanding of the physical processes governing what we observe in the data. In the hope of facilitating this exploration and promoting reproducibility of geophysical simulations and inversions, we have developed the open source software package, SimPEG. The software has been designed to be modular and extensible with the goal of allowing researchers to interrogate all of the components and to facilitate the exploration of new inversion strategies. We present an overview of the software in its application to airborne EM and demonstrate its use for visualizing fields and fluxes in a forward simulation as well as its flexibility in formulating and solving the inverse problem. We invert a line of airborne TDEM data over a conductive vertical plate using a 1D voxel-inversion, a 2D voxel inversion and a parametric inversion, where all of the forward modelling is done on a 3D grid. The results in this paper can be reproduced by using the provided Jupyter notebooks. The Python software can also be modified to allow users to experiment with parameters and explore the physics of the electromagnetics and intricacies of inversion.


title: Examples of (non)-braided tensor categories

abstract: Are introduced six examples of non-braidable tensor categories which are extensions of the category Comod(H), for H a super-group algebra; and two examples of braided categories where the only possible braiding is the trivial braiding.


title: Iterative oscillation tests for differential equations with several non-monotone arguments

abstract: Sufficient oscillation conditions involving $\limsup $ and $\liminf $ for first-order differential equations with several non-monotone deviating arguments and nonnegative coefficients are obtained. The results are based on the iterative application of the Gr\"{o}nwall inequality. Examples illustrating the significance of the results are also given.


title: A high precision calibration of the nonlinear energy response at Daya Bay

abstract: A high precision calibration of the nonlinearity in the energy response of the Daya Bay Reactor Neutrino Experiment's antineutrino detectors is presented in detail. The energy nonlinearity originates from the particle-dependent light yield of the scintillator and charge-dependent electronics response. The nonlinearity model is constrained by $\gamma$ calibration points from deployed and naturally occurring radioactive sources, the $\beta$ spectrum from $^{12}$B decays, and a direct measurement of the electronics nonlinearity with a new flash analog-to-digital converter readout system. Less than 0.5% uncertainty in the energy nonlinearity calibration is achieved for positrons of kinetic energies greater than 1 MeV.


title: Jamming transition in non-spherical particle systems: pentagons vs. disks

abstract: We investigate the jamming transition in a quasi-2D granular material composed of regular pentagons or disks subjected to quasistatic uniaxial compression. We report six major findings based on experiments with monodisperse photoelastic particles with static friction coefficient $\mu\approx 1$. (1) For both pentagons and disks, the onset of rigidity occurs when the average coordination number of non-rattlers, $Z_{nr}$ , reaches 3, and the dependence of $Z_{nr}$ on the packing fraction $\phi$ changes again when $Z_{nr}$ reaches 4. (2) Though the packing fractions $\phi_{c1}$ and $\phi_{c2}$ at these transitions differ from run to run, for both shapes the data from all runs with different initial configurations collapses when plotted as a function of the non-rattler fraction. (3) The averaged values of $\phi_{c1}$ and $\phi_{c2}$ for pentagons are around 1% smaller than those for disks. (4) Both jammed pentagons and disks show Gamma distribution of the Voronoi cell area with same parameters. (5) The jammed pentagons have similar translational order for particle centers but slightly less orientational order for contacting pairs comparing to jammed disks. (6) For jammed pentagons, the angle between edges at a face-to-vertex contact point shows a uniform distribution and the size of a cluster connected by face-to-face contacts shows a power-law distribution.


title: Critical jammed phase of the linear perceptron

abstract: Criticality in statistical physics naturally emerges at isolated points in the phase diagram. Jamming of spheres is not an exception: varying density, it is the critical point that separates the unjammed phase where spheres do not overlap and the jammed phase where they cannot be arranged without overlaps. The same remains true in more general constraint satisfaction problems with continuous variables (CCSP) where jamming coincides with the (protocol dependent) satisfiability transition point. In this work we show that by carefully choosing the cost function to be minimized, the region of criticality extends to occupy a whole region of the jammed phase. As a working example, we consider the spherical perceptron with a linear cost function in the unsatisfiable (UNSAT) jammed phase and we perform numerical simulations which show critical power laws emerging in the configurations obtained minimizing the linear cost function. We develop a scaling theory to compute the emerging critical exponents.


title: Optimizing Sparse RFI Prediction using Deep Learning

abstract: Radio Frequency Interference (RFI) is an ever-present limiting factor among radio telescopes even in the most remote observing locations. When looking to retain the maximum amount of sensitivity and reduce contamination for Epoch of Reionization studies, the identification and removal of RFI is especially important. In addition to improved RFI identification, we must also take into account computational efficiency of the RFI-Identification algorithm as radio interferometer arrays such as the Hydrogen Epoch of Reionization Array grow larger in number of receivers. To address this, we present a Deep Fully Convolutional Neural Network (DFCN) that is comprehensive in its use of interferometric data, where both amplitude and phase information are used jointly for identifying RFI. We train the network using simulated HERA visibilities containing mock RFI, yielding a known "ground truth" dataset for evaluating the accuracy of various RFI algorithms. Evaluation of the DFCN model is performed on observations from the 67 dish build-out, HERA-67, and achieves a data throughput of 1.6$\times 10^{5}$ HERA time-ordered 1024 channeled visibilities per hour per GPU. We determine that relative to an amplitude only network including visibility phase adds important adjacent time-frequency context which increases discrimination between RFI and Non-RFI. The inclusion of phase when predicting achieves a Recall of 0.81, Precision of 0.58, and $F_{2}$ score of 0.75 as applied to our HERA-67 observations.


title: 3D electromagnetic modelling and inversion: A case for open-source

abstract: Electromagnetics has an important role to play in solving the next generation of geoscience problems. These problems are multidisciplinary, complex, and require collaboration. This is especially true at the base scientific level where the underlying physical equations need to be solved, and data, associated with physical experiments, need to be inverted. In this paper, we present arguments for adopting an open-source methodology for geophysics and provide some background about open-source software for electromagnetics. Immediate benefits are the reduced time required to carry out research, being able to collaborate, having reproducible results, and being able to disseminate results quickly. To illustrate the use of an open-source methodology in electromagnetics, we present two challenges. The first is to simulate data from a time domain airborne system over a conductive plate buried in a more resistive earth. The second is to jointly invert airborne TDEM and FDEM data with ground TDEM. SimPEG, Simulation and Parameter Estimation in Geophysics, (https://simpeg.xyz) is used for the open-source software. The figures in this paper can be reproduced by downloading the Jupyter Notebooks we provide with this paper (https://github.com/simpeg-research/oldenburg-2018-AEM). Access to the source code allows the researcher to explore the simulations and inversions by changing model and inversion parameters, plot fields and fluxes to gain further insight about the EM phenomena, and solve a new research problem by using open-source software as a base. By providing results in a manner that allows others to reproduce, further explore, and even extend them, we hope to demonstrate that an open-source paradigm has the potential to enable more rapid progress of the geophysics community as a whole.