title: $\mathrm{P}$-Optimal Proof Systems for Each $\mathrm{coNP}$-Complete Set and no Complete Problems in $\mathrm{NP}\cap\mathrm{coNP}$ Relative to an Oracle

abstract: We build on a working program initiated by Pudl\'ak [Pud17] and construct an oracle relative to which each $\mathrm{coNP}$-complete set has $\mathrm{P}$-optimal proof systems and $\mathrm{NP}\cap\mathrm{coNP}$ does not have complete problems.


title: \'Equir\'epartition de sommes exponentielles (travaux de Katz)

abstract: Many exponential sums over finite fields, including Gauss sums and Kloosterman sums, arise as the Fourier transform with respect to a character of the trace function of an $\ell$-adic sheaf on a commutative algebraic group. We study the equidistribution of these sums when the sheaf is fixed but the character varies over larger and larger extensions of the finite field. For the additive group, monodromy governs equidistribution by a theorem of Deligne. A few years ago, Katz solved the multiplicative variant of the question in a work where Tannakian ideas play an essential role.


title: Coupling Oceanic Observation Systems to Study Mesoscale Ocean Dynamics

abstract: Understanding local currents in the North Atlantic region of the ocean is a key part of modelling heat transfer and global climate patterns. Satellites provide a surface signature of the temperature of the ocean with a high horizontal resolution while in situ autonomous probes supply high vertical resolution, but horizontally sparse, knowledge of the ocean interior thermal structure. The objective of this paper is to develop a methodology to combine these complementary ocean observing systems measurements to obtain a three-dimensional time series of ocean temperatures with high horizontal and vertical resolution. Within an observation-driven framework, we investigate the extent to which mesoscale ocean dynamics in the North Atlantic region may be decomposed into a mixture of dynamical modes, characterized by different local regressions between Sea Surface Temperature (SST), Sea Level Anomalies (SLA) and Vertical Temperature fields. Ultimately we propose a Latent-class regression method to improve prediction of vertical ocean temperature.


title: Enlarging the scope of resonant di-Higgs searches: Hunting for Higgs-to-Higgs cascades in $4b$ final states at the LHC and future colliders

abstract: We extend the coverage of resonant di-Higgs searches in the $b \bar{b} b \bar{b}$ final state to the process $p p \to H_1 \to H_2 H_2 \to b \bar{b} b \bar{b}$, where both $H_{1,2}$ are spin-$0$ states beyond the Standard Model. Such a process constitutes a joint discovery mode for the new states $H_1$ and $H_2$. We present the first sensitivity study of this channel, using public LHC data to validate our analysis. We also provide a first estimate of the sensitivity of the search for the HL-LHC and future facilities like the HE-LHC and FCC-hh. We analyze the discovery potential of this search for several non-minimal scalar sector scenarios: an extension of the SM with two extra singlet scalar fields, the two-Higgs-doublet model and a two-Higgs doublet model plus a singlet, which captures the scalar potential features of the NMSSM. We find that this channel represents a novel, very powerful probe for extended Higgs sectors, offering complementary sensitivity to existing analyses.


title: Modelling Online Comment Threads from their Start

abstract: The social Web is a widely used platform for online discussion. Across social media, users can start discussions by posting a topical image, url, or message. Upon seeing this initial post, other users may add their own comments to the post, or to another user's comment. The resulting online discourse produces a comment thread, which constitutes an enormous portion of modern online communication. Comment threads are often viewed as trees: nodes represent the post and its comments, while directed edges represent reply-to relationships. The goal of the present work is to predict the size and shape of these comment threads. Existing models do this by observing the first several comments and then fitting a predictive model. However, most comment threads are relatively small, and waiting for data to materialize runs counter to the goal of the prediction task. We therefore introduce the Comment Thread Prediction Model (CTPM) that accurately predicts the size and shape of a comment thread using only the text of the initial post, allowing for the prediction of new posts without observable comments. We find that the CTPM significantly outperforms existing models and competitive baselines on thousands of Reddit discussions from nine varied subreddits, particularly for new posts.


title: Resistive contribution in electrical switching experiments with antiferromagnets

abstract: Recent research demonstrated the electrical switching of antiferromagnets via intrinsic spin-orbit torque or the spin Hall effect of an adjacent heavy metal layer. The electrical readout is typically realized by measuring the transverse anisotropic magnetoresistance at planar cross- or star-shaped devices with four or eight arms, respectively. Depending on the material, the current density necessary to switch the magnetic state can be large, often close to the destruction threshold of the device. We demonstrate that the resulting electrical stress changes the film resistivity locally and thereby breaks the fourfold rotational symmetry of the conductor. This symmetry breaking due to film inhomogeneity produces signals, that resemble the anisotropic magnetoresistance and is experimentally seen as a "saw-tooth"-shaped transverse resistivity. This artifact can persist over many repeats of the switching experiment and is not easily separable from the magnetic contribution. We discuss the origin of the artifact, elucidate the role of the film crystallinity, and propose approaches how to separate the resistive contribution from the magnetic contribution.


title: Weak-field induced nonmagnetic state in a Co-based honeycomb

abstract: Layered honeycomb magnets are of interest as potential realizations of the Kitaev quantum spin liquid (KQSL), a quantum state with long-range spin entanglement and an exactly solvable Hamiltonian. Conventional magnetically ordered states are present for all currently known candidate materials, however, because non-Kitaev terms in the Hamiltonians obscure the Kitaev physics. Current experimental studies of the KQSL are focused on 4d- or 5d-transition-metal-based honeycombs, in which strong spin-orbit coupling can be expected, yielding Kitaev interaction that dominate in an applied magnetic field. In contrast, for 3d-based layered honeycomb magnets, spin orbit coupling is weak and thus Kitaev-physics should be substantially less accessible. Here we report our studies on BaCo2(AsO4)2, for which we find that the magnetic order associated with the non-Kitaev interactions can be fully suppressed by a relatively low magnetic field, yielding a non-magnetic material and implying the presence of strong magnetic frustration and weak non-Kitaev interactions.


title: Bootstrapping the $\mathcal{N}=1$ Wess-Zumino models in three dimensions

abstract: Using numerical bootstrap method, we determine the critical exponents of the minimal three-dimensional N = 1 Wess-Zumino models with cubic superpotetential $W\sim d_{ijk}\Phi_i\Phi_j\Phi_k$. The tensor $d_{ijk}$ is taken to be the invariant tensor of either permutation group $S_N$, special unitary group $SU(N)$, or a series of groups called F4 family of Lie groups. Due to the equation of motion, at the Wess-Zumino fixed point, the operator $d_{ijk}\Phi_i\Phi_j\Phi_k$ is a (super)descendant of $\Phi_i$ . We observe such super-multiplet recombination in numerical bootstrap, which allows us to determine the scaling dimension of the super-field $\Phi_i$.


title: Towards Interpretable Graph Modeling with Vertex Replacement Grammars

abstract: An enormous amount of real-world data exists in the form of graphs. Oftentimes, interesting patterns that describe the complex dynamics of these graphs are captured in the form of frequently reoccurring substructures. Recent work at the intersection of formal language theory and graph theory has explored the use of graph grammars for graph modeling and pattern mining. However, existing formulations do not extract meaningful and easily interpretable patterns from the data. The present work addresses this limitation by extracting a special type of vertex replacement grammar, which we call a KT grammar, according to the Minimum Description Length (MDL) heuristic. In experiments on synthetic and real-world datasets, we show that KT-grammars can be efficiently extracted from a graph and that these grammars encode meaningful patterns that represent the dynamics of the real-world system.


title: An empirical model of Energetic Neutral Atom imaging of the heliosphere and its implications for future heliospheric missions at great heliocentric distances

abstract: Several concepts for heliospheric missions operating at heliocentric distances far beyond Earth orbit are currently investigated by the scientific community. The mission concept of the Interstellar Probe (McNutt et al. 2018), e.g., aims at reaching a distance of 1000 au away from the Sun within this century. This would allow the coming generation to obtain a global view of our heliosphere from an outside vantage point by measuring the Energetic Neutral Atoms (ENAs) originating from the various plasma regions. It would also allow for direct sampling of unperturbed interstellar medium, and for many observation opportunities beyond heliospheric science, such as visits to Kuiper Belt Objects, a comprehensive view on the interplanetary dust populations, and infrared astronomy free from the foreground emission of the Zodiacal cloud. In this study, we present a simple empirical model of ENAs from the heliosphere and derive basic requirements for ENA instrumentation onboard a spacecraft at great heliocentric distances. We consider the full energy range of heliospheric ENAs from 10 eV to 100 keV because each part of the energy spectrum has its own merits for heliospheric science. To cover the full ENA energy range, two or three different ENA instruments are needed. Thanks to parallax observations, some insights about the nature of the IBEX Ribbon and the dimensions of the heliosphere can already be gained by ENA imaging from a few au heliocentric distance. To directly reveal the global shape of the heliosphere, measurements from outside the heliosphere are, of course, the best option.


title: Towards Quantifying Intrinsic Generalization of Deep ReLU Networks

abstract: Understanding the underlying mechanisms that enable the empirical successes of deep neural networks is essential for further improving their performance and explaining such networks. Towards this goal, a specific question is how to explain the "surprising" behavior of the same over-parametrized deep neural networks that can generalize well on real datasets and at the same time "memorize" training samples when the labels are randomized. In this paper, we demonstrate that deep ReLU networks generalize from training samples to new points via piece-wise linear interpolation. We provide a quantified analysis on the generalization ability of a deep ReLU network: Given a fixed point $\mathbf{x}$ and a fixed direction in the input space $\mathcal{S}$, there is always a segment such that any point on the segment will be classified the same as the fixed point $\mathbf{x}$. We call this segment the $generalization \ interval$. We show that the generalization intervals of a ReLU network behave similarly along pairwise directions between samples of the same label in both real and random cases on the MNIST and CIFAR-10 datasets. This result suggests that the same interpolation mechanism is used in both cases. Additionally, for datasets using real labels, such networks provide a good approximation of the underlying manifold in the data, where the changes are much smaller along tangent directions than along normal directions. On the other hand, however, for datasets with random labels, generalization intervals along mid-lines of triangles with the same label are much smaller than those on the datasets with real labels, suggesting different behaviors along other directions. Our systematic experiments demonstrate for the first time that such deep neural networks generalize through the same interpolation and explain the differences between their performance on datasets with real and random labels.


