title: Improved tracking of particles with highly correlated motion

abstract: Despite significant advances in particle imaging technologies over the past two decades, few advances have been made in particle tracking, i.e. linking individual particle positions across time series data. The state-of-the-art tracking algorithm is highly effective for systems in which the particles behave mostly independently. However, these algorithms become inaccurate when particle motion is highly correlated, such as in dense or strongly interacting systems. Accurate particle tracking is essential in the study of the physics of dense colloids, such as the study of dislocation formation, nucleation, and shear transformations. Here, we present a new method for particle tracking that incorporates information about the correlated motion of the particles. We demonstrate significant improvement over the state-of-the-art tracking algorithm in simulated data on highly correlated systems.


title: Evaluation of data imputation strategies in complex, deeply-phenotyped data sets: the case of the EU-AIMS Longitudinal European Autism Project

abstract: An increasing number of large-scale multi-modal research initiatives has been conducted in the typically developing population, as well as in psychiatric cohorts. Missing data is a common problem in such datasets due to the difficulty of assessing multiple measures on a large number of participants. The consequences of missing data accumulate when researchers aim to explore relationships between multiple measures. Here we aim to evaluate different imputation strategies to fill in missing values in clinical data from a large (total N=764) and deeply characterised (i.e. range of clinical and cognitive instruments administered) sample of N=453 autistic individuals and N=311 control individuals recruited as part of the EU-AIMS Longitudinal European Autism Project (LEAP) consortium. In particular we consider a total of 160 clinical measures divided in 15 overlapping subsets of participants. We use two simple but common univariate strategies, mean and median imputation, as well as a Round Robin regression approach involving four independent multivariate regression models including a linear model, Bayesian Ridge regression, as well as several non-linear models, Decision Trees, Extra Trees and K-Neighbours regression. We evaluate the models using the traditional mean square error towards removed available data, and consider in addition the KL divergence between the observed and the imputed distributions. We show that all of the multivariate approaches tested provide a substantial improvement compared to typical univariate approaches. Further, our analyses reveal that across all 15 data-subsets tested, an Extra Trees regression approach provided the best global results. This allows the selection of a unique model to impute missing data for the LEAP project and deliver a fixed set of imputed clinical data to be used by researchers working with the LEAP dataset in the future.


title: Deep Reinforcement Learning with Spiking Q-learning

abstract: With the help of special neuromorphic hardware, spiking neural networks (SNNs) are expected to realize artificial intelligence (AI) with less energy consumption. It provides a promising energy-efficient way for realistic control tasks by combining SNNs with deep reinforcement learning (RL). There are only a few existing SNN-based RL methods at present. Most of them either lack generalization ability or employ Artificial Neural Networks (ANNs) to estimate value function in training. The former needs to tune numerous hyper-parameters for each scenario, and the latter limits the application of different types of RL algorithm and ignores the large energy consumption in training. To develop a robust spike-based RL method, we draw inspiration from non-spiking interneurons found in insects and propose the deep spiking Q-network (DSQN), using the membrane voltage of non-spiking neurons as the representation of Q-value, which can directly learn robust policies from high-dimensional sensory inputs using end-to-end RL. Experiments conducted on 17 Atari games demonstrate the DSQN is effective and even outperforms the ANN-based deep Q-network (DQN) in most games. Moreover, the experiments show superior learning stability and robustness to adversarial attacks of DSQN.


title: Pneumatic Computers for Embedded Control of Microfluidics

abstract: Alternative computing approaches that interface readily with physical systems are well suited for embedded control of those systems. We demonstrate finite state machines implemented as pneumatic circuits of microfluidic valves, and we employ these controllers to direct microfluidic liquid handling procedures such as serial dilution on the same chip. These monolithic integrated systems require only power to be supplied externally, in the form of a vacuum source. User input can be provided directly to the chip by covering pneumatic ports with a finger. State machines with up to four bits of state memory are demonstrated, and next-state combinational logic can be fully reprogrammed by changing the hole-punch pattern on a membrane in the chip. These pneumatic computers demonstrate a new framework for the embedded control of physical systems.


title: The Price of Symmetric Line Plans in the Parametric City

abstract: We consider the line planning problem in public transport in the Parametric City, an idealized model that captures typical scenarios by a (small) number of parameters. The Parametric City is rotation symmetric, but optimal line plans are not always symmetric. This raises the question to quantify the symmetry gap between the best symmetric and the overall best solution. For our analysis, we formulate the line planning problem as a mixed integer linear program, that can be solved in polynomial time if the solutions are forced to be symmetric. The symmetry gap is provably small when a specific Parametric City parameter is fixed, and we give an approximation algorithm for line planning in the Parametric City in this case. While the symmetry gap can be arbitrarily large in general, we show that symmetric line plans are a good choice in most practical situations.


title: Frames via Unilateral Iterations of a Bounded Operator

abstract: Motivated by recent work in Dynamical Sampling, we prove a necessary and sufficient condition for a frame in a separable and infinite-dimensional Hilbert space to admit the form $\{T^{n} \varphi \}_{n \geq 0}$ with $T \in B(H)$. Also, a characterization of all the vectors $\varphi$ for which $\{T^{n} \varphi \}_{n \geq 0}$ is a frame for some $T \in B(H)$ is provided. Some auxiliary results on operator representations of Riesz frames are given as well.


title: Almost prime ideals in noncommutative rings

abstract: A proper ideal $P$ of a commutative ring with identity is an almost prime ideal if $ab \in P{\setminus}P^2$ implies $a \in P$ or $b \in P$. In this paper we define almost prime ideals of a noncommutative ring, and provide some equivalent definitions. We also examine some cases such that all right ideals of a noncommutative ring are almost prime right ideals.


title: Exploration of Hyperdimensional Computing Strategies for Enhanced Learning on Epileptic Seizure Detection

abstract: Wearable and unobtrusive monitoring and prediction of epileptic seizures has the potential to significantly increase the life quality of patients, but is still an unreached goal due to challenges of real-time detection and wearable devices design. Hyperdimensional (HD) computing has evolved in recent years as a new promising machine learning approach, especially when talking about wearable applications. But in the case of epilepsy detection, standard HD computing is not performing at the level of other state-of-the-art algorithms. This could be due to the inherent complexity of the seizures and their signatures in different biosignals, such as the electroencephalogram (EEG), the highly personalized nature, and the disbalance of seizure and non-seizure instances. In the literature, different strategies for improved learning of HD computing have been proposed, such as iterative (multi-pass) learning, multi-centroid learning and learning with sample weight ("OnlineHD"). Yet, most of them have not been tested on the challenging task of epileptic seizure detection, and it stays unclear whether they can increase the HD computing performance to the level of the current state-of-the-art algorithms, such as random forests. Thus, in this paper, we implement different learning strategies and assess their performance on an individual basis, or in combination, regarding detection performance and memory and computational requirements. Results show that the best-performing algorithm, which is a combination of multi-centroid and multi-pass, can indeed reach the performance of the random forest model on a highly unbalanced dataset imitating a real-life epileptic seizure detection application.


title: Multi-Graph Fusion Networks for Urban Region Embedding

abstract: Learning the embeddings for urban regions from human mobility data can reveal the functionality of regions, and then enables the correlated but distinct tasks such as crime prediction. Human mobility data contains rich but abundant information, which yields to the comprehensive region embeddings for cross domain tasks. In this paper, we propose multi-graph fusion networks (MGFN) to enable the cross domain prediction tasks. First, we integrate the graphs with spatio-temporal similarity as mobility patterns through a mobility graph fusion module. Then, in the mobility pattern joint learning module, we design the multi-level cross-attention mechanism to learn the comprehensive embeddings from multiple mobility patterns based on intra-pattern and inter-pattern messages. Finally, we conduct extensive experiments on real-world urban datasets. Experimental results demonstrate that the proposed MGFN outperforms the state-of-the-art methods by up to 12.35% improvement.


title: Modelling preventive measures and their effect on generation times in emerging epidemics

abstract: We present a stochastic epidemic model to study the effect of various preventive measures, such as uniform reduction of contacts and transmission, vaccination, isolation, screening and contact tracing, on a disease outbreak in a homogeneously mixing community. The model is based on an infectivity process, which we define through stochastic contact and infectiousness processes, so that each individual has an independent infectivity profile. In particular, we monitor variations of the reproduction number and of the distribution of generation times. We show that some interventions, i.e. uniform reduction and vaccination, affect the former while leaving the latter unchanged, whereas other interventions, i.e. isolation, screening and contact tracing, affect both quantities. We provide a theoretical analysis of the variation of these quantities, and we show that, in practice, the variation of the generation time distribution can be significant and that it can cause biases in the estimation of reproduction numbers. The framework, because of its general nature, captures the properties of many infectious diseases, but particular emphasis is on COVID-19, for which numerical results are provided.


title: Symmetry results for compactly supported steady solutions of the 2D Euler equations

abstract: In this paper we prove symmetry of compactly supported steady solutions of the 2D Euler equations. Assuming that $\Omega = \{x \in \mathbb{R}^2:\ u(x) \neq 0\}$ is an annular domain, we prove that the streamlines of the flow are circular. We are also able to remove the topological condition on $\Omega$ if we impose regularity and nondegeneracy assumptions on $u$ at $\partial \Omega$. The proof uses that the corresponding stream function solves an elliptic semilinear problem $-\Delta \phi = f(\phi)$ with $\nabla \phi=0$ at the boundary. One of the main difficulties in our study is that $f$ is not Lipschitz continuous near the boundary values. However, $f(\phi)$ vanishes at the boundary values and then we can apply a local symmetry result of F. Brock to conclude. In the case $\partial_{\nu} u \neq 0$ at $\partial \Omega$ this argument is not possible. In this case we are able to use the moving plane scheme to show symmetry, despite the possible lack of regularity of $f$. We think that such result is interesting in its own right and will be stated and proved also for higher dimensions. The proof requires the study of maximum principles, Hopf lemma and Serrin corner lemma for elliptic linear operators with singular coefficients.


