title: Properties of the Sheath Regions of Coronal Mass Ejections with or without Shocks from STEREO in situ Observations near 1 AU

abstract: We examine 188 coronal mass ejections (CMEs) measured by the twin STEREO spacecraft during 2007-2016 to investigate the generic features of the CME sheath and the magnetic ejecta (ME) and dependencies of average physical parameters of the sheath on the ME. We classify the MEs into three categories, focusing on whether a ME drives both a shock and sheath, or only a sheath, or neither, near 1 AU. We also reevaluate our initial classification through an automated algorithm and visual inspection. We observe that even for leading edge speeds greater than 500 km/s, 1 out of 4 MEs do not drive shocks near 1 AU. MEs driving both shocks and sheaths are the fastest and propagate in high magnetosonic solar wind, whereas MEs driving only sheaths are the slowest and propagate in low magnetosonic solar wind. Our statistical and superposed epoch analyses indicate that all physical parameters are more enhanced in the sheath regions following shocks than in sheaths without shocks. However, differences within sheaths become statistically less significant for similar driving MEs. We also find that the radial thickness of ME-driven sheaths apparently has no clear linear correlation with the speed profile and associated Mach numbers of the driver.


title: Design, upgrade and characterization of the silicon photomultiplier front-end for the AMIGA detector at the Pierre Auger Observatory

abstract: AMIGA (Auger Muons and Infill for the Ground Array) is an upgrade of the Pierre Auger Observatory to complement the study of ultra-high-energy cosmic rays (UHECR) by measuring the muon content of extensive air showers (EAS). It consists of an array of 61 water Cherenkov detectors on a denser spacing in combination with underground scintillation detectors used for muon density measurement. Each detector is composed of three scintillation modules, with 10 m$^2$ detection area per module, buried at 2.3 m depth, resulting in a total detection area of 30 m$^2$. Silicon photomultiplier sensors (SiPM) measure the amount of scintillation light generated by charged particles traversing the modules. In this paper, the design of the front-end electronics to process the signals of those SiPMs and test results from the laboratory and from the Pierre Auger Observatory are described. Compared to our previous prototype, the new electronics shows a higher performance, higher efficiency and lower power consumption, and it has a new acquisition system with increased dynamic range that allows measurements closer to the shower core. The new acquisition system is based on the measurement of the total charge signal that the muonic component of the cosmic ray shower generates in the detector.


title: Structure of the thin disk in the background of a distorted naked singularity versus a distorted static black hole

abstract: We consider the standard relativistic geometrically thin and optically thick accretion disc around a distorted static black hole and a distorted naked singularity. The distortion is the result of the existence of a static axisymmetric external distribution of matter via exercising quadrupole moments. Our main purpose of this work is to investigate whether the naked singularity models have their own observational fingerprint if they do exist. In fact, understanding the astrophysical behavior of naked singularities seems to be extremely important in the perspective of general relativity and understanding of black hole physics today.


title: Phase transitions and critical behavior in hadronic transport with a relativistic density functional equation of state

abstract: We develop a flexible, relativistically covariant parameterization of dense nuclear matter equation of state suited for inclusion in computationally demanding hadronic transport simulations. Within an implementation in the hadronic transport code SMASH, we show that effects due to bulk thermodynamic behavior are reproduced in dynamic hadronic systems, demonstrating that hadronic transport can be used to study critical behavior in dense nuclear matter, both at and away from equilibrium. We also show that two-particle correlations calculated from hadronic transport simulation data follow theoretical expectations based on the second order cumulant ratio, and constitute a clear signature of crossing the phase diagram above the critical point.


title: A data driven heuristic for rapid convergence of Scheduled Relaxation Jacobi schemes

abstract: The Scheduled Relaxation Jacobi (SRJ) method is a viable candidate as a high performance linear solver for elliptic partial differential equations (PDEs). The method greatly improves the convergence of the standard Jacobi iteration by applying a sequence of $M$ well-chosen overrelaxation and underrelaxation factors in each cycle of the algorithm to effectively attenuate the solution error. In previous work, optimal SRJ schemes (sets of relaxation factors) have been derived to accelerate convergence for specific discretizations of elliptic PDEs. In this work, we develop a family of SRJ schemes which can be applied to solve elliptic PDEs regardless of the specific discretization employed. To achieve favorable convergence, we train an algorithm to select which scheme in this family to apply at each cycle of the linear solve process, based on convergence data collected from applying these schemes to the one-dimensional Poisson equation. The automatic selection heuristic that is developed based on this limited data is found to provide good convergence for a wide range of problems.


title: On approximation of maps into real algebraic homogeneous spaces

abstract: Let X be a compact (resp. compact and nonsingular) real algebraic variety and let Y be a homogeneous space for some linear real algebraic group. We prove that a continuous (resp. C^infinity) map f:X-->Y can be approximated by regular maps in the C^o (resp. C^infinity) topology if and only if it is homotopic to a regular map. Taking Y=S^p, the unit p-dimensional sphere, we obtain solutions of several problems that have been open since the 1980's and which concern approximation of maps with values in the unit spheres. This has several consequences for approximation of maps between unit spheres. For example, we prove that for every positive integer n every C^infinity map from S^n into S^n can be approximated by regular maps in the C^infinity topology. Up to now such a result has only been known for five special values of n, namely, n=1,2,3,4 or 7.


title: Resolving nonuniform temperature distributions with single-beam absorption spectroscopy. Part II: Implementation from broadband spectra

abstract: Several past studies have described how absorption spectroscopy can be used to determine spatial temperature variations along the optical path by measuring the unique, nonlinear response to temperature of many molecular absorption transitions and performing an inversion. New laser absorption spectroscopy techniques are well-suited to this nonuniformity measurement, yet present analysis approaches use only isolated features rather than a full broadband spectral measurement. In this work, we develop a constrained spectral fitting technique called E"-binning to fit an absorption spectrum arising from a nonuniform environment. The information extracted from E"-binning is then input to the inversion approach from the previous paper in this series (Malarich and Rieker, JQSRT 107455) to determine the temperature distribution. We demonstrate this approach by using dual frequency comb laser measurements to resolve convection cells in a tube furnace. The recovered temperature distributions at each measurement height agree with an existing natural convection model. Finally, we show that for real-world measurements with noise and absorption model error, increasing the bandwidth and the number of measured absorption transitions may improve the temperature distribution precision. We make the fitting code publicly available for use with any broadband absorption measurement.


title: Empirical Performance Analysis of Conventional Deep Learning Models for Recognition of Objects in 2-D Images

abstract: Artificial Neural Networks, an essential part of Deep Learning, are derived from the structure and functionality of the human brain. It has a broad range of applications ranging from medical analysis to automated driving. Over the past few years, deep learning techniques have improved drastically - models can now be customized to a much greater extent by varying the network architecture, network parameters, among others. We have varied parameters like learning rate, filter size, the number of hidden layers, stride size and the activation function among others to analyze the performance of the model and thus produce a model with the highest performance. The model classifies images into 3 categories, namely, cars, faces and aeroplanes.


title: Invariants of Self-Intersected N-Periodics in the Elliptic Billiard

abstract: We study self-intersected N-periodics in the elliptic billiard, describing new facts about their geometry (e.g., self-intersected 4-periodics have vertices concyclic with the foci). We also check if some invariants listed in "Eighty New Invariants of N-Periodics in the Elliptic Billiard" (2020), arXiv:2004.12497, remain invariant in the self-intersected case. Toward that end, we derive explicit expressions for many low-N simple and self-intersected cases. We identify two special cases (one simple, one self-intersected) where a quantity prescribed to be invariant is actually variable.


title: Constraining velocity-dependent Lorentz/CPT-violations using Lunar Laser Ranging

abstract: The possibility for Lorentz/CPT-breaking, which is motivated by unification theories, can be systematically tested within the standard-model extension framework. In the pure gravity sector, the mass dimension 5 operators produce new Lorentz and CPT-breaking terms in the 2-body equations of motion that depend on the relative velocity of the bodies. In this Letter, we report new constraints on 15 independent SME coefficients for Lorentz/CPT-violations with mass dimension 5 using lunar laser ranging. We perform a global analysis of lunar ranging data within the SME framework using more than 26,000 normal points between 1969 and 2018. We also perform a jackknife analysis in order to provide realistic estimates of the systematic uncertainties. No deviation from Lorentz/CPT symmetries is reported. In addition, when fitting simultaneously for the 15 canonical SME coefficients for Lorentz/CPT-violations, we improve up to three orders of magnitude previous post-fit constraints from radio pulsars.


title: Context-aware Stand-alone Neural Spelling Correction

abstract: Existing natural language processing systems are vulnerable to noisy inputs resulting from misspellings. On the contrary, humans can easily infer the corresponding correct words from their misspellings and surrounding context. Inspired by this, we address the stand-alone spelling correction problem, which only corrects the spelling of each token without additional token insertion or deletion, by utilizing both spelling information and global context representations. We present a simple yet powerful solution that jointly detects and corrects misspellings as a sequence labeling task by fine-turning a pre-trained language model. Our solution outperforms the previous state-of-the-art result by 12.8% absolute F0.5 score.


title: Quantum Blackbody Thermometry

abstract: Blackbody radiation (BBR) sources are calculable radiation sources that are frequently used in radiometry, temperature dissemination, and remote sensing. Despite their ubiquity, blackbody sources and radiometers have a plethora of systematics. We envision a new, primary route to measuring blackbody radiation using ensembles of polarizable quantum systems, such as Rydberg atoms and diatomic molecules. Quantum measurements with these exquisite electric field sensors could enable active feedback, improved design, and, ultimately, lower radiometric and thermal uncertainties of blackbody standards. A portable, calibration-free Rydberg-atom physics package could also complement a variety of classical radiation detector and thermometers. The successful merger of quantum and blackbody-based measurements provides a new, fundamental paradigm for blackbody physics.