title: Are white dwarf magnetic fields in close binaries generated during common-envelope evolution?

abstract: Understanding the origin of the magnetic fields in white dwarfs (WDs) has been a puzzle for decades. A scenario that has gained considerable attention in the past years assumes that such magnetic fields are generated through a dynamo process during common-envelope evolution. We performed binary population models using an up-to-date version of the BSE code to confront the predictions of this model with observational results. We found that this hypothesis can only explain the observed distribution of WD magnetic fields in polars and pre-polars and the low-temperature WDs in pre-polars if it is re-scaled to fit the observational data. Furthermore, in its present version, the model fails to explain the absence of young close detached WD+M-dwarf binaries harbouring hot magnetic WDs and predicts that the overwhelming majority of WDs in close binaries should be strongly magnetic, which is also in serious conflict with the observations. We conclude that either the common-envelope dynamo scenario needs to be substantially revised or a different mechanism is responsible for the generation of strong WD magnetic fields in close binaries.

categories: astro-ph.SR astro-ph.EP astro-ph.GA astro-ph.HE


title: The C-Band All-Sky Survey (C-BASS): Total intensity point-source detection over the northern sky

abstract: We present a point-source detection algorithm that employs the second order Spherical Mexican Hat wavelet filter (SMHW2), and use it on C-BASS northern intensity data to produce a catalogue of point-sources. This catalogue allows us to cross-check the C-BASS flux-density scale against existing source surveys, and provides the basis for a source mask which will be used in subsequent C-BASS and cosmic microwave background (CMB) analyses. The SMHW2 allows us to filter the entire sky at once, avoiding complications from edge effects arising when filtering small sky patches. The algorithm is validated against a set of Monte Carlo simulations, consisting of diffuse emission, instrumental noise, and various point-source populations. The simulated source populations are successfully recovered. The SMHW2 detection algorithm is used to produce a $4.76\,\mathrm{GHz}$ northern sky source catalogue in total intensity, containing 1784 sources and covering declinations $\delta\geq-10^{\circ}$. The C-BASS catalogue is matched with the Green Bank 6\,cm (GB6) and Parkes-MIT-NRAO (PMN) catalogues over their areas of common sky coverage. From this we estimate the $90$ per cent completeness level to be approximately $610\,\mathrm{mJy}$, with a corresponding reliability of $98$ per cent, when masking the brightest $30$ per cent of the diffuse emission in the C-BASS northern sky map. We find the C-BASS and GB6 flux-density scales to be consistent with one another to within approximately $4$ per cent.


title: Anisotropy-driven thermal conductivity switching and thermal hysteresis in a ferroelectric

abstract: We present a theoretical proposal for the design of a thermal switch based on the anisotropy of the thermal conductivity of PbTiO3 and of the possibility to rotate the ferroelectric polarization with an external electric field. Our calculations are based on an iterative solution of the phonon Boltzmann Transport Equation and rely on interatomic force constants computed within an efficient second-principles density functional theory scheme. We also characterize the hysteresis cycle of the thermal conductivity in presence of an applied electric field and show that the response time would be limited by speed of the ferroelectric switch itself and thus can operate in the high-frequency regime.


title: Tropical floor plans and enumeration of complex and real multi-nodal surfaces

abstract: The family of complex projective surfaces in projective three space of degree $d$ having precisely $\delta$ nodes as their only singularities has codimension $\delta$ in the linear system of surfaces of degree $d$ for sufficiently large $d$ and is of degree $N_{\delta,complex}(d)=(4(d-1)^3)^\delta/\delta!+O(d^{3\delta-3})$. In particular, this number is polynomial in $d$. By means of tropical geometry, we explicitly describe $(4d^3)^\delta/\delta!+O(d^{3\delta-1})$ surfaces passing through a suitable generic configuration of $n=\binom{d+3}{3}-\delta-1$ points in projective three space. These surfaces are close to tropical limits which we characterize combinatorially, introducing the concept of floor plans for multinodal tropical surfaces. The concept of floor plans is similar to the well-known floor diagrams (a combinatorial tool for tropical curve counts): with it, we keep the combinatorial essentials of a multinodal tropical surface which are sufficient to reconstruct the surface. In the real case, we estimate the range for possible numbers of real multi-nodal surfaces satisfying point conditions. We show that, for a special configuration $w$ of real points, the number $N_{\delta,real}(d,w)$ of real surfaces of degree $d$ having $\delta$ real nodes and passing through $w$ is bounded from below by $(\frac{3}{2}d^3)^\delta/\delta! +O(d^{3\delta-1})$. We prove analogous statements for counts of multinodal surfaces in $P^1\times P^2$ and $P^1\times P^1\times P^1$.


title: Dispersion Properties, Nonlinear Waves and Birefringence in Classical Nonlinear Electrodynamics

abstract: Using the very basic physics principles, we have studied the implications of quantum corrections to classical electrodynamics and the propagation of electromagnetic waves and pulses. The initial nonlinear wave equation for the electromagnetic vector potential is solved perturbatively about the known exact plane wave solution in both the free vacuum case, as well as when a constant magnetic field is applied. A nonlinear wave equation with nonzero convective part for the (relatively) slowly varying amplitude of the first-order perturbation has been derived. This equation governs the propagation of electromagnetic waves with a reduced speed of light, where the reduction is roughly proportional to the intensity of the initial pumping plane wave. A system of coupled nonlinear wave equations for the two slowly varying amplitudes of the first-order perturbation, which describe the two polarization states, has been obtained for the case of constant magnetic field background. Further, the slowly varying wave amplitude behavior is shown to be similar to that of a cnoidal wave, known to describe surface gravity waves in shallow water. It has been demonstrated that the two wave modes describing the two polarization states are independent, and they propagate at different wave frequencies. This effect is usually called nonlinear birefringence.

categories: physics.plasm-ph astro-ph.HE hep-th nlin.PS quant-ph


title: Frame Matroids, Toric Ideals, and a Conjecture of White

abstract: Blasiak verified a conjecture of White for graphic matroids by showing that the toric ideal of a graphic matroid is generated by quadrics. In this paper, we extend this result to frame matroids satisfying a linearity condition. Such classes of matroids include graphic matroids, bicircular matroids, signed graphic matroids, and more generally frame matroids obtained from group-labelled graphs.


title: Improved photometric redshifts with colour-constrained galaxy templates for future wide-area surveys

abstract: Cosmology and galaxy evolution studies with LSST, \Euclid, and {\it Roman}, will require accurate redshifts for the detected galaxies. In this study, we present improved photometric redshift estimates for galaxies using a template library that populates three-color space and is constrained by \HST/CANDELS photometry. For the training sample, we use a sample of galaxies having photometric redshifts which allows us to train on a large, unbiased galaxy sample having deep, unconfused photometry at optical-to-mid infrared wavelengths. Galaxies in the training sample are assigned to cubes in three-dimensional color space, $V-H$, $I-J$, and $z-H$. We then derive the best-fit spectral energy distributions of the training sample at the fixed CANDELS median photometric redshifts to construct the new template library for each individual color cube (i.e. color-cube-based template library). We derive photometric redshifts (photo-z) of our target galaxies using our new color-cube-based template library and with photometry in only a limited set of bands, as expected for the aforementioned surveys. As a result, our method yields $\sigma_{NMAD}$ of 0.026 and an outlier fraction of 6\% using only photometry in the LSST and \Euclid/{\it Roman} bands. This is an improvement of $\sim$10\% on $\sigma_{NMAD}$ and a reduction in outlier fraction of $\sim$13\% compared to other techniques. In particular, we improve the photo-z precision by about 30\% at $2 < z < 3$. We also assess photo-z improvements by including $K$ or mid-infrared bands to the $ugrizYJH$ photometry. Our color-cube-based template library is a powerful tool to constrain photometric redshifts for future large surveys.


title: Decoupling feature propagation from the design of graph auto-encoders

abstract: We present two instances, L-GAE and L-VGAE, of the variational graph auto-encoding family (VGAE) based on separating feature propagation operations from graph convolution layers typically found in graph learning methods to a single linear matrix computation made prior to input in standard auto-encoder architectures. This decoupling enables the independent and fixed design of the auto-encoder without requiring additional GCN layers for every desired increase in the size of a node's local receptive field. Fixing the auto-encoder enables a fairer assessment on the size of a nodes receptive field in building representations. Furthermore a by-product of fixing the auto-encoder design often results in substantially smaller networks than their VGAE counterparts especially as we increase the number of feature propagations. A comparative downstream evaluation on link prediction tasks show comparable state of the art performance to similar VGAE arrangements despite considerable simplification. We also show the simple application of our methodology to more challenging representation learning scenarios such as spatio-temporal graph representation learning.


title: Anderson Acceleration of Proximal Gradient Methods

abstract: Anderson acceleration is a well-established and simple technique for speeding up fixed-point computations with countless applications. Previous studies of Anderson acceleration in optimization have only been able to provide convergence guarantees for unconstrained and smooth problems. This work introduces novel methods for adapting Anderson acceleration to (non-smooth and constrained) proximal gradient algorithms. Under some technical conditions, we extend the existing local convergence results of Anderson acceleration for smooth fixed-point mappings to the proposed scheme. We also prove analytically that it is not, in general, possible to guarantee global convergence of native Anderson acceleration. We therefore propose a simple scheme for stabilization that combines the global worst-case guarantees of proximal gradient methods with the local adaptation and practical speed-up of Anderson acceleration.


