title: Possible Production of Solar Spicules by Microfilament Eruptions

abstract: We examine Big Bear Solar Observatory (BBSO) Goode Solar Telescope (GST) high-spatial resolution (0''.06), high-cadence (3.45 s), H-alpha-0.8 Angstrom images of central-disk solar spicules, using data of Samanta et al. (2019). We compare with coronal-jet chromospheric-component observations of Sterling et al. (2010a). Morphologically, bursts of spicules, referred to as "enhanced spicular activities" by Samanta et al. (2019), appear as scaled-down versions of the jet's chromospheric component. Both the jet and the enhanced spicular activities appear as chromospheric-material strands, undergoing twisting-type motions of ~20---50 km/s in the jet and ~20---30 km/s in the enhanced spicular activities. Presumably, the jet resulted from a minifilament-carrying magnetic eruption. For two enhanced spicular activities that we examine in detail, we find tentative candidates for corresponding erupting microfilaments, but not expected corresponding base brightenings. Nonetheless, the enhanced-spicular-activities' interacting mixed-polarity base fields, frequent-apparent-twisting motions, and morphological similarities to the coronal jet's chromospheric-temperature component, suggest that erupting microfilaments might drive the enhanced spicular activities but be hard to detect, perhaps due to H-alpha opacity. Degrading the BBSO/GST-image resolution with a 1''.0-FWHM smoothing function yields enhanced spicular activities resembling the "classical spicules" described by, e.g., Beckers (1968). Thus, a microfilament eruption might be the fundamental driver of many spicules, just as a minifilament eruption is the fundamental driver of many coronal jets. Similarly, a 0".5-FWHM smoothing renders some enhanced spicular activities to resemble previously-reported "twinned" spicules, while the full-resolution features might account for spicules sometimes appearing as 2D-sheet-like structures.


title: A fast and effective MIP-based heuristic for a selective and periodic inventory routing problem in reverse logistics

abstract: We consider an NP-hard selective and periodic inventory routing problem (SPIRP) in a waste vegetable oil collection environment. This SPIRP arises in the context of reverse logistics where a biodiesel company has daily requirements of oil to be used as raw material in its production process. These requirements can be fulfilled by using the available inventory, collecting waste vegetable oil or purchasing virgin oil. The problem consists in determining a period (cyclic) planning for the collection and purchasing of oil such that the total collection, inventory and purchasing costs are minimized, while meeting the company's oil requirements and all the operational constraints. We propose a MIP-based heuristic which solves a relaxed model without routing, constructs routes taking into account the relaxation's solution and then improves these routes by solving the capacitated vehicle routing problem associated to each period. Following this approach, an a posteriori performance guarantee is ensured, as the approach provides both a lower bound and a feasible solution. The performed computational experiments show that the MIP-based heuristic is very fast and effective as it is able to encounter near optimal solutions with low gaps within seconds, improving several of the best known results using just a fraction of the time spent by a state-of-the-art heuristic. A remarkable fact is that the proposed MIP-based heuristic improves over the best known results for all the large instances available in the literature.


title: Disjoint type graphs with no short odd cycles

abstract: In this note, we provide a proof of a technical result of Erd\H{o}s and Hajnal about the existence of disjoint type graphs with no odd cycles. We also prove that this result is sharp in a certain sense.


title: State-dependent regulation of cortical processing speed via gain modulation

abstract: To thrive in dynamic environments, animals must be capable of rapidly and flexibly adapting behavioral responses to a changing context and internal state. Examples of behavioral flexibility include faster stimulus responses when attentive and slower responses when distracted. Contextual or state-dependent modulations may occur early in the cortical hierarchy and may be implemented via top-down projections from cortico-cortical or neuromodulatory pathways. However, the computational mechanisms mediating the effects of such projections are not known. Here, we introduce a theoretical framework to classify the effects of cell-type specific top-down perturbations on the information processing speed of cortical circuits. Our theory demonstrates that perturbation effects on stimulus processing can be predicted by intrinsic gain modulation, which controls the timescale of the circuit dynamics. Our theory leads to counter-intuitive effects such as improved performance with increased input variance. We tested the model predictions using large-scale electrophysiological recordings from the visual hierarchy in freely running mice, where we found that a decrease in single-cell intrinsic gain during locomotion led to an acceleration of visual processing. Our results establish a novel theory of cell-type specific perturbations, applicable to top-down modulation as well as optogenetic and pharmacological manipulations. Our theory links connectivity, dynamics, and information processing via gain modulation.


title: Infinitesimal Poisson Algebras and Linearization of Hamiltonian Systems

abstract: Using the notion of a contravariant derivative, we give some algebraic and geometric characterizations of Poisson algebras associated to the infinitesimal data of Poisson submanifolds. We show that such a class of Poisson algebras provides a suitable framework for the study of the Hamiltonization problem for the linearized dynamics along Poisson submanifolds.


title: The GeoLifeCLEF 2020 Dataset

abstract: Understanding the geographic distribution of species is a key concern in conservation. By pairing species occurrences with environmental features, researchers can model the relationship between an environment and the species which may be found there. To facilitate research in this area, we present the GeoLifeCLEF 2020 dataset, which consists of 1.9 million species observations paired with high-resolution remote sensing imagery, land cover data, and altitude, in addition to traditional low-resolution climate and soil variables. We also discuss the GeoLifeCLEF 2020 competition, which aims to use this dataset to advance the state-of-the-art in location-based species recommendation.


title: Convergence rates and approximation results for SGD and its continuous-time counterpart

abstract: This paper proposes a thorough theoretical analysis of Stochastic Gradient Descent (SGD) with non-increasing step sizes. First, we show that the recursion defining SGD can be provably approximated by solutions of a time inhomogeneous Stochastic Differential Equation (SDE) using an appropriate coupling. In the specific case of a batch noise we refine our results using recent advances in Stein's method. Then, motivated by recent analyses of deterministic and stochastic optimization methods by their continuous counterpart, we study the long-time behavior of the continuous processes at hand and establish non-asymptotic bounds. To that purpose, we develop new comparison techniques which are of independent interest. Adapting these techniques to the discrete setting, we show that the same results hold for the corresponding SGD sequences. In our analysis, we notably improve non-asymptotic bounds in the convex setting for SGD under weaker assumptions than the ones considered in previous works. Finally, we also establish finite-time convergence results under various conditions, including relaxations of the famous {\L}ojasiewicz inequality, which can be applied to a class of non-convex functions.


title: Stochastic quantization of Liouville conformal field theory

abstract: We study a nonlinear stochastic heat equation forced by a space-time white noise on closed surfaces, with nonlinearity $e^{\beta u}$. This equation corresponds to the stochastic quantization of the Liouville quantum gravity (LQG) measure. (i) We first revisit the construction of the LQG measure in Liouville conformal field theory (LCFT) in the $L^2$ regime $0<\beta<\sqrt{2}$. This uniformizes in this regime the approaches of David-Kupiainen-Rhodes-Vargas (2016), David-Rhodes-Vargas (2016) and Guillarmou-Rhodes-Vargas (2019) which treated the case of a closed surface with genus 0, 1 and $> 1$ respectively. Moreover, our argument shows that this measure is independent of the approximation procedure for a large class of smooth approximations. (ii) We prove almost sure global well-posedness of the parabolic stochastic dynamics, and invariance of the measure under this stochastic flow. In particular, our results improve previous results obtained by Garban (2020) in the cases of the sphere and the torus with their canonical metric, and are new in the case of closed surfaces with higher genus.


title: Thermodynamically consistent equation of state for an accreted neutron star crust

abstract: We study equation of state (EOS) of an accreting neutron star crust. Usually, such EOS is obtained assuming (implicitly) that the free (unbound) neutrons and nuclei in the inner crust move together. We argue, that this assumption violates the condition $\mu_n^\infty={\rm const}$, required for hydrostatic (and diffusion) equilibrium of unbound neutrons ($\mu^\infty_n$ is the redshifted neutron chemical potential). We construct a new EOS respecting this condition, working in the compressible liquid-drop approximation. We demonstrate that it is close to the catalyzed EOS in most part of the inner crust, being very different from EOSs of accreted crust discussed in the literature. In particular, the pressure at the outer-inner crust interface does not coincide with the neutron drip pressure, usually calculated in the literature, and is determined by hydrostatic (and diffusion) equilibrium conditions within the star. We also find an instability at the bottom of fully accreted crust that transforms nuclei into homogeneous nuclear matter. It guarantees that the structure of fully accreted crust remains self-similar during accretion.


title: Second-order cone representation for convex subsets of the plane

abstract: Semidefinite programming (SDP) is the task of optimizing a linear function over the common solution set of finitely many linear matrix inequalities (LMIs). For the running time of SDP solvers, the maximal matrix size of these LMIs is usually more critical than their number. The semidefinite extension degree $\text{sxdeg}(K)$ of a convex set $K\subseteq\mathbb R^n$ is the smallest number $d$ such that $K$ is a linear image of a finite intersection $S_1\cap\dots\cap S_N$, where each $S_i$ is a spectrahedron defined by a linear matrix inequality of size $\le d$. Thus $\text{sxdeg}(K)$ can be seen as a measure for the complexity of performing semidefinite programs over the set $K$. We give several equivalent characterizations of $\text{sxdeg}(K)$, and use them to prove our main result: $\text{sxdeg}(K)\le2$ holds for any closed convex semialgebraic set $K\subseteq\mathbb R^2$. In other words, such $K$ can be represented using the second-order cone.


title: Quantum Approximate Optimization of Non-Planar Graph Problems on a Planar Superconducting Processor