title: A Ti/Pt/Co multilayer stack for transfer function based magnetic force microscopy calibrations

abstract: Magnetic force microscopy is a widespread technique for imaging magnetic structures with a resolution of some 10 nanometers. MFM can be calibrated to obtain quantitative spatially resolved magnetization data in units of A/m by determining the calibrated point spread function of the instrument, its instrument calibration function (ICF), from a measurement of a well-known reference sample. Beyond quantifying the MFM data, a deconvolution of the MFM image data with the ICF also corrects the smearing caused by the finite width of the MFM tip stray field distribution. However, the quality of the calibration depends critically on the calculability of the magnetization distribution of the reference sample. Here, we discuss a Ti/Pt/Co multilayer stack that shows a stripe domain pattern as a suitable reference material. A precise control of the fabrication process, combined with a characterization of the sample micromagnetic parameters, allows reliable calculation of the sample's magnetic stray field, proven by a very good agreement between micromagnetic simulations and qMFM measurements. A calibrated qMFM measurement using the Ti/Pt/Co stack as a reference sample is shown and validated, and the application area for quantitative MFM measurements calibrated with the Ti/Pt/Co stack is discussed


title: QED Meson Description of the X17 and Other Anomalous Particles

abstract: The X17 particle, the E38 particle, and the anomalous soft photons are anomalous particles because they do not appear to belong to any known Standard Model families. We propose a QED meson description of the anomalous particles as composite systems of a light quark and a light antiquark bound and confined by the compact QED interaction, by combining Polyakov's transverse confinement of opposite electric charges in compact QED in (2+1)D and Schwinger's longitudinal confinement for massless opposite electric charges in QED in (1+1)D. With predicted QED meson masses close to the observed X17 and E38 masses, QED mesons may be good candidates for the description of the anomalous particles.


title: Generative Planning for Temporally Coordinated Exploration in Reinforcement Learning

abstract: Standard model-free reinforcement learning algorithms optimize a policy that generates the action to be taken in the current time step in order to maximize expected future return. While flexible, it faces difficulties arising from the inefficient exploration due to its single step nature. In this work, we present Generative Planning method (GPM), which can generate actions not only for the current step, but also for a number of future steps (thus termed as generative planning). This brings several benefits to GPM. Firstly, since GPM is trained by maximizing value, the plans generated from it can be regarded as intentional action sequences for reaching high value regions. GPM can therefore leverage its generated multi-step plans for temporally coordinated exploration towards high value regions, which is potentially more effective than a sequence of actions generated by perturbing each action at single step level, whose consistent movement decays exponentially with the number of exploration steps. Secondly, starting from a crude initial plan generator, GPM can refine it to be adaptive to the task, which, in return, benefits future explorations. This is potentially more effective than commonly used action-repeat strategy, which is non-adaptive in its form of plans. Additionally, since the multi-step plan can be interpreted as the intent of the agent from now to a span of time period into the future, it offers a more informative and intuitive signal for interpretation. Experiments are conducted on several benchmark environments and the results demonstrated its effectiveness compared with several baseline methods.


title: Design Strategies and Approximation Methods for High-Performance Computing Variability Management

abstract: Performance variability management is an active research area in high-performance computing (HPC). We focus on input/output (I/O) variability. To study the performance variability, computer scientists often use grid-based designs (GBDs) to collect I/O variability data, and use mathematical approximation methods to build a prediction model. Mathematical approximation models could be biased particularly if extrapolations are needed. Space-filling designs (SFDs) and surrogate models such as Gaussian process (GP) are popular for data collection and building predictive models. The applicability of SFDs and surrogates in the HPC variability needs investigation. We investigate their applicability in the HPC setting in terms of design efficiency, prediction accuracy, and scalability. We first customize the existing SFDs so that they can be applied in the HPC setting. We conduct a comprehensive investigation of design strategies and the prediction ability of approximation methods. We use both synthetic data simulated from three test functions and the real data from the HPC setting. We then compare different methods in terms of design efficiency, prediction accuracy, and scalability. In synthetic and real data analysis, GP with SFDs outperforms in most scenarios. With respect to approximation models, GP is recommended if the data are collected by SFDs. If data are collected using GBDs, both GP and Delaunay can be considered. With the best choice of approximation method, the performance of SFDs and GBD depends on the property of the underlying surface. For the cases in which SFDs perform better, the number of design points needed for SFDs is about half of or less than that of the GBD to achieve the same prediction accuracy. SFDs that can be tailored to high dimension and non-smooth surface are recommended especially when large numbers of input factors need to be considered in the model.


title: Mid-infrared Variability of Low-redshift Active Galactic Nuclei : Constraints on a Hot Dust Component with a Variable Covering Factor

abstract: We utilize mid-infrared multi-epoch data from Wide-field Infrared Survey Explorer (WISE) over a $\sim10$-year period in W1 (3.4$\mu$m) and W2 (4.6$\mu$m) bands to investigate the structure of dusty torus in low-redshift ($0.15<z\le 0.4$) active galactic nuclei (AGN). We calculate a Spearman correlation coefficient ($r_{12}$) between W1 magnitude and W1$-$W2 color based on the light curve in individual objects. Interestingly, $r_{12}$ spans a broad range from $-1$ to $1$, and is detected to be correlated with mean W1$-$W2 color and AGN bolometric luminosity, in the sense that objects with blue W1$-$W2 color and low AGN luminosity tend to become redder (bluer) with increasing (decreasing) W1 brightness in the light curve (i.e. $r_{12} < 0$), although the correlation of $r_{12}$ with the bolometric luminosity is relatively weak. The fit for the spectral energy distribution reveals a significant contribution from the host galaxy in the W1 and W2 bands. However, the dependencies of $r_{12}$ on W1$-$W2 color and AGN luminosity still persist even after careful elimination of the host light contribution. We propose that this result can be explained if covering factor of the hot dust component decreases as AGN luminosity increases.


title: Study of voltage cycling conditions on Pt oxidation and dissolution in polymer electrolyte fuel cells

abstract: This paper is devoted to study the electrochemical behavior of Pt catalyst in a polymer electrolyte fuel cell at various operating conditions and at different electric potential difference (also known as voltage) cycling applied in accelerated stress tests. The degradation of platinum is considered with respect to the Pt ion dissolution and the Pt oxide coverage of catalyst described by a one-dimensional model. In the model, degradation rate increases with temperature and decreasing particle diameter of Pt nano-particles. The theoretical study of the underlying diffusion system with the nonlinear reactions is presented by analytical methods and gives explicit solutions through a first integral of the ODE system. Numerical tests are obtained using a second order implicit-explicit scheme. The computer simulation shows that the lifetime of the catalyst depends on the voltage profile and the upper potential level. By this Pt mass loss is more significant at the membrane surface than at the gas diffusion layer.


title: A Sorted Datalog Hammer for Supervisor Verification Conditions Modulo Simple Linear Arithmetic

abstract: In a previous paper, we have shown that clause sets belonging to the Horn Bernays-Sch\"onfinkel fragment over simple linear real arithmetic (HBS(SLR)) can be translated into HBS clause sets over a finite set of first-order constants. The translation preserves validity and satisfiability and it is still applicable if we extend our input with positive universally or existentially quantified verification conditions (conjectures). We call this translation a Datalog hammer. The combination of its implementation in SPASS-SPL with the Datalog reasoner VLog establishes an effective way of deciding verification conditions in the Horn fragment. We verify supervisor code for two examples: a lane change assistant in a car and an electronic control unit of a supercharged combustion engine. In this paper, we improve our Datalog hammer in several ways: we generalize it to mixed real-integer arithmetic and finite first-order sorts; we extend the class of acceptable inequalities beyond variable bounds and positively grounded inequalities; and we significantly reduce the size of the hammer output by a soft typing discipline. We call the result the sorted Datalog hammer. It not only allows us to handle more complex supervisor code and to model already considered supervisor code more concisely, but it also improves our performance on real world benchmark examples. Finally, we replace the before file-based interface between SPASS-SPL and VLog by a close coupling resulting in a single executable binary.


title: On a Finite Group Generated by Subnormal Supersoluble Subgroups

abstract: Supersolubility of a finite group $G=\langle A,B\rangle$ with the nilpotent derived subgroup $G^\prime$ is established under the condition that the subgroups $A$ and $B$ are both subnormal and supersoluble.


title: Higher differentiability results in the scale of Besov spaces to a class of double-phase obstacle problems

abstract: We study the higher fractional differentiability properties of the gradient of the solutions to variational obstacle problems of the form \begin{gather*} \min \biggl\{ \int_{\Omega} F(x,w,Dw) d x \ : \ w \in \mathcal{K}_{\psi}(\Omega) \biggr\}, \end{gather*} with $F$ double phase functional of the form \begin{equation*} F(x,w,z)=b(x,w)(|z|^p+a(x)|z|^q), \end{equation*} where $\Omega$ is a bounded open subset of $\mathbb{R}^n$, $\psi \in W^{1,p}(\Omega)$ is a fixed function called \textit{obstacle} and $\mathcal{K}_{\psi}(\Omega)= \{ w \in W^{1,p}(\Omega) : w \geq \psi \ \text{a.e. in} \ \Omega \}$ is the class of admissible functions. Assuming that the gradient of the obstacle belongs to a suitable Besov space, we are able to prove that the gradient of the solution preserves some fractional differentiability property.


title: In Defence of Visual Analytics Systems: Replies to Critics