title: Solving the measurement problem within standard quantum theory

abstract: A misunderstanding of entangled states has spawned decades of concern about quantum measurements and a plethora of quantum interpretations. The "measurement state" or "Schrodinger's cat state" of a superposed quantum system and its detector is nonlocally entangled, suggesting that we turn to nonlocality experiments for insight into measurements. By studying the full range of superposition phases, these experiments show precisely what the measurement state does and does not superpose. These experiments reveal that the measurement state is not, as had been supposed, a paradoxical superposition of detector states. It is instead a nonparadoxical superposition of two correlations between detector states and system states. In this way, the experimental results resolve the problem of definite outcomes ("Schrodinger's cat"), leading to a resolution of the measurement problem. However, this argument does not yet resolve the measurement problem because it is based on the results of experiments, while measurement is a theoretical problem: How can standard quantum theory explain the definite outcomes seen experimentally? Thus, we summarize the nonlocality experiments' supporting theory, which rigorously predicts the experimental results directly from optical paths. Several previous theoretical analyses of the measurement problem have relied on the reduced density operators derived from the measurement state, but these solutions have been rejected due to criticism of reduced density operators. Because it avoids reduced density operators, the optical-path analysis is immune to such criticism.


title: Automatic Post-Editing for Machine Translation

abstract: Automatic Post-Editing (APE) aims to correct systematic errors in a machine translated text. This is primarily useful when the machine translation (MT) system is not accessible for improvement, leaving APE as a viable option to improve translation quality as a downstream task - which is the focus of this thesis. This field has received less attention compared to MT due to several reasons, which include: the limited availability of data to perform a sound research, contrasting views reported by different researchers about the effectiveness of APE, and limited attention from the industry to use APE in current production pipelines. In this thesis, we perform a thorough investigation of APE as a downstream task in order to: i) understand its potential to improve translation quality; ii) advance the core technology - starting from classical methods to recent deep-learning based solutions; iii) cope with limited and sparse data; iv) better leverage multiple input sources; v) mitigate the task-specific problem of over-correction; vi) enhance neural decoding to leverage external knowledge; and vii) establish an online learning framework to handle data diversity in real-time. All the above contributions are discussed across several chapters, and most of them are evaluated in the APE shared task organized each year at the Conference on Machine Translation. Our efforts in improving the technology resulted in the best system at the 2017 APE shared task, and our work on online learning received a distinguished paper award at the Italian Conference on Computational Linguistics. Overall, outcomes and findings of our work have boost interest among researchers and attracted industries to examine this technology to solve real-word problems.


title: Generative Adversarial Networks And Domain Adaptation For Training Data Independent Image Registration

abstract: Medical image registration is an important task in automated analysis of multi-modal images and temporal data involving multiple patient visits. Conventional approaches, although useful for different image types, are time consuming. Of late, deep learning (DL) based image registration methods have been proposed that outperform traditional methods in terms of accuracy and time. However,DL based methods are heavily dependent on training data and do not generalize well when presented with images of different scanners or anatomies. We present a DL based approach that can perform medical image registration of one image type despite being trained with images of a different type. This is achieved by unsupervised domain adaptation in the registration process and allows for easier application to different datasets without extensive retraining.To achieve our objective we train a network that transforms the given input image pair to a latent feature space vector using autoencoders. The resultant encoded feature space is used to generate the registered images with the help of generative adversarial networks (GANs). This feature transformation ensures greater invariance to the input image type. Experiments on chest Xray, retinal and brain MR images show that our method, trained on one dataset gives better registration performance for other datasets, outperforming conventional methods that do not incorporate domain adaptation


title: Acceleration techniques for semiclassical Maxwell-Bloch systems: An application to discrete quantum dot ensembles

abstract: The solution to Maxwell-Bloch systems using an integral-equation-based framework has proven effective at capturing collective features of laser-driven and radiation-coupled quantum dots, such as light localization and modifications of Rabi oscillations. Importantly, it enables observation of the dynamics of each quantum dot in large ensembles in a rigorous, error-controlled, and self-consistent way without resorting to spatial averaging. Indeed, this approach has demonstrated convergence in ensembles containing up to $10^4$ interacting quantum dots. Scaling beyond $10^4$ quantum dots tests the limit of computational horsepower, however, due to the $\mathcal{O}(N_t N_s^2)$ scaling (where $N_t$ and $N_s$ denote the number of temporal and spatial degrees of freedom). In this work, we present an algorithm that reduces the cost of analysis to $\mathcal{O}(N_t N_s \log^2 N_s)$. While the foundations of this approach rely on well-known particle-particle/particle-mesh and adaptive integral methods, we add refinements specific to transient systems and systems with multiple spatial and temporal derivatives. Accordingly, we offer numerical results that validate the accuracy, effectiveness and utility of this approach in analyzing the dynamics of large ensembles of quantum dots.


title: Identifying the Most Explainable Classifier

abstract: We introduce the notion of pointwise coverage to measure the explainability properties of machine learning classifiers. An explanation for a prediction is a definably simple region of the feature space sharing the same label as the prediction, and the coverage of an explanation measures its size or generalizability. With this notion of explanation, we investigate whether or not there is a natural characterization of the most explainable classifier. According with our intuitions, we prove that the binary linear classifier is uniquely the most explainable classifier up to negligible sets.


title: Wellposedness, Spectral Analysis and Asymptotic Stability of a Multilayered Heat-Wave-Wave System

abstract: n this work we consider a multilayered heat-wave system where a 3-D heat equation is coupled with a 3-D wave equation via a 2-D interface whose dynamics is described by a 2-D wave equation. This system can be viewed as a simplification of a certain fluid-structure interaction (FSI) PDE model where the structure is of composite-type; namely it consists of a \textquotedblleft thin\textquotedblright\ layer and a \textquotedblleft thick\textquotedblright\ layer. We associate the wellposedness of the system with a strongly continuous semigroup and establish its asymptotic decay. Our first result is semigroup well-posedness for the (FSI) PDE dynamics. Utilizing here a Lumer-Phillips approach, we show that the fluid-structure system generates a $C_0$-semigroup on a chosen finite energy space of data. As our second result, we prove that the solution to the (FSI) dynamics generated by the $C_0$-semigroup tends asymptotically to the zero state for all initial data. That is, the semigroup of the (FSI) system is strongly stable. For this stability work, we analyze the spectrum of the generator $\mathbf A$ and show that the spectrum of $\mathbf A$ does not intersect the imaginary axis. \vskip.3cm \noindent \textbf{Key terms:} Fluid-structure interaction, heat-wave system, well-posedness, semigroup, strong stability.


title: Robust Learning Rate Selection for Stochastic Optimization via Splitting Diagnostic

abstract: This paper proposes SplitSGD, a new dynamic learning rate schedule for stochastic optimization. This method decreases the learning rate for better adaptation to the local geometry of the objective function whenever a stationary phase is detected, that is, the iterates are likely to bounce at around a vicinity of a local minimum. The detection is performed by splitting the single thread into two and using the inner product of the gradients from the two threads as a measure of stationarity. Owing to this simple yet provably valid stationarity detection, SplitSGD is easy-to-implement and essentially does not incur additional computational cost than standard SGD. Through a series of extensive experiments, we show that this method is appropriate for both convex problems and training (non-convex) neural networks, with performance compared favorably to other stochastic optimization methods. Importantly, this method is observed to be very robust with a set of default parameters for a wide range of problems and, moreover, can yield better generalization performance than other adaptive gradient methods such as Adam.


title: Neutron spin dynamics in polarized targets

abstract: We present neutron elastic scattering amplitude for arbitrary polarized target in irreducible spherical tensor representation. The general approach for the description of neutron spin dynamics for the propagation trough the medium with an arbitrary polarization is discussed in a relation to the search for time reversal invariance violation in neutron scattering.


title: Noncrossing structured additive multiple-output Bayesian quantile regression models

abstract: Quantile regression models are a powerful tool for studying different points of the conditional distribution of univariate response variables. Their multivariate counterpart extension though is not straightforward, starting with the definition of multivariate quantiles. We propose here a flexible Bayesian quantile regression model when the response variable is multivariate, where we are able to define a structured additive framework for all predictor variables. We build on previous ideas considering a directional approach to define the quantiles of a response variable with multiple-outputs and we define noncrossing quantiles in every directional quantile model. We define a Markov Chain Monte Carlo (MCMC) procedure for model estimation, where the noncrossing property is obtained considering a Gaussian process design to model the correlation between several quantile regression models. We illustrate the results of these models using two data sets: one on dimensions of inequality in the population, such as income and health; the second on scores of students in the Brazilian High School National Exam, considering three dimensions for the response variable.


title: Global Existence for the N Body Euler-Poisson System

abstract: In this paper we investigate the problem of multiple expanding Newtonian stars that interact via their gravitational effect on each other. It is clear physically that if two stars at rest are separated initially, and start expanding as well as moving according to the laws of Newtonian gravity, they may eventually collide. Thus, one can ask whether each star can be given an initial position and velocity such that they can keep expanding without touching. We show that even with gravitational interaction between the bodies, a large class of initial positions and velocities give global-in-time solutions to the N Body Euler-Poisson system. To do this we use a scaling mechanism present in the compressible Euler system shown in [63] and a careful analysis of how the gravitational interaction between stars affects their dynamics.


title: Single-Shot Characterization of High Transformer Ratio Wakefields in Nonlinear Plasma Acceleration