title: Functorial Semantics for Partial Theories

abstract: We provide a Lawvere-style definition for partial theories, extending the classical notion of equational theory by allowing partially defined operations. As in the classical case, our definition is syntactic: we use an appropriate class of string diagrams as terms. This allows for equational reasoning about the class of models defined by a partial theory. We demonstrate the expressivity of such equational theories by considering a number of examples, including partial combinatory algebras and cartesian closed categories. Moreover, despite the increase in expressivity of the syntax we retain a well-behaved notion of semantics: we show that our categories of models are precisely locally finitely presentable categories, and that free models exist.


title: A priori bounds for rough differential equations with a non-linear damping term

abstract: We consider a rough differential equation with a non-linear damping drift term: \begin{align*} dY(t) = - |Y|^{m-1} Y(t) dt + \sigma(Y(t)) dX(t), \end{align*} where $X$ is a branched rough path of arbitrary regularity $\alpha >0$, $m>1$ and where $\sigma$ is smooth and satisfies an $m$ and $\alpha$-dependent growth property. We show a strong a priori bound for $Y$, which includes the "coming down from infinity" property, i.e. the bound on $Y(t)$ for a fixed $t>0$ holds uniformly over all choices of initial datum $Y(0)$. The method of proof builds on recent work by Chandra, Moinat and Weber on a priori bounds for the $\phi^4$ SPDE in arbitrary subcritical dimension. A key new ingredient is an extension of the algebraic framework which permits to derive an estimate on higher order conditions of a coherent controlled rough path in terms of the regularity condition at lowest level.


title: Investigation of Edge States in Artificial Graphene Nano-Flakes

abstract: Graphene nano-flakes (GNFs) are predicted to host spin-polarized metallic edge states, which are envisioned for exploration of spintronics at the nanometer scale. To date, experimental realization of GNFs is only in its infancy because of the limitation of precise cutting or synthesizing methods at the nanometer scale. Here, we use low temperature scanning tunneling microscope (STM) to manipulate coronene molecules on a Cu(111) surface to build artificial triangular and hexagonal GNFs with either zigzag or armchair type of edges. We observe that the metallic edge states only exist in the GNFs with zigzag edge and localize at the most outside one type of the sublattice. The experimental results agree well with the tight-binding calculations. To our knowledge, our work renders the first systematic experimental confirmation of the predicated electronic properties of the GNFs.


title: Captured Small Solar System Bodies in the Ice Giant Region

abstract: This white paper advocates for the inclusion of small, captured Outer Solar system objects, found in the Ice Giant region in the next Decadal Survey. These objects include the Trojans and Irregular satellite populations of Uranus and Neptune. The captured small bodies provide vital clues as to the formation of our Solar system. They have unique dynamical situations, which any model of Solar system formation needs to explain. The major issue is that so few of these objects have been discovered, with very little information known about them. The purpose of this document is to prioritize further discovery and characterization of these objects. This will require the use of NASA and NSF facilities over the 2023 2032 decade, including additional support for analysis. This is in preparation for potential future insitu missions in the following decades.


title: Impact of a cost functional on the optimal control and the cost-effectiveness: control of a spreading infection as a case study

abstract: In applications of the optimal control theory to problems in medicine and biology, the dependency of the objective functional on the control itself is often a matter of controversy. In this paper, we explore the impact of the dependency using reasonably simple \emph{SIR} and \emph{SEIRS} epidemic models. To qualitatively compare the outcomes for different objective functionals, we apply the cost-effectiveness analysis. Our result shows that, at least for the comparatively inexpensive controls, the variation of the power at the controls in a biologically feasible range does not significantly affect the forms of the optimal controls and the corresponding optimal state solutions. Moreover, the costs and effectiveness are affected even less. At the same time, the dependency of the cost on the state variables can be very significant.


title: Hyperbolicity of maximal entropy measures for certain maps isotopic to Anosov

abstract: We prove the hyperbolicity of ergodic maximal entropy measures for a class of partially hyperbolic diffeomorphisms of $\mathbb{T}^{d}$, which have a compact two-dimensional center foliation.


title: BeyondPlanck VI. Noise characterization and modelling

abstract: We present a Bayesian method for estimating instrumental noise parameters and propagating noise uncertainties within the global BeyondPlanck Gibbs sampling framework, and apply this to Planck Low Frequency Instrument (LFI) time-ordered data. Following previous literature, we initially adopt a $1/f$ model for the noise power spectral density (PSD), but find the need for an additional lognormal component in the noise model for the 30 and 44\,GHz bands. We implement an optimal Wiener-filter (or constrained realization) gap-filling procedure to account for masked data. We then use this procedure to both estimate the gapless correlated noise in the time-domain, $n_\mathrm{corr}$, and to sample the noise PSD parameters, $\xi^n = \{\sigma_0, f_\mathrm{knee}, \alpha, A_\mathrm{p}\}$. In contrast to previous \textit{Planck} analyses, we assume piecewise stationary noise only within each pointing period (PID), not throughout the full mission, but we adopt the LFI Data Processing Center (DPC) results as priors on $\alpha$ and $f_\mathrm{knee}$. On average, we find best-fit correlated noise parameters that are mostly consistent with previous results, with a few notable exceptions. However, a detailed inspection of the time-dependent results reveals many important findings. First and foremost, we find strong evidence for statistically significant temporal variations in all noise PSD parameters, many of which are directly correlated with satellite housekeeping data. Second, while the simple $1/f$ model appears to be an excellent fit for the LFI 70 GHz channel, there is evidence for additional correlated noise not described by a $1/f$ model in the 30 and 44 GHz channels, including within the primary science frequency range of 0.1--1 Hz. (Abridged)


title: NuSTAR Observation of Energy Release in Eleven Solar Microflares

abstract: Solar flares are explosive releases of magnetic energy. Hard X-ray (HXR) flare emission originates from both hot (millions of Kelvin) plasma and nonthermal accelerated particles, giving insight into flare energy release. The Nuclear Spectroscopic Telescope ARray (NuSTAR) utilizes direct focusing optics to attain much higher sensitivity in the HXR range than that of previous indirect imagers. This paper presents eleven NuSTAR microflares from two active regions (AR 12671 on 2017 August 21, and AR 12712 on 2018 May 29). The temporal, spatial, and energetic properties of each are discussed in context with previously published HXR brightenings. They are seen to display several 'large-flare' properties, such as impulsive time profiles and earlier peaktimes in higher energy HXRs. For two events where active region background could be removed, microflare emission did not display spatial complexity: differing NuSTAR energy ranges had equivalent emission centroids. Finally, spectral fitting showed a high energy excess over a single thermal model in all events. This excess was consistent with additional higher-temperature plasma volumes in 10/11 microflares, and consistent only with an accelerated particle distribution in the last. Previous NuSTAR studies focused on one or a few microflares at a time, making this the first to collectively examine a sizable number of events. Additionally, this paper introduces an observed variation in the NuSTAR gain unique to the extremely low-livetime (<1%) regime, and establishes a correction method to be used in future NuSTAR solar spectral analysis.


title: A modeling framework for coupling plasticity with species diffusion

abstract: This paper presents a modeling framework---mathematical model and computational framework---to study the response of a plastic material due to the presence and transport of a chemical species in the host material. Such a modeling framework is important to a wide variety of problems ranging from Li-ion batteries, moisture diffusion in cementitious materials, hydrogen diffusion in metals, to consolidation of soils under severe loading-unloading regimes. The mathematical model incorporates experimental observations reported in the literature on how (elastic and plastic) material properties change because of the presence and transport of a chemical species. Also, the model accounts for one-way (transport affects the deformation but not vice versa) and two-way couplings between deformation and transport subproblems. The resulting coupled equations are not amenable to analytical solutions; so, we present a robust computational framework for obtaining numerical solutions. Given that popular numerical formulations do not produce nonnegative solutions, the computational framework uses an optimized-based nonnegative formulation that respects physical constraints (e.g., nonnegative concentrations). For completeness, we will also show the effect and propagation of the negative concentrations, often produced by contemporary transport solvers, into the overall predictions of deformation and concentration fields. Notably, anisotropy of the diffusion process exacerbates these unphysical violations. Using representative numerical examples, we will discuss how the concentration field affects plastic deformations of a degrading solid. Based on these numerical examples, we also discuss how plastic zones spread because of material degradation. To illustrate how the proposed computational framework performs, we report various performance metrics such as optimization iterations and time-to-solution.


title: Ultrabroad-Band Direct Digital Refractive Index Imaging Based on Suspended Graphene Plasmon Cavities

abstract: Mid-infrared spectroscopy is essential for chemical identification and compositional analysis, due to the existence of characteristic molecular absorption fingerprints. However, it is very challenging to determine the refractive index of an analyte at low concentrations using current photonic systems in a broad mid-infrared spectral range. We propose an imaging-based nanophotonic technique for refractive index determination. The technique is based on deeply subwavelength graphene plasmon cavities and allows for the retrieval of molecular concentration. This method features a two-dimensional array of suspended graphene plasmon cavities, in which the extremely high field enhancement and extraordinary compression of graphene plasmons can be realized simultaneously by combining shallow and deep cavities. This enables resonant unit cells to be read out in the spatial absorption pattern of the array at multiple spectral points, and the resulting information is then translated into the refractive index of the analytes. The proposed technique gives complementary information compared with the current nanophotonic techniques based on molecular absorption, including the ability of the refractive index measurement, the ultra-broadband measurable spectral range, and the small volume of the analyte required, thereby pushing the potential for refractometric sensing technologies into the infrared frequencies.