abstract: The last decade has witnessed many visual analytics (VA) systems that make successful applications to wide-ranging domains like urban analytics and explainable AI. However, their research rigor and contributions have been extensively challenged within the visualization community. We come in defence of VA systems by contributing two interview studies for gathering critics and responses to those criticisms. First, we interview 24 researchers to collect criticisms the review comments on their VA work. Through an iterative coding and refinement process, the interview feedback is summarized into a list of 36 common criticisms. Second, we interview 17 researchers to validate our list and collect their responses, thereby discussing implications for defending and improving the scientific values and rigor of VA systems. We highlight that the presented knowledge is deep, extensive, but also imperfect, provocative, and controversial, and thus recommend reading with an inclusive and critical eye. We hope our work can provide thoughts and foundations for conducting VA research and spark discussions to promote the research field forward more rigorously and vibrantly.


title: Investigations on a Fuzzy Process: Effect of Diffusion on Calibration and Particle Identification in Liquid Argon Time Projection Chambers

abstract: Ionization electron diffusion in Liquid Argon Time Projection Chambers (LArTPCs) has typically been considered at the detector design stage, but little attention has been given to its effects on calibration and particle identification. We use a GEANT4-based simulation to study how diffusion impacts these techniques, and give consideration to how this effect is simulated. We find that diffusion can cause a drift-dependent bias to both the median and Most Probable Value (MPV) of charge deposition per unit length distributions. The bias is estimated to be ~2.5% (median) and ~5.0% (MPV) for typical maximum drift times in currently running LArTPCs before adding detector specific considerations such as electric field non-uniformities. This indicates that these metrics should not be used for calibration without care, contrary to the conventional wisdom. The impact of diffusion on the ability of LArTPCs to separate muons and protons is small, and not expected to pose any problems in future detectors. Diffusion may however be a significant source of systematic uncertainty when separating particles of more similar masses (muons and pions, kaons and protons). Separation of such populations may be improved by implementation of a drift-time dependent particle identification.


title: Hiding Behind Backdoors: Self-Obfuscation Against Generative Models

abstract: Attack vectors that compromise machine learning pipelines in the physical world have been demonstrated in recent research, from perturbations to architectural components. Building on this work, we illustrate the self-obfuscation attack: attackers target a pre-processing model in the system, and poison the training set of generative models to obfuscate a specific class during inference. Our contribution is to describe, implement and evaluate a generalized attack, in the hope of raising awareness regarding the challenge of architectural robustness within the machine learning community.


title: Spin noise in birefringent and inhomogeneous media

abstract: It is known that linear birefringence of the medium essentially hinders measuring the Faraday effect. For this reason, optically anisotropic materials have never been considered as objects of the Faraday-rotation-based spin noise spectroscopy (SNS). We show, both theoretically and experimentally, that strong optical anisotropy that may badly suppress the regular Faraday rotation of the medium, practically does not affect the measurement of the spatially uncorrelated spin fluctuations. An important consequence of this result is that the Faraday-rotation noise should be also insensitive to spatially nonuniform birefringence, which makes the SNS applicable to a wide class of optically anisotropic and inhomogeneous materials. We also show that the birefringent media provide additional opportunity to measure spatial spin correlations. Results of the experimental measurements of the spin-noise spectra performed on 4$f$--4$f$ transitions of Nd$^{3+}$ ions in the CaWO$_4$ and LiYF$_4$ crystals well agree with the theory.


title: Plasma Cleaning of LCLS-II-HE verification cryomodule cavities

abstract: Plasma cleaning is a technique that can be applied in superconducting radio-frequency (SRF) cavities in situ in cryomodules in order to decrease their level of field emission. We developed the technique for the Linac Coherent Light Source II (LCLS-II) cavities and we present in this paper the full development and application of plasma processing to the LCLS-II High Energy (HE) verification cryomodule (vCM). We validated our plasma processing procedure on the vCM, fully processing four out of eight cavities of this CM, demonstrating that cavities performance were preserved in terms of both accelerating field and quality factor. Applying plasma processing to this clean, record breaking cryomodule also showed that no contaminants were introduced in the string, maintaining the vCM field emission-free up to the maximum field reached by each cavity. We also found that plasma processing eliminates multipacting (MP) induced quenches that are typically observed frequently within the MP band field range. This suggests that plasma processing could be employed in situ in CMs to mitigate both field emission and multipacting, significantly decreasing the testing time of cryomodules, the linac commissioning time and cost and increasing the accelerator reliability.


title: RISING a new framework for few-view tomographic image reconstruction with deep learning

abstract: This paper proposes a new two-step procedure for sparse-view tomographic image reconstruction. It is called RISING, since it combines an early-stopped Rapid Iterative Solver with a subsequent Iteration Network-based Gaining step. So far, regularized iterative methods have widely been used for X-ray computed tomography image reconstruction from low-sampled data, since they converge to a sparse solution in a suitable domain, as upheld by the Compressed Sensing theory. Unfortunately, their use is practically limited by their high computational cost which imposes to perform only a few iterations in the available time for clinical exams. Data-driven methods, using neural networks to post-process a coarse and noisy image obtained from geometrical algorithms, have been recently studied and appreciated for both their computational speed and accurate reconstructions. However, there is no evidence, neither theoretically nor numerically, that neural networks based algorithms solve the mathematical inverse problem modeling the tomographic reconstruction process. In our two-step approach, the first phase executes very few iterations of a regularized model-based algorithm whereas the second step completes the missing iterations by means of a neural network. The resulting hybrid deep-variational framework preserves the convergence properties of the iterative method and, at the same time, it exploits the computational speed and flexibility of a data-driven approach. Experiments performed on a simulated and a real data set confirm the numerical and visual accuracy of the reconstructed RISING images in short computational times.


title: Traceable Dosimetry for MeV Ion Beams

abstract: Standard dosimetry protocols exist for highly penetrating photon and particle beams used in the clinic and in research. However, these protocols cannot be directly applied to shallow penetration MeV-range ion beams. The Radiological Research Accelerator Facility has been using such beams for almost 50 years to irradiate cell monolayers, using self-developed dosimetry, based on tissue equivalent ionization chambers. To better align with the internationally accepted standards, we describe implementation of a commercial, NIST-traceable, air-filled ionization chamber for measurement of absorbed dose to water from low energy ions, using radiation quality correction factors calculated using TRS-398 recommendations. The reported dose does not depend on the ionization density in the range of 10-150 keV/micron.


title: Operating in a deep underground facility improves the locking of gradiometric fluxonium qubits at the sweet spots

abstract: We demonstrate flux-bias locking and operation of a gradiometric fluxonium artificial atom using two symmetric granular aluminum (grAl) loops to implement the superinductor. The gradiometric fluxonium shows two orders of magnitude suppression of sensitivity to homogeneous magnetic fields, which can be an asset for hybrid quantum systems requiring strong magnetic field biasing. By cooling down the device in an external magnetic field while crossing the metal-to-superconductor transition, the gradiometric fluxonium can be locked either at $0$ or $\Phi_0/2$ effective flux bias, corresponding to an even or odd number of trapped fluxons, respectively. At mK temperatures, the fluxon parity prepared during initialization survives to magnetic field bias exceeding $100 \,\Phi_0$. However, even for states biased in the vicinity of $1 \,\Phi_0$, we observe unexpectedly short fluxon lifetimes of a few hours, which cannot be explained by thermal or quantum phase slips. When operating in a deep-underground cryostat of the Gran Sasso laboratory, the fluxon lifetimes increase to days, indicating that ionizing events activate phase slips in the grAl superinductor.


title: Well-posedness, ill-posedness, and traveling waves for models of pulsatile flow in viscoelastic vessels

abstract: We study dispersive models of fluid flow in viscoelastic vessels, derived in the study of blood flow. The unknowns in the models are the velocity of the fluid in the axial direction and the displacement of the vessel wall from rest. We prove that one such model has a well-posed initial value problem, while we argue that a related model instead has an ill-posed initial value problem; in the second case, we still prove the existence of solutions in analytic function spaces. Finally we prove the existence of some periodic traveling waves.


title: Uncertainty principles with error term in Gelfand-Shilov spaces

abstract: In this note, an alternative approach to establish observability for semigroups based on their smoothing properties is presented. The results discussed here are closely related to those recently obtained in [arXiv:2112.01788], but the current proof allows to get rid of several technical assumptions by following the standard complex analytic approach established by Kovrijkine combined with an idea from [arXiv:2201.02370].


title: Inferring taxonomic placement from DNA barcoding allowing discovery of new taxa

abstract: In ecology it has become common to apply DNA barcoding to biological samples leading to datasets containing a large number of nucleotide sequences. The focus is then on inferring the taxonomic placement of each of these sequences by leveraging on existing databases containing reference sequences having known taxa. This is highly challenging because i) sequencing is typically only available for a relatively small region of the genome due to cost considerations; ii) many of the sequences are from organisms that are either unknown to science or for which there are no reference sequences available. These issues can lead to substantial classification uncertainty, particularly in inferring new taxa. To address these challenges, we propose a new class of Bayesian nonparametric taxonomic classifiers, BayesANT, which use species sampling model priors to allow new taxa to be discovered at each taxonomic rank. Using a simple product multinomial likelihood with conjugate Dirichlet priors at the lowest rank, a highly efficient algorithm is developed to provide a probabilistic prediction of the taxa placement of each sequence at each rank. BayesANT is shown to have excellent performance in real data, including when many sequences in the test set belong to taxa unobserved in training.


title: Multi-multifractality and dynamic scaling in stochastic porous lattice

abstract: In this article, we extend the idea of stochastic dyadic Cantor set to weighted planar stochastic lattice that leads to a stochastic porous lattice. The process starts with an initiator which we choose to be a square of unit area for convenience. We then define a generator that divides the initiator or one of the blocks, picked preferentially with respect to their areas, to divide it either horizontally or vertically into two rectangles of which one of them is removed with probability $q=1-p$. We find that the remaining number of blocks and their mass varies with time as $t^{p}$ and $t^{-q}$ respectively. Analytical solution shows that the dynamics of this process is governed by infinitely many hidden conserved quantities each of which is a multifractal measure with porous structure as it contains missing blocks of various different sizes. The support where these measures are distributed is fractal with fractal dimension $2p$ provided $0<p<1$. We find that if the remaining blocks are characterized by their respective area then the corresponding block size distribution function obeys dynamic scaling.


