title: Internal $\infty$-Categorical Models of Dependent Type Theory: Towards 2LTT Eating HoTT

abstract: Using dependent type theory to formalise the syntax of dependent type theory is a very active topic of study and goes under the name of "type theory eating itself" or "type theory in type theory." Most approaches are at least loosely based on Dybjer's categories with families (CwF's) and come with a type CON of contexts, a type family TY indexed over it modelling types, and so on. This works well in versions of type theory where the principle of unique identity proofs (UIP) holds. In homotopy type theory (HoTT) however, it is a long-standing and frequently discussed open problem whether the type theory "eats itself" and can serve as its own interpreter. The fundamental underlying difficulty seems to be that categories are not suitable to capture a type theory in the absence of UIP. In this paper, we develop a notion of $\infty$-categories with families ($\infty$-CwF's). The approach to higher categories used relies on the previously suggested semi-Segal types, with a new construction of identity substitutions that allow for both univalent and non-univalent variations. The type-theoretic universe as well as the internalised syntax are models, although it remains a conjecture that the latter is initial. To circumvent the known unsolved problem of constructing semisimplicial types, the definition is presented in two-level type theory (2LTT). Apart from introducing $\infty$-CwF's, the paper explains the shortcomings of 1-categories in type theory without UIP as well as the difficulties of and approaches to internal higher-dimensional categories.


title: Model extraction from counterfactual explanations

abstract: Post-hoc explanation techniques refer to a posteriori methods that can be used to explain how black-box machine learning models produce their outcomes. Among post-hoc explanation techniques, counterfactual explanations are becoming one of the most popular methods to achieve this objective. In particular, in addition to highlighting the most important features used by the black-box model, they provide users with actionable explanations in the form of data instances that would have received a different outcome. Nonetheless, by doing so, they also leak non-trivial information about the model itself, which raises privacy issues. In this work, we demonstrate how an adversary can leverage the information provided by counterfactual explanations to build high-fidelity and high-accuracy model extraction attacks. More precisely, our attack enables the adversary to build a faithful copy of a target model by accessing its counterfactual explanations. The empirical evaluation of the proposed attack on black-box models trained on real-world datasets demonstrates that they can achieve high-fidelity and high-accuracy extraction even under low query budgets.


title: Linear Optical Approach to Supersymmetric Dynamics

abstract: The concept of supersymmetry developed in particle physics has been applied to various fields of modern physics. In quantum mechanics, the supersymmetric systems refer to the systems involving two supersymmetric partner Hamiltonians, whose energy levels are degeneracy except one of the systems has an extra ground state possibly, and the eigenstates of the partner systems can be mapped onto each other. Recently, an interferometric scheme has been proposed to show this relationship in ultracold atoms [Phys. Rev.A 96, 043624 (2017)]. Here this approach is generalized to linear optics for observing the supersymmetric dynamics with photons. The time evolution operator is simulated approximately via Suzuki-Trotter expansion with considering the realization of the kinetic and potential terms separately. The former is realized through the diffraction nature of light and the later is implemented using phase plate. Additionally, we propose an interferometric approach which can be implemented perfectly using amplitude alternator to realize the non-unitary operator. The numerical results show that our scheme is universal and can be realized with current technologies.


title: Laver Trees in the Generalized Baire Space

abstract: We prove that any suitable generalization of Laver forcing to the space $ \kappa^\kappa$, for uncountable regular $\kappa$, necessarily adds a Cohen $\kappa$-real. We also study a dichotomy and an ideal naturally related to generalized Laver forcing. Using this dichotomy, we prove the following stronger result: if $ \kappa^{<\kappa}=\kappa$, then every $<\kappa$-distributive tree forcing on $\kappa^\kappa$ adding a dominating $\kappa$-real which is the image of the generic under a continuous function in the ground model, adds a Cohen $\kappa$-real. This is a contribution to the study of generalized Baire spaces and answers a question from arXiv:1611.08140


title: Robust Homomorphic Video Hashing

abstract: The Internet has been weaponized to carry out cybercriminal activities at an unprecedented pace. The rising concerns for preserving the privacy of personal data while availing modern tools and technologies is alarming. End-to-end encrypted solutions are in demand for almost all commercial platforms. On one side, it seems imperative to provide such solutions and give people trust to reliably use these platforms. On the other side, this creates a huge opportunity to carry out unchecked cybercrimes. This paper proposes a robust video hashing technique, scalable and efficient in chalking out matches from an enormous bulk of videos floating on these commercial platforms. The video hash is validated to be robust to common manipulations like scaling, corruptions by noise, compression, and contrast changes that are most probable to happen during transmission. It can also be transformed into the encrypted domain and work on top of encrypted videos without deciphering. Thus, it can serve as a potential forensic tool that can trace the illegal sharing of videos without knowing the underlying content. Hence, it can help preserve privacy and combat cybercrimes such as revenge porn, hateful content, child abuse, or illegal material propagated in a video.


title: The effect of luminosity outbursts on the protoplanetary disk dynamics

abstract: Aims: Response of a protoplanetary disk to luminosity bursts of various duration is studied with the purpose to determine the effect of the bursts on the strength and sustainability of gravitational instability in the disk. A special emphasis is paid to the spatial distribution of gas and grown dust (from 1 mm to a few cm) during and after the burst. Methods: Numerical hydrodynamics simulations using the FEOSAD code were employed to study the dynamics of gas and dust in the thin-disk limit. Dust-to-gas friction including back reaction and dust growth were also considered. Bursts of various duration (from 100 to 500 yr) were initiated in accordance with a thermally ignited magnetorotational instability. Luminosity curves for constant- and declining-magnitude bursts were adopted to represent two typical limiting cases for FU-Orionis-type eruptions. Results: The short-term effect of the burst is to reduce the strength of gravitational instability by heating and expanding the disk. The longest bursts with duration comparable to the revolution period of the spiral can completely dissolve the original two-armed spiral pattern in the gas disk by the end of the burst, while the shortest bursts only weaken the spiral pattern. The reaction of grown dust to the burst is somewhat different. The spiral-like initial distribution with deep cavities in the inter-armed regions transforms into a ring-like distribution with deep gaps. This transformation is most expressed for the longest-duration bursts. The long-term effect of the burst depends on the initial disk conditions at the onset of the burst. In some cases, vigorous disk fragmentation sets in several thousand years after the burst, which was absent in the model without the bursts (abridged).


title: Compactness of a restricted X-ray transform

abstract: We show that the X-ray transform restricted along the moment curve possesses extremizers and that $L^p$-normalized extremizing sequences are pre-compact modulo symmetry.


title: Odd-spin glueballs, AdS/QCD and information entropy

abstract: Odd-spin glueballs in the dynamical AdS/QCD model are scrutinized, in the paradigm of the configurational entropy (CE). Configurational-entropic Regge trajectories, that relate the CE underlying odd-spin glueballs to their mass spectra and spin, are then engendered. They predict the mass spectra of odd-spin glueballs, besides pointing towards the configurational stability of odd-spin glueball resonances. The exponential modified dilaton with logarithmic anomalous dimensions comprises the most suitable choice to derive the mass spectra of odd-spin glueballs, compatible to lattice QCD. It is then used in a hybrid paradigm that takes both lattice QCD and the AdS/QCD correspondence into account.


title: Improving the Usability of Virtual Reality Neuron Tracing with Topological Elements

abstract: Researchers in the field of connectomics are working to reconstruct a map of neural connections in the brain in order to understand at a fundamental level how the brain processes information. Constructing this wiring diagram is done by tracing neurons through high-resolution image stacks acquired with fluorescence microscopy imaging techniques. While a large number of automatic tracing algorithms have been proposed, these frequently rely on local features in the data and fail on noisy data or ambiguous cases, requiring time-consuming manual correction. As a result, manual and semi-automatic tracing methods remain the state-of-the-art for creating accurate neuron reconstructions. We propose a new semi-automatic method that uses topological features to guide users in tracing neurons and integrate this method within a virtual reality (VR) framework previously used for manual tracing. Our approach augments both visualization and interaction with topological elements, allowing rapid understanding and tracing of complex morphologies. In our pilot study, neuroscientists demonstrated a strong preference for using our tool over prior approaches, reported less fatigue during tracing, and commended the ability to better understand possible paths and alternatives. Quantitative evaluation of the traces reveals that users' tracing speed increased, while retaining similar accuracy compared to a fully manual approach.


title: First principles analysis of precipitation in Mg-Zn alloys

abstract: Precipitation in Mg-Zn alloys was analyzed by means of first principles calculations. Formation energies of symmetrically distinct hcp Mg1-xZnx (0 < x < 1) configurations were determined and potential candidates for Guinier-Preston zones coherent with the matrix were identified from the convex hull. The most likely structures were ranked depending on the interface energy along the basal plane. In addition, the formation energy and vibrational entropic contributions of several phases reported experimentally (Mg4Zn7, MgZn2 cubic, MgZn2 hexagonal, Mg21Zn25 and Mg2Zn11) were calculated. The formation energies of Mg4Zn7, MgZn2 cubic, and MgZn2 hexagonal Laves phases were very close because they were formed by different arrangements of rhombohedral and hexagonal lattice units. It was concluded that \beta_1^' precipitates were formed by a mixture of all of them. Nevertheless, the differences in the geometrical arrangements led to variations in the entropic energy contributions which determined the high temperature stability. It was found that the MgZn2 hexagonal Laves phase is the most stable phase at high temperature and, thus, \beta_1^' precipitates tend to transform into the \beta_2^' (MgZn2 hexagonal) precipitates with higher aging temperature or longer aging times. Finally, the equilibrium \beta phase (Mg21Zn25) was found to be a long-range order that precipitates the last one on account of the kinetic processes necessary to trigger the transformation from a short-range order phase \beta_2^' to \beta .


title: On characteristic rank for matrix and tensor completion

abstract: In this lecture note, we discuss a fundamental concept, referred to as the {\it characteristic rank}, which suggests a general framework for characterizing the basic properties of various low-dimensional models used in signal processing. Below, we illustrate this framework using two examples: matrix and three-way tensor completion problems, and consider basic properties include identifiability of a matrix or tensor, given partial observations. In this note, we consider cases without observation noise to illustrate the principle.


title: PolSIRD: Modeling Epidemic Spread under Intervention Policies

abstract: Epidemic spread in a population is traditionally modeled via compartmentalized models which represent the free evolution of disease in absence of any intervention policies. In addition, these models assume full observability of disease cases and do not account for under-reporting. We present a mathematical model, namely PolSIRD, which accounts for the under-reporting by introducing an observation mechanism. It also captures the effects of intervention policies on the disease spread parameters by leveraging intervention policy data along with the reported disease cases. Furthermore, we allow our recurrent model to learn the initial hidden state of all compartments end-to-end along with other parameters via gradient-based training. We apply our model to the spread of the recent global outbreak of COVID-19 in the United States, where our model outperforms the methods employed by the CDC in predicting the spread. We also provide counterfactual simulations from our model to analyze the effect of lifting the intervention policies prematurely and our model correctly predicts the second wave of the epidemic.


title: The Dynamics, Destruction, and Survival of Supernova-Formed Dust Grains

abstract: Observations have demonstrated that supernovae efficiently produce dust. This is consistent with the hypothesis that supernovae and asymptotic giant branch stars are the primary producers of dust in the Universe. However, there has been a longstanding question of how much of the dust detected in the interiors of young supernova remnants can escape into the interstellar medium. We present new hydrodynamical calculations of the evolution of dust grains that were formed in dense ejecta clumps within a Cas A-like remnant. We follow the dynamics of the grains as they decouple from the gas after their clump is hit by the reverse shock. They are subsequently subject to destruction by thermal and kinetic sputtering as they traverse the remnant. Grains that are large enough ($\sim 0.25\,\mu$m for silicates and $\sim 0.1\,\mu$m for carbonaceous grains) escape into the interstellar medium while smaller grains get trapped and destroyed. However, grains that reach the interstellar medium still have high velocities, and are subject to further destruction as they are slowed down. We find that for initial grain size distributions that include large ($\sim 0.25 - 0.5\,\mu$m) grains, 10--20\% of silicate grains can survive, while 30--50\% of carbonaceous grains survive even when the initial size distribution cuts off at smaller ($0.25\,\mu$m) sizes. For a 19 M$_{\odot}$ star similar to the progenitor of Cas A, up to 0.1 M$_{\odot}$ of dust can survive if the dust grains formed are large. Thus we show that supernovae under the right conditions can be significant sources of interstellar dust.


title: Author Mentions in Science News Reveal Widespread Disparities Across Name-inferred Ethnicities

abstract: Media outlets play a key role in spreading scientific knowledge to the general public and raising the profile of researchers among their peers. Yet, how journalists choose to present researchers in their stories is poorly understood. Using a comprehensive dataset of 223,587 news stories from 288 U.S. outlets reporting on 100,486 research papers across all areas of science, we investigate if the authors' ethnicities, as inferred from names, are associated with whether journalists explicitly mention them by name. By focusing on research papers news outlets chose to cover, our analysis reduces concerns that differences in name mentions are driven by differences in research quality or newsworthiness. We find substantial disparities in name mention rates across ethnically-distinctive names. Researchers with non-Anglo names, especially those with East Asian and African names, are significantly less likely to be mentioned in news stories covering their research, even when comparing stories from a particular news outlet reporting on publications in a particular scientific venue on a particular research topic. The disparities are not fully explained by authors' affiliation locations, suggesting that pragmatic factors such as difficulties in scheduling interviews play only a partial role. Furthermore, among U.S.-based authors, journalists more often use authors' institutions instead of names when referring to non-Anglo-named authors, suggesting that journalists' rhetorical choices are also key. Overall, this study finds evidence of ethnic disparities in how researchers are described in the media coverage of their research, likely affecting thousands of non-Anglo-named scholars in our data alone.


title: Estimation of marriage incidence rates by combining two cross-sectional retrospective designs: Event history analysis of two dependent processes

abstract: The aim of this work is to develop methods for studying the determinants of marriage incidence using marriage histories collected under two different types of retrospective cross-sectional study designs. These designs are: sampling of ever married women before the cross-section, a prevalent cohort, and sampling of women irrespective of marital status, a general cross-sectional cohort. While retrospective histories from a prevalent cohort do not identify incidence rates without parametric modelling assumptions, the rates can be identified when combined with data from a general cohort. Moreover, education, a strong endogenous covariate, and marriage processes are correlated. Hence, they need to be modelled jointly in order to estimate the marriage incidence. For this purpose, we specify a multi-state model and propose a likelihood-based estimation method. We outline the assumptions under which a likelihood expression involving only marriage incidence parameters can be derived. This is of particular interest when either retrospective education histories are not available or related parameters are not of interest. Our simulation results confirm the gain in efficiency by combining data from the two designs, while demonstrating how the parameter estimates are affected by violations of the assumptions used in deriving the simplified likelihood expressions. Two Indian National Family Health Surveys are used as motivation for the methodological development and to demonstrate the application of the methods.


title: Chui's conjecture in Bergman spaces

abstract: We solve Chui's conjecture on the simplest fractions (i.e., sums of Cauchy kernels with unit coefficients) in weighted (Hilbert) Bergman spaces. Namely, for a wide class of weights, we prove that for every $N$, the simplest fractions with $N$ poles on the unit circle have minimal norm if and only if the poles are equispaced on the circle. We find sharp asymptotics of these norms. Furthermore, we describe the closure of the simplest fractions in weighted Bergman spaces, using an $L^2$ version of Thompson's theorem on dominated approximation by simplest fractions.


title: Limit groups over coherent right-angled Artin groups

abstract: A new class of groups $\mathcal{C}$, containing all coherent RAAGs and all toral relatively hyperbolic groups, is defined. It is shown that, for a group $G$ in the class $\mathcal{C}$, the $\mathbb{Z}[t]$-exponential group $G^{\mathbb{Z}[t]}$ may be constructed as an iterated centraliser extension. Using this fact, it is proved that $G^{\mathbb{Z}[t]}$ is fully residually $G$ (i.e. it has the same universal theory as $G$) and so its finitely generated subgroups are limit groups over $G$. If $\mathbb{G}$ is a coherent RAAG, then the converse also holds - any limit group over $\mathbb{G}$ embeds into $\mathbb{G}^{\mathbb{Z}[t]}$. Moreover, it is proved that limit groups over $\mathbb{G}$ are finitely presented, coherent and CAT$(0)$, so in particular have solvable word and conjugacy problems.


title: Effect of interactions in the interference pattern of Bose Einstein condensates

abstract: Understanding the effect of interactions in the phase evolution of expanding atomic Bose Einstein condensates is fundamental to describe the basic phenomenon of matter wave interference. Many theoretical and experimental works tackled this problem, always with the implicit assumption that the mutual interaction between two expanding condensates rigidly modifies the phase evolution through an effective force. In this paper, we present a combined experimental and theoretical investigation of the interference profile of expanding $^{87}$Rb condensates, with a specific focus on the effect of interactions. We come to the different conclusion that the mutual interaction produces local modifications of the condensate phase only in the region where the wavepackets overlap.


title: Phytoplankton temporal strategies increase entropy production in a marine food web model

abstract: We develop a trait-based model founded on the hypothesis that biological systems evolve and organize to maximize entropy production by dissipating chemical and electromagnetic potentials over longer time scales than abiotic processes by implementing temporal strategies. A marine food web consisting of phytoplankton, bacteria and consumer functional groups is used to explore how temporal strategies, or the lack there of, change entropy production in a shallow pond that receives a continuous flow of reduced organic carbon plus inorganic nitrogen and illumination from solar radiation with diel and seasonal dynamics. Results show that a temporal strategy that employs an explicit circadian clock produces more entropy than a passive strategy that uses internal carbon storage or a balanced growth strategy that requires phytoplankton to grow with fixed stoichiometry. When the community is forced to operate at high specific growth rates near 2 d-1, the optimization-guided model selects for phytoplankton ecotypes that exhibit complementary for winter versus summer environmental conditions to increase entropy production. We also present a new type of trait-based modeling where trait values are determined by maximizing entropy production rather than by random selection.


title: IoT-based Contact Tracing Systems for Infectious Diseases: Architecture and Analysis

abstract: The recent COVID-19 pandemic has become a major threat to human health and well-being. Non-pharmaceutical interventions such as contact tracing solutions are important to contain the spreads of COVID-19-like infectious diseases. However, current contact tracing solutions are fragmented with limited use of sensing technologies and centered on monitoring the interactions between individuals without an analytical framework for evaluating effectiveness. Therefore, we need to first explore generic architecture for contact tracing in the context of today's Internet of Things (IoT) technologies based on a broad range of applicable sensors. A new architecture for IoT based solutions to contact tracing is proposed and its overall effectiveness for disease containment is analyzed based on the traditional epidemiological models with the simulation results. The proposed work aims to provide a framework for assisting future designs and evaluation of IoT-based contact tracing solutions and to enable data-driven collective efforts on combating current and future infectious diseases.


title: Gravitational wave signatures from discrete flavor symmetries

abstract: Non-Abelian discrete symmetries have been widely used to explain the patterns of lepton masses and flavor mixing. In these models, a given symmetry is assumed at a high scale and then is spontaneously broken by scalars (the flavons), which acquire vacuum expectation values. Typically, the resulting leading order predictions for the oscillation parameters require corrections in order to comply with neutrino oscillation data. We introduce such corrections through an explicit small breaking of the symmetry. This has the advantage of solving the cosmological problems of these models without resorting to inflation. The explicit breaking induces an energy difference or "bias" between different vacua and drives the evolution of the domain walls, unavoidably produced after the symmetry breaking, towards their annihilation. Importantly, the wall annihilation leads to gravitational waves which may be observed in current and/or future experiments. We show that a distinctive pattern of gravitational waves with multiple overlapped peaks is generated when walls annihilate, which is within the reach of future detectors. We also show that cosmic walls from discrete flavor symmetries can be cosmologically safe for any spontaneous breaking scale between 1 and $10^{18}$ GeV, if the bias is chosen adequately, without the need to inflate the walls away. We use as an example a particular $A_4$ model in which an explicit breaking is included in right-handed neutrino mass terms.


title: Bayesian model selection on Scalar $\epsilon$-Field Dark Energy

abstract: The main aim of this paper is to analyse minimally-coupled scalar-fields -- quintessence and phantom -- as the main candidates to explain the accelerated expansion of the universe and compare its observables to current cosmological observations; as a byproduct we present its python module. This work includes a parameter $\epsilon$ which allows to incorporate both quintessence and phantom fields within the same analysis. Examples of the potentials, so far included, are $V(\phi)=V_0\phi^{\mu}e^{\beta \phi^\alpha}$ and $V(\phi)=V_0(\cosh(\alpha \phi)+\beta)$ with $\alpha$, $\mu$ and $\beta$ being free parameters, but the analysis can be easily extended to any other scalar field potential. Additional to the field component and the standard content of matter, the study also incorporates the contribution from spatial curvature ($\Omega_k$), as it has been the focus in recent studies. The analysis contains the most up-to-date datasets along with a nested sampler to produce posterior distributions along with the Bayesian evidence, that allows to perform a model selection. In this work we constrain the parameter-space describing the two generic potentials, and among several combinations, we found that the best-fit to current datasets is given by a model slightly favouring the quintessence field with potential $V(\phi)=V_0\phi^\mu e^{\beta \phi}$ with $\beta=0.22\pm 1.56$, $\mu = -0.41\pm 1.90$, and slightly negative curvature $\Omega_{k,0}=-0.0016\pm0.0018$, which presents deviations of $1.6\sigma$ from the standard $\Lambda$CDM model. Even though this potential contains three extra parameters, the Bayesian evidence $\mathcal{B}_{\Lambda, \phi} =2.0$ is unable to distinguish this model compared to the $\Lambda$CDM with curvature ($\Omega_{k,0}=0.0013\pm0.0018$). The potential that provides the minimal Bayesian evidence corresponds to $V(\phi)=V_0 \cosh(\alpha \phi)$ with $\alpha=-0.61\pm 1.36$.


title: Frequency-Dependent Material Motion Benchmarks for Radiative Transfer

abstract: We present a general solution for the radiation intensity in front of a purely absorbing slab moving toward an observer at constant speed and with a constant temperature. The solution is obtained by integrating the lab-frame radiation transport equation through the slab to the observer. We present comparisons between our benchmark and results from the Kull simulation code for an aluminum slab moving toward the observer at 2% the speed-of-light. We demonstrate that ignoring certain material motion correction terms in the transport equation can lead to 20-80% errors with the error magnitude growing as the frequency resolution is improved. Our results also indicate that our benchmark can identify potential errors in the implementation of material motion corrections.


title: Optoelectronic characterisation of twisted germanium sulfide nanowires with experimental observation of intrinsic ferroelectricity

abstract: We report the optoelectrical characterisation of Eshelby twisted Germanium sulfide (GeS) nanowires with first experimental observation of ferroelectric order at room temperature in GeS which is an otherwise centrosymmetric molecule. The chemical composition and structure of these nanowires were confirmed by various spectroscopic, microscopic and diffractive techniques. In addition, the nanowires were found to be stable over time. From the optoelectronic measurements we found that these p-type semiconducting GeS nanowires have up to two orders higher charge carrier mobility than GeS nanosheets and have sustainable, robust photo-switching property. The existence of room temperature ferroelectricity is confirmed by piezoresponse force microscopy which showed hysteresis and butterfly loop, characteristics of a ferroelectric material. Our observations reveal that the properties of twisted GeS nanowires can be harnessed in making efficient electric sensors, photodetectors, data memories and flexible electronics.


title: The Little W-Net That Could: State-of-the-Art Retinal Vessel Segmentation with Minimalistic Models

abstract: The segmentation of the retinal vasculature from eye fundus images represents one of the most fundamental tasks in retinal image analysis. Over recent years, increasingly complex approaches based on sophisticated Convolutional Neural Network architectures have been slowly pushing performance on well-established benchmark datasets. In this paper, we take a step back and analyze the real need of such complexity. Specifically, we demonstrate that a minimalistic version of a standard U-Net with several orders of magnitude less parameters, carefully trained and rigorously evaluated, closely approximates the performance of current best techniques. In addition, we propose a simple extension, dubbed W-Net, which reaches outstanding performance on several popular datasets, still using orders of magnitude less learnable weights than any previously published approach. Furthermore, we provide the most comprehensive cross-dataset performance analysis to date, involving up to 10 different databases. Our analysis demonstrates that the retinal vessel segmentation problem is far from solved when considering test images that differ substantially from the training data, and that this task represents an ideal scenario for the exploration of domain adaptation techniques. In this context, we experiment with a simple self-labeling strategy that allows us to moderately enhance cross-dataset performance, indicating that there is still much room for improvement in this area. Finally, we also test our approach on the Artery/Vein segmentation problem, where we again achieve results well-aligned with the state-of-the-art, at a fraction of the model complexity in recent literature. All the code to reproduce the results in this paper is released.


title: Progress in Solar Cycle Predictions: Sunspot Cycles 24-25 in Perspective

abstract: The dynamic activity of the Sun -- sustained by a magnetohydrodynamic dynamo mechanism working in its interior -- modulates the electromagnetic, particulate and radiative environment in space. While solar activity variations on short timescale create space weather, slow long-term modulation forms the basis of space climate. Space weather impacts diverse space-reliant technologies while space climate influences planetary atmospheres and climate. Having prior knowledge of the Sun's activity is important in these contexts. However, forecasting solar-stellar magnetic activity has remained an outstanding challenge. In this review, predictions for sunspot cycle 24 and the upcoming cycle 25 are summarized, and critically assessed. The analysis demonstrates that while predictions based on diverse techniques disagree across solar cycles 24--25, physics-based predictions for solar cycle 25 have converged and indicates a weak sunspot cycle 25. It is argued that this convergence in physics-based predictions is indicative of progress in the fundamental understanding of solar cycle predictability. Based on this understanding, resolutions to several outstanding questions related to solar cycle predictions are discussed.


title: Characterizing the Evolved Stellar Population in the Galactic Foreground I: Bolometric Magnitudes, Spatial Distribution and P-L Relations

abstract: Radio campaigns using maser stellar beacons have provided crucial information to characterize Galactic stellar populations. Currently, the Bulge Asymmetries and Dynamical Evolution (BAaDE) project is surveying infrared (IR) color-selected targets for SiO masers. This provides a sample of evolved stars that can be used to study the inner, optically obscured Galaxy using line of sight velocities and possibly VLBI proper motions. In order to use the BAaDE sample for kinematic studies, the stellar population should be characterized. In this study, the BAaDE targets have been cross-matched with IR (2MASS) and optical Gaia samples. By exploring the synergies of this cross-match together with Gaia parallaxes and extinction maps, the local ($d < 2$ kpc) AGB stars can be characterized. We have defined a \textit{BAaDE-Gaia} sample of 20,111 sources resulting from cross-matching BAaDE targets with IR and optical surveys. From this sample, a~{\local} sample of 1,812 evolved stars with accurate parallax measurements, confirmed evolved stellar evolution stage, and within 2 kpc distance around the Sun was selected, for which absolute (bolometric) magnitudes are estimated. The evolved stellar population with Gaia counterparts that are variable seems to be predominantly associated with AGB stars with moderate luminosity ($1,500^{+3,000}_{-500} \ L_\odot$) and periods between 250 and 1,250 days.


title: TB2J: a python package for computing magnetic interaction parameters

abstract: We present TB2J, a Python package for the automatic computation of magnetic interactions, including exchange and Dzyaloshinskii-Moriya interactions, between atoms of magnetic crystals from the results of density functional calculations. The program is based on the Green's function method with the local rigid spin rotation treated as a perturbation. As input,the package uses the output of either Wannier90, which is interfaced with many density functional theory packages,or of codes based on localized orbitals. A minimal user input is needed, which allows for easy integration into high-throughput workflows. The package is open source under BSD 2-Clause license, available at https://github.com/mailhexu/TB2J.


title: Numerical differentiation of noisy data: A unifying multi-objective optimization framework

abstract: Computing derivatives of noisy measurement data is ubiquitous in the physical, engineering, and biological sciences, and it is often a critical step in developing dynamic models or designing control. Unfortunately, the mathematical formulation of numerical differentiation is typically ill-posed, and researchers often resort to an \textit{ad hoc} process for choosing one of many computational methods and its parameters. In this work, we take a principled approach and propose a multi-objective optimization framework for choosing parameters that minimize a loss function to balance the faithfulness and smoothness of the derivative estimate. Our framework has three significant advantages. First, the task of selecting multiple parameters is reduced to choosing a single hyper-parameter. Second, where ground-truth data is unknown, we provide a heuristic for automatically selecting this hyper-parameter based on the power spectrum and temporal resolution of the data. Third, the optimal value of the hyper-parameter is consistent across different differentiation methods, thus our approach unifies vastly different numerical differentiation methods and facilitates unbiased comparison of their results. Finally, we provide an extensive open-source Python library \texttt{pynumdiff} to facilitate easy application to diverse datasets (https://github.com/florisvb/PyNumDiff).


title: Nonlinear AC and DC Conductivities in a Two-Subband n-GaAs/AlAs Heterostructure

abstract: The DC and AC conductivities of the n-GaAs/AlAs heterostructure with two filled size quantization levels are studied within a wide magnetic field range. The electron spectrum of such heterostructure is characterized by two subbands (symmetric $S$ and antisymmetric $AS$), separated by the band gap $\Delta_{12}=15.5$ meV. It is shown that, in the linear regime at the applied magnetic field $B >3$ T, the system exhibits oscillations corresponding to the integer quantum Hall effect. A quite complicated pattern of such oscillations is well interpreted in terms of transitions between Landau levels related to different subbands. At $B <1$ T, magneto-intersubband resistance oscillations (MISOs) are observed. An increase in the conductivity with the electric current flowing across the sample or in the intensity of the surface acoustic wave (SAW) in the regime of the integer quantum Hall effect is determined by an increase in the electron gas temperature. In the case of intersubband transitions, it is found that nonlinearity cannot be explained by heating. At the same time, the decrease in the AC conductivity with increasing SAW electric field is independent of frequency, but the corresponding behavior does not coincide with that corresponding to the dependence of the DC conductivity on the Hall voltage $E_y$.


title: Critical percolation on the kagome hypergraph

abstract: We study the percolation critical surface of the kagome lattice in which each triangle is allowed an arbitrary connectivity. Using the method of critical polynomials, we find points along this critical surface to high precision. This kagome hypergraph contains many unsolved problems as special cases, including bond percolation on the kagome and $(3,12^2)$ lattices, and site percolation on the hexagonal, or honeycomb, lattice, as well as a single point for which there is an exact solution. We are able to compute enough points along the critical surface to find a very accurate fit, essentially a Taylor series about the exact point, that allows estimations of the critical point of any system that lies on the surface to precision rivaling Monte Carlo and traditional techniques of similar accuracy. We find also that this system sheds light on some of the surprising aspects of the method of critical polynomials, such as why it is so accurate for certain problems, like the kagome and $(3,12^2)$ lattices. The bond percolation critical points of these lattices can be found to 17 and 18 digits, respectively, because they are in close proximity, in a sense that can be made quantitative, to the exact point on the critical surface. We also discuss in detail a parallel implementation of the method which we use here for a few calculations.


title: Inverse Compton Cooling in the Coronae of Simulated Black Hole Accretion Flows

abstract: We present a formulation for a local cooling function to be employed in the diffuse, hot corona region of 3D GRMHD simulations of accreting black holes. This new cooling function calculates the cooling rate due to inverse Compton scattering by considering the relevant microphysics in each cell in the corona and approximating the radiation energy density and Compton temperature there by integrating over the thermal seed photon flux from the disk surface. The method either assumes ion and electron temperatures are equal (1T), or calculates them separately (2T) using an instantaneous equilibrium approach predicated on the actual relevant rate equations (Coulomb and Compton). The method is shown to be consistent with a more detailed ray-tracing calculation where the bulk of the cooling occurs, but is substantially less costly to perform. As an example, we apply these methods to a \textsc{harm3d} simulation of a $10 M_\odot$, non-spinning black hole, accreting at nominally 1\% the Eddington value. Both 1T and 2T approaches lead to increased radiative efficiency and a larger fraction of total cooling in the corona as compared to the original target-temperature cooling function used by \textsc{harm3d}, especially in the 1T case. Time-averaged post-processing reveals that the continuum spectral observations predicted from these simulations are qualitatively similar to actual X-ray binary data, especially so for the 1T approach which yields a harder power-law component ($\Gamma = 2.25$) compared to the 2T version ($\Gamma = 2.53$)


title: Deuteron production in AuAu collisions at $\sqrt{s_{NN}} = 7-200$ GeV via pion catalysis

abstract: We study deuteron production using no-coalescence hydrodynamic + transport simulations of central AuAu collisions at $\sqrt{s_{NN}} = 7 - 200$ GeV. Deuterons are sampled thermally at the transition from hydrodynamics to transport, and interact in transport dominantly via $\pi p n \leftrightarrow \pi d$ reactions. The measured proton, Lambda, and deuteron transverse momentum spectra and yields are reproduced well for all collision energies considered. We further provide a possible explanation for the measured minimum in the energy dependence of the coalescence parameter, $B_2(\sqrt{s_{NN}})$ as well as for the difference between $B_2(d)$ for deuterons and that for anti-deuterons, $B_2(\bar{d})$.


title: Dual-Polarized All-Metallic Metagratings for Perfect Anomalous Reflection

abstract: We theoretically formulate and experimentally demonstrate the design of metagratings (MGs) composed of periodic rectangular grooves in a metallic medium, intended for perfect anomalous reflection. Using mode matching, a semianalytical scheme for analysis and synthesis of such MGs, containing multiple, arbitrarily arranged, grooves per period, is derived. Following the typical MG design approach, we use this formalism to identify the relevant Floquet-Bloch (FB) modes and conveniently formulate constraints for suppression of spurious scattering, directly tying the structure's geometrical degrees of freedom (DOFs) to the desired functionality. Solving this set of constraints, in turn, yields a detailed fabrication-ready MG design, without any full-wave optimization. Besides providing means to realize highly-efficient beam deflection with all-metallic formations, we show that the rectangular (two-dimensional) groove configuration enables \emph{simultaneous} manipulation of both transverse electric (TE) and transverse magnetic (TM) polarized fields, unavailable to date with common, printed-circuit-board-based, microwave MGs. In addition, we highlight a physical limitation on the TE-polarization performance, preventing the ability to achieve perfect anomalous reflection in any desired angle. These capabilities are verified using three MG prototypes, produced with standard computer numerical control (CNC) machines, demonstrating both single- and dual-polarized control of multiple diffraction modes. These results enable the use of MGs for a broader range of applications, where dual-polarized control is required, or all-metallic devices are preferable (e.g., spaceborne systems or at high operating frequencies).


title: Quantum Boltzmann equation for fermions: An attempt to calculate the NMR relaxation and decoherence times using quantum field theory techniques

abstract: Extracting macroscopic properties of a system from microscopic interactions has always been an interesting topic with the most diverse applications. Here, we use the quantum Boltzmann equation to investigate the density matrix evolution of a system of nucleons. Using the quantum field theory tools for constructing the density matrix operators and calculating the interactions is the main advantage of this equation. The right-hand side of this equation involves forward scattering and usual collision terms. As examples of application, we calculate the standard Bloch equations for the nucleon system in the presence of a constant and an oscillating magnetic field from the forward scattering term. We find the longitudinal and transverse (decoherence) relaxation times from the collision term by considering the nucleon-nucleon scattering.


title: Mobility of <c+a> dislocations in zirconium

abstract: Plasticity in hexagonal close-packed zirconium is mainly controlled by the glide of dislocations with 1/3<1-210> Burgers vectors. As these dislocations cannot accommodate deformation in the [0001] direction , twinning or glide of <c+a> dislocations, i.e. dislocations with 1/3<1-213> Burgers vector, have to be activated. We have performed in situ straining experiments in a transmission electron microscope to study the glide of <c+a> dislocations in two different zirconium samples, pure zirconium and Zircaloy-4, at room temperature. These experiments show that <c+a> dislocations exclusively glide in first-order pyramidal planes with cross-slip being activated. A much stronger lattice friction is opposing the glide of <c+a> dislocations when their orientation corresponds to the <a> direction defined by the intersection of their glide plane with the basal plane. This results in long dislocations straightened along <a> which glide either viscously or jerkily. This <a> direction governs the motion of segments with other orientations, whose shape is merely driven by the minimization of the line tension. The friction due to solute atoms is also discussed.


title: Neutrino telescopes and high-energy cosmic neutrinos

abstract: In this review paper, we present the main aspects of high-energy cosmic neutrino astrophysics. We begin by describing the generic expectations for cosmic neutrinos, including the effects of propagation from their sources to the detectors. Then we introduce the operating principles of current neutrino telescopes, and examine the main features (topologies) of the observable events. After a discussion of the main background processes, due to the concomitant presence of secondary particles produced in the terrestrial atmosphere by cosmic rays, we summarize the current status of the observations with astrophysical relevance that have been greatly contributed by IceCube detector. Then, we examine various interpretations of these findings, trying to assess the best candidate sources of cosmic neutrinos. We conclude with a brief perspective on how the field could evolve within a few years.


title: RAMSES: A full-stack application for detecting seizures and reducing data during continuous EEG monitoring

abstract: Objective: Continuous EEG (cEEG) monitoring is associated with lower mortality in critically ill patients, however it is underutilized due to the difficulty of manually interpreting prolonged streams of cEEG data. Here we present a novel real-time, machine learning-based alerting and monitoring system for epilepsy and seizures (RAMSES) that dramatically reduces the amount of manual EEG review. Methods: We developed a custom data reduction algorithm using a random forest, and deployed it within an online cloud-based platform which streams data and communicates interactively with caregivers via a web interface to display algorithm results. We validate RAMSES on cEEG recordings from 77 patients undergoing routine scalp ICU EEG monitoring. Results: On subjects with seizures we achieved >80% overall data reduction, while detecting a mean of 84% of seizures across all validation patients, with 19/27 patients achieving 100% seizure detection. On seizure free-patients, the majority of cEEG records, we reduced data requiring manual review by >83%. Conclusion: This study validates a platform for machine-learning assisted data reduction. Significance: This work represents a meaningful step towards improving utility and decreasing cost for cEEG monitoring We also make our high-quality annotated dataset of 77 ICU cEEG recordings public for others to validate and improve upon our methods.


title: A Visual Analytics Approach to Debugging Cooperative, Autonomous Multi-Robot Systems' Worldviews

abstract: Autonomous multi-robot systems, where a team of robots shares information to perform tasks that are beyond an individual robot's abilities, hold great promise for a number of applications, such as planetary exploration missions. Each robot in a multi-robot system that uses the shared-world coordination paradigm autonomously schedules which robot should perform a given task, and when, using its worldview--the robot's internal representation of its belief about both its own state, and other robots' states. A key problem for operators is that robots' worldviews can fall out of sync (often due to weak communication links), leading to desynchronization of the robots' scheduling decisions and inconsistent emergent behavior (e.g., tasks not performed, or performed by multiple robots). Operators face the time-consuming and difficult task of making sense of the robots' scheduling decisions, detecting de-synchronizations, and pinpointing the cause by comparing every robot's worldview. To address these challenges, we introduce MOSAIC Viewer, a visual analytics system that helps operators (i) make sense of the robots' schedules and (ii) detect and conduct a root cause analysis of the robots' desynchronized worldviews. Over a year-long partnership with roboticists at the NASA Jet Propulsion Laboratory, we conduct a formative study to identify the necessary system design requirements and a qualitative evaluation with 12 roboticists. We find that MOSAIC Viewer is faster- and easier-to-use than the users' current approaches, and it allows them to stitch low-level details to formulate a high-level understanding of the robots' schedules and detect and pinpoint the cause of the desynchronized worldviews.


title: The mixed affine quermassintegrals

abstract: In this paper, we introduce first the mixed affine quermassintegrals. The Aleksandrov-Fenchel inequality for the mixed affine quermassintegrals is established. As an application, the Minkowski, Brunn-Minkowski inequalities for the mixed affine quermassintegrals are also derived.


title: Heterogeneity and Superspreading Effect on Herd Immunity

abstract: We model and calculate the fraction of infected population necessary to reach herd immunity, taking into account the heterogeneity in infectiousness and susceptibility, as well as the correlation between those two parameters. We show that these cause the effective reproduction number to decrease more rapidly, and consequently have a drastic effect on the estimate of the necessary percentage of the population that has to contract the disease for herd immunity to be reached. We quantify the difference between the size of the infected population when the effective reproduction number decreases below 1 vs. the ultimate fraction of population that had contracted the disease. This sheds light on an important distinction between herd immunity and the end of the disease and highlights the importance of limiting the spread of the disease even if we plan to naturally reach herd immunity. We analyze the effect of various lock-down scenarios on the resulting final fraction of infected population. We discuss implications to COVID-19 and other pandemics and compare our theoretical results to population-based simulations. We consider the dependence of the disease spread on the architecture of the infectiousness graph and analyze different graph architectures and the limitations of the graph models.


title: Introduction to Medical Image Registration with DeepReg, Between Old and New

abstract: This document outlines a tutorial to get started with medical image registration using the open-source package DeepReg. The basic concepts of medical image registration are discussed, linking classical methods to newer methods using deep learning. Two iterative, classical algorithms using optimisation and one learning-based algorithm using deep learning are coded step-by-step using DeepReg utilities, all with real, open-accessible, medical data.


title: The mass-energy relation and the Doppler shift of a relativistic light source

abstract: This work considers the cospectral and arbitrary light emission of a moving source. The observed wavelengths of the emitted photons are described in term of kinematic and dynamical Doppler shifts in which the mass-energy relation plays a fundamental role. The presentation is an alternative way of emphasizing the importance of the concept of proper mass as a conserved quantity and the implications of the mass-energy relation when a source emits radiation. The physical contexts in which the source changes velocity after emission are discussed and a set of additional problems is presented.


title: On the Various Extensions of the BMS Group

abstract: The Bondi-Metzner-Sachs-van der Burg (BMS) group is the asymptotic symmetry group of radiating asymptotically flat spacetimes. It has recently received renewed interest in the context of the flat holography and the infrared structure of gravity. In this thesis, we investigate the consequences of considering extensions of the BMS group in four dimensions with superrotations. In particular, we apply the covariant phase space methods on a class of first order gauge theories that includes the Cartan formulation of general relativity and specify this analysis to gravity in asymptotically flat spacetime. Furthermore, we renormalize the symplectic structure at null infinity to obtain the generalized BMS charge algebra associated with smooth superrotations. We then study the vacuum structure of the gravitational field, which allows us to relate the so-called superboost transformations to the velocity kick/refraction memory effect. Afterward, we propose a new set of boundary conditions in asymptotically locally (A)dS spacetime that leads to a version of the BMS group in the presence of a non-vanishing cosmological constant, called the $\Lambda$-BMS asymptotic symmetry group. Using the holographic renormalization procedure and a diffeomorphism between Bondi and Fefferman-Graham gauges, we construct the phase space of $\Lambda$-BMS and show that it reduces to the one of the generalized BMS group in the flat limit.


title: The engulfment of aqueous droplets on perfectly wetting oil layers

abstract: Place a droplet of mineral oil on water and the oil will spread to cover the water surface in a thin film -- a phenomenon familiar to many, owing to the rainbow-faced puddles left behind leaking buses on rainy days. In this paper we study the everted problem: an aqueous droplet deposited onto a deep layer of silicone oil. As it is energetically favourable for the oil phase to spread to cover the droplet surface completely, the droplet is ultimately engulfed in the oil layer. We present a detailed study of engulfment dynamics, from the instant the droplet first impacts the oil surface until it finally sediments into the less dense oil. We study a broad range of droplet sizes (micrometric to millimetric) and oil kinematic viscosities ($10^2$ to $10^5$ cSt), corresponding to a viscosity-dominated parameter regime with relevance to oil spills. Our investigation primarily examines droplet engulfment dynamics over two distinct stages: a rapid earlier stage in which the droplet is almost entirely submerged, driven by capillary forces in the oil surface, and cloaked by a thin layer of oil; and a much slower later stage in which gravity pulls on the drop adhered to the oil surface, thus driving a peeling flow. This means that gravitational effects are essential to complete the engulfmet of the droplet, even for micrometric droplets. We observe the longest engulfment times for droplets of intermediate size. Experiments at fixed droplet size reveal a power law dependence of engulfment time on oil kinematic viscosity.


title: Efficient Algorithms to Mine Maximal Span-Trusses From Temporal Graphs

abstract: Over the last decade, there has been an increasing interest in temporal graphs, pushed by a growing availability of temporally-annotated network data coming from social, biological and financial networks. Despite the importance of analyzing complex temporal networks, there is a huge gap between the set of definitions, algorithms and tools available to study large static graphs and the ones available for temporal graphs. An important task in temporal graph analysis is mining dense structures, i.e., identifying high-density subgraphs together with the span in which this high density is observed. In this paper, we introduce the concept of $(k, \Delta)$-truss (span-truss) in temporal graphs, a temporal generalization of the $k$-truss, in which $k$ captures the information about the density and $\Delta$ captures the time span in which this density holds. We then propose novel and efficient algorithms to identify maximal span-trusses, namely the ones not dominated by any other span-truss neither in the order $k$ nor in the interval $\Delta$, and evaluate them on a number of public available datasets.


title: Stellar scattering and the formation of exponential discs in self-gravitating systems

abstract: We show, using the N-body code GADGET-2, that stellar scattering by massive clumps can produce exponential discs, and the effectiveness of the process depends on the mass of scattering centres, as well as the stability of the galactic disc. Heavy, dense scattering centres in a less stable disc generate an exponential profile quickly, with a timescale shorter than 1 Gyr. The profile evolution due to scattering can make a near-exponential disc under various initial stellar distributions. This result supports analytic theories that predict the scattering processes always favour the zero entropy gradient solution to the Jeans/Poisson equations, whose profile is a near-exponential. Profile changes are accompanied by disc thickening, and a power-law increase in stellar velocity dispersion in both vertical and radial directions is also observed through the evolution. Close encounters between stars and clumps can produce abrupt changes in stellar orbits and shift stars radially. These events can make trajectories more eccentric, but many leave eccentricities little changed. On average, orbital eccentricities of stars increase moderately with time.


title: Sparse Sensing and Optimal Precision: Robust $\mathcal{H}_{\infty}$ Optimal Observer Design with Model Uncertainty

abstract: We present a framework which incorporates three aspects of the estimation problem, namely, sparse sensor configuration, optimal precision, and robustness in the presence of model uncertainty. The problem is formulated in the $\mathcal{H}_{\infty}$ optimal observer design framework. We consider two types of uncertainties in the system, i.e. structured affine and unstructured uncertainties. The objective is to design an observer with a given $\mathcal{H}_{\infty}$ performance index with minimal number of sensors and minimal precision values, while guaranteeing the performance for all admissible uncertainties. The problem is posed as a convex optimization problem subject to linear matrix inequalities. Numerical simulations demonstrate the application of the theoretical results presented in this work.


title: Network Coding-Based Post-Quantum Cryptography

abstract: We propose a novel hybrid universal network-coding cryptosystem (HUNCC) to obtain secure post-quantum cryptography at high communication rates. The secure network-coding scheme we offer is hybrid in the sense that it combines information-theory security with public-key cryptography. In addition, the scheme is general and can be applied to any communication network, and to any public-key cryptosystem. Our hybrid scheme is based on the information theoretic notion of individual secrecy, which traditionally relies on the assumption that an eavesdropper can only observe a subset of the communication links between the trusted parties - an assumption that is often challenging to enforce. For this setting, several code constructions have been developed, where the messages are linearly mixed before transmission over each of the paths in a way that guarantees that an adversary which observes only a subset has sufficient uncertainty about each individual message. Instead, in this paper, we take a computational viewpoint, and construct a coding scheme in which an arbitrary secure cryptosystem is utilized on a subset of the links, while a pre-processing similar to the one in individual security is utilized. Under this scheme, we demonstrate 1) a computational security guarantee for an adversary which observes the entirety of the links 2) an information theoretic security guarantee for an adversary which observes a subset of the links, and 3) information rates which approach the capacity of the network and greatly improve upon the current solutions. A perhaps surprising consequence of our scheme is that, to guarantee a computational security level b, it is sufficient to encrypt a single link using a computational post-quantum scheme. In addition, the information rate approaches 1 as the number of communication links increases.