title: Utilizing Ensemble Learning for Performance and Power Modeling and Improvement of Parallel Cancer Deep Learning CANDLE Benchmarks

abstract: Machine learning (ML) continues to grow in importance across nearly all domains and is a natural tool in modeling to learn from data. Often a tradeoff exists between a model's ability to minimize bias and variance. In this paper, we utilize ensemble learning to combine linear, nonlinear, and tree-/rule-based ML methods to cope with the bias-variance tradeoff and result in more accurate models. Hardware performance counter values are correlated with properties of applications that impact performance and power on the underlying system. We use the datasets collected for two parallel cancer deep learning CANDLE benchmarks, NT3 (weak scaling) and P1B2 (strong scaling), to build performance and power models based on hardware performance counters using single-object and multiple-objects ensemble learning to identify the most important counters for improvement. Based on the insights from these models, we improve the performance and energy of P1B2 and NT3 by optimizing the deep learning environments TensorFlow, Keras, Horovod, and Python under the huge page size of 8 MB on the Cray XC40 Theta at Argonne National Laboratory. Experimental results show that ensemble learning not only produces more accurate models but also provides more robust performance counter ranking. We achieve up to 61.15% performance improvement and up to 62.58% energy saving for P1B2 and up to 55.81% performance improvement and up to 52.60% energy saving for NT3 on up to 24,576 cores.


title: Performance and Power Modeling and Prediction Using MuMMI and Ten Machine Learning Methods

abstract: In this paper, we use modeling and prediction tool MuMMI (Multiple Metrics Modeling Infrastructure) and ten machine learning methods to model and predict performance and power and compare their prediction error rates. We use a fault-tolerant linear algebra code and a fault-tolerant heat distribution code to conduct our modeling and prediction study on the Cray XC40 Theta and IBM BG/Q Mira at Argonne National Laboratory and the Intel Haswell cluster Shepard at Sandia National Laboratories. Our experiment results show that the prediction error rates in performance and power using MuMMI are less than 10% for most cases. Based on the models for runtime, node power, CPU power, and memory power, we identify the most significant performance counters for potential optimization efforts associated with the application characteristics and the target architectures, and we predict theoretical outcomes of the potential optimizations. When we compare the prediction accuracy using MuMMI with that using 10 machine learning methods, we observe that MuMMI not only results in more accurate prediction in both performance and power but also presents how performance counters impact the performance and power models. This provides some insights about how to fine-tune the applications and/or systems for energy efficiency.


title: 2020 Nobel Prize for Physics: Black holes and the Milky Way's darkest secret

abstract: This article was written at the invitation of Current Science to explain the history and Science behind this year's Nobel prize in Physics. The article is aimed at a general audience and provides a popular account and perspective on the subject of black holes.


title: Application of phenomenological Landau approach to the description of magnetic phase transitions in the ferromagnetic superconductor at ambient pressure

abstract: We study the possibility to apply phenomenological approach to the description of magnetic transitions in UGe2 with the help of Landau free energy expanded to 8-th order in magnetisation. The analysis shows that for certain values of parameters in front of M^4, and M^6 terms in the free energy there id possibility for the appearance of two successive phase transitions between the low-magnetisation and high-magnetisation phase with the same structure. We establish the relation of the parameters in Landau energy, for which the phase transition from the disordered to low-magnetisation phase is of second order as the experimental data shows. Within our approximation the transition between two magnetically ordered phases is of first order


title: Fast Convergence Algorithm for Analog Federated Learning

abstract: In this paper, we consider federated learning (FL) over a noisy fading multiple access channel (MAC), where an edge server aggregates the local models transmitted by multiple end devices through over-the-air computation (AirComp). To realize efficient analog federated learning over wireless channels, we propose an AirComp-based FedSplit algorithm, where a threshold-based device selection scheme is adopted to achieve reliable local model uploading. In particular, we analyze the performance of the proposed algorithm and prove that the proposed algorithm linearly converges to the optimal solutions under the assumption that the objective function is strongly convex and smooth. We also characterize the robustness of proposed algorithm to the ill-conditioned problems, thereby achieving fast convergence rates and reducing communication rounds. A finite error bound is further provided to reveal the relationship between the convergence behavior and the channel fading and noise. Our algorithm is theoretically and experimentally verified to be much more robust to the ill-conditioned problems with faster convergence compared with other benchmark FL algorithms.


title: A comprehensive mathematical model of a low-friction servopneumatic actuator

abstract: This paper presents a comprehensive mathematical model of a servopneumatic system, aimed at its consolidation in literature. The work exploits system's friction forces, temperature and pressure evolution, heat transfer, leakage between chambers and environment, equilibrium of cylinder forces, resistance of the pipes, mass flow rate in valve output, varying area of valve orifices and equilibrium of valve forces. Numerical simulations are performed, where system's elastic behaviour due to the air compressibility is thoroughly exploited.


title: Communication complexity of Nash equilibrium in potential games

abstract: We prove communication complexity lower bounds for (possibly mixed) Nash equilibrium in potential games. In particular, we show that finding a Nash equilibrium requires $poly(N)$ communication in two-player $N \times N$ potential games, and $2^{poly(n)}$ communication in $n$-player two-action games. To the best of our knowledge, these are the first results to demonstrate hardness in any model of (possibly mixed) Nash equilibrium in potential games.


title: Stabilization of the fluidic pinball with gradient-enriched machine learning control

abstract: We stabilize the flow past a cluster of three rotating cylinders, the fluidic pinball, with automated gradient-enriched machine learning algorithms. The control laws command the rotation speed of each cylinder in an open- and closed-loop manner. These laws are optimized with respect to the average distance from the target steady solution in three successively richer search spaces. First, stabilization is pursued with steady symmetric forcing. Second, we allow for asymmetric steady forcing. And third, we determine an optimal feedback controller employing nine velocity probes downstream. As expected, the control performance increases with every generalization of the search space. Surprisingly, both open- and closed-loop optimal controllers include an asymmetric forcing, which surpasses symmetric forcing. Intriguingly, the best performance is achieved by a combination of phasor control and asymmetric steady forcing. We hypothesize that asymmetric forcing is typical for pitchfork bifurcated dynamics of nominally symmetric configurations. Key enablers are automated machine learning algorithms augmented with gradient search: explorative gradient method for the open-loop parameter optimization and a gradient-enriched machine learning control (gMLC) for the feedback optimization. gMLC learns the control law significantly faster than previously employed genetic programming control. The gMLC source code is freely available online.


title: Mathematical Model of Hippocampal Microdialysis: Validation of in vivo Methodology

abstract: Microdialysis is a well-established method for in vivo neurochemical measurements of small molecules, with implanted concentric-design probes offering minimized tissue damage and good temporal and spatial resolution. However, the large majority of measurements do not allow the perfusate to reach equilibrium with the brain, so that inferential methods of sample concentration correction such as zero-net-flux must be used to determine actual brain extracellular fluid glucose concentrations. In order for such methods to be valid, steady-state transfer of the analyte of interest within the brain is required, but this situation has not previously been confirmed. A first-principles mathematical model of fluid flow and analyte diffusion around an implanted microdialysis probe was developed and implemented in COMSOL in order to validate the zero-net-flux approach, using measurement of extracellular brain glucose levels as a well-explored example system against which to compare the model. Results from the model accurately reproduced and predicted results from in vivo experiments. Importantly, the model predicts that the time for an implanted probe to achieve steady-state equilibrium with the surrounding extracellular fluid is on the order of one to two minutes, supporting the validity of this technique for quantitative measurement of in vivo neurochemistry.


title: Patient Recruitment Using Electronic Health Records Under Selection Bias: a Two-phase Sampling Framework

abstract: Electronic health records (EHRs) are increasingly recognized as a cost-effective resource for patient recruitment in clinical research. However, how to optimally select a cohort from millions of individuals to answer a scientific question of interest remains unclear. Consider a study to estimate the mean or mean difference of an expensive outcome. Inexpensive auxiliary covariates predictive of the outcome may often be available in patients' health records, presenting an opportunity to recruit patients selectively which may improve efficiency in downstream analyses. In this paper, we propose a two-phase sampling design that leverages available information on auxiliary covariates in EHR data. A key challenge in using EHR data for multi-phase sampling is the potential selection bias, because EHR data are not necessarily representative of the target population. Extending existing literature on two-phase sampling design, we derive an optimal two-phase sampling method that improves efficiency over random sampling while accounting for the potential selection bias in EHR data. We demonstrate the efficiency gain from our sampling design via simulation studies and an application to evaluating the prevalence of hypertension among US adults leveraging data from the Michigan Genomics Initiative, a longitudinal biorepository in Michigan Medicine.


title: Path regularity of coupled McKean-Vlasov FBSDEs

abstract: This paper establishes H\"{o}lder time regularity of solutions to coupled McKean-Vlasov forward-backward stochastic differential equations (MV-FBSDEs). This is not only of fundamental mathematical interest, but also essential for their numerical approximations. We show that a solution triple to a MV-FBSDE with Lipschitz coefficients is 1/2-H\"{o}lder continuous in time in the $L^p$-norm provided that it admits a Lipschitz decoupling field. Special examples include decoupled MV-FBSDEs, coupled MV-FBSDEs with a small time horizon and coupled stochastic Pontryagin systems arsing from mean field control problems.


title: Domain-Level Explainability -- A Challenge for Creating Trust in Superhuman AI Strategies