title: Time distance based computation of the state space of preemptive real time systems

abstract: We explore in this paper a novel approach that builds an overapproximation of the state space of preemptive real time systems. Our graph construction extends the expression of a class to the time distance system that encodes the quantitative properties of past fired subsequences. This makes it possible to restore relevant time information that is used to tighten still more the DBM overapproximation of reachable classes. We succeed thereby to build efficiently tighter approximated graphs which are more appropriate to restore the quantitative properties of the model. The simulation results show that the computed graphs are of the same size as the exact graphs while improving by far the times needed for their computation.


title: Unifying and Boosting Gradient-Based Training-Free Neural Architecture Search

abstract: Neural architecture search (NAS) has gained immense popularity owing to its ability to automate neural architecture design. A number of training-free metrics are recently proposed to realize NAS without training, hence making NAS more scalable. Despite their competitive empirical performances, a unified theoretical understanding of these training-free metrics is lacking. As a consequence, (a) the relationships among these metrics are unclear, (b) there is no theoretical interpretation for their empirical performances, and (c) there may exist untapped potential in existing training-free NAS, which probably can be unveiled through a unified theoretical understanding. To this end, this paper presents a unified theoretical analysis of gradient-based training-free NAS, which allows us to (a) theoretically study their relationships, (b) theoretically guarantee their generalization performances, and (c) exploit our unified theoretical understanding to develop a novel framework named hybrid NAS (HNAS) which consistently boosts training-free NAS in a principled way. Remarkably, HNAS can enjoy the advantages of both training-free (i.e., the superior search efficiency) and training-based (i.e., the remarkable search effectiveness) NAS, which we have demonstrated through extensive experiments.


title: Aerial Energy Provisioning for Massive Energy-Constrained IoT by UAVs

abstract: Autonomy of devices is a major challenge in many Internet of Things (IoT) applications, in particular when the nodes are deployed remotely or difficult to assess places. In this paper we present an approach to provide energy to these devices by Unmanned Aerial Vehicles (UAVs). Therefore, the two major challenges, finding and charging the node are presented. We propose a model to give the energy constrained node an unlimited autonomy by taken the Wireless Power Transfer (WPT) link and battery capacity into account. Selecting the most suitable battery technology allows a reduction in battery capacity and waste. Moreover, an upgrade of existing IoT nodes is feasible with a limited impact on the design and form factor.


title: Using Computational Grounded Theory to Understand Tutors' Experiences in the Gig Economy

abstract: The introduction of online marketplace platforms has led to the advent of new forms of flexible, on-demand (or 'gig') work. Yet, most prior research concerning the experience of gig workers examines delivery or crowdsourcing platforms, while the experience of the large numbers of workers who undertake educational labour in the form of tutoring gigs remains understudied. To address this, we use a computational grounded theory approach to analyse tutors' discussions on Reddit. This approach consists of three phases including data exploration, modelling and human-centred interpretation. We use both validation and human evaluation to increase the trustworthiness and reliability of the computational methods. This paper is a work in progress and reports on the first of the three phases of this approach.


title: Tidal stripping of dark matter subhalos by baryons from analytical perspectives: disk shocking and encounters with stars

abstract: The cold dark matter (CDM) scenario predicts that galactic halos should host a huge amount of subhalos possibly lighter than planets, depending on the nature of dark matter. Predicting their abundance and distribution has important implications for dark matter searches and searches for subhalos themselves, as they could provide a decisive test of the CDM paradigm. A major difficulty in subhalo population model building is to account for the gravitational stripping induced by baryons, which strongly impact on the overall dynamics inside galaxies. In this paper, we focus on these "baryonic" tides from analytical perspectives, summarizing previous work on galactic disk shocking, and thoroughly revisiting the impact of individual encounters with stars. For the latter, we go beyond the reference calculation of Gerhard and Fall (1983) to deal with penetrative encounters, and provide new analytical results. Based upon a full statistical analysis of subhalo energy change during multiple stellar encounters possibly occurring during disk crossing, we show that subhalos lighter than $\sim 1$~M$_\odot$ are very efficiently pruned by stellar encounters. This modifies their mass function in a stellar environment. In contrast, disk shocking is more efficient at pruning massive subhalos. In short, if reasonably resilient, subhalos surviving disk crossing have lost all their mass but an inner cuspy part, with a tidal mass function strongly departing from the cosmological one. If fragile, stellar encounters make their number density drop by an additional order of magnitude with respect to disk-shocking effects only (e.g., at the solar position in the Milky Way). Our results can be incorporated to any analytical or numerical subhalo population model, as we show for illustration. This study complements those based on cosmological simulations, which cannot resolve dark matter subhalos on such small scales.


title: Theoretical and Observational Aspecs in Metric-Affine Gravity: A field theoretic perspective

abstract: In this PhD thesis we deal with several theoretical and phenomenological apsects of metric-affine theories of gravity. Concretely, we first give a broad introduction to the necessary tools to understand the framework and elaborate on some subtleties of the minimal coupling prescription between geometry and matter in presence of torsion and nonmetricity. Then we dedicate the central part of the thesis to study the structure of Ricci Based gravity (RBG) theories, which will be of later use to understand generic properties of metric-affine theories. We begin by analysing the structure of the RBG field equations and nontrivial aspects of their solution space. We then analyse the abrosption spectra of some spherically symmetric solutions. Then, we show that, if the projective symmetry in these theories is explicitly broken, then there arise ghost degrees of freedom, and we argue that this will be a generic feature of metric-affine gravity theories. Having done this, we analyse metricafine theories through the EFT lens, showing how the nonmetricity tkes a particular form in generic theories where the symmetrised Ricci tensor appears in the action beyond the Einstein-Hilbert term. This sources effective interactions that we use to place tight constraints to these theories. In the third part of the thesis we present a miscelanea of works which are not so related to the structure of RBG theories. We begin by studying a model for spontaneous breaking of Lorentz symmetry, namely the bumblebee model, in the metric-affine approach. In the following chapter we generalise a conformal invariant definition of proper time given by Perlick to the case with general nonmetricity. Finally, we present arguments that show that the recently proposed D4EGB theory is not well defined in its original form. We finish with a brief outlook.


title: Linear Laws of Markov Chains with an Application for Anomaly Detection in Bitcoin Prices

abstract: The goals of this paper are twofold: (1) to present a new method that is able to find linear laws governing the time evolution of Markov chains and (2) to apply this method for anomaly detection in Bitcoin prices. To accomplish these goals, first, the linear laws of Markov chains are derived by using the time embedding of their (categorical) autocorrelation function. Then, a binary series is generated from the first difference of Bitcoin exchange rate (against the United States Dollar). Finally, the minimum number of parameters describing the linear laws of this series is identified through stepped time windows. Based on the results, linear laws typically became more complex (containing an additional third parameter that indicates hidden Markov property) in two periods: before the crash of cryptocurrency markets inducted by the COVID-19 pandemic (12 March 2020), and before the record-breaking surge in the price of Bitcoin (Q4 2020 - Q1 2021). In addition, the locally high values of this third parameter are often related to short-term price peaks, which suggests price manipulation.


title: Nonconventional Quantized Hall Resistances Obtained with $\nu = 2$ Equilibration in Epitaxial Graphene $p-n$ Junctions

abstract: We have demonstrated the millimeter-scale fabrication of monolayer epitaxial graphene $p-n$ junction devices using simple ultraviolet photolithography, thereby significantly reducing device processing time compared to that of electron beam lithography typically used for obtaining sharp junctions. This work presents measurements yielding nonconventional, fractional multiples of the typical quantized Hall resistance at $\nu=2$ ($R_H\approx 12906 {\Omega}$) that take the form: $\frac{a}{b}R_H$. Here, $a$ and $b$ have been observed to take on values such 1, 2, 3, and 5 to form various coefficients of $R_H$. Additionally, we provide a framework for exploring future device configurations using the LTspice circuit simulator as a guide to understand the abundance of available fractions one may be able to measure. These results support the potential for drastically simplifying device processing time and may be used for many other two-dimensional materials.


title: Patches Are All You Need?

abstract: Although convolutional networks have been the dominant architecture for vision tasks for many years, recent experiments have shown that Transformer-based models, most notably the Vision Transformer (ViT), may exceed their performance in some settings. However, due to the quadratic runtime of the self-attention layers in Transformers, ViTs require the use of patch embeddings, which group together small regions of the image into single input features, in order to be applied to larger image sizes. This raises a question: Is the performance of ViTs due to the inherently-more-powerful Transformer architecture, or is it at least partly due to using patches as the input representation? In this paper, we present some evidence for the latter: specifically, we propose the ConvMixer, an extremely simple model that is similar in spirit to the ViT and the even-more-basic MLP-Mixer in that it operates directly on patches as input, separates the mixing of spatial and channel dimensions, and maintains equal size and resolution throughout the network. In contrast, however, the ConvMixer uses only standard convolutions to achieve the mixing steps. Despite its simplicity, we show that the ConvMixer outperforms the ViT, MLP-Mixer, and some of their variants for similar parameter counts and data set sizes, in addition to outperforming classical vision models such as the ResNet. Our code is available at https://github.com/locuslab/convmixer.


title: Interpolation of curves on Fano hypersurfaces

abstract: On a general hypersurface of degree $d\leq n$ in $\mathbb P^n$ or $\mathbb P^n$ itself, we prove the existence of curves of any genus and high enough degree depending on the genus passing through the expected number $t$ of general points or incident to a general collection of subvarieties of suitable codimensions. In some cases we also show that the family of curves through $t$ fixed points has general moduli as family of $t$-pointed curves. These results imply positivity of certain intersection numbers on Kontsevich spaces of stable maps. An arithmetical appendix by M. C. Chang descibes the set of numerical characters ($n, d$, curve degree, genus) to which our results apply.