title: Unified Picture of Electromagnetic and Gravitational Forces in Two-Spinor Language

abstract: The well known Geodesic Equation of General Relativity is newly formulated in Weyl two-spinor language in a convenient way susceptible of being combined with a set of two-spinor equations, equivalent to the Lorentz Force of Electrodynamics, obtained in previous studies. General spinor equations of motion for a test charged spin 1/2 particle in arbitrary electromagnetic and gravitational fields fields are obtained describing not only spacetime trajectories but also spin precession. For the case of the Schwarzschild metric and radial trajectories, spinor equations unifying gravitational and electric forces are obtained and solved for a weak gravitational field (or long distance from the source).


title: How Hummingbirds Hum: Oscillating Aerodynamic Forces Explain Timbre of the Humming Sound

abstract: The source of the hummingbirds distinctive hum is not well understood, but there are clues to its origin in the acoustic nearfield and farfield. To unravel this mystery, we recorded the acoustic nearfield generated by six freely hovering Annas hummingbirds using a 2176 microphone array. We also directly measured the 3D aerodynamic forces generated by the hummingbird in vivo using a new aerodynamic force platform. To determine the degree to which the aerodynamic forces cause the hum, we developed a simple first-principles model to predict the acoustic field radiated by the 3D oscillating forces. The correspondence between the predicted and measured acoustic field shows the primary acoustic sources of the hum are the lift and drag forces that oscillate as the flapping wings move back and forth. The model also shows how the aerodynamic force profile of the flapping wing determines the hums timbre and sound pressure level: profiles that support bodyweight with higher harmonics radiate more acoustic power. Further, extending this model across birds and flying insects, we show how the radiated acoustic power scales with body mass to the second power - with allometric deviation making larger birds quieter and elongated flies louder. The model's ability to predict the acoustic signature of flapping wings suggests it can be applied for such diverse uses as differentiating wing and feather sounds in bird display behaviors to interpreting how insects generate courtship songs with their wings to making flapping robot wings more silent.