abstract: For strategic problems, intelligent systems based on Deep Reinforcement Learning (DRL) have demonstrated an impressive ability to learn advanced solutions that can go far beyond human capabilities, especially when dealing with complex scenarios. While this creates new opportunities for the development of intelligent assistance systems with groundbreaking functionalities, applying this technology to real-world problems carries significant risks and therefore requires trust in their transparency and reliability. With superhuman strategies being non-intuitive and complex by definition and real-world scenarios prohibiting a reliable performance evaluation, the key components for trust in these systems are difficult to achieve. Explainable AI (XAI) has successfully increased transparency for modern AI systems through a variety of measures, however, XAI research has not yet provided approaches enabling domain level insights for expert users in strategic situations. In this paper, we discuss the existence of superhuman DRL-based strategies, their properties, the requirements and challenges for transforming them into real-world environments, and the implications for trust through explainability as a key technology.


title: Integration of Fully-Actuated Multirotors into Real-World Applications

abstract: The introduction of fully-actuated multirotors has opened the door to new possibilities and more efficient solutions to many real-world applications. However, their integration had been slower than expected, partly due to the need for new tools to take full advantage of these robots. As far as we know, all the groups currently working on the fully-actuated multirotors develop new full-pose (6-D) tools and methods to use their robots, which is inefficient, time-consuming, and requires many resources. We propose a way of bridging the gap between the tools already available for underactuated robots and the new fully-actuated vehicles. The approach can extend the existing underactuated flight controllers to support the fully-actuated robots, or enhance the existing fully-actuated controllers to support existing underactuated flight stacks. We introduce attitude strategies that work with the underactuated controllers, tools, planners and remote control interfaces, all while allowing taking advantage of the full actuation. Moreover, new methods are proposed that can properly handle the limited lateral thrust suffered by many fully-actuated UAV designs. The strategies are lightweight, simple, and allow rapid integration of the available tools with these new vehicles for the fast development of new real-world applications. The real experiments on our robots and simulations on several UAV architectures with different underlying controller methods show how these strategies can be utilized to extend existing flight controllers for fully-actuated applications. We have provided the source code for the PX4 firmware enhanced with our proposed methods to showcase an example flight controller for underactuated multirotors that can be modified to seamlessly support fully-actuated vehicles while retaining the rest of the flight stack unchanged.


title: Nonlocal superconducting quantum interference device

abstract: Superconducting quantum interference devices (SQUIDs) that incorporate two superconductor/insulator/superconductor (SIS) Josephson junctions in a closed loop form the core of some of the most sensitive detectors of magnetic and electric fields currently available. SQUIDs in these applications are typically operated with a finite voltage which generates microwave radiation through the ac Josephson effect. This radiation may impact the system being measured. We describe here a SQUID in which the Josephson junctions are formed from strips of normal metal (N) in good electrical contact with the superconductor (S). Such SNS SQUIDs can be operated under a finite voltage bias with performance comparable or potentially better than conventional SIS SQUIDs. However, they also permit a novel mode of operation that is based on the unusual interplay of quasiparticle currents and supercurrents in the normal metal of the Josephson junction. The new method allows measurements of the flux dependence of the critical current of the SNS SQUID without applying a finite voltage bias across the SNS junction, enabling sensitive flux detection without generating microwave radiation.


title: Mean perimeter and area of the convex hull of a planar Brownian motion in the presence of resetting

abstract: We compute exactly the mean perimeter and the mean area of the convex hull of a $2$-d Brownian motion of duration $t$ and diffusion constant $D$, in the presence of resetting to the origin at a constant rate $r$. We show that for any $t$, the mean perimeter is given by $\langle L(t)\rangle= 2 \pi \sqrt{\frac{D}{r}}\, f_1(rt)$ and the mean area is given by $\langle A(t) \rangle= 2\pi\frac{D}{r}\, f_2(rt)$ where the scaling functions $f_1(z)$ and $f_2(z)$ are computed explicitly. For large $t\gg 1/r$, the mean perimeter grows extremely slowly as $\langle L(t)\rangle \propto \ln (rt)$ with time. Likewise, the mean area also grows slowly as $\langle A(t)\rangle \propto \ln^2(rt)$ for $t\gg 1/r$. Our exact results indicate that the convex hull, in the presence of resetting, approaches a circular shape at late times. Numerical simulations are in perfect agreement with our analytical predictions.

categories: cond-mat.stat-mech math-ph math.MP math.PR


title: Predicting laboratory earthquakes with machine learning

abstract: The successful prediction of earthquakes is one of the holy grails in Earth Sciences. Traditional predictions use statistical information on recurrence intervals, but those predictions are not accurate enough. In a recent paper, a machine learning approach was proposed and applied to data of laboratory earthquakes. The machine learning algorithm utilizes continuous measurements of radiated energy through acoustic emissions and the authors were able to successfully predict the timing of laboratory earthquakes. Here, we reproduced their model which was applied to a gouge layer of glass beads and applied it to a data set obtained using a gouge layer of salt. In this salt experiment different load point velocities were set, leading to variable recurrence times. The machine learning technique we use is called random forest and uses the acoustic emissions during the interseismic period. The random forest model succeeds in making a relatively reliable prediction for both materials, also long before the earthquake. Apparently there is information in the data on the timing of the next earthquake throughout the experiment. For glass beads energy is gradually and increasingly released whereas for salt energy is only released during precursor activity, therefore the important features used in the prediction are different. We interpret the difference in results to be due to the different micromechanics of slip. The research shows that a machine learning approach can reveal the presence of information in the data on the timing of unstable slip events (earthquakes). Further research is needed to identify the responsible micromechanical processes which might be then be used to extrapolate to natural conditions.


title: A Formal Model for Adaptive Free Choice in Complex Systems

abstract: In this article, I develop a formal model of free will for complex systems based on emergent properties and adaptive selection. The model is based on a process ontology in which a free choice is a singular process that takes a system from one macrostate to another. I quantify the model by introducing a formal measure of the `freedom' of a singular choice. The `free will' of a system, then, is emergent from the aggregate freedom of the choice processes carried out by the system. The focus in this model is on the actual choices themselves viewed in the context of processes. That is, the nature of the system making the choices is not considered. Nevertheless, my model does not necessarily conflict with models that are based on internal properties of the system. Rather it takes a behavioral approach by focusing on the externalities of the choice process.


title: Disassemblable Fieldwork CT Scanner Using a 3D-printed Calibration Phantom

abstract: The use of computed tomography (CT) imaging has become of increasing interest to academic areas outside of the field of medical imaging and industrial inspection, e.g., to biology and cultural heritage research. The pecularities of these fields, however, sometimes require that objects need to be imaged on-site, e.g., in field-work conditions or in museum collections. Under these circumstances, it is often not possible to use a commercial device and a custom solution is the only viable option. In order to achieve high image quality under adverse conditions, reliable calibration and trajectory reproduction are usually key requirements for any custom CT scanning system. Here, we introduce the construction of a low-cost disassemblable CT scanner that allows calibration even when trajectory reproduction is not possible due to the limitations imposed by the project conditions. Using 3D-printed in-image calibration phantoms, we compute a projection matrix directly from each captured X-ray projection. We describe our method in detail and show successful tomographic reconstructions of several specimen as proof of concept.


title: Why is the universe comprehensible?

abstract: Why is the universe comprehensible? How is it that we can come to know its regularities well-enough to exploit them for our own gain? In this essay I argue that the nature of our comprehension lies in the mutually agreed upon methodology we use to attain it and on the basic stability of the universe. But I also argue that the very act of comprehension itself places constraints on what we can comprehend by forcing us to establish a context for our knowledge. In this way the universe has managed to conspire to make itself objectively comprehensible to subjective observers.


title: Symbolically Solving Partial Differential Equations using Deep Learning

abstract: We describe a neural-based method for generating exact or approximate solutions to differential equations in the form of mathematical expressions. Unlike other neural methods, our system returns symbolic expressions that can be interpreted directly. Our method uses a neural architecture for learning mathematical expressions to optimize a customizable objective, and is scalable, compact, and easily adaptable for a variety of tasks and configurations. The system has been shown to effectively find exact or approximate symbolic solutions to various differential equations with applications in natural sciences. In this work, we highlight how our method applies to partial differential equations over multiple variables and more complex boundary and initial value conditions.


title: Capillary and Viscous Fracturing During Drainage in Porous Media

abstract: Detailed understanding of the couplings between fluid flow and solid deformation in porous media is crucial for the development of novel technologies relating to a wide range of geological and biological processes. A particularly challenging phenomenon that emerges from these couplings is the transition from fluid invasion to fracturing during multiphase flow. Previous studies have shown that this transition is highly sensitive to fluid flow rate, capillarity, and the structural properties of the porous medium. However, a comprehensive characterization of the relevant fluid flow and material failure regimes does not exist. Here, we used our newly developed Multiphase Darcy-Brinkman-Biot framework to examine the transition from drainage to material failure during viscously-stable multiphase flow in soft porous media in a broad range of flow, wettability, and solid rheology conditions. We demonstrate the existence of three distinct material failure regimes controlled by non-dimensional numbers that quantify the balance of viscous, capillary, and structural forces in the porous medium.


title: A theoretical scenario for Galactic RR Lyrae in the Gaia database: constraints on the parallax offset

abstract: On the basis of an extended set of nonlinear convective RR Lyrae pulsation models we derive the first theoretical light curves in the Gaia bands and the corresponding intensity-weighted mean magnitudes and pulsation amplitudes. The effects of chemical com-position on the derived Bailey diagrams in the Gaia filters are discussed for both Fundamental and First Overtone mode pulsators. The inferred mean magnitudes and colors are used to derive the first theoretical Period-Wesenheit relations for RR Lyrae in the Gaia filters. The application of the theoretical Period-Wesenheit relations for both the Fundamental and FirstOvertone mode to Galactic RR Lyrae in the Gaia Data Release 2 database and complementary information on individual metal abundances, allows us to derive theoretical estimates of their individual parallaxes. These results are compared with the astrometric solutions to conclude that a very small offset, consistent with zero, is required in order to reconcile the predicted distances with Gaia results.