title: A note on increasing paths in countable hypergraphs

abstract: An old result of M\"uller and R\"odl states that a countable graph $G$ has a subgraph whose vertices all have infinite degree if and only if for any vertex labeling of $G$ by positive integers, an infinite increasing path can be found. They asked whether an analogous equivalence holds for edge labelings, which Reiterman answered in the affirmative. Recently, Arman, Elliott, and R\"odl extended this problem to linear $k$-uniform hypergraphs $H$ and generalized the original equivalence for vertex labelings. They asked whether Reiterman's result for edge labelings can similarly be extended. We confirm this for the case where $H$ admits only finitely many Berge cycles.


title: Driven toroidal helix as a generalization of Kapitzas pendulum

abstract: We explore a model system consisting of a particle confined to move along a toroidal helix while being exposed to a static potential as well as a driving force due to a harmonically oscillating electric field. It is shown that in the limit of a vanishing helix radius the governing equations of motion coincide with those of the well-known Kapitza pendulum - a classical pendulum with oscillating pivot - implying that the driven toroidal helix represents a corresponding generalization. It is shown that the two dominant static fixed points present in the Kapitza pendulum are also present for a finite helix radius. The dependence of the stability of these two fixed points on the helix radius, the driving amplitude, and the static potential are analyzed both analytically and numerically. Additionally, the most prominent deviations of the driven helix from Kapitzas pendulum with respect to the resulting phase space are investigated and analyzed in some detail. These effects include an unusual transition to chaos and an effective directed transport due to the simultaneous presence of multiple chaotic phase space regions.


title: Inhomogeneous Jets from Neutron Star Mergers: One Jet to Rule them all

abstract: Using the resultant profiles from 3D hydrodynamic simulations of relativistic jets interacting with neutron star merger wind ejecta, we show how the inhomogeneity of energy and velocity {across the jet surface profile} can alter the observed afterglow lightcurve. We find that the peak afterglow flux depends sensitively on the observer's line-of-sight, not only via the jet inclination but also through the jet rotation: for an observer viewing the afterglow within the GRB-bright jet core, we find a peak flux variability on the order $<0.5$ dex through rotational orientation and $<1.3$ dex for the polar inclination. An observed afterglow's peak flux can be used to infer the jet kinetic energy, and where a top-hat jet is assumed, we find the range of inferred jet kinetic energies for our various model afterglow lightcurves (with fixed model parameters), covers $\sim 1/3$ of the observed short GRB population. Additionally, we present an analytic jet structure function that includes physically motivated parameter uncertainties due to variability through the rotation of the source. % An approximation for the change in collimation due to the merger ejecta mass is included and we show that by considering the observed range of merger ejecta masses from short GRB kilonova candidates, a population of merger jets with a fixed intrinsic jet energy is capable of explaining the observed broad diversity seen in short GRB afterglows.


title: ExoCAM: A 3D Climate Model for Exoplanet Atmospheres

abstract: The TRAPPIST-1 Habitable Atmosphere Intercomparison (THAI) project was initiated to compare 3D climate models that are commonly used for predicting theoretical climates of habitable zone extrasolar planets. One of the core models studied as part of THAI is ExoCAM, an independently curated exoplanet branch of the National Center for Atmospheric Research (NCAR) Community Earth System Model (CESM) version 1.2.1. ExoCAM has been used for studying atmospheres of terrestrial extrasolar planets around a variety of stars. To accompany the THAI project and provide a primary reference, here we describe ExoCAM and what makes it unique from standard configurations of CESM. Furthermore, we also conduct a series of intramodel sensitivity tests of relevant moist physical tuning parameters while using the THAI protocol as our starting point. A common criticism of 3D climate models used for exoplanet modeling is that cloud and convection routines often contain free parameters that are tuned to the modern Earth, and thus may be a source of uncertainty in evaluating exoplanet climates. Here, we explore sensitivities to numerous configuration and parameter selections, including a recently updated radiation scheme, a different cloud and convection physics package, different cloud and precipitation tuning parameters, and a different sea ice albedo. Improvements to our radiation scheme and the modification of cloud particle sizes have the largest effect on global mean temperatures, with variations up to ~10 K, highlighting the requirement for accurate radiative transfer and the importance of cloud microphysics for simulating exoplanetary climates. However for the vast majority of sensitivity tests, climate differences are small. For all cases studied, intramodel differences do not bias general conclusions regarding climate states and habitability.


title: IMO$^3$: Interactive Multi-Objective Off-Policy Optimization

abstract: Most real-world optimization problems have multiple objectives. A system designer needs to find a policy that trades off these objectives to reach a desired operating point. This problem has been studied extensively in the setting of known objective functions. We consider a more practical but challenging setting of unknown objective functions. In industry, this problem is mostly approached with online A/B testing, which is often costly and inefficient. As an alternative, we propose interactive multi-objective off-policy optimization (IMO$^3$). The key idea in our approach is to interact with a system designer using policies evaluated in an off-policy fashion to uncover which policy maximizes her unknown utility function. We theoretically show that IMO$^3$ identifies a near-optimal policy with high probability, depending on the amount of feedback from the designer and training data for off-policy estimation. We demonstrate its effectiveness empirically on multiple multi-objective optimization problems.


title: Neural Architecture Searching for Facial Attributes-based Depression Recognition

abstract: Recent studies show that depression can be partially reflected from human facial attributes. Since facial attributes have various data structure and carry different information, existing approaches fail to specifically consider the optimal way to extract depression-related features from each of them, as well as investigates the best fusion strategy. In this paper, we propose to extend Neural Architecture Search (NAS) technique for designing an optimal model for multiple facial attributes-based depression recognition, which can be efficiently and robustly implemented in a small dataset. Our approach first conducts a warmer up step to the feature extractor of each facial attribute, aiming to largely reduce the search space and providing customized architecture, where each feature extractor can be either a Convolution Neural Networks (CNN) or Graph Neural Networks (GNN). Then, we conduct an end-to-end architecture search for all feature extractors and the fusion network, allowing the complementary depression cues to be optimally combined with less redundancy. The experimental results on AVEC 2016 dataset show that the model explored by our approach achieves breakthrough performance with 27\% and 30\% RMSE and MAE improvements over the existing state-of-the-art. In light of these findings, this paper provides solid evidences and a strong baseline for applying NAS to time-series data-based mental health analysis.


title: Low temperature specific heat and thermal conductivity in superconducting UTe2

abstract: The measurements (Phys.Rev.B {\bf 100}, 220504(R) (2019)) do not detect noticeable thermal conductivity in superconducting UTe$_2$ in the T=0 limit. At the same time the same crystals exhibit a large residual density of states comparable with its normal state value. The improvement of samples quality leads to the augmentation of critical temperature of transition to the superconducting state accompanied by the decreasing of residual specific heat ratio $C(T)/T$ at $T\to 0$. There is presented an analytic derivation of this inverse correlation property in concrete cases of triplet superconducting states with node-less and point-nodes order parameters allowed by symmetry in UTe$_2$. The obtained explicit formulas for the residual thermal conductivity are in reasonable correspondence with observations.


title: How to Solve "The Hardest Logic Puzzle Ever" and Its Generalization

abstract: Raymond Smullyan came up with a puzzle that George Boolos called "The Hardest Logic Puzzle Ever".[1] The puzzle has truthful, lying, and random gods who answer yes or no questions with words that we don't know the meaning of. The challenge is to figure out which type each god is. The puzzle has attracted some general attention -- for example, one popular presentation of the puzzle has been viewed 10 million times.[2] Various "top-down" solutions to the puzzle have been developed.[1,3] Here a systematic bottom-up approach to the puzzle and its generalization is presented. We prove that an n gods puzzle is solvable if and only if the random gods are less than the non-random gods, for arbitrary cardinals. We develop a solution using 4.15 questions to the 5 gods variant with 2 random and 3 lying gods.


title: Constrained Policy Optimization via Bayesian World Models

abstract: Improving sample-efficiency and safety are crucial challenges when deploying reinforcement learning in high-stakes real world applications. We propose LAMBDA, a novel model-based approach for policy optimization in safety critical tasks modeled via constrained Markov decision processes. Our approach utilizes Bayesian world models, and harnesses the resulting uncertainty to maximize optimistic upper bounds on the task objective, as well as pessimistic upper bounds on the safety constraints. We demonstrate LAMBDA's state of the art performance on the Safety-Gym benchmark suite in terms of sample efficiency and constraint violation.


title: Deciphering dynamics of recent COVID-19 outbreak in India: An age-structured modeling

abstract: Infectious disease transmission dynamics are particularly sensitive to social contact patterns, and the precautions people take to limit disease transmission. It depends on the age distribution of the community. Thus, knowing the age$-$specific prevalence and incidence of infectious diseases is critical for predicting future disease burden and the efficacy of interventions like immunization. This study uses an SEIR age$-$structured multi-group epidemic model to understand how social contact affects disease control. We created location$-$specific social contact matrices in the community to see how social mixing has affected illness spread. We estimated the basic reproduction number $(R_0)$ of the system and plotted its global behavior in terms of $R_0$. Optimal control for the problem has also been established quantitatively. The proposed model's transmission rate for India from September 1, 2020, to December 31, 2020, has also been estimated. We replicated the lifting of non-pharmaceutical therapies by allowing participants to return to work in phases and studied the impact of this. Our findings imply that identifying symptomatic sick people aged $20-49$ can help lower the number of infected people when schools are closed. When some schools are partially open, awareness of symptomatic infected persons helps reduce cases. The simulation results also suggest that limiting contact at school and other meeting areas could significantly lower the number of instances. Using the least square approach, it was discovered that the time$-$dependent transmission rate is more realistic than the constant spread rate for COVID$-19$ in India. To reduce the COVID$-19$ in burden, we found that gradually loosening control methods could flatten and lower other peaks. Our findings may help health policymakers decide on timely age$-$based immunization distribution strategies and hence control the disease.