title: Detection of AI-Synthesized Speech Using Cepstral & Bispectral Statistics

abstract: Digital technology has made possible unimaginable applications come true. It seems exciting to have a handful of tools for easy editing and manipulation, but it raises alarming concerns that can propagate as speech clones, duplicates, or maybe deep fakes. Validating the authenticity of a speech is one of the primary problems of digital audio forensics. We propose an approach to distinguish human speech from AI synthesized speech exploiting the Bi-spectral and Cepstral analysis. Higher-order statistics have less correlation for human speech in comparison to a synthesized speech. Also, Cepstral analysis revealed a durable power component in human speech that is missing for a synthesized speech. We integrate both these analyses and propose a machine learning model to detect AI synthesized speech.


title: Gas-phase Metallicity as a Diagnostic of the Drivers of Star-formation on Different Spatial Scales

abstract: We examine the correlations of star formation rate (SFR) and gas-phase metallicity $Z$. We first predict how the SFR, cold gas mass and $Z$ will change with variations in inflow rate or in star-formation efficiency (SFE) in a simple gas-regulator framework. The changes $\Delta {\rm log}$SFR and $\Delta {\rm log} Z$, are found to be negatively (positively) correlated when driving the gas-regulator with time-varying inflow rate (SFE). We then study the correlation of $\Delta {\rm log}$sSFR (specific SFR) and $\Delta {\rm log}$(O/H) from observations, at both $\sim$100 pc and galactic scales, based on two 2-dimensional spectroscopic surveys with different spatial resolutions, MAD and MaNGA. After taking out the overall mass and radial dependences, which may reflect changes in inflow gas metallicity and/or outflow mass-loading, we find that $\Delta {\rm log}$sSFR and $\Delta {\rm log}$(O/H) on galactic are found to be negatively correlated, but $\Delta {\rm log}$sSFR and $\Delta {\rm log}$(O/H) are positively correlated on $\sim$100 pc scales within galaxies. If we assume that the variations across the population reflect temporal variations in individual objects, we conclude that variations in the star formation rate are primarily driven by time-varying inflow at galactic scales, and driven by time-varying SFE at $\sim$100 pc scales. We build a theoretical framework to understand the correlation between SFR, gas mass and metallicity, as well as their variability, which potentially uncovers the relevant physical processes of star formation at different scales.