title: Phonons and thermal transport in Si/SiO$_2$ multishell nanotubes: Atomistic study

abstract: Thermal transport in the Si/SiO$_2$ multishell nanotubes is investigated theoretically. The phonon energy spectra are obtained using the atomistic Lattice Dynamics approach. Thermal conductivity is calculated using the Boltzmann transport equation within the relaxation time approximation. Redistribution of the vibrational spectra in multishell nanotubes leads to a decrease of the phonon group velocity and the thermal conductivity as compared to homogeneous Si nanowires. Phonon scattering on the Si/SiO$_2$ interfaces is another key factor of strong reduction of the thermal conductivity in these structures (down to 0.2 W/mK at room temperature). We demonstrate that phonon thermal transport in Si/SiO$_2$ nanotubes can be efficiently suppressed by a proper choice of nanotube's geometrical parameters: lateral cross-section, thickness and number of shells.


title: Gauge field theory without groups

abstract: Non-standard topics underlying a partly original approach to gauge field theory are concisely introduced, expressing ideas that were broached in several papers and, eventually, exposed in an organized form in a recently published book. By proposing a change of perspective about the roles and relative importance of several notions, this approach seeks to obtain an overall clarification of foundational matters. In particular, by consistently relying on natural differential geometry, the role of groups is shown to be downgradable to secondary.


title: Comment on Demographic gaps or preparation gaps?: The large impact of incoming preparation on performance of students in introductory physics

abstract: A recent paper by Salehi et al. claims that the differences found between major demographic groups on scores in introductory college physics tests are due to differences in pre-college preparation. No evidence is produced, however, to show that preparation differences are more causally important than any other differences. In one case, the male/female difference, the paper actually provides evidence that preparation gaps are unimportant.


title: Trajectory Prediction in Autonomous Driving with a Lane Heading Auxiliary Loss

abstract: Predicting a vehicle's trajectory is an essential ability for autonomous vehicles navigating through complex urban traffic scenes. Bird's-eye-view roadmap information provides valuable information for making trajectory predictions, and while state-of-the-art models extract this information via image convolution, auxiliary loss functions can augment patterns inferred from deep learning by further encoding common knowledge of social and legal driving behaviors. Since human driving behavior is inherently multimodal, models which allow for multimodal output tend to outperform single-prediction models on standard metrics. We propose a loss function which enhances such models by enforcing expected driving rules on all predicted modes. Our contribution to trajectory prediction is twofold; we propose a new metric which addresses failure cases of the off-road rate metric by penalizing trajectories that oppose the ascribed heading (flow direction) of a driving lane, and we show this metric to be differentiable and therefore suitable as an auxiliary loss function. We then use this auxiliary loss to extend the the standard multiple trajectory prediction (MTP) and MultiPath models, achieving improved results on the nuScenes prediction benchmark by predicting trajectories which better conform to the lane-following rules of the road.


title: Cognitive RF-FSO Fronthaul Assignment in Cell-Free and User-Centric mMIMO Networks

abstract: Cell-free massive MIMO (CF-mMIMO) network and its user-centric (UC) version are considered as promising techniques for the next generations of wireless networks. However, fronthaul and backhaul assignments are challenging issues in these networks. In this paper, energy efficiencies of uplink transmission for the CF- and UC-mMIMO networks are studied, wherein access points (APs) are connected to aggregation nodes (ANs) through radio frequency (RF) and/or free-space optic (FSO) fronthauls, and the ANs are connected to a central processing unit via fiber backhauls. The achievable data rates are derived by taking into account the effects of hardware non-ideality at the APs and ANs, FSO alignment and weather conditions. To have a robust and energy-efficient network, especially in the presence of FSO misalignment and adverse weather conditions, firstly, a cognitive RF--FSO fronthaul assignment algorithm is proposed at the cost of sharing the available RF bandwidth between the access and fronthaul links. Then, optimal power allocations at the users and APs are investigated, and two analytical approaches are proposed to solve the non-convex optimization problem. Through numerical results, we have discussed how utilizing the cognitive RF--FSO fronthaul assignment achieves higher energy efficiency compared to that of FSO-only, RF-only, or simultaneously using RF and FSO fronthaul links, e.g., achieving up to $198\%$ higher energy efficiency under unfavorable conditions. Moreover, the effects of FSO misalignment, weather conditions, and power allocations on the performances of the CF- and UC-mMIMO networks are discussed.


title: Discrete Cosine Transform Based Causal Convolutional Neural Network for Drift Compensation in Chemical Sensors

abstract: Sensor drift is a major problem in chemical sensors that requires addressing for reliable and accurate detection of chemical analytes. In this paper, we develop a causal convolutional neural network (CNN) with a Discrete Cosine Transform (DCT) layer to estimate the drift signal. In the DCT module, we apply soft-thresholding nonlinearity in the transform domain to denoise the data and obtain a sparse representation of the drift signal. The soft-threshold values are learned during training. Our results show that DCT layer-based CNNs are able to produce a slowly varying baseline drift signal. We train the CNN on synthetic data and test it on real chemical sensor data. Our results show that we can have an accurate and smooth drift estimate even when the observed sensor signal is very noisy.


title: Clustering microbiome data using mixtures of logistic normal multinomial models

abstract: Discrete data such as counts of microbiome taxa resulting from next-generation sequencing are routinely encountered in bioinformatics. Taxa count data in microbiome studies are typically high-dimensional, over-dispersed, and can only reveal relative abundance therefore being treated as compositional. Analyzing compositional data presents many challenges because they are restricted on a simplex. In a logistic normal multinomial model, the relative abundance is mapped from a simplex to a latent variable that exists on the real Euclidean space using the additive log-ratio transformation. While a logistic normal multinomial approach brings in flexibility for modeling the data, it comes with a heavy computational cost as the parameter estimation typically relies on Bayesian techniques. In this paper, we develop a novel mixture of logistic normal multinomial models for clustering microbiome data. Additionally, we utilize an efficient framework for parameter estimation using variational Gaussian approximations (VGA). Adopting a variational Gaussian approximation for the posterior of the latent variable reduces the computational overhead substantially. The proposed method is illustrated on simulated and real datasets.


title: Waring's Problem For Locally Nilpotent Groups: The Case of Discrete Heisenberg Groups

abstract: Kamke \cite{Kamke1921} solved an analog of Waring's problem with $n$th powers replaced by integer-valued polynomials. Larsen and Nguyen \cite{LN2019} explored the view of algebraic groups as a natural setting for Waring's problem. This paper applies the theory of polynomial maps and polynomial sequences in locally nilpotent groups developed in previous work \cite{Hu2020} to solve an analog of Waring's problem for the general discrete Heisenberg groups $H_{2n+1}(\mathbb{Z})$ for any integer $n\ge1$.


title: Fibers are not (P)Threads: The Case for Loose Coupling of Asynchronous Programming Models and MPI Through Continuations

abstract: Asynchronous programming models (APM) are gaining more and more traction, allowing applications to expose the available concurrency to a runtime system tasked with coordinating the execution. While MPI has long provided support for multi-threaded communication and non-blocking operations, it falls short of adequately supporting APMs as correctly and efficiently handling MPI communication in different models is still a challenge. Meanwhile, new low-level implementations of light-weight, cooperatively scheduled execution contexts (fibers, aka user-level threads (ULT)) are meant to serve as a basis for higher-level APMs and their integration in MPI implementations has been proposed as a replacement for traditional POSIX thread support to alleviate these challenges. In this paper, we first establish a taxonomy in an attempt to clearly distinguish different concepts in the parallel software stack. We argue that the proposed tight integration of fiber implementations with MPI is neither warranted nor beneficial and instead is detrimental to the goal of MPI being a portable communication abstraction. We propose MPI Continuations as an extension to the MPI standard to provide callback-based notifications on completed operations, leading to a clear separation of concerns by providing a loose coupling mechanism between MPI and APMs. We show that this interface is flexible and interacts well with different APMs, namely OpenMP detached tasks, OmpSs-2, and Argobots.


title: Second Generation Readout For Large Format Photon Counting Microwave Kinetic Inductance Detectors

abstract: We present the development of a second generation digital readout system for photon counting microwave kinetic inductance detector (MKID) arrays operating in the optical and near-IR wavelength bands. Our system retains much of the core signal processing architecture from the first generation system, but with a significantly higher bandwidth, enabling readout of kilopixel MKID arrays. Each set of readout boards is capable of reading out 1024 MKID pixels multiplexed over 2 GHz of bandwidth; two such units can be placed in parallel to read out a full 2048 pixel microwave feedline over a 4 -- 8 GHz band. As in the first generation readout, our system is capable of identifying, analyzing, and recording photon detection events in real time with a time resolution of order a few microseconds. Here, we describe the hardware and firmware, and present an analysis of the noise properties of the system. We also present a novel algorithm for efficiently suppressing IQ mixer sidebands to below -30 dBc.


title: Nuclearity for partial crossed products by exact discrete groups

abstract: We study partial actions of exact discrete groups on C*-algebras. We show that the partial crossed product of a commutative C*-algebra by an exact discrete group is nuclear whenever the full and reduced partial crossed products coincide. This generalises a result by Matsumura in the context of global actions. In general, we prove that a partial action of an exact discrete group on a C*-algebra $A$ has Exel's approximation property if and only if the full and reduced partial crossed products associated to the diagonal partial action on $A\otimes_{\max} A^\mathrm{op}$ coincide. We apply our results to show that the reduced semigroup C*-algebra $\mathrm{C}^*_{\lambda}(P)$ of a submonoid of an exact discrete group is nuclear if the left regular representation on $\ell^2(P)$ is an isomorphism between the full and reduced C*-algebras. We also show that nuclearity is equivalent to the weak containment property in the case of C*-algebras associated to separated graphs.