abstract: Plasma wakefields can enable very high accelerating gradients for frontier high energy particle accelerators, in excess of 10 GeV/m. To overcome limits on total acceleration achievable, specially shaped drive beams can be used in both linear and nonlinear plasma wakefield accelerators (PWFA), to increase the transformer ratio, implying that the drive beam deceleration is minimized relative to acceleration obtained in the wake. In this Letter, we report the results of a nonlinear PWFA, high transformer ratio experiment using high-charge, longitudinally asymmetric drive beams in a plasma cell. An emittance exchange process is used to generate variable drive current profiles, in conjunction with a long (multiple plasma wavelength) witness beam. The witness beam is energy-modulated by the wakefield, yielding a response that contains detailed spectral information in a single-shot measurement. Using these methods, we generate a variety of beam profiles and characterize the wakefields, directly observing beam-loaded transformer ratios up to R=7.8. Furthermore, a spectrally-based reconstruction technique, validated by 3D particle-in-cell simulations, is introduced to obtain the drive beam current profile from the decelerating wake data.


title: Wolf-Rayet stars in young massive star clusters as potential sources of Galactic cosmic rays

abstract: For most elements, the isotopic ratios seen in cosmic rays (CRs) are similar to those in the solar wind. The most important exception to this is $^{22}$Ne/$^{20}$Ne where the CR value is $\sim 5$ times that of the solar wind. According to most recent models of nucleosynthesis, a large amount of $^{22}\mathrm{Ne}$ is generated in Wolf-Rayet (WR) stars. In the winds of carbon sequence WR stars, i.e., WC stars, the isotopic ratio $^{22}\mathrm{Ne}$/$^{20}\mathrm{Ne}$ can be much larger than in the solar wind. Here, we consider CRs produced by $^{22}\mathrm{Ne}$-enriched WR winds in compact clusters of young massive stars like Westerlund~1. We assume that efficient CR acceleration in clusters occurs over the WR lifetime in an ensemble of shock waves from multiple massive star winds. We estimate the fraction of all Galactic CRs such sources may produce for a given set of parameters.


title: Magnetic field control of topological magnon-polaron bands in two-dimensional ferromagnets

abstract: We theoretically study magnon-phonon hybrid excitations in a square lattice ferromagnet subjected to a magnetic field by varying the field direction. The bulk bands of hybrid excitations, which are referred to as magnon-polarons, are investigated by considering all three phonon modes: vertical phonon, transverse phonon, and longitudinal phonon. We show that the topological proprieties of three hybridizations are different in terms of the Berry curvature and the Chern numbers. We also find that the topological properties of the bands can be controlled by changing the direction of the magnetic field, exhibiting one or more topological phase transitions. The dependence of thermal Hall conductivity as a function of magnetic field direction is proposed as an experiment probe of our theoretical results.


title: Digital Democracy: Episode IV -- A New Hope, How a Corporation for Public Software Could Transform Digital Engagement for Government and Civil Society

abstract: Though successive generations of digital technology have become increasingly powerful in the past twenty years, digital democracy has yet to realize its potential for deliberative transformation. The undemocratic exploitation of massive social media systems continued this trend, but it only worsened an existing problem of modern democracies, which were already struggling to develop deliberative infrastructure independent of digital technologies. There have been many creative conceptions of civic tech, but implementation has lagged behind innovation. This essay argues for implementing one such vision of digital democracy through the establishment of a public corporation. Modeled on the Corporation for Public Broadcasting in the U.S., this entity would foster the creation of new digital technology by providing a stable source of funding to nonprofit technologists, interest groups, civic organizations, government, researchers, private companies, and the public. Funded entities would produce and maintain software infrastructure for public benefit. The concluding sections identify what circumstances might create and sustain such an entity.


title: Machine learning Calabi-Yau metrics

abstract: We apply machine learning to the problem of finding numerical Calabi-Yau metrics. Building on Donaldson's algorithm for calculating balanced metrics on K\"ahler manifolds, we combine conventional curve fitting and machine-learning techniques to numerically approximate Ricci-flat metrics. We show that machine learning is able to predict the Calabi-Yau metric and quantities associated with it, such as its determinant, having seen only a small sample of training data. Using this in conjunction with a straightforward curve fitting routine, we demonstrate that it is possible to find highly accurate numerical metrics much more quickly than by using Donaldson's algorithm alone, with our new machine-learning algorithm decreasing the time required by between one and two orders of magnitude.


title: AI Safety for High Energy Physics

abstract: The field of high-energy physics (HEP), along with many scientific disciplines, is currently experiencing a dramatic influx of new methodologies powered by modern machine learning techniques. Over the last few years, a growing body of HEP literature has focused on identifying promising applications of deep learning in particular, and more recently these techniques are starting to be realized in an increasing number of experimental measurements. The overall conclusion from this impressive and extensive set of studies is that rarer and more complex physics signatures can be identified with the new set of powerful tools from deep learning. However, there is an unstudied systematic risk associated with combining the traditional HEP workflow and deep learning with high-dimensional data. In particular, calibrating and validating the response of deep neural networks is in general not experimentally feasible, and therefore current methods may be biased in ways that are not covered by current uncertainty estimates. By borrowing ideas from AI safety, we illustrate these potential issues and propose a method to bound the size of unaccounted for uncertainty. In addition to providing a pragmatic diagnostic, this work will hopefully begin a dialogue within the community about the robust application of deep learning to experimental analyses.


title: Exorcising Spectres with Secure Compilers

abstract: Attackers can access sensitive information of programs by exploiting the side-effects of speculatively-executed instructions using Spectre attacks. To mitigate theses attacks, popular compilers deployed a wide range of countermeasures. The security of these countermeasures, however, has not been ascertained: while some of them are believed to be secure, others are known to be insecure and result in vulnerable programs. To reason about the security guarantees of these compiler-inserted countermeasures, this paper presents a framework comprising several secure compilation criteria characterizing when compilers produce code resistant against Spectre attacks. With this framework, we perform a comprehensive security analysis of compiler-level countermeasures against Spectre attacks implemented in major compilers. This work provides sound foundations to formally reason about the security of compiler-level countermeasures against Spectre attacks as well as the first proofs of security and insecurity of said countermeasures.


title: Detecting Gravitational Scattering of Interstellar Objects Using Pulsar Timing

abstract: Gravitational scattering events, in which the path of an interstellar object is deflected by a pulsar or the solar system, give rise to reflex motion which can potentially be detected using pulsar timing. We determine the form of the timing signal expected from a gravitational scattering event, which is ramp-like and resembles the signal produced by a glitch or a gravitational wave burst with memory (BWM), and investigate the prospects for detecting such a signal using a pulsar timing array. The level of timing precision currently achieved for some millisecond pulsars makes it possible to detect objects as small as $10^{-10}$ solar masses, less than the mass of the dwarf planet Ceres, at impact parameters as large as 1 AU. The signals produced by gravitational scattering could provide independent constraints on models of dark matter involving asteroid-mass objects or subhalos, and should be considered as potential false positives in searches for BWMs.


title: Asymptotic Behavior of Polynomially Bounded Solutions of Linear Fractional Differential Equations

abstract: In this paper we study the asymptotic behavior of solutions of fractional differential equations of the form $D^{\alpha}_Cu(t)=Au(t)+f(t)$ on the half line, where $D^{\alpha}_Cu(t)$ is the derivative of the function $u$ in Caputo's sense, $A$ is generally an unbounded closed operator, $f$ is polynomially bounded. To this end we develop a spectral theory for functions of polynomial growth on the half line. Our main result claims that if $u$ is mild solution of the Cauchy problem such that $\lim_{h\downarrow 0} \sup_{t\ge 0} \| u(t+h)-u(t)\|/(1+t)^n=0$, and $\sup_{t\ge 0} \| u(t)\| /(1+t)^n <\infty$, then, $\lim_{t\to\infty} u(t)/(1+t)^n =0$ provided that the spectral set $\Sigma (A,\alpha )\cap i\R$ is countable, where $\Sigma (A,\alpha )$ is defined to be the set of complex numbers $\xi$ such that $\lambda^{\alpha -1} (\lambda^\alpha -A)^{-1}$ is analytic in a neighborhood of $\xi$, and $u$ satisfies some ergodic


title: Probing quantum spin liquids in equilibrium using the inverse spin Hall effect