title: Optimal Control of Convection-Cooling and Numerical Implementation

abstract: This paper is concerned with the problem of enhancing convection-cooling via active control of the incompressible velocity field, described by a stationary diffusion-convection model. This essentially leads to a bilinear optimal control problem. A rigorous proof of the existence of an optimal control is presented and the first order optimality conditions are derived for solving the control using a variational inequality. Moreover, the second order sufficient conditions are established to characterize the local minimizer. Finally, numerical experiments are conducted utilizing finite elements methods together with nonlinear iterative schemes, to demonstrate and validate the effectiveness of our control design.


title: The Information Catastrophe

abstract: Currently we produce 10 to power 21 digital bits of information annually on Earth. Assuming 20 percent annual growth rate, we estimate that 350 years from now, the number of bits produced will exceed the number of all atoms on Earth, or 10 to power 50. After 250 years, the power required to sustain this digital production will exceed 18.5 TW, or the total planetary power consumption today, and 500 years from now the digital content will account for more than half of the Earths mass, according to the mass energy information equivalence principle. Besides the existing global challenges such as climate, environment, population, food, health, energy and security, our estimates here point to another singularity event for our planet, called the Information Catastrophe.


title: Multi-Perspective Semantic Information Retrieval

abstract: Information Retrieval (IR) is the task of obtaining pieces of data (such as documents or snippets of text) that are relevant to a particular query or need from a large repository of information. While a combination of traditional keyword- and modern BERT-based approaches have been shown to be effective in recent work, there are often nuances in identifying what information is "relevant" to a particular query, which can be difficult to properly capture using these systems. This work introduces the concept of a Multi-Perspective IR system, a novel methodology that combines multiple deep learning and traditional IR models to better predict the relevance of a query-sentence pair, along with a standardized framework for tuning this system. This work is evaluated on the BioASQ Biomedical IR + QA challenges.