title: Greedy Motzkin-Kaczmarz methods for solving linear systems

abstract: The famous greedy randomized Kaczmarz (GRK) method uses the greedy selection rule on maximum distance to determine a subset of the indices of working rows. In this paper, with the greedy selection rule on maximum residual, we propose the greedy randomized Motzkin-Kaczmarz (GRMK) method for linear systems. The block version of the new method is also presented. We analyze the convergence of the two methods and provide the corresponding convergence factors. Extensive numerical experiments show that the GRMK method has almost the same performance as the GRK method for dense matrices and the former performs better in computing time for some sparse matrices, and the block versions of the GRMK and GRK methods always have almost the same performance.


title: Block sampling Kaczmarz-Motzkin methods for consistent linear systems

abstract: The sampling Kaczmarz-Motzkin (SKM) method is a generalization of the randomized Kaczmarz and Motzkin methods. It first samples some rows of coefficient matrix randomly to build a set and then makes use of the maximum violation criterion within this set to determine a constraint. Finally, it makes progress by enforcing this single constraint. In this paper, on the basis of the framework of the SKM method and considering the greedy strategies, we present two block sampling Kaczmarz-Motzkin methods for consistent linear systems. Specifically, we also first sample a subset of rows of coefficient matrix and then determine an index in this set using the maximum violation criterion. Unlike the SKM method, in the rest of the block methods, we devise different greedy strategies to build index sets. Then, the new methods make progress by enforcing the corresponding multiple constraints simultaneously. Theoretical analyses demonstrate that these block methods converge at least as quickly as the SKM method, and numerical experiments show that, for the same accuracy, our methods outperform the SKM method in terms of the number of iterations and computing time.


title: An exploratory assessment of a multidimensional healthcare and economic data on COVID-19 in Nigeria

abstract: The coronavirus disease of 2019 (COVID-19) is a pandemic that is ravaging Nigeria and the world at large. This data article provides a dataset of daily updates of COVID-19 as reported online by the Nigeria Centre for Disease Control (NCDC) from February 27, 2020 to September 29, 2020. The data were obtained through web scraping from different sources and it includes some economic variables such as the Nigeria budget for each state in 2020, population estimate, healthcare facilities, and the COVID-19 laboratories in Nigeria. The dataset has been processed using the standard of the FAIR data principle which encourages its findability, accessibility, interoperability, and reusability and will be relevant to researchers in different fields such as Data Science, Epidemiology, Earth Modelling, and Health Informatics.


title: Adversarial Image Color Transformations in Explicit Color Filter Space

abstract: Deep Neural Networks have been shown to be vulnerable to adversarial images. Conventional attacks strive for indistinguishable adversarial images with strictly restricted perturbations. Recently, researchers have moved to explore distinguishable yet non-suspicious adversarial images and demonstrated that color transformation attacks are effective. In this work, we propose Adversarial Color Filter (AdvCF), a novel color transformation attack that is optimized with gradient information in the parameter space of a simple color filter. In particular, our color filter space is explicitly specified so that we are able to provide a systematic analysis of model robustness against adversarial color transformations, from both the attack and defense perspectives. In contrast, existing color transformation attacks do not offer the opportunity for systematic analysis due to the lack of such an explicit space. We further demonstrate the effectiveness of our AdvCF in fooling image classifiers and also compare it with other color transformation attacks regarding their robustness to defenses and image acceptability through an extensive user study. We also highlight the human-interpretability of AdvCF and show its superiority over the state-of-the-art human-interpretable color transformation attack on both image acceptability and efficiency. Additional results provide interesting new insights into model robustness against AdvCF in another three visual tasks.


title: CNN-based driving of block partitioning for intra slices encoding

abstract: This paper provides a technical overview of a deep-learning-based encoder method aiming at optimizing next generation hybrid video encoders for driving the block partitioning in intra slices. An encoding approach based on Convolutional Neural Networks is explored to partly substitute classical heuristics-based encoder speed-ups by a systematic and automatic process. The solution allows controlling the trade-off between complexity and coding gains, in intra slices, with one single parameter. This algorithm was proposed at the Call for Proposals of the Joint Video Exploration Team (JVET) on video compression with capability beyond HEVC. In All Intra configuration, for a given allowed topology of splits, a speed-up of $\times 2$ is obtained without BD-rate loss, or a speed-up above $\times 4$ with a loss below 1\% in BD-rate.


title: Demonstration of a MOT in a Sub-Millimeter Membrane Hole

abstract: We demonstrate the generation of a cold-atom ensemble within a sub-millimeter diameter hole in a transparent membrane, a so-called "membrane MOT". With a sub-Doppler cooling process, the atoms trapped by the membrane MOT are cooled down to 10 uK. The atom number inside the unbridged/bridged membrane hole is about 10^4 to 10^5, and the 1/e^2-diameter of the MOT cloud is about 180 um for a 400 um-diameter membrane hole. Such a membrane device can, in principle, efficiently load cold atoms into the evanescent-field optical trap generated by the suspended membrane waveguide for strong atom-light interaction and provide the capability of sufficient heat dissipation at the waveguide. This represents a key step toward the photonic atom trap integrated platform (ATIP).


title: The Uncertain Shape of Grey Swans: Extreme Value Theory with Uncertain Threshold

abstract: Extreme Value Theory (EVT) is one of the most commonly used approaches in finance for measuring the downside risk of investment portfolios, especially during financial crises. In this paper, we propose a novel approach based on EVT called Uncertain EVT to improve its forecast accuracy and capture the statistical characteristics of risk beyond the EVT threshold. In our framework, the extreme risk threshold, which is commonly assumed a constant, is a dynamic random variable. More precisely, we model and calibrate the EVT threshold by a state-dependent hidden variable, called Break-Even Risk Threshold (BRT), as a function of both risk and ambiguity. We will show that when EVT approach is combined with the unobservable BRT process, the Uncertain EVT's predicted VaR can foresee the risk of large financial losses, outperforms the original EVT approach out-of-sample, and is competitive to well-known VaR models when back-tested for validity and predictability.


title: Study of genetic mutations and dynamic spread of SARS-CoV-2 pandemic and prediction of its evolution according to the SIR model

abstract: In this work, we aim to study that the dynamics behavior for cumulative number of SARS-CoV-2 pandemic can provide information on the overall behavior of the spread over daily time.The cumulative data can be synthesized in an empirical form obtained from a Susceptible-Infected-Recovered (SIR) model previously studied on a Euclidean network. From the study we carried out, we can conclude that the SIR model on the Euclidean network can reproduce data from several countries with a deviation of precision for given parameter values. This gives an insight into the different agents that influence the behavior of SARS-CoV-2 especially during the virus mutation period. We are thus trying to analyze the effect of genetic mutations in different countries, and how a specific mutation can make the virus more contagious.


title: When Should We (Not) Interpret Linear IV Estimands as LATE?

abstract: In this paper I revisit the interpretation of the linear instrumental variables (IV) estimand as a weighted average of conditional local average treatment effects (LATEs). I focus on a practically relevant situation in which additional covariates are required for identification while the reduced-form and first-stage regressions implicitly restrict the effects of the instrument to be homogeneous, and are thus possibly misspecified. I show that the weights on some conditional LATEs are negative and the IV estimand is no longer interpretable as a causal effect under a weaker version of monotonicity, i.e. when there are compliers but no defiers at some covariate values and defiers but no compliers elsewhere. The problem of negative weights disappears in the overidentified specification of Angrist and Imbens (1995) and in an alternative method, termed "reordered IV," that I also outline. Even if all weights are positive, the IV estimand in the just identified specification is not interpretable as the unconditional LATE parameter unless the groups with different values of the instrument are roughly equal sized. I illustrate my findings in an application to causal effects of college education using the college proximity instrument. The benchmark estimates suggest that college attendance yields earnings gains of about 60 log points, which is well outside the range of estimates in the recent literature. I demonstrate, however, that this result is driven by the presence of negative weights. Corrected estimates indicate that attending college causes earnings to be roughly 20% higher.


title: Data-driven Analysis of Relight variability of Jet Fuels induced by Turbulence

abstract: For safety purposes, reliable reignition of aircraft engines in the event of flame blow-out is a critical requirement. Typically, an external ignition source in the form of a spark is used to achieve a stable flame in the combustor. However, such forced turbulent ignition may not always successfully relight the combustor, mainly because the state of the combustor cannot be precisely determined. Uncertainty in the turbulent flow inside the combustor, inflow conditions, and spark discharge characteristics can lead to variability in sparking outcomes even for nominally identical operating conditions. Prior studies have shown that of all the uncertain parameters, turbulence is often dominant and can drastically alter ignition behavior. For instance, even when different fuels have similar ignition delay times, their ignition behavior in practical systems can be completely different. In practical operating conditions, it is challenging to understand why ignition fails and how much variation in outcomes can be expected. The focus of this work is to understand relight variability induced by turbulence for two different aircraft fuels, namely Jet-A and a variant named C1. A detailed, previously developed simulation approach is used to generate a large number of successful and failed ignition events. Using this data, the cause of misfire is evaluated based on a discriminant analysis that delineates the difference between turbulent initial conditions that lead to ignition or failure. It was determined that the cause of ignition failure is different for the two fuels. While it was found that Jet-A is influenced by fuel entrainment, C1 was found to be more sensitive to small scale turbulence features. A larger variability is found in the ignition modes of C1, which can be subject to extreme events induced by kernel breakdown.


title: Jet measurements at CMS