abstract: We demonstrate the application of the Google Sycamore superconducting qubit quantum processor to combinatorial optimization problems with the quantum approximate optimization algorithm (QAOA). Like past QAOA experiments, we study performance for problems defined on the (planar) connectivity graph of our hardware; however, we also apply the QAOA to the Sherrington-Kirkpatrick model and MaxCut, both high dimensional graph problems for which the QAOA requires significant compilation. Experimental scans of the QAOA energy landscape show good agreement with theory across even the largest instances studied (23 qubits) and we are able to perform variational optimization successfully. For problems defined on our hardware graph we obtain an approximation ratio that is independent of problem size and observe, for the first time, that performance increases with circuit depth. For problems requiring compilation, performance decreases with problem size but still provides an advantage over random guessing for circuits involving several thousand gates. This behavior highlights the challenge of using near-term quantum computers to optimize problems on graphs differing from hardware connectivity. As these graphs are more representative of real world instances, our results advocate for more emphasis on such problems in the developing tradition of using the QAOA as a holistic, device-level benchmark of quantum processors.


title: Bayesian Interpolants as Explanations for Neural Inferences

abstract: The notion of Craig interpolant, used as a form of explanation in automated reasoning, is adapted from logical inference to statistical inference and used to explain inferences made by neural networks. The method produces explanations that are at the same time concise, understandable and precise.


title: Transferable, Controllable, and Inconspicuous Adversarial Attacks on Person Re-identification With Deep Mis-Ranking

abstract: The success of DNNs has driven the extensive applications of person re-identification (ReID) into a new era. However, whether ReID inherits the vulnerability of DNNs remains unexplored. To examine the robustness of ReID systems is rather important because the insecurity of ReID systems may cause severe losses, e.g., the criminals may use the adversarial perturbations to cheat the CCTV systems. In this work, we examine the insecurity of current best-performing ReID models by proposing a learning-to-mis-rank formulation to perturb the ranking of the system output. As the cross-dataset transferability is crucial in the ReID domain, we also perform a back-box attack by developing a novel multi-stage network architecture that pyramids the features of different levels to extract general and transferable features for the adversarial perturbations. Our method can control the number of malicious pixels by using differentiable multi-shot sampling. To guarantee the inconspicuousness of the attack, we also propose a new perception loss to achieve better visual quality. Extensive experiments on four of the largest ReID benchmarks (i.e., Market1501 [45], CUHK03 [18], DukeMTMC [33], and MSMT17 [40]) not only show the effectiveness of our method, but also provides directions of the future improvement in the robustness of ReID systems. For example, the accuracy of one of the best-performing ReID systems drops sharply from 91.8% to 1.4% after being attacked by our method. Some attack results are shown in Fig. 1. The code is available at https://github.com/whj363636/Adversarial-attack-on-Person-ReID-With-Deep-Mis-Ranking.


title: Creating solitons with controllable and near zero velocity in Bose-Einstein condensates

abstract: Established techniques for deterministically creating dark solitons in repulsively interacting atomic Bose-Einstein condensates (BECs) can only access a narrow range of soliton velocities. Because velocity affects the stability of individual solitons and the properties of soliton-soliton interactions, this technical limitation has hindered experimental progress. Here we create dark solitons in highly anisotropic cigar-shaped BECs with arbitrary position and velocity by simultaneously engineering the amplitude and phase of the condensate wavefunction, improving upon previous techniques which only explicitly manipulated the condensate phase. The single dark soliton solution present in true 1D systems corresponds to the kink soliton in anisotropic 3D systems and is joined by a host of additional dark solitons including vortex ring and solitonic vortex solutions. We readily create dark solitons with speeds from zero to half the sound speed. The observed soliton oscillation frequency suggests that we imprinted solitonic vortices, which for our cigar-shaped system are the only stable solitons expected for these velocities. Our numerical simulations of 1D BECs show this technique to be equally effective for creating kink solitons when they are stable. We demonstrate the utility of this technique by deterministically colliding dark solitons with domain walls in two-component spinor BECs.


title: Finite Element Approximation of Finite Deformation Dislocation Mechanics

abstract: We develop and demonstrate the first general computational tool for finite deformation static and dynamic dislocation mechanics. A finite element formulation of finite deformation (Mesoscale) Field Dislocation Mechanics theory is presented. The model is a minimal enhancement of classical crystal/$J_2$ plasticity that fundamentally accounts for polar/excess dislocations at the mesoscale. It has the ability to compute the static and dynamic finite deformation stress fields of arbitrary (evolving) dislocation distributions in finite bodies of arbitrary shape and elastic anisotropy under general boundary conditions. This capability is used to present a comparison of the static stress fields, at finite and small deformations, for screw and edge dislocations, revealing heretofore unexpected differences. The computational framework is verified against the sharply contrasting predictions of geometrically linear and nonlinear theories for the stress field of a spatially homogeneous dislocation distribution in the body, as well as against other exact results of the theory. Verification tests of the time-dependent numerics are also presented. Size effects in crystal and isotropic versions of the theory are shown to be a natural consequence of the model and are validated against available experimental data. With inertial effects incorporated, the development of an (asymmetric) propagating Mach cone is demonstrated in the finite deformation theory when a dislocation moves at speeds greater than the linear elastic shear wave speed of the material.


title: The ELM Survey South. I. An Effective Search for Extremely Low Mass White Dwarfs

abstract: We begin the search for extremely-low mass ($M\leq0.3M_{\odot}$, ELM) white dwarfs (WDs) in the southern sky based on photometry from the VST ATLAS and SkyMapper surveys. We use a similar color-selection method as the Hypervelocity star survey. We switched to an astrometric selection once Gaia Data Release 2 became available. We use the previously known sample of ELM white dwarfs to demonstrate that these objects occupy a unique parameter space in parallax and magnitude. We use the SOAR 4.1m telescope to test the Gaia-based selection, and identify more than two dozen low-mass white dwarfs, including 6 new ELM white dwarf binaries with periods as short as 2 h. The better efficiency of the Gaia-based selection enables us to extend the ELM Survey footprint to the southern sky. We confirm one of our candidates, J0500$-$0930, to become the brightest ($G=12.6$ mag) and closest ($d=72$ pc) ELM white dwarf binary currently known. Remarkably, the Transiting Exoplanet Survey Satellite (TESS) full-frame imaging data on this system reveals low-level ($<0.1$%) but significant variability at the orbital period of this system ($P=9.5$ h), likely from the relativistic beaming effect. TESS data on another system, J0642$-$5605, reveals ellipsoidal variations due to a tidally distorted ELM WD. These demonstrate the power of TESS full-frame images in confirming the orbital periods of relatively bright compact object binaries.


title: On recursively defined combinatorial classes and labelled trees

abstract: We define and prove isomorphisms between three combinatorial classes involving labeled trees. We also give an alternative proof by means of generating functions.


title: Embedding spheres in knot traces

abstract: The trace of $n$-framed surgery on a knot in $S^3$ is a 4-manifold homotopy equivalent to the 2-sphere. We characterise when a generator of the second homotopy group of such a manifold can be realised by a locally flat embedded 2-sphere whose complement has abelian fundamental group. Our characterisation is in terms of classical and computable 3-dimensional knot invariants. For each $n$, this provides conditions that imply a knot is topologically $n$-shake slice, directly analogous to the result of Freedman and Quinn that a knot with trivial Alexander polynomial is topologically slice.


title: Optimal Orbital Selection for Full Configuration Interaction (OptOrbFCI): Pursuing the Basis Set Limit under a Budget

abstract: Full configuration interaction (FCI) solvers are limited to small basis sets due to their expensive computational costs. An optimal orbital selection for FCI (OptOrbFCI) is proposed to boost the power of existing FCI solvers to pursue the basis set limit under a computational budget. The optimization problem coincides with that of the complete active space SCF method (CASSCF), while OptOrbFCI is algorithmically quite different. OptOrbFCI effectively finds an optimal rotation matrix via solving a constrained optimization problem directly to compress the orbitals of large basis sets to one with a manageable size, conducts FCI calculations only on rotated orbital sets, and produces a variational ground-state energy and its wave function. Coupled with coordinate descent full configuration interaction (CDFCI), we demonstrate the efficiency and accuracy of the method on the carbon dimer and nitrogen dimer under basis sets up to cc-pV5Z. We also benchmark the binding curve of the nitrogen dimer under the cc-pVQZ basis set with 28 selected orbitals, which provide consistently lower ground-state energies than the FCI results under the cc-pVDZ basis set. The dissociation energy in this case is found to be of higher accuracy.


title: C++11/14 Mutation Operators Based on Common Fault Patterns

abstract: The C++11/14 standard offers a wealth of features aimed at helping programmers write better code. Unfortunately, some of these features may cause subtle programming faults, likely to go unnoticed during code reviews. In this paper we propose four new mutation operators for C++11/14 based on common fault patterns, which allow to verify whether a unit test suite is capable of testing against such faults. We validate the relevance of the proposed mutation operators by performing a case study on seven real-life software systems.


title: Computing Convex Partitions for Point Sets in the Plane: The CG:SHOP Challenge 2020

abstract: We give an overview of the 2020 Computational Geometry Challenge, which targeted the problem of partitioning the convex hull of a given planar point set P into the smallest number of convex faces, such that no point of P is contained in the interior of a face.


title: How Crisp is the Crease? A Subjective Study on Web Browsing Perception of Above-The-Fold

abstract: Quality of Experience (QoE) for various types of websites has gained significant attention in recent years. In order to design and evaluate websites, a metric that can estimate a user's experienced quality robustly for diverse content is necessary. SpeedIndex (SI) has been widely adopted to estimate perceived web page loading progress. It measures the speed of rendering pixels for the webpage that is visible in the browser window. This is termed Above-The-Fold (ATF). The influence of animated content on the perception of ATF has been less comprehensively explored. In this paper, we present an experimental design and methodology to measure ATF perception for websites with and without animated elements for various page content categories. We found that pages with animated elements caused people to have more varied perceptions of ATF under different network conditions. Animated content also impacts the page load estimation accuracy of SI for websites. We discuss how the difference in the perception of ATF will impact the QoE management of web applications. We explain the necessity of revisiting the visual assessment of ATF to include the animated contents and improve the robustness of metrics like SI.


title: A single image deep learning approach to restoration of corrupted remote sensing products

abstract: Remote sensing images are used for a variety of analyses, from agricultural monitoring, to disaster relief, to resource planning, among others. The images can be corrupted due to a number of reasons, including instrument errors and natural obstacles such as clouds. We present here a novel approach for reconstruction of missing information in such cases using only the corrupted image as the input. The Deep Image Prior methodology eliminates the need for a pre-trained network or an image database. It is shown that the approach easily beats the performance of traditional single-image methods.


