title: Secure Information Flow Connections

abstract: Denning's lattice model provided secure information flow analyses with an intuitive mathematical foundation: the lattice ordering determines permitted flows. We examine how this framework may be extended to support the flow of information between autonomous organisations, each employing possibly quite different security lattices and information flow policies. We propose a connection framework that permits different organisations to exchange information while maintaining both security of information flow as well as their autonomy in formulating and maintaining security policies. Our prescriptive framework is based on the rigorous mathematical framework of Lagois connections proposed by Melton, together with a simple operational model for transferring object data between domains. The merit of this formulation is that it is simple, minimal, adaptable and intuitive. We show that our framework is semantically sound, by proving that the connections proposed preserve standard correctness notions such as non-interference. We then illustrate how Lagois theory also provides a robust framework and methodology for negotiating and maintaining secure agreements on information flow between autonomous organisations, even when either or both organisations change their security lattices. Composition and decomposition properties indicate support for a modular approach to secure flow frameworks in complex organisations. We next show that this framework extends naturally and conservatively to the Decentralised Labels Model of Myers et al. - a Lagois connection between the hierarchies of principals in two organisations naturally induces a Lagois connection between the corresponding security label lattices, thus extending the security guarantees ensured by the decentralised model to encompass bidirectional inter-organisational flows.


title: Kernel Dependence Network

abstract: We propose a greedy strategy to spectrally train a deep network for multi-class classification. Each layer is defined as a composition of linear weights with the feature map of a Gaussian kernel acting as the activation function. At each layer, the linear weights are learned by maximizing the dependence between the layer output and the labels using the Hilbert Schmidt Independence Criterion (HSIC). By constraining the solution space on the Stiefel Manifold, we demonstrate how our network construct (Kernel Dependence Network or KNet) can be solved spectrally while leveraging the eigenvalues to automatically find the width and the depth of the network. We theoretically guarantee the existence of a solution for the global optimum while providing insight into our network's ability to generalize.


title: Understanding Double Descent Requires a Fine-Grained Bias-Variance Decomposition

abstract: Classical learning theory suggests that the optimal generalization performance of a machine learning model should occur at an intermediate model complexity, with simpler models exhibiting high bias and more complex models exhibiting high variance of the predictive function. However, such a simple trade-off does not adequately describe deep learning models that simultaneously attain low bias and variance in the heavily overparameterized regime. A primary obstacle in explaining this behavior is that deep learning algorithms typically involve multiple sources of randomness whose individual contributions are not visible in the total variance. To enable fine-grained analysis, we describe an interpretable, symmetric decomposition of the variance into terms associated with the randomness from sampling, initialization, and the labels. Moreover, we compute the high-dimensional asymptotic behavior of this decomposition for random feature kernel regression, and analyze the strikingly rich phenomenology that arises. We find that the bias decreases monotonically with the network width, but the variance terms exhibit non-monotonic behavior and can diverge at the interpolation boundary, even in the absence of label noise. The divergence is caused by the \emph{interaction} between sampling and initialization and can therefore be eliminated by marginalizing over samples (i.e. bagging) \emph{or} over the initial parameters (i.e. ensemble learning).


title: Learning to Respond with Your Favorite Stickers: A Framework of Unifying Multi-Modality and User Preference in Multi-Turn Dialog

abstract: Stickers with vivid and engaging expressions are becoming increasingly popular in online messaging apps, and some works are dedicated to automatically select sticker response by matching the stickers image with previous utterances. However, existing methods usually focus on measuring the matching degree between the dialog context and sticker image, which ignores the user preference of using stickers. Hence, in this paper, we propose to recommend an appropriate sticker to user based on multi-turn dialog context and sticker using history of user. Two main challenges are confronted in this task. One is to model the sticker preference of user based on the previous sticker selection history. Another challenge is to jointly fuse the user preference and the matching between dialog context and candidate sticker into final prediction making. To tackle these challenges, we propose a \emph{Preference Enhanced Sticker Response Selector} (PESRS) model. Specifically, PESRS first employs a convolutional based sticker image encoder and a self-attention based multi-turn dialog encoder to obtain the representation of stickers and utterances. Next, deep interaction network is proposed to conduct deep matching between the sticker and each utterance. Then, we model the user preference by using the recently selected stickers as input, and use a key-value memory network to store the preference representation. PESRS then learns the short-term and long-term dependency between all interaction results by a fusion network, and dynamically fuse the user preference representation into the final sticker selection prediction. Extensive experiments conducted on a large-scale real-world dialog dataset show that our model achieves the state-of-the-art performance for all commonly-used metrics. Experiments also verify the effectiveness of each component of PESRS.


title: Blue-tilted inflationary tensor spectrum and reheating in the light of NANOGrav results

abstract: We discuss the possibility of explaining the recent NANOGrav results by inflationary gravitational waves (IGWs) with a blue-tilted primordial spectrum. Although such IGWs can account for the NANOGrav signal without contradicting the upper bound on the tensor-to-scalar ratio at the cosmic microwave background scale, the predicted spectrum is in strong tension with the upper bound on the amplitude of the stochastic gravitational wave background by big-bang nucleosynthesis (BBN) and the second LIGO-Virgo observation run. However, the thermal history of the Universe, such as reheating and late-time entropy production, affects the spectral shape of IGWs at high frequencies and permits evading the upper bounds. We show that, for the standard reheating scenario, when the reheating temperature is relatively low, a blue tensor spectrum can explain the recent NANOGrav signal without contradicting the BBN and the LIGO-Virgo constraints. We further find that, when one considers a late-time entropy production, the NANOGrav signal can be explained even for an instant reheating scenario.


title: Search for the brightest stars in galaxies outside the Local Group

abstract: This paper shows a technique for searching for bright massive stars in galaxies beyond the Local Group. To search for massive stars, we used the results of stellar photometry of the Hubble Space Telescope images using the DAOPHOT and DOLPHOT packages. The results of such searches are shown on the example of the galaxies DDO68, M94 and NGC1672. In the galaxy DDO68 the LBV star changes its brightness, and in M94 massive stars can be identified by the excess in the H${\alpha}$ band. For the galaxy NGC1672, we measured the distance for the first time by the TRGB method, which made it possible to determine the luminosities of the brightest stars, likely hypergiants, in the young star formation region. So far we have performed stellar photometry of HST images of 320 northern sky galaxies located at a distance below 12Mpc. This allowed us to identify 53 galaxies with probable hypergiants. Further photometric and spectral observations of these galaxies are planned to search for massive stars.


title: DNN-based Detectors for Massive MIMO Systems with Low-Resolution ADCs

abstract: Low-resolution analog-to-digital converters (ADCs) have been considered as a practical and promising solution for reducing cost and power consumption in massive Multiple-Input-Multiple-Output (MIMO) systems. Unfortunately, low-resolution ADCs significantly distort the received signals, and thus make data detection much more challenging. In this paper, we develop a new deep neural network (DNN) framework for efficient and low-complexity data detection in low-resolution massive MIMO systems. Based on reformulated maximum likelihood detection problems, we propose two model-driven DNN-based detectors, namely OBMNet and FBMNet, for one-bit and few-bit massive MIMO systems, respectively. The proposed OBMNet and FBMNet detectors have unique and simple structures designed for low-resolution MIMO receivers and thus can be efficiently trained and implemented. Numerical results also show that OBMNet and FBMNet significantly outperform existing detection methods.


title: On members of Lucas sequences which are either products of factorials or product of middle binomial coefficients and Catalan numbers

abstract: Let $\{U_n\}_{n\geq 0}$ be a Lucas sequence. Then the equation $$|U_n|=m_1!m_2!\cdots m_k!$$ with $1<m_1\leq m_2\leq \cdots\leq m_k$ implies $n\in \{1,2, 3, 4, 6, 8, 12\}$. Further the equation $$|U_n|=D_{m_1}D_{m_2}\cdots D_{m_k}, \qquad D_{m_i}\in \{B_{m_i}, C_{m_i}\}$$ with $1<m_1\leq m_2\leq \cdots\leq m_k$ implies $n\in \{1,2, 3, 4, 6, 8, 12, 16\}$. Here $B_m$ is the middle binomial coefficient $\binom{2m}{m}$ and $C_m$ is the Catalan number $\frac{1}{m+1}\binom{2m}{m}$.


title: Fighting an Infodemic: COVID-19 Fake News Dataset