abstract: We propose an experimental method utilizing a strongly spin-orbit coupled metal to quantum magnet bilayer that will probe quantum magnets lacking long range magnetic order, e.g., quantum spin liquids, via examination of the voltage noise spectrum in the metal layer. The bilayer is held in thermal and chemical equilibrium, and spin fluctuations arising across the single interface are converted into voltage fluctuations in the metal as a result of the inverse spin Hall effect. We elucidate the theoretical workings of the proposed bilayer system, and provide precise predictions for the frequency characteristics of the enhancement to the ac electrical resistance measured in the metal layer for three candidate quantum spin liquid models. Application to the Heisenberg spin-$1/2$ kagom{\'e} lattice model should allow for the extraction of any spinon gap present. A quantum spin liquid consisting of fermionic spinons coupled to a $U(1)$ gauge field should cause subdominant $\W^{4/3}$ scaling of the resistance of the coupled metal. Finally, if the magnet is well-captured by the Kitaev model in the gapless spin liquid phase, then the proposed bilayer can extract the two-flux gap which arises in spite of the gapless spectrum of the fermions. We therefore show that spectral analysis of the ac resistance in the metal in a single interface, equilibrium bilayer can test the relevance of a quantum spin liquid model to a given candidate material.


title: A multilevel analysis to systemic exposure: insights from local and system-wide information

abstract: In the aftermath of the financial crisis, the growing literature on financial networks has widely documented the predictive power of topological characteristics (e.g. degree centrality measures) to explain the systemic impact or systemic vulnerability of financial institutions. In this work, we show that considering alternative topological measures based on local sub-network environment improves our ability to identify systemic institutions. To provide empirical evidence, we apply a two-step procedure. First, we recover network communities (i.e. close-peer environment) on a spillover network of financial institutions. Second, we regress alternative measures of vulnerability on three levels of topological measures: the global level (i.e. firm topological characteristics computed over the whole system), local level (i.e. firm topological characteristics computed over the community) and aggregated level by averaging individual characteristics over the community. The sample includes $46$ financial institutions (banks, broker-dealers, insurance and real-estate companies) listed in the Standard \& Poor's 500 index. Our results confirm the informational content of topological metrics based on close-peer environment. Such information is different from the one embeds in traditional system wide topological metrics and is proved to be predictor of distress for financial institutions in time of crisis.


title: Trajectory Design for Energy Minimization in UAV-enabled Wireless Communications with Latency Constraints

abstract: This paper studies energy-efficient unmanned aerial vehicle (UAV)-enabled wireless communications, where the UAV acts as a flying base station (BS) to serve the ground users (GUs) within some predetermined latency constraints, e.g., requested timeout (RT). Our goal is to design the UAV trajectory to minimize the total energy consumption while satisfying the RT requirement and energy budget, which is accomplished via jointly optimizing the trajectory and UAV's velocities along subsequent hops. The corresponding optimization problem is difficult to solve due to its non-convexity and combinatorial nature. To overcome this difficulty, we solve the original problem via two consecutive steps. Firstly, we propose two algorithms, namely heuristic search, and dynamic programming (DP) to obtain a feasible set of trajectories without violating the GU's RT requirements based on the traveling salesman problem with time window (TSPTW). Then, they are compared with exhaustive search and traveling salesman problem (TSP) used as reference methods. While the exhaustive algorithm achieves the best performance at a high computation cost, the heuristic algorithm exhibits poorer performance with low complexity. As a result, the DP is proposed as a practical trade-off between the exhaustive and heuristic algorithms. Specifically, the DP algorithm results in near-optimal performance at a much lower complexity. Secondly, for given feasible trajectories, we propose an energy minimization problem via a joint optimization of the UAV's velocities along subsequent hops. Finally, numerical results are presented to demonstrate the effectiveness of our proposed algorithms. ...


title: Poisson CNN: Convolutional neural networks for the solution of the Poisson equation on a Cartesian mesh

abstract: The Poisson equation is commonly encountered in engineering, for instance in computational fluid dynamics (CFD) where it is needed to compute corrections to the pressure field to ensure the incompressibility of the velocity field. In the present work, we propose a novel fully convolutional neural network (CNN) architecture to infer the solution of the Poisson equation on a 2D Cartesian grid with different resolutions given the right hand side term, arbitrary boundary conditions and grid parameters. It provides unprecedented versatility for a CNN approach dealing with partial differential equations. The boundary conditions are handled using a novel approach by decomposing the original Poisson problem into a homogeneous Poisson problem plus four inhomogeneous Laplace sub-problems. The model is trained using a novel loss function approximating the continuous $L^p$ norm between the prediction and the target. Even when predicting on grids denser than previously encountered, our model demonstrates encouraging capacity to reproduce the correct solution profile. The proposed model, which outperforms well-known neural network models, can be included in a CFD solver to help with solving the Poisson equation. Analytical test cases indicate that our CNN architecture is capable of predicting the correct solution of a Poisson problem with mean percentage errors below 10%, an improvement by comparison to the first step of conventional iterative methods. Predictions from our model, used as the initial guess to iterative algorithms like Multigrid, can reduce the RMS error after a single iteration by more than 90% compared to a zero initial guess.


title: A geometric framework to compare classical field theories and to transfer solutions between PDEs

abstract: In this contribution, a mathematical framework is constructed to relate and compare non-linear partial differential equations (PDEs) in the category of smooth manifolds. In particular, it can be used to compare those aspects of field theories (e.g. of classical (Newtonian) mechanics, hydrodynamics, electrodynamics, relativity theory, classical Yang-Mills theory and so on) that are described by such equations. Employing a geometric (jet space) approach, a suitable notion of shared structure of two systems of PDEs is identified. It is proven that this shared structure can serve to transfer solutions from one theory to another and a generalization of so-called B\"acklund transformations is derived that can be used to generate non-trivial solutions of some non-linear PDEs. A procedure (based on formal integrability) is introduced with which one can explicitly compute the minimal consistency conditions that two systems of PDEs need to fulfill in order to share structure under a given correspondence. Furthermore, it is shown how symmetry groups can be used to identify useful correspondences and structure that is shared up to symmetries. Thereby, the role that B\"acklund transformations play in the theory of quotient equations is clarified. Explicit examples illustrate the general ideas throughout the text and in the last chapter, the framework is applied to systems related to electrodynamics and hydrodynamics.


title: Fitting a Kalman Smoother to Data

abstract: This paper considers the problem of fitting the parameters of a Kalman smoother to data. We formulate the Kalman smoothing problem with missing measurements as a constrained least squares problem and provide an efficient method to solve it based on sparse linear algebra. We then introduce the Kalman smoother tuning problem, which seeks to find parameters that achieve low prediction error on held out measurements. We derive a Kalman smoother auto-tuning algorithm, which is based on the proximal gradient method, that finds good, if not the best, parameters for a given dataset. Central to our method is the computation of the gradient of the prediction error with respect to the parameters of the Kalman smoother; we describe how to compute this at little to no additional cost. We demonstrate the method on population migration within the United States as well as data collected from an IMU+GPS system while driving. The paper is accompanied by an open-source implementation.


title: An exactly mass conserving space-time embedded-hybridized discontinuous Galerkin method for the Navier-Stokes equations on moving domains

abstract: This paper presents a space-time embedded-hybridized discontinuous Galerkin (EHDG) method for the Navier--Stokes equations on moving domains. This method uses a different hybridization compared to the space-time hybridized discontinuous Galerkin (HDG) method we presented previously in (Int. J. Numer. Meth. Fluids 89: 519--532, 2019). In the space-time EHDG method the velocity trace unknown is continuous while the pressure trace unknown is discontinuous across facets. In the space-time HDG method, all trace unknowns are discontinuous across facets. Alternatively, we present also a space-time embedded discontinuous Galerkin (EDG) method in which all trace unknowns are continuous across facets. The advantage of continuous trace unknowns is that the formulation has fewer global degrees-of-freedom for a given mesh than when using discontinuous trace unknowns. Nevertheless, the discrete velocity field obtained by the space-time EHDG and EDG methods, like the space-time HDG method, is exactly divergence-free, even on moving domains. However, only the space-time EHDG and HDG methods result in divergence-conforming velocity fields. An immediate consequence of this is that the space-time EHDG and HDG discretizations of the conservative form of the Navier--Stokes equations are energy stable. The space-time EDG method, on the other hand, requires a skew-symmetric formulation of the momentum advection term to be energy-stable. Numerical examples will demonstrate the differences in solution obtained by the space-time EHDG, EDG, and HDG methods.


title: A Bid-Validity Mechanism for Sequential Heat and Electricity Market Clearing

abstract: Coordinating the operation of units at the interface between heat and electricity systems, such as combined heat and power plants and heat pumps, is essential to reduce inefficiencies in each system and help achieve a cost-effective and efficient operation of the overall energy system. These energy systems are currently operated by sequential markets, which interface the technical and economic aspects of the systems. In that context, this study introduces an electricity-aware heat unit commitment model, which seeks to optimize the operation of the heat system while accounting for the techno-economic interdependencies between heat and electricity markets. These interdependencies are represented by bid-validity constraints, which model the linkage between the heat and electricity outputs and costs of combined heat and power plants and heat pumps. This approach also constitutes a novel market mechanism for the coordination of heat and electricity systems, which defines heat bids conditionally on electricity prices. Additionally, a tractable reformulation of the resulting trilevel optimization problem as a mixed integer linear program is proposed. Finally, it is shown on a case study that the proposed model yields a 23% reduction in total operating cost and a 6% reduction in wind curtailment compared to a traditional decoupled unit commitment model.


title: Electric-field control of spin transitions in molecular compounds

abstract: We present a theoretical model of spin transitions in stacks of molecular layers. Our model captures the already established physics of these systems (thermal hysteretic transitions and crossovers) and suggests a way towards in situ control of this physics by means of an external electric field. Our results pave the way toward both temperature and voltage controllable organic memory.


title: Atomically Controlled Tunable Doping in High Performance WSe2 Devices

abstract: Two-dimensional transitional metal dichalcogenide (TMD) field-effect transistors (FETs) are promising candidates for future electronic applications, owing to their excellent transport properties and potential for ultimate device scaling. However, it is widely acknowledged that substantial contact resistance associated with the contact-TMD interface has impeded device performance to a large extent. It has been discovered that O2 plasma treatment can convert WSe2 into WO3-x and substantially improve contact resistances of p-type WSe2 devices by strong doping induced thinner depletion width. In this paper, we carefully study the temperature dependence of this conversion, demonstrating an oxidation process with a precise monolayer control at room temperature and multilayer conversion at elevated temperatures. Furthermore, the lateral oxidation of WSe2 under the contact revealed by HR-STEM leads to potential unpinning of the metal Fermi level and Schottky barrier lowering, resulting in lower contact resistances. The p-doping effect is attributed to the high electron affinity of the formed WO3-x layer on top of the remaining WSe2 channel, and the doping level is found to be dependent on the WO3-x thickness that is controlled by the temperature. Comprehensive materials and electrical characterizations are presented, with a low contact resistance of ~528 ohm-um and record high on-state current of 320 uA/um at -1V bias being reported.


title: Observation of all pre- and post-reconnection structures involved in three-dimensional reconnection geometries in solar eruptions

abstract: We report on observations of the two newly-identified reconnection geometries involving erupting flux ropes. In 3D, a flux rope can reconnect either with a surrounding coronal arcade (recently named "ar-rf" reconnection) or with itself ("rr-rf" reconnection), and both kinds of reconnection create a new flux rope field line and a flare loop. For the first time, we identify all four constituents of both reconnections in a solar eruptive event, the filament eruption of 2011 June 07 observed by SDO/AIA. The ar-rf reconnection manifests itself as shift of one leg of the filament by more than 25" northward. At its previous location, a flare arcade is formed, while the new location of the filament leg previously corresponded to a footpoint of a coronal loop in 171 A. In addition, the evolution of the flare ribbon hooks is also consistent with the occurrence of ar--rf reconnection as predicted by MHD simulations. Specifically, the growing hook sweeps footpoints of preeruptive coronal arcades, and these locations become inside the hook. Furthermore, the rr-rf reconnection occurs during the peak phase above the flare arcade, in an apparently X-type geometry involving a pair of converging bright filament strands in the erupting filament. A new flare loop forms near the leg of one of the strands, while a bright blob, representing a remnant of the same strand, is seen ascending into the erupting filament. All together, these observations vindicate recent predictions of the 3D standard solar flare model.


title: Abelian group actions and hypersmooth equivalence relations

abstract: We show that a Borel action of a standard Borel group which is isomorphic to a sum of a countable abelian group with a countable sum of real lines and circles induces an orbit equivalence relation which is hypersmooth, i.e., Borel reducible to eventual agreement on sequences of reals, and it follows from this result along with the structure theory for locally compact abelian groups that Borel actions of Polish LCA groups induce orbit equivalence relations which are essentially hyperfinite, extending a result of Gao and Jackson and answering a question of Ding and Gao.


title: Probing primordial $^3$He from hyperfine line afterglows around supercritical black holes

abstract: We consider the possibility of the detection of $^3$HeII hyperfine line (rest frequency, $8.67 \, \rm GHz$) emission from ionized zones around accreting black holes (BHs) formed at high redshifts, $z=15\hbox{--}30$. We show that the brightness temperature in 8.67GHz line increases and reaches a peak value after the accretion onto the BH exhausts and HeIII recombines into HeII. This period of brightening last up to 40~million years. We find that during this period the maximum brightness temperature reaches $\simeq 0.2\hbox{--}0.5 \mu$K, depending on the epoch when such a black hole starts growing. The maximum angular size of the region emitting in the hyperfine line is around $0.5'$. The flux from such a region ($\simeq 0.3 \, \rm nJy$) is too small to be detected by SKA1-MID. The RMS of the collective flux from many emitting regions from a volume bounded by the synthesized beam and the band-width of SKA1-MID might reach 100~nJy, which is potentially detectable by SKA1-MID.


title: A Fast Saddle-Point Dynamical System Approach to Robust Deep Learning

abstract: Recent focus on robustness to adversarial attacks for deep neural networks produced a large variety of algorithms for training robust models. Most of the effective algorithms involve solving the min-max optimization problem for training robust models (min step) under worst-case attacks (max step). However, they often suffer from high computational cost from running several inner maximization iterations (to find an optimal attack) inside every outer minimization iteration. Therefore, it becomes difficult to readily apply such algorithms for moderate to large size real world data sets. To alleviate this, we explore the effectiveness of iterative descent-ascent algorithms where the maximization and minimization steps are executed in an alternate fashion to simultaneously obtain the worst-case attack and the corresponding robust model. Specifically, we propose a novel discrete-time dynamical system-based algorithm that aims to find the saddle point of a min-max optimization problem in the presence of uncertainties. Under the assumptions that the cost function is convex and uncertainties enter concavely in the robust learning problem, we analytically show that our algorithm converges asymptotically to the robust optimal solution under a general adversarial budget constraints as induced by $\ell_p$ norm, for $1\leq p\leq \infty$. Based on our proposed analysis, we devise a fast robust training algorithm for deep neural networks. Although such training involves highly non-convex robust optimization problems, empirical results show that the algorithm can achieve significant robustness compared to other state-of-the-art robust models on benchmark data sets.


title: A lattice model of heavy-light three-body system

abstract: We present a study of a $1+1$ dimensional heavy-light three-body system in finite volume. The heavy-light system is simulated by a coupled-channel $\phi^4$ type lattice model, and both ground state and excited states of multiparticle energy spectra are measured on various lattices. The lattice simulation data analysis is performed based on variational approach.


title: Comparing Greedy Constructive Heuristic Subtour Elimination Methods for the Traveling Salesman Problem

abstract: This paper further defines the class of fragment constructive heuristics used to compute feasible solutions for the Traveling Salesman Problem into arc-greedy and node-greedy subclasses. Since these subclasses of heuristics can create subtours, two known methodologies for subtour elimination on symmetric instances are reviewed and are expanded to cover asymmetric problem instances. This paper introduces a third novel methodology, the Greedy Tracker, and compares it to both known methodologies. Computational results are generated across multiple symmetric and asymmetric instances. The results demonstrate the Greedy Tracker is the fastest method for preventing subtours for instances below 400 nodes. A distinction between fragment constructive heuristics and the subtour elimination methodology used to ensure the feasibility of resulting solutions enables the introduction of a new node-greedy fragment heuristic called Ordered Greedy.


title: Collection of Historical Weather Data: Issues with Missing Values

abstract: Weather data collected from automated weather stations have become a crucial component for making decisions in agriculture and in forestry. Over time, weather stations may become out-of-order or stopped for maintenance, and therefore, during those periods, the data values will be missing. Unfortunately, this will cause huge problems when analysing the data. The main aim of this study is to create high-quality historical weather datasets by dealing efficiently with missing values. In this paper, we present a set of missing data imputation methods and study their effectiveness. These methods were used based on different types of missing values. The experimental results show that two the proposed methods are very promising and can be used at larger scale.


title: On the quantum behavior and clustering properties of correlated financial portfolios

abstract: We investigate 17 digital currencies making an analogy with quantum systems and develop the concept of eigenportfolios. We show that the density of states of the correlation matrix of these assets shows a behavior between that of the Wishart ensemble and one whose elements are Cauchy distributed. A metric for the participation matrix based on superposition of Gaussian functions is proposed and we show that small eigenvalues correspond to localized states. Nonetheless, some level of localization is also present for bigger eigenvalues probably caused by the fat tails of the distribution of returns of these assets. We also show through a clustering study that the digital currencies tend to stagger together. We conclude the paper showing that this correlation structure leads to an Epps effect.


title: Sector Neutral Portfolios: Long memory motifs persistence in market structure dynamics

abstract: We study soft persistence (existence in subsequent temporal layers of motifs from the initial layer) of motif structures in Triangulated Maximally Filtered Graphs (TMFG) generated from time-varying Kendall correlation matrices computed from stock prices log-returns over rolling windows with exponential smoothing. We observe long-memory processes in these structures in the form of power law decays in the number of persistent motifs. The decays then transition to a plateau regime with a power-law decay with smaller exponent. We demonstrate that identifying persistent motifs allows for forecasting and applications to portfolio diversification. Balanced portfolios are often constructed from the analysis of historic correlations, however not all past correlations are persistently reflected into the future. Sector neutrality has also been a central theme in portfolio diversification and systemic risk. We present an unsupervised technique to identify persistently correlated sets of stocks. These are empirically found to identify sectors driven by strong fundamentals. Applications of these findings are tested in two distinct ways on four different markets, resulting in significant reduction in portfolio volatility. A persistence-based measure for portfolio allocation is proposed and shown to outperform volatility weighting when tested out of sample.


title: Neural Logic Networks

abstract: Recent years have witnessed the great success of deep neural networks in many research areas. The fundamental idea behind the design of most neural networks is to learn similarity patterns from data for prediction and inference, which lacks the ability of logical reasoning. However, the concrete ability of logical reasoning is critical to many theoretical and practical problems. In this paper, we propose Neural Logic Network (NLN), which is a dynamic neural architecture that builds the computational graph according to input logical expressions. It learns basic logical operations as neural modules, and conducts propositional logical reasoning through the network for inference. Experiments on simulated data show that NLN achieves significant performance on solving logical equations. Further experiments on real-world data show that NLN significantly outperforms state-of-the-art models on collaborative filtering and personalized recommendation tasks.


title: Betatron frequency and the Poincare rotation number

abstract: Symplectic maps are routinely used to describe single-particle dynamics in circular accelerators. In the case of a linear accelerator map, the rotation number (the betatron frequency) can be easily calculated from the map itself. In the case of a nonlinear map, the rotation number is normally obtained numerically, by iterating the map for given initial conditions, or through a normal form analysis, a type of a perturbation theory for maps. Integrable maps, a subclass of symplectic maps, allow for an analytic evaluation of their rotation numbers. In this paper we propose an analytic expression to determine the rotation number for integrable symplectic maps of the plane and present several examples, relevant to accelerators. These new results can be used to analyze the topology of the accelerator Hamiltonians as well as to serve as the starting point for a perturbation theory for maps.


title: Extension of a residually finite group by a residually finite group is weakly sofic

abstract: We show that residually finite by residually finite extensions are weakly sofic.


title: 28 GHz mmWave Channel Measurements and Modeling in a Library Environment

abstract: To fully exploit the millimeter-wave bands for the fifth generation cellular systems, an accurate understanding of the channel propagation characteristics is required, and hence extensive measurement campaigns in different environments are needed. In this paper, we use a rotated directional antenna-based channel sounder for measurements at 28 GHz in large indoor environments at a library setting. We present models for power angular-delay profile and large-scale path loss based on the measurements over distances ranging from 10 m to 50 m. In total, nineteen different line-of-sight (LOS) and non-line-of-sight (NLOS) scenarios are considered, including the cases where the transmitter and the receiver are placed on different floors. Results show that the close-in free space reference distance and the floating intercept path loss models both perform well in fitting the empirical data. The path loss exponent obtained for the LOS scenarios is found to be very close to that of the free space path loss model.


title: Focusing on the dynamics of the entanglement in spin junction

abstract: We study the dynamics of entanglement in the one-dimensional spin-1/2 XY model in the presence of a transverse magnetic field. A pair of spins are considered as an open quantum system, while the rest of the chain plays the role of the environment. Our study focuses on the pair of spins in the system, the edge spins, and the environment. It is observed that the entanglement between the pair of spins in the system decreases and it can transfer to the rest of the spins. For a value of anisotropy leading to the Ising model, the entanglement is completely back to the system by passing time. On the other hand, the entanglement can only be seen under certain conditions between edge spins of the system and the environment. The pair of spins on the edge will be entangled very quickly and it will disappear after a very short time. A pair of spins far from the system was chosen to examine the behavior of entanglement in the environment. As expected, the transmission of entanglement from the system to the environment takes notable time.


title: Universal Composability is Robust Compilation

abstract: This paper discusses the relationship between two frameworks: universal composability (UC) and robust compilation (RC). In cryptography, UC is a framework for the specification and analysis of cryptographic protocols with a strong compositionality guarantee: UC protocols remain secure even when composed with other protocols. In programming language security, RC is a novel framework for determining secure compilation by proving whether compiled programs are as secure as their source-level counterparts no matter what target-level code they interact with. Presently, these disciplines are studied in isolation, though we argue that there is a deep connection between them and exploring this connection will benefit both research fields. This paper formally proves the connection between UC and RC and then it explores the benefits of this connection. For this, this paper first identifies which conditions must programming languages fulfil in order to possibly attain UC-like composition. Then, it proves UC of both an existing and a new commitment protocol as a corollary of the related compilers attaining RC. Finally, it mechanises these proofs in Deepsec, obtaining symbolic guarantees that the protocol is indeed UC. Our connection lays the groundwork towards a better and deeper understanding of both UC and RC, and the benefits we showcase from this connection provide first evidence of scalable mechanised proofs for UC.


title: Tree-based Intelligent Intrusion Detection System in Internet of Vehicles

abstract: The use of autonomous vehicles (AVs) is a promising technology in Intelligent Transportation Systems (ITSs) to improve safety and driving efficiency. Vehicle-to-everything (V2X) technology enables communication among vehicles and other infrastructures. However, AVs and Internet of Vehicles (IoV) are vulnerable to different types of cyber-attacks such as denial of service, spoofing, and sniffing attacks. In this paper, an intelligent intrusion detection system (IDS) is proposed based on tree-structure machine learning models. The results from the implementation of the proposed intrusion detection system on standard data sets indicate that the system has the ability to identify various cyber-attacks in the AV networks. Furthermore, the proposed ensemble learning and feature selection approaches enable the proposed system to achieve high detection rate and low computational cost simultaneously.


title: The TCGA Meta-Dataset Clinical Benchmark

abstract: Machine learning is bringing a paradigm shift to healthcare by changing the process of disease diagnosis and prognosis in clinics and hospitals. This development equips doctors and medical staff with tools to evaluate their hypotheses and hence make more precise decisions. Although most current research in the literature seeks to develop techniques and methods for predicting one particular clinical outcome, this approach is far from the reality of clinical decision making in which you have to consider several factors simultaneously. In addition, it is difficult to follow the recent progress concretely as there is a lack of consistency in benchmark datasets and task definitions in the field of Genomics. To address the aforementioned issues, we provide a clinical Meta-Dataset derived from the publicly available data hub called The Cancer Genome Atlas Program (TCGA) that contains 174 tasks. We believe those tasks could be good proxy tasks to develop methods which can work on a few samples of gene expression data. Also, learning to predict multiple clinical variables using gene-expression data is an important task due to the variety of phenotypes in clinical problems and lack of samples for some of the rare variables. The defined tasks cover a wide range of clinical problems including predicting tumor tissue site, white cell count, histological type, family history of cancer, gender, and many others which we explain later in the paper. Each task represents an independent dataset. We use regression and neural network baselines for all the tasks using only 150 samples and compare their performance.


title: Markovian and Non-Markovian dynamics in the one-dimensional transverse-field XY model

abstract: We consider an anisotropic spin-1/2 XY Heisenberg chain in the presence of a transverse magnetic field. Selecting the nearest neighbor pair spins as an open quantum system, the rest of the chain plays the role of the structured environment. In fact, the aforementioned system is used as a quantum probe signifying nontrivial features of the environment with which is interacting. We use a general measure that is based on the trace distance for the degree of non-Markovian behavior in open quantum systems. The witness of non-Markovianity takes on nonzero values whenever there is a flow of information from the environment back to the open system. We have shown that the dynamics of the system with isotropic Heisenberg interaction is Markovian. A dynamical transition into the non-Markovian regime is observed as soon as the anisotropy, {\gamma}, is applied. At the Ising value of the anisotropy {\gamma} = 1.0, all the information flows back from the environment to the system. The additional dynamical transition from the non-Markovian to the Markovian is obtained by applying the transverse magnetic field. In addition, we have focused on the time evolution of the Loschmidtecho return rate function. It is found that a non-analyticity can be seen in the time evolution of the Loschmidt-echo return rate function exactly at the critical points where a dynamical transition from the Markovian to the non-Markovian occurs.


title: Extended Search for the Invisible Axion with the Axion Dark Matter Experiment

abstract: This paper reports on a cavity haloscope search for dark matter axions in the galactic halo in the mass range $2.81$-$3.31$ ${\mu}eV$. This search excludes the full range of axion-photon coupling values predicted in benchmark models of the invisible axion that solve the strong CP problem of quantum chromodynamics, and marks the first time a haloscope search has been able to search for axions at mode crossings using an alternate cavity configuration. Unprecedented sensitivity in this higher mass range is achieved by deploying an ultra low-noise Josephson parametric amplifier as the first stage signal amplifier.


title: OffWorld Gym: open-access physical robotics environment for real-world reinforcement learning benchmark and research

abstract: Success stories of applied machine learning can be traced back to the datasets and environments that were put forward as challenges for the community. The challenge that the community sets as a benchmark is usually the challenge that the community eventually solves. The ultimate challenge of reinforcement learning research is to train real agents to operate in the real environment, but until now there has not been a common real-world RL benchmark. In this work, we present a prototype real-world environment from OffWorld Gym -- a collection of real-world environments for reinforcement learning in robotics with free public remote access. Close integration into existing ecosystem allows the community to start using OffWorld Gym without any prior experience in robotics and takes away the burden of managing a physical robotics system, abstracting it under a familiar API. We introduce a navigation task, where a robot has to reach a visual beacon on an uneven terrain using only the camera input and provide baseline results in both the real environment and the simulated replica. To start training, visit https://gym.offworld.ai


title: Are Perceptually-Aligned Gradients a General Property of Robust Classifiers?

abstract: For a standard convolutional neural network, optimizing over the input pixels to maximize the score of some target class will generally produce a grainy-looking version of the original image. However, Santurkar et al. (2019) demonstrated that for adversarially-trained neural networks, this optimization produces images that uncannily resemble the target class. In this paper, we show that these "perceptually-aligned gradients" also occur under randomized smoothing, an alternative means of constructing adversarially-robust classifiers. Our finding supports the hypothesis that perceptually-aligned gradients may be a general property of robust classifiers. We hope that our results will inspire research aimed at explaining this link between perceptually-aligned gradients and adversarial robustness.


title: Nonhedgeable risk and Credit Risk Pricing

abstract: We introduce a new model for pricing corporate bonds, which is a modification of the classical model of Merton. In this new model, we drop the liquidity assumption of the firm's asset value process, and assume that there is a liquidly traded asset in the market whose value is correlated with the firm's asset value, and all portfolios can be constructed using solely this asset and the money market account. We formulate the market price of the corporate bond as the product of the price of an optimal replicating portfolio and exp(- kappa x replication error), where kappa is a positive constant. The interpretation is that the representative investor accepts the price of the optimal replicating portfolio as a benchmark, however, requests compensation for the non-hedgeable risk. We show that if the replication error is measured relative to the firm's value, the resulting formula is arbitrage free with mild restrictions on the parameters.


title: Cocycle superrigidity for profinite actions of irreducible lattices

abstract: Let $\Gamma$ be an irreducible lattice in a product of two locally compact groups and assume that $\Gamma$ is densely embedded in a profinite group $K$. We give necessary conditions which imply that the left translation action $\Gamma\curvearrowright K$ is "virtually" cocycle superrigid: any cocycle $w:\Gamma\times K\rightarrow\Delta$ with values in a countable group $\Delta$ is cohomologous to a cocycle which factors through the map $\Gamma\times K\rightarrow\Gamma\times K_0$, for some finite quotient group $K_0$ of $K$. As a corollary, we deduce that any ergodic profinite action of $\Gamma=\text{SL}_2(\mathbb Z[S^{-1}])$ is virtually cocycle superrigid and virtually W$^*$-superrigid, for any finite nonempty set of primes $S$.


title: Intracranial Hemorrhage Segmentation Using Deep Convolutional Model

abstract: Traumatic brain injuries could cause intracranial hemorrhage (ICH). ICH could lead to disability or death if it is not accurately diagnosed and treated in a time-sensitive procedure. The current clinical protocol to diagnose ICH is examining Computerized Tomography (CT) scans by radiologists to detect ICH and localize its regions. However, this process relies heavily on the availability of an experienced radiologist. In this paper, we designed a study protocol to collect a dataset of 82 CT scans of subjects with traumatic brain injury. Later, the ICH regions were manually delineated in each slice by a consensus decision of two radiologists. Recently, fully convolutional networks (FCN) have shown to be successful in medical image segmentation. We developed a deep FCN, called U-Net, to segment the ICH regions from the CT scans in a fully automated manner. The method achieved a Dice coefficient of 0.31 for the ICH segmentation based on 5-fold cross-validation. The dataset is publicly available online at PhysioNet repository for future analysis and comparison.


title: Initialization methods for optimum average silhouette width clustering

abstract: A unified clustering approach that can estimate number of clusters and produce clustering against this number simultaneously is proposed. Average silhouette width (ASW) is a widely used standard cluster quality index. A distance based objective function that optimizes ASW for clustering is defined. The proposed algorithm named as OSil, only, needs data observations as an input without any prior knowledge of the number of clusters. This work is about thorough investigation of the proposed methodology, its usefulness and limitations. A vast spectrum of clustering structures were generated, and several well-known clustering methods including partitioning, hierarchical, density based, and spatial methods were consider as the competitor of the proposed methodology. Simulation reveals that OSil algorithm has shown superior performance in terms of clustering quality than all clustering methods included in the study. OSil can find well separated, compact clusters and have shown better performance for the estimation of number of clusters as compared to several methods. Apart from the proposal of the new methodology and it's investigation the paper offers a systematic analysis on the estimation of cluster indices, some of which never appeared together in comparative simulation setup before. The study offers many insightful findings useful for the selection of the clustering methods and indices for clustering quality assessment.


title: Conserved currents for electromagnetic fields in the Kerr spacetime

abstract: We construct a variety of conserved currents for test electromagnetic fields on a Kerr background. Our procedure, which involves the symplectic product for electromagnetism and symmetry operators, generates the conserved currents given by Andersson, B\"ackdahl, and Blue (2015), as well as a new conserved current. These currents reduce to the sum of (positive powers of) the Carter constants of photons in the geometric optics limit, and generalize the current for scalar fields discovered by Carter (1977) involving the Killing tensor. We furthermore show that the fluxes of our new current through null infinity and the horizon are finite.


title: The Bitwise Hashing Trick for Personalized Search

abstract: Many real world problems require fast and efficient lexical comparison of large numbers of short text strings. Search personalization is one such domain. We introduce the use of feature bit vectors using the hashing trick for improving relevance in personalized search and other personalization applications. We present results of several lexical hashing and comparison methods. These methods are applied to a user's historical behavior and are used to predict future behavior. Using a single bit per dimension instead of floating point results in an order of magnitude decrease in data structure size, while preserving or even improving quality. We use real data to simulate a search personalization task. A simple method for combining bit vectors demonstrates an order of magnitude improvement in compute time on the task with only a small decrease in accuracy.


title: Blameworthiness in Security Games

abstract: Security games are an example of a successful real-world application of game theory. The paper defines blameworthiness of the defender and the attacker in security games using the principle of alternative possibilities and provides a sound and complete logical system for reasoning about blameworthiness in such games. Two of the axioms of this system capture the asymmetry of information in security games.


title: n-m-Variant Systems: Adversarial-Resistant Software Rejuvenation for Cloud-Based Web Applications

abstract: Web servers are a popular target for adversaries as they are publicly accessible and often vulnerable to compromise. Compromises can go unnoticed for months, if not years, and recovery often involves a complete system rebuild. In this paper, we propose n-m-Variant Systems, an adversarial-resistant software rejuvenation framework for cloud-based web applications. We improve the state-of-the-art by introducing a variable m that provides a knob for administrators to tune an environment to balance resource usage, performance overhead, and security guarantees. Using m, security guarantees can be tuned for seconds, minutes, days, or complete resistance. We design and implement an n-m-Variant System prototype to protect a Mediawiki PHP application serving dynamic content from an external SQL persistent storage. Our performance evaluation shows a throughput reduction of 65% for 108 seconds of resistance and 83% for 12 days of resistance to sophisticated adversaries, given appropriate resource allocation. Furthermore, we use theoretical analysis and simulation to characterize the impact of system parameters on resilience to adversaries. Through these efforts, our work demonstrates how properties of cloud-based servers can enhance the integrity of Web servers.


title: Poisson stochastic master equation unravellings and the measurement problem: a quantum stochastic calculus perspective

abstract: The paper studies a class of quantum stochastic differential equations, modeling an interaction of a system with its environment in the quantum noise approximation. The space representing quantum noise is the symmetric Fock space over L^2(R_+). Using the isomorphism of this space with the space of square-integrable functionals of the Poisson process, the equations can be represented as classical stochastic differential equations, driven by Poisson processes. This leads to a discontinuous dynamical state reduction which we compare to the Ghirardi-Rimini-Weber model. A purely quantum object, the norm process, is found which plays the role of an observer (in the sense of Everett [H. Everett III, Reviews of modern physics, 29.3, 454, (1957)]), encoding all events occurring in the system space. An algorithm introduced by Dalibard et al [J. Dalibard, Y. Castin, and K. M{\o}lmer, Physical review letters, 68.5, 580 (1992)] to numerically solve quantum master equations is interpreted in the context of unravellings and the trajectories of expected values of system observables are calculated.


title: Toward Metrics for Differentiating Out-of-Distribution Sets

abstract: Vanilla CNNs, as uncalibrated classifiers, suffer from classifying out-of-distribution (OOD) samples nearly as confidently as in-distribution samples. To tackle this challenge, some recent works have demonstrated the gains of leveraging available OOD sets for training end-to-end calibrated CNNs. However, a critical question remains unanswered in these works: how to differentiate OOD sets for selecting the most effective one(s) that induce training such CNNs with high detection rates on unseen OOD sets? To address this pivotal question, we provide a criterion based on generalization errors of Augmented-CNN, a vanilla CNN with an added extra class employed for rejection, on in-distribution and unseen OOD sets. However, selecting the most effective OOD set by directly optimizing this criterion incurs a huge computational cost. Instead, we propose three novel computationally-efficient metrics for differentiating between OOD sets according to their "protection" level of in-distribution sub-manifolds. We empirically verify that the most protective OOD sets -- selected according to our metrics -- lead to A-CNNs with significantly lower generalization errors than the A-CNNs trained on the least protective ones. We also empirically show the effectiveness of a protective OOD set for training well-generalized confidence-calibrated vanilla CNNs. These results confirm that 1) all OOD sets are not equally effective for training well-performing end-to-end models (i.e., A-CNNs and calibrated CNNs) for OOD detection tasks and 2) the protection level of OOD sets is a viable factor for recognizing the most effective one. Finally, across the image classification tasks, we exhibit A-CNN trained on the most protective OOD set can also detect black-box FGS adversarial examples as their distance (measured by our metrics) is becoming larger from the protected sub-manifolds.


