title: HDR-cGAN: Single LDR to HDR Image Translation using Conditional GAN

abstract: The prime goal of digital imaging techniques is to reproduce the realistic appearance of a scene. Low Dynamic Range (LDR) cameras are incapable of representing the wide dynamic range of the real-world scene. The captured images turn out to be either too dark (underexposed) or too bright (overexposed). Specifically, saturation in overexposed regions makes the task of reconstructing a High Dynamic Range (HDR) image from single LDR image challenging. In this paper, we propose a deep learning based approach to recover details in the saturated areas while reconstructing the HDR image. We formulate this problem as an image-to-image (I2I) translation task. To this end, we present a novel conditional GAN (cGAN) based framework trained in an end-to-end fashion over the HDR-REAL and HDR-SYNTH datasets. Our framework uses an overexposed mask obtained from a pre-trained segmentation model to facilitate the hallucination task of adding details in the saturated regions. We demonstrate the effectiveness of the proposed method by performing an extensive quantitative and qualitative comparison with several state-of-the-art single-image HDR reconstruction techniques.


title: Rerunning OCR: A Machine Learning Approach to Quality Assessment and Enhancement Prediction

abstract: Iterating with new and improved OCR solutions enforces decision making when it comes to targeting the right candidates for reprocessing. This especially applies when the underlying data collection is of considerable size and rather diverse in terms of fonts, languages, periods of publication and consequently OCR quality. This article captures the efforts of the National Library of Luxembourg to support those targeting decisions. They are crucial in order to guarantee low computational overhead and reduced quality degradation risks, combined with a more quantifiable OCR improvement. In particular, this work explains the methodology of the library with respect to text block level quality assessment. Through extension of this technique, a regression model, that is able to take into account the enhancement potential of a new OCR engine, is also presented. They both mark promising approaches, especially for cultural institutions dealing with historical data of lower quality.


title: Theory Motivation: What measurements are needed?

abstract: I give a theory motivation for future measurements in quark flavour physics, trying to identify observables, which are less familiar, but nevertheless interesting and promising.


title: Global Convergence and Stability of Stochastic Gradient Descent

abstract: In machine learning, stochastic gradient descent (SGD) is widely deployed to train models using highly non-convex objectives with equally complex noise models. Unfortunately, SGD theory often makes restrictive assumptions that fail to capture the non-convexity of real problems, and almost entirely ignore the complex noise models that exist in practice. In this work, we make substantial progress on this shortcoming. First, we establish that SGD's iterates will either globally converge to a stationary point or diverge under nearly arbitrary nonconvexity and noise models. Under a slightly more restrictive assumption on the joint behavior of the non-convexity and noise model that generalizes current assumptions in the literature, we show that the objective function cannot diverge, even if the iterates diverge. As a consequence of our results, SGD can be applied to a greater range of stochastic optimization problems with confidence about its global convergence behavior and stability.


title: Estimating Potential Outcome Distributions with Collaborating Causal Networks

abstract: Traditional causal inference approaches leverage observational study data to estimate the difference in observed and unobserved outcomes for a potential treatment, known as the Conditional Average Treatment Effect (CATE). However, CATE corresponds to the comparison on the first moment alone, and as such may be insufficient in reflecting the full picture of treatment effects. As an alternative, estimating the full potential outcome distributions could provide greater insights. However, existing methods for estimating treatment effect potential outcome distributions often impose restrictive or simplistic assumptions about these distributions. Here, we propose Collaborating Causal Networks (CCN), a novel methodology which goes beyond the estimation of CATE alone by learning the full potential outcome distributions. Estimation of outcome distributions via the CCN framework does not require restrictive assumptions of the underlying data generating process. Additionally, CCN facilitates estimation of the utility of each possible treatment and permits individual-specific variation through utility functions. CCN not only extends outcome estimation beyond traditional risk difference, but also enables a more comprehensive decision-making process through definition of flexible comparisons. Under assumptions commonly made in the causal literature, we show that CCN learns distributions that asymptotically capture the true potential outcome distributions. Furthermore, we propose an adjustment approach that is empirically effective in alleviating sample imbalance between treatment groups in observational data. Finally, we evaluate the performance of CCN in multiple synthetic and semi-synthetic experiments. We demonstrate that CCN learns improved distribution estimates compared to existing Bayesian and deep generative methods as well as improved decisions with respects to a variety of utility functions.


title: Adinkras: Graphs of Clifford Algebra Representations, Supersymmetry, and Codes

abstract: An Adinkra is a graph from the study of supersymmetry in particle physics, but it can be adapted to study Clifford algebra representations. The graph in this context is called a Cliffordinkra, and puts some standard ideas in Clifford algebra representations in a geometric and visual context. In the past few years there have been developments in Adinkras that have shown how they are connected to error correcting codes, algebraic topology, algebraic geometry, and combinatorics. These connections also arise for Cliffordinkras. This paper introduces Cliffordinkras and describes the relationship to these subjects in that context. No previous knowledge of Adinkras and supersymmetry is assumed.


title: Magnetized Ablative Rayleigh-Taylor Instability in 3-D

abstract: 3-D extended-MHD simulations of the magnetized ablative Rayleigh-Taylor instability are presented for the first time. Previous 2-D simulations claiming perturbation suppression by magnetic tension are shown to be misleading, as they do not include the most unstable dimension. For perturbation modes along the applied field direction, the magnetic field simultaneously reduces ablative stabilization and adds magnetic tension stabilization; the stabilizing term is found to dominate for applied fields >5T, with both effects increasing in importance at short wavelengths. For modes perpendicular to the applied field, magnetic tension does not directly stabilize the perturbation, but can result in moderately slower growth due to the perturbation appearing 2-D (albeit in a different orientation to 2-D ICF simulations). In cases where thermal ablative stabilization is dominant the applied field increases the peak bubble-spike height. Resistive diffusion is shown to be important for short wavelengths and long timescales, reducing the effectiveness of tension stabilization.


title: Wide field beamformed observation with MeerKAT

abstract: Large-scale beamforming with radio interferometers has the potential to revolutionize the science done with pulsars and fast radio bursts by improving the survey efficiency for these sources. We describe a wide-field beamformer for the MeerKAT radio telescope and outline strategies to optimally design such surveys. A software implementation of these techniques, ${\rm M{\small OSAIC}}$ is introduced and its application in the MeerKAT telescope is presented. We show initial results using the beamformer by observing a globular cluster to track several pulsars simultaneously and demonstrate the source localization capability of this observation.


title: Learning to shortcut and shortlist order fulfillment deciding

abstract: With the increase of order fulfillment options and business objectives taken into consideration in the deciding process, order fulfillment deciding is becoming more and more complex. For example, with the advent of ship from store retailers now have many more fulfillment nodes to consider, and it is now common to take into account many and varied business goals in making fulfillment decisions. With increasing complexity, efficiency of the deciding process can become a real concern. Finding the optimal fulfillment assignments among all possible ones may be too costly to do for every order especially during peak times. In this work, we explore the possibility of exploiting regularity in the fulfillment decision process to reduce the burden on the deciding system. By using data mining we aim to find patterns in past fulfillment decisions that can be used to efficiently predict most likely assignments for future decisions. Essentially, those assignments that can be predicted with high confidence can be used to shortcut, or bypass, the expensive deciding process, or else a set of most likely assignments can be used for shortlisting -- sending a much smaller set of candidates for consideration by the fulfillment deciding system.


title: Solving realistic security-constrained optimal power flow problems

abstract: We present a decomposition approach for obtaining good feasible solutions for the security-constrained alternating-current optimal power flow (SCACOPF) problem at an industrial scale and under real-world time and computational limits. The approach aims at complementing the existing body of literature on bounding the problem via convex relaxations. It was designed for the participation in ARPA-E's Grid Optimization (GO) Competition Challenge 1. The challenge focused on a near-real-time version of the SCACOPF problem where a base case operating point is optimized taking into account possible single-element contingencies, after which the system adapts its operating point following the response of automatic frequency drop controllers and voltage regulators. Our solution approach for this problem relies on state-of-the-art nonlinear programming algorithms and employs nonconvex relaxations for complementarity constraints, a specialized two-stage decomposition technique with sparse approximations of recourse terms, and contingency ranking and pre-screening. The paper also outlines the salient features of our implementation, such as fast model functions and derivatives evaluation, warm-starting strategies, and asynchronous parallelism. We discuss the results of the independent benchmark of our approach done by ARPA-E's GO team in Challenge 1, which found that our methodology consistently produces high quality solutions across a wide range of network sizes and difficulty. Finally, we conclude by outlining potential extensions and improvements of our methodology.


title: A manifold learning approach for gesture recognition from micro-Doppler radar measurements

abstract: A recent paper (Neural Networks, {\bf 132} (2020), 253-268) introduces a straightforward and simple kernel based approximation for manifold learning that does not require the knowledge of anything about the manifold, except for its dimension. In this paper, we examine how the pointwise error in approximation using least squares optimization based on similarly localized kernels depends upon the data characteristics and deteriorates as one goes away from the training data. The theory is presented with an abstract localized kernel, which can utilize any prior knowledge about the data being located on an unknown sub-manifold of a known manifold. We demonstrate the performance of our approach using a publicly available micro-Doppler data set, and investigate the use of different preprocessing measures, kernels, and manifold dimensions. Specifically, it is shown that the localized kernel introduced in the above mentioned paper when used with PCA components leads to a near-competitive performance to deep neural networks, and offers significant improvements in training speed and memory requirements. To demonstrate the fact that our methods are agnostic to the domain knowledge, we examine the classification problem in a simple video data set.


title: Open quantum system dynamics and the mean force Gibbs state