abstract: Recent measurements based on jet production in high energy pp collisions at the CERN LHC with the CMS detector are reported. Specifically, the measurement of the inclusive jet production cross section as a function of the anti-$k_t$ distance parameter $R$ divided by the jet cross section at $R = 0.4$, is discussed. The cross section ratio is sensitive to various perturbative and non-perturbative treatments of the jet formation process in quantum chromodynamics (QCD). In the second study, the regimes of validity of the parton shower and matrix element approaches are tested on multijet events in pp collisions at $\sqrt{s} = 8$ and $13$ TeV. Special attention is given to the second and third leading $p_\text{T}$ jets and their momenta and rapidity-azimuth correlations. Finally, studies of dijet production where the two leading $p_T$ jets are separated by a large pseudorapidity interval void of charged particles is presented. This signature is expected from hard color singlet exchange (two-gluon exchange in perturbative QCD). The latter can be treated with perturbative QCD techniques based on the Balitsky-Fadin-Kuraev-Lipatov (BFKL) evolution equations.


title: Robust Policies via Mid-Level Visual Representations: An Experimental Study in Manipulation and Navigation

abstract: Vision-based robotics often separates the control loop into one module for perception and a separate module for control. It is possible to train the whole system end-to-end (e.g. with deep RL), but doing it "from scratch" comes with a high sample complexity cost and the final result is often brittle, failing unexpectedly if the test environment differs from that of training. We study the effects of using mid-level visual representations (features learned asynchronously for traditional computer vision objectives), as a generic and easy-to-decode perceptual state in an end-to-end RL framework. Mid-level representations encode invariances about the world, and we show that they aid generalization, improve sample complexity, and lead to a higher final performance. Compared to other approaches for incorporating invariances, such as domain randomization, asynchronously trained mid-level representations scale better: both to harder problems and to larger domain shifts. In practice, this means that mid-level representations could be used to successfully train policies for tasks where domain randomization and learning-from-scratch failed. We report results on both manipulation and navigation tasks, and for navigation include zero-shot sim-to-real experiments on real robots.


title: Mid-infrared single photon detector with superconductor Mo$_{80}$Si$_{20}$ nanowire

abstract: A mid-infrared single photon detector (MIR-SNSPD) was reported based on 30 nm-wide superconductor molybdenum silicide nanowires in this work. Saturated quantum efficiencies (QEs) were achieved at the wavelength ranging from 1.55 to 5.07 micrometer in experiments. At the same time, the intrinsic dark count rate (DCR) was below 100 cps. Thus, this device produced a noise equivalent power (NEP) of 4.5 * 10-19 W/sqrt(Hz). The results provide the foundation of developing 10 micrometer-SNSPD for the applications of infrared astronomy observation.

categories: physics.app-ph astro-ph.IM physics.ins-det quant-ph


title: Deep Reinforcement Learning of Transition States

abstract: Combining reinforcement learning (RL) and molecular dynamics (MD) simulations, we propose a machine-learning approach (RL$^\ddag$) to automatically unravel chemical reaction mechanisms. In RL$^\ddag$, locating the transition state of a chemical reaction is formulated as a game, where a virtual player is trained to shoot simulation trajectories connecting the reactant and product. The player utilizes two functions, one for value estimation and the other for policy making, to iteratively improve the chance of winning this game. We can directly interpret the reaction mechanism according to the value function. Meanwhile, the policy function enables efficient sampling of the transition paths, which can be further used to analyze the reaction dynamics and kinetics. Through multiple experiments, we show that RL{\ddag} can be trained tabula rasa hence allows us to reveal chemical reaction mechanisms with minimal subjective biases.


title: Calculable microscopic theory for $^{12}$C($\alpha, \gamma$)$^{16}$O cross section near Gamow window

abstract: $^{12}$C$(\alpha, \gamma)^{16}$O radiative-capture process is a key reaction to produce the element of oxygen in stars. Measuring the cross section near the Gamow window is extremely hard because it is too small. To make a theoretical contribution towards resolving the long-standing problem, I present a microscopic formulation that aims at providing all materials needed to calculate the cross section. The states of $^{12}$C and $^{16}$O relevant to the reaction are respectively described with fully microscopic 3 $\alpha$-particle and 4 $\alpha$-particle configurations, in which the relative motion among the $\alpha$ particles is expanded in terms of correlated Gaussian basis functions. The configuration space has the advantage of being able to well describe the reduced $\alpha$-width amplitudes of the states of $^{16}$O. Both electric dipole and electric quadrupole transitions are responsible for the radiative-capture process. The $\alpha$ particle is described with a $(0s)^4$ configuration admixed with a small amount of an isospin $T=1$ impurity component, which is crucially important to account for the isovector electric dipole transition. The isoscalar electric dipole operators are also taken into account up to the first order beyond the long-wavelength approximation. All the necessary ingredients are provided to make the paper self-contained and ready for numerical computations.


title: Neural Network Training Techniques Regularize Optimization Trajectory: An Empirical Study

abstract: Modern deep neural network (DNN) trainings utilize various training techniques, e.g., nonlinear activation functions, batch normalization, skip-connections, etc. Despite their effectiveness, it is still mysterious how they help accelerate DNN trainings in practice. In this paper, we provide an empirical study of the regularization effect of these training techniques on DNN optimization. Specifically, we find that the optimization trajectories of successful DNN trainings consistently obey a certain regularity principle that regularizes the model update direction to be aligned with the trajectory direction. Theoretically, we show that such a regularity principle leads to a convergence guarantee in nonconvex optimization and the convergence rate depends on a regularization parameter. Empirically, we find that DNN trainings that apply the training techniques achieve a fast convergence and obey the regularity principle with a large regularization parameter, implying that the model updates are well aligned with the trajectory. On the other hand, DNN trainings without the training techniques have slow convergence and obey the regularity principle with a small regularization parameter, implying that the model updates are not well aligned with the trajectory. Therefore, different training techniques regularize the model update direction via the regularity principle to facilitate the convergence.


title: Quasianalytic functionals and ultradistributions as boundary values of harmonic functions

abstract: We study boundary values of harmonic functions in spaces of quasianalytic functionals and spaces of ultradistributions of non-quasianalytic type. As an application, we provide a new approach to H\"ormander's support theorem for quasianalytic functionals. Our main technical tool is a description of ultradifferentiable functions by almost harmonic functions, a concept that we introduce in this article. We work in the setting of ultradifferentiable classes defined via weight matrices. In particular, our results simultaneously apply to the two standard classes defined via weight sequences and via weight functions.


title: Diffusion models for Handwriting Generation

abstract: In this paper, we propose a diffusion probabilistic model for handwriting generation. Diffusion models are a class of generative models where samples start from Gaussian noise and are gradually denoised to produce output. Our method of handwriting generation does not require using any text-recognition based, writer-style based, or adversarial loss functions, nor does it require training of auxiliary networks. Our model is able to incorporate writer stylistic features directly from image data, eliminating the need for user interaction during sampling. Experiments reveal that our model is able to generate realistic , high quality images of handwritten text in a similar style to a given writer. Our implementation can be found at https://github.com/tcl9876/Diffusion-Handwriting-Generation


title: Towards observing the neutron star collapse with gravitational wave detectors

abstract: Gravitational waves from binary neutron star inspirals have been detected along with the electromagnetic transients coming from the aftermath of the merger in GW170817. However, much is still unknown about the post-merger dynamics that connects these two sets of observables. This includes if, and when, the post-merger remnant star collapses to a black hole, and what are the necessary conditions to power a short gamma-ray burst, and other observed electromagnetic counterparts. Observing the collapse of the post-merger neutron star would shed led on these questions, constraining models for the short gamma-ray burst engine and the hot neutron star equation of state. In this work, we explore the scope of using gravitational wave detectors to measure the timing of the collapse either indirectly, by establishing the shut-off of the post-merger gravitational emission, or---more challengingly---directly, by detecting the collapse signal. For the indirect approach, we consider a kilohertz high-frequency detector design that utilises a previously studied coupled arm cavity and signal recycling cavity resonance. This design would give a signal-to-noise ratio of 0.5\,-\,8.6 (depending on the variation of waveform parameters) for a collapse gravitational wave signal occurring at 10\,ms post-merger of a binary at 50\,Mpc and with total mass $2.7 M_\odot$. For the direct approach, we propose a narrow-band detector design, utilising the sensitivity around the frequency of the arm cavity free spectral range. The proposed detector achieves a signal-to-noise ratio of 0.3\,-\,1.9, independent of the collapse time. This detector is limited by both the fundamental classical and quantum noise with the arm cavity power chosen as 10\,MW.


title: Regression Trees for Cumulative Incidence Functions

abstract: The use of cumulative incidence functions for characterizing the risk of one type of event in the presence of others has become increasingly popular over the past decade. The problems of modeling, estimation and inference have been treated using parametric, nonparametric and semi-parametric methods. Efforts to develop suitable extensions of machine learning methods, such as regression trees and related ensemble methods, have begun only recently. In this paper, we develop a novel approach to building regression trees for estimating cumulative incidence curves in a competing risks setting. The proposed methods employ augmented estimators of the Brier score risk as the primary basis for building and pruning trees. The proposed methods are easily implemented using the R statistical software package. Simulation studies demonstrate the utility of our approach in the competing risks setting. Data from the Radiation Therapy Oncology Group (trial 9410) is used to illustrate these new methods.


title: Network Topology Invariant Stability Certificates for DC Microgrids with Arbitrary Load Dynamics

abstract: DC microgrids are prone to small-signal instabilities due to the presence of tightly regulated loads. This paper develops a decentralized stability certificate which is capable of certifying the small-signal stability of an islanded DC network containing such loads. Utilizing a novel homotopy approach, the proposed standards ensure that no system eigenmodes are able to cross into the unstable right half plane for a continuous range of controller gain levels. The resulting "standards" can be applied to variety of grid components which meet the specified, but non-unique, criteria. These standards thus take a step towards offering plug-and-play operability of DC microgrids. The proposed theorems are explicitly illustrated and numerically validated on multiple DC microgrid test-cases containing both buck and boost converter dynamics.