title: Studying the Landau mass parameter of the extended $\sigma$-$\omega$ model for neutron star matter

abstract: We present a Bayesian analysis of the Landau mass within the extended $\sigma$-$\omega$ model for neutron star matter. To this purpose, we consider the mass measurement of the object PSR 0740+6620, the tidal deformability estimation from the GW170817 and the mass-radius estimate of PSR J0030+0451 by NICER. Using Landau mass as free parameter of the theory, we rely on the prediction power of the Bayesian method to find the best value for this nuclear quantity.


title: Do Null-Type Mutation Operators Help Prevent Null-Type Faults?

abstract: The null-type is a major source of faults in Java programs, and its overuse has a severe impact on software maintenance. Unfortunately traditional mutation testing operators do not cover null-type faults by default, hence cannot be used as a preventive measure. We address this problem by designing four new mutation operators which model null-type faults explicitly. We show how these mutation operators are capable of revealing the missing tests, and we demonstrate that these mutation operators are useful in practice. For the latter, we analyze the test suites of 15 open-source projects to describe the trade-offs related to the adoption of these operators to strengthen the test suite.


title: A characterization of singular Schr\"odinger operators on the half-line

abstract: We study a class of delta-like perturbations of the Laplacian on the half-line, characterized by Robin boundary conditions at the origin. Using the formalism of nonstandard analysis, we derive a simple connection with a suitable family of Schr\"{o}dinger operators with potentials of very large (infinite) magnitude and very short (infinitesimal) range. As a consequence, we also derive a similar result for point interactions in the Euclidean space $\mathbb{R}^3$, in the case of radial potentials. Moreover, we discuss explicitly our results in the case of potentials that are linear in a neighbourhood of the origin.


title: $\Lambda_c^{\pm}$ production in pp collisions with a new fragmentation function

abstract: We study inclusive $\Lambda_c^{\pm}$-baryon production in $pp$ collisions in the general-mass variable-flavor number scheme and compare with data from the LHCb, ALICE and CMS collaborations. We perform a new fit of the $c \to \Lambda_c^+$ fragmentation function combining $e^+e^-$ data from OPAL and Belle. The agreement with LHC data is slightly worse compared with a calculation using an older fragmentation function, and the tension between different determinations of $\Lambda_c^{\pm}$ production cross sections from the LHC experimental collaborations is not resolved. The ratio of data for $\Lambda_c^+$-baryon and $D^0$-meson production seems to violate the universality of $c$-charm quark to $c$-hadron fragmentation.


title: Optimal Runtime Verification of Finite State Properties over Lossy Event Streams

abstract: Monitoring programs for finite state properties is challenging due to high memory and execution time overheads it incurs. Some events if skipped or lost naturally can reduce both overheads, but lead to uncertainty about the current monitor state. In this work, we present a theoretical framework to model these lossy event streams and provide a construction for a monitor which observes them without producing false positives. The constructed monitor is optimally sound among all complete monitors. We model several loss types of practical relevance using our framework and provide construction of smaller approximate monitors for properties with a large number of states.


title: Generating functions for a lattice path model introduced by Deutsch

abstract: The lattice path model suggested by E. Deutsch is derived from ordinary Dyck paths, but with additional down-steps of size -3,-5,-7,... . For such paths, we find the generating functions of them, according to length, ending at level $i$, both, when considering them from left to right and from right to left. The generating functions are intrinsically cubic, and thus (for $i=0$) in bijection to various objects, like even trees, ternary trees, etc.


title: Generating Counter Narratives against Online Hate Speech: Data and Strategies

abstract: Recently research has started focusing on avoiding undesired effects that come with content moderation, such as censorship and overblocking, when dealing with hatred online. The core idea is to directly intervene in the discussion with textual responses that are meant to counter the hate content and prevent it from further spreading. Accordingly, automation strategies, such as natural language generation, are beginning to be investigated. Still, they suffer from the lack of sufficient amount of quality data and tend to produce generic/repetitive responses. Being aware of the aforementioned limitations, we present a study on how to collect responses to hate effectively, employing large scale unsupervised language models such as GPT-2 for the generation of silver data, and the best annotation strategies/neural architectures that can be used for data filtering before expert validation/post-editing.


title: Exact solution and the multidimensional Godunov scheme for the acoustic equations

abstract: The acoustic equations derived as a linearization of the Euler equations are a valuable system for studies of multi-dimensional solutions. Additionally they possess a low Mach number limit analogous to that of the Euler equations. Aiming at understanding the behaviour of the multi-dimensional Godunov scheme in this limit, first the exact solution of the corresponding Cauchy problem in three spatial dimensions is derived. The appearance of logarithmic singularities in the exact solution of the 4-quadrant Riemann Problem in two dimensions is discussed. The solution formulae are then used to obtain the multidimensional Godunov finite volume scheme in two dimensions. It is shown to be superior to the dimensionally split upwind/Roe scheme concerning its domain of stability and ability to resolve multi-dimensional Riemann problems. It is shown experimentally and theoretically that despite taking into account multi-dimensional information it is, however, not able to resolve the low Mach number limit.


title: Spin dynamics and a nearly continuous magnetic phase transition in an entropy-stabilized oxide antiferromagnet

abstract: The magnetic order and the spin dynamics in the antiferromagnetic entropy-stabilized oxide (Mg$_{0.2}$Co$_{0.2}$Ni$_{0.2}$Cu$_{0.2}$Zn$_{0.2}$)O (MgO-ESO) have been studied using muon spin relaxation ($\mu$SR) and inelastic neutron scattering. We find that antiferromagnetic order develops gradually in the sample volume as it is cooled below 140 K, becoming fully ordered around 100 K. The spin dynamics show a critical slowing down in the vicinity of the transition, and the magnetic order parameter grows continuously in the ordered state. These results indicate that the antiferromagnetic transition is continuous but proceeds with a Gaussian distribution of ordering temperatures. The magnetic contribution to the specific heat determined from inelastic neutron scattering likewise shows a broad feature centered around 120 K. High-resolution inelastic neutron scattering further reveals an initially gapped spectrum at low temperature which sees an increase in a quasielastic contribution upon heating until the ordering temperature.


title: Dehn fillings of knot manifolds containing essential twice-punctured tori

abstract: We show that if a hyperbolic knot manifold $M$ contains an essential twice-punctured torus $F$ with boundary slope $\beta$ and admits a filling with slope $\alpha$ producing a Seifert fibred space, then the distance between the slopes $\alpha$ and $\beta$ is less than or equal to $5$ unless $M$ is the exterior of the figure eight knot. The result is sharp; the bound of $5$ can be realized on infinitely many hyperbolic knot manifolds. We also determine distance bounds in the case that the fundamental group of the $\alpha$-filling contains no non-abelian free group. The proofs are divided into the four cases $F$ is a semi-fibre, $F$ is a fibre, $F$ is non-separating but not a fibre, and $F$ is separating but not a semi-fibre, and we obtain refined bounds in each case.


title: Determination of spatial configuration of an underwater swarm with minimum data

abstract: The subject is the localization problem of an underwater swarm of autonomous underwater robots (AUV), in the frame of the HARNESS project; by localization, we mean the relative swarm configuration, i.e., the geometrical shape of the group. The result is achieved by using the signals that the robots exchange. The swarm is organized by rules and conceived to perform tasks, ranging from environmental monitoring to terrorism-attack surveillance. Two methods of determining the shape of the swarm, both based on trilateration calculation, are proposed. The first method focuses on the robot's speed. In this case, we use our knowledge of the speeds and distances between the machines, while the second method considers only distances and the orientation angles of the robots. Unlike a trilateration problem, we do not know the position of the beacons and this renders the problem a difficult one. Moreover, we have very few data. More than one step of motion is needed to resolve the multiple solutions found, owing to the symmetries of the system and optimization process of one or more objective functions leading to the final configuration. We subsequently checked our algorithm using a simulator taking into account random errors affecting the measurements


title: Saliency-based Weighted Multi-label Linear Discriminant Analysis

abstract: In this paper, we propose a new variant of Linear Discriminant Analysis (LDA) to solve multi-label classification tasks. The proposed method is based on a probabilistic model for defining the weights of individual samples in a weighted multi-label LDA approach. Linear Discriminant Analysis is a classical statistical machine learning method, which aims to find a linear data transformation increasing class discrimination in an optimal discriminant subspace. Traditional LDA sets assumptions related to Gaussian class distributions and single-label data annotations. To employ the LDA technique in multi-label classification problems, we exploit intuitions coming from a probabilistic interpretation of class saliency to redefine the between-class and within-class scatter matrices. The saliency-based weights obtained based on various kinds of affinity encoding prior information are used to reveal the probability of each instance to be salient for each of its classes in the multi-label problem at hand. The proposed Saliency-based weighted Multi-label LDA approach is shown to lead to performance improvements in various multi-label classification problems.


title: Centralized and decentralized isolation strategies and their impact on the COVID-19 pandemic dynamics

abstract: The infectious diseases are spreading due to human interactions enabled by various social networks. Therefore, when a new pathogen such as SARS-CoV-2 causes an outbreak, the non-pharmaceutical isolation strategies (e.g., social distancing) are the only possible response to disrupt its spreading. To this end, we introduce the new epidemic model (SICARS) and compare the centralized (C), decentralized (D), and combined (C+D) social distancing strategies, and analyze their efficiency to control the dynamics of COVID-19 on heterogeneous complex networks. Our analysis shows that the centralized social distancing is necessary to minimize the pandemic spreading. The decentralized strategy is insufficient when used alone, but offers the best results when combined with the centralized one. Indeed, the (C+D) is the most efficient isolation strategy at mitigating the network superspreaders and reducing the highest node degrees to less than 10% of their initial values. Our results also indicate that stronger social distancing, e.g., cutting 75% of social ties, can reduce the outbreak by 75% for the C isolation, by 33% for the D isolation, and by 87% for the (C+D) isolation strategy. Finally, we study the impact of proactive versus reactive isolation strategies, as well as their delayed enforcement. We find that the reactive response to the pandemic is less efficient, and delaying the adoption of isolation measures by over one month (since the outbreak onset in a region) can have alarming effects; thus, our study contributes to an understanding of the COVID-19 pandemic both in space and time. We believe our investigations have a high social relevance as they provide insights into understanding how different degrees of social distancing can reduce the peak infection ratio substantially; this can make the COVID-19 pandemic easier to understand and control over an extended period of time.