abstract: The dynamical convergence of a system to the thermal distribution, or Gibbs state, is a standard assumption across all of the physical sciences. The Gibbs state is determined just by temperature and the system's energies alone. But at decreasing system sizes, i.e. for nanoscale and quantum systems, the interaction with their environments is not negligible. The question then arises: Is the system's steady state still the Gibbs state? And if not, how may the steady state depend on the interaction details? Here we provide an overview of recent progress on answering these questions. We expand on the state-of-the-art along two general avenues: First we take the static point-of-view which postulates the so-called mean force Gibbs state. This view is commonly adopted in the field of strong coupling thermodynamics, where modified laws of thermodynamics and non-equilibrium fluctuation relations are established on the basis of this modified state. Second, we take the dynamical point-of-view, originating from the field of open quantum systems, which examines the time-asymptotic steady state within two paradigms. We describe the mathematical paradigm which proves return to equilibrium, i.e. convergence to the mean force Gibbs state, and then discuss a number of microscopic physical methods, particularly master equations. We conclude with a summary of established links between statics and equilibration dynamics, and provide an extensive list of open problems. This comprehensive overview will be of interest to researchers in the wider fields of quantum thermodynamics, open quantum systems, mesoscopic physics, statistical physics and quantum optics, and will find applications whenever energy is exchanged on the nanoscale, from quantum chemistry and biology, to magnetism and nanoscale heat management.


title: Longest paths related to Steenrod length of real projective spaces

abstract: We answer an open problem in homotopy theory: Find a combinatorial interpretation for a function $f(n)$ related to the Steenrod length of real projective spaces.


title: Model-completeness and decidability of the additive structure of integers expanded with a function for a Beatty sequence

abstract: We introduce a model-complete theory which completely axiomatizes the structure $Z_{\alpha}=(Z, +, 0, 1, f)$ where $f : x \to \lfloor{\alpha} x \rfloor $ is a unary function with $\alpha$ a fixed transcendental number. When $\alpha$ is computable, our theory is recursively enumerable, and hence decidable as a result of completeness. Therefore, this result fits into the more general theme of adding traces of multiplication to integers without losing decidability.


title: Operads in Unstable Global Homotopy Theory

abstract: We study operads in unstable global homotopy theory, which is the homotopy theory of spaces with compatible actions by all compact Lie groups. We show that the theory of these operads works remarkably well, as for example it is possible to give a model structure for the category of algebras over any such operad. We define global $E_\infty$-operads, a good generalization of $E_\infty$-operads to the global setting, and we give a rectification result for algebras over them.


title: Gas phase Elemental abundances in Molecular cloudS (GEMS) V. Methanol in Taurus

abstract: Methanol, one of the simplest complex organic molecules in the Interstellar Medium (ISM), has been shown to be present and extended in cold environments such as starless cores. We aim at studying methanol emission across several starless cores and investigate the physical conditions at which methanol starts to be efficiently formed, as well as how the physical structure of the cores and their surrounding environment affect its distribution. Methanol and C$^{18}$O emission lines at 3 mm have been observed with the IRAM 30m telescope within the large program "Gas phase Elemental abundances in Molecular CloudS" (GEMS) towards 66 positions across 12 starless cores in the Taurus Molecular Cloud. A non-LTE radiative transfer code was used to compute the column densities in all positions. We then used state-of-the-art chemical models to reproduce our observations. We have computed N(CH$_3$OH)/N(C$^{18}$O) column density ratios for all the observed offsets, and two different behaviours can be recognised: the cores where the ratio peaks at the dust peak, and the cores where the ratio peaks with a slight offset with respect to the dust peak ($\sim $10000 AU). We suggest that the cause of this behaviour is the irradiation on the cores due to protostars nearby which accelerate energetic particles along their outflows. The chemical models, which do not take into account irradiation variations, can reproduce fairly well the overall observed column density of methanol, but cannot reproduce the two different radial profiles observed. We confirm the substantial effect of the environment onto the distribution of methanol in starless cores. We suggest that the clumpy medium generated by protostellar outflows might cause a more efficient penetration of the interstellar radiation field in the molecular cloud and have an impact on the distribution of methanol in starless cores.


title: Deep Learning Approach Protecting Privacy in Camera-Based Critical Applications

abstract: Many critical applications rely on cameras to capture video footage for analytical purposes. This has led to concerns about these cameras accidentally capturing more information than is necessary. In this paper, we propose a deep learning approach towards protecting privacy in camera-based systems. Instead of specifying specific objects (e.g. faces) are privacy sensitive, our technique distinguishes between salient (visually prominent) and non-salient objects based on the intuition that the latter is unlikely to be needed by the application.


title: Inductive learning for product assortment graph completion

abstract: Global retailers have assortments that contain hundreds of thousands of products that can be linked by several types of relationships like style compatibility, "bought together", "watched together", etc. Graphs are a natural representation for assortments, where products are nodes and relations are edges. Relations like style compatibility are often produced by a manual process and therefore do not cover uniformly the whole graph. We propose to use inductive learning to enhance a graph encoding style compatibility of a fashion assortment, leveraging rich node information comprising textual descriptions and visual data. Then, we show how the proposed graph enhancement improves substantially the performance on transductive tasks with a minor impact on graph sparsity.


title: Energy Full Counting Statistics and Return to Equilibrium

abstract: We consider a finite dimensional quantum system $S$ in an arbitrary initial state coupled to an infinitely extended quantum thermal reservoir $R$ in equilibrium at inverse temperature $\beta$. The coupling is given by a bounded perturbation of the dynamics and the coupling strength is controlled by a parameter $\lambda$. We assume the system $S + R$ has the property of return to equilibrium, which means that after sufficiently long time, the joint system will have reached equilibrium at inverse temperature $\beta$. In this context, we prove a refinement to the first law of thermodynamics, which states that the total energy of the system and reservoir is conserved. Specifically, we define two measures which encode all the information about the fluctuations of the system and reservoir energy when two measurements are made at time $0$ and time $t$. These measures are called the full counting statistics (FCS). We prove weak convergence of the system and reservoir FCS in the double limit $t \rightarrow \infty$ and $\lambda \rightarrow 0$.


title: Outer billiards in the spaces of oriented geodesics of the three dimensional space forms

abstract: Let $M_{\kappa }$ be the three-dimensional space form of constant curvature $\kappa =0,1,-1$, that is, Euclidean space $\mathbb{R}^{3}$, the sphere $S^{3} $, or hyperbolic space $H^{3}$. Let $S$ be a smooth, closed, strictly convex surface in $M_{\kappa }$. We define an outer billiard map $B$ on the four dimensional space $\mathcal{G}_{\kappa }$ of oriented complete geodesics of $M_{\kappa }$, for which the billiard table is the subset of $\mathcal{G}_{\kappa }$ consisting of all oriented geodesics not intersecting $S$. We show that $B$ is a diffeomorphism when $S$ is quadratically convex. For $\kappa =1,-1$, $\mathcal{G}_{\kappa }$ has a K\"{a}hler structure associated with the Killing form of $\operatorname{Iso}(M_{\kappa })$. We prove that $B$ is a symplectomorphism with respect to its fundamental form and that $B$ can be obtained as an analogue to the construction of Tabachnikov of the outer billiard in $\mathbb{R}^{2n}$ defined in terms of the standard symplectic structure. We show that $B$ does not preserve the fundamental symplectic form on $\mathcal{G}_{\kappa }$ associated with the cross product on $M_{\kappa }$, for $\kappa =0,1,-1$. We initiate the dynamical study of this outer billiard in the hyperbolic case by introducing and discussing a notion of holonomy for periodic points.


title: How You Move Your Head Tells What You Do: Self-supervised Video Representation Learning with Egocentric Cameras and IMU Sensors

abstract: Understanding users' activities from head-mounted cameras is a fundamental task for Augmented and Virtual Reality (AR/VR) applications. A typical approach is to train a classifier in a supervised manner using data labeled by humans. This approach has limitations due to the expensive annotation cost and the closed coverage of activity labels. A potential way to address these limitations is to use self-supervised learning (SSL). Instead of relying on human annotations, SSL leverages intrinsic properties of data to learn representations. We are particularly interested in learning egocentric video representations benefiting from the head-motion generated by users' daily activities, which can be easily obtained from IMU sensors embedded in AR/VR devices. Towards this goal, we propose a simple but effective approach to learn video representation by learning to tell the corresponding pairs of video clip and head-motion. We demonstrate the effectiveness of our learned representation for recognizing egocentric activities of people and dogs.


title: Computable limits of optical multiple-access communications

abstract: Communication rates over quantum channels can be boosted by entanglement, via superadditivity phenomena or entanglement assistance. Superadditivity refers to the capacity improvement from entangling inputs across multiple channel uses. Nevertheless, when unlimited entanglement assistance is available, the entanglement between channel uses becomes unnecessary -- the entanglement-assisted (EA) capacity of a single-sender and single-receiver channel is additive. We generalize the additivity of EA capacity to general multiple-access channels (MACs) for the total communication rate. Furthermore, for optical communication modelled as phase-insensitive bosonic Gaussian MACs, we prove that the optimal total rate is achieved by Gaussian entanglement and therefore can be efficiently evaluated. To benchmark entanglement's advantage, we propose computable outer bounds for the capacity region without entanglement assistance. Finally, we formulate an EA version of minimum entropy conjecture, which leads to the additivity of the capacity region of phase-insensitive bosonic Gaussian MACs if it is true. The computable limits confirm entanglement's boosts in optical multiple-access communications.


title: Microlocal analysis of borehole seismic data

abstract: Borehole seismic data is obtained by receivers located in a well, with sources located on the surface or in another well. Using microlocal analysis, we study possible approximate reconstruction via linearized, filtered backprojection of an isotropic sound speed in the subsurface for three types of data sets. The sources may form a dense array on the surface, or be located along a line on the surface (walkaway geometry) or in another borehole (crosswell). We show that for the dense array, reconstruction is feasible, with no artifacts in the absence of caustics in the background ray geometry, and mild artifacts in the presence of fold caustics in a sense that we define. In contrast, the walkaway and crosswell data sets both give rise to strong, nonremovable artifacts.