title: Accurate TLS Fingerprinting using Destination Context and Knowledge Bases

abstract: Network fingerprinting is used to identify applications, provide insight into network traffic, and detect malicious activity. With the broad adoption of TLS, traditional fingerprinting techniques that rely on clear-text data are no longer viable. TLS-specific techniques have been introduced that create a fingerprint string from carefully selected data features in the client_hello to facilitate process identification before data is exchanged. Unfortunately, this approach fails in practice because hundreds of processes can map to the same fingerprint string. We solve this problem by presenting a TLS fingerprinting system that makes use of the destination address, port, and server name in addition to a carefully constructed fingerprint string. The destination context is used to disambiguate the set of processes that match a fingerprint string by applying a weighted naive Bayes classifier, resulting in far greater performance.


title: COVID-19 Policy Impact Evaluation: A guide to common design issues

abstract: Policy responses to COVID-19, particularly those related to non-pharmaceutical interventions, are unprecedented in scale and scope. Epidemiologists are more involved in policy decisions and evidence generation than ever before. However, policy impact evaluations always require a complex combination of circumstance, study design, data, statistics, and analysis. Beyond the issues that are faced for any policy, evaluation of COVID-19 policies is complicated by additional challenges related to infectious disease dynamics and lags, lack of direct observation of key outcomes, and a multiplicity of interventions occurring on an accelerated time scale. The methods needed for policy-level impact evaluation are not often used or taught in epidemiology, and differ in important ways that may not be obvious. The volume and speed, and methodological complications of policy evaluations can make it difficult for decision-makers and researchers to synthesize and evaluate strength of evidence in COVID-19 health policy papers. In this paper, we (1) introduce the basic suite of policy impact evaluation designs for observational data, including cross-sectional analyses, pre/post, interrupted time-series, and difference-in-differences analysis, (2) demonstrate key ways in which the requirements and assumptions underlying these designs are often violated in the context of COVID-19, and (3) provide decision-makers and reviewers a conceptual and graphical guide to identifying these key violations. The overall goal of this paper is to help epidemiologists, policy-makers, journal editors, journalists, researchers, and other research consumers understand and weigh the strengths and limitations of evidence that is essential to decision-making.