title: Estimates of the change rate of solar mass and gravitational constant based on the dynamics of the Solar System

abstract: The estimate of the change rate of the solar gravitational parameter $\mathrm{d}(GM_{\odot})/\mathrm{d}t$ is obtained from processing modern positional observations of planets and spacecraft. Observations were processed and parameters were determined basing on the numerical planetary ephemeris EPM2019. The obtained annual decrease in solar mass $M_{\odot}$ accounts for the loss through radiation ${\dot M}_{{\odot}\mathrm{rad}}$, through the outgoing solar wind ${\dot M}_{{\odot}{\mathrm{wind}}}$, and for the material falling on the Sun ${\dot M}_{{\odot}\mathrm{fall}}$. The estimated relative value is within $-13.4\cdot 10^{-14} < (\dot M_{\odot} /M_{\odot})_{\mathrm{rad}+\mathrm{wind}+ \mathrm{fall}} < -8.7\cdot 10^{-14}$ per year. The following range for the change rate of the gravitational constant $G$ was obtained: $-2.9 \cdot 10^{-14} < \dot G /G < + 4.6 \cdot 10^{-14}$ per year $(3\sigma$). The new result reduces the interval for the change in $G$ and narrows the limits of possible deviations for alternative gravitational theories from general relativity.

categories: astro-ph.EP astro-ph.IM astro-ph.SR gr-qc


title: Non-Linear Dynamic Inversion with Actuator Dynamics: an Incremental Control Perspective

abstract: In this paper, we derive a sensor based Nonlinear Dynamic Inversion (NDI) control law for a nonlinear system with first-order linear actuators, and compare it to Incremental Nonlinear Dynamic Inversion (INDI), which has gained popularity in recent years. It is shown that for first-order actuator dynamics, INDI approximates the corresponding NDI control law arbitrarily well under the condition of sufficiently fast actuators. If the actuator bandwidth is low compared to changes in the states, the derived NDI control law has the following advantages compared to INDI: 1) compensation of state derivative terms 2) well defined error dynamics and 3) exact tracking of a reference model, independent of error controller gains in nominal conditions. The comparison of the INDI control law with the well established control design method NDI adds to the understanding of incremental control. It is additionally shown how to quantify the deficiency of the INDI control law with respect to the exact NDI law for actuators with finite bandwidth. The results are confirmed through simulation results of the rolling motion of a fixed wing aircraft.


title: Infinite Growth: A Curse or a Blessing?

abstract: The article discusses the controversy over infinite growth. It analyzes the two dominant perspectives on economic growth and finds them based on a limited and subjective view of reality. An examination of the principal aspects of economic growth (economic activity, value, and value creation) helps to understand what fuels economic growth. The article also discusses the correlations between production, consumption, as well as resources and population dynamics. The article finds that infinite and exponential economic growth is essential for the survival of our civilization.


title: Transformation-optics modeling of 3D-printed freeform waveguides

abstract: Multi-photon lithography allows to complement planar photonic integrated circuits (PIC) by in-situ 3D-printed freeform waveguide structures. However, design and optimization of such freeform waveguides using time-domain Maxwell's equations solvers often requires comparatively large computational volumes, within which the structure of interest only occupies a small fraction, thus leading to poor computational efficiency. In this paper, we present a solver-independent transformation-optics-(TO-) based technique that allows to greatly reduce the computational effort related to modeling of 3D freeform waveguides. The concept relies on transforming freeform waveguides with curved trajectories into equivalent waveguide structures with modified material properties but geometrically straight trajectories, that can be efficiently fit into rather small cuboid-shaped computational volumes. We demonstrate the viability of the technique and benchmark its performance using a series of different freeform waveguides, achieving a reduction of the simulation time by a factor of 3-6 with a significant potential for further improvement. We also fabricate and experimentally test the simulated waveguides by 3D-printing on a silicon photonic chip, and we find good agreement between the simulated and the measured transmission at $\lambda = 1550 \textrm{ nm}$.


title: Hund's coupling and electronic anisotropy in the spin-density wave state of iron pnictides

abstract: In the multiband systems, Hund's coupling ($J$) plays a significant role in the spin and charge excitations. We study the dependence of electronic anisotropy on $J$ in terms of Drude-weight along different directions as well as the orbital order in the four-fold symmetry broken spin-density wave state of iron pnictides. A robust behavior of the Drude-weight anisotropy within a small window around $J \sim 0.25U$ with $U$ as intraorbital Coulomb interaction is described in terms of orbital-weight distribution along the reconstructed Fermi surfaces. We also find that the ferro-orbital order increases with $J$ in the widely accepted regime for the latter, which is explained as a consequence of rising exchange field with an increase in magnetization.