title: Hypergeometric Series Representations of Feynman Integrals by GKZ Hypergeometric Systems

abstract: We show that almost all Feynman integrals as well as their coefficients in a Laurent series in dimensional regularization can be written in terms of Horn hypergeometric functions. By applying the results of Gelfand-Kapranov-Zelevinsky (GKZ) we derive a formula for a class of hypergeometric series representations of Feynman integrals, which can be obtained by triangulations of the Newton polytope $\Delta_G$ corresponding to the Lee-Pomeransky polynomial $G$. Those series can be of higher dimension, but converge fast for convenient kinematics, which also allows numerical applications. Further, we discuss possible difficulties which can arise in a practical usage of this approach and give strategies to solve them.


title: On the shift-invert Lanczos method for the buckling eigenvalue problem

abstract: We consider the problem of extracting a few desired eigenpairs of the buckling eigenvalue problem $Kx = \lambda K_Gx$, where $K$ is symmetric positive semi-definite, $K_G$ is symmetric indefinite, and the pencil $K - \lambda K_G$ is singular, namely, $K$ and $K_G$ share a non-trivial common nullspace. Moreover, in practical buckling analysis of structures, bases for the nullspace of $K$ and the common nullspace of $K$ and $K_G$ are available. There are two open issues for developing an industrial strength shift-invert Lanczos method: (1) the shift-invert operator $(K - \sigma K_G)^{-1}$ does not exist or is extremely ill-conditioned, and (2) the use of the semi-inner product induced by $K$ drives the Lanczos vectors rapidly towards the nullspace of $K$, which leads to a rapid growth of the Lanczos vectors in norms and cause permanent loss of information and the failure of the method. In this paper, we address these two issues by proposing a generalized buckling spectral transformation of the singular pencil $K - \lambda K_G$ and a regularization of the inner product via a low-rank updating of the semi-positive definiteness of $K$. The efficacy of our approach is demonstrated by numerical examples, including one from industrial buckling analysis.