title: Dense CNN with Self-Attention for Time-Domain Speech Enhancement

abstract: Speech enhancement in the time domain is becoming increasingly popular in recent years, due to its capability to jointly enhance both the magnitude and the phase of speech. In this work, we propose a dense convolutional network (DCN) with self-attention for speech enhancement in the time domain. DCN is an encoder and decoder based architecture with skip connections. Each layer in the encoder and the decoder comprises a dense block and an attention module. Dense blocks and attention modules help in feature extraction using a combination of feature reuse, increased network depth, and maximum context aggregation. Furthermore, we reveal previously unknown problems with a loss based on the spectral magnitude of enhanced speech. To alleviate these problems, we propose a novel loss based on magnitudes of enhanced speech and a predicted noise. Even though the proposed loss is based on magnitudes only, a constraint imposed by noise prediction ensures that the loss enhances both magnitude and phase. Experimental results demonstrate that DCN trained with the proposed loss substantially outperforms other state-of-the-art approaches to causal and non-causal speech enhancement.


title: On system-wide safety staffing of large-scale parallel server networks

abstract: We introduce a "system-wide safety staffing" (SWSS) parameter for multiclass multi-pool networks of any tree topology, Markovian or non-Markovian, in the Halfin-Whitt regime. This parameter can be regarded as the optimal reallocation of the capacity fluctuations (positive or negative) of order $\sqrt{n}$ when each server pool employs a square-root staffing rule. We provide an explicit form of the SWSS as a function of the system parameters, which is derived using a graph theoretic approach based on Gaussian elimination. For Markovian networks, we give an equivalent characterization of the SWSS parameter via the drift parameters of the limiting diffusion. We show that if the SWSS parameter is negative, the limiting diffusion and the diffusion-scaled queueing processes are transient under any Markov control, and cannot have a stationary distribution when this parameter is zero. If it is positive, we show that the diffusion-scaled queueing processes are uniformly stabilizable, that is, there exists a scheduling policy under which the stationary distributions of the controlled processes are tight over the size of the network. In addition, there exists a control under which the limiting controlled diffusion is exponentially ergodic. Thus we have identified a necessary and sufficient condition for the uniform stabilizability of such networks in the Halfin-Whitt regime. We use a constant control resulting from the leaf elimination algorithm to stabilize the limiting controlled diffusion, while a family of Markov scheduling policies which are easy to compute are used to stabilize the diffusion-scaled processes. Finally, we show that under these controls the processes are exponentially ergodic and the stationary distributions have exponential tails.