title: Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

abstract: Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available.


title: Irreversibility of Structure Tensors of Modules

abstract: Determining the matrix multiplication exponent $\omega$ is one of the greatest open problems in theoretical computer science. We show that it is impossible to prove $\omega = 2$ by starting with structure tensors of modules of fixed degree and using arbitrary restrictions. It implies that the same is impossible by starting with $1_A$-generic non-diagonal tensors of fixed size with minimal border rank. This generalizes the work of Bl\"aser and Lysikov [3]. Our methods come from both commutative algebra and complexity theory.


title: Boomerang webs up to three-loop order

abstract: Webs are sets of Feynman diagrams which manifest soft gluon exponentiation in gauge theory scattering amplitudes: individual webs contribute to the logarithm of the amplitude and their ultraviolet renormalization encodes its infrared structure. In this paper, we consider the particular class of boomerang webs, consisting of multiple gluon exchanges, but where at least one gluon has both of its endpoints on the same Wilson line. First, we use the replica trick to prove that diagrams involving self-energy insertions along the Wilson line do not contribute to the web, i.e. their exponentiated colour factor vanishes. Consequently boomerang webs effectively involve only integrals where boomerang gluons straddle one or more gluons that connect to other Wilson lines. Next we classify and calculate all boomerang webs involving semi-infinite non-lightlike Wilson lines up to three-loop order, including a detailed discussion of how to regulate and renormalize them. Furthermore, we show that they can be written using a basis of specific harmonic polylogarithms, that has been conjectured to be sufficient for expressing all multiple gluon exchange webs. However, boomerang webs differ from other gluon-exchange webs by featuring a lower and non-uniform transcendental weight. We cross-check our results by showing how certain boomerang webs can be determined by the so-called collinear reduction of previously calculated webs. Our results are a necessary ingredient of the soft anomalous dimension for non-lightlike Wilson lines at three loops.


title: Learning, Computing, and Trustworthiness in Intelligent IoT Environments: Performance-Energy Tradeoffs

abstract: An Intelligent IoT Environment (iIoTe) is comprised of heterogeneous devices that can collaboratively execute semi-autonomous IoT applications, examples of which include highly automated manufacturing cells or autonomously interacting harvesting machines. Energy efficiency is key in such edge environments, since they are often based on an infrastructure that consists of wireless and battery-run devices, e.g., e-tractors, drones, Automated Guided Vehicle (AGV)s and robots. The total energy consumption draws contributions from multipleiIoTe technologies that enable edge computing and communication, distributed learning, as well as distributed ledgers and smart contracts. This paper provides a state-of-the-art overview of these technologies and illustrates their functionality and performance, with special attention to the tradeoff among resources, latency, privacy and energy consumption. Finally, the paper provides a vision for integrating these enabling technologies in energy-efficient iIoTe and a roadmap to address the open research challenges


title: The Observation of High-Energy Neutrinos from the Cosmos: Lessons Learned for Multimessenger Astronomy

abstract: The IceCube neutrino telescope discovered PeV-energy neutrinos originating beyond our Galaxy with an energy flux that is comparable to that of GeV-energy gamma rays and EeV-energy cosmic rays. These neutrinos provide the only unobstructed view of the cosmic accelerators that power the highest energy radiation reaching us from the universe. We will review the results from IceCube's first decade of operations, emphasizing the measurement of the diffuse multiflavored neutrino flux from the universe and the identification of the supermassive black hole TXS 0506+056 as a source of cosmic neutrinos and, therefore, cosmic rays. We will speculate on the lessons learned for multimessenger astronomy, among them that extragalactic neutrino sources may be a relatively small subset of the cosmic accelerators observed in high-energy gamma rays and that these may be gamma-ray-obscured at the times that they emit neutrinos.


title: When can relative risks provide causal estimates?

abstract: It is emphasised that for epidemiological studies where disease incidence is rare, results from conventional proportional hazards models can often correctly estimate causal associations. The well-known "backdoor criteria" from causal-inference is applied to the common epidemiological study of rare diseases with a proportional hazards model, providing an example of when and how estimates from conventional proportional hazards studies can be used. A similar study with the "frontdoor criteria", that allows studies with unmeasured confounders, finds similar results to conventional mediation analysis with measured confounders. Reasons for this are discussed.


title: Motion Control of Redundant Robots with Generalised Inequality Constraints

abstract: We present an improved version of the Saturation in the Null Space (SNS) algorithm for redundancy resolution at the velocity level. In addition to hard bounds on joint space motion, we consider also Cartesian box constraints that cannot be violated at any time. The modified algorithm combines all bounds into a single augmented generalised vector and gives equal, highest priority to all inequality constraints. When needed, feasibility of the original task is enforced by the SNS task scaling procedure. Simulation results are reported for a 6R planar robot.


title: Simultaneous Constraints on the Star Formation History and Nucleosynthesis of Sculptor dSph

abstract: We demonstrate that using up to seven stellar abundance ratios can place observational constraints on the star formation histories (SFHs) of Local Group dSphs, using the Sculptor dSph as a test case. We use a one-zone chemical evolution model to fit the overall abundance patterns of $\alpha$ elements (which probe the core-collapse supernovae that occur shortly after star formation), $s$-process elements (which probe AGB nucleosynthesis at intermediate delay times), and iron-peak elements (which probe delayed Type Ia supernovae). Our best-fit model indicates that Sculptor dSph has an ancient SFH, consistent with previous estimates from deep photometry. However, we derive a total star formation duration of $\sim0.9$ Gyr, which is shorter than photometrically-derived SFHs. We explore the effect of various model assumptions on our measurement and find that modifications to these assumptions still produce relatively short SFHs of duration $\lesssim1.4$ Gyr. Our model is also able to compare sets of predicted nucleosynthetic yields for supernovae and AGB stars, and can provide insight into the nucleosynthesis of individual elements in Sculptor dSph. We find that observed [Mn/Fe] and [Ni/Fe] trends are most consistent with sub-$M_{\mathrm{Ch}}$ Type Ia supernova models, and that a combination of "prompt" (delay times similar to core-collapse supernovae) and "delayed" (minimum delay times $\gtrsim50$ Myr) $r$-process events may be required to reproduce observed [Ba/Mg] and [Eu/Mg] trends.


title: AI Chains: Transparent and Controllable Human-AI Interaction by Chaining Large Language Model Prompts

abstract: Although large language models (LLMs) have demonstrated impressive potential on simple tasks, their breadth of scope, lack of transparency, and insufficient controllability can make them less effective when assisting humans on more complex tasks. In response, we introduce the concept of Chaining LLM steps together, where the output of one step becomes the input for the next, thus aggregating the gains per step. We first define a set of LLM primitive operations useful for Chain construction, then present an interactive system where users can modify these Chains, along with their intermediate results, in a modular way. In a 20-person user study, we found that Chaining not only improved the quality of task outcomes, but also significantly enhanced system transparency, controllability, and sense of collaboration. Additionally, we saw that users developed new ways of interacting with LLMs through Chains: they leveraged sub-tasks to calibrate model expectations, compared and contrasted alternative strategies by observing parallel downstream effects, and debugged unexpected model outputs by "unit-testing" sub-components of a Chain. In two case studies, we further explore how LLM Chains may be used in future applications


title: Developing Adaptive Secondary Mirror Concepts for the APF and W.M. Keck Observatory Based on HVR Technology

abstract: An Adaptive secondary mirror (ASM) allows for the integration of adaptive optics (AO) into the telescope itself. Adaptive secondary mirrors, based on hybrid variable reluctance (HVR) actuator technology, developed by TNO, provide a promising path to telescope-integrated AO. HVR actuators have the advantage of allowing mirrors that are stiffer, more power efficient, and potentially less complex than similar, voice-coil based ASM's. We are exploring the application of this technology via a laboratory testbed that will validate the technical approach. In parallel, we are developing conceptual designs for ASMs at several telescopes including the Automated Planet Finder Telescope (APF) and for Keck Observatory. An ASM for APF has the potential to double the light through the slit for radial velocity measurements, and dramatically improved the image stability. An ASM for WMKO enables ground layer AO correction and lower background infrared AO observations, and provides for more flexible deployment of instruments via the ability to adjust the location of the Cassegrain focus.


title: Performance of Large-Format Deformable Mirrors Constructed with Hybrid Variable Reluctance Actuators II: Initial Lab Results from FLASH

abstract: Advancements in high-efficiency hybrid variable reluctance (HVR) actuators are an enabling technology for building the next generation of large-format deformable mirrors, including adaptive secondary mirrors. The Netherlands Organization for Applied Scientific Research (TNO) has developed a new style of hybrid variable reluctance actuator that requires approximately seventy-five times less power to operate as compared to the traditional style of voice-coil actuators. We present the initial performance results from laboratory testing of TNO's latest 19-actuator prototype deformable mirror, FLASH. We report the actuator cross-coupling, linearity, hysteresis, natural shape flattening, and drift as measured with a Zygo interferometer and a set of four capacitive sensors. We also present results of the dynamic performance of the FLASH on sub-millisecond timescales to estimate the limits of this technology for high-contrast imaging adaptive optics. We confirm that this technology has strong potential for use in on-sky adaptive secondary mirrors without the need for active cooling.


title: The weak Ramsey property and extreme amenability

abstract: We extend the Kechris--Pestov--Todor\v{c}evi\'c correspondence to weak Fra\"{\i}ss\'{e} categories and automorphism groups of generic objects. The new ingredient is the weak Ramsey property. We demonstrate the theory on several examples including monoid categories, the category of almost linear orders, and categories of strong embeddings of trees.


title: Cosmological constraints on slow-roll inflation: an update