title: Particle Trapping in Axisymmetric Electron Holes

abstract: Electron orbits are calculated in solitary two-dimensional axisymmetric electrostatic potential structures, typical of plasma electron holes, in order to establish the conditions for the particles to remain trapped. Analytic calculations of the evolution of the parallel energy caused by the perturbing radial electric field are shown to agree well with full numerical orbit integration Poincar\'e plots. The predominant mechanism of detrapping is resonance between the gyro frequency in the parallel magnetic field and harmonics of the parallel bounce frequency. A region of phase-space adjacent to the trapped-passing boundary in parallel energy is generally stochastic because of island overlap of different harmonics, but except for very strong radial electric field perturbation, more deeply trapped orbits have well-defined islands and are permanently confined. A simple universal quantitative algorithm is given, and its results plotted as a function of magnetic field strength and hole radial scale-length, determining the phase space volume available to sustain the electron hole by depression of the permanently trapped distribution function.


title: An effective version of Belyi's theorem in positive characteristic

abstract: The purpose of the present paper is to give an effective version of the noncritical $p$-tame Belyi theorem. That is to say, we compute explicitly an upper bound of the minimal degree of tamely ramified Belyi maps in positive characteristic which are unramified at a prescribed finite set of points.


title: Measuring Emotions in the COVID-19 Real World Worry Dataset

abstract: The COVID-19 pandemic is having a dramatic impact on societies and economies around the world. With various measures of lockdowns and social distancing in place, it becomes important to understand emotional responses on a large scale. In this paper, we present the first ground truth dataset of emotional responses to COVID-19. We asked participants to indicate their emotions and express these in text. This resulted in the Real World Worry Dataset of 5,000 texts (2,500 short + 2,500 long texts). Our analyses suggest that emotional responses correlated with linguistic measures. Topic modeling further revealed that people in the UK worry about their family and the economic situation. Tweet-sized texts functioned as a call for solidarity, while longer texts shed light on worries and concerns. Using predictive modeling approaches, we were able to approximate the emotional responses of participants from text within 14% of their actual value. We encourage others to use the dataset and improve how we can use automated methods to learn about emotional responses and worries about an urgent problem.


title: Schmidt rank constraints in Quantum Information Theory

abstract: Can vectors with low Schmidt rank form mutually unbiased bases? Can vectors with high Schmidt rank form positive under partial transpose states? In this work, we address these questions by presenting several new results related to Schmidt rank constraints and their compatibility with other properties. We provide an upper bound on the number of mutually unbiased bases of $\mathbb{C}^m\otimes\mathbb{C}^n$ $(m\leq n)$ formed by vectors with low Schmidt rank. In particular, the number of mutually unbiased product bases of $\mathbb{C}^m\otimes\mathbb{C}^n$ cannot exceed $m+1$, which solves a conjecture proposed by McNulty et al. Then we show how to create a positive under partial transpose entangled state from any state supported on the antisymmetric space and how their Schmidt numbers are exactly related. Finally, we show that the Schmidt number of operator Schmidt rank 3 states of $\mathcal{M}_m\otimes \mathcal{M}_n\ (m\leq n)$ that are invariant under left partial transpose cannot exceed $m-2$.


title: Formal Test Synthesis for Safety-Critical Autonomous Systems based on Control Barrier Functions

abstract: The prolific rise in autonomous systems has led to questions regarding their safe instantiation in real-world scenarios. Failures in safety-critical contexts such as human-robot interactions or even autonomous driving can ultimately lead to loss of life. In this context, this paper aims to provide a method by which one can algorithmically test and evaluate an autonomous system. Given a black-box autonomous system with some operational specifications, we construct a minimax problem based on control barrier functions to generate a family of test parameters designed to optimally evaluate whether the system can satisfy the specifications. To illustrate our results, we utilize the Robotarium as a case study for an autonomous system that claims to satisfy waypoint navigation and obstacle avoidance simultaneously. We demonstrate that the proposed test synthesis framework systematically finds those sequences of events (tests) that identify points of system failure.


title: Asking and Answering Questions to Evaluate the Factual Consistency of Summaries

abstract: Practical applications of abstractive summarization models are limited by frequent factual inconsistencies with respect to their input. Existing automatic evaluation metrics for summarization are largely insensitive to such errors. We propose an automatic evaluation protocol called QAGS (pronounced "kags") that is designed to identify factual inconsistencies in a generated summary. QAGS is based on the intuition that if we ask questions about a summary and its source, we will receive similar answers if the summary is factually consistent with the source. To evaluate QAGS, we collect human judgments of factual consistency on model-generated summaries for the CNN/DailyMail (Hermann et al., 2015) and XSUM (Narayan et al., 2018) summarization datasets. QAGS has substantially higher correlations with these judgments than other automatic evaluation metrics. Also, QAGS offers a natural form of interpretability: The answers and questions generated while computing QAGS indicate which tokens of a summary are inconsistent and why. We believe QAGS is a promising tool in automatically generating usable and factually consistent text.


title: Integrative dynamic structural biology unveils conformers essential for the oligomerization of a large GTPase

abstract: Guanylate binding proteins (GBPs) are soluble dynamin-like proteins with structured domains that undergo a conformational transition for GTP-controlled oligomerization to exert their function as part of the innate immune system of mammalian cells - attacking intra-cellular parasites by disrupting their membranes. The structural basis and mechanism of this process is unknown. Therefore, we apply neutron spin echo, X-ray scattering, fluorescence, and EPR spectroscopy as techniques for integrative dynamic structural biology to human GBP1 (hGBP1). We mapped hGBP1's essential dynamics from nanoseconds to milliseconds by motional spectra of sub-domains. We find a GTP-independent flexibility of the C-terminal effector domain in the $\mu$s-regime and structurally characterize conformers being essential that hGBP1 can open like a pocketknife for oligomerization. This unveils the intrinsic flexibility, a GTP-triggered association of the GTPase-domains and assembly-dependent GTP-hydrolysis as functional design principles of hGBP1 that control its reversible oligomerization in polar assemblies and the subsequent formation of condensates.


title: Impossible moons -- Transit timing effects that cannot be due to an exomoon

abstract: Exomoons are predicted to produce transit timing variations (TTVs) upon their host planet. Unfortunately, so are many other astrophysical phenomena - most notably other planets in the system. In this work, an argument of reductio ad absurdum is invoked, by deriving the transit timing effects that are impossible for a single exomoon to produce. Our work derives three key analytic tests. First, one may exploit the fact that a TTV signal from an exomoon should be accompanied by transit duration variations (TDVs), and that one can derive a TDV floor as a minimum expected level of variability. Cases for which the TDV upper limit is below this floor can thus be killed as exomoon candidates. Second, formulae are provided for estimating whether moons are expected to be 'killable' when no TDVs presently exist, thus enabling the community to estimate whether it's even worth deriving TDVs in the first place. Third, a TTV ceiling is derived, above which exomoons should never be able to produce TTV amplitudes. These tools are applied to a catalog of TTVs and TDVs for two and half thousand Kepler Objects Interest, revealing over two hundred cases that cannot be due to a moon. These tests are also applied to the exomoon candidate Kepler-1625b i, which comfortably passes the criteria. These simple analytic results should provide a means of rapidly rejecting putative exomoons and streamlining the search for satellites.


title: Stars at infinity in Teichm\"uller space

abstract: We investigate a metric structure on the Thurston boundary of Teichm\"uller space. To do this, we develop tools in sup metrics and apply Minsky's theorem.


title: Skew braces of size $p^2q$ II: non-abelian type

abstract: In this paper we enumerate the skew braces of non-abelian type of size $p^2q$ for $p,q$ primes with $q>2$ by the classification of regular subgroups of the holomorph of the non-abelian groups of the same order. Since Crespo dealt with the case $q=2$, this paper completes the enumeration of skew braces of size $p^2q$ started in a previous work by the authors. In some cases, we provide also a structural description of the skew braces. As an application, we prove a conjecture posed by V. Bardakov, M. Neshchadim and M. Yadav.


title: Symplectic Rigidity of Fibers in Cotangent Bundles of Open Riemann Surfaces

abstract: We study symplectic rigidity phenomena for fibers in cotangent bundles of Riemann surfaces. Our main result can be seen as a generalization to open Riemann surfaces of arbitrary genus of work of Eliashberg and Polterovich on the Nearby Lagrangian Conjecture for $T^* \mathbb{R}^2$. As a corollary, we answer a strong version in dimension $2n=4$ of a question of Eliashberg about linking of Lagrangian disks in $T^* \mathbb{R}^n$, which was previously answered by Ekholm and Smith in dimensions $2n \geq 8$.


title: Dispersion relations applied to double-folding potentials from chiral EFT

abstract: We present a determination of optical potentials using the double-folding method based on chiral effective field theory nucleon-nucleon interactions at next-to-next-to-leading order combined with dispersion relations to constrain the imaginary part. This approach is benchmarked on $^{16}$O--$^{16}$O collisions, and extended to the $^{12}$C--$^{12}$C and $^{12}$C--$^{16}$O cases. Predictions derived from these potentials are compared to data for elastic scattering at energies up to 1000 MeV, as well as for fusion at low energy. Without adjusting parameters, excellent agreement with experiment is found. In addition, we study the sensitivity of the corresponding cross sections to the nucleon-nucleon interactions and nuclear densities used.


title: Application of Jordan Decomposition to Non-Hermitian Lattice Models with Spectrally-Isolated Lower Dimensional States