title: Measurements of multiplicity fluctuations of identified hadrons in inelastic proton-proton interactions at the CERN Super Proton Synchrotron

abstract: Measurements of multiplicity fluctuations of identified hadrons produced in inelastic p+p interactions at 31, 40, 80, and 158~\GeVc beam momentum are presented. Three different measures of multiplicity fluctuations are used: the scaled variance $\omega$ and strongly intensive measures $\Sigma$ and $\Delta$. These fluctuation measures involve second and first moments of joint multiplicity distributions. Data analysis is performed using the Identity method which corrects for incomplete particle identification. Strongly intensive quantities are calculated in order to allow for a direct comparison to corresponding results on nucleus-nucleus collisions. The results for different hadron types are shown as a function of collision energy. A comparison with predictions of string-resonance Monte-Carlo models: Epos, Smash and Venus, is also presented.


title: Coherent Manipulation of the Internal State of Ultracold $^{87}$Rb$^{133}$Cs Molecules with Multiple Microwave Fields

abstract: We explore coherent multi-photon processes in $^{87}$Rb$^{133}$Cs molecules using 3-level lambda and ladder configurations of rotational and hyperfine states, and discuss their relevance to future applications in quantum computation and quantum simulation. In the lambda configuration, we demonstrate the driving of population between two hyperfine levels of the rotational ground state via a two-photon Raman transition. Such pairs of states may be used in the future as a quantum memory, and we measure a Ramsey coherence time for a superposition of these states of 58(9) ms. In the ladder configuration, we show that we can generate and coherently populate microwave dressed states via the observation of an Autler-Townes doublet. We demonstrate that we can control the strength of this dressing by varying the intensity of the microwave coupling field. Finally, we perform spectroscopy of the rotational states of $^{87}$Rb$^{133}$Cs up to $N=6$, highlighting the potential of ultracold molecules for quantum simulation in synthetic dimensions. By fitting the measured transition frequencies we determine a new value of the centrifugal distortion coefficient $D_v=h\times207.3(2)~$Hz.


title: Three Years of HARPS-N High-Resolution Spectroscopy and Precise Radial Velocity Data for the Sun

abstract: The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such high-cadence, long-baseline data set is crucial for understanding spurious radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub-\ms\,systematics that needed to be corrected for. The goal of this manuscript is to i) present a new data reduction software for HARPS-N, ii) demonstrate the improvement brought by this new software on the first three years of the HARPS-N solar data set, and iii) release all the obtained solar products, from extracted spectra to precise radial velocities. To correct for the instrumental systematics observed in the data reduced with the current version of the HARPS-N data reduction software (DRS version 3.7), we adapted the newly available ESPRESSO DRS (version 2.2.3) to HARPS-N and developed new optimized recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS. The most significant improvement brought by the new DRS is a strong decrease in the day-to-day radial-velocity scatter, from 1.27 to 1.07\ms; this is thanks to a more robust method to derive wavelength solutions, but also to the use of calibrations closer in time. The newly derived solar radial-velocities are also better correlated with the chromospheric activity level of the Sun on the long-term, with a Pearson correlation coefficient of 0.93 compared to 0.77 before, which is expected from our understanding of stellar signals. Finally, we also discuss how HARPS-N spectral ghosts contaminate the measurement of the calcium activity index, and present an efficient technique to derive an index free of instrumental systematics. This paper presents a new data reduction software for HARPS-N, and demonstrates its improvements [...]


title: Triangle conics and cubics

abstract: This is a paper about triangle cubics and conics in classical geometry with elements of projective geometry. In recent years, N.J. Wildberger has actively dealt with this topic using an algebraic perspective. Triangle conics were also studied in detail by H.M. Cundy and C.F. Parry recently. The main task of the article was to develop an algorithm for creating curves, which pass through triangle centers. During the research, it was noticed that some different triangle centers in distinct triangles coincide. The simplest example: an incenter in a base triangle is an orthocenter in an excentral triangle. This was the key for creating an algorithm. Indeed, we can match points belonging to one curve (base curve) with other points of another triangle. Therefore, we get a new intersting geometrical object. During the research were derived number of new triangle conics and cubics, were considered their properties in Euclidian space. In addition, was discussed corollaries of the obtained theorems in projective geometry, what proves that all of the descovered results could be transfered to the projeticve plane.


title: Practical and Parallelizable Algorithms for Non-Monotone Submodular Maximization with Size Constraint

abstract: We present combinatorial and parallelizable algorithms for maximization of a submodular function, not necessarily monotone, with respect to a size constraint. We improve the best approximation factor achieved by an algorithm that has optimal adaptivity and nearly optimal query complexity to $0.193 - \varepsilon$. The conference version of this work mistakenly employed a subroutine that does not work for non-monotone, submodular functions. In this version, we propose a fixed and improved subroutine to add a set with high average marginal gain, ThreshSeq, which returns a solution in $O( \log(n) )$ adaptive rounds with high probability. Moreover, we provide two approximation algorithms. The first has approximation ratio $1/6 - \varepsilon$, adaptivity $O( \log (n) )$, and query complexity $O( n \log (k) )$, while the second has approximation ratio $0.193 - \varepsilon$, adaptivity $O( \log^2 (n) )$, and query complexity $O(n \log (k))$. Our algorithms are empirically validated to use a low number of adaptive rounds and total queries while obtaining solutions with high objective value in comparison with state-of-the-art approximation algorithms, including continuous algorithms that use the multilinear extension.


title: On the intrinsic torsion of spacetime structures

abstract: We briefly review the notion of the intrinsic torsion of a $G$-structure and then go on to classify the intrinsic torsion of the $G$-structures associated with spacetimes: namely, galilean (or Newton-Cartan), carrollian, aristotelian and bargmannian. In the case of galilean structures, the intrinsic torsion classification agrees with the well-known classification into torsionless, twistless torsional and torsional Newton-Cartan geometries. In the case of carrollian structures, we find that intrinsic torsion allows us to classify Carroll manifolds into four classes, depending on the action of the Carroll vector field on the spatial metric, or equivalently in terms of the nature of the null hypersurfaces of a lorentzian manifold into which a carrollian geometry may embed. By a small refinement of the results for galilean and carrollian structures, we show that there are sixteen classes of aristotelian structures, which we characterise geometrically. Finally, the bulk of the paper is devoted to the case of bargmannian structures, where we find twenty-seven classes which we also characterise geometrically while simultaneously relating some of them to the galilean and carrollian structures. This paper is dedicated to Dmitri Vladimirovich Alekseevsky on his 80th birthday.


title: Strong-Field Terahertz Control of Plasmon Induced Opacity in Photoexcited Metamaterial

abstract: A terahertz metamaterial consisting of radiative slot antennas and subradiant complementary split-ring resonators exhibits plasmon induced opacity in a narrow spectral range due to the destructive interference between the bright and dark modes of the coupled oscillators. Femtosecond optical excitations instantly quench the mode coupling and plasmon oscillations, injecting photocarriers into the metamaterial. The plasmon resonances in the coupled metamaterial are transiently restored by intense terahertz pulses. The strong terahertz fields induce intervalley scattering and interband tunneling of the photocarries, and achieve significant reduction of the photocarrier mobility. The ultrafast dynamics of the nonlinear THz interactions reveals intricate interplay between photocarriers and plasmon oscillations. The high-field THz control of the plasmon oscillations implies potential applications to ultrahigh-speed plasmonics.


title: A DFT Investigation of the Electronic and Optical Properties of Pentadiamond

abstract: Recently, a new carbon 3D carbon allotrope named pentadiamond was proposed. Pentadiamond is composed of carbon atoms in mixed sp$^2$ and sp$^3$-like hybridization. In this work, we have carried out a detailed investigation of the electronic and optical properties of pentadiamond structure using first-principles (DFT) methods. Our results show that pentadiamond has an indirect bandgap semiconductor of $2.50$ eV with GGA-PBE and $3.31$ eV with HSE06. Its static dielectric constant is $4.70$ and the static refractive index is $2.16$. Pentadiamond presents low reflectivity, almost 40$\%$, for all-optical spectrum, making it a good structure to be used as a UV collector. Also, pentadiamond exhibits optical activity in the UV range where other carbon allotropes, such as diamond and 8-tetra(2,2) tubulane show no activity.


title: Zero products of Toeplitz operators on Reinhardt domains

abstract: Let $\Omega$ be a bounded Reinhardt domain in $\mathbb{C}^n$ and $\phi_1,\ldots,\phi_m$ be finite sums of bounded quasi-homogeneous functions. We show that if the product of Toeplitz operators $T_{\phi_m}\cdots T_{\phi_1}=0$ on the Bergman space on $\Omega$, then $\phi_j=0$ for some $j$.


title: Augmenting machine learning photometric redshifts with Gaussian mixture models