abstract: In light of the most recent cosmological observations, we provide new updated constraints on the slow-roll inflation in different extended scenarios beyond the $\Lambda\rm{CDM}$ cosmological model. Along with the usual six parameters, we simultaneously vary different combinations of additional parameters, including the running of the scalar spectral index $\alpha_s$, its running of running $\beta_s$, the tensor amplitude $r$ and the spatial curvature $\Omega_k$. From the Planck 2018 data, we find no evidence for a scalar running or a running of running, while analyzing the Atacama Cosmology Telescope data combined with WMAP 9-years observations data we find a preference for non-zero $\alpha_s$ and $\beta_s$ at the level of 2.9$\sigma$ and 2.7$\sigma$, respectively. Anyway, this preference is reduced when the tensor amplitude can vary in the model or $\beta_s$ is fixed to zero. The upper bound on $r$ is only slightly affected by the additional parameters while the differences in the datasets can remarkably change the compatibility among the different inflationary models, sometimes leading to discordant conclusions.


title: Critical dynamics of relativistic diffusion

abstract: We study the dynamics of self-interacting scalar fields with $Z_2$ symmetry governed by a relativistic Israel-Stuart type diffusion equation in the vicinity of a critical point. We calculate spectral functions of the order parameter in mean-field approximation as well as using first-principles classical-statistical lattice simulations in real-time. We observe that the spectral functions are well-described by single Breit-Wigner shapes. Away from criticality, the dispersion matches the expectations from the mean-field approach. At the critical point, the spectral functions largely keep their Breit-Wigner shape, albeit with non-trivial power-law dispersion relations. We extract the characteristic time-scales as well as the dynamic critical exponent $z$, verifying the existence of a dynamic scaling regime. In addition, we derive the universal scaling functions implied by the Breit-Wigner shape with critical power-law dispersion and show that they match the data. Considering equations of motion for a system coupled to a heat bath as well as an isolated system, we perform this study for two different dynamic universality classes, both in two and three spatial dimensions.

categories: hep-lat cond-mat.stat-mech hep-ph nucl-th


title: Bilevel hyperparameter optimization for support vector classification: theoretical analysis and a solution method

abstract: Support vector classification (SVC) is a classical and well-performed learning method for classification problems. A regularization parameter, which significantly affects the classification performance, has to be chosen and this is usually done by the cross-validation procedure. In this paper, we reformulate the hyperparameter selection problem for support vector classification as a bilevel optimization problem in which the upper-level problem minimizes the average number of misclassified data points over all the cross-validation folds, and the lower-level problems are the l1-loss SVC problems, with each one for each fold in T-fold cross-validation. The resulting bilevel optimization model is then converted to a mathematical program with equilibrium constraints (MPEC). To solve this MPEC, we propose a global relaxation cross-validation algorithm (GR-CV) based on the well-know Sholtes-type global relaxation method (GRM). It is proven to converge to a C-stationary point. Moreover, we prove that the MPEC-tailored version of the Mangasarian-Fromovitz constraint qualification (MFCQ), which is a key property to guarantee the convergence of the GRM, automatically holds at each feasible point of this MPEC. Extensive numerical results verify the efficiency of the proposed approach. In particular, compared with other methods, our algorithm enjoys superior generalization performance over almost all the data sets used in this paper.


title: HYPPO: A Surrogate-Based Multi-Level Parallelism Tool for Hyperparameter Optimization

abstract: We present a new software, HYPPO, that enables the automatic tuning of hyperparameters of various deep learning (DL) models. Unlike other hyperparameter optimization (HPO) methods, HYPPO uses adaptive surrogate models and directly accounts for uncertainty in model predictions to find accurate and reliable models that make robust predictions. Using asynchronous nested parallelism, we are able to significantly alleviate the computational burden of training complex architectures and quantifying the uncertainty. HYPPO is implemented in Python and can be used with both TensorFlow and PyTorch libraries. We demonstrate various software features on time-series prediction and image classification problems as well as a scientific application in computed tomography image reconstruction. Finally, we show that (1) we can reduce by an order of magnitude the number of evaluations necessary to find the most optimal region in the hyperparameter space and (2) we can reduce by two orders of magnitude the throughput for such HPO process to complete.


title: Mediated interactions beyond the nearest neighbor in an array of superconducting qubits

abstract: We consider mediated interactions in an array of floating transmons, where each qubit capacitor consists of two superconducting pads galvanically isolated from ground. Each such pair contributes two quantum degrees of freedom, one of which is used as a qubit, while the other remains fixed. However, these extraneous modes can generate coupling between the qubit modes that extends beyond the nearest neighbor. We present a general formalism describing the formation of this coupling and calculate it for a one-dimensional chain of transmons. We show that the strength of coupling and its range (that is, the exponential falloff) can be tuned independently via circuit design to realize a continuum from nearest-neighbor-only interactions to interactions that extend across the length of the chain. We present designs with capacitance and microwave simulations showing that various interaction configurations can be achieved in realistic circuits. Such coupling could be used in analog simulation of different quantum regimes or to increase connectivity in digital quantum systems. Thus mechanism must also be taken into account in other types of qubits with extraneous modes.


title: On the Maximum Achievable Sum-rate of the RIS-aided MIMO Broadcast Channel

abstract: Reconfigurable intelligent surfaces (RISs) represent a new technology that can shape the radio wave propagation and thus offers a great variety of possible performance and implementation gains. Motivated by this, we investigate the achievable sum-rate optimization in a broadcast channel (BC) in the presence of RISs. We solve this problem by exploiting the well-known duality between the Gaussian multiple-input multiple-output (MIMO) BC and the multiple-access channel (MAC), and we correspondingly derive three algorithms which optimize the users' covariance matrices and the RIS phase shifts in the dual MAC. The users' covariance matrices are optimized by a dual decomposition method with block coordinate maximization (BCM), or by a gradient-based method. The RIS phase shifts are either optimized sequentially by using a closed-form expression, or are computed in parallel by using a gradient-based method. We present a computational complexity analysis for the proposed algorithms. Simulation results show that the proposed algorithms tend to converge to the same achievable sum-rate overall, but may produce different sum-rate performance for some specific situations, due to the non-convexity of the considered problem. Also, the gradient-based optimization methods are generally more time efficient. In addition, we demonstrate that the proposed algorithms can provide a significant gain in the RIS-assisted BC assisted by multiple RISs and that the gain depends on the placement of the RISs.


title: Multiscale Mechanical Consequences of Ocean Acidification for Cold-Water Corals

abstract: Ocean acidification is a threat to deep-sea corals and could lead to dramatic and rapid loss of the reef framework habitat they build. Weakening of structurally critical parts of the coral reef framework can lead to physical habitat collapse on an ecosystem scale, reducing the potential for biodiversity support. The mechanism underpinning crumbling and collapse of corals can be described via a combination of laboratory-scale experiments and mathematical and computational models. We synthesise data from electron back-scatter diffraction, micro-computed tomography, and micromechanical experiments, supplemented by molecular dynamics and continuum micromechanics simulations to predict failure of coral structures under increasing porosity and dissolution. Results reveal remarkable mechanical properties of cold-water coral skeletons of 462 MPa compressive strength and 45-67 GPa stiffness. This is 10 times stronger than concrete, twice as strong than ultrahigh performance fibre reinforced concrete, or nacre. Contrary to what would be expected, CWCs skeletons retain their strength despite a loss of stiffness and even when synthesised under future oceanic conditions. Our models capture the impact of corrosive waters on exposed skeletons and illustrate how small modifications in their skeleton lead to significantly increased risk of crumbling coral habitat. This new understanding, combined with projections of how seawater chemistry will change over the coming decades, will help support future conservation and management efforts of these vulnerable marine ecosystems by identifying which ecosystems are at risk and when they will be at risk, allowing assessment of the impact upon associated biodiversity.


title: Using C++ to Calculate $\mathsf{SO(10)}$ Tensor Couplings

abstract: \noindent Model building in $\mathsf{SO(10)}$, which is the leading grand unification framework, often involves large Higgs representations and their couplings. Explicit calculations of such couplings is a multi-step process that involves laborious calculations that are time-consuming and error-prone, an issue which only grows as the complexity of the coupling increases. Therefore, there exists an opportunity to leverage the abilities of computer software in order to algorithmically perform these calculations on-demand. This paper outlines the details of such a software, implemented in C++ using in-built libraries. The software is capable of accepting invariant couplings involving an arbitrary number of $\mathsf{SO(10)}$ Higgs tensors, each having up to 5 indices. The output is then produced in \LaTeX, so that it is universally readable and sufficiently expressive. Through the use of this software, $\mathsf{SO(10)}$ coupling analysis can be performed in a way that minimizes calculation time, eliminates errors, and allows for experimentation with couplings that have not been computed before in the literature. Furthermore, this software can be expanded in the future to account for similar Higgs-Spinor coupling analysis, or extended to include further $\mathsf{SO(N)}$ invariant couplings.


title: Affine dimers from characteristic polygons

abstract: Recent work by Forsg{\aa}rd indicates that not every convex lattice polygon arises as the characteristic polygon of an affine dimer or, equivalently, an admissible oriented line arrangement on the torus in general position. We begin the classification of convex lattice polygons arising as characteristic polygons of affine dimers. We present several general constructions of new affine dimers from old, and an algorithm for finding affine dimers with prescribed polygon. With these tools we prove that all lattice triangles, generalised parallelograms, and polygons of genus at most two admit an affine dimer.


title: Predicting far-infrared maps of galaxies via machine learning techniques