abstract: When analyzing non-Hermitian lattice systems, the standard eigenmode decomposition utilized for the analysis of Hermitian systems must be replaced by Jordan decomposition. This approach enables us to identify the correct number of the left and right eigenstates of a large finite-sized lattice system, and to form a complete basis for calculating the resonant excitation of the system. Specifically, we derive the procedure for applying Jordan decomposition to a system with spectrally-isolated states. We use a non-Hermitian quadrupole insulator with zero-energy corner states as an example of a large system whose dimensionality can be drastically reduced to derive a low-dimensional "defective" Hamiltonian describing such localized states. Counter-intuitive and non-local properties of the resonant response of the system near zero energy are explained using the Jordan decomposition approach. Depending on the excitation properties of the corner states, we classify our non-Hermitian quadrupolar insulator into three categories: trivial, near-Hermitian, and non-local.

categories: cond-mat.mes-hall cond-mat.quant-gas cond-mat.str-el physics.optics


title: Rotation of Low-Mass Stars in Taurus with K2

abstract: We present an analysis of K2 light curves (LCs) from Campaigns 4 and 13 for members of the young ($\sim$3 Myr) Taurus association, in addition to an older ($\sim$30 Myr) population of stars that is largely in the foreground of the Taurus molecular clouds. Out of 156 of the highest-confidence Taurus members, we find that 81\% are periodic. Our sample of young foreground stars is biased and incomplete, but nearly all (37/38) are periodic. The overall distribution of rotation rates as a function of color (a proxy for mass) is similar to that found in other clusters: the slowest rotators are among the early M spectral types, with faster rotation towards both earlier FGK and later M types. The relationship between period and color/mass exhibited by older clusters such as the Pleiades is already in place by Taurus age. The foreground population has very few stars, but is consistent with the USco and Pleiades period distributions. As found in other young clusters, stars with disks rotate on average slower, and few with disks are found rotating faster than $\sim$2 d. The overall amplitude of the light curves decreases with age and higher mass stars have generally lower amplitudes than lower mass stars. Stars with disks have on average larger amplitudes than stars without disks, though the physical mechanisms driving the variability and the resulting light curve morphologies are also different between these two classes.


title: Hitting the Ground Running: Computational physics education to prepare students for computational physics research

abstract: Momentum exists in the physics community for integrating computation into the undergraduate curriculum. One of many benefits would be preparation for computational research. Our investigation poses the question of which computational skills might be best learned in the curriculum (prior to research) versus during research. Based on a survey of computational physicists, we present evidence that many relevant skills are developed naturally in a research context while others stand out as best learned in advance.


title: Control and Trajectory Optimization for Soft Aerial Manipulation

abstract: Manipulation and grasping with unmanned aerial vehicles (UAVs) currently require accurate positioning and are often executed at reduced speed to ensure successful grasps. This is due to the fact that typical UAVs can only accommodate rigid manipulators with few degrees of freedom, which limits their capability to compensate for disturbances caused by the vehicle positioning errors. Moreover, UAVs have to minimize external contact forces in order to maintain stability. Biological systems, on the other hand, exploit softness to overcome similar limitations, and leverage compliance to enable aggressive grasping. This paper investigates control and trajectory optimization for a soft aerial manipulator, consisting of a quadrotor and a tendon-actuated soft gripper, in which the advantages of softness can be fully exploited. To the best of our knowledge, this is the first work at the intersection between soft manipulation and UAV control. We present a decoupled approach for the quadrotor and the soft gripper, combining (i) a geometric controller and a minimum-snap trajectory optimization for the quadrotor (rigid) base, with (ii) a quasi-static finite element model and control-space interpolation for the soft gripper. We prove that the geometric controller asymptotically stabilizes the quadrotor velocity and attitude despite the addition of the soft load. Finally, we evaluate the proposed system in a realistic soft dynamics simulator, and show that: (i) the geometric controller is fairly insensitive to the soft payload, (ii) the platform can reliably grasp unknown objects despite inaccurate positioning and initial conditions, and (iii) the decoupled controller is amenable for real-time execution.


title: On the universality of the KRK factorization scheme

abstract: The factorisation scheme (FS) abbreviated as KRK FS including a new definition of the PDFs for initial hadrons was formulated while developing KrkNLO scheme of matching QCD NLO corrections for the hard process with the parton shower heavy boson production in hadron-hadron collision and for deep inelastic lepton-hadron scattering. KRK FS (originally called Monte Carlo FS) can be regarded as a variant of the $\overline{MS}$ system. It is therefore trivially universal, that is process independent. The question of its universality is formulated differently: As the basic role of KRK FS is to simplify drastically NLO corrections, the question is now whether the same {\em single} variant of PDFs in the KRK FS is able to achieve the same maximal simplification of the NLO corrections for {\em all processes} with one or two initial hadrons and any number of the final hadrons? Our answer is positive and the proof is elaborated in the present note within the Catani-Seymour subtraction methodology. KRK FS is mandatory in the KrkNLO method of matching NLO calculation and parton shower -- a much simpler alternative of POWHEG and/or MC\@ NLO. However, the use of KRK FS and the corresponding PDFs simplifies NLO calculations for any other method of calculating NLO corrections and for arbitrary processes as well.


title: Higgs self-coupling measurements using deep learning in the $b\bar{b}b\bar{b}$ final state

abstract: Measuring the Higgs trilinear self-coupling $\lambda_{hhh}$ is experimentally demanding but fundamental for understanding the shape of the Higgs potential. We present a comprehensive analysis strategy for the HL-LHC using di-Higgs events in the four $b$-quark channel ($hh \to 4b$), extending current methods in several directions. We perform deep learning to suppress the formidable multijet background with dedicated optimisation for BSM $\lambda_{hhh}$ scenarios. We compare the $\lambda_{hhh}$ constraining power of events using different multiplicities of large radius jets with a two-prong structure that reconstruct boosted $h \to bb$ decays. We show that current uncertainties in the SM top Yukawa coupling $y_t$ can modify $\lambda_{hhh}$ constraints by $\sim 20\%$. For SM $y_t$, we find prospects of $-0.8 < \lambda_{hhh} / \lambda_{hhh}^\text{SM} < 6.6$ at 68% CL under simplified assumptions for 3000~fb$^{-1}$ of HL-LHC data. Our results provide a careful assessment of di-Higgs identification and machine learning techniques for all-hadronic measurements of the Higgs self-coupling and sharpens the requirements for future improvement.


title: Recent developments in determining the evolution of magnetic OB stars

abstract: We review recent developments in determining the evolution of magnetic massive OB stars. One of the important scientific questions is the completeness and the detection limits of contemporaneous spectropolarimetric surveys across the HR diagram. We present the characteristics of the MiMeS survey of O-type stars, and how the limits of the current available observations warrant the design of new high-precision surveys that target the older O-type star population. Another important question is whether the presence of the magnetic fields changes stellar evolutionary tracks in a significant way, hence leading to a wrong determination of stellar parameters. We review new evolution models that include the effect of magnetic wind quenching, and suppression of convection in the iron opacity peak zone.


title: Deep Manifold Prior

abstract: We present a prior for manifold structured data, such as surfaces of 3D shapes, where deep neural networks are adopted to reconstruct a target shape using gradient descent starting from a random initialization. We show that surfaces generated this way are smooth, with limiting behavior characterized by Gaussian processes, and we mathematically derive such properties for fully-connected as well as convolutional networks. We demonstrate our method in a variety of manifold reconstruction applications, such as point cloud denoising and interpolation, achieving considerably better results against competitive baselines while requiring no training data. We also show that when training data is available, our method allows developing alternate parametrizations of surfaces under the framework of AtlasNet, leading to a compact network architecture and better reconstruction results on standard image to shape reconstruction benchmarks.


title: Error correction and extraction in request dialogs

abstract: We propose a dialog system utility component that gets the last two utterances of a user and can detect whether the last utterance is an error correction of the second last utterance. If yes, it corrects the second last utterance according to the error correction in the last utterance and outputs the extracted pairs of reparandum and repair entity. This component offers two advantages, learning the concept of corrections to avoid collecting corrections for every new domain and extracting reparandum and repair pairs, which offers the possibility to learn out of it. For the error correction one sequence labeling and two sequence to sequence approaches are presented. For the error correction detection these three error correction approaches can also be used and in addition, we present a sequence classification approach. One error correction detection and one error correction approach can be combined to a pipeline or the error correction approaches can be trained and used end-to-end to avoid two components. We modified the EPIC-KITCHENS-100 dataset to evaluate the approaches for correcting entity phrases in request dialogs. For error correction detection and correction, we got an accuracy of 96.40 % on synthetic validation data and an accuracy of 77.81 % on human-created real-world test data.


title: Variable Rate Video Compression using a Hybrid Recurrent Convolutional Learning Framework

abstract: In recent years, neural network-based image compression techniques have been able to outperform traditional codecs and have opened the gates for the development of learning-based video codecs. However, to take advantage of the high temporal correlation in videos, more sophisticated architectures need to be employed. This paper presents PredEncoder, a hybrid video compression framework based on the concept of predictive auto-encoding that models the temporal correlations between consecutive video frames using a prediction network which is then combined with a progressive encoder network to exploit the spatial redundancies. A variable-rate block encoding scheme has been proposed in the paper that leads to remarkably high quality to bit-rate ratios. By joint training and fine-tuning of this hybrid architecture, PredEncoder has been able to gain significant improvement over the MPEG-4 codec and has achieved bit-rate savings over the H.264 codec in the low to medium bit-rate range for HD videos and comparable results over most bit-rates for non-HD videos. This paper serves to demonstrate how neural architectures can be leveraged to perform at par with the highly optimized traditional methodologies in the video compression domain.


title: Folding of Hitchin systems and crepant resolutions

abstract: Folding of ADE-Dynkin diagrams according to graph automorphisms yields irreducible Dynkin diagrams of ABCDEFG-types. This folding procedure allows to trace back the properties of the corresponding simple Lie algebras or groups to those of ADE-type. In this article, we implement the techniques of folding by graph automorphisms for Hitchin integrable systems. We show that the fixed point loci of these automorphisms are isomorphic as algebraic integrable systems to the Hitchin systems of the folded groups away from singular fibers. The latter Hitchin systems are isomorphic to the intermediate Jacobian fibrations of Calabi--Yau orbifold stacks constructed by the first author. We construct simultaneous crepant resolutions of the associated singular quasi-projective Calabi--Yau threefolds and compare the resulting intermediate Jacobian fibrations to the corresponding Hitchin systems.