abstract: Along with COVID-19 pandemic we are also fighting an `infodemic'. Fake news and rumors are rampant on social media. Believing in rumors can cause significant harm. This is further exacerbated at the time of a pandemic. To tackle this, we curate and release a manually annotated dataset of 10,700 social media posts and articles of real and fake news on COVID-19. We benchmark the annotated dataset with four machine learning baselines - Decision Tree, Logistic Regression, Gradient Boost, and Support Vector Machine (SVM). We obtain the best performance of 93.46% F1-score with SVM. The data and code is available at: https://github.com/parthpatwa/covid19-fake-news-dectection


title: Understanding Work Practices of Autonomous Agile Teams: A Social-psychological Review

abstract: The purpose of this paper is to suggest additional aspects of social psychology that could help when making sense of autonomous agile teams. To make use of well-tested theories in social psychology and instead see how they replicated and differ in the autonomous agile team context would avoid reinventing the wheel. This was done, as an initial step, through looking at some very common agile practices and relate them to existing findings in social-psychological research. The two theories found that I argue could be more applied to the software engineering context are social identity theory and group socialization theory. The results show that literature provides social-psychological reasons for the popularity of some agile practices, but that scientific studies are needed to gather empirical evidence on these under-researched topics. Understanding deeper psychological theories could provide a better understanding of the psychological processes when building autonomous agile team, which could then lead to better predictability and intervention in relation to human factors.


title: Gigantic magnetochiral anisotropy in the topological semimetal ZrTe5

abstract: Topological materials with broken inversion symmetry can give rise to nonreciprocal responses, such as the current rectification controlled by magnetic fields via magnetochiral anisotropy. Bulk nonreciprocal responses usually stem from relativistic corrections and are always very small. Here we report our discovery that ZrTe5 crystals in proximity to a topological quantum phase transition present gigantic magnetochiral anisotropy, which is the largest ever observed to date. We argue that a very low carrier density, inhomogeneities, and a torus-shaped Fermi surface induced by breaking of inversion symmetry in a Dirac material are central to explain this extraordinary property.


title: Safe trajectory of a piece moved by a robot

abstract: In this work, we propose a mathematical model for a physical problem based on the movement of a metal piece held by a robot. Using the principles of Kirchoff plate theory, a set of equations determining stresses and deformations caused during the motion, have been provided. We also discuss possible numerical treatment of these equations and finally, a solution to the one-dimensional analog of the problem has been presented.


title: Scalable Unsupervised Multi-Criteria Trajectory Segmentation and Driving Preference Mining

abstract: We present analysis techniques for large trajectory data sets that aim to provide a semantic understanding of trajectories reaching beyond them being point sequences in time and space. The presented techniques use a driving preference model w.r.t. road segment traversal costs, e.g., travel time and distance, to analyze and explain trajectories. In particular, we present trajectory mining techniques that can (a) find interesting points within a trajectory indicating, e.g., a via-point, and (b) recover the driving preferences of a driver based on their chosen trajectory. We evaluate our techniques on the tasks of via-point identification and personalized routing using a data set of more than 1 million vehicle trajectories collected throughout Denmark during a 3-year period. Our techniques can be implemented efficiently and are highly parallelizable, allowing them to scale to millions or billions of trajectories.


title: The fractal dimension of the Appalachian Trail

abstract: The Appalachian Trail (AT) is a 2193-mile-long hiking trail in the eastern United States. The trail has many bends and turns at different length scales, which gives it a nontrivial fractal dimension. Here I use GPS data from the Appalachian Trail Conservancy to estimate the fractal dimension of the AT. I find that, at length scales between $\sim 20$ m and $\sim 100$ km, the trail has a well-defined "divider dimension" of $\approx 1.08$. This dimension can be used to estimate the true hiking distance between two points, given the distance as estimated from a map with finite spatial resolution (e.g., Google Maps).


title: Contribution of Magnetic Reconnection Events to Energy Dissipation in Magnetosheath Turbulence

abstract: By analyzing the magnetosheath measurements from MMS, we obtain the statistical results for the contribution of magnetic reconnection (MR) events at electron scales to the energy dissipation of coherent structures. The Partial Variance of Increments (PVI) method is employed to find coherent structures in the magnetic field data. The current sheet structures with reversal of magnetic field components are further selected. We consider the following criteria to identify the MR events, such as current sheet with magnetic field reversal, significant energy dissipation, and evident electron outflow velocity. Statistically, for most MR events, their PVI values are larger than that of other types of coherent structures, and their energy dissipations are also stronger than that of others. However, due to the relatively small proportion of MR events, their contribution to coherent structures' energy dissipation is relatively trivial. If taken into account the dissipation of non-coherent structures, the MR's contribution to energy dissipation would be less. Hence, we suggest that MR events, though have strong dissipation locally, are not the major contributor to the energy dissipation in the magnetosheath. After analyzing the features of non-MR current sheets, we propose that non-MR current sheets are mainly coherent structures inherent to kinetic Alfv\'en fluctuations.


title: Occlusion-Aware Search for Object Retrieval in Clutter

abstract: We address the manipulation task of retrieving a target object from a cluttered shelf. When the target object is hidden, the robot must search through the clutter for retrieving it. Solving this task requires reasoning over the likely locations of the target object. It also requires physics reasoning over multi-object interactions and future occlusions. In this work, we present a data-driven hybrid planner for generating occlusion-aware actions in closed-loop. The hybrid planner explores likely locations of the occluded target object as predicted by a learned distribution from the observation stream. The search is guided by a heuristic trained with reinforcement learning to act on observations with occlusions. We evaluate our approach in different simulation and real-world settings (video available on https://youtu.be/dY7YQ3LUVQg). The results validate that our approach can search and retrieve a target object in near real time in the real world while only being trained in simulation.


title: Automatic Differentiation in PCF

abstract: We study the correctness of automatic differentiation (AD) in the context of a higher-order, Turing-complete language (PCF with real numbers), both in forward and reverse mode. Our main result is that, under mild hypotheses on the primitive functions included in the language, AD is almost everywhere correct, that is, it computes the derivative or gradient of the program under consideration except for a set of Lebesgue measure zero. Stated otherwise, there are inputs on which AD is incorrect, but the probability of randomly choosing one such input is zero. Our result is in fact more precise, in that the set of failure points admits a more explicit description: for example, in case the primitive functions are just constants, addition and multiplication, the set of points where AD fails is contained in a countable union of zero sets of non-identically-zero polynomials.


title: Fast state estimation under sensor attacks: a senor categorization approach

abstract: In a sensor network, some sensors usually provide the same or equivalent measurement information, which is not taken into account by the existing secure state estimation methods against sparse sensor attacks such that the computational efficiency of these methods needs to be further improved. In this paper, by considering the observation equivalence of sensor measurement, a concept of analytic sensor types is introduced based on the equivalent class to develop a fast state estimation algorithm. By testing the measurement data of a sensor type, the attack location information can be extracted to exclude some mismatching search candidates, without loss of estimation correctness. This confirms high speed performance of the proposed algorithm, since the number of sensor types is usually far less than the number of sensors.


title: Quantum corrections to tidal Love number for Schwarzschild black holes

abstract: A sum rule for tidal Love number is derived from quantum field theory computations, which relates tidal susceptibility of a spinless body to transition rates of its single graviton emission processes. An analogous sum rule for electromagnetism is given as an example, which is substantiated by comparing to the solved problem of the hydrogen atom. Based on the semiclassical Hawking radiation spectrum, a finite nonvanishing value for quantum corrections to the Love number of Schwarzschild black holes in general relativity is computed using the sum rule, which is known to classically vanish.


title: A Sub-Element Adaptive Shock Capturing Approach for Discontinuous Galerkin Methods

abstract: In this paper, a new strategy for a sub-element based shock capturing for discontinuous Galerkin (DG) approximations is presented. The idea is to interpret a DG element as a collection of data and construct a hierarchy of low to high order discretizations on this set of data, including a first order finite volume scheme up to the full order DG scheme. The different DG discretizations are then blended according to sub-element troubled cell indicators, resulting in a final discretization that adaptively blends from low to high order within a single DG element. The goal is to retain as much high order accuracy as possible, even in simulations with very strong shocks, as e.g. presented in the Sedov test. The framework retains the locality of the standard DG scheme and is hence well suited for a combination with adaptive mesh refinement and parallel computing. The numerical tests demonstrate the sub-element adaptive behavior of the new shock capturing approach and its high accuracy.


title: Effect of Short-Term Debt on Financial Growth of Non-Financial Firms Listed at Nairobi Securities Exchange

abstract: A significant number of the non-financial firms listed at Nairobi Securities Exchange (NSE) have been experiencing declining financial performance which deter investors from investing in such firms. The lenders are also not willing to lend to such firms. As such, the firms struggle to raise funds for their operations. Prudent financing decisions can lead to financial growth of the firm. The purpose of this study is to assess the effect of short-term debt on financial growth of non-financial firms listed at Nairobi Securities Exchange for a period of ten years from 2008 to 2017. Financial firms were excluded because of their specific sector characteristics and stringent regulatory framework. The study is guided by Agency Theory and Theory of Growth of the Firm. Explanatory research design was adopted. The target population of the study comprised of 45 non-financial firms listed at the NSE for a period of ten years from 2008 to 2017. The study conducted both descriptive statistics analysis and panel data analysis. The result indicates that, short term debt explains 45.99% and 25.6% of variations in financial growth as measured by growth in earnings per share and growth in market capitalization respectively. Short term debt positively and significantly influences financial growth measured using both growth in earnings per share and growth in market capitalization. The study recommends that, the management of non-financial firms listed at Nairobi Securities Exchange to employ financing means that can improve the earnings per share, market capitalization and enhance the value of the firm for the benefit of its stakeholders.


title: Testing Short Distance Anisotropy in Space

abstract: The isotropy of space is not a logical requirement but rather is an empirical question; indeed there is suggestive evidence that universe might be anisotropic. A plausible source of these anisotropies could be quantum gravity corrections. If these corrections happen to be between the electroweak scale and the Planck scale, then these anisotropies can have measurable consequences at short distances and their effects can be measured using ultra sensitive condensed matter systems. We investigate how such anisotropic quantum gravity corrections modify low energy physics through an anisotropic deformation of the Heisenberg algebra. We discuss how such anisotropies might be observed using a scanning tunneling microscope.


title: Evolution of cosmological perturbations in the models with interacting dynamical dark energy

abstract: Evolution of cosmological perturbations is considered in the model with dynamical dark energy which non-gravitationally interacts with dark matter. The dark energy equation of state parameter (EoS) is varying in time and is parameterized by its adiabatic sound speed. Such model of interacting dark energy has advantages over model with constant EoS, because it avoids non-adiabatic instabilities of dark energy at radiation dominated epoch for certain types of interaction in dark sector. The stability conditions for solutions of equations of dark energy perturbations were derived. The impact of strength of additional interaction between dark components on the evolution of density and velocity perturbations in them is analysed for quintessence and phantom types of dark energy.


title: On a conjecture regarding quantum hypothesis testing

abstract: In this MSc thesis I consider the asymptotic behaviour of the symmetric error in composite hypothesis testing. In the classical case, when the null and alternative hypothesis are finite sets of states, the best achievable symmetric error exponent is simply the the best achievable symmetric error exponent corresponding to the "worst case pair". The conjecture -- raised several years ago -- is that this remains true in the quantum case, too. This is known to be true in some special case. However, all of the known special cases are in some sense "too nice", e.g., have certain symmetries. Hoping to find a counter-example, in my thesis I consider a new special case, which on one hand is as "asymmetrical" as possible, yet still analytically computable. However, as it turns out from some involved computation, the conjecture is also true in this case, and thus gives further evidence to this conjecture.


title: Temporal nonlinear dynamics of plasmon-solitons, a Duffing oscillator-based approach

abstract: This paper deals with the temporal nonlinear dynamics of plasmon-solitons in a plasmonic waveguide. Duffing equation is recognized as the temporal part of the nonlinear amplitude equation governing the plasmonic waveguide. It is shown that Duffing oscillator waveforms stand for the temporal nonlinear dynamics of plasmon-soliton waves. The exchange of Lorentz-type bright and dark solitons energies gives rise to a Fano resonance. It is thus shown that the interaction of solitons and the formation of plasmon-solitons is inherently nonlinear. It is accordingly indicated that the nonlinear modulation of the plasmon-solitons is achievable via tuning the nonlinearity of the plasmonic waveguide.


title: Can We Detect Mastitis earlier than Farmers?

abstract: The aim of this study was to build a modelling framework that would allow us to be able to detect mastitis infections before they would normally be found by farmers through the introduction of machine learning techniques. In the making of this we created two different modelling framework's, one that works on the premise of detecting Sub Clinical mastitis infections at one Somatic Cell Count recording in advance called SMA and the other tries to detect both Sub Clinical mastitis infections aswell as Clinical mastitis infections at any time the cow is milked called AMA. We also introduce the idea of two different feature sets for our study, these represent different characteristics that should be taken into account when detecting infections, these were the idea of a cow differing to a farm mean and also trends in the lactation. We reported that the results for SMA are better than those created by AMA for Sub Clinical infections yet it has the significant disadvantage of only being able to classify Sub Clinical infections due to how we recorded Sub Clinical infections as being any time a Somatic Cell Count measurement went above a certain threshold where as CM could appear at any stage of lactation. Thus in some cases the lower accuracy values for AMA might in fact be more beneficial to farmers.


title: Two-point evolution equations for incompressible variable density turbulence

abstract: A derivation of the "exact" two-point equations analogous to those used as a basis for one-point Reynolds-Averaged Navier-Stokes turbulence model for variable density, incompressible turbulence. The purpose is to present the statistical correlations that must be dealt with to yield a two-point, or alternatively, a spectral model of variable density turbulence.


title: DeepDFT: Neural Message Passing Network for Accurate Charge Density Prediction

abstract: We introduce DeepDFT, a deep learning model for predicting the electronic charge density around atoms, the fundamental variable in electronic structure simulations from which all ground state properties can be calculated. The model is formulated as neural message passing on a graph, consisting of interacting atom vertices and special query point vertices for which the charge density is predicted. The accuracy and scalability of the model are demonstrated for molecules, solids and liquids. The trained model achieves lower average prediction errors than the observed variations in charge density obtained from density functional theory simulations using different exchange correlation functionals.


title: Template-Free Preparation of Thermoresponsive Mag-netic Cilia Compatible with Biological Conditions

abstract: Bio-inspired materials are commonly used in the development of functional devices. The fabrication of artificial cilia mimicking the biological functions has emerged as a promising strategy for fluid manipulation in miniaturized systems. In this study, we propose a different physicochemical insight for the preparation of magnetic cilia based on the temperature-triggered reversible assembly of coated iron oxide nanoparticles in a bio-compatible template-free approach. The length of the prepared cilia could be tuned between 10 and 100 microns reaching aspect ratios up to 100 in a very dense array of flexible structures with persistence lengths around 8 microns. Magnetic actuation of the cilia revealed robust structures (over several hours of actuation) with a wide range of bending amplitudes resulting from high susceptibility of the filaments. The results demonstrate that the proposed strategy is an efficient and versatile alternative for templated fabrication methods and producing cilia with remarkable characteristics and dimensions within the template-free approaches.


title: Drone Positioning for Visible Light Communication with Drone-Mounted LED and Camera

abstract: The world is often stricken by catastrophic disasters. On-demand drone-mounted visible light communication (VLC) networks are suitable for monitoring disaster-stricken areas for leveraging disaster-response operations. The concept of an image sensor-based VLC has also attracted attention in the recent past for establishing stable links using unstably moving drones. However, existing works did not sufficiently consider the one-to-many image sensor-based VLC system. Thus, this paper proposes the concept of a one-to-many image sensor-based VLC between a camera and multiple drone-mounted LED lights with a drone-positioning algorithm to avoid interference among VLC links. Multiple drones are deployed on-demand in a disaster-stricken area to monitor the ground and continuously send image data to a camera with image sensor-based visible light communication (VLC) links. The proposed idea is demonstrated with the proof-of-concept (PoC) implemented with drones that are equipped with LED panels and a 4K camera. As a result, we confirmed the feasibility of the proposed system.


title: DeepM&Mnet for hypersonics: Predicting the coupled flow and finite-rate chemistry behind a normal shock using neural-network approximation of operators

abstract: In high-speed flow past a normal shock, the fluid temperature rises rapidly triggering downstream chemical dissociation reactions. The chemical changes lead to appreciable changes in fluid properties, and these coupled multiphysics and the resulting multiscale dynamics are challenging to resolve numerically. Using conventional computational fluid dynamics (CFD) requires excessive computing cost. Here, we propose a totally new efficient approach, assuming that some sparse measurements of the state variables are available that can be seamlessly integrated in the simulation algorithm. We employ a special neural network for approximating nonlinear operators, the DeepONet, which is used to predict separately each individual field, given inputs from the rest of the fields of the coupled multiphysics system. We demonstrate the effectiveness of DeepONet by predicting five species in the non-equilibrium chemistry downstream of a normal shock at high Mach numbers as well as the velocity and temperature fields. We show that upon training, DeepONets can be over five orders of magnitude faster than the CFD solver employed to generate the training data and yield good accuracy for unseen Mach numbers within the range of training. Outside this range, DeepONet can still predict accurately and fast if a few sparse measurements are available. We then propose a composite supervised neural network, DeepM&Mnet, that uses multiple pre-trained DeepONets as building blocks and scattered measurements to infer the set of all seven fields in the entire domain of interest. Two DeepM&Mnet architectures are tested, and we demonstrate the accuracy and capacity for efficient data assimilation. DeepM&Mnet is simple and general: it can be employed to construct complex multiphysics and multiscale models and assimilate sparse measurements using pre-trained DeepONets in a "plug-and-play" mode.


title: Intra-Domain Task-Adaptive Transfer Learning to Determine Acute Ischemic Stroke Onset Time

abstract: Treatment of acute ischemic strokes (AIS) is largely contingent upon the time since stroke onset (TSS). However, TSS may not be readily available in up to 25% of patients with unwitnessed AIS. Current clinical guidelines for patients with unknown TSS recommend the use of MRI to determine eligibility for thrombolysis, but radiology assessments have high inter-reader variability. In this work, we present deep learning models that leverage MRI diffusion series to classify TSS based on clinically validated thresholds. We propose an intra-domain task-adaptive transfer learning method, which involves training a model on an easier clinical task (stroke detection) and then refining the model with different binary thresholds of TSS. We apply this approach to both 2D and 3D CNN architectures with our top model achieving an ROC-AUC value of 0.74, with a sensitivity of 0.70 and a specificity of 0.81 for classifying TSS < 4.5 hours. Our pretrained models achieve better classification metrics than the models trained from scratch, and these metrics exceed those of previously published models applied to our dataset. Furthermore, our pipeline accommodates a more inclusive patient cohort than previous work, as we did not exclude imaging studies based on clinical, demographic, or image processing criteria. When applied to this broad spectrum of patients, our deep learning model achieves an overall accuracy of 75.78% when classifying TSS < 4.5 hours, carrying potential therapeutic implications for patients with unknown TSS.


title: Affine Invariant Analysis of Frank-Wolfe on Strongly Convex Sets

abstract: It is known that the Frank-Wolfe (FW) algorithm, which is affine-covariant, enjoys accelerated convergence rates when the constraint set is strongly convex. However, these results rely on norm-dependent assumptions, usually incurring non-affine invariant bounds, in contradiction with FW's affine-covariant property. In this work, we introduce new structural assumptions on the problem (such as the directional smoothness) and derive an affine invariant, norm-independent analysis of Frank-Wolfe. Based on our analysis, we propose an affine invariant backtracking line-search. Interestingly, we show that typical backtracking line-searches using smoothness of the objective function surprisingly converge to an affine invariant step size, despite using affine-dependent norms in the step size's computation. This indicates that we do not necessarily need to know the set's structure in advance to enjoy the affine-invariant accelerated rate.


title: Deep learning architectures for inference of AC-OPF solutions

abstract: We present a systematic comparison between neural network (NN) architectures for inference of AC-OPF solutions. Using fully connected NNs as a baseline we demonstrate the efficacy of leveraging network topology in the models by constructing abstract representations of electrical grids in the graph domain, for both convolutional and graph NNs. The performance of the NN architectures is compared for regression (predicting optimal generator set-points) and classification (predicting the active set of constraints) settings. Computational gains for obtaining optimal solutions are also presented.

categories: cs.LG cs.SY eess.SP eess.SY physics.data-an


title: Self-Supervised Learning for Biological Sample Localization in 3D Tomographic Images

abstract: In synchrotron-based Computed Tomography (CT) there is a trade-off between spatial resolution, field of view and speed of positioning and alignment of samples. The problem is even more prominent for high-throughput tomography--an automated setup, capable of scanning large batches of samples without human interaction. As a result, in many applications, only 20-30% of the reconstructed volume contains the actual sample. Such data redundancy clutters the storage and increases processing time. Hence, an automated sample localization becomes an important practical problem. In this work, we describe two self-supervised losses designed for biological CT. We further demonstrate how to employ the uncertainty estimation for sample localization. This approach shows the ability to localize a sample with less than 1.5\% relative error and reduce the used storage by a factor of four. We also show that one of the proposed losses works reasonably well as a pre-training task for the semantic segmentation.


title: Vertex Fault-Tolerant Geometric Spanners for Weighted Points

abstract: Given a set S of n points, a weight function w to associate a non-negative weight to each point in S, a positive integer k \ge 1, and a real number \epsilon > 0, we present algorithms for computing a spanner network G(S, E) for the metric space (S, d_w) induced by the weighted points in S. The weighted distance function d_w on the set S of points is defined as follows: for any p, q \in S, d_w(p, q) is equal to w(p) + d_\pi(p, q) + w(q) if p \ne q, otherwise, d_w(p, q) is 0. Here, d_\pi(p, q) is the Euclidean distance between p and q if points in S are in \mathbb{R}^d, otherwise, it is the geodesic (Euclidean) distance between p and q. The following are our results: (1) When the weighted points in S are located in \mathbb{R}^d, we compute a k-vertex fault-tolerant (4+\epsilon)-spanner network of size O(k n). (2) When the weighted points in S are located in the relative interior of the free space of a polygonal domain \cal P, we detail an algorithm to compute a k-vertex fault-tolerant (4+\epsilon)-spanner network with O(\frac{kn\sqrt{h+1}}{\epsilon^2} \lg{n}) edges. Here, h is the number of simple polygonal holes in \cal P. (3) When the weighted points in S are located on a polyhedral terrain \cal T, we propose an algorithm to compute a k-vertex fault-tolerant (4+\epsilon)-spanner network, and the number of edges in this network is O(\frac{kn}{\epsilon^2} \lg{n}).


title: About the Teaching of Plane Motion of Rigid Bodies

abstract: The study of the motion of a rigid body on a plane (RBP motion) is usually one of the most challenging topics that students face in introductory physics courses. In this paper, we discuss a couple of problems which are typically used in basic physics courses, in order to highlight some aspects related to RBP motion which are not usually well understood by physics students. The first problem is a pendulum composed of a rod and disk. The angular frequency of the pendulum is calculated in two situations: disk fixed to the rod and disk free to spin. A detailed explanation of the change in the angular frequency from one case to another is given. The second problem is a ladder which slides touching a frictionless surface. We use this problem to highlight the fact that the contact forces applied by the surface perform translational and rotational work despite that the total mechanical energy of the ladder is conserved.


title: Gamma Radiation Source Localization for Micro Aerial Vehicles with a Miniature Single-Detector Compton Event Camera

abstract: A novel method for localization and estimation of compact sources of gamma radiation for Micro Aerial Vehicles (MAVs) is presented in this paper. The method is developed for a novel single-detector Compton camera, developed by the authors. The detector is extremely small and weighs only 40 g, which opens the possibility for use on sub-1 kg class of drones. The Compton camera uses the MiniPIX TPX3 CdTe event camera to measure Compton scattering products of incoming high-energy gamma photons. The 3D position and the sub-nanosecond time delay of the measured scattering products are used to reconstruct sets of possible directions to the source. The proposed method utilizes a filter for fusing the measurements and estimating the radiation source position during the flight. The computations are executed in real-time onboard and allow integration of the detector info into a fully-autonomous system. Moreover, the real-time nature of the estimator potentially allows estimating states of a moving radiation source. The proposed method was validated in simulations and demonstrated in a real-world experiment with a Cs137 radiation source. The approach can localize a gamma source without estimating the gradient or contours of radiation intensity, which opens possibilities for localizing sources in a cluttered and urban environment.


title: A Precedence-Driven Approach for Concurrent Model Synchronization Scenarios using Triple Graph Grammars

abstract: Concurrent model synchronization is the task of restoring consistency between two correlated models after they have been changed concurrently and independently. To determine whether such concurrent model changes conflict with each other and to resolve these conflicts taking domain- or user-specific preferences into account is highly challenging. In this paper, we present a framework for concurrent model synchronization algorithms based on Triple Graph Grammars (TGGs). TGGs specify the consistency of correlated models using grammar rules; these rules can be used to derive different consistency restoration operations. Using TGGs, we infer a causal dependency relation for model elements that enables us to detect conflicts non-invasively. Different kinds of conflicts are detected first and resolved by the subsequent conflict resolution process. Users configure the overall synchronization process by orchestrating the application of consistency restoration fragments according to several conflict resolution strategies to achieve individual synchronization goals. As proof of concept, we have implemented this framework in the model transformation tool eMoflon. Our initial evaluation shows that the runtime of our presented approach scales with the size of model changes and conflicts, rather than model size.


title: Generalization of Quasi-Newton Methods: Application to Robust Symmetric Multisecant Updates

abstract: Quasi-Newton techniques approximate the Newton step by estimating the Hessian using the so-called secant equations. Some of these methods compute the Hessian using several secant equations but produce non-symmetric updates. Other quasi-Newton schemes, such as BFGS, enforce symmetry but cannot satisfy more than one secant equation. We propose a new type of quasi-Newton symmetric update using several secant equations in a least-squares sense. Our approach generalizes and unifies the design of quasi-Newton updates and satisfies provable robustness guarantees.


title: A New Inference algorithm of Dynamic Uncertain Causality Graph based on Conditional Sampling Method for Complex Cases

abstract: Dynamic Uncertain Causality Graph(DUCG) is a recently proposed model for diagnoses of complex systems. It performs well for industry system such as nuclear power plants, chemical system and spacecrafts. However, the variable state combination explosion in some cases is still a problem that may result in inefficiency or even disability in DUCG inference. In the situation of clinical diagnoses, when a lot of intermediate causes are unknown while the downstream results are known in a DUCG graph, the combination explosion may appear during the inference computation. Monte Carlo sampling is a typical algorithm to solve this problem. However, we are facing the case that the occurrence rate of the case is very small, e.g. $10^{-20}$, which means a huge number of samplings are needed. This paper proposes a new scheme based on conditional stochastic simulation which obtains the final result from the expectation of the conditional probability in sampling loops instead of counting the sampling frequency, and thus overcomes the problem. As a result, the proposed algorithm requires much less time than the DUCG recursive inference algorithm presented earlier. Moreover, a simple analysis of convergence rate based on a designed example is given to show the advantage of the proposed method. % In addition, supports for logic gate, logic cycles, and parallelization, which exist in DUCG, are also addressed in this paper. The new algorithm reduces the time consumption a lot and performs 3 times faster than old one with 2.7% error ratio in a practical graph for Viral Hepatitis B.


title: A Gleason-Kahane-\.Zelazko theorem for reproducing kernel Hilbert spaces

abstract: We establish the following Hilbert-space analogue of the Gleason-Kahane-\.Zelazko theorem. If $\mathcal{H}$ is a reproducing kernel Hilbert space with a normalized complete Pick kernel, and if $\Lambda$ is a linear functional on $\mathcal{H}$ such that $\Lambda(1)=1$ and $\Lambda(f)\ne0$ for all cyclic functions $f\in\mathcal{H}$, then $\Lambda$ is multiplicative, in the sense that $\Lambda(fg)=\Lambda(f)\Lambda(g)$ for all $f,g\in\mathcal{H}$ such that $fg\in\mathcal{H}$. Moreover $\Lambda$ is automatically continuous. We give examples to show that the theorem fails if the hypothesis of a complete Pick kernel is omitted. We also discuss conditions under which $\Lambda$ has to be a point evaluation.


title: Hadamard multipliers on weighted Dirichlet spaces

abstract: The Hadamard product of two power series is obtained by multiplying them coefficientwise. In this paper we characterize those power series that act as Hadamard multipliers on all weighted Dirichlet spaces on the disk with superharmonic weights, and we obtain sharp estimates on the corresponding multiplier norms. Applications include an analogue of Fej\'er's theorem in these spaces, and a new estimate for the weighted Dirichlet integrals of dilates.


title: Linear polynomial approximation schemes in Banach holomorphic function spaces

abstract: Let $X$ be a Banach holomorphic function space on the unit disk. A linear polynomial approximation scheme for $X$ is a sequence of bounded linear operators $T_n:X\to X$ with the property that, for each $f\in X$, the functions $T_n(f)$ are polynomials converging to $f$ in the norm of the space. We completely characterize those spaces $X$ that admit a linear polynomial approximation scheme. In particular, we show that it is NOT sufficient merely that polynomials be dense in $X$.


title: Domain Adaptive Person Re-Identification via Coupling Optimization

abstract: Domain adaptive person Re-Identification (ReID) is challenging owing to the domain gap and shortage of annotations on target scenarios. To handle those two challenges, this paper proposes a coupling optimization method including the Domain-Invariant Mapping (DIM) method and the Global-Local distance Optimization (GLO), respectively. Different from previous methods that transfer knowledge in two stages, the DIM achieves a more efficient one-stage knowledge transfer by mapping images in labeled and unlabeled datasets to a shared feature space. GLO is designed to train the ReID model with unsupervised setting on the target domain. Instead of relying on existing optimization strategies designed for supervised training, GLO involves more images in distance optimization, and achieves better robustness to noisy label prediction. GLO also integrates distance optimizations in both the global dataset and local training batch, thus exhibits better training efficiency. Extensive experiments on three large-scale datasets, i.e., Market-1501, DukeMTMC-reID, and MSMT17, show that our coupling optimization outperforms state-of-the-art methods by a large margin. Our method also works well in unsupervised training, and even outperforms several recent domain adaptive methods.


title: Critical point for Bose-Einstein condensation of excitons in graphite

abstract: An exciton is an electron-hole pair bound by attractive Coulomb interaction. Short-lived excitons have been detected by a variety of experimental probes in numerous contexts. An excitonic insulator, a collective state of such excitons, has been more elusive. Here, thanks to Nernst measurements in pulsed magnetic fields, we show that in graphite there is a critical temperature (T = 9.2 K) and a critical magnetic field (B = 47 T) for Bose-Einstein condensation of excitons. At this critical field, hole and electron Landau sub-bands simultaneously cross the Fermi level and allow exciton formation. By quantifying the effective mass and the spatial separation of the excitons in the basal plane, we show that the degeneracy temperature of the excitonic fluid corresponds to this critical temperature. This identification would explain why the field-induced transition observed in graphite is not a universal feature of three-dimensional electron systems pushed beyond the quantum limit.


title: Parallel optically detected magnetic resonance spectrometer for dozens of single nitrogen-vacancy centers using laser-spot lattice

abstract: We develop a parallel optically detected magnetic resonance (PODMR) spectrometer to address, manipulate and read out an array of single nitrogen-vacancy (NV) centers in diamond in parallel. In this spectrometer, we use an array of micro-lens to generate 20 * 20 laser-spot lattice (LSL) on the objective focal plane, and then align the LSL with an array of single NV centers. The quantum states of NV centers are manipulated by a uniform microwave field from a {\Omega}-shape coplanar coil. As an experimental demonstration, we observe 80 NV centers in the field of view. Among them, magnetic resonance (MR) spectrums and Rabi oscillations of 18 NV centers along the external magnetic field are measured in parallel. These results can be directly used to realize parallel quantum sensing and multiple times speedup compared with the confocal technique. Regarding the nanoscale MR technique, PODMR will be crucial for high throughput single molecular MR spectrum and imaging.


title: Frenet-Serret analysis of helical Bloch modes in N-fold rotationally symmetric rings of coupled spiralling optical waveguides

abstract: The behavior of electromagnetic waves in chirally twisted structures is a topic of enduring interest, dating back at least to the invention in the 1940s of the microwave travelling wave tube amplifier and culminating in contemporary studies of chiral metamaterials, metasurfaces, and photonic crystal fibers (PCFs). Optical fibers with chiral microstructures, drawn from a spinning preform, have many useful properties, exhibiting for example circular birefringence and circular dichroism. It has recently been shown that chiral fibers with N fold rotationally symmetric (symmetry group CN) transverse microstructures support families of helical Bloch modes (HBMs), each of which consists of a superposition of azimuthal Bloch harmonics (or optical vortices). An example is a fiber with N coupled cores arranged in a ring around its central axis (N core single ring fiber). Although this type of fiber can be readily modelled using scalar coupled mode theory, a full description of its optical properties requires a vectorial analysis that takes account of the polarization state of the light particularly important in studies of circular and vortical birefringence. In this paper we develop, using an orthogonal two dimensional helicoidal coordinate system embedded in a cylindrical surface at constant radius, a rigorous vector coupled mode description of the fields using local Frenet Serret frames that rotate and twist with each of the N cores. The analysis places on a firm theoretical footing a previous HBM theory in which a heuristic approach was taken, based on physical intuition of the properties of Bloch waves. We believe this study provides clarity in what can sometimes be a rather difficult field, and will facilitate further exploration of real-world applications of these fascinating waveguiding systems.


title: Disentangling 3D Prototypical Networks For Few-Shot Concept Learning

abstract: We present neural architectures that disentangle RGB-D images into objects' shapes and styles and a map of the background scene, and explore their applications for few-shot 3D object detection and few-shot concept classification. Our networks incorporate architectural biases that reflect the image formation process, 3D geometry of the world scene, and shape-style interplay. They are trained end-to-end self-supervised by predicting views in static scenes, alongside a small number of 3D object boxes. Objects and scenes are represented in terms of 3D feature grids in the bottleneck of the network. We show that the proposed 3D neural representations are compositional: they can generate novel 3D scene feature maps by mixing object shapes and styles, resizing and adding the resulting object 3D feature maps over background scene feature maps. We show that classifiers for object categories, color, materials, and spatial relationships trained over the disentangled 3D feature sub-spaces generalize better with dramatically fewer examples than the current state-of-the-art, and enable a visual question answering system that uses them as its modules to generalize one-shot to novel objects in the scene.


title: Randomized Quaternion Singular Value Decomposition for Low-Rank Approximation

abstract: This paper presents a randomized quaternion singular value decomposition (QSVD) algorithm for low-rank matrix approximation problems, which are widely used in color face recognition, video compression, and signal processing problems. With quaternion normal distribution based random sampling, the randomized QSVD algorithm projects a high-dimensional data to a low-dimensional subspace and then identifies an approximate range subspace of the quaternion matrix. The key statistical properties of quaternion Wishart distribution are proposed and used to perform the approximation error analysis of the algorithm. Theoretical results show that the randomized QSVD algorithm can trace dominant singular value decomposition triplets of a quaternion matrix with acceptable accuracy. Numerical experiments also indicate the rationality of proposed theories. Applied to color face recognition problems, the randomized QSVD algorithm obtains higher recognition accuracies and behaves more efficient than the known Lanczos-based partial QSVD and a quaternion version of fast frequent directions algorithm.


title: Asteroid Resource Utilization: Ethical Concerns and Progress

abstract: As asteroid mining moves toward reality, the high bar to entering the business may limit participation and increase inequality, reducing or eliminating any benefit gained by marginalized people or developing nations. Consideration of ethical issues is urgently needed, as well as participation in international, not merely multilateral, solutions.


title: Hydro-mechanical simulation and analysis of induced seismicity for a hydraulic stimulation test at the Reykjanes geothermal field, Iceland

abstract: The hydraulic stimulation of the well RN-34 at the Reykjanes geothermal field in Iceland caused increased seismic activity near the well. Here, we use this as a case study for investigation on how seismic analysis can be combined with physics-based simulation studies to further understand injection-induced fault reactivation. The work presents new analysis of the seismic data combined with application of a recent simulation software for modeling of coupled hydromechanical processes and fault deformation caused by fluid injection. The simulation model incorporates an explicit model of the fault network based on geological characterization combined with insights from seismic analysis. The 3D faulted reservoir model is then calibrated based on injection data. Despite limited data, the work shows how seismic interpretations can be used in developing simulation models and, reciprocally, how the modeling can add to the seismic interpretations in analysis of dynamics.


title: Explaining Differences in Classes of Discrete Sequences

abstract: While there are many machine learning methods to classify and cluster sequences, they fail to explain what are the differences in groups of sequences that make them distinguishable. Although in some cases having a black box model is sufficient, there is a need for increased explainability in research areas focused on human behaviors. For example, psychologists are less interested in having a model that predicts human behavior with high accuracy and more concerned with identifying differences between actions that lead to divergent human behavior. This paper presents techniques for understanding differences between classes of discrete sequences. Approaches introduced in this paper can be utilized to interpret black box machine learning models on sequences. The first approach compares k-gram representations of sequences using the silhouette score. The second method characterizes differences by analyzing the distance matrix of subsequences. As a case study, we trained black box supervised learning methods to classify sequences of GitHub teams and then utilized our sequence analysis techniques to measure and characterize differences between event sequences of teams with bots and teams without bots. In our second case study, we classified Minecraft event sequences to infer their high-level actions and analyzed differences between low-level event sequences of actions.


title: FDNAS: Improving Data Privacy and Model Diversity in AutoML

abstract: To prevent the leakage of private information while enabling automated machine intelligence, there is an emerging trend to integrate federated learning and Neural Architecture Search (NAS). Although promising as it may seem, the coupling of difficulties from both two tenets makes the algorithm development quite challenging. In particular, how to efficiently search the optimal neural architecture directly from massive non-iid data of clients in a federated manner remains to be a hard nut to crack. To tackle this challenge, in this paper, by leveraging the advances in proxy-less NAS, we propose a Federated Direct Neural Architecture Search (FDNAS) framework that allows hardware-aware NAS from decentralized non-iid data of clients. To further adapt for various data distributions of clients, inspired by meta-learning, a cluster Federated Direct Neural Architecture Search (CFDNAS) framework is proposed to achieve client-aware NAS, in the sense that each client can learn a tailored deep learning model for its particular data distribution. Extensive experiments on real-world non-iid datasets show state-of-the-art accuracy-efficiency trade-offs for various hardware and data distributions of clients. Our codes will be released publicly upon paper acceptance.


title: Linear fractional group as Galois group

abstract: We compute all signatures of $PSL_2(\mathbb{F}_7)$, and $PSL_2(\mathbb{F}_{11})$ which classify all orientation preserving actions of the groups $PSL_2(\mathbb{F}_7)$, and $PSL_2(\mathbb{F}_{11})$ on compact, connected, orientable surfaces with orbifold genus $\geq 0$. This classification is well-grounded in the other branches of Mathematics like topology, smooth, and conformal geometry, algebraic categories, and it is also directly related to the inverse Galois problem.


title: Covariant formulation of gluon self-energy in presence of ellipsoidal anisotropy

abstract: In this work, a covariant formulation of the gluon self-energy in presence of ellipsoidal anisotropy is considered. It is shown that the general structure of the gluon self-energy can be written in terms of six linearly independent projection tensors. Similar to the spheroidal anisotropy, mass scales can be introduced for each of the collective modes considering the static limits. With a simplified ellipsoidal generalization of the Romatschke-Strickland form, the angular dependencies of the mass scales are studied. It is observed that, compared to the spheroidal case, additional unstable mode may appear in presence of ellipsoidal anisotropy depending upon the choice of the parameters.


title: A Scalable MIP-based Method for Learning Optimal Multivariate Decision Trees

abstract: Several recent publications report advances in training optimal decision trees (ODT) using mixed-integer programs (MIP), due to algorithmic advances in integer programming and a growing interest in addressing the inherent suboptimality of heuristic approaches such as CART. In this paper, we propose a novel MIP formulation, based on a 1-norm support vector machine model, to train a multivariate ODT for classification problems. We provide cutting plane techniques that tighten the linear relaxation of the MIP formulation, in order to improve run times to reach optimality. Using 36 data-sets from the University of California Irvine Machine Learning Repository, we demonstrate that our formulation outperforms its counterparts in the literature by an average of about 10% in terms of mean out-of-sample testing accuracy across the data-sets. We provide a scalable framework to train multivariate ODT on large data-sets by introducing a novel linear programming (LP) based data selection method to choose a subset of the data for training. Our method is able to routinely handle large data-sets with more than 7,000 sample points and outperform heuristics methods and other MIP based techniques. We present results on data-sets containing up to 245,000 samples. Existing MIP-based methods do not scale well on training data-sets beyond 5,500 samples.


title: Epidemic Plateau: A Phenomenon under Adaptive Prevention Strategies

abstract: Since the beginning of the COVID-19 spreading, the number of studies on the epidemic models increased dramatically. It is important for policymakers to know how the disease will spread and what are the effects of the policies and environment on the spreading. In this paper, we propose two extensions to the standard SIR model: (a) we consider the prevention measures adopted based on the current severity of the infection. Those measures are adaptive and change over time; (b) multiple cities and regions are considered, with population movements between those cities and regions, while taking into account that each region may have different prevention measures. Although the adaptive measures and mobility of the population were often observed during the pandemic, these effects are rarely explicitly modeled and studied in the classical epidemic models. The model we propose gives rise to a plateau phenomenon: the number of people infected by the disease stays at the same level during an extended period of time. We show what are conditions need to be met in order for the spreading to exhibit a plateau period in a single city. In addition, this phenomenon is interdependent when considering multiple cities. We verify from the real-world data that the plateau phenomenon does exist in many regions of the world in the current COVID-19 development. Finally, we provide theoretical analysis on the plateau phenomenon for the single-city model and derive a series of results on the emergence and the ending of the plateau, as well as on the height and length of the plateau. Our theoretical results match well with our experimental findings.


title: On flatness and coherence with respect to modules of flat dimension at most one

abstract: This paper introduces and studies homological properties of new classes of modules, namely, the $\mathcal F_1$-flat modules and the $\mathcal F_1^{\fp}$-flat modules, where $\mathcal F_1$ stands for the class of right modules of flat dimension at most one and $\mathcal F_1^{\fp}$ its subclass consisting of finitely presented elements. This leads us to introduce a new class of rings that we term $\mathcal F_1^{\fp}$-coherent rings as they behave nicely with respect to $\mathcal F_1^{\fp}$-flat modules as do coherent rings with respect to flat modules. The new class of $\mathcal F_1^{\fp}$-coherent rings turns out to be a large one and it includes coherent rings, perfect rings, semi-hereditary rings and all rings $R$ such that $\lim\limits_{\lr}\mathcal P_1=\mathcal F_1$. As a particular case of rings satisfying $\lim\limits_{\lr}\mathcal P_1=\mathcal F_1$ figures the important class of integral domains.


title: Fractional Fibonacci groups with an odd number of generators

abstract: The Fibonacci groups $F(n)$ are known to exhibit significantly different behaviour depending on the parity of $n$. We extend known results for $F(n)$ for odd $n$ to the family of Fractional Fibonacci groups $F^{k/l}(n)$. We show that for odd $n$ the group $F^{k/l}(n)$ is not the fundamental group of an orientable hyperbolic 3-orbifold of finite volume. We obtain results concerning the existence of torsion in the groups $F^{k/l}(n)$ (where $n$ is odd) paying particular attention to the groups $F^k(n)$ and $F^{k/l}(3)$, and observe consequences concerning asphericity of relative presentations of their shift extensions. We show that if $F^{k}(n)$ (where $n$ is odd) and $F^{k/l}(3)$ are non-cyclic 3-manifold groups then they are isomorphic to the direct product of the quaternion group $Q_8$ and a finite cyclic group.


title: Joint Sensing and Communication over Memoryless Broadcast Channels

abstract: A memoryless state-dependent broadcast channel (BC) is considered, where the transmitter wishes to convey two private messages to two receivers while simultaneously estimating the respective states via generalized feedback. The model at hand is motivated by a joint radar and communication system where radar and data applications share the same frequency band. For physically degraded BCs with i.i.d. state sequences, we characterize the capacity-distortion tradeoff region. For general BCs, we provide inner and outer bounds on the capacitydistortion region, as well as a sufficient condition when it is equal to the product of the capacity region and the set of achievable distortion. Interestingly, the proposed synergetic design significantly outperforms a conventional approach that splits the resource either for sensing or communication.


title: Outbreak diversity in epidemic waves propagating through distinct geographical scales

abstract: A central feature of an emerging infectious disease in a pandemic scenario is the spread through geographical scales and the impacts on different locations according to the adopted mitigation protocols. We investigated a stochastic epidemic model with the metapopulation approach in which patches represent municipalities. Contagion follows a stochastic compartmental model for municipalities; the latter, in turn, interact with each other through recurrent mobility. As a case of study, we consider the epidemic of COVID-19 in Brazil performing data-driven simulations. Properties of the simulated epidemic curves have very broad distributions across different geographical locations and scales, from states, passing through intermediate and immediate regions down to municipality levels. Correlations between delay of the epidemic outbreak and distance from the respective capital cities were predicted to be strong in several states and weak in others, signaling influences of multiple epidemic foci propagating towards the inland cities. Responses of different regions to a same mitigation protocol can vary enormously implying that the policies of combating the epidemics must be engineered according to the region' specificity but integrated with the overall situation. Real series of reported cases confirm the qualitative scenarios predicted in simulations. Even though we restricted our study to Brazil, the prospects and model can be extended to other geographical organizations with heterogeneous demographic distributions.


title: Deep Learning-based Cattle Activity Classification Using Joint Time-frequency Data Representation

abstract: Automated cattle activity classification allows herders to continuously monitor the health and well-being of livestock, resulting in increased quality and quantity of beef and dairy products. In this paper, a sequential deep neural network is used to develop a behavioural model and to classify cattle behaviour and activities. The key focus of this paper is the exploration of a joint time-frequency domain representation of the sensor data, which is provided as the input to the neural network classifier. Our exploration is based on a real-world data set with over 3 million samples, collected from sensors with a tri-axial accelerometer, magnetometer and gyroscope, attached to collar tags of 10 dairy cows and collected over a one month period. The key results of this paper is that the joint time-frequency data representation, even when used in conjunction with a relatively basic neural network classifier, can outperform the best cattle activity classifiers reported in the literature. With a more systematic exploration of neural network classifier architectures and hyper-parameters, there is potential for even further improvements. Finally, we demonstrate that the time-frequency domain data representation allows us to efficiently trade-off a large reduction of model size and computational complexity for a very minor reduction in classification accuracy. This shows the potential for our classification approach to run on resource-constrained embedded and IoT devices.


title: Gravity Field Mapping Using Laser Coupled Quantum Accelerometers in Space

abstract: The emergence of quantum technologies, including cold atom based accelerometers, offers an opportunity to improve the performances of space geodesy missions. In this context, CNES initiated an assessment study called GRICE (GRadiom\'etrie \`a Interf\'erom\`etres quantiques Corr\'el\'es pour l'Espace) in order to evaluate the contribution of cold atom technologies to space geodesy and to the end users of geodetic data. In this paper, we present mission scenario for gravity field mapping based on a long baseline gradiometer. The mission is based on a constellation of two satellites, flying at an altitude of 373 km, each equipped with a cold atom accelerometer with a sensitivity of $6 \times 10^{-10}$~m.s$^{-2}$.$\mathrm{\tau}^{-1/2}$. A laser link measures the distance between the two satellites and couples these two instruments in order to produce a correlated differential acceleration measurement. The main parameters, determining the performances of the payload, have been investigated. We carried out a general study of satellite architecture and simulations of the mission performances in terms of restitution of the gravity field. The simulations show that this concept would give its best performance in terms of monthly gravity fields recovery under 1000~km resolution. In the resolution band between 1000 and 222~km, the improvement of the GRICE gradient approach over the traditional range-rate approach is globally in the order of 10 to 25\%.


title: Adversarial Skill Learning for Robust Manipulation

abstract: Deep reinforcement learning has made significant progress in robotic manipulation tasks and it works well in the ideal disturbance-free environment. However, in a real-world environment, both internal and external disturbances are inevitable, thus the performance of the trained policy will dramatically drop. To improve the robustness of the policy, we introduce the adversarial training mechanism to the robotic manipulation tasks in this paper, and an adversarial skill learning algorithm based on soft actor-critic (SAC) is proposed for robust manipulation. Extensive experiments are conducted to demonstrate that the learned policy is robust to internal and external disturbances. Additionally, the proposed algorithm is evaluated in both the simulation environment and on the real robotic platform.


title: Suppression of Correlated Noise with Similarity-based Unsupervised Deep Learning

abstract: Image denoising is a prerequisite for downstream tasks in many fields. Low-dose and photon-counting computed tomography (CT) denoising can optimize diagnostic performance at minimized radiation dose. Supervised deep denoising methods are popular but require paired clean or noisy samples that are often unavailable in practice. Limited by the independent noise assumption, current unsupervised denoising methods cannot process correlated noises as in CT images. Here we propose the first-of-its-kind similarity-based unsupervised deep denoising approach, referred to as Noise2Sim, that works in a nonlocal and nonlinear fashion to suppress not only independent but also correlated noises. Theoretically, Noise2Sim is asymptotically equivalent to supervised learning methods under mild conditions. Experimentally, Nosie2Sim recovers intrinsic features from noisy low-dose CT and photon-counting CT images as effectively as or even better than supervised learning methods on practical datasets visually, quantitatively and statistically. Noise2Sim is a general unsupervised denoising approach and has great potential in diverse applications.


title: Regularized Stochastic Team Problems

abstract: In this paper, we introduce regularized stochastic team problems. Under mild assumptions, we prove that there exists an unique fixed point of the best response operator, where this unique fixed point is the optimal regularized team decision rule. Then, we establish an asynchronous distributed algorithm to compute this optimal strategy. We also provide a bound that shows how the optimal regularized team decision rule performs in the original stochastic team problem.


title: Computing sets from all infinite subsets

abstract: A set is introreducible if it can be computed by every infinite subset of itself. Such a set can be thought of as coding information very robustly. We investigate introreducible sets and related notions. Our two main results are that the collection of introreducible sets is $\Pi^1_1$-complete, so that there is no simple characterization of the introreducible sets; and that every introenumerable set has an introreducible subset.


title: Experimental quantification of atomically-resolved HAADF-STEMimages using EDX

abstract: Atomically-resolved mappings of the indium composition in InGaN/GaN multi-quantum wellstructures have been obtained by quantifying the contrast in HAADF-STEM. The quantificationprocedure presented here does not rely on computation-intensive simulations, but rather uses EDXmeasurements to calibrate the HAADF-STEM contrast. The histogram of indium compositionsobtained from the mapping provides unique insights into the growth of InGaN: the transitionfrom GaN to InGaN and vice versa occurs in discreet increments of composition; each incrementcorresponds to one monolayer of the interface, indicating that nucleation takes longer than thelateral growth of the step. Strain-state analysis is also performed by applying Peak-Pair Analysisto the positions of the atomic columns identified the quantification of the contrast. The strainmappings yield an estimate of the composition in good agreement with the one obtained fromquantified HAADF-STEM, albeit with a lower precision. Possible improvements to increase theprecision of the strain mappings are discussed, opening potential pathways for the quantification ofarbitrary quaternary alloys at atomic scales.


title: Singular Boundary Conditions for Sturm--Liouville Operators via Perturbation Theory

abstract: We show that all self-adjoint extensions of semi-bounded Sturm--Liouville operators with general limit-circle endpoint(s) can be obtained via an additive singular form bounded self-adjoint perturbation of rank equal to the deficiency indices, say $d\in\{1,2\}$. This characterization generalizes the well-known analog for semi-bounded Sturm--Liouville operators with regular endpoints. Explicitly, every self-adjoint extension of the minimal operator can be written as \begin{align*} \boldsymbol{A}_\Theta=\boldsymbol{A}_0+{\bf B}\Theta{\bf B}^*, \end{align*} where $\boldsymbol{A}_0$ is a distinguished self-adjoint extension and $\Theta$ is a self-adjoint linear relation in $\mathbb{C}^d$. The perturbation is singular in the sense that it does not belong to the underlying Hilbert space but is form bounded with respect to $\boldsymbol{A}_0$, i.e. it belongs to $\mathcal{H}_{-1}(\boldsymbol{A}_0)$. The construction of a boundary triple and compatible boundary pair for the symmetric operator ensure that the perturbation is well-defined and self-adjoint extensions are in a one-to-one correspondence with self-adjoint relations $\Theta$. As an example, self-adjoint extensions of the classical symmetric Jacobi differential equation (which has two limit-circle endpoints) are obtained and their spectra are analyzed with tools both from the theory of boundary triples and perturbation theory.


title: Multi-modal excitation to model the Quasi-Biennial Oscillation

abstract: The Quasi-Biennial Oscillation (QBO) of stratospheric winds is the most striking example of mean-flow generation and reversal by the non-linear interactions of internal waves. Previous studies have used an idealized monochromatic forcing to investigate the QBO. Here we instead force a more realistic continuous wave spectrum. Unexpectedly, spreading the wave energy across a wide frequency range leads to more regular oscillations. We also find that different forcing spectra can yield the same QBO. Multi-modal wave forcing is thus essential for understanding wave/mean-flow interactions in nature.


title: Phase boundary near a magnetic percolation transition

abstract: Motivated by recent experimental observations [Phys. Rev. 96, 020407 (2017)] on hexagonal ferrites, we revisit the phase diagrams of diluted magnets close to the lattice percolation threshold. We perform large-scale Monte Carlo simulations of XY and Heisenberg models on both simple cubic lattices and lattices representing the crystal structure of the hexagonal ferrites. Close to the percolation threshold $p_c$, we find that the magnetic ordering temperature $T_c$ depends on the dilution $p$ via the power law $T_c \sim |p-p_c|^\phi$ with exponent $\phi=1.09$, in agreement with classical percolation theory. However, this asymptotic critical region is very narrow, $|p-p_c| \lesssim 0.04$. Outside of it, the shape of the phase boundary is well described, over a wide range of dilutions, by a nonuniversal power law with an exponent somewhat below unity. Nonetheless, the percolation scenario does not reproduce the experimentally observed relation $T_c \sim (x_c -x)^{2/3}$ in PbFe$_{12-x}$Ga$_x$O$_{19}$. We discuss the generality of our findings as well as implications for the physics of diluted hexagonal ferrites.


title: Measuring and characterizing the line profile of HARPS with a laser frequency comb

abstract: Aims. We study the 2D spectral line profile of HARPS (High Accuracy Radial Velocity Planet Searcher), measuring its variation with position across the detector and with changing line intensity. The characterization of the line profile and its variations are important for achieving the precision of the wavelength scales of 10^{-10} or 3.0 cm/s necessary to detect Earth-twins in the habitable zone around solar-like stars. Methods. We used a laser frequency comb (LFC) with unresolved and unblended lines to probe the instrument line profile. We injected the LFC light (attenuated by various neutral density filters) into both the object and the reference fibres of HARPS, and we studied the variations of the line profiles with the line intensities. We applied moment analysis to measure the line positions, widths, and skewness as well as to characterize the line profile distortions induced by the spectrograph and detectors. Based on this, we established a model to correct for point spread function distortions by tracking the beam profiles in both fibres. Results. We demonstrate that the line profile varies with the position on the detector and as a function of line intensities. This is consistent with a charge transfer inefficiency (CTI) effect on the HARPS detector. The estimate of the line position depends critically on the line profile, and therefore a change in the line amplitude effectively changes the measured position of the lines, affecting the stability of the wavelength scale of the instrument. We deduce and apply the correcting functions to re-calibrate and mitigate this effect, reducing it to a level consistent with photon noise.


title: Did Hurricane Katrina Reduce Mortality?

abstract: In a recent article in the American Economic Review, Tatyana Deryugina and David Molitor (DM) analyzed the effect of Hurricane Katrina on the mortality of elderly and disabled residents of New Orleans. The authors concluded that Hurricane Katrina improved the eight-year survival rate of elderly and disabled residents of New Orleans by 3% and that most of this decline in mortality was due to declines in mortality among those who moved to places with lower mortality. In this article, I provide a critical assessment of the evidence provided by DM to support their conclusions. There are three main problems. First, DM generally fail to account for the fact that people of different ages, races or sex will have different probabilities of dying as time goes by, and when they do allow for this, results change markedly. Second, DM do not account for the fact that residents in New Orleans are likely to be selected non-randomly on the basis of health because of the relatively high mortality rate in New Orleans compared to the rest of the country. Third, there is considerable evidence that among those who moved from New Orleans, the destination chosen was non-random. Finally, DM never directly assessed changes in mortality of those who moved, or stayed, in New Orleans before and after Hurricane Katrina. These problems lead me to conclude that the evidence presented by DM does not support their inferences.


title: Input-to-state stability of a scalar conservation law with nonlocal velocity

abstract: In this paper, we study input-to-state stability (ISS) of an equilibrium for a scalar conservation law with nonlocal velocity and measurement error arising in a highly re-entrant manufacturing system. By using a suitable Lyapunov function, we prove sufficient and necessary conditions on ISS. We also analyze the numerical discretization of ISS for a discrete scalar conservation law with nonlocal velocity and measurement error. A suitable discretized Lyapunov function is also analyzed to provide ISS of an equilibrium for the numerical approximation. Finally, we show numerical simulations to validate the theoretical findings.


title: Fluid dynamics simulations show that facial masks can suppress the spread of COVID-19 in indoor environments

abstract: The Coronavirus disease outbreak of 2019 has been causing significant loss of life and unprecedented economical loss throughout the world. Social distancing and face masks are widely recommended around the globe in order to protect others and prevent the spread of the virus through breathing, coughing, and sneezing. To expand the scientific underpinnings of such recommendations, we carry out high-fidelity computational fluid dynamics simulations of unprecedented resolution and realism to elucidate the underlying physics of saliva particulate transport during human cough with and without facial masks. Our simulations: (a) are carried out under both a stagnant ambient flow (indoor) and a mild unidirectional breeze (outdoor); (b) incorporate the effect of human anatomy on the flow; (c) account for both medical and non-medical grade masks; and (d) consider a wide spectrum of particulate sizes, ranging from 10 micro m to 300 micro m. We show that during indoor coughing some saliva particulates could travel up to 0.48 m, 0.73 m, and 2.62 m for the cases with medical-grade, non-medical grade, and without facial masks, respectively. Thus, in indoor environments either medical or non-medical grade facial masks can successfully limit the spreading of saliva particulates to others. Under outdoor conditions with a unidirectional mild breeze, however, leakage flow through the mask can cause saliva particulates to be entrained into the energetic shear layers around the body and transported very fast at large distances by the turbulent flow, thus, limiting the effectiveness of facial masks.


title: Underspecification Presents Challenges for Credibility in Modern Machine Learning

abstract: ML models often exhibit unexpectedly poor behavior when they are deployed in real-world domains. We identify underspecification as a key reason for these failures. An ML pipeline is underspecified when it can return many predictors with equivalently strong held-out performance in the training domain. Underspecification is common in modern ML pipelines, such as those based on deep learning. Predictors returned by underspecified pipelines are often treated as equivalent based on their training domain performance, but we show here that such predictors can behave very differently in deployment domains. This ambiguity can lead to instability and poor model behavior in practice, and is a distinct failure mode from previously identified issues arising from structural mismatch between training and deployment domains. We show that this problem appears in a wide variety of practical ML pipelines, using examples from computer vision, medical imaging, natural language processing, clinical risk prediction based on electronic health records, and medical genomics. Our results show the need to explicitly account for underspecification in modeling pipelines that are intended for real-world deployment in any domain.


title: Bayes Security: A Not So Average Metric

abstract: Security system designers favor worst-case security metrics, such as those derived from differential privacy (DP), due to the strong guarantees they provide. On the downside, these guarantees result in a high penalty on the system's performance. In this paper, we study Bayes security, a security metric inspired by the cryptographic advantage. Similarly to DP, Bayes security i) is independent of an adversary's prior knowledge, ii) it captures the worst-case scenario for the two most vulnerable secrets (e.g., data records); and iii) it is easy to compose, facilitating security analyses. Additionally, Bayes security iv) can be consistently estimated in a black-box manner, contrary to DP, which is useful when a formal analysis is not feasible; and v) provides a better utility-security trade-off in high-security regimes because it quantifies the risk for a specific threat model as opposed to threat-agnostic metrics such as DP. We formulate a theory around Bayes security, and we provide a thorough comparison with respect to well-known metrics, identifying the scenarios where Bayes Security is advantageous for designers.


title: Molecular Structure, Quantum Coherence and Solvent Effects on the Ultrafast Electron Transport in BODIPY--C$_{60}$ Derivatives

abstract: Molecular systems containing donor-bridge-acceptor sites or molecular antennas constitute promising candidates for organic photovoltaic device implementation. Photo-induced electron transfer in multi-chromophore molecular systems is defined by a subtle interaction between the donor and the molecular bridge, and by the system-solvent coupling. Here, we address the computation of quantum properties such as population inversion and electron transfer in molecular photo-systems composed of fulleroisoxazoline, fulleropyrrolidine, BODIPY and Zn-porphyrin, as well as their system-solvent ultrafast dynamics. The molecular complexes are modelled as two- and three-site systems, and we use the density functional theory (DFT) for obtaining the site energies required in the construction of the open system diabatic Hamiltonians relevant to the computation of the electron transfer. The site energies and electronic couplings are calculated by using a continuous polarizable model that allow for the analysis of different solvent environments, and the site-to-site couplings are computed by means of the generalized Mulliken Hush method at the DFT level. We find that the stabilization energy of the charge transfer states exhibit a significant variation for a compound embedded in different polar environments, and thus the effect due to the solvent has been analyzed for the specific cases of Methanol, THF and Toluene. We show that the incorporation of a molecular bridge generates the creation of an intermediate state that plays a crucial role in the charge transfer process by defining $\Lambda$ or cascaded-type energy schemes; this affects the asymptotic value of the transfer rate and favors the cascaded-type configuration.


title: Structure and dynamics of responsive colloids with dynamical polydispersity

abstract: Dynamical polydispersity in single-particle properties, for example a fluctuating particle size, shape, charge density, etc., is intrinsic to responsive colloids (RCs), such as biomacromolecules or microgels, but is typically not resolved in coarse-grained mesoscale simulations. Here, we present Brownian Dynamics (BD) simulations of suspensions of RCs modeling soft hydrogel colloids, for which the size of the individual particles is an explicitly resolved (Gaussian) degree of freedom and dynamically responds to the local interacting environment. We calculate the liquid structure, emergent size distributions, long-time diffusion, and property (size) relaxation kinetics for a wide range of densities and intrinsic property relaxation times in the canonical ensemble. Comparison to interesting reference cases, such as conventional polydisperse suspensions with a frozen parent distribution, or conventional monodisperse systems interacting with an effective pair potential for one fixed size, shows a significant spread in the structure and dynamics. The differences, most apparent in the high density regimes, are due to many-body correlations and the dynamical coupling between property and translation in RC systems, not explicitly accounted for in the conventional treatments. In particular, the translational diffusion in the RC systems is surprisingly close to the free (single RC) diffusion, mainly due to a cancellation of crowding and size compression effects. We show that an effective monodisperse pair potential can be constructed that describes the many-body correlations reasonably well by convoluting the RC pair potential with the density-dependent emergent size distributions and using a mean effective diffusion constant.


title: Quantum simulation of three-body interactions in weakly driven quantum systems

abstract: The realization of effective Hamiltonians featuring many-body interactions beyond pairwise coupling would enable the quantum simulation of central models underpinning topological physics and quantum computation. We overcome crucial limitations of perturbative Floquet engineering and discuss the highly accurate realization of a purely three-body Hamiltonian in superconducting circuits and molecular nanomagnets.


title: Ap\'ery Limits: Experiments and Proofs

abstract: An important component of Ap\'ery's proof that $\zeta (3)$ is irrational involves representing $\zeta (3)$ as the limit of the quotient of two rational solutions to a three-term recurrence. We present various approaches to such Ap\'ery limits and highlight connections to continued fractions as well as the famous theorems of Poincar\'e and Perron on difference equations. In the spirit of Jon Borwein, we advertise an experimental-mathematics approach by first exploring in detail a simple but instructive motivating example. We conclude with various open problems.


title: Upper bounds for Betti numbers from constraints on the Hilbert function

abstract: We describe an algorithm for finding sharp upper bounds for the total Betti numbers of a saturated ideal given certain constraints on its Hilbert function. This algorithm is implemented in the Macaulay2 package, MaxBettiNumbers, along with variations that produce ideals with maximal total Betti numbers.


title: The journey of the German expedition to observe the transit of Venus on December 9, 1874 at the Kerguelen Islands and its sojourn there

abstract: Astronomers are the surveyors of the sky. They have always dedicated much time and resources to determining the scale of distances in the universe. Today, the Hubble constant; yesterday, the solar parallar. The following account, attributed to Ladislas Weinek, a modest actor of this scientific epic, describes one of its numerous chapters.


title: Beating classical heuristics for the binary paint shop problem with the quantum approximate optimization algorithm

abstract: The binary paint shop problem (BPSP) is an APX-hard optimization problem of the automotive industry. In this work, we show how to use the Quantum Approximate Optimization Algorithm (QAOA) to find solutions of the BPSP and demonstrate that QAOA with constant depth is able to beat classical heuristics on average in the infinite size limit $n\rightarrow\infty$. For the BPSP, it is known that no classical algorithm can exist which approximates the problem in polynomial runtime. We introduce a BPSP instance which is hard to solve with QAOA, and numerically investigate its performance and discuss QAOA's ability to generate approximate solutions. We complete our studies by running first experiments of small-sized instances on a trapped-ion quantum computer through AWS Braket.


title: Heat current anticorrelation effects leading to thermal conductivity reduction in nanoporous Si

abstract: Prevailing nanostructuring strategies focus on increasing phonon scattering and reducing the mean-free-path of phonons across the spectrum. In nanoporous Si materials, for example, boundary scattering reduces thermal conductivity drastically. In this work, we identify an unusual anticorrelated specular phonon scattering effect which can result in additional reductions in thermal conductivity of up to ~ 80% for specific nanoporous geometries. We further find evidence that this effect has its origin in heat trapping between large pores with narrow necks. As the heat becomes trapped between the pores, phonons undergo multiple specular reflections such that their contribution to the thermal conductivity is partly undone. We find this effect to be wave-vector dependent at low temperatures. We use large-scale molecular dynamics simulations, wave packet analysis, as well as an analytical model to illustrate the anticorrelation effect, evaluate its impact on thermal conductivity, and detail how it can be controlled to manipulate phonon transport in nanoporous materials.


title: Multiscale Control of Stackelberg Games

abstract: We present a linear--quadratic Stackelberg game with a large number of followers and we also derive the mean field limit of infinitely many followers. The relation between optimization and mean-field limit is studied and conditions for consistency are established. Finally, we propose a numerical method based on the derived models and present numerical results.


title: Early-Dynastic tables from Southern Mesopotamia, or the multiple facets of the quantification of surfaces

abstract: How were surfaces evaluated before the invention of the sexagesimal place value notation in Mesopotamia? This chapter examines a group of five tablets containing tables for surfaces of squares and rectangles dated to the Early Dynastic period (ca. 2600-2350 BCE) and unearthed in southern Mesopotamia. In order to capture the methods used by ancient scribes to quantify surfaces, special attention is paid to the layout and organization of the tables, as well as to the way in which measurement values are written down. It is argued that these methods vary according to the dimensions of the squares or rectangles concerned: the quantification of small surfaces does not use the same mathematical tools as the quantification of large parcels of land. The chapter shows a reciprocal influence between the metrological systems adopted by the ancient scribes and the methods of calculation of surface they implemented. Some methods may reflect ancient land-surveying practices, and others may testify the emergence of new mathematical concepts applied to all kinds of surfaces, large or small. Ultimately, several different conceptualizations of the notion of surface emerge from the examination of these tables.


title: Dark Energy Survey Year 3 Results: Photometric Data Set for Cosmology

abstract: We describe the Dark Energy Survey (DES) photometric data set assembled from the first three years of science operations to support DES Year 3 cosmology analyses, and provide usage notes aimed at the broad astrophysics community. Y3 Gold improves on previous releases from DES, Y1 Gold and Data Release 1 (DES DR1), presenting an expanded and curated data set that incorporates algorithmic developments in image detrending and processing, photometric calibration, and object classification. Y3 Gold comprises nearly 5000 square degrees of grizY imaging in the south Galactic cap, including nearly 390 million objects, with depth reaching S/N ~ 10 for extended objects up to $i_{AB}\sim 23.0$, and top-of-the-atmosphere photometric uniformity $< 3$ mmag. Compared to DR1, photometric residuals with respect to Gaia are reduced by $50\%$, and per-object chromatic corrections are introduced. Y3 Gold augments DES DR1 with simultaneous fits to multi-epoch photometry for more robust galaxy color measurements and corresponding photometric redshift estimates. Y3 Gold features improved morphological star-galaxy classification with efficiency $>98\%$ and purity $>99\%$ for galaxies with $19 < i_{AB} < 22.5$. Additionally, it includes per-object quality information, and accompanying maps of the footprint coverage, masked regions, imaging depth, survey conditions, and astrophysical foregrounds that are used to select the cosmology analysis samples. This paper will be complemented by online resources.