abstract: The ultraviolet (UV) to sub-millimetre (submm) spectral energy distribution of galaxies can be roughly divided into two sections: the stellar emission (attenuated by dust) at UV to near-infrared wavelengths and dust emission at longer wavelengths. In Dobbels et al. (2020), we show that these two sections are strongly related, and we can predict the global dust properties from the integrated UV to mid-infrared emission with the help of machine learning techniques. We investigate if these machine learning techniques can also be extended to resolved scales. Our aim is to predict resolved maps of the specific dust luminosity, specific dust mass, and dust temperature starting from a set of surface brightness images from UV to mid-infrared wavelengths. We used a selection of nearby galaxies retrieved from the DustPedia sample, in addition to M31 and M33. These were convolved and resampled to a range of pixel sizes, ranging from 150 pc to 3 kpc. We trained a random forest model which considers each pixel individually. We find that the predictions work well on resolved scales, with the dust mass and temperature having a similar root mean square error as on global scales (0.32 dex and 3.15 K on 18" scales respectively), and the dust luminosity being noticeably better (0.11 dex). We find no significant dependence on the pixel scale. Predictions on individual galaxies can be biased, and we find that about two-thirds of the scatter can be attributed to scatter between galaxies (rather than within galaxies). A machine learning approach can be used to create dust maps, with its resolution being only limited to the input bands, thus achieving a higher resolution than Herschel. These dust maps can be used to improve global estimates of dust properties, they can lead to a better estimate of dust attenuation, and they can be used as a constraint on cosmological simulations that trace dust.


title: Let there be a clock on the beach: Reducing Object Hallucination in Image Captioning

abstract: Explaining an image with missing or non-existent objects is known as object bias (hallucination) in image captioning. This behaviour is quite common in the state-of-the-art captioning models which is not desirable by humans. To decrease the object hallucination in captioning, we propose three simple yet efficient training augmentation method for sentences which requires no new training data or increase in the model size. By extensive analysis, we show that the proposed methods can significantly diminish our models' object bias on hallucination metrics. Moreover, we experimentally demonstrate that our methods decrease the dependency on the visual features. All of our code, configuration files and model weights will be made public.


title: Event-based hyperspectral EELS: towards nanosecond temporal resolution

abstract: The acquisition of a hyperspectral image is nowadays a standard technique used in the scanning transmission electron microscope. It relates the spatial position of the electron probe to the spectral data associated with it. In the case of electron energy loss spectroscopy (EELS), frame-based hyperspectral acquisition is much slower than the achievable rastering time of the scan unit (SU), which sometimes leads to undesirable effects in the sample, such as electron irradiation damage, that goes unperceived during frame acquisition. In this work, we have developed an event-based hyperspectral EELS by using a Timepix3 application-specific integrated circuit detector with two supplementary time-to-digital (TDC) lines embedded. In such a system, electron events are characterized by their positional and temporal coordinates, but TDC events only by temporal ones. By sending reference signals from the SU to the TDC line, it is possible to reconstruct the entire spectral image with SU-limited scanning pixel dwell time and thus acquire, with no additional cost, a hyperspectral image at the same rate as that of a single channel detector, such as annular dark-field. To exemplify the possibilities behind event-based hyperspectral EELS, we have studied the decomposition of calcite (CaCO$_3$) into calcium oxide (CaO) and carbon dioxide (CO$_2$) under the electron beam irradiation.


title: The Limits of Optimal Pricing in the Dark

abstract: A ubiquitous learning problem in today's digital market is, during repeated interactions between a seller and a buyer, how a seller can gradually learn optimal pricing decisions based on the buyer's past purchase responses. A fundamental challenge of learning in such a strategic setup is that the buyer will naturally have incentives to manipulate his responses in order to induce more favorable learning outcomes for him. To understand the limits of the seller's learning when facing such a strategic and possibly manipulative buyer, we study a natural yet powerful buyer manipulation strategy. That is, before the pricing game starts, the buyer simply commits to "imitate" a different value function by pretending to always react optimally according to this imitative value function. We fully characterize the optimal imitative value function that the buyer should imitate as well as the resultant seller revenue and buyer surplus under this optimal buyer manipulation. Our characterizations reveal many useful insights about what happens at equilibrium. For example, a seller with concave production cost will obtain essentially 0 revenue at equilibrium whereas the revenue for a seller with convex production cost is the Bregman divergence of her cost function between no production and certain production. Finally, and importantly, we show that a more powerful class of pricing schemes does not necessarily increase, in fact, may be harmful to, the seller's revenue. Our results not only lead to an effective prescriptive way for buyers to manipulate learning algorithms but also shed lights on the limits of what a seller can really achieve when pricing in the dark.


title: A DEIM driven reduced basis method for the diffuse Stokes/Darcy model coupled at parametric phase-field interfaces

abstract: In this article, we develop a reduced basis method for efficiently solving the coupled Stokes/Darcy equations with parametric internal geometry. To accommodate possible changes in topology, we define the Stokes and Darcy domains implicitly via a phase-field indicator function. In our reduced order model, we approximate the parameter-dependent phase-field function with a discrete empirical interpolation method (DEIM) that enables affine decomposition of the associated linear and bilinear forms. In addition, we introduce a modification of DEIM that leads to non-negativity preserving approximations, thus guaranteeing positive-semidefiniteness of the system matrix. We also present a strategy for determining the required number of DEIM modes for a given number of reduced basis functions. We couple reduced basis functions on neighboring patches to enable the efficient simulation of large-scale problems that consist of repetitive subdomains. We apply our reduced basis framework to efficiently solve the inverse problem of characterizing the subsurface damage state of a complete in-situ leach mining site.


title: Benchmarking Memory-Centric Computing Systems: Analysis of Real Processing-in-Memory Hardware

abstract: Many modern workloads such as neural network inference and graph processing are fundamentally memory-bound. For such workloads, data movement between memory and CPU cores imposes a significant overhead in terms of both latency and energy. A major reason is that this communication happens through a narrow bus with high latency and limited bandwidth, and the low data reuse in memory-bound workloads is insufficient to amortize the cost of memory access. Fundamentally addressing this data movement bottleneck requires a paradigm where the memory system assumes an active role in computing by integrating processing capabilities. This paradigm is known as processing-in-memory (PIM). Recent research explores different forms of PIM architectures, motivated by the emergence of new technologies that integrate memory with a logic layer, where processing elements can be easily placed. Past works evaluate these architectures in simulation or, at best, with simplified hardware prototypes. In contrast, the UPMEM company has designed and manufactured the first publicly-available real-world PIM architecture. The UPMEM PIM architecture combines traditional DRAM memory arrays with general-purpose in-order cores, called DRAM Processing Units (DPUs), integrated in the same chip. This paper presents key takeaways from the first comprehensive analysis of the first publicly-available real-world PIM architecture. We provide four key takeaways about the UPMEM PIM architecture, which stem from our study. More insights about suitability of different workloads to the PIM system, programming recommendations for software designers, and suggestions and hints for hardware and architecture designers of future PIM systems are available in arXiv:2105.03814


title: PyTorrent: A Python Library Corpus for Large-scale Language Models

abstract: A large scale collection of both semantic and natural language resources is essential to leverage active Software Engineering research areas such as code reuse and code comprehensibility. Existing machine learning models ingest data from Open Source repositories (like GitHub projects) and forum discussions (like Stackoverflow.com), whereas, in this showcase, we took a step backward to orchestrate a corpus titled PyTorrent that contains 218,814 Python package libraries from PyPI and Anaconda environment. This is because earlier studies have shown that much of the code is redundant and Python packages from these environments are better in quality and are well-documented. PyTorrent enables users (such as data scientists, students, etc.) to build off the shelf machine learning models directly without spending months of effort on large infrastructure. The dataset, schema and a pretrained language model is available at: https://github.com/fla-sil/PyTorrent


title: LazySets.jl: Scalable Symbolic-Numeric Set Computations

abstract: LazySets.jl is a Julia library that provides ways to symbolically represent sets of points as geometric shapes, with a special focus on convex sets and polyhedral approximations. LazySets provides methods to apply common set operations, convert between different set representations, and efficiently compute with sets in high dimensions using specialized algorithms based on the set types. LazySets is the core library of JuliaReach, a cutting-edge software addressing the fundamental problem of reachability analysis: computing the set of states that are reachable by a dynamical system from all initial states and for all admissible inputs and parameters. While the library was originally designed for reachability and formal verification, its scope goes beyond such topics. LazySets is an easy-to-use, general-purpose and scalable library for computations that mix symbolics and numerics. In this article we showcase the basic functionality, highlighting some of the key design choices.


title: Feedback control of social distancing for COVID-19 via elementary formulae

abstract: Social distancing has been enacted in order to mitigate the spread of COVID-19. Like many authors, we adopt the classic epidemic SIR model, where the infection rate is the control variable. Its differential flatness property yields ele mentary closed-form formulae for open-loop social distancing scenarios, where, for instance, the increase of the number of uninfected people may be taken into account. Those formulae might therefore be useful to decision makers. A feedback loop stemming from model-free control leads to a remarkable robustness with respect to severe uncertainties and mismatches. Although an identification procedure is presented, a good knowledge of the recovery rate is not necessary for our control strategy.


title: Density functional Bogoliubov-de Gennes analysis of superconducting Nb and Nb(110) surfaces

abstract: We report on the implementation of the Bogoliubov-de Gennes method into the JuKKR code [https: //jukkr.fz-juelich.de], an implementation of the relativistic all-electron, full-potential Korringa-Kohn-Rostoker Green function method, which allows a material-specific description of inhomogeneous super-conductors and heterostructures on the basis of density functional theory. We describe the formalism and report on calculations for the s-wave superconductor Nb, a potential component of the materials platform enabling the realization of the Majorana zero modes in the field of topological quantum computing. We compare the properties of the superconducting state both in the bulk and for (110) surfaces. We compare slab calculations for different thicknesses and comment on the importance of spin-orbit coupling, the effect of surface relaxations and the influence of a softening of phonon modes on the surface for the resulting superconducting gap.


title: Anomalous temperature dependence of phonon pumping by ferromagnetic resonance in Co/Pd multilayers with perpendicular anisotropy

abstract: We demonstrate the pumping of phonons by ferromagnetic resonance in a series of [Co(0.8 nm)/Pd(1.5 nm)]$_n$ multilayers ($n =$ 6, 11, 15, and 20) with large magnetostriction and perpendicular magnetic anisotropy. The effect is shown using broadband ferromagnetic resonance over a range of temperatures (10 to 300 K), where a resonant damping enhancement is observed at frequencies corresponding to standing wave phonons across the multilayer. The strength of this effect is enhanced by approximately a factor of 4 at 10 K compared to room temperature, which is anomalous in the sense that the temperature dependence of the magnetostriction predicts an enhancement that is less than a factor of 2. Lastly, we demonstrate that the damping enhancement is correlated with a shift in the ferromagnetic resonance field as predicted quantitatively from linear response theory.


title: Epidemiological parameter sensitivity in Covid-19 dynamics and estimation

abstract: Covid-19 is one of the most dreaded pandemics/epidemics in the world threatening the human population. The dynamics of this pandemic is quite complicated and prediction of pandemic states often fails. In this work, we study and correlate the SIR epidemiological model with the ongoing pandemic and found that pandemic dynamics and states are quite sensitively dependent on model parameters. The analysis of the exact parametric solution of the deterministic SIR model shows that the fixed points ($g^*$) depend on the SIR parameters, where, if $g^*\langle\frac{\gamma}{\beta x_0}$ then $g^*$ is stable and the pandemic can be controlled, whereas, if $g^*\rangle\frac{\gamma}{\beta x_0}$ then $g^*$ is unstable indicating active pandemic state, and $g^*=0$ corresponds to an endemic state. The dynamics show asymptotic stability bifurcation. The analytical solution of the stochastic SIR model allows to predict endemic time $t_E\propto\frac{1}{\gamma}$ and mean infected population becomes constant for small values of $\gamma$. The noise in the system even can modulate pandemic dynamics. Hence, we estimated the values of these parameters during lockdown periods using the ABC SSA algorithm and found strong dependence on control strategy, namely, lockdown and follows power-law behaviors. The numerical results show strong agreement with the real data and have the positive effect of social distancing or lockdown in the pandemic control.


title: Highly Optimized Full-Core Reactor Simulations on Summit

abstract: Nek5000/RS is a highly-performant open-source spectral element code for simulation of incompressible and low-Mach fluid flow, heat transfer, and combustion with a particular focus on turbulent flows in complex domains. It is based on high-order discretizations that realize the same (or lower) cost per gridpoint as traditional low-order methods. State-of-the-art multilevel preconditioners, efficient high-order time-splitting methods, and runtime-adaptive communication strategies are built on a fast OCCA-based kernel library, libParanumal, to provide scalability and portability across the spectrum of current and future high-performance computing platforms. On Summit, Nek5000/RS has recently achieved an milestone in the simulation of nuclear reactors: the first full-core computational fluid dynamics simulations of reactor cores, including pebble beds with > 350,000 pebbles and 98M elements advanced in less than 0.25 seconds per Navier-Stokes timestep. With carefully tuned algorithms, it is possible to simulate a single flow-through time for a full reactor core in less than six hours on all of Summit.


title: Molecule3D: A Benchmark for Predicting 3D Geometries from Molecular Graphs

abstract: Graph neural networks are emerging as promising methods for modeling molecular graphs, in which nodes and edges correspond to atoms and chemical bonds, respectively. Recent studies show that when 3D molecular geometries, such as bond lengths and angles, are available, molecular property prediction tasks can be made more accurate. However, computing of 3D molecular geometries requires quantum calculations that are computationally prohibitive. For example, accurate calculation of 3D geometries of a small molecule requires hours of computing time using density functional theory (DFT). Here, we propose to predict the ground-state 3D geometries from molecular graphs using machine learning methods. To make this feasible, we develop a benchmark, known as Molecule3D, that includes a dataset with precise ground-state geometries of approximately 4 million molecules derived from DFT. We also provide a set of software tools for data processing, splitting, training, and evaluation, etc. Specifically, we propose to assess the error and validity of predicted geometries using four metrics. We implement two baseline methods that either predict the pairwise distance between atoms or atom coordinates in 3D space. Experimental results show that, compared with generating 3D geometries with RDKit, our method can achieve comparable prediction accuracy but with much smaller computational costs. Our Molecule3D is available as a module of the MoleculeX software library (https://github.com/divelab/MoleculeX).


title: Randomized Projection Learning Method forDynamic Mode Decomposition

abstract: A data-driven analysis method known as dynamic mode decomposition (DMD) approximates the linear Koopman operator on projected space. In the spirit of Johnson-Lindenstrauss Lemma, we will use random projection to estimate the DMD modes in reduced dimensional space. In practical applications, snapshots are in high dimensional observable space and the DMD operator matrix is massive. Hence, computing DMD with the full spectrum is infeasible, so our main computational goal is estimating the eigenvalue and eigenvectors of the DMD operator in a projected domain. We will generalize the current algorithm to estimate a projected DMD operator. We focus on a powerful and simple random projection algorithm that will reduce the computational and storage cost. While clearly, a random projection simplifies the algorithmic complexity of a detailed optimal projection, as we will show, generally the results can be excellent nonetheless, and quality understood through a well-developed theory of random projections. We will demonstrate that modes can be calculated for a low cost by the projected data with sufficient dimension. Keyword: Koopman Operator, Dynamic Mode Decomposition(DMD), Johnson-Lindenstrauss Lemma, Random Projection, Data-driven method.


title: Semiclassical gravity with a conformally covariant field in globally hyperbolic spacetimes

abstract: We prove that semiclassical gravity in conformally static, globally hyperbolic spacetimes with a massless, conformally coupled Klein-Gordon field is well posed, when viewed as a coupled theory for the dynamical conformal factor of the metric and the Klein-Gordon theory. Namely, it admits unique and stable solutions whenever constrained fourth-order initial data for the conformal factor and suitably defined Hadamard initial data for the Klein-Gordon state are provided on a spacelike Cauchy surface. As no spacetime symmetries are imposed on the conformal factor, the present result implies that, provided constrained initial data exists, there also exist exact solutions to the semiclassical gravity equations beyond the isotropic, homogeneous or static cases.


title: Separation of Artifacts from Spin-Torque Ferromagnetic Resonance Measurements of Spin-Orbit Torque for the Low-Symmetry van der Waals Semi-Metal ZrTe$_\textbf{3}$

abstract: We measure spin-orbit torque generated by exfoliated layers of the low-symmetry semi-metal ZrTe$_3$ using the spin-torque ferromagnetic resonance (ST-FMR) technique. When the ZrTe$_3$ has a thickness greater than about 10 nm, artifacts due to spin pumping and/or resonant heating can cause the standard ST-FMR analysis to overestimate the true magnitude of the torque efficiency by as much as a factor of 30, and to indicate incorrectly that the spin-orbit torque depends strongly on the ZrTe$_3$ layer thickness. Artifact-free measurements can still be achieved over a substantial thickness range by the method developed recently to detect ST-FMR signals in the Hall geometry as well as the longitudinal geometry. ZrTe$_3$/Permalloy samples generate a conventional in-plane anti-damping spin torque efficiency $\xi_{||}^{\text{DL}}$ = 0.014 $\pm$ 0.004, and an unconventional in-plane field-like torque efficiency $|\xi_{||}^{\text{FL}}|$ = 0.003 $\pm$ 0.001. The out-of-plane anti-damping torque is negligible. We suggest that artifacts similarly interfere with the standard ST-FMR analysis for other van der Waals samples thicker than about 10 nm.


title: Image reconstruction from tissue scattered events for ${\beta}^{+}\gamma$ coincidences in Compton-PET

abstract: For long time non-pure beta emitters are avoided from PET imaging due to extra dose and increase in background from Compton scattering. But advent of high-resolution Compton camera system opens up new domain of imaging. Various non-pure beta emitters are formed as beam irradiation byproduct in therapy which can be used in online beam range verification. In this case, the number of usable counts for imaging is generally 1-3 order lesser than normal PET scan. On the other hand, we know that in human PET scanner, 30-60\% can be tissue scattered coincidences in 3D case containing 80\% single scattered events. In this work, we have investigated feasibility of imaging using only single scattered coincidences for non-pure beta emitters in a Compton-PET system. The locus of tissue scatter point can be reduced to in generally two points after using Compton cone from both ends of 511 keV detections. Finally, annihilation point is estimated using Compton cone of 1157 keV gamma and time-of-flight information for the 511 keV. We believe independent assessment of underlying activity from single scattered data sets will increase confidence in image interpretation.


title: Wireless Link Scheduling via Graph Representation Learning: A Comparative Study of Different Supervision Levels

abstract: We consider the problem of binary power control, or link scheduling, in wireless interference networks, where the power control policy is trained using graph representation learning. We leverage the interference graph of the wireless network as an underlying topology for a graph neural network (GNN) backbone, which converts the channel matrix to a set of node embeddings for all transmitter-receiver pairs. We show how the node embeddings can be trained in several ways, including via supervised, unsupervised, and self-supervised learning, and we compare the impact of different supervision levels on the performance of these methods in terms of the system-level throughput, convergence behavior, sample efficiency, and generalization capability.


title: Density maximizers of layered permutations

abstract: A permutation is layered if it contains neither 231 nor 312 as a pattern. It is known that, if $\sigma$ is a layered permutation, then the density of $\sigma$ in a permutation of order $n$ is maximized by a layered permutation. Albert, Atkinson, Handley, Holton and Stromquist [Electron. J. Combin. 9 (2002), R#5] claimed that the density of a layered permutation with layers of sizes $(a,1,b)$ where $a,b\geq2$ is asymptotically maximized by layered permutations with a bounded number of layers, and conjectured that the same holds if a layered permutation has no consecutive layers of size one and its first and last layers are of size at least two. We show that, if $\sigma$ is a layered permutation whose first layer is sufficiently large and second layer is of size one, then the number of layers tends to infinity in every sequence of layered permutations asymptotically maximizing the density of $\sigma$. This disproves the conjecture and the claim of Albert et al. We complement this result by giving sufficient conditions on a layered permutation to have asymptotic or exact maximizers with a bounded number of layers.


title: Optimization of the resonator-induced phase gate for superconducting qubits

abstract: The resonator-induced phase gate is a two-qubit operation in which driving a bus resonator induces a state-dependent phase shift on the qubits equivalent to an effective $ZZ$ interaction. In principle, the dispersive nature of the gate offers flexibility for qubit parameters. However, the drive can cause resonator and qubit leakage, the physics of which cannot be fully captured using either the existing Jaynes-Cummings or Kerr models. In this paper, we adopt an ab-initio model based on Josephson nonlinearity for transmon qubits. The ab-initio analysis agrees well with the Kerr model in terms of capturing the effective $ZZ$ interaction in the weak-drive dispersive regime. In addition, however, it reveals numerous leakage transitions involving high-excitation qubit states. We analyze the physics behind such novel leakage channels, demonstrate the connection with specific qubits-resonator frequency collisions, and lay out a plan towards device parameter optimization. We show this type of leakage can be substantially suppressed using very weakly anharmonic transmons. In particular, weaker qubit anharmonicity mitigates both collision density and leakage amplitude, while larger qubit frequency moves the collisions to occur only at large anharmonicity not relevant to experiment. Our work is broadly applicable to the physics of weakly anharmonic transmon qubits coupled to linear resonators. In particular, our analysis confirms and generalizes the measurement-induced state transitions noted in Sank et al. (Phys. Rev. Lett. 117, 190503) and lays the groundwork for both strong-drive resonator-induced phase gate implementation and strong-drive dispersive qubit measurement.


title: LLOL: Low-Latency Odometry for Spinning Lidars

abstract: In this paper, we present a low-latency odometry system designed for spinning lidars. Many existing lidar odometry methods wait for an entire sweep from the lidar before processing the data. This introduces a large delay between the first laser firing and its pose estimate. To reduce this latency, we treat the spinning lidar as a streaming sensor and process packets as they arrive. This effectively distributes expensive operations across time, resulting in a very fast and lightweight system with much higher throughput and lower latency. Our open-source implementation is available at \url{https://github.com/versatran01/llol}.


title: On-Demand Networking for Ubiquitous Connectivity and Network Resilience: A Network-in-a-Box Solution

abstract: Recently, the wireless community has initiated research on the sixth generation (6G) cellular network for the next decade. The 6G visions are still under development but are converging toward ubiquitous, sustainable, and automated digital society. A network-in-a-box (NIB) is a portable and fully-fledged networking solution that has many potentials to stimulate 6G visions, especially for ubiquitous and resilient network connectivity. In this article, we highlight how NIB features suit 6G use cases and requirements and how it can be used for 6G communications. In addition, we discuss the challenges of the potential enabling technologies of 6G that can reinforce the NIB performance.


title: Multimodal Driver State Modeling through Unsupervised Learning

abstract: Naturalistic driving data (NDD) can help understand drivers' reactions to each driving scenario and provide personalized context to driving behavior. However, NDD requires a high amount of manual labor to label certain driver's state and behavioral patterns. Unsupervised analysis of NDD can be used to automatically detect different patterns from the driver and vehicle data. In this paper, we propose a methodology to understand changes in driver's physiological responses within different driving patterns. Our methodology first decomposes a driving scenario by using a Bayesian Change Point detection model. We then apply the Latent Dirichlet Allocation method on both driver state and behavior data to detect patterns. We present two case studies in which vehicles were equipped to collect exterior, interior, and driver behavioral data. Four patterns of driving behaviors (i.e., harsh brake, normal brake, curved driving, and highway driving), as well as two patterns of driver's heart rate (HR) (i.e., normal vs. abnormal high HR), and gaze entropy (i.e., low versus high), were detected in these two case studies. The findings of these case studies indicated that among our participants, the drivers' HR had a higher fraction of abnormal patterns during harsh brakes, accelerating and curved driving. Additionally, free-flow driving with close to zero accelerations on the highway was accompanied by more fraction of normal HR as well as a lower gaze entropy pattern. With the proposed methodology we can better understand variations in driver's psychophysiological states within different driving scenarios. The findings of this work, has the potential to guide future autonomous vehicles to take actions that are fit to each specific driver.


title: A note on the two dimensional Lagrangian mean curvature equation

abstract: In this note, we use Warren-Yuan's super isoperimetric inequality on the level sets of subharmonic functions, which is available only in two dimensions, to derive a modified Hessian bound for solutions of the two dimensional Lagrangian mean curvature equation. We assume the Lagrangian phase to be supercritical with bounded second derivatives. Unlike the previous approach, the simplified approach in this proof does not require the Michael-Simon mean value and Sobolev inequalities on generalized submanifolds of $\mathbb{R}^n$.


title: Multilevel orthogonal Bochner function subspaces with applications to robust machine learning

abstract: In our approach, we consider the data as instances of a random field within a relevant Bochner space. Our key observation is that the classes can predominantly reside in two distinct subspaces. To uncover the separation between these classes, we employ the Karhunen-Loeve expansion and construct the appropriate subspaces. This allows us to effectively reveal the distinction between the classes. The novel features forming the above bases are constructed by applying a coordinate transformation based on the recent Functional Data Analysis theory for anomaly detection. The associated signal decomposition is an exact hierarchical tensor product expansion with known optimality properties for approximating stochastic processes (random fields) with finite dimensional function spaces. Using a hierarchical finite dimensional expansion of the nominal class, a series of orthogonal nested subspaces is constructed for detecting anomalous signal components. Projection coefficients of input data in these subspaces are then used to train a Machine Learning (ML classifier. However, due to the split of the signal into nominal and anomalous projection components, clearer separation surfaces for the classes arise. In fact we show that with a sufficiently accurate estimation of the covariance structure of the nominal class, a sharp classification can be obtained. This is particularly advantageous for large unbalanced datasets. We demonstrate it on a number of high-dimensional datasets. This approach yields significant increases in accuracy of ML methods compared to using the same ML algorithm with the original feature data. Our tests on the Alzheimer's Disease ADNI dataset shows a dramatic increase in accuracy (from 48% to 89% accuracy). Furthermore, tests using unbalanced semi-synthetic datasets created from the benchmark GCM dataset confirm increased accuracy as the dataset becomes more unbalanced.


title: Pre-Quantized Deep Learning Models Codified in ONNX to Enable Hardware/Software Co-Design

abstract: This paper presents a methodology to separate the quantization process from the hardware-specific model compilation stage via a pre-quantized deep learning model description in standard ONNX format. Separating the quantization process from the model compilation stage enables independent development. The methodology is expressive to convey hardware-specific operations and to embed key quantization parameters into a ONNX model which enables hardware/software co-design. Detailed examples are given for both MLP and CNN based networks, which can be extended to other networks in a straightforward fashion.


title: Heritable Nongenetic Information in the Form of the DNA-Autonomous Tissue Spatial Code that Governs Organismal Development, Tissue Regeneration, and Tumor Architecture

abstract: Numerous studies of the tumor microenvironment, interspecies xenografting, and limb regeneration suggest the existence of a tissue spatial code [TSC] that controls tissue structure in a quasi "epigenetic" fashion. Epigenetic is an inadequate label, because this information does not act directly upon DNA molecules like methylation. A broader term is needed to capture the diversity of three-dimensional spatial codes in biology. One such term is Heritable Nongenetic Information [HNI], which encompasses the TSC. The term heritable is appropriate because this information is passed onto offspring, otherwise it would have disappeared during evolution. Another reason for the heritability of HNI is that the spatial information observed in tissues is not reducible to the laws of physics, meaning structures like epithelial tubes or neural circuits do not spontaneously form in an aqueous solution. Pre-existing physiological information in the microenvironment is necessary.


title: A faster algorithm for counting the integer points number in $\Delta$-modular polyhedra (corrected version)

abstract: Let a polytope $P$ be defined by a system $A x \leq b$. We consider the problem of counting the number of integer points inside $P$, assuming that $P$ is $\Delta$-modular, where the polytope $P$ is called $\Delta$-modular if all the rank sub-determinants of $A$ are bounded by $\Delta$ in the absolute value. We present a new FPT-algorithm, parameterized by $\Delta$ and by the maximal number of vertices in $P$, where the maximum is taken by all r.h.s. vectors $b$. We show that our algorithm is more efficient for $\Delta$-modular problems than the approach of A. Barvinok et al. To this end, we do not directly compute the short rational generating function for $P \cap Z^n$, which is commonly used for the considered problem. Instead, we use the dynamic programming principle to compute its particular representation in the form of exponential series that depends on a single variable. We completely do not rely to the Barvinok's unimodular sign decomposition technique. Using our new complexity bound, we consider different special cases that may be of independent interest. For example, we give FPT-algorithms for counting the integer points number in $\Delta$-modular simplices and similar polytopes that have $n + O(1)$ facets. As a special case, for any fixed $m$, we give an FPT-algorithm to count solutions of the unbounded $m$-dimensional $\Delta$-modular subset-sum problem.


title: Dust dissipation timescales in the intermediate and outer regions of protoplanetary disks

abstract: Dust and gas in protoplanetary disks dissipate as central stars evolve. In order to estimate the dust dissipation timescales in the protoplanetary disks, we stacked the WISE 12, 22, and the AKARI 90 $\mu$m survey images of known T Tauri stars and derived the average fluxes, well below the survey flux limit in the 90 $\mu$m band. We classified 4,783 T Tauri stars into three age groups, which are young ($<$2 Myr), mid-age (2-6 Myr), and old ($>$6 Myr) groups, and stacked the WISE 12 and 22 and the AKARI 90 $\mu$m images in each group. The photometry of the stacked image shows the flux decay timescales of 1.4$\pm$0.2, 1.38$\pm$0.05, and 1.4$^{+0.6}_{-0.5}$ Myr in the 12, 22, and 90 $\mu$m bands, respectively. In optically thin disks with one-solar luminosity central stars, the 12 and 22 $\mu$m fluxes are attributed to the emission from the intermediate ($\sim$ 1 au) region and the 90 $\mu$m flux corresponds to that from the outer ($\sim$ 10 au) region in the disk. We hence conclude the dust dissipation timescale is $\tau_{\rm med,dust}\sim$1.4 Myr in the intermediate disks and is $\tau_{\rm outer, dust}=$1.4$^{+0.6}_{-0.5}$ Myr in the outer disks. The dust-dissipation time difference between the outer and intermediate disks is $\Delta\tau_{\rm dust} = \tau_{\rm outer,dust} -\tau_{\rm med,dust}=0.0^{+0.6}_{-0.5}$ Myr, indicating that the dust in the intermediate and outer disks dissipates on almost the same timescale.


title: Distributed Model Predictive Control of Buildings and Energy Hubs

abstract: Model predictive control (MPC) strategies can be applied to the coordination of energy hubs to reduce their energy consumption. Despite the effectiveness of these techniques, their potential for energy savings are potentially underutilized due to the fact that energy demands are often assumed to be fixed quantities rather than controlled dynamic variables. The joint optimization of energy hubs and buildings' energy management systems can result in higher energy savings. This paper investigates how different MPC strategies perform on energy management systems in buildings and energy hubs. We first discuss two MPC approaches; centralized and decentralized. While the centralized control strategy offers optimal performance, its implementation is computationally prohibitive and raises privacy concerns. On the other hand, the decentralized control approach, which offers ease of implementation, displays significantly lower performance. We propose a third strategy, distributed control based on dual decomposition, which has the advantages of both approaches. Numerical case studies and comparisons demonstrate that the performance of distributed control is close to the performance of the centralized case, while maintaining a significantly lower computational burden, especially in large-scale scenarios with many agents. Finally, we validate and verify the reliability of the proposed method through an experiment on a full-scale energy hub system in the NEST demonstrator in D\"{u}bendorf, Switzerland.


title: On partially hyperbolic diffeomorphisms in dimension three via a notion of autonomous dynamics

abstract: We introduce a notion of autonomous dynamical systems and apply it to prove rigidity of partially hyperbolic diffeomorphisms on closed compact three-manifolds under some smoothness hypothesis of their associated framing.


title: AdjointBackMapV2: Precise Reconstruction of Arbitrary CNN Unit's Activation via Adjoint Operators

abstract: Adjoint operators have been found to be effective in the exploration of CNN's inner workings [1]. However, the previous no-bias assumption restricted its generalization. We overcome the restriction via embedding input images into an extended normed space that includes bias in all CNN layers as part of the extended space and propose an adjoint-operator-based algorithm that maps high-level weights back to the extended input space for reconstructing an effective hypersurface. Such hypersurface can be computed for an arbitrary unit in the CNN, and we prove that this reconstructed hypersurface, when multiplied by the original input (through an inner product), will precisely replicate the output value of each unit. We show experimental results based on the CIFAR-10 and CIFAR-100 data sets where the proposed approach achieves near 0 activation value reconstruction error.


title: A Decomposition Algorithm for Large-Scale Security-Constrained AC Optimal Power Flow

abstract: A decomposition algorithm for solving large-scale security-constrained AC optimal power flow problems is presented. The formulation considered is the one used in the ARPA-E Grid Optimization (GO) Competition, Challenge 1, held from November 2018 through October 2019. The techniques found to be most effective in terms of performance in the challenge are presented, including strategies for contingency selection, fast contingency evaluation, handling complementarity constraints, avoiding issues related to degeneracy, and exploiting parallelism. The results of numerical experiments are provided to demonstrate the effectiveness of the proposed techniques as compared to alternative strategies.


title: Implementing commercial inverse design tools for compact, phase-encoded, plasmonic digital logic devices

abstract: Numerical simulations have become an essential design tool in the field of photonics, especially for nanophotonics. In particular, 3D finite-difference-time-domain (FDTD) simulations are popular for their powerful design capabilities. Increasingly, researchers are developing or using inverse design tools to improve device footprints and performance. These tools often make use of 3D FDTD simulations and the adjoint optimization method. In this work, we implement a commercial inverse design tool with these features for several plasmonic devices that push the boundaries of the tool. We design a logic gate with complex design requirements, as well as a y-splitter and waveguide crossing. With minimal code changes, we implement a design that incorporates phase encoded inputs in a dielectric-loaded surface plasmon polariton waveguide. The complexity of the requirements in conjunction with limitations in the inverse design tool force us to make concessions regarding the density of encoding and to use on-off keying to encode the outputs. We compare the performance of the inverse-designed devices to conventionally-designed devices with the same operational behavior. Finally, we discuss the limitations of the inverse design tools for realizing complex device designs and comment on what is possible at present and where improvements can be made.


title: Two-loop gauge coupling $\beta$-function in a four-dimensional framework: the Standard Model case

abstract: At present, the gauge coupling $\beta$-function in the Standard Model (SM) is known up to four-loop order. As most SM calculations, dimensional regularization was employed. Despite its striking success, other regularization schemes have emerged, which aim to stay in the physical dimension as far as possible. In this contribution, we apply Implicit Regularization, a scheme defined in momentum space that complies with BPHZ, to obtain the two-loop gauge coupling $\beta$-function in the Standard Model. We reproduce the same result obtained when applying dimensional regularization.


title: Error Correction for FrodoKEM Using the Gosset Lattice

abstract: We consider FrodoKEM, a lattice-based cryptosystem based on LWE, and propose a new error correction mechanism to improve its performance. Our encoder maps the secret key block-wise into the Gosset lattice $E_8$. We propose two sets of parameters for our modified implementation. Thanks to the improved error correction, the first implementation outperforms FrodoKEM in terms of concrete security by $10$ to $13$ bits by increasing the error variance; the second allows to reduce the bandwidth by $7\%$ by halving the modulus $q$. In both cases, the decryption failure probability is improved compared to the original FrodoKEM. Unlike some previous works on error correction for lattice-based protocols, we provide a rigorous error probability bound by decomposing the error matrix into blocks with independent error coefficients.


title: Beware the Gini Index! A New Inequality Measure

abstract: The Gini index underestimates inequality for heavy-tailed distributions: for example, a Pareto distribution with exponent 1.5 (which has infinite variance) has the same Gini index as any exponential distribution (a mere 0.5). This is because the Gini index is relatively robust to extreme observations: while a statistic's robustness to extremes is desirable for data potentially distorted by outliers, it is misleading for heavy-tailed distributions, which inherently exhibit extremes. We propose an alternative inequality index: the variance normalized by the second moment. This ratio is more stable (hence more reliable) for large samples from an infinite-variance distribution than the Gini index paradoxically. Moreover, the new index satisfies the normative axioms of inequality measurement; notably, it is decomposable into inequality within and between subgroups, unlike the Gini index.


title: Epileptic Seizure Classification Using Combined Labels and a Genetic Algorithm

abstract: Epilepsy affects 50 million people worldwide and is one of the most common serious neurological disorders. Seizure detection and classification is a valuable tool for diagnosing and maintaining the condition. An automated classification algorithm will allow for accurate diagnosis. Utilising the Temple University Hospital (TUH) Seizure Corpus, six seizure types are compared; absence, complex partial, myoclonic, simple partial, tonic and tonic- clonic models. This study proposes a method that utilises unique features with a novel parallel classifier - Parallel Genetic Naive Bayes (NB) Seizure Classifier (PGNBSC). The PGNBSC algorithm searches through the features and by reclassifying the data each time, the algorithm will create a matrix for optimum search criteria. Ictal states from the EEGs are segmented into 1.8 s windows, where the epochs are then further decomposed into 13 different features from the first intrinsic mode function (IMF). The features are compared using an original NB classifier in the first model. This is improved upon in a second model by using a genetic algorithm (Binary Grey Wolf Optimisation, Option 1) with a NB classifier. The third model uses a combination of the simple partial and complex partial seizures to provide the highest classification accuracy for each of the six seizures amongst the three models (20%, 53%, and 85% for first, second, and third model, respectively).


title: Design and Characterization of a 3D-printed Pneumatically-driven Bistable Valve with Tunable Characteristics

abstract: Although research studies in pneumatic soft robots develop rapidly, most pneumatic actuators are still controlled by rigid valves and conventional electronics. The existence of these rigid, electronic components sacrifices the compliance and adaptability of soft robots.} Current electronics-free valve designs based on soft materials are facing challenges in behaviour consistency, design flexibility, and fabrication complexity. Taking advantages of soft material 3D printing, this paper presents a new design of a bi-stable pneumatic valve, which utilises two soft, pneumatically-driven, and symmetrically-oriented conical shells with structural bistability to stabilise and regulate the airflow. The critical pressure required to operate the valve can be adjusted by changing the design features of the soft bi-stable structure. Multi-material printing simplifies the valve fabrication, enhances the flexibility in design feature optimisations, and improves the system repeatability. In this work, both a theoretical model and physical experiments are introduced to examine the relationships between the critical operating pressure and the key design features. Results with valve characteristic tuning via material stiffness changing show better effectiveness compared to the change of geometry design features (demonstrated largest tunable critical pressure range from 15.3 to 65.2 kPa and fastest response time $\leq$ 1.8 s.