title: A generalized Fermat-Torricelli tree that has acquired a subconscious on a surface

abstract: We study a generalized Fermat-Torricelli (S.FT) problem for infinitesimal geodesic triangles on a C^2 complete surface M with variable Gaussian curvature a < K < b, for a, b in R, such that the intersection point (generalized Fermat-Torricelli point) of the three geodesics acquires a positive real number (subconscious). The solution of the S.FT problem is a generalized Fermat-Torricelli tree with one node that has acquired a subconscious. This solution is based on a new variational method of the length of a geodesic arc with respect to arc length, which coincides with the first variational formula for geodesics on a surface with K < 0, or 0 < K < c. The 'plasticity' solution of the inverse S.FT problem gives a connection of the absolute value of the Gaussian curvature K(F) at the generalized Fermat-Torricelli point F with the absolute value of the Aleksandrov curvature of the geodesic triangle by acquiring both of them the subconscious of the g.FT point.


title: Applications of the Coase Theorem

abstract: The Coase Theorem has a central place in the theory of environmental economics and regulation. But its applicability for solving real-world externality problems remains debated. In this paper, we first place this seminal contribution in its historical context. We then survey the experimental literature that has tested the importance of the many, often tacit assumptions in the Coase Theorem in the laboratory. We discuss a selection of applications of the Coase Theorem to actual environmental problems, distinguishing between situations in which the polluter or the pollutee pays. While limited in scope, Coasian bargaining over externalities offers a pragmatic solution to problems that are difficult to solve in any other way.


title: The intermediate nebular phase of SN 2014J: onset of clumping as the source of recombination

abstract: At the age of about 1 year, the spectra of most Type Ia supernovae are dominated by strong forbidden nebular emission lines of FeII and FeIII. Later observations (at about 2 years) of the nearby SN 2011fe showed an unexpected shift of ionization to FeI and FeII. Spectra of the very nearby SN Ia 2014J at an intermediate phase (1 - 1.5 years) that are presented here show a progressive decline of FeIII emission, while FeI is not yet strong. The decrease in ionization can be explained if the degree of clumping in the ejecta increases significantly at ~1.5 years, at least in the Fe-dominated zone. Models suggest that clumps remain coherent after about one year, behaving like shrapnel. The high density in the clumps, combined with the decreasing heating rate, would then cause recombination. These data may witness the phase of transition from relatively smooth ejecta to the very clumpy morphology that is typical of SN remnants. The origin of the increased clumping may be the development of local magnetic fields.


title: GeneCAI: Genetic Evolution for Acquiring Compact AI

abstract: In the contemporary big data realm, Deep Neural Networks (DNNs) are evolving towards more complex architectures to achieve higher inference accuracy. Model compression techniques can be leveraged to efficiently deploy such compute-intensive architectures on resource-limited mobile devices. Such methods comprise various hyper-parameters that require per-layer customization to ensure high accuracy. Choosing such hyper-parameters is cumbersome as the pertinent search space grows exponentially with model layers. This paper introduces GeneCAI, a novel optimization method that automatically learns how to tune per-layer compression hyper-parameters. We devise a bijective translation scheme that encodes compressed DNNs to the genotype space. The optimality of each genotype is measured using a multi-objective score based on accuracy and number of floating point operations. We develop customized genetic operations to iteratively evolve the non-dominated solutions towards the optimal Pareto front, thus, capturing the optimal trade-off between model accuracy and complexity. GeneCAI optimization method is highly scalable and can achieve a near-linear performance boost on distributed multi-GPU platforms. Our extensive evaluations demonstrate that GeneCAI outperforms existing rule-based and reinforcement learning methods in DNN compression by finding models that lie on a better accuracy-complexity Pareto curve.


title: An Improved Cutting Plane Method for Convex Optimization, Convex-Concave Games and its Applications

abstract: Given a separation oracle for a convex set $K \subset \mathbb{R}^n$ that is contained in a box of radius $R$, the goal is to either compute a point in $K$ or prove that $K$ does not contain a ball of radius $\epsilon$. We propose a new cutting plane algorithm that uses an optimal $O(n \log (\kappa))$ evaluations of the oracle and an additional $O(n^2)$ time per evaluation, where $\kappa = nR/\epsilon$. $\bullet$ This improves upon Vaidya's $O( \text{SO} \cdot n \log (\kappa) + n^{\omega+1} \log (\kappa))$ time algorithm [Vaidya, FOCS 1989a] in terms of polynomial dependence on $n$, where $\omega < 2.373$ is the exponent of matrix multiplication and $\text{SO}$ is the time for oracle evaluation. $\bullet$ This improves upon Lee-Sidford-Wong's $O( \text{SO} \cdot n \log (\kappa) + n^3 \log^{O(1)} (\kappa))$ time algorithm [Lee, Sidford and Wong, FOCS 2015] in terms of dependence on $\kappa$. For many important applications in economics, $\kappa = \Omega(\exp(n))$ and this leads to a significant difference between $\log(\kappa)$ and $\mathrm{poly}(\log (\kappa))$. We also provide evidence that the $n^2$ time per evaluation cannot be improved and thus our running time is optimal. A bottleneck of previous cutting plane methods is to compute leverage scores, a measure of the relative importance of past constraints. Our result is achieved by a novel multi-layered data structure for leverage score maintenance, which is a sophisticated combination of diverse techniques such as random projection, batched low-rank update, inverse maintenance, polynomial interpolation, and fast rectangular matrix multiplication. Interestingly, our method requires a combination of different fast rectangular matrix multiplication algorithms.


title: DAG With Omitted Objects Displayed (DAGWOOD): A framework for revealing causal assumptions in DAGs

abstract: Directed acyclic graphs (DAGs) are frequently used in epidemiology as a method to encode causal inference assumptions. We propose the DAGWOOD framework to bring many of those encoded assumptions to the forefront. DAGWOOD combines a root DAG (the DAG in the proposed analysis) and a set of branch DAGs (alternative hidden assumptions to the root DAG). All branch DAGs share a common ruleset, and must 1) change the root DAG, 2) be a valid DAG, and either 3a) change the minimally sufficient adjustment set or 3b) change the number of frontdoor paths. Branch DAGs comprise a list of assumptions which must be justified as negligible. We define two types of branch DAGs: exclusion branch DAGs add a single- or bidirectional pathway between two nodes in the root DAG (e.g. direct pathways and colliders), while misdirection branch DAGs represent alternative pathways that could be drawn between objects (e.g., creating a collider by reversing the direction of causation for a controlled confounder). The DAGWOOD framework 1) organizes causal model assumptions, 2) reinforces best DAG practices, 3) provides a framework for evaluation of causal models, and 4) can be used for generating causal models.


title: Projective modules and the homotopy classification of $(G,n)$-complexes

abstract: A $(G,n)$-complex is an $n$-dimensional CW-complex with fundamental group $G$ and whose universal cover is $(n-1)$-connected. If $G$ has periodic cohomology then, for appropriate $n$, we show that there is a one-to-one correspondence between the homotopy types of finite $(G,n)$-complexes and the orbits of the stable class of a certain projective $\mathbb{Z} G$-module under the action of $\text{Aut}(G)$. We develop techniques to compute this action explicitly and use this to give an example where the action is non-trivial.


title: Extragalactic cosmic rays diffusing from two populations of sources

abstract: We consider the possibility of explaining the observed spectrum and composition of the cosmic rays with energies above $10^{17}$ eV in terms of two different extragalactic populations of sources in the presence of a turbulent intergalactic magnetic field (including also a fading Galactic cosmic-ray component). The populations are considered to be the superposition of different nuclear species having rigidity dependent spectra. The first extragalactic population is dominant in the energy range $10^{17}$ to $10^{18}$ eV and consists of sources having a relatively large density ($> 10^{-3}$ Mpc$^{-3}$) and a steep spectrum. The second extragalactic population dominates the cosmic ray flux above few EeV, it has a harder spectral slope and has a high-energy cutoff at few $Z$ EeV (where $eZ$ is the associated cosmic ray charge). This population has a lower density of sources ($<10^{-4}$ Mpc$^{-3}$), so that the typical intersource separation is larger than few tens of Mpc, being significantly affected by a magnetic horizon effect that strongly suppresses its flux for energies below $\sim Z$ EeV. We discuss how this scenario could be reconciled with the values of the cosmic-ray source spectral indices that are expected to result from the diffusive shock acceleration mechanism.


title: Posterior computation with the Gibbs zig-zag sampler

abstract: An intriguing new class of piecewise deterministic Markov processes (PDMPs) has recently been proposed as an alternative to Markov chain Monte Carlo (MCMC). In order to facilitate the application to a larger class of problems, we propose a new class of PDMPs termed Gibbs zig-zag samplers, which allow parameters to be updated in blocks with a zig-zag sampler applied to certain parameters and traditional MCMC-style updates to others. We demonstrate the flexibility of this framework on posterior sampling for logistic models with shrinkage priors for high-dimensional regression and random effects and provide conditions for geometric ergodicity and the validity of a central limit theorem.


title: Qubit transformations on Rashba ring with periodic potential

abstract: A spin-qubit transformation protocol is proposed for an electron in a mesoscopic quantum ring with tunable Rashba interaction controlled by the external electric field. The dynamics of an electron driven around the ring by a series of Landau-Zenner-like transitions between a finite number of local voltage gates is determined analytically. General single-qubit transformations are demonstrated to be feasible in a dynamical basis of localized pseudo-spin states. It is also demonstrated that by the use of suitable protocols based on changes of the Rashba interaction full Bloch sphere can be covered. The challenges of a possible realization of the proposed system in semiconductor heterostructures are discussed.


title: Federated Multi-view Matrix Factorization for Personalized Recommendations

