title: Towards Models for Availability and Security Evaluation of Cloud Computing with Moving Target Defense

abstract: Security is one of the most relevant concerns in cloud computing. With the evolution of cyber-security threats, developing innovative techniques to thwart attacks is of utmost importance. One recent method to improve cloud computing security is Moving Target Defense (MTD). MTD makes use of dynamic reconfiguration in virtualized environments to "confuse" attackers or to nullify their knowledge about the system state. However, there is still no consolidated mechanism to evaluate the trade-offs between availability and security when using MTD on cloud computing. The evaluation through measurements is complex as one needs to deal with unexpected events as failures and attacks. To overcome this challenge, we intend to propose a set of models to evaluate the availability and security of MTD in cloud computing environments. The expected results include the quantification of availability and security levels under different conditions (e.g., different software aging rates, varying workloads, different attack intensities).


title: Effects of spatial dispersion on Self--induced transparency in two--level media

abstract: We study the effects of dispersion in carrier waves on the properties of soliton self--induced transparency (SIT) in two level media. We found substantial impact of dispersion effects on typical SIT soliton features. For example, the degree of SIT pulse velocity slowing down (acceleration) is determined by the ratio of the incoming pulse frequency over atomic transition frequency - $x=\omega/\omega_0$. Specifically, an immediate pulse stopping is predicted for absorbing media when pulse duration time exceeds some critical value. In the sharp line limit stopping may emerge only for frequency ratio above unity, while for the inhomogeneously broadened systems it appears irrespective of the value of $x$. Analysis performed on the basis of Mcall\& Hahn \textit{Area theorem} implies that pulse stopping is achieved when Ber's absorption coefficient approaches infinity, that is, pulse energy is fully absorbed in the medium. In the case of amplifying media super-luminal motion is predicted as in the case of resonance. However, there is a lowest value in the frequency ratio below which the pulse velocity tends to the sub-luminal region. These new features of the SIT phenomenon open novel ways on how it may be exploited for the control of electromagnetic wave radiation in two-level media. This may be achieved by varying frequency ratio.


title: A Novel Loss Function Incorporating Imaging Acquisition Physics for PET Attenuation Map Generation using Deep Learning

abstract: In PET/CT imaging, CT is used for PET attenuation correction (AC). Mismatch between CT and PET due to patient body motion results in AC artifacts. In addition, artifact caused by metal, beam-hardening and count-starving in CT itself also introduces inaccurate AC for PET. Maximum likelihood reconstruction of activity and attenuation (MLAA) was proposed to solve those issues by simultaneously reconstructing tracer activity ($\lambda$-MLAA) and attenuation map ($\mu$-MLAA) based on the PET raw data only. However, $\mu$-MLAA suffers from high noise and $\lambda$-MLAA suffers from large bias as compared to the reconstruction using the CT-based attenuation map ($\mu$-CT). Recently, a convolutional neural network (CNN) was applied to predict the CT attenuation map ($\mu$-CNN) from $\lambda$-MLAA and $\mu$-MLAA, in which an image-domain loss (IM-loss) function between the $\mu$-CNN and the ground truth $\mu$-CT was used. However, IM-loss does not directly measure the AC errors according to the PET attenuation physics, where the line-integral projection of the attenuation map ($\mu$) along the path of the two annihilation events, instead of the $\mu$ itself, is used for AC. Therefore, a network trained with the IM-loss may yield suboptimal performance in the $\mu$ generation. Here, we propose a novel line-integral projection loss (LIP-loss) function that incorporates the PET attenuation physics for $\mu$ generation. Eighty training and twenty testing datasets of whole-body 18F-FDG PET and paired ground truth $\mu$-CT were used. Quantitative evaluations showed that the model trained with the additional LIP-loss was able to significantly outperform the model trained solely based on the IM-loss function.


title: Modelling the AGN broad line region using single-epoch spectra II. Nearby AGN

abstract: The structure of the broad line region (BLR) is an essential ingredient in the determination of active galactic nuclei (AGN) virial black hole masses, which in turn are important to study the role of black holes in galaxy evolution. Constraints on the BLR geometry and dynamics can be obtained from velocity-resolved studies using reverberation mapping data (i.e. monitoring data). However, monitoring data are observationally expensive and only available for a limited sample of AGN, mostly confined to the local Universe. Here we explore a new version of a Bayesian inference, physical model of the BLR which uses an individual spectrum and prior information on the BLR size from the radius-luminosity relation, to model the AGN BLR geometry and dynamics. We apply our model to a sample of 11 AGN, which have been previously modelled using monitoring data. Our single-epoch BLR model is able to constrain some of the BLR parameters with inferred parameter values that agree within the uncertainties with those determined from the modelling of monitoring data. We find that our model is able to derive stronger constraints on the BLR for AGN with broad emission lines that qualitatively have more substructure and more asymmetry, presumably as they contain more information to constrain the physical model. The performance of this model makes it a practical and cost-effective tool to determine some of the BLR properties of a large sample of low and high redshift AGN, for which monitoring data are not available.


title: Streaming controlled by meniscus shape

abstract: Surface waves called meniscus waves often appear in the systems that are close to the capillary length scale. Since the meniscus shape determines the form of the meniscus waves, the resulting streaming circulation has a structure distinct from that caused by other capillary-gravity waves recently reported in the literature. In the present study, we produce symmetric and antisymmetric meniscus shapes by controlling boundary wettability and excite meniscus waves by oscillating the meniscus vertically. The symmetric and antisymmetric configurations produce different surface capillary-gravity wave modes and streaming flow structures. The energy density of the streaming circulation increases at the rate of the fourth power of the forcing amplitude in both configurations. The flow symmetry of streaming circulation is retained under the symmetric meniscus, while it is lost under the antisymmetric meniscus. In our experiments, the streaming circulation primarily originates from the Stokes boundary layer beneath the meniscus and can be successfully explained using the existing streaming theory.


title: Global Optima is not Limit Computable

abstract: We study the limit computability of finding a global optimum of a continuous function. We give a short proof to show that the problem of checking whether a point is a global minimum is not limit computable. Thereby showing the same for the problem of finding a global minimum. In the next part, we give an algorithm that converges to the global minima when a lower bound on the size of the basin of attraction of the global minima is known. We prove the convergence of this algorithm and provide some numerical experiments.


title: Probing thermalization through spectral analysis with matrix product operators

abstract: We combine matrix product operator techniques with Chebyshev polynomial expansions and present a method that is able to explore spectral properties of quantum many-body Hamiltonians. In particular, we show how this method can be used to probe thermalization of large spin chains without explicitly simulating their time evolution, as well as to compute full and local densities of states. The performance is illustrated with the examples of the Ising and PXP spin chains. For the non-integrable Ising chain, our findings corroborate the presence of thermalization for several initial states, well beyond what direct time-dependent simulations have been able to achieve so far.


title: Uniqueness of solution of an inverse source problem for ultrahyperbolic equations

abstract: The aim of this article is to investigate the uniqueness of solution of an inverse problem for ultrahyperbolic equations. We first reduce the inverse problem to a Cauchy problem for an integro-differential equation and then by using a pointwise Carleman type inequality we prove the uniqueness.


title: Quasar X-ray and UV flux, baryon acoustic oscillation, and Hubble parameter measurement constraints on cosmological model parameters

abstract: We use the Risaliti & Lusso (2015) compilation of 808 X-ray and UV flux measurements of quasars (QSOs) in the redshift range $0.061 \leq z \leq 6.28$, alone and in conjuction with baryon acoustic oscillation (BAO) and Hubble parameter [$H(z)$] measurements, to constrain cosmological parameters in six cosmological models. The QSO data constraints are significantly weaker than, but consistent with, those from the $H(z)$ + BAO data. A joint analysis of the QSO + $H(z)$ + BAO data is consistent with the current standard model, spatially-flat $\Lambda$CDM, but mildly favors closed spatial hypersurfaces and dynamical dark energy.


title: Brain2Char: A Deep Architecture for Decoding Text from Brain Recordings

abstract: Decoding language representations directly from the brain can enable new Brain-Computer Interfaces (BCI) for high bandwidth human-human and human-machine communication. Clinically, such technologies can restore communication in people with neurological conditions affecting their ability to speak. In this study, we propose a novel deep network architecture Brain2Char, for directly decoding text (specifically character sequences) from direct brain recordings (called Electrocorticography, ECoG). Brain2Char framework combines state-of-the-art deep learning modules --- 3D Inception layers for multiband spatiotemporal feature extraction from neural data and bidirectional recurrent layers, dilated convolution layers followed by language model weighted beam search to decode character sequences, optimizing a connectionist temporal classification (CTC) loss. Additionally, given the highly non-linear transformations that underlie the conversion of cortical function to character sequences, we perform regularizations on the network's latent representations motivated by insights into cortical encoding of speech production and artifactual aspects specific to ECoG data acquisition. To do this, we impose auxiliary losses on latent representations for articulatory movements, speech acoustics and session specific non-linearities. In 3 participants tested here, Brain2Char achieves 10.6\%, 8.5\% and 7.0\% Word Error Rates (WER) respectively on vocabulary sizes ranging from 1200 to 1900 words. Brain2Char also performs well when 2 participants silently mimed sentences. These results set a new state-of-the-art on decoding text from brain and demonstrate the potential of Brain2Char as a high-performance communication BCI.


title: Stability of dynamical quantum phase transitions in quenched topological insulators: From multiband to disordered systems