title: Link-homotopy classes of 4-component links and claspers

abstract: Two links are link-homotopic if they are transformed into each other by a sequence of self-crossing changes and ambient isotopies. The link-homotopy classes of 4-component links were classified by Levine with enormous algebraic computations. We modify the results by using Habiro's clasper theory. The new classification gives more symmetrical and schematic points of view to the link-homotopy classes of 4-component links. As applications, we give several new subsets of the link-homotopy classes of 4-component links which are classified by comparable invariants and give an algorithm which determines whether given two links are link-homotopic or not.


title: PyTorchPipe: a framework for rapid prototyping of pipelines combining language and vision

abstract: Access to vast amounts of data along with affordable computational power stimulated the reincarnation of neural networks. The progress could not be achieved without adequate software tools, lowering the entry bar for the next generations of researchers and developers. The paper introduces PyTorchPipe (PTP), a framework built on top of PyTorch. Answering the recent needs and trends in machine learning, PTP facilitates building and training of complex, multi-modal models combining language and vision (but is not limited to those two modalities). At its core, PTP employs a component-oriented approach and relies on the concept of a pipeline, defined as a directed acyclic graph of loosely coupled components. A user defines a pipeline using yaml-based (thus human-readable) configuration files, whereas PTP provides generic workers for their loading, training, and testing using all the computational power (CPUs and GPUs) that is available to the user. The paper covers the main concepts of PyTorchPipe, discusses its key features and briefly presents the currently implemented tasks, models and components.