abstract: We introduce the federated multi-view matrix factorization method that extends the federated learning framework to matrix factorization with multiple data sources. Our method is able to learn the multi-view model without transferring the user's personal data to a central server. As far as we are aware this is the first federated model to provide recommendations using multi-view matrix factorization. The model is rigorously evaluated on three datasets on production settings. Empirical validation confirms that federated multi-view matrix factorization outperforms simpler methods that do not take into account the multi-view structure of the data, in addition, it demonstrates the usefulness of the proposed method for the challenging prediction tasks of cold-start federated recommendations.


title: Incidence weighting estimation under bipartite incidence graph sampling

abstract: Bipartite incidence graph sampling provides a unified representation of many sampling situations for the purpose of estimation, including the existing unconventional sampling methods, such as indirect, network or adaptive cluster sampling, which are not originally described as graph problems. We develop a large class of linear estimators based on the edges in the sample bipartite incidence graph, subjected to a general condition of design unbiasedness. The class contains as special cases the classic Horvitz-Thompson estimator, as well as the other unbiased estimators in the literature of unconventional sampling, which can be traced back to Birnbaum and Sirken (1965). Our generalisation allows one to devise other unbiased estimators, thereby providing a potential of efficiency gains in applications. Illustrations are given for adaptive cluster sampling, line-intercept sampling and simulated graphs.


title: Estimating Fiber Orientation Distribution through Blockwise Adaptive Thresholding with Application to HCP Young Adults Data

abstract: Due to recent technological advances, large brain imaging data sets can now be collected. Such data are highly complex so extraction of meaningful information from them remains challenging. Thus, there is an urgent need for statistical procedures that are computationally scalable and can provide accurate estimates that capture the neuronal structures and their functionalities. We propose a fast method for estimating the fiber orientation distribution(FOD) based on diffusion MRI data. This method models the observed dMRI signal at any voxel as a convolved and noisy version of the underlying FOD, and utilizes the spherical harmonics basis for representing the FOD, where the spherical harmonic coefficients are adaptively and nonlinearly shrunk by using a James-Stein type estimator. To further improve the estimation accuracy by enhancing the localized peaks of the FOD, as a second step a super-resolution sharpening process is then applied. The resulting estimated FODs can be fed to a fiber tracking algorithm to reconstruct the white matter fiber tracts. We illustrate the overall methodology using both synthetic data and data from the Human Connectome Project.


title: A new Topology of a High-Voltage-Gain DC-DC Converter Based on Modified Greinacher Voltage Multiplier

abstract: This paper proposes a new topology of a non-isolated two-phase interleaved boost DC-DC converters with modified Greinacher voltage multiplier cells (VMCs) in continuous conduction mode (CCM). The proposed topology can quickly achieve a high gain with low stress on components and proper duty cycle. Moreover, this topology requires a small inductance size in CCM while maintaining the continuity of the input current due to the interleaved technique used in this converter. Because of these features, this topology can be used in many renewable applications such as PV, DC microgrids, etc. The design procedures and component selection are presented in this paper. Furthermore, the proposed converter was simulated to convert an input voltage of 20 V to output voltage of 400 V with an output power of 200 W. Finally, experimental results have been provided to confirm the analysis and simulation results.


title: Ultra-sensitive refractive index gas sensor with functionalized silicon nitride photonic circuits

abstract: Portable and cost-effective gas sensors are gaining demand for a number of environmental, biomedical and industrial applications, yet current devices are confined into specialized labs and cannot be extended to general use. Here, we demonstrate a part-per-billion-sensitive refractive index gas sensor on a photonic chip based on silicon nitride waveguides functionalized with a mesoporous silica top-cladding layer. Low-concentration chemical vapors are detected by monitoring the output spectral pattern of an integrated Mach-Zehnder interferometer having one coated arm exposed to the gas vapors. We retrieved a limit of detection of 65 ppb, 247 ppb and 1.6 ppb for acetone, isopropyl alcohol and ethanol, respectively. Our on-chip refractive index sensor provides, to the best of our knowledge, an unprecedented sensitivity for low gas concentrations based on photonic integrated circuits. As such, our results herald the implementation of compact, portable and inexpensive devices for on-site and real-time environmental monitoring and medical diagnostics.


title: Ti-alloying of BaZrS3 chalcogenide perovskite for photovoltaics

abstract: BaZrS3, a prototypical chalcogenide perovskite, has been shown to possess a direct band gap, an exceptionally strong near band edge light absorption, and good carrier transport. Coupled with its great stability, non-toxicity with earth abundant elements, it is thus a promising candidate for thin film solar cells. However, its reported band gap in the range of 1.7-1.8 eV is larger than the optimal value required to reach the Shockley-Queisser limit of a single junction solar cell. Here we report the synthesis of Ba(Zr1-xTix)S3 perovskite compounds with a reduced band gap. It is found that Ti alloying is extremely effective in band gap reduction of BaZrS3: a mere 4 at% alloying decreases the band gap from 1.78 to 1.51 eV, resulting in a theoretical maximum power conversion efficiency of 32%. Higher Ti-alloying concentration is found to destabilize the distorted chalcogenide perovskite phase.


title: How and why non smooth solutions of the 3D Navier-Stokes equations could possibly develop

abstract: Fluid configurations in three-dimensions, displaying a plausible decay of regularity in a finite time, are suitably built and examined. Vortex rings are the primary ingredients in this study. The full Navier-Stokes system is converted into a 3D scalar problem, where appropriate numerical methods are implemented in order to figure out the behavior of the solutions. Further simplifications in 2D and 1D provide interesting toy problems, that may be used as a starting platform for a better understanding of blowup phenomena.


title: Strong chemical tagging with APOGEE: 21 candidate star clusters that have dissolved across the Milky Way disc

abstract: Chemically tagging groups of stars born in the same birth cluster is a major goal of spectroscopic surveys. To investigate the feasibility of such strong chemical tagging, we perform a blind chemical tagging experiment on abundances measured from APOGEE survey spectra. We apply a density-based clustering algorithm to the eight dimensional chemical space defined by [Mg/Fe], [Al/Fe], [Si/Fe], [K/Fe], [Ti/Fe], [Mn/Fe], [Fe/H], and [Ni/Fe], abundances ratios which together span multiple nucleosynthetic channels. In a high quality sample of 182,538 giant stars, we detect twenty-one candidate clusters with more than fifteen members. Our candidate clusters are more chemically homogeneous than a population of non-member stars with similar [Mg/Fe] and [Fe/H], even in abundances not used for tagging. Group members are consistent with having the same age and fall along a single stellar-population track in logg vs. Teff space. Each group's members are distributed over multiple kpc, and the spread in their radial and azimuthal actions increases with age. We qualitatively reproduce this increase using N-body simulations of cluster dissolution in Galactic potentials that include transient winding spiral arms. Observing our candidate birth clusters with high-resolution spectroscopy in other wavebands to investigate their chemical homogeneity in other nucleosynthetic groups will be essential to confirming the efficacy of strong chemical tagging. Our initially spatially-compact but now widely dispersed candidate clusters will provide novel limits on chemical evolution and orbital diffusion in the Galactic disc, and constraints on star formation in loosely-bound groups.


title: A geometric model of gravitating mass formation and baryon mass spectrum

abstract: The operator of elementary particle mass spectrum is constructed on the basis of the postulate of a six-dimensional space-time manifold having a Riemann-Cartan space structure, proceeding only from its geometric properties. Special solutions found for this operator well describe the mass spectrum of the known quark triplet and hence of some baryon supermultiplets.


title: Extremist ideology as a complex contagion: the spread of far-right radicalization in the United States between 2005-2017

abstract: Increasing levels of far-right extremist violence have generated public concern about the spread of radicalization in the United States. Previous research suggests that radicalized individuals are destabilized by various environmental (or endemic) factors, exposed to extremist ideology, and subsequently reinforced by members of their community. As such, the spread of radicalization may proceed through a social contagion process, in which extremist ideologies behave like complex contagions that require multiple exposures for adoption. In this study, I applied an epidemiological method called two-component spatio-temporal intensity modeling to data from 416 far-right extremists exposed in the United States between 2005 and 2017. The results indicate that patterns of far-right radicalization in the United States are consistent with a complex contagion process, in which reinforcement is required for transmission. Both social media usage and group membership enhance the spread of extremist ideology, suggesting that online and physical organizing remain primary recruitment tools of the far-right movement. Additionally, I identified several endemic factors, such as poverty, that increase the probability of radicalization in particular regions. Future research should investigate how specific interventions, such as online counter-narratives to battle propaganda, may be effectively implemented to mitigate the spread of far-right extremism in the United States.


title: Interplay between phonon downconversion efficiency, density of states at Fermi energy, and intrinsic energy resolution for microwave kinectic inductance detectors

abstract: Microwave Kinetic Inductance detectors (MKIDs) have been recognized as a powerful new tool for single photon detection. These highly multiplexed superconducting devices give timing and energy measurement for every detected photon. However, the full potential of MKID single photon spectroscopy has not been reached , the achieved energy resolution is lower than expected from first principles. Here, we study the efficiency in the phonon downconversion process following the absorption of energetic photons in MKIDs. Assuming previously published material properties, we measure an average downconversion efficiency for three TiN resonators is $\eta$=0.14. We discuss how this efficiency can impact the intrinsic energy resolution of MKID, and how any uncertainty in the unknown density of electron states at the Fermi energy directly affects the efficiency estimations.


title: On Weighted Generalized Entropy for Double Truncated Distribution

abstract: The notion of weighted Renyi's entropy for truncated random variables has recently been proposed in the information-theoretic literature. In this paper, we introduce a generalized measure of it for double truncated distribution, namely weighted generalized interval entropy (WGIE), and study it in the context of reliability analysis. Several properties, including monotonicity, bounds and uniqueness of WGIE are investigated. Moreover, a simulation study is carried out to demonstrate the performance of the estimates of the proposed measure using simulated and real data sets. The role of WGIE in reliability modeling has also been investigated for a real-life problem.


title: Enhancing quantum vacuum signatures with tailored laser beams