abstract: Dynamical quantum phase transitions (DQPTs) represent a counterpart in non-equilibrium quantum time evolution of thermal phase transitions at equilibrium, where real time becomes analogous to a control parameter such as temperature. In quenched quantum systems, recently the occurrence of DQPTs has been demonstrated, both with theory and experiment, to be intimately connected to changes of topological properties. Here, we contribute to broadening the systematic understanding of this relation between topology and DQPTs to multi-orbital and disordered systems. Specifically, we provide a detailed ergodicity analysis to derive criteria for DQPTs in all spatial dimensions, and construct basic counter-examples to the occurrence of DQPTs in multi-band topological insulator models. As a numerical case study illustrating our results, we report on microscopic simulations of the quench dynamics in the Harper-Hofstadter model. Furthermore, going gradually from multi-band to disordered systems, we approach random disorder by increasing the (super) unit cell within which random perturbations are switched on adiabatically. This leads to an intriguing order of limits problem which we address by extensive numerical calculations on quenched one-dimensional topological insulators and superconductors with disorder.

categories: cond-mat.quant-gas cond-mat.mes-hall cond-mat.stat-mech quant-ph


title: On ordered $k$-paths and rims for certain families of Kazhdan-Lusztig cells of $S_n$

abstract: For a composition $\lambda$ of $n$ we consider the Kazhdan-Lusztig cell in the symmetric group $S_n$ containing the longest element of the standard parabolic subgroup of $S_n$ associated to $\lambda$. In this paper we extend some of the ideas and results in [{Beitr{\"a}ge} zur Algebra und Geometrie, \textbf{59} (2018), no.~3, 523--547]. In particular, by introducing the notion of an ordered $k$-path, we are able to obtain alternative explicit descriptions for some additional families of cells associated to compositions. This is achieved by first determining the rim of the cell, from which reduced forms for all the elements of the cell are easily obtained.


title: Multiple Populations in Integrated Light Spectroscopy of Intermediate Age Clusters

abstract: The presence of star-to-star light-element abundance variations (a.k.a. multiple populations, MPs) appears to be ubiquitous within old and massive clusters in the Milky Way and all studied nearby galaxies. Most previous studies have focussed on resolved images or spectroscopy of individual stars, although there has been significant effort in the past few years to look for multiple population signatures in integrated light spectroscopy. If proven feasible, integrated light studies offer a potential way to vastly open parameter space, as clusters out to tens of Mpc can be studied. We use the NaD lines in the integrated spectra of two clusters with similar ages ($2-3$ Gyr) but very different masses, NGC 1978 ($\sim3\times10^5$ Msun) in the LMC and G114 ($1.7\times10^7$ Msun) in NGC 1316. For NGC 1978, our findings agree with resolved studies of individual stars which did not find evidence for Na spreads. However, for G114, we find clear evidence for the presence of multiple populations. The fact that the same anomalous abundance patterns are found in both the intermediate age and ancient GCs lends further support to the notion that young massive clusters are effectively the same as the ancient globular clusters, only separated in age.


title: Studying black holes on horizon scales with space-VLBI

abstract: The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.


title: Seasonal variation of atmospheric muons in IceCube

abstract: After more than seven years of data taking with the full IceCube detector triggering at an average rate of 2.15 kHz, a sample of half a trillion muon events is available for analysis. The extreme temperature variations in the stratosphere together with the high data rate reveal features on both long and short time scales with unprecedented precision. In this paper we report an analysis in terms of the atmospheric profile for production of muons from decay of charged pions and kaons. We comment on the implications for seasonal variations of neutrinos, which are presented in a separate paper at this conference.


title: Polarization resolved radiation angular patterns of orientationally ordered nanorods

abstract: We employ the transfer matrix approach combined with the Green's function method to theoretically study polarization resolved far-field angular distributions of photoluminescence from quantum nanorods (NRs) embedded in an anisotropic polymer film. The emission and excitation properties of NRs are described by the emission and excitation anisotropy tensors. These tensors and the solution of the emission problem expressed in terms of the evolution operators are used to derive the orientationally averaged coherency matrix of the emitted wavefield. For the case of in-plane alignment and unpolarized excitation, we estimate the emission anisotropy parameter and compute the angular profiles for the photoluminescence polarization parameter such as the degree of linear polarization, the Stokes parameter $s_1$, the ellipticity and the polarization azimuth. We show that the alignment order parameter has a profound effect on the angular profiles.


title: Extremely long baseline interferometry with Origins Space Telescope

abstract: Operating 1.5 million km from Earth at the Sun-Earth L2 Lagrange point, the Origins Space Telescope equipped with a slightly modified version of its HERO heterodyne instrument could function as a uniquely valuable node in a VLBI network. The unprecedented angular resolution resulting from the combination of Origins with existing ground-based millimeter/submillimeter telescope arrays would increase the number of spatially resolvable black holes by a factor of a million, permit the study of these black holes across all of cosmic history, and enable new tests of general relativity by unveiling the photon ring substructure in the nearest black holes.


title: Use of a controlled experiment and computational models to measure the impact of sequential peer exposures on decision making

abstract: It is widely believed that one's peers influence product adoption behaviors. This relationship has been linked to the number of signals a decision-maker receives in a social network. But it is unclear if these same principles hold when the pattern by which it receives these signals vary and when peer influence is directed towards choices which are not optimal. To investigate that, we manipulate social signal exposure in an online controlled experiment using a game with human participants. Each participant in the game makes a decision among choices with differing utilities. We observe the following: (1) even in the presence of monetary risks and previously acquired knowledge of the choices, decision-makers tend to deviate from the obvious optimal decision when their peers make similar decision which we call the influence decision, (2) when the quantity of social signals vary over time, the forwarding probability of the influence decision and therefore being responsive to social influence does not necessarily correlate proportionally to the absolute quantity of signals. To better understand how these rules of peer influence could be used in modeling applications of real world diffusion and in networked environments, we use our behavioral findings to simulate spreading dynamics in real world case studies. We specifically try to see how cumulative influence plays out in the presence of user uncertainty and measure its outcome on rumor diffusion, which we model as an example of sub-optimal choice diffusion. Together, our simulation results indicate that sequential peer effects from the influence decision overcomes individual uncertainty to guide faster rumor diffusion over time. However, when the rate of diffusion is slow in the beginning, user uncertainty can have a substantial role compared to peer influence in deciding the adoption trajectory of a piece of questionable information.


title: Oblivious Sketching of High-Degree Polynomial Kernels

abstract: Kernel methods are fundamental tools in machine learning that allow detection of non-linear dependencies between data without explicitly constructing feature vectors in high dimensional spaces. A major disadvantage of kernel methods is their poor scalability: primitives such as kernel PCA or kernel ridge regression generally take prohibitively large quadratic space and (at least) quadratic time, as kernel matrices are usually dense. Some methods for speeding up kernel linear algebra are known, but they all invariably take time exponential in either the dimension of the input point set (e.g., fast multipole methods suffer from the curse of dimensionality) or in the degree of the kernel function. Oblivious sketching has emerged as a powerful approach to speeding up numerical linear algebra over the past decade, but our understanding of oblivious sketching solutions for kernel matrices has remained quite limited, suffering from the aforementioned exponential dependence on input parameters. Our main contribution is a general method for applying sketching solutions developed in numerical linear algebra over the past decade to a tensoring of data points without forming the tensoring explicitly. This leads to the first oblivious sketch for the polynomial kernel with a target dimension that is only polynomially dependent on the degree of the kernel function, as well as the first oblivious sketch for the Gaussian kernel on bounded datasets that does not suffer from an exponential dependence on the dimensionality of input data points.


title: Studying Black Holes on Horizon Scales with VLBI Ground Arrays

abstract: High-resolution imaging of supermassive black holes is now possible, with new applications to testing general relativity and horizon-scale accretion and relativistic jet formation processes. Over the coming decade, the EHT will propose to add new strategically placed VLBI elements operating at 1.3mm and 0.87mm wavelength. In parallel, development of next-generation backend instrumentation, coupled with high throughput correlation architectures, will boost sensitivity, allowing the new stations to be of modest collecting area while still improving imaging fidelity and angular resolution. The goal of these efforts is to move from imaging static horizon scale structure to dynamic reconstructions that capture the processes of accretion and jet launching in near real time.


title: Mixture Probabilistic Principal Geodesic Analysis

abstract: Dimensionality reduction on Riemannian manifolds is challenging due to the complex nonlinear data structures. While probabilistic principal geodesic analysis~(PPGA) has been proposed to generalize conventional principal component analysis (PCA) onto manifolds, its effectiveness is limited to data with a single modality. In this paper, we present a novel Gaussian latent variable model that provides a unique way to integrate multiple PGA models into a maximum-likelihood framework. This leads to a well-defined mixture model of probabilistic principal geodesic analysis (MPPGA) on sub-populations, where parameters of the principal subspaces are automatically estimated by employing an Expectation Maximization algorithm. We further develop a mixture Bayesian PGA (MBPGA) model that automatically reduces data dimensionality by suppressing irrelevant principal geodesics. We demonstrate the advantages of our model in the contexts of clustering and statistical shape analysis, using synthetic sphere data, real corpus callosum, and mandible data from human brain magnetic resonance~(MR) and CT images.


title: An analytical perturbative solution to the Merton Garman model using symmetries

abstract: In this paper, we introduce an analytical perturbative solution to the Merton Garman model. It is obtained by doing perturbation theory around the exact analytical solution of a model which possesses a two-dimensional Galilean symmetry. We compare our perturbative solution of the Merton Garman model to Monte Carlo simulations and find that our solutions performs surprisingly well for a wide range of parameters. We also show how to use symmetries to build option pricing models. Our results demonstrate that the concept of symmetry is important in mathematical finance.


title: From type theory to setoids and back

abstract: A model of Martin-L\"of extensional type theory with universes is formalized in Agda, an interactive proof system based on Martin-L\"of intensional type theory. This may be understood, we claim, as a solution to the old problem of modelling the full extensional theory in the intensional theory. Types are interpreted as setoids, and the model is therefore a setoid model. We solve the problem of intepreting type universes by utilizing Aczel's type of iterative sets, and show how it can be made into a setoid of small setoids containing the necessary setoid constructions. In addition we interpret the bracket types of Awodey and Bauer. Further quotient types should be interpretable.


title: Reliability Bounds for Dependent Fading Wireless Channels

abstract: Unreliable fading wireless channels are the main challenge for strict performance guarantees in mobile communications. Diversity schemes including massive number of antennas, huge spectrum bands and multi-connectivity links are applied to improve the outage performance. The success of these approaches relies heavily on the joint distribution of the underlying fading channels. In this work, we consider the $\varepsilon$-outage capacity of slowly fading wireless diversity channels and provide lower and upper bounds for fixed marginal distributions of the individual channels. This answers the question about the best and worst case outage probability achievable over $n$ fading channels with a given distribution, e.g., Rayleigh fading, but not necessarily statistically independent. Interestingly, the best-case joint distribution enables achieving a zero-outage capacity greater than zero without channel state information at the transmitter for $n \geq 2$. Furthermore, the results are applied to characterize the worst- and best-case joint distribution for zero-outage capacity with perfect channel state information everywhere. All results are specialized to Rayleigh fading and compared to the standard assumption of independent and identically distributed fading component channels. The results show a significant impact of the joint distribution and the gap between worst- and best-case can be arbitrarily large.


title: Thing/Machine-s (Thimacs) Applied to Structural Description in Software Engineering

abstract: We are pursuing a modeling methodology that views the world as a realm of things. A thing is defined as something that can be created, processed, released, transferred, and received. Additionally, in this modeling approach, a thing is a five-dimensional structure referred to as a thinging (abstract) machine. On the other hand, machines are things that are operated on; that is, they are created, processed, released, transferred, and received. The intertwining with the world is accomplished by integrating these two modes of an entity s being: being a thing that flows through machines and being a machine that processes things. This paper further enriches these notions of things and machines. We present further exploration of the thinging machine model through introducing a new notion called the thing/machine (thimac) as a label of the unity of things/machines. Thimacs replace traditional categorization, properties, and behavior with creating, processing, releasing, transferring, and receiving, as well as the two linking notions of flow and triggering. The paper discusses the concept of thimacs with examples and focuses on the notion of structure as it applies to various diagrammatic modeling methodologies.


title: Multi-level Attention network using text, audio and video for Depression Prediction

abstract: Depression has been the leading cause of mental-health illness worldwide. Major depressive disorder (MDD), is a common mental health disorder that affects both psychologically as well as physically which could lead to loss of lives. Due to the lack of diagnostic tests and subjectivity involved in detecting depression, there is a growing interest in using behavioural cues to automate depression diagnosis and stage prediction. The absence of labelled behavioural datasets for such problems and the huge amount of variations possible in behaviour makes the problem more challenging. This paper presents a novel multi-level attention based network for multi-modal depression prediction that fuses features from audio, video and text modalities while learning the intra and inter modality relevance. The multi-level attention reinforces overall learning by selecting the most influential features within each modality for the decision making. We perform exhaustive experimentation to create different regression models for audio, video and text modalities. Several fusions models with different configurations are constructed to understand the impact of each feature and modality. We outperform the current baseline by 17.52% in terms of root mean squared error.


title: A categorification of the Malvenuto--Reutenauer algebra via a tower of groups

abstract: There is a long tradition of categorifying combinatorial Hopf algebras by the modules of a tower of algebras (or even better via the representation theory of a tower of groups). From the point of view of combinatorics, such a categorification supplies canonical bases, inner products, and a natural avenue to prove positivity results. Recent ideas in supercharacter theory have made fashioning the representation theory of a tower of groups into a Hopf structure more tractable. This paper applies such a program to the Malvenuto--Reutenauer Hopf algebra. In particular, we design functors on the representation theory of a tower of $p$-groups that realize the Hopf structure of the Malvenuto--Reutenauer algebra in such a way that its well-known fundamental basis corresponds to a supercharacter basis.


title: Learning Koopman Eigenfunctions and Invariant Subspaces from Data: Symmetric Subspace Decomposition

abstract: This paper develops data-driven methods to identify eigenfunctions of the Koopman operator associated to a dynamical system and subspaces that are invariant under the operator. We build on Extended Dynamic Mode Decomposition (EDMD), a data-driven method that finds a finite-dimensional approximation of the Koopman operator on the span of a predefined dictionary of functions. We propose a necessary and sufficient condition to identify Koopman eigenfunctions based on the application of EDMD forward and backward in time. Moreover, we propose the Symmetric Subspace Decomposition (SSD) algorithm, an iterative method which provably identifies the maximal Koopman-invariant subspace and the Koopman eigenfunctions in the span of the dictionary. We also introduce the Streaming Symmetric Subspace Decomposition (SSSD) algorithm, an online extension of SSD that only requires a small, fixed memory and incorporates new data as is received. Finally, we propose an extension of SSD that approximates Koopman eigenfunctions and invariant subspaces when the dictionary does not contain sufficient informative eigenfunctions.


title: The interplay of magnetically-dominated turbulence and magnetic reconnection in producing nonthermal particles

abstract: Magnetized turbulence and magnetic reconnection are often invoked to explain the nonthermal emission observed from a wide variety of astrophysical sources. By means of fully-kinetic 2D and 3D PIC simulations, we investigate the interplay between turbulence and reconnection in generating nonthermal particles in magnetically-dominated pair plasmas. A generic by-product of the turbulence evolution is the generation of a nonthermal particle spectrum with a power-law energy range. The power-law slope $p$ is harder for larger magnetizations and stronger turbulence fluctuations, and it can be as hard as $p < 2$. The Larmor radius of particles at the high-energy cutoff is comparable to the size $l$ of the largest eddies. Plasmoid-mediated reconnection, which self-consistently occurs in the turbulent plasma, controls the physics of particle injection. Then, particles are further accelerated by stochastic scattering off turbulent fluctuations. The work done by parallel electric fields - naturally expected in reconnection layers - is responsible for most of the initial energy increase, and is proportional to the magnetization $\sigma$ of the system, while the subsequent energy gain, which dominates the overall energization of high-energy particles, is powered by the perpendicular electric fields of turbulent fluctuations. The two-stage acceleration process leaves an imprint in the particle pitch-angle distribution: low-energy particles are aligned with the field, while the highest energy particles move preferentially orthogonal to it. The energy diffusion coefficient of stochastic acceleration scales as $D_\gamma\sim 0.1\sigma(c/l)\gamma^2$, where $\gamma$ is the particle Lorentz factor. This results in fast acceleration timescales $t_{acc}\sim (3/\sigma)\,l/c$. Our findings have important implications for understanding the generation of nonthermal particles in high-energy astrophysical sources.


title: Mining Insights from Weakly-Structured Event Data

abstract: This thesis focuses on process mining on event data where such a normative specification is absent and, as a result, the event data is less structured. The thesis puts special emphasis on one application domain that fits this description: the analysis of smart home data where sequences of daily activities are recorded. In this thesis we propose a set of techniques to analyze such data, which can be grouped into two categories of techniques. The first category of methods focuses on preprocessing event logs in order to enable process discovery techniques to extract insights into unstructured event data. In this category we have developed the following techniques: - An unsupervised approach to refine event labels based on the time at which the event took place, allowing for example to distinguish recorded eating events into breakfast, lunch, and dinner. - An approach to detect and filter from event logs so-called chaotic activities, which are activities that cause process discovery methods to overgeneralize. - A supervised approach to abstract low-level events into more high-level events, where we show that there exist situations where process discovery approaches overgeneralize on the low-level event data but are able to find precise models on the high-level event data. The second category focuses on mining local process models, i.e., collections of process model patterns that each describe some frequent pattern, in contrast to the single global process model that is obtained with existing process discovery techniques. Several techniques are introduced in the area of local process model mining, including a basic method, fast but approximate heuristic methods, and constraint-based techniques.


title: Optimization with Equality and Inequality Constraints Using Parameter Continuation

abstract: We generalize the successive continuation paradigm introduced by Kern\'evez and Doedel [16] for locating locally optimal solutions of constrained optimization problems to the case of simultaneous equality and inequality constraints. The analysis shows that potential optima may be found at the end of a sequence of easily-initialized separate stages of continuation, without the need to seed the first stage of continuation with nonzero values for the corresponding Lagrange multipliers. A key enabler of the proposed generalization is the use of complementarity functions to define relaxed complementary conditions, followed by the use of continuation to arrive at the limit required by the Karush-Kuhn-Tucker theory. As a result, a successful search for optima is found to be possible also from an infeasible initial solution guess. The discussion shows that the proposed paradigm is compatible with the staged construction approach of the COCO software package. This is evidenced by a modified form of the COCO core used to produce the numerical results reported here. These illustrate the efficacy of the continuation approach in locating stationary solutions of an objective function along families of two-point boundary value problems and in optimal control problems.


title: Exploration Without Global Consistency Using Local Volume Consolidation

abstract: In exploration, the goal is to build a map of an unknown environment. Most state-of-the-art approaches use map representations that require drift-free state estimates to function properly. Real-world state estimators, however, exhibit drift. In this paper, we present a 2D map representation for exploration that is robust to drift. Rather than a global map, it uses local metric volumes connected by relative pose estimates. This pose-graph does not need to be globally consistent. Overlaps between the volumes are resolved locally, rather than on the faulty estimate of space. We demonstrate our representation with a frontier-based exploration approach, evaluate it under different conditions and compare it with a commonly-used grid-based representation. We show that, at the cost of longer exploration time, using the proposed representation allows full coverage of space even for very large drift in the state estimate, contrary to the grid-based representation. The system is validated in a real world experiment and we discuss its extension to 3D.


title: Runaway O-star Bow Shocks as Particle Accelerators? The Case of AE Aur revisited

abstract: We present results of our Chandra/ACIS observations of the field centered on the fast, runaway O star AE Aur and its bow shock. Previous XMM-Newton observations revealed an X-ray "blob" near the IR arc tracing the bow shock, possibly a nonthermal source consistent with models of Inverse Compton scattering of dust IR photons by electrons accelerated at the shock. The new, subarcsecond resolution Chandra data, while confirming the presence of the XMM-Newton source, clearly indicate that the latter is neither extended nor coincident with the IR arc and strongly suggest it is a background AGN. Motivated by results published for the bow shock of BD+43 3654, we extended our study to the radio domain, by analyzing archival EVLA data. We find no radio emission from the AE Aur bow shock either. The corresponding upper limits for the absorbed (unabsorbed) X-ray flux of 5.9(7.8)x10^-15 erg/cm^2/s (3 sigma) and, in the radio range, of 2 mJy (1.4 GHz), and 0.4 mJy (5.0 GHz), are used to put constraints on model predictions for particle acceleration within the bow shock. In the "classical" framework of Diffusive Shock Acceleration, we find that the predicted X-ray and radio emission by the bow shock is at least two orders of magnitude below the current upper limits, consistent with the systematic non-detections of up to 60 stellar bow shocks. The only exception so far remains that of BD+43 3654, probably the result of its very large mass-loss rate among runaway O stars.


title: Topological Transition on the Conformal Manifold

abstract: Despite great successes in the study of gapped phases, a comprehensive understanding of the gapless phases and their transitions is still under developments. In this paper, we study a general phenomenon in the space of (1+1)$d$ critical phases with fermionic degrees of freedom described by a continuous family of conformal field theories (CFT), a.k.a. the conformal manifold. Along a one-dimensional locus on the conformal manifold, there can be a transition point, across which the fermionic CFTs on the two sides differ by stacking an invertible fermionic topological order (IFTO), point-by-point along the locus. At every point on the conformal manifold, the order and disorder operators have power-law two-point functions, but their critical exponents cross over with each other at the transition point, where stacking the IFTO leaves the fermionic CFT unchanged. We call this continuous transition on the fermionic conformal manifold a topological transition. By gauging the fermion parity, the IFTO stacking becomes a Kramers-Wannier duality between the corresponding bosonic CFTs. Both the IFTO stacking and the Kramers-Wannier duality are induced by the electromagnetic duality of the (2+1)$d$ $\mathbb{Z}_2$ topological order. We provide several examples of topological transitions, including the familiar Luttinger model of spinless fermions (i.e. the $c=1$ massless Dirac fermion with the Thirring interaction), and a new class of $c=2$ examples describing $U(1)\times SU(2)$-protected gapless phases.


title: Clean, Robust Alkali Sources by Intercalation within Highly-Oriented Pyrolytic Graphite

abstract: We report the fabrication, characterization, and use of rubidium vapor dispensers based on highly-oriented pyrolytic graphite (HOPG) intercalated with metallic rubidium. Compared to commercial chromate salt dispensers, these intercalated HOPG (IHOPG) dispensers hold an order of magnitude more rubidium in a similar volume, require less than one-fourth the heating power, and emit less than one-half as many impurities. Appropriate processing permits exposure of the IHOPG to atmosphere for over ninety minutes without any adverse effects. Intercalation of cesium and potassium into HOPG have also been demonstrated in the literature, which suggests that IHOPG dispensers may also be be made for those metals.


title: Separating subgroups of mapping class groups in homological representations

abstract: Let $\Gamma$ be either the mapping class group of a closed surface of genus $\geq 2$, or the automorphism group of a free group of rank $\geq 3$. Given any homological representation $\rho$ of $\Gamma$ corresponding to a finite cover, and any term $\mathcal{I}_k$ of the Johnson filtration, we show that $\rho(\mathcal{I}_k)$ has finite index in $\rho(\mathcal{I})$, the Torelli subgroup of $\Gamma$. Since $[\mathcal{I}: \mathcal{I}_k] = \infty$ for $k > 1$, this implies for instance that no such representation is faithful.


title: Map plasticity

abstract: With the arrival of digital maps, the ubiquity of maps has increased sharply and new map functionalities have become available such as changing the scale on the fly or displaying/hiding layers. Users can now interact with maps on multiple devices (e.g. smartphones, desktop computers, large-scale displays, head-mounted displays) using different means of interaction such as touch, voice or gestures. However, ensuring map functionalities and good user experience across these devices and modalities frequently entails dedicated development efforts for each combination. In this paper, we argue that introducing an abstract representation of what a map contains and affords can unlock new opportunities. For this purpose, we propose the concept of map plasticity, the capability of a map-based system to support different contexts of use while preserving usability and functionality. Based on this definition, we discuss core components and an example. We also propose a research agenda for realising map plasticity and its benefits.


title: Statistical physics and mesoscopic modeling to interpret tethered particle motion experiments

abstract: Tethered particle motion experiments are versatile single-molecule techniques enabling one to address in vitro the molecular properties of DNA and its interactions with various partners involved in genetic regulations. These techniques provide raw data such as the tracked particle amplitude of movement, from which relevant information about DNA conformations or states must be recovered. Solving this inverse problem appeals to specific theoretical tools that have been designed in the two last decades, together with the data pre-processing procedures that ought to be implemented to avoid biases inherent to these experimental techniques. These statistical tools and models are reviewed in this paper.


title: Detecting Depinning and Nonequilibrium Transitions with Unsupervised Machine Learning

abstract: Using numerical simulations of a model disk system, we demonstrate that a machine learning generated order parameter can detect depinning transitions and different dynamic flow phases in systems driven far from equilibrium. We specifically consider monodisperse passive disks with short range interactions undergoing a depinning phase transition when driven over quenched disorder. The machine learning derived order parameter identifies the depinning transition as well as different dynamical regimes, such as the transition from a flowing liquid to a phase separated liquid-solid state that is not readily distinguished with traditional measures such as velocity-force curves or Voronoi tessellation. The order parameter also shows markedly distinct behavior in the limit of high density where jamming effects occur. Our results should be general to the broad class of particle-based systems that exhibit depinning transitions and nonequilibrium phase transitions.


title: Trapped state at a dislocation in a weak magnetomechanical topological insulator

abstract: Topological insulators (TIs) are characterized by an insulating bulk and symmetry protected bound state on their boundaries. A "strong" topological insulator is characterized by robust conducting states on \emph{all} boundaries protected by certain internal symmetries. A "weak" topological insulator (WTI) however, requires lattice translation symmetry, making it more sensitive to disorder. However, this sensitivity gives rise to interesting characteristics such as anisotropic edge modes, quantized charge polarization, and bound states appearing at dislocation defects. Despite hosting interesting features, the sensitivity of WTIs to disorder poses an experimental confirmation challenge. Here we realize a 2D magneto-mechanical metamaterial and demonstrate experimentally the unique features of a WTI. Specifically, we show that the 2D WTI is anisotropic and hosts edge modes only on certain edges, as well as hosting a bound state at a dislocation defect. We construct the 2D WTI from stacked 1D SSH chains for which we show experimentally the different gapped phases of the 1D model.


title: Adversarial Robustness of Similarity-Based Link Prediction

abstract: Link prediction is one of the fundamental problems in social network analysis. A common set of techniques for link prediction rely on similarity metrics which use the topology of the observed subnetwork to quantify the likelihood of unobserved links. Recently, similarity metrics for link prediction have been shown to be vulnerable to attacks whereby observations about the network are adversarially modified to hide target links. We propose a novel approach for increasing robustness of similarity-based link prediction by endowing the analyst with a restricted set of reliable queries which accurately measure the existence of queried links. The analyst aims to robustly predict a collection of possible links by optimally allocating the reliable queries. We formalize the analyst problem as a Bayesian Stackelberg game in which they first choose the reliable queries, followed by an adversary who deletes a subset of links among the remaining (unreliable) queries by the analyst. The analyst in our model is uncertain about the particular target link the adversary attempts to hide, whereas the adversary has full information about the analyst and the network. Focusing on similarity metrics using only local information, we show that the problem is NP-Hard for both players, and devise two principled and efficient approaches for solving it approximately. Extensive experiments with real and synthetic networks demonstrate the effectiveness of our approach.


title: Simultaneous Cubic and Quadratic Diagonal Equations In 12 Prime Variables

abstract: The system of equations \[ u_1p_1^2 + \ldots + u_sp_s^2 = 0 \] \[ v_1p_1^3 + \ldots + v_sp_s^3 = 0 \] has prime solutions $(p_1, \ldots, p_s)$ for $s \geq 12$, assuming that the system has solutions modulo each prime $p$. This is proved via the Hardy-Littlewood circle method, building on Wooley's work on the corresponding system over the integers and recent results on Vinogradov's mean value theorem. Additionally, a set of sufficient conditions for local solvability is given: If both equations are solvable modulo 2, the quadratic equation is solvable modulo 3, and for each prime $p$ at least 7 of each of the $u_i$, $v_i$ are not zero modulo $p$, then the system has solutions modulo each prime $p$.


title: Towards Hilbert's Tenth Problem for rings of integers through Iwasawa theory and Heegner points

abstract: For a positive proportion of primes $p$ and $q$, we prove that $\mathbb{Z}$ is Diophantine in the ring of integers of $\mathbb{Q}(\sqrt[3]{p},\sqrt{-q})$. This provides a new and explicit infinite family of number fields $K$ such that Hilbert's tenth problem for $O_K$ is unsolvable. Our methods use Iwasawa theory and congruences of Heegner points in order to obtain suitable rank stability properties for elliptic curves.


title: Submillimeter radiation as the thermal component of the Neupert Effect

abstract: The Neupert effect is the empirical observation that the time evolution of non-thermal emission (e.g. hard X-rays) is frequently proportional to the time derivative of the thermal emission flux (soft X-rays), or, vice versa, that time integrated non-thermal flux is proportional to thermal flux. We analyzed the GOES M2.2 event SOL2011-02-14T17:25, and found that the 212 GHz emission plays quite well the role of the thermal component of the Neupert effect. We show that the maximum of the hard X-ray flux for energies above 50 keV is coincident in time with the time-derivative of the 212 GHz flux, within the uncertainties. The microwave flux density at 15.4 GHz, produced by optically thin gyrosynchrotron mechanism, and hard-X rays above 25 keV mark the typical impulsive phase, and have similar time evolution. On the other hand, the 12 GHz emission is delayed by about 25 seconds with respect of the microwave and hard X-ray peak. We argue that this delay cannot be explained by magnetic trapping of non-thermal electrons. With all the observational evidence, we suggest that the 212 GHz emission is produced by thermal bremsstrahlung, initially in the chromosphere, and shifting to optically thin emission from thehot coronal loops at the end of the gradual phase.


title: Discriminative Topic Modeling with Logistic LDA

abstract: Despite many years of research into latent Dirichlet allocation (LDA), applying LDA to collections of non-categorical items is still challenging. Yet many problems with much richer data share a similar structure and could benefit from the vast literature on LDA. We propose logistic LDA, a novel discriminative variant of latent Dirichlet allocation which is easy to apply to arbitrary inputs. In particular, our model can easily be applied to groups of images, arbitrary text embeddings, and integrates well with deep neural networks. Although it is a discriminative model, we show that logistic LDA can learn from unlabeled data in an unsupervised manner by exploiting the group structure present in the data. In contrast to other recent topic models designed to handle arbitrary inputs, our model does not sacrifice the interpretability and principled motivation of LDA.


title: Small-worldness favours network inference

abstract: A main goal in the analysis of a complex system is to infer its underlying network structure from time-series observations of its behaviour. The inference process is often done by using bi-variate similarity measures, such as the cross-correlation (CC), however, the main factors favouring or hindering its success are still puzzling. Here, we use synthetic neuron models in order to reveal the main topological properties that frustrate or facilitate inferring the underlying network from CC measurements. Specifically, we use pulse-coupled Izhikevich neurons connected as in the Caenorhabditis elegans neural networks as well as in networks with similar randomness and small-worldness. We analyse the effectiveness and robustness of the inference process under different observations and collective dynamics, contrasting the results obtained from using membrane potentials and inter-spike interval time-series. We find that overall, small-worldness favours network inference and degree heterogeneity hinders it. In particular, success rates in C. elegans networks -- that combine small-world properties with degree heterogeneity -- are closer to success rates in Erd\"os-R\'enyi network models rather than those in Watts-Strogatz network models. These results are relevant to understand better the relationship between topological properties and function in different neural networks.


title: Opaque Lowermost Mantle

abstract: Earth/s lowermost mantle displays complex geological structures that likely result from heterogeneous thermal and electromagnetic interaction with the core. Geophysical models of the core-mantle boundary (CMB) region rely on the thermal and electrical conductivities of appropriate geomaterials which, however, have never been probed at representative pressure and temperature (P-T) conditions. Here we report on the opacity of single crystalline bridgmanite and ferropericlase, which is linked to both their radiative and electrical conductivity, measured in dynamically- and statically-heated diamond anvil cells as well as computed from first-principles at CMB conditions. Our results show that light absorption in the visible spectral range is enhanced upon heating in both minerals but the rate of change in opacity with temperature is a factor of six higher in ferropericlase. As a result, bridgmanite in the lowermost mantle is moderately transparent while ferropericlase is highly opaque. Our measurements suggest a very low (< 1 W/m/K) and largely temperature-independent radiative conductivity in the lowermost mantle, at odds with previous studies. This implies that the radiative mechanism has not contributed significantly to cooling the Earth/s core throughout the geologic time and points to a present-day CMB heat flow of 9-11 TW. Opaque ferropericlase is electrically conducting and mediates strong core-mantle electromagnetic coupling, explaining the intradecadal oscillations in the length of day, low secular geomagnetic variations in Central Pacific, and the preferred paths of geomagnetic pole reversals.


title: Universality of clone dynamics during tissue development

abstract: The emergence of complex organs is driven by the coordinated proliferation, migration and differentiation of precursor cells. The fate behaviour of these cells is reflected in the time evolution their progeny, termed clones, which serve as a key experimental observable. In adult tissues, where cell dynamics is constrained by the condition of homeostasis, clonal tracing studies based on transgenic animal models have advanced our understanding of cell fate behaviour and its dysregulation in disease. But what can be learned from clonal dynamics in development, where the spatial cohesiveness of clones is impaired by tissue deformations during tissue growth? Drawing on the results of clonal tracing studies, we show that, despite the complexity of organ development, clonal dynamics may converge to a critical state characterized by universal scaling behaviour of clone sizes. By mapping clonal dynamics onto a generalization of the classical theory of aerosols, we elucidate the origin and range of scaling behaviours and show how the identification of universal scaling dependences may allow lineage-specific information to be distilled from experiments. Our study shows the emergence of core concepts of statistical physics in an unexpected context, identifying cellular systems as a laboratory to study non-equilibrium statistical physics.


title: LCA: Loss Change Allocation for Neural Network Training

abstract: Neural networks enjoy widespread use, but many aspects of their training, representation, and operation are poorly understood. In particular, our view into the training process is limited, with a single scalar loss being the most common viewport into this high-dimensional, dynamic process. We propose a new window into training called Loss Change Allocation (LCA), in which credit for changes to the network loss is conservatively partitioned to the parameters. This measurement is accomplished by decomposing the components of an approximate path integral along the training trajectory using a Runge-Kutta integrator. This rich view shows which parameters are responsible for decreasing or increasing the loss during training, or which parameters "help" or "hurt" the network's learning, respectively. LCA may be summed over training iterations and/or over neurons, channels, or layers for increasingly coarse views. This new measurement device produces several insights into training. (1) We find that barely over 50% of parameters help during any given iteration. (2) Some entire layers hurt overall, moving on average against the training gradient, a phenomenon we hypothesize may be due to phase lag in an oscillatory training process. (3) Finally, increments in learning proceed in a synchronized manner across layers, often peaking on identical iterations.


title: CrossWeigh: Training Named Entity Tagger from Imperfect Annotations

abstract: Everyone makes mistakes. So do human annotators when curating labels for named entity recognition (NER). Such label mistakes might hurt model training and interfere model comparison. In this study, we dive deep into one of the widely-adopted NER benchmark datasets, CoNLL03 NER. We are able to identify label mistakes in about 5.38% test sentences, which is a significant ratio considering that the state-of-the-art test F1 score is already around 93%. Therefore, we manually correct these label mistakes and form a cleaner test set. Our re-evaluation of popular models on this corrected test set leads to more accurate assessments, compared to those on the original test set. More importantly, we propose a simple yet effective framework, CrossWeigh, to handle label mistakes during NER model training. Specifically, it partitions the training data into several folds and train independent NER models to identify potential mistakes in each fold. Then it adjusts the weights of training data accordingly to train the final NER model. Extensive experiments demonstrate significant improvements of plugging various NER models into our proposed framework on three datasets. All implementations and corrected test set are available at our Github repo: https://github.com/ZihanWangKi/CrossWeigh.


title: Dispersive Landau levels and valley currents in strained graphene nanoribbons

abstract: We describe a simple setup generating pure valley currents -- valley transport without charge transport -- in strained graphene nanoribbons with zigzag edges. The crucial ingredient is a uniaxial strain pattern which couples to the low-energy Dirac electrons as a uniform pseudomagnetic field. Remarkably, the resulting pseudo-Landau levels are not flat but disperse linearly from the Dirac points, with an opposite slope in the two valleys. We show how this is a natural consequence of an inhomogeneous Fermi velocity arising in the low-energy theory describing the system, which maps to an exactly-solvable singular Sturm-Liouville problem. The velocity of the valley currents can be controlled by tuning the magnitude of strain and by applying bias voltages across the ribbon. Furthermore, applying an electric field along the ribbon leads to pumping of charge carriers between the two valleys, realizing a valley analog of the chiral anomaly in one spatial dimension. These effects produce unique signatures that can be observed experimentally by performing ordinary charge transport measurements and spectroscopy.


title: Self-consistent approach to the description of relaxation processes in classical multiparticle systems

abstract: The concept of time correlation functions is a very convenient theoretical tool in describing relaxation processes in multiparticle systems because, on one hand, correlation functions are directly related to experimentally measured quantities (for example, intensities in spectroscopic studies and kinetic coefficients via the Kubo-Green relation) and, on the other hand, the concept is also applicable beyond the equilibrium case. We show that the formalism of memory functions and the method of recurrence relations allow formulating a self-consistent approach for describing relaxation processes in classical multiparticle systems without needing a priori approximations of time correlation functions by model dependencies and with the satisfaction of sum rules and other physical conditions guaranteed. We also demonstrate that the approach can be used to treat the simplest relaxation scenarios and to develop microscopic theories of transport phenomena in liquids, the propagation of density fluctuations in equilibrium simple liquids, and structure relaxation in supercooled liquids. This approach generalizes the mode-coupling approximation in the G\"{o}tze-Leutheusser realization and the Yulmetyev-Shurygin correlation approximations.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.mtrl-sci physics.chem-ph physics.class-ph


title: Gamma-ray Emitting Narrow-Line Seyfert 1 Galaxies: Past, Present, and Future

abstract: This article reviews our current understanding about $\gamma$-ray detected narrow-line Seyfert 1 ($\gamma$-NLSy1) galaxies. The detection with the Large Area Telescope onboard {\it Fermi}~Gamma-ray Space Telescope has provided the strongest evidence for the presence of closely aligned relativistic jet in these intriguing active galactic nuclei (AGN) and opened up a realm to explore the physical conditions needed to launch the jet in a different central engine and host galaxy environment than that is known for blazars. Promising results acquired from various multi-wavelength campaigns are converging to a scenario in which the $\gamma$-NLSy1 galaxies can be considered as `young' blazars. These enigmatic sources hold the key to unravel the jet triggering mechanism and evolution of the AGN phase of a galaxy, in general. As such, $\gamma$-NLSy1s should be considered as one of the top priority targets for next generation observational facilities.


title: Zero-sum Stochastic Games with Asymmetric Information

abstract: A general model for zero-sum stochastic games with asymmetric information is considered. In this model, each player's information at each time can be divided into a common information part and a private information part. Under certain conditions on the evolution of the common and private information, a dynamic programming characterization of the value of the game (if it exists) is presented. If the value of the zero-sum game does not exist, then the dynamic program provides bounds on the upper and lower values of the game. This dynamic program is then used for a class of zero-sum stochastic games with complete information on one side and partial information on the other, that is, games where one player has complete information about state, actions and observation history while the other player may only have partial information about the state and action history. For such games, it is shown that the value exists and can be characterized using the dynamic program. It is further shown that for this class of games, the dynamic program can be used to compute an equilibrium strategy for the more informed player in which the player selects its action using its private information and the common information belief.


title: Eruptions and Flaring Activity in Emerging Quadrupolar Regions

abstract: Context. Some of the most dynamic active regions are associated with complex photospheric magnetic configurations such as quadrupolar regions, and especially ones with a $\delta$-spot configuration and a strong Polarity Inversion Line (PIL). Aims. We study the formation and eruption of magnetic flux ropes in quadrupolar regions. Methods. 3D MHD simulations of the partial emergence of a highly twisted flux tube from the solar interior into a non-magnetized stratified atmosphere. Results. We emerge two $\Omega$-shaped loops forming a quadrupolar region. The emerging flux forms two initially separated bipoles, that later come in contact creating a $\delta$-spot central region. Above the two bipoles, two magnetic lobes expand and interact through a series of current sheets at the interface between them. Two recurrent confined eruptions are produced. In both cases, reconnection between sheared low-lying field lines forms a flux rope. Reconnection between the two lobes high in the atmosphere forms retracting field lines that push against the flux rope, creating a current sheet between them. It also forms a third magnetic lobe between the two emerged ones, that later acts as a strapping field. The flux rope eruptions are triggered when the reconnection between the flux ropes and the field above them becomes efficient enough to remove the tension of the overlying field. These reconnection events occur internally in the system. The first erupting flux rope almost fully reconnects with the overlying field. The second eruption is confined by the overlying field. During the confined eruption, the flux rope is enhanced in size, flux and twist, similar to confined-flare-to-flux-rope observations. Proxies of the emission reveal the two erupting filaments channels. A flare arcade is formed only in the second eruption due to the longer-lasting and more efficient reconnection below the flux rope.


title: Weak Analytic Geometry and a Trace Formula for Families of $p$-adic Representations

abstract: The eigencurve is a powerful tool introduced by Coleman and Mazur to study $p$-adic families of overconvergent modular forms. In this article, we introduce an analogous set of tools for understanding families of "overconvergent" $p$-adic representations of $\pi_1(X)$, where $X$ is a smooth affine variety over a finite field of characteristic $p$. Our main theorem is a trace formula relating the $L$-function of such a family to the geometry of a sequence of associated eigenvarieties. In the case of a single $p$-adic representation, our result reduces to the well known trace formula of Monsky. We apply our theory to the study of $T$-adic exponential sums attached to $\mathbb{Z}_p$-towers over $X$. Special cases of this theory have been applied by Davis, Wan, and Xiao to prove a spectral halo decomposition of the eigencurve attached to $\mathbb{Z}_p$-towers over $X=\mathbb{A}^1$.


title: Reflective prolate-spheroidal operators and the KP/KdV equations

abstract: Commuting integral and differential operators connect the topics of Signal Processing, Random Matrix Theory, and Integrable Systems. Previously, the construction of such pairs was based on direct calculation and concerned concrete special cases, leaving behind important families such as the operators associated to the rational solutions of the KdV equation. We prove a general theorem that the integral operator associated to every wave function in the infinite dimensional Adelic Grassmannian Gr ad of Wilson always reflects a differential operator (in the sense of Definition 1 below). This intrinsic property is shown to follow from the symmetries of Grassmannians of KP wave functions, where the direct commutativity property holds for operators associated to wave functions fixed by Wilson's sign involution but is violated in general. Based on this result, we prove a second main theorem that the integral operators in the computation of the singular values of the truncated generalized Laplace transforms associated to all bispectral wave functions of rank 1 reflect a differential operator. A 90$\deg$ rotation argument is used to prove a third main theorem that the integral operators in the computation of the singular values of the truncated generalized Fourier transforms associated to all such KP wave functions commute with a differential operator. These methods produce vast collections of integral operators with prolate-spheroidal properties, including as special cases the integral operators associated to all rational solutions of the KdV and KP hierarchies considered by Airault-McKean-Moser and Krichever, respectively, in the late 70's. Many novel examples are presented.


title: Boron Arsenide Heterostructures: Lattice-Matched Heterointerfaces, and Strain Effects on Band Alignments and Mobility

abstract: BAs is III-V semiconductor with ultra-high thermal conductivity, but many of its electronic properties are unknown. This work applies predictive atomistic calculations to investigate the properties of BAs heterostructures, such as strain effects on band alignments and carrier mobility, considering BAs as both a thin film and a substrate for lattice-matched materials. The results show that strain decreases the band gap independent of sign or direction. In addition, biaxial tensile strain increases the in-plane electron and hole mobilities by more than 60% compared to the unstrained values due to a reduction of the electron effective mass and of hole interband scattering. Moreover, BAs is shown to be nearly lattice-matched with InGaN and ZnSnN2, two important optoelectronic semiconductors with tunable band gaps by alloying and cation disorder, respectively. The results predict type-II band alignments and determine the absolute band offsets of these two materials with BAs. The combination of the ultra-high thermal conductivity and intrinsic p-type character of BAs, with its high electron and hole mobilities that can be further increased by tensile strain, as well as the lattice-match and the type-II band alignment with intrinsically n-type InGaN and ZnSnN2 demonstrate the potential of BAs heterostructures for electronic and optoelectronic devices.


title: Defects in Cd3As2 Epilayers Via Molecular Beam Epitaxy and Strategies for Reducing Them

abstract: Molecular beam epitaxy offers an exciting avenue for investigating the behavior of topological semimetal Cd3As2, by providing routes for doping, alloying, strain engineering, and heterostructure formation. To date, however, minimal exploration has been devoted to the impact of defects that are incorporated into epilayers due to contraints imposed by the substrate and narrow growth window. Here, we use a combination of lattice-matched ZnxCd1-xTe buffer layers, miscut substrates and broadband illumination to study how dislocations, twins and point defects influence the electron mobility of Cd3As2. A combination of defect suppression approaches produces Cd3As2 epilayers with electron mobilities upwards of 15,000 cm2/V-s at room temperature.


title: Renormalization of bi-cubic circle maps

abstract: We develop a renormalization theory for analytic homeomorphisms of the circle with two cubic critical points. We prove a renormalization hyperbolicity theorem. As a basis for the proofs, we develop complex a priori bounds for multi-critical circle maps.


title: Using causal networks to represent the targets of resource coordination

abstract: The resources framework emphasizes the potential productivity of student intuitions for constructing a canonical understanding of physics. It models learning as the progressive coordination and refinement of these resources. Yet, there is a lack of theoretical clarity about how resources should be coordinated and refined to align with canonical physics. We present causal network diagrams as a tool for representing the targets of research coordination. As an example, we compare student reasoning about projectile motion to the causal network describing that motion. We argue that the causal networks make manifest and explicit two types of resource coordination required to achieve a correct physical understanding: (i) integrating additional causal influences and mediators and (ii) using qualitative logic to draw valid inferences.


title: Towards Megajoule X-ray Lasers via Relativistic 4-Photon Cascade in Plasma

abstract: An extraordinarily efficient mechanism, operating at high laser intensities and powers, is identified for spectral transferring huge laser energies to shorter ultraviolet and x-ray wavelengths. With megajoule laser energies currently available at near-optical wavelengths, this transfer would enable megajoule x-ray lasers, a huge advance over the millijoules x-ray pulses produced now.


title: Asymptotic Giant Branch Stars in the Nearby Dwarf Galaxy Leo P

abstract: We have conducted a highly sensitive census of the evolved-star population in the metal-poor dwarf galaxy Leo P and detected four asymptotic giant branch (AGB) star candidates. Leo P is one of the best examples of a nearby analog of high-redshift galaxies because of its primitive metal content (2% of the solar value), proximity, and isolated nature, ensuring a less complicated history. Using medium-band optical photometry from the Hubble Space Telescope (HST), we have classified the AGB candidates by their chemical type. We have identified one oxygen-rich source which appears to be dusty in both the HST and Spitzer observations. Its brightness, however, suggests it may be a planetary nebula or post-AGB object. We have also identified three carbon-rich candidates, one of which may be dusty. Follow-up observations are needed to confirm the nature of these sources and to study the composition of any dust that they produce. If dust is confirmed, these stars would likely be among the most metal-poor examples of dust-producing stars known and will provide valuable insight into our understanding of dust formation at high redshift.


title: A new approach for uncovering student resources with multiple-choice questions

abstract: The traditional approach to studying student understanding presents a question and uses the student answers to make inferences about their knowledge. However, this method does not capture the range of possible alternative ideas available to students. We use a new approach, asking students to generate a plausible explanation for every choice of a multiple-choice question, to capture a range of explanations that students can generate in answering physics questions. Asking 16 students to provide explanations in this way revealed alternative possibilities for student thinking that would not have been captured if they only provided one solution. The findings show two ways these alternatives can be productive for learning physics: (i) even students who ultimately chose the wrong answer could often generate the correct explanation and (ii) many incorrect explanations contained elements of correct physical reasoning. We discuss the instructional implications of this multiple-choice questioning approach and of student alternative ideas.


title: Topologically-Guided Color Image Enhancement

abstract: Enhancement is an important step in post-processing digital images for personal use, in medical imaging, and for object recognition. Most existing manual techniques rely on region selection, similarity, and/or thresholding for editing, never really considering the topological structure of the image. In this paper, we leverage the contour tree to extract a hierarchical representation of the topology of an image. We propose 4 topology-aware transfer functions for editing features of the image using local topological properties, instead of global image properties. Finally, we evaluate our approach with grayscale and color images.


title: Freeze-in Production of Dark Matter Prior to Early Matter Domination

abstract: Freeze-out or freeze-in during a period of early matter domination can yield the correct dark matter abundance for small values of the velocity-averaged annihilation cross section, $\langle \sigma_{\rm ann} v \rangle_{\rm f} < 3 \times 10^{-26}$ cm$^3$ s$^{-1}$. However, in a generic non-standard thermal history, such a period is typically preceded by other phases. Here, we study production of dark matter in a simple post-inflationary history where a radiation-dominated phase after reheating is followed by an epoch of early matter domination. Focusing on the freeze-in regime, we show that dark matter production prior to early matter domination can dominate the relic abundance in large parts of the parameter space, including weak scale dark matter masses, and the allowed regions are highly dependent on the entire post-inflationary history. Moreover, for a very broad range of $\langle \sigma_{\rm ann} v \rangle_{\rm f}$ spanning over several decades, dark matter particles can start in chemical equilibrium early on and decouple during early matter domination, thereby rendering the relic abundance essentially independent of $\langle \sigma_{\rm ann} v \rangle_{\rm f}$. We briefly discuss connections to different observables as a possible means to test the elusive freeze-in scenario in this case.


title: New four-color optical photometry of the young stars V521 Cyg and FHO 27

abstract: Results from new BVRI photometric observations of the pre-MS stars V521 Cyg and FHO 27 collected during the period from February 2015 to July 2019 are presented. These stars are embedded in the LDN 935, dubbed "The Gulf of Mexico" in the active star-forming complex NGC 7000/IC 5070. Previous studies revealed them as variables with large amplitude brightness variations. Our new observations show that the stars continue to exhibit strong photometric variability and fading events in the brightness. On the basis of the new data received and the data available in the literature, we were able to specify the periodicity in the brightness variations of V521 Cyg. In the case of FHO 27, we found a significant peak in its periodogram corresponding to 893 days period.


title: Interpretable Word Embeddings via Informative Priors

abstract: Word embeddings have demonstrated strong performance on NLP tasks. However, lack of interpretability and the unsupervised nature of word embeddings have limited their use within computational social science and digital humanities. We propose the use of informative priors to create interpretable and domain-informed dimensions for probabilistic word embeddings. Experimental results show that sensible priors can capture latent semantic concepts better than or on-par with the current state of the art, while retaining the simplicity and generalizability of using priors.


title: Measurement and simulation of a temperature field

abstract: Measurements of surface temperature fields are used to determine the heat transfer by conduction and convection from an inhomogeneously heated metallic tube into environment. For most of the here reported measurements we use a low-cost infra-red camera Seek Thermal, but the results are also compared to measurements done with a professional Fluke Ti20 Thermal imager. The results are interpreted using two simple theoretical models which give estimates of the thermal conductivity of the tube and the heat transfer coefficient between the tube and the environment.


title: A note on pseudorandom Ramsey graphs

abstract: For fixed $s \ge 3$, we prove that if optimal $K_s$-free pseudorandom graphs exist, then the Ramsey number $r(s,t) = t^{s-1+o(1)}$ as $t \rightarrow \infty$. Our method also improves the best lower bounds for $r(C_{\ell},t)$ obtained by Bohman and Keevash from the random $C_{\ell}$-free process by polylogarithmic factors for all odd $\ell \geq 5$ and $\ell \in \{6,10\}$. For $\ell = 4$ it matches their lower bound from the $C_4$-free process. We also prove, via a different approach, that $r(C_5, t)> (1+o(1))t^{11/8}$ and $r(C_7, t)> (1+o(1))t^{11/9}$. These improve the exponent of $t$ in the previous best results and appear to be the first examples of graphs $F$ with cycles for which such an improvement of the exponent for $r(F, t)$ is shown over the bounds given by the random $F$-free process and random graphs.


title: Predicting Specificity in Classroom Discussion

abstract: High quality classroom discussion is important to student development, enhancing abilities to express claims, reason about other students' claims, and retain information for longer periods of time. Previous small-scale studies have shown that one indicator of classroom discussion quality is specificity. In this paper we tackle the problem of predicting specificity for classroom discussions. We propose several methods and feature sets capable of outperforming the state of the art in specificity prediction. Additionally, we provide a set of meaningful, interpretable features that can be used to analyze classroom discussions at a pedagogical level.


title: Prospect Theory Based Crowdsourcing for Classification in the Presence of Spammers

abstract: We consider the $M$-ary classification problem via crowdsourcing, where crowd workers respond to simple binary questions and the answers are aggregated via decision fusion. The workers have a reject option to skip answering a question when they do not have the expertise, or when the confidence of answering that question correctly is low. We further consider that there are spammers in the crowd who respond to the questions with random guesses. Under the payment mechanism that encourages the reject option, we study the behavior of honest workers and spammers, whose objectives are to maximize their monetary rewards. To accurately characterize human behavioral aspects, we employ prospect theory to model the rationality of the crowd workers, whose perception of costs and probabilities are distorted based on some value and weight functions, respectively. Moreover, we estimate the number of spammers and employ a weighted majority voting decision rule, where we assign an optimal weight for every worker to maximize the system performance. The probability of correct classification and asymptotic system performance are derived. We also provide simulation results to demonstrate the effectiveness of our approach.


title: Rates of Convergence for Large-scale Nearest Neighbor Classification

abstract: Nearest neighbor is a popular class of classification methods with many desirable properties. For a large data set which cannot be loaded into the memory of a single machine due to computation, communication, privacy, or ownership limitations, we consider the divide and conquer scheme: the entire data set is divided into small subsamples, on which nearest neighbor predictions are made, and then a final decision is reached by aggregating the predictions on subsamples by majority voting. We name this method the big Nearest Neighbor (bigNN) classifier, and provide its rates of convergence under minimal assumptions, in terms of both the excess risk and the classification instability, which are proven to be the same rates as the oracle nearest neighbor classifier and cannot be improved. To significantly reduce the prediction time that is required for achieving the optimal rate, we also consider the pre-training acceleration technique applied to the bigNN method, with proven convergence rate. We find that in the distributed setting, the optimal choice of the neighbor $k$ should scale with both the total sample size and the number of partitions, and there is a theoretical upper limit for the latter. Numerical studies have verified the theoretical findings.


title: Towards Gradual Checking of Reference Capabilities

abstract: Concurrent and parallel programming is difficult due to the presence of memory side-effects, which may introduce data races. Type qualifiers, such as reference capabilities, can remove data races by restricting sharing of mutable data. Unfortunately, reference capability languages are an all-in or nothing game, i.e., all the types must be annotated with reference capabilities. In this work in progress, we propose to mix the ideas from the reference capability literature with gradual typing, leading to gradual reference capabilities.


title: Polariton dynamics in strongly interacting quantum many-body systems

abstract: We develop a theory for light propagating in an atomic Bose-Einstein condensate in the presence of strong interactions. The resulting many-body correlations are shown to have profound effects on the optical properties of this interacting medium. For weak atom-light coupling, there is a well-defined quasiparticle, the polaron-polariton, supporting light propagation with spectral features differing significantly from the noninteracting case. The damping of the polaron-polariton depends nonmonotonically on the light-matter coupling strength, initially increasing and then decreasing. This gives rise to an interesting crossover between two quasiparticles: a bare polariton and a polaron-polariton, separated by a complex and lossy mixture of light and matter.


title: L-Sweeps: A scalable, parallel preconditioner for the high-frequency Helmholtz equation

abstract: We present the first fast solver for the high-frequency Helmholtz equation that scales optimally in parallel, for a single right-hand side. The L-sweeps approach achieves this scalability by departing from the usual propagation pattern, in which information flows in a 180 degree cone from interfaces in a layered decomposition. Instead, with L-sweeps, information propagates in 90 degree cones induced by a checkerboard domain decomposition (CDD). We extend the notion of accurate transmission conditions to CDDs and introduce a new sweeping strategy to efficiently track the wave fronts as they propagate through the CDD. The new approach decouples the subdomains at each wave front, so that they can be processed in parallel, resulting in better parallel scalability than previously demonstrated in the literature. The method has an overall O((N/p) log w) empirical run-time for N=n^d total degrees-of-freedom in a d-dimensional problem, frequency w, and p=O(n) processors. We introduce the algorithm and provide a complexity analysis for our parallel implementation of the solver. We corroborate all claims in several two- and three-dimensional numerical examples involving constant, smooth, and discontinuous wave speeds.


title: Organizing genome engineering for the gigabase scale

abstract: Engineering the entire genome of an organism enables large-scale changes in organization, function, and external interactions, with significant implications for industry, medicine, and the environment. Improvements to DNA synthesis and organism engineering are already enabling substantial changes to organisms with megabase genomes, such as Escherichia coli and Saccharomyces cerevisiae. Simultaneously, recent advances in genome-scale modeling are increasingly informing the design of metabolic networks. However, major challenges remain for integrating these and other relevant technologies into workflows that can scale to the engineering of gigabase genomes. In particular, we find that a major under-recognized challenge is coordinating the flow of models, designs, constructs, and measurements across the large teams and complex technological systems that will likely be required for gigabase genome engineering. We recommend that the community address these challenges by 1) adopting and extending existing standards and technologies for representing and exchanging information at the gigabase genomic scale, 2) developing new technologies to address major open questions around data curation and quality control, 3) conducting fundamental research on the integration of modeling and design at the genomic scale, and 4) developing new legal and contractual infrastructure to better enable collaboration across multiple institutions.


title: Generalized chi-squared detector for LTI systems with non-Gaussian noise

abstract: Previously, we derived exact relationships between the properties of a linear time-invariant control system and properties of an anomaly detector that quantified the impact an attacker can have on the system if that attacker aims to remain stealthy to the detector. A necessary first step in this process is to be able to precisely tune the detector to a desired level of performance (false alarm rate) under normal operation, typically through the selection of a threshold parameter. To-date efforts have only considered Gaussian noises. Here we generalize the approach to tune a chi-squared anomaly detector for noises with non-Gaussian distributions. Our method leverages a Gaussian Mixture Model to represent the arbitrary noise distributions, which preserves analytic tractability and provides an informative interpretation in terms of a collection of chi-squared detectors and multiple Gaussian disturbances.


title: Electro-optic entanglement source for microwave to telecom quantum state transfer

abstract: We propose an efficient microwave-photonic modulator as a resource for stationary entangled microwave-optical fields and develop the theory for deterministic entanglement generation and quantum state transfer in multi-resonant electro-optic systems. The device is based on a single crystal whispering gallery mode resonator integrated into a 3D microwave cavity. The specific design relies on a new combination of thin-film technology and conventional machining that is optimized for the lowest dissipation rates in the microwave, optical and mechanical domains. We extract important device properties from finite element simulations and predict continuous variable entanglement generation rates on the order of a Mebit/s for optical pump powers of only a few tens of microwatt. We compare the quantum state transfer fidelities of coherent, squeezed and non-Gaussian cat-states for both teleportation and direct conversion protocols under realistic conditions. Combining the unique capabilities of circuit quantum electrodynamics with the resilience of fiber optic communication could facilitate long distance solid-state qubit networks, new methods for quantum signal synthesis, quantum key distribution, and quantum enhanced detection, as well as more power-efficient classical sensing and modulation.


title: Search for Neutron Star Binaries in the Local Group Galaxies Using LISA

abstract: We discuss the prospects of LISA for detecting neutron star binaries (NSBs) in the local group galaxies such as LMC and M31. Using the recently estimated merger rate $ {\rm 1540 \, Gpc^{-3}yr^{-1}}$ and inversely applying the conventional arguments based on the $B$-band galaxy luminosities, we estimate the frequency distributions of NSBs in the local galaxies. We find that, after 10 year observation with its current design sensitivity, LISA might detect $\sim 5$ NSBs both in LMC and M31 with signal-to-noise ratios larger than 10. Some of the NSBs might be three-dimensionally localized well within LMC. These binaries will be useful for studying various topics including the origin of r-process elements.


title: Further Results on an Abstract Model for Branching and its Application to Mixed-Integer Programming

abstract: A key ingredient in branch and bound (B&B) solvers for mixed-integer programming (MIP) is the selection of branching variables since poor or arbitrary selection can affect the size of the resulting search trees by orders of magnitude. A recent article by Le Bodic and Nemhauser [Mathematical Programming, (2017)] investigated variable selection rules by developing a theoretical model of B&B trees from which they developed some new, effective scoring functions for MIP solvers. In their work, Le Bodic and Nemhauser left several open theoretical problems, solutions to which could guide the future design of variable selection rules. In this article, we first solve many of these open theoretical problems. We then implement an improved version of the model-based branching rules in SCIP 6.0, a state-of-the-art academic MIP solver, in which we observe an 11% geometric average time and node reduction on instances of the MIPLIB 2017 Benchmark Set that require large B&B trees.


title: Asynchronous Time-Parallel Method based on Laplace Transform

abstract: Laplace transform method has proved to be very efficient and easy to parallelize for the solution of time-dependent problems. However, the synchronization delay among processors implies an upper bound on the expectable acceleration factor, which leads to a lot of wasted time. In this paper, we propose an original asynchronous Laplace transform method formalized for quasilinear problems based on the well-known Gaver-Stehfest algorithm. Parallel experiments show the convergence of our new method, as well as several interesting properties compared with the classical algorithms.


title: Effective Strategies for Using Hashtags in Online Communication

abstract: The features of use of hashtags among students of Lviv were investigated. The list of optimal strategies for using these communicative tools for personal branding is determined. The effective strategy for using hashtags in online communication for the personal and company branding is considered. The results of calculation of effectiveness of hashtags related to #education is calculated. The reports of using hashtag #education in social networks is presented.


title: Critical role of the sign structure in the doped Mott insulator: Superconductivity vs. Fermi liquid

abstract: Mechanism of superconductivity (SC) in a purely interacting electron system has been one of the most challenging issues in condensed matter physics. In the BCS theory, the Landau's Fermi liquid is a normal state against which an SC instability occurs once an \emph{additional} pairing force is added. We show that in the doped Mott insulator an \emph{intrinsic} SC ground state (specifically a Luther-Emery state in a finite-doping two-leg $t$-$J$ ladder) can be directly turned into a Fermi-gas-like state by merely switching off the hidden statistical sign structure via two schemes. It points to a new pairing paradigm, %beyond the resonating-valence-bond scheme, which is an "Amperean-like pairing" with a "stringlike" pairing force as shown by an adiabatic continuity to a strong anisotropic limit of the model.


title: How much research shared on Facebook happens outside of public pages and groups? A comparison of public and private online activity around PLOS ONE papers

abstract: Despite its undisputed position as the biggest social media platform, Facebook has never entered the main stage of altmetrics research. In this study, we argue that the lack of attention by altmetrics researchers is due, in part, to the challenges in collecting Facebook data regarding activity that takes place outside of public pages and groups. We present a new method of collecting aggregate counts of shares, reactions, and comments across the platform-including users' personal timelines-and use it to gather data for all articles published between 2015 to 2017 in the journal PLOS ONE. We compare the gathered data with altmetrics collected and aggregated by Altmetric. The results show that 58.7% of papers shared on Facebook happen outside of public spaces and that, when collecting all shares, the volume of activity approximates patterns of engagement previously only observed for Twitter. Both results suggest that the role and impact of Facebook as a medium for science and scholarly communication has been underestimated. Furthermore, they emphasise the importance of openness and transparency around the collection and aggregation of altmetrics.


title: Filtering Approaches for Dealing with Noise in Anomaly Detection

abstract: The leading workhorse of anomaly (and attack) detection in the literature has been residual-based detectors, where the residual is the discrepancy between the observed output provided by the sensors (inclusive of any tampering along the way) and the estimated output provided by an observer. These techniques calculate some statistic of the residual and apply a threshold to determine whether or not to raise an alarm. To date, these methods have not leveraged the frequency content of the residual signal in making the detection problem easier, specifically dealing with the case of (e.g., measurement) noise. Here we demonstrate some opportunities to combine filtering to enhance the performance of residual-based detectors. We also demonstrate how filtering can provide a compelling alternative to residual-based methods when paired with a robust observer. In this process, we consider the class of attacks that are stealthy, or undetectable, by such filtered detection methods and the impact they can have on the system.


title: Trans-Neptunian objects found in the first four years of the Dark Energy Survey

abstract: We present a catalog of 316 trans-Neptunian bodies detected by the Dark Energy Survey (DES). These objects include 245 discoveries by DES (139 not previously published) detected in $\approx 60,000$ exposures from the first four seasons of the survey ("Y4" data). The survey covers a contiguous 5000 deg$^2$ of the southern sky in the $grizY$ optical/NIR filter set, with a typical TNO in this part of the sky being targeted by $25-30$ Y4 exposures. We describe the processes for detection of transient sources and the linkage into TNO orbits, which are made challenging by the absence of the few-hour repeat observations employed by TNO-optimized surveys. We also describe the procedures for determining detection efficiencies vs. magnitude and estimating rates of false-positive linkages. This work presents all TNOs which were detected on $\ge 6$ unique nights in the Y4 data and pass a "sub-threshold confirmation" test wherein we demand the the object be detectable in a stack of the individual images in which the orbit indicates an object should be present, but was not detected. This eliminates false positives and yields TNO detections complete to $r\lesssim 23.3$ mag with virtually no dependence on orbital properties for bound TNOs at distance $30\,{\rm AU}<d<2500\,{\rm AU}.$ The final DES TNO catalog is expected to yield $>0.3$ mag more depth, and arcs of $>4$ years for nearly all detections.


title: Fast Gradient Methods with Alignment for Symmetric Linear Systems without Using Cauchy Step

abstract: The performance of gradient methods has been considerably improved by the introduction of delayed parameters. After two and a half decades, the revealing of second-order information has recently given rise to the Cauchy-based methods with alignment, which reduce asymptotically the search spaces in smaller and smaller dimensions. They are generally considered as the state of the art of gradient methods. This paper reveals the spectral properties of minimal gradient and asymptotically optimal steps, and then suggests three fast methods with alignment without using the Cauchy step. The convergence results are provided, and numerical experiments show that the new methods provide competitive and more stable alternatives to the classical Cauchy-based methods. In particular, alignment gradient methods present advantages over the Krylov subspace methods in some situations, which makes them attractive in practice.