abstract: Wide-area imaging surveys are one of the key ways of advancing our understanding of cosmology, galaxy formation physics, and the large-scale structure of the Universe in the coming years. These surveys typically require calculating redshifts for huge numbers (hundreds of millions to billions) of galaxies - almost all of which must be derived from photometry rather than spectroscopy. In this paper we investigate how using statistical models to understand the populations that make up the colour-magnitude distribution of galaxies can be combined with machine learning photometric redshift codes to improve redshift estimates. In particular we combine the use of Gaussian Mixture Models with the high performing machine learning photo-z algorithm GPz and show that modelling and accounting for the different colour-magnitude distributions of training and test data separately can give improved redshift estimates, reduce the bias on estimates by up to a half, and speed up the run-time of the algorithm. These methods are illustrated using data from deep optical and near infrared data in two separate deep fields, where training and test data of different colour-magnitude distributions are constructed from the galaxies with known spectroscopic redshifts, derived from several heterogeneous surveys.


title: Why should I not follow you? Reasons For and Reasons Against in Responsible Recommender Systems

abstract: A few Recommender Systems (RS) resort to explanations so as to enhance trust in recommendations. However, current techniques for explanation generation tend to strongly uphold the recommended products instead of presenting both reasons for and reasons against them. We argue that an RS can better enhance overall trust and transparency by frankly displaying both kinds of reasons to users.We have developed such an RS by exploiting knowledge graphs and by applying Snedegar's theory of practical reasoning. We show that our implemented RS has excellent performance and we report on an experiment with human subjects that shows the value of presenting both reasons for and against, with significant improvements in trust, engagement, and persuasion.


title: Analysis on quasidisks; a unified approach through transmission and jump problems

abstract: We give an exposition of results from a crossroad between geometric function theory, harmonic analysis, boundary value problems and approximation theory, which characterize quasicircles. We will specifically expose the interplay between the jump decomposition, singular integral operators and approximation by Faber series. Our unified point of view is made possible by the the concept of transmission.


title: Development of an Extrudable Paste to build Mycelium-bound Composites

abstract: Mycelium-bound composites are promising materials for sustainable packaging, insulation, fashion, and architecture. However, moulding is the main fabrication process explored to date, strongly limiting the ability to design the complex shapes that could widen the range of applications. Extrusion is a facile and low energy-cost process that has not been explored yet for mycelium-bound composites with design freedom and structural properties. In this study, we combine cheap, easily and commonly available agricultural waste materials, bamboo microfibres, chitosan, and mycelium from Ganoderma Lucidum, to establish a composite mixture that is workable, extrudable and buildable. We study the impact of bamboo fibre size, chitosan concentration, pH and weight ratio of bamboo to chitosan to determine the optimum growth condition for the mycelium as well as highest mechanical stiffness. The resulting materials have thus low energy costs, are sustainable and can be shaped in diverse forms easily. The developed composition is promising to further explore the use of mycelium-bound materials for structural applications using agricultural waste.


title: Compression-aware Continual Learning using Singular Value Decomposition

abstract: We propose a compression based continual task learning method that can dynamically grow a neural network. Inspired from the recent model compression techniques, we employ compression-aware training and perform low-rank weight approximations using singular value decomposition (SVD) to achieve network compaction. By encouraging the network to learn low-rank weight filters, our method achieves compressed representations with minimal performance degradation without the need for costly fine-tuning. Specifically, we decompose the weight filters using SVD and train the network on incremental tasks in its factorized form. Such a factorization allows us to directly impose sparsity-inducing regularizers over the singular values and allows us to use fewer number of parameters for each task. We further introduce a novel shared representational space based learning between tasks. This promotes the incoming tasks to only learn residual task-specific information on top of the previously learnt weight filters and greatly helps in learning under fixed capacity constraints. Our method significantly outperforms prior continual learning approaches on three benchmark datasets, demonstrating accuracy improvements of 10.3%, 12.3%, 15.6% on 20-split CIFAR-100, miniImageNet and a 5-sequence dataset, respectively, over state-of-the-art. Further, our method yields compressed models that have ~3.64x, 2.88x, 5.91x fewer number of parameters respectively, on the above mentioned datasets in comparison to baseline individual task models. Our source code is available at https://github.com/pavanteja295/CACL.


title: Interpolation in Model Spaces

abstract: In this paper we consider interpolation in model spaces, $H^2 \ominus B H^2$ with $B$ a Blaschke product. We study unions of interpolating sequences for two sequences that are far from each other in the pseudohyperbolic metric as well as two sequences that are close to each other in the pseudohyperbolic metric. The paper concludes with a discussion of the behavior of Frostman sequences under perturbations.


title: The relative emission from chromospheres and coronae: dependence on spectral type and age

abstract: Extreme-ultraviolet and X-ray emissions from stellar coronae drive mass loss from exoplanet atmospheres, and ultraviolet emission from stellar chromospheres drives photo-chemistry in exoplanet atmospheres. Comparisons of the spectral energy distributions of host stars are, therefore, essential for understanding the evolution and habitability of exoplanets. The large number of stars observed with the MUSCLES, Mega-MUSCLES, and other recent HST observing programs has provided for the first time a large sample (79 stars) of reconstructed Lyman-alpha fluxes that we compare with X-ray fluxes to identify significant patterns in the relative emission from these two atmospheric regions as a function of stellar age and effective temperature. We find that as stars age on the main sequence, the emissions from their chromospheres and coronae follow a pattern in response to the amount of magnetic heating in these atmospheric layers. A single trendline slope describes the pattern of X-ray vs. Lyman-alpha emission for G and K dwarfs, but the different trendlines for M dwarf stars show that the Lyman-alpha fluxes of M stars are significantly smaller than warmer stars with the same X-ray flux. The X-ray and Lyman-alpha luminosities divided by the stellar bolometric luminosities show different patterns depending on stellar age. The L(Lyman-alpha)/L(bol) ratios increase smoothly to cooler stars of all ages, but the L(X)/L(bol) ratios show different trends. For older stars, the increase in coronal emission with decreasing T(eff) is much steeper than chromospheric emission. We suggest a fundamental link between atmospheric properties and trendlines relating coronal and chromospheric heating,


title: CoNCRA: A Convolutional Neural Network Code Retrieval Approach

abstract: Software developers routinely search for code using general-purpose search engines. However, these search engines cannot find code semantically unless it has an accompanying description. We propose a technique for semantic code search: A Convolutional Neural Network approach to code retrieval (CoNCRA). Our technique aims to find the code snippet that most closely matches the developer's intent, expressed in natural language. We evaluated our approach's efficacy on a dataset composed of questions and code snippets collected from Stack Overflow. Our preliminary results showed that our technique, which prioritizes local interactions (words nearby), improved the state-of-the-art (SOTA) by 5% on average, retrieving the most relevant code snippets in the top 3 (three) positions by almost 80% of the time. Therefore, our technique is promising and can improve the efficacy of semantic code retrieval.


title: Multimodal Autonomous Last Mile Delivery System Design and Application

abstract: With the rapid increase in congestion, alternative solutions are needed to efficiently use the capacity of our existing networks. This paper focuses on exploring the emerging autonomous technologies for on-demand food delivery in congested urban cities. Three different last mile food delivery systems are proposed in this study employing aerial and ground autonomous vehicles technologies. The three proposed systems are: robot delivery system, drone delivery system and a hybrid delivery system. In the hybrid system the concept of hub-and-spoke network is explored in order to consolidate orders and reach more destinations in less time. To investigate the performance of the three proposed delivery systems, they are applied to the city of Mississauga network, in an in-house agent-based simulation in MATLAB. 18 Scenarios are tested differing in terms of demand and fleet size. The results show that the hybrid robot-drone delivery system performs the best with a fleet side of 25 robots and 15 drones and with an average preparation and delivery time less than the individual robot and drone system by 48% and 42% respectively.


title: Augmented Gaussian Random Field: Theory and Computation

abstract: We propose the novel augmented Gaussian random field (AGRF), which is a universal framework incorporating the data of observable and derivatives of any order. Rigorous theory is established. We prove that under certain conditions, the observable and its derivatives of any order are governed by a single Gaussian random field, which is the aforementioned AGRF. As a corollary, the statement ``the derivative of a Gaussian process remains a Gaussian process'' is validated, since the derivative is represented by a part of the AGRF. Moreover, a computational method corresponding to the universal AGRF framework is constructed. Both noiseless and noisy scenarios are considered. Formulas of the posterior distributions are deduced in a nice closed form. A significant advantage of our computational method is that the universal AGRF framework provides a natural way to incorporate arbitrary order derivatives and deal with missing data. We use four numerical examples to demonstrate the effectiveness of the computational method. The numerical examples are composite function, damped harmonic oscillator, Korteweg-De Vries equation, and Burgers' equation.


title: Uniformization and Constructive Analytic Continuation of Taylor Series

abstract: We analyze the problem of global reconstruction of functions as accurately as possible, based on partial information in the form of a truncated power series at some point, and additional analyticity properties. This situation occurs frequently in applications. The question of the optimal procedure was open, and we formulate it as a well-posed mathematical problem. Its solution leads to a practical method which provides dramatic accuracy improvements over existing techniques. Our procedure is based on uniformization of Riemann surfaces. As an application, we show that our procedure can be implemented for solutions of a wide class of nonlinear ODEs. We find a new uniformization method, which we use to construct the uniformizing maps needed for special functions, including solution of the Painlev'e equations P_I-P_V. We also introduce a new rigorous and constructive method of regularization, elimination of singularities whose position and type are known. If these are unknown, the same procedure enables a highly sensitive resonance method to determine the position and type of a singularity. In applications where less explicit information is available about the Riemann surface, our approach and techniques lead to new approximate, but still much more precise reconstruction methods than existing ones, especially in the vicinity of singularities, which are the points of greatest interest.