abstract: We demonstrate that tailored laser beams provide a powerful means to make quantum vacuum signatures in strong electromagnetic fields accessible in experiment. Typical scenarios aiming at the detection of quantum vacuum nonlinearities at the high-intensity frontier envision the collision of focused laser pulses. The effective interaction of the driving fields mediated by vacuum fluctuations gives rise to signal photons encoding the signature of quantum vacuum nonlinearity. Isolating a small number of signal photons from the large background of the driving laser photons poses a major experimental challenge. The main idea of the present work is to modify the far-field properties of a driving laser beam to exhibit a field-free hole in its center, thereby allowing for an essentially background free measurement of the signal scattered in the forward direction. Our explicit construction makes use of a peculiar far-field/focus duality.


title: Scaling laws for the geometry of an impact-induced magma ocean

abstract: Here, we develop scaling laws for (1) the distribution of impact-induced heat within the mantle and (2) shape of the impact-induced melt based on more than 100 smoothed particle hydrodynamic (SPH) simulations. We use Legendre polynomials to describe these scaling laws and determine their coefficients by linear regression, minimizing the error between our model and SPH simulations. The input parameters are the impact angle $\theta$ ($0^{\circ}, 30^{\circ}, 60^{\circ}$, and $90^{\circ}$), total mass $M_T$ ($1M_{\rm Mars}-53M_{\rm Mars}$, where $M_{\rm Mars}$ is the mass of Mars), impact velocity $v_{\rm imp}$ ($v_{\rm esc} - 2v_{\rm esc}$, where $v_{\rm esc}$ is the mutual escape velocity), and impactor-to-total mass ratio $\gamma$ ($0.03-0.5$). We find that the equilibrium pressure at the base of a melt pool can be higher (up to $\approx 80 \%$) than those of radially-uniform global magma ocean models. This could have a significant impact on element partitioning. These melt scaling laws are publicly available on GitHub ($\href{https://github.com/mikinakajima/MeltScalingLaw}{https://github.com/mikinakajima/MeltScalingLaw}$).


title: The Spotify Podcast Dataset

abstract: Podcasts are a relatively new form of audio media. Episodes appear on a regular cadence, and come in many different formats and levels of formality. They can be formal news journalism or conversational chat; fiction or non-fiction. They are rapidly growing in popularity and yet have been relatively little studied. As an audio format, podcasts are more varied in style and production types than, say, broadcast news, and contain many more genres than typically studied in video research. The medium is therefore a rich domain with many research avenues for the IR and NLP communities. We present the Spotify Podcast Dataset, a set of approximately 100K podcast episodes comprised of raw audio files along with accompanying ASR transcripts. This represents over 47,000 hours of transcribed audio, and is an order of magnitude larger than previous speech-to-text corpora.


title: Simple Formulas for a Two-Tier Strategy to Flatten the Curve

abstract: Some basic facts that enable or constrain efforts to "flatten the curve" by non-pharmaceutical interventions (NPI) are summarized and placed in the current context of the COVID-19 epidemic. Analytical formulas are presented for a simple two-tier NPI strategy that places different social distancing targets for high- and low-risk groups. The results aim to facilitate rapid what-if analysis rather than replace detailed simulation, which remains indispensable for informed decision making. The appearance of the oft neglected Lambert function in the resulting formulas is noted.


title: Deep Learning and Open Set Malware Classification: A Survey

abstract: As the Internet is growing rapidly these years, the variant of malicious software, which often referred to as malware, has become one of the major and serious threats to Internet users. The dramatic increase of malware has led to a research area of not only using cutting edge machine learning techniques classify malware into their known families, moreover, recognize the unknown ones, which can be related to Open Set Recognition (OSR) problem in machine learning. Recent machine learning works have shed light on Open Set Recognition (OSR) from different scenarios. Under the situation of missing unknown training samples, the OSR system should not only correctly classify the known classes, but also recognize the unknown class. This survey provides an overview of different deep learning techniques, a discussion of OSR and graph representation solutions and an introduction of malware classification systems.


title: Effect of wave versus particle phonon nature in thermal transport through nanostructures

abstract: Comprehensive understanding of thermal transport in nanostructured materials needs large scale simulations bridging length scales dictated by different physics related to the wave versus particle nature of phonons. Yet, available computational approaches implicitly treat phonons as either just waves or as particles. In this work, using a full wave-based Non-Equilibrium Green's Function (NEGF) method, and a particle-based ray-tracing Monte Carlo (MC) approach, we investigate the qualitative differences in the wave and particle-based phonon transport at the vicinity of nanoscale features. For the simple example of a nanoporous geometry, we show that phonon transmission agrees very well for both methods with an error margin of approximately 15%, across phonon wavelengths even for features with sizes down to 3-4 nm. For cases where phonons need to squeeze in smaller regions to propagate, we find that MC underestimates the transmission of long wavelength phonons whereas wave treatment within NEGF indicates that those long wavelength phonons can propagate more easily. We also find that particle-based simulation methods are somewhat more sensitive to structural variations compared to the wave-based NEGF method. The insight extracted from comparing wave and particle methods can be used to provide a better and more complete understanding of phonon transport in nanomaterials.


title: Convergence Results for Implicit--Explicit General Linear Methods

abstract: This paper studies fixed-step convergence of implicit-explicit general linear methods. We focus on a subclass of schemes that is internally consistent, has high stage order, and favorable stability properties. Classical, index-1 differential algebraic equation, and singular perturbation convergence analyses results are given. For all these problems IMEX GLMs from the class of interest converge with the full theoretical orders under general assumptions. The convergence results require the time steps to be sufficiently small, with upper bounds that are independent on the stiffness of the problem.


title: Ensemble Kalman Filter with perturbed observations in weather forecasting and data assimilation

abstract: Data assimilation provides algorithms for widespread applications in various fields. It is of practical use to deal with a large amount of information in the complex system that is hard to estimate. Weather forecasting is one of the applications, where the prediction of meteorological data are corrected given the observations. Numerous approaches are contained in data assimilation. One specific sequential method is the Kalman Filter. The core is to estimate unknown information with the new data that is measured and the prior data that is predicted. As a matter of fact, there are different improved methods in the Kalman Filter. In this project, the Ensemble Kalman Filter with perturbed observations is considered. It is achieved by Monte Carlo simulation. In this method, the ensemble is involved in the calculation instead of the state vectors. In addition, the measurement with perturbation is viewed as the suitable observation. These changes compared with the Linear Kalman Filter make it more advantageous in that applications are not restricted in linear systems any more and less time is taken when the data are calculated by computers. The thesis seeks to develop the Ensemble Kalman Filter with perturbed observation gradually. With the Mathematical preliminaries including the introduction of dynamical systems, the Linear Kalman Filter is built. Meanwhile, the prediction and analysis processes are derived. After that, we use the analogy thoughts to lead in the non-linear Ensemble Kalman Filter with perturbed observations. Lastly, a classic Lorenz 63 model is illustrated by MATLAB. In the example, we experiment on the number of ensemble members and seek to investigate the relationships between the error of variance and the number of ensemble members. We reach the conclusion that on a limited scale the larger number of ensemble members indicates the smaller error of prediction.


title: nPINNs: nonlocal Physics-Informed Neural Networks for a parametrized nonlocal universal Laplacian operator. Algorithms and Applications

abstract: Physics-informed neural networks (PINNs) are effective in solving inverse problems based on differential and integral equations with sparse, noisy, unstructured, and multi-fidelity data. PINNs incorporate all available information into a loss function, thus recasting the original problem into an optimization problem. In this paper, we extend PINNs to parameter and function inference for integral equations such as nonlocal Poisson and nonlocal turbulence models, and we refer to them as nonlocal PINNs (nPINNs). The contribution of the paper is three-fold. First, we propose a unified nonlocal operator, which converges to the classical Laplacian as one of the operator parameters, the nonlocal interaction radius $\delta$ goes to zero, and to the fractional Laplacian as $\delta$ goes to infinity. This universal operator forms a super-set of classical Laplacian and fractional Laplacian operators and, thus, has the potential to fit a broad spectrum of data sets. We provide theoretical convergence rates with respect to $\delta$ and verify them via numerical experiments. Second, we use nPINNs to estimate the two parameters, $\delta$ and $\alpha$. The strong non-convexity of the loss function yielding multiple (good) local minima reveals the occurrence of the operator mimicking phenomenon: different pairs of estimated parameters could produce multiple solutions of comparable accuracy. Third, we propose another nonlocal operator with spatially variable order $\alpha(y)$, which is more suitable for modeling turbulent Couette flow. Our results show that nPINNs can jointly infer this function as well as $\delta$. Also, these parameters exhibit a universal behavior with respect to the Reynolds number, a finding that contributes to our understanding of nonlocal interactions in wall-bounded turbulence.


title: GRANDMA Observations of Advanced LIGO's and Advanced Virgo's Third Observational Campaign

abstract: GRANDMA is a network of 25 telescopes of different sizes, including both photometric and spectroscopic facilities. The network aims to coordinate follow-up observations of gravitational-wave candidate alerts, especially those with large localisation uncertainties, to reduce the delay between the initial detection and the optical confirmation. In this paper, we detail GRANDMA's observational performance during Advanced LIGO/Advanced Virgo Observing Run 3 (O3), focusing on the second part of O3; this includes summary statistics pertaining to coverage and possible astrophysical origin of the candidates. To do so, we quantify our observation efficiency in terms of delay between gravitational-wave candidate trigger time, observations, and the total coverage. Using an optimised and robust coordination system, GRANDMA followed-up about 90 % of the gravitational-wave candidate alerts, i.e. 49 out of 56 candidates. This led to coverage of over 9000 deg2 during O3. The delay between the gravitational-wave candidate trigger and the first observation was below 1.5 hour for 50 % of the alerts. We did not detect any electromagnetic counterparts to the gravitational-wave candidates during O3, likely due to the very large localisation areas (on average thousands of degrees squares) and relatively large distance of the candidates (above 200 Mpc for 60 % of BNS candidates). We derive constraints on potential kilonova properties for two potential binary neutron star coalescences (GW190425 and S200213t), assuming that the events' locations were imaged.