title: Empirical evidence for a global atmospheric temperature control system: physical structure

abstract: There is evidence that a natural control system influences global atmospheric surface temperature (Leggett and Ball, 2020). The present paper sets up and tests a hypothesis concerning the physical makeup of the sequential elements of the control system and its outcomes. The final outcome the control system influences is defined as global atmospheric surface temperature. The terms used for the control system element types in the hypothesis are, in sequence, leading element, controller and actuator. Actuators are hypothesised to affect in turn the final outcome either directly or via penultimate outcomes. The existence of the control system is evidenced by demonstration of statistically significant one way Granger causality across each step of the hypothesised control system sequence. Evidence is presented that the leading element of the control system, represented by the Normalized Difference Vegetation Index, is the global biosphere. The biosphere as a control system has previously been referred to as Gaia (Lovelock and Margulis, 1974). A fourth (second derivative) term is found to enhance the Proportional, Integral and Derivative process terms of the control system shown in Leggett and Ball (2020). The main actuators of the control system found are shown to be wind speed and cloud cover. Cloud cover is shown to influence the final outcome, global surface temperature, directly. It and wind speed also influence the penultimate outcomes found, those of enhanced ocean heat uptake and enhanced outgoing longwave radiation. These together lead to control system output to the final outcome, global atmospheric temperature. Overall, evidence for the activity of the control system is present across many major physical dimensions of the global atmosphere.


title: Active Reinforcement Learning: Observing Rewards at a Cost

abstract: Active reinforcement learning (ARL) is a variant on reinforcement learning where the agent does not observe the reward unless it chooses to pay a query cost c > 0. The central question of ARL is how to quantify the long-term value of reward information. Even in multi-armed bandits, computing the value of this information is intractable and we have to rely on heuristics. We propose and evaluate several heuristic approaches for ARL in multi-armed bandits and (tabular) Markov decision processes, and discuss and illustrate some challenging aspects of the ARL problem.


title: Comments on Efficient Singular Value Thresholding Computation

abstract: We discuss how to evaluate the proximal operator of a convex and increasing function of a nuclear norm, which forms the key computational step in several first-order optimization algorithms such as (accelerated) proximal gradient descent and ADMM. Various special cases of the problem arise in low-rank matrix completion, dropout training in deep learning and high-order low-rank tensor recovery, although they have all been solved on a case-by-case basis. We provide an unified and efficiently computable procedure for solving this problem.


title: Morris-Thorne Wormhole in the Vector-Tensor theories with Abelian gauge symmetry breaking

abstract: We construct an asymptotically flat Morris-Thorne wormhole solution supported by anisotropic matter fluid and a vector field which is coupled to gravity in a non-minimal way with broken Abelian gauge symmetry. In this paper, a specific shape function is considered. We find that the ansatz of vector field plays a significant role in determining the spacetime geometry of the wormhole. If there exists the electrostatic potential only, the redshift function could be considered as a constant value, implying the vanishing tidal force. However, when the vector potential in radial-direction is involved, the r-component of extended Maxwell equations at the wormhole's throat is invalid. To solve this issue, a thin shell is introduced near the throat, dividing the spacetime into two parts. Furthermore, it is proved that the spacetime geometry of wormhole could be smooth at junction position if the expressions of redshift function and vector potential are given appropriately. Finally, the energy conditions and the volume integral quantifer are explored.


title: Long-Range Surface-Assisted Molecule-Molecule Hybridization

abstract: Metalated phthalocyanines (Pc's) are robust and versatile molecular complexes, whose properties can be tuned by changing their functional groups and central metal atom. The electronic structure of magnesium Pc (MgPc) - structurally and electronically similar to chlorophyll - adsorbed on the Ag(100) surface is investigated by low-temperature scanning tunneling microscopy (STM) and spectroscopy (STS), non-contact atomic force microscopy (ncAFM) and density functional theory (DFT). Single, isolated MgPc's exhibit a flat, four-fold rotationally symmetric morphology, with doubly degenerate, partially populated (due to surface-to-molecule electron transfer) lowest unoccupied molecular orbitals (LUMOs). In contrast, MgPc's with neighbouring molecules in proximity undergo a lift of LUMOs degeneracy, with a near-Fermi local density of states with reduced two-fold rotational symmetry, indicative of a long-range attractive intermolecular interaction. The latter is assigned to a surface-mediated two-step electronic hybridization process. First, LUMOs interact with Ag(100) conduction electrons, forming hybrid molecule-surface orbitals with enhanced spatial extension. Then, these delocalized molecule-surface states further hybridize with those of neighbouring molecules. This work highlights how the electronic structure of molecular adsorbates - including orbital degeneracies and symmetries - can be significantly altered via surface-mediated intermolecular hybridization, over extended distances (beyond 3 nm), having important implications for prospects of molecule-based solid-state technologies.


title: iHorology: Lowering the Barrier to Microsecond-level Internet Time

abstract: High precision, synchronized clocks are essential to a growing number of Internet applications. Standard protocols and their associated server infrastructure have been shown to typically enable client clocks to synchronize on the order of tens of milliseconds. We address one of the key challenges to high precision Internet timekeeping - the intrinsic contribution to clock error of path asymmetry between client and time server, a fundamental barrier to microsecond level accuracy. We first exploit results of a measurement study to quantify asymmetry and its effect on timing. We then describe three approaches to addressing the path asymmetry problem: LBBE, SBBE and K-SBBE, each based on timestamp exchange with multiple servers, with the goal of tightening bounds on asymmetry for each client. We explore their capabilities and limitations through simulation and argument. We show that substantial improvements are possible, and discuss whether, and how, the goal of microsecond accuracy might be attained.


title: Manipulation of Spin Transport in Graphene/Transition Metal Dichalcogenide Heterobilayers upon Twisting

abstract: Proximity effects are one of the pillars of exotic phenomena and technological applications of two dimensional materials. However, the interactions nature depends strongly on the materials involved, their crystalline symmetries, and interfacial properties. Here we used large-scale first-principle calculations to demonstrate that strain and twist-angle are efficient knobs to tailor the spin-orbit coupling in graphene transition metal dichalcogenide heterobilayers. We found that by choosing a twist-angle of 30 degrees, the spin relaxation times increase by two orders of magnitude, opening a path to improve these heterostructures spin transport capability. Moreover, we demonstrate that strain and twist angle will modify the relative values of valley-Zeeman and Rashba spin-orbit coupling, allowing to tune the system into an ideal Dirac-Rashba regime. These results enable us to envision an answer for the variability of spin-orbit coupling found in different experiments and have significant consequences for applications that depend on polycrystallinity, where grains form at different orientations.


title: An Efficient High-Order Meshless Method for Advection-Diffusion Equations on Time-Varying Irregular Domains

abstract: We present a high-order radial basis function finite difference (RBF-FD) framework for the solution of advection-diffusion equations on time-varying domains. Our framework is based on a generalization of the recently developed Overlapped RBF-FD method that utilizes a novel automatic procedure for computing RBF-FD weights on stencils in variable-sized regions around stencil centers. This procedure eliminates the overlap parameter $\delta$, thereby enabling tuning-free assembly of RBF-FD differentiation matrices on moving domains. In addition, our framework utilizes a simple and efficient procedure for updating differentiation matrices on moving domains tiled by node sets of time-varying cardinality. Finally, advection-diffusion in time-varying domains is handled through a combination of rapid node set modification, a new high-order semi-Lagrangian method that utilizes the new tuning-free overlapped RBF-FD method, and a high-order time-integration method. The resulting framework has no tuning parameters and has $O(N \log N)$ time complexity. We demonstrate high-orders of convergence for advection-diffusion equations on time-varying 2D and 3D domains for both small and large Peclet numbers. We also present timings that verify our complexity estimates. Finally, we utilize our method to solve a coupled 3D problem motivated by models of platelet aggregation and coagulation, once again demonstrating high-order convergence rates on a moving domain.


title: Dependency-based Anomaly Detection: a General Framework and Comprehensive Evaluation

abstract: Anomaly detection is crucial for understanding unusual behaviors in data, as anomalies offer valuable insights. This paper introduces Dependency-based Anomaly Detection (DepAD), a general framework that utilizes variable dependencies to uncover meaningful anomalies with better interpretability. DepAD reframes unsupervised anomaly detection as supervised feature selection and prediction tasks, which allows users to tailor anomaly detection algorithms to their specific problems and data. We extensively evaluate representative off-the-shelf techniques for the DepAD framework. Two DepAD algorithms emerge as all-rounders and superior performers in handling a wide range of datasets compared to nine state-of-the-art anomaly detection methods. Additionally, we demonstrate that DepAD algorithms provide new and insightful interpretations for detected anomalies.


title: Coordinated Motion Control and Event-based Obstacle-crossing for Four Wheel-leg Independent Motor-driven Robotic System via MPC

abstract: This work presents the coordinated motion control and obstacle-crossing problem for the four wheel-leg independent motor-driven robotic systems via a model predictive control (MPC) approach based on an event-triggering mechanism. The modeling of a wheel-leg robotic control system with a dynamic supporting polygon is organized. The system dynamic model is 3 degrees of freedom (DOF) ignoring the pitch, roll and vertical motions. The single wheel dynamic is analyzed considering the characteristics of motor-driven and the Burckhardt nonlinear tire model. As a result, an over-actuated predictive model is proposed with the motor torques as inputs and the system states as outputs. As the supporting polygon is only adjusted at certain conditions, an event-based triggering mechanism is designed to save hardware resources and energy. The MPC controller is evaluated on a virtual prototype as well as a physical prototype. The simulation results guide the parameter tuning for the controller implementation in the physical prototype. The experimental results on these two prototypes verify the efficiency of the proposed approach.