title: Reconstruction for the time-dependent coefficients of a quasilinear dynamical Schr{\"o}dinger equation

abstract: We study an inverse problem related to the dynamical Schr{\"o}dinger equation in a bounded domain of $\Rb^n,n\geq 2$. Since the concerned non-linear Schr\"odinger equation possesses a trivial solution, we linearize the equation around the trivial solution. Demonstrating the well-posedness of the direct problem under appropriate conditions on initial and boundary data, it is observed that the solution admits $\eps$-expansion. By taking into account the fact that the terms $\Oh(|\nabla u(t,x)|^3)$ are negligible in this context, we shall reconstruct the time-dependent coefficients such as electric potential and vector-valued function associated with quadratic nonlinearity from the knowledge of input-output map using the geometric optics solution and Fourier inversion.


title: Dispersion of activity at an active-passive nematic interface

abstract: Efficient nutrient mixing is crucial for the survival of bacterial colonies and other living systems. This raises the question of whether the optimization of mixing through the emergence of active turbulent motion in bacterial swarms played a role in the evolution of bacterial shapes. Here, to address this question, we solve the hydrodynamic equation for active nematics coupled with an advection-diffusion equation for the nutrients. The latter models a conserved activity field and mimics the conservation of nutrients in bacterial swarms. At the interface between active and passive nematic phases, in addition to diffusion, the activity is transported by interfacial flows and in turn modifies them through active stresses. We find that the interfacial dispersion of the conserved activity is subdiffusive due to the emergence of a barrier of negative defects at the active-passive interface, which hinders the propagation of the motile positive defects. Furthermore, we observe a non-monotonic dependence of the generalized diffusion coefficient on the aligning parameter, which is related to the shape of the particles. Our simulations suggest that there is an optimal shape that maximizes the dispersion of conserved activity at the active-passive nematic interface.


title: Imputing Missing Values in the Occupational Requirements Survey

abstract: The U.S. Bureau of Labor Statistics allows public access to much of the data acquired through its Occupational Requirements Survey (ORS). This data can be used to draw inferences about the requirements of various jobs and job classes within the United States workforce. However, the dataset contains a multitude of missing observations and estimates, which somewhat limits its utility. Here, we propose a method by which to impute these missing values that leverages many of the inherent features present in the survey data, such as known population limit and correlations between occupations and tasks. An iterative regression fit, implemented with a recent version of XGBoost and executed across a set of simulated values drawn from the distribution described by the known values and their standard deviations reported in the survey, is the approach used to arrive at a distribution of predicted values for each missing estimate. This allows us to calculate a mean prediction and bound said estimate with a 95% confidence interval. We discuss the use of our method and how the resulting imputations can be utilized to inform and pursue future areas of study stemming from the data collected in the ORS. Finally, we conclude with an outline of WIGEM, a generalized version of our weighted, iterative imputation algorithm that could be applied to other contexts.


title: Natural Selection and Random Matrix Theory

abstract: We will study the relationship between two well-known theories, genetic evolution and random matrix theory in the context of many-body systems. We show that the time evolution of certain quantum mechanical toy models is similar to that of a living cell. It is also suggested that genetic evolution can be described by a random matrix theory with statistical distribution in which natural selection acts as a Gross-Witten-Wadia phase transition.


title: Towards a single scale-dependent Pomeron in holographic light-front QCD

abstract: The Pomeron Regge trajectory underlies the dynamics dependence of hadronic total cross sections and diffractive reactions at high energies. The physics of the Pomeron is closely related to the gluon distribution function and the gluon gravitational form factor of the target hadron. In this article we examine the scale dependence of the nonperturbative gluon distribution in the nucleon and the pion which was derived in a previous article [Phys. Rev. D 104, 114005 (2021)] in the framework of holographic light-front QCD and the Veneziano model. We argue that the QCD evolution of the gluon distribution function $g(x,\mu)$ to large $\mu^2$ leads to a single scale-dependent Pomeron. The resulting Pomeron trajectory $\alpha_P(t, \mu)$ not only depends on the momentum transfer squared $t$, but also on the physical scale $\mu$ of the amplitude, such as the virtuality $Q^2$ of the interacting photon in inclusive diffractive electroproduction. This can explain not only the $Q^2$ evolution of the proton structure function $F_2(x,Q^2)$ at small $x$, but also the observed energy and $Q^2$ dependence of high energy diffractive processes involving virtual photons up to LHC energies.


title: Pattern Formation in Thermal Convective Systems: Spatio-temporal Thermal Statistics, Emergent Flux, and Local Equilibrium

abstract: We discuss spatio-temporal pattern formation in two separate thermal convective systems. In the first system, hydrothermal waves (HTW) are modeled numerically in an annular channel. A temperature difference is imposed across the channel, which induces a surface tension gradient on the free surface of the fluid, leading to a surface flow towards the cold side. The flow pattern is axially symmetric along the temperature gradient with an internal circulation for a small temperature difference. This axially symmetric flow (ASF) becomes unstable beyond a given temperature difference threshold, and subsequently, symmetry-breaking flow, i.e., rotational oscillating waves or HTW, appears. For the second system, Rayleigh-B\'enard convection (RBC) is experimentally studied in the non-turbulent regime. When a thin film of liquid is heated, the competing forces of viscosity and buoyancy give rise to convective instabilities. This convective instability creates a spatio-temporal non-uniform temperature distribution on the surface of the fluid film. The surface temperature statistics are studied in both these systems as `order' and `disorder' phase separates. Although the mechanisms that give rise to convective instabilities are different in both cases, we find an agreement on the macroscopic nature of the thermal distributions in these emergent structures.

categories: physics.flu-dyn cond-mat.soft cond-mat.stat-mech nlin.AO nlin.PS


title: Analytic Mutual Information in Bayesian Neural Networks

abstract: Bayesian neural networks have successfully designed and optimized a robust neural network model in many application problems, including uncertainty quantification. However, with its recent success, information-theoretic understanding about the Bayesian neural network is still at an early stage. Mutual information is an example of an uncertainty measure in a Bayesian neural network to quantify epistemic uncertainty. Still, no analytic formula is known to describe it, one of the fundamental information measures to understand the Bayesian deep learning framework. In this paper, we derive the analytical formula of the mutual information between model parameters and the predictive output by leveraging the notion of the point process entropy. Then, as an application, we discuss the parameter estimation of the Dirichlet distribution and show its practical application in the active learning uncertainty measures by demonstrating that our analytical formula can improve the performance of active learning further in practice.


title: Pseudospectral continuation for aeroelastic stability analysis

abstract: This technical note is concerned with aeroelastic flutter problems: the analysis of aeroelastic systems undergoing airspeed-dependent dynamic instability. Existing continuation methods for parametric stability analysis are based on marching along an airspeed parameter until the flutter point is found - an approach which may waste computational effort on low-airspeed system behavior, before a flutter point is located and characterized. Here, we describe a pseudospectral continuation approach which instead marches outwards from the system's flutter points, from points of instability to points of increasing damping, allowing efficient characterization of the subcritical and supercritical behavior of the system. This approach ties together aeroelastic stability analysis and abstract linear algebra, and provides efficient methods for computing practical aeroelastic stability properties - for instance, flight envelopes based on maximum modal damping, and the location of borderline-stable zones.


title: Dark Matter Halos in Interacting Dark Energy Models: Formation History, Density Profile, Spin and Shape

abstract: The interacting dark energy (IDE) model, which considers the interaction between dark energy and dark matter, provides a natural mechanism to alleviate the coincidence problem and can also relieve the observational tensions under the $\Lambda$CDM model. Previous studies have put constraints on IDE models by observations of cosmic expansion history, cosmic microwave background and large-scale structures. However, these data are not yet enough to distinguish IDE models from $\Lambda$CDM effectively. Because the non-linear structure formation contains rich cosmological information, it can provide additional means to differentiate alternative models. In this paper, based on a set of $N$-body simulations for IDE models, we investigate the formation histories and properties of dark matter halos, and compare with their $\Lambda$CDM counterparts. For the model with dark matter decaying into dark energy and the parameters being the best-fit values from previous constraints, the structure formation is markedly slowed down, and the halos have systematically lower mass, looser internal structure, higher spin and anisotropy. This is inconsistent with the observed structure formation, and thus this model can be safely ruled out from the perspective of non-linear structure formation. Moreover, we find that the ratio of halo concentrations between IDE and $\Lambda$CDM counterparts depends sensitively on the interaction parameter and is independent of halo mass. This can act as a powerful probe to constrain IDE models. Our results concretely demonstrate that the interaction of the two dark components can affect the halo formation considerably, and therefore the constraints from non-linear structures are indispensable.


title: Optimal SQ Lower Bounds for Learning Halfspaces with Massart Noise

abstract: We give tight statistical query (SQ) lower bounds for learnining halfspaces in the presence of Massart noise. In particular, suppose that all labels are corrupted with probability at most $\eta$. We show that for arbitrary $\eta \in [0,1/2]$ every SQ algorithm achieving misclassification error better than $\eta$ requires queries of superpolynomial accuracy or at least a superpolynomial number of queries. Further, this continues to hold even if the information-theoretically optimal error $\mathrm{OPT}$ is as small as $\exp\left(-\log^c(d)\right)$, where $d$ is the dimension and $0 < c < 1$ is an arbitrary absolute constant, and an overwhelming fraction of examples are noiseless. Our lower bound matches known polynomial time algorithms, which are also implementable in the SQ framework. Previously, such lower bounds only ruled out algorithms achieving error $\mathrm{OPT} + \epsilon$ or error better than $\Omega(\eta)$ or, if $\eta$ is close to $1/2$, error $\eta - o_\eta(1)$, where the term $o_\eta(1)$ is constant in $d$ but going to 0 for $\eta$ approaching $1/2$. As a consequence, we also show that achieving misclassification error better than $1/2$ in the $(A,\alpha)$-Tsybakov model is SQ-hard for $A$ constant and $\alpha$ bounded away from 1.


title: Stone Duality for Topological Convexity Spaces

abstract: A convexity space is a set X with a chosen family of subsets (called convex subsets) that is closed under arbitrary intersections and directed unions. There is a lot of interest in spaces that have both a convexity space and a topological space structure. In this paper, we study the category of topological convexity spaces and extend the Stone duality between coframes and topological spaces to an adjunction between topological convexity spaces and sup-lattices. We factor this adjunction through the category of preconvexity spaces (somtimes called closure spaces).


title: Optically Active and Nonlinear Chiral Metamaterials

abstract: Advanced photonic nanostructures have enabled the maximization of synthetic chiroptic activities. The unique structuring of these building blocks has empowered chiral selective interactions with electromagnetic waves in plasmonic structures and dielectric media. Given the repertoire of optimized chiral surfaces in the literature and the ubiquity of chirality in the organic realm, the natural direction to consider is the operation of these devices in larger optical systems much like their chiral organic counterparts. In this review, we recapitulate advances in active and nonlinear chiral metamaterials. Many of the results, such as the magneto-chiral anisotropy and third-harmonic Rayleigh scattering optical activity, are relatively unknown members of the more conventional family of tuning methodology and nonlinear processes. We believe they are poised to play an instrumental role in designing advanced chiroptic systems for applications in biochemistry, valleytronics, spintronics, and chiral quantum optics.


title: Heralded single-photon source based on ensemble of Raman active molecules

abstract: Light with high mutual correlations at different frequencies can be used to create heralded single-photon sources, which may serve as the basic elements of existing quantum cryptography and quantum teleportation schemes. One of the important examples in natural systems of light with high mutual correlations is the light produced by spontaneous Raman scattering on an ensemble of molecules. In this paper, we investigate the possibility of using Raman light to create a heralded single-photon source. We show that when using Stokes scattered light for postselection of an anti-Stokes scattered light, the latter may posses the single-photon properties. We analyze the influence of the various negative factors on the characteristics of such a heralded single-photon source, which include a time delay between Stokes and anti-Stokes photons, the finiteness of the correlation radius of an external source, and background radiation. We show that the high purity of the single-photon source is preserved even when the flow of uncorrelated photons exceeds the flow of correlated photons in the scattered Raman light by an order of magnitude.


title: Spectral-PQ: A Novel Spectral Sensitivity-Orientated Perceptual Compression Technique for RGB 4:4:4 Video Data

abstract: There exists an intrinsic relationship between the spectral sensitivity of the Human Visual System (HVS) and colour perception; these intertwined phenomena are often overlooked in perceptual compression research. In general, most previously proposed visually lossless compression techniques exploit luminance (luma) masking including luma spatiotemporal masking, luma contrast masking and luma texture/edge masking. The perceptual relevance of color in a picture is often overlooked, which constitutes a gap in the literature. With regard to the spectral sensitivity phenomenon of the HVS, the color channels of raw RGB 4:4:4 data contain significant color-based psychovisual redundancies. These perceptual redundancies can be quantized via color channel-level perceptual quantization. In this paper, we propose a novel spatiotemporal visually lossless coding method named Spectral Perceptual Quantization (Spectral-PQ). With application for RGB 4:4:4 video data, Spectral-PQ exploits HVS spectral sensitivity-related color masking in addition to spatial masking and temporal masking; the proposed method operates at the Coding Block (CB) level and the Prediction Unit (PU) level in the HEVC standard. Spectral-PQ perceptually adjusts the Quantization Step Size (QStep) at the CB level if high variance spatial data in G, B and R CBs is detected and also if high motion vector magnitudes in PUs are detected. Compared with anchor 1 (HEVC HM 16.17 RExt), Spectral-PQ considerably reduces bitrates with a maximum reduction of approximately 81%. The Mean Opinion Score (MOS) in the subjective evaluations show that Spectral-PQ successfully achieves perceptually lossless quality.


title: Exponential ergodicity for a stochastic two-layer quasi-geostrophic model

abstract: Ergodic properties of a stochastic medium complexity model for atmosphere and ocean dynamics are analysed. More specifically, a two-layer quasi-geostrophic model for geophysical flows is studied, with the upper layer being perturbed by additive noise. This model is popular in the geosciences, for instance to study the effects of a stochastic wind forcing on the ocean. A rigorous mathematical analysis however meets with the challenge that in the model under study, the noise configuration is spatially degenerate as the stochastic forcing acts only on the top layer. Exponential convergence of solutions laws to the invariant measure is established, implying a spectral gap of the associated Markov semigroup on a space of H\"older continuous functions. The approach provides a general framework for generalised coupling techniques suitable for applications to dissipative SPDEs. In case of the two-layer quasi-geostrophic model, the results require the second layer to obey a certain passivity condition.


title: The Design, Construction, Operation and Performance of the Belle II Silicon Vertex Detector

abstract: The Silicon Vertex Detector of Belle II is a state-of-the-art tracking and vertexing system based on double-sided silicon strip sensors, designed and fabricated by a large international collaboration in the period 2012--2018. Since 2019 it has been in operation providing high quality data with a small number of defective channels (<1%), a large hit-finding efficiency (>99%), a good signal-to-noise ratio (well in excess of 10 for all sensor configurations and tracks). Together with the good control over the alignment, these are all essential factors to achieve good tracking reconstruction and physics performance. In this extended paper we try to document all the aspects of the SVD challenges and achievements, in the spirit of providing information to the broader community and help the development of high quality detector systems, which are fundamental tools to carry out physics research.


title: Supported Sets -- A New Foundation For Nominal Sets And Automata

abstract: The present work proposes and discusses the category of supported sets which provides a uniform foundation for nominal sets of various kinds, such as those for equality symmetry, for the order symmetry, and renaming sets. We show that all these differently flavoured categories of nominal sets are monadic over supported sets. Thus, supported sets provide a canonical finite way to represent nominal sets and the automata therein, e.g. register automata. Name binding in supported sets is modelled by a functor following the idea of de Bruijn indices. This functor lifts to the well-known abstraction functor in nominal sets. Together with the monadicity result, this gives rise to a transformation process that takes the finite representation of a register automaton in supported sets and transforms it into its configuration automaton in nominal sets.


title: Tidal disruption of star clusters in galaxy formation simulations

abstract: We investigate the evolution of the tidal field experienced by massive star clusters using cosmological simulations of Milky Way-sized galaxies. Clusters in our simulations experience the strongest tidal force in the first few hundred Myr after formation, when the maximum eigenvalue of the tidal tensor reaches several times $10^4$ Gyr$^{-2}$. After about 1 Gyr the tidal field plateaus at a lower value, with the median $\lambda_{\rm m} \sim 3 \times 10^3$ Gyr$^{-2}$. The fraction of time clusters spend in high tidal strength ($\lambda_{\rm m} > 3 \times 10^4$ Gyr$^{-2}$) regions also decreases with their age from $\sim$20% immediately after formation to less than 1% after 1 Gyr. At early ages both the in situ and ex situ clusters experience similar tidal fields, while at older ages the in situ clusters in general experience stronger tidal field due to their lower orbits in host galaxy. This difference is reflected in the survival of clusters: we looked into cluster disruption calculated in simulation runtime and found that ex situ star clusters of the same initial mass typically end up with higher bound fraction at the last available simulation snapshot than the in situ ones.


title: Mixed Precision GMRES-based Iterative Refinement with Recycling

abstract: With the emergence of mixed precision capabilities in hardware, iterative refinement schemes for solving linear systems $Ax=b$ have recently been revisited and reanalyzed in the context of three or more precisions. These new analyses show that under certain constraints on condition number, the LU factorization of the matrix can be computed in low precision without affecting the final accuracy. Another promising technique is GMRES-based iterative refinement, which, in contrast to the standard approach, use GMRES preconditioned by the low-precision triangular factors to solve for the approximate solution update in each refinement step. This more accurate solution method extends the range of problems which can be solved with a given combination of precisions. However, in certain settings, GMRES may require too many iterations per refinement step, making it potentially more expensive than simply recomputing the LU factors in a higher precision. Krylov subspace recycling is a well-known technique for reusing information across sequential invocations of a Krylov subspace method on systems with the same or a slowly changing coefficient matrix. In this work, we incorporate the idea of Krylov subspace recycling into a mixed precision GMRES-based iterative refinement solver. The insight is that in each refinement step, we call preconditioned GMRES on a linear system with the same coefficient matrix $A$, with only the right-hand side changing. In this way, the GMRES solves in subsequent refinement steps can be accelerated by recycling information obtained from the first step. We perform extensive numerical experiments on various random dense problems, Toeplitz problems (prolate matrices), and problems from real applications, which confirm the benefits of the recycling approach.


title: MMLatch: Bottom-up Top-down Fusion for Multimodal Sentiment Analysis

abstract: Current deep learning approaches for multimodal fusion rely on bottom-up fusion of high and mid-level latent modality representations (late/mid fusion) or low level sensory inputs (early fusion). Models of human perception highlight the importance of top-down fusion, where high-level representations affect the way sensory inputs are perceived, i.e. cognition affects perception. These top-down interactions are not captured in current deep learning models. In this work we propose a neural architecture that captures top-down cross-modal interactions, using a feedback mechanism in the forward pass during network training. The proposed mechanism extracts high-level representations for each modality and uses these representations to mask the sensory inputs, allowing the model to perform top-down feature masking. We apply the proposed model for multimodal sentiment recognition on CMU-MOSEI. Our method shows consistent improvements over the well established MulT and over our strong late fusion baseline, achieving state-of-the-art results.


title: Training Data Selection for Accuracy and Transferability of Interatomic Potentials

abstract: Advances in machine learning (ML) techniques have enabled the development of interatomic potentials that promise both the accuracy of first principles methods and the low-cost, linear scaling, and parallel efficiency of empirical potentials. Despite rapid progress in the last few years, ML-based potentials often struggle to achieve transferability, that is, to provide consistent accuracy across configurations that significantly differ from those used to train the model. In order to truly realize the promise of ML-based interatomic potentials, it is therefore imperative to develop systematic and scalable approaches for the generation of diverse training sets that ensure broad coverage of the space of atomic environments. This work explores a diverse-by-construction approach that leverages the optimization of the entropy of atomic descriptors to create a very large ($>2\cdot10^{5}$ configurations, $>7\cdot10^{6}$ atomic environments) training set for tungsten in an automated manner, i.e., without any human intervention. This dataset is used to train polynomial as well as multiple neural network potentials with different architectures. For comparison, a corresponding family of potentials were also trained on an expert-curated dataset for tungsten. The models trained to entropy-optimized data exhibited vastly superior transferability compared to the expert-curated models. Furthermore, while the models trained with heavy user input (i.e., domain expertise) yield the lowest errors when tested on similar configurations, out-sample predictions are dramatically more robust when the models are trained on a deliberately diverse set of training data. Herein we demonstrate the development of both accurate and transferable ML potentials using automated and data-driven approaches for generating large and diverse training sets.


abstract: Devising augmentations for graph contrastive learning is challenging due to their irregular structure, drastic distribution shifts, and nonequivalent feature spaces across datasets. We introduce LG2AR, Learning Graph Augmentations to Learn Graph Representations, which is an end-to-end automatic graph augmentation framework that helps encoders learn generalizable representations on both node and graph levels. LG2AR consists of a probabilistic policy that learns a distribution over augmentations and a set of probabilistic augmentation heads that learn distributions over augmentation parameters. We show that LG2AR achieves state-of-the-art results on 18 out of 20 graph-level and node-level benchmarks compared to previous unsupervised models under both linear and semi-supervised evaluation protocols. The source code will be released here: https://github.com/kavehhassani/lg2ar


title: The Image Deblurring Problem: Matrices, Wavelets, and Multilevel Methods

abstract: The image deblurring problem consists of reconstructing images from blur and noise contaminated available data. In this AMS Notices article, we provide an overview of some well known numerical linear algebra techniques that are use for solving this problem. In particular, we start by carefully describing how to represent images, the process of blurring an image and modeling different kind of added noise. Then, we present regularization methods such as Tikhonov (on the standard and general form), Total Variation and other variations with sparse and edge preserving properties. Additionally, we briefly overview some of the main matrix structures for the blurring operator and finalize presenting multilevel methods that preserve such structures. Numerical examples are used to illustrate the techniques described.


title: Imaging cloaked objects: diffraction tomography of realistic invisibility devices

abstract: Invisibility cloaks have become one of the most outstanding developments among the wide range of applications in the field of metamaterials. So far, most efforts in invisibility science have been devoted to achieving practically realizable cloak designs and to improving the effectiveness of these devices in reducing their scattering cross-section (SCS), a scalar quantity accounting for the total electromagnetic energy scattered by an object. In contrast, little attention has been paid to the opposite side of the technology: the development of more efficient techniques for the detection of invisibility devices. For instance, the SCS ignores the phase change introduced by the cloak, as well as the angular dependence of the incident and scattered waves. Here, we propose to take advantage of the smarter way in which diffraction tomography processes all this overlooked information to improve the efficiency in unveiling the presence of invisibility devices. We show that this approach not only results in a considerable sensitivity enhancement in the detection of different kinds of cloaks based on both scattering cancellation and transformation optics, but also enables us to even obtain images depicting the approximate shape and size of the cloak. As a side result, we find that scattering cancellation coatings may be used to enhance the contrast of tomographic images of closely-packed particles. Our work could be easily extended to the detection of sound cloaks.


title: Cache Allocation in Multi-Tenant Edge Computing via online Reinforcement Learning

abstract: We consider in this work Edge Computing (EC) in a multi-tenant environment: the resource owner, i.e., the Network Operator (NO), virtualizes the resources and lets third party Service Providers (SPs - tenants) run their services, which can be diverse and with heterogeneous requirements. Due to confidentiality guarantees, the NO cannot observe the nature of the traffic of SPs, which is encrypted. This makes resource allocation decisions challenging, since they must be taken based solely on observed monitoring information. We focus on one specific resource, i.e., cache space, deployed in some edge node, e.g., a base station. We study the decision of the NO about how to partition cache among several SPs in order to minimize the upstream traffic. Our goal is to optimize cache allocation using purely data-driven, model-free Reinforcement Learning (RL). Differently from most applications of RL, in which the decision policy is learned offline on a simulator, we assume no previous knowledge is available to build such a simulator. We thus apply RL in an \emph{online} fashion, i.e., the policy is learned by directly perturbing the actual system and monitoring how its performance changes. Since perturbations generate spurious traffic, we also limit them. We show in simulation that our method rapidly converges toward the theoretical optimum, we study its fairness, its sensitivity to several scenario characteristics and compare it with a method from the state-of-the-art.