title: A new non-negative distribution with both finite and infinite support

abstract: The Tukey-$\lambda$ distribution has interesting properties including (i) for some parameters values it has finite support, and for others infinite support, and (ii) it can mimic several other distributions such that parameter estimation for the Tukey distribution is a method for identifying an appropriate class of distribution to model a set of data. The Tukey-$\lambda$ is, however, symmetric. Here we define a new class of {\em non-negative} distribution with similar properties to the Tukey-$\lambda$ distribution. As with the Tukey-$\lambda$ distribution, our distribution is defined in terms of its quantile function, which in this case is given by the polylogarithm function. We show the support of the distribution to be the Riemann zeta function (when finite), and we provide a closed form for the expectation, provide simple means to calculate the CDF and PDF, and show that it has relationships to the uniform, exponential, inverse beta and extreme-value distributions.


title: Uncertainty-Aware Deep Calibrated Salient Object Detection

abstract: Existing deep neural network based salient object detection (SOD) methods mainly focus on pursuing high network accuracy. However, those methods overlook the gap between network accuracy and prediction confidence, known as the confidence uncalibration problem. Thus, state-of-the-art SOD networks are prone to be overconfident. In other words, the predicted confidence of the networks does not reflect the real probability of correctness of salient object detection, which significantly hinder their real-world applicability. In this paper, we introduce an uncertaintyaware deep SOD network, and propose two strategies from different perspectives to prevent deep SOD networks from being overconfident. The first strategy, namely Boundary Distribution Smoothing (BDS), generates continuous labels by smoothing the original binary ground-truth with respect to pixel-wise uncertainty. The second strategy, namely Uncertainty-Aware Temperature Scaling (UATS), exploits a relaxed Sigmoid function during both training and testing with spatially-variant temperature scaling to produce softened output. Both strategies can be incorporated into existing deep SOD networks with minimal efforts. Moreover, we propose a new saliency evaluation metric, namely dense calibration measure C, to measure how the model is calibrated on a given dataset. Extensive experimental results on seven benchmark datasets demonstrate that our solutions can not only better calibrate SOD models, but also improve the network accuracy.


title: Descriptive and Predictive Analysis of Aggregating Functions in Serverless Clouds: the Case of Video Streaming

abstract: Serverless clouds allocate multiple tasks (e.g., micro-services) from multiple users on a shared pool of computing resources. This enables serverless cloud providers to reduce their resource usage by transparently aggregate similar tasks of a certain context (e.g., video processing) that share the whole or part of their computation. To this end, it is crucial to know the amount of time-saving achieved by aggregating the tasks. Lack of such knowledge can lead to uninformed merging and scheduling decisions that, in turn, can cause deadline violation of either the merged tasks or other following tasks. Accordingly, in this paper, we study the problem of estimating execution-time saving resulted from merging tasks with the example in the context of video processing. To learn the execution-time saving in different forms of merging, we first establish a set of benchmarking videos and examine a wide variety of video processing tasks -- with and without merging in place. We observed that although merging can save up to 44% in the execution-time, the number of possible merging cases is intractable. Hence, in the second part, we leverage the benchmarking results and develop a method based on Gradient Boosting Decision Tree (GBDT) to estimate the time-saving for any given task merging case. Experimental results show that the method can estimate the time-saving with the error rate of 0.04, measured based on Root Mean Square Error (RMSE).


title: Motion and Force Planning for Manipulating Heavy Objects by Pivoting

abstract: Manipulation of objects by exploiting their contact with the environment can enhance both the dexterity and payload capability of robotic manipulators. A common way to manipulate heavy objects beyond the payload capability of a robot is to use a sequence of pivoting motions, wherein, an object is moved while some contact points between the object and a support surface are kept fixed. The goal of this paper is to develop an algorithmic approach for automated plan generation for object manipulation with a sequence of pivoting motions. A plan for manipulating a heavy object consists of a sequence of joint angles of the manipulator, the corresponding object poses, as well as the joint torques required to move the object. The constraint of maintaining object contact with the ground during manipulation results in nonlinear constraints in the configuration space of the robot, which is challenging for motion planning algorithms. Exploiting the fact that pivoting motion corresponds to movements in a subgroup of the group of rigid body motions, SE(3), we present a novel task-space based planning approach for computing a motion plan for both the manipulator and the object while satisfying contact constraints. We also combine our motion planning algorithm with a grasping force synthesis algorithm to ensure that friction constraints at the contacts and actuator torque constraints are satisfied. We present simulation results with a dual-armed Baxter robot to demonstrate our approach.


title: Cost-to-Go Function Generating Networks for High Dimensional Motion Planning

abstract: This paper presents c2g-HOF networks which learn to generate cost-to-go functions for manipulator motion planning. The c2g-HOF architecture consists of a cost-to-go function over the configuration space represented as a neural network (c2g-network) as well as a Higher Order Function (HOF) network which outputs the weights of the c2g-network for a given input workspace. Both networks are trained end-to-end in a supervised fashion using costs computed from traditional motion planners. Once trained, c2g-HOF can generate a smooth and continuous cost-to-go function directly from workspace sensor inputs (represented as a point cloud in 3D or an image in 2D). At inference time, the weights of the c2g-network are computed very efficiently and near-optimal trajectories are generated by simply following the gradient of the cost-to-go function. We compare c2g-HOF with traditional planning algorithms for various robots and planning scenarios. The experimental results indicate that planning with c2g-HOF is significantly faster than other motion planning algorithms, resulting in orders of magnitude improvement when including collision checking. Furthermore, despite being trained from sparsely sampled trajectories in configuration space, c2g-HOF generalizes to generate smoother, and often lower cost, trajectories. We demonstrate cost-to-go based planning on a 7 DoF manipulator arm where motion planning in a complex workspace requires only 0.13 seconds for the entire trajectory.


title: Robustness and Transferability of Universal Attacks on Compressed Models

abstract: Neural network compression methods like pruning and quantization are very effective at efficiently deploying Deep Neural Networks (DNNs) on edge devices. However, DNNs remain vulnerable to adversarial examples-inconspicuous inputs that are specifically designed to fool these models. In particular, Universal Adversarial Perturbations (UAPs), are a powerful class of adversarial attacks which create adversarial perturbations that can generalize across a large set of inputs. In this work, we analyze the effect of various compression techniques to UAP attacks, including different forms of pruning and quantization. We test the robustness of compressed models to white-box and transfer attacks, comparing them with their uncompressed counterparts on CIFAR-10 and SVHN datasets. Our evaluations reveal clear differences between pruning methods, including Soft Filter and Post-training Pruning. We observe that UAP transfer attacks between pruned and full models are limited, suggesting that the systemic vulnerabilities across these models are different. This finding has practical implications as using different compression techniques can blunt the effectiveness of black-box transfer attacks. We show that, in some scenarios, quantization can produce gradient-masking, giving a false sense of security. Finally, our results suggest that conclusions about the robustness of compressed models to UAP attacks is application dependent, observing different phenomena in the two datasets used in our experiments.


title: Exploring Deep Neural Networks and Transfer Learning for Analyzing Emotions in Tweets

abstract: In this paper, we present an experiment on using deep learning and transfer learning techniques for emotion analysis in tweets and suggest a method to interpret our deep learning models. The proposed approach for emotion analysis combines a Long Short Term Memory (LSTM) network with a Convolutional Neural Network (CNN). Then we extend this approach for emotion intensity prediction using transfer learning technique. Furthermore, we propose a technique to visualize the importance of each word in a tweet to get a better understanding of the model. Experimentally, we show in our analysis that the proposed models outperform the state-of-the-art in emotion classification while maintaining competitive results in predicting emotion intensity.


title: Superabsorption in an organic microcavity: towards a quantum battery

abstract: The rate at which matter emits or absorbs light can be modified by its environment, as dramatically exemplified by the widely-studied phenomenon of superradiance. The reverse process, superabsorption, is harder to demonstrate due to the challenges of probing ultrafast processes, and has only been seen for small numbers of atoms. Its central idea - superextensive scaling of absorption meaning larger systems absorb faster - is also the key idea underpinning quantum batteries. Here we implement experimentally a paradigmatic model of a quantum battery, constructed of a microcavity enclosing a molecular dye. Ultrafast optical spectroscopy allows us to observe charging dynamics at femtosecond resolution to demonstrate superextensive charging rates and storage capacity, in agreement with our theoretical modelling. We find that decoherence plays an important role in stabilising energy storage. Our work opens new opportunities for harnessing collective effects in light-matter coupling for nanoscale energy capture, storage, and transport technologies.


title: Comparison of Update and Genetic Training Algorithms in a Memristor Crossbar Perceptron

abstract: Memristor-based computer architectures are becoming more attractive as a possible choice of hardware for the implementation of neural networks. However, at present, memristor technologies are susceptible to a variety of failure modes, a serious concern in any application where regular access to the hardware may not be expected or even possible. In this study, we investigate whether certain training algorithms may be more resilient to particular hardware failure modes, and therefore more suitable for use in those applications. We implement two training algorithms -- a local update scheme and a genetic algorithm -- in a simulated memristor crossbar, and compare their ability to train for a simple image classification task as an increasing number of memristors fail to adjust their conductance. We demonstrate that there is a clear distinction between the two algorithms in several measures of the rate of failure to train.


title: Exactly unitary discrete representations of the metaplectic transform for linear-time algorithms

abstract: The metaplectic transform (MT), a generalization of the Fourier transform sometimes called the linear canonical transform, is a tool used ubiquitously in modern optics, for example, when calculating the transformations of light beams in paraxial optical systems. The MT is also an essential ingredient of the geometrical-optics modeling of caustics that was recently proposed by the authors. In particular, this application relies on the near-identity MT (NIMT); however, the NIMT approximation used so far is not exactly unitary and leads to numerical instability. Here, we develop a discrete MT that is exactly unitary, and approximate it to obtain a discrete NIMT that is also unitary and can be computed in linear time. We prove that the discrete NIMT converges to the discrete MT when iterated, thereby allowing the NIMT to compute MTs that are not necessarily near-identity. We then demonstrate the new algorithms with a series of examples.


title: Correlated Charge Noise and Relaxation Errors in Superconducting Qubits

abstract: The central challenge in building a quantum computer is error correction. Unlike classical bits, which are susceptible to only one type of error, quantum bits ("qubits") are susceptible to two types of error, corresponding to flips of the qubit state about the $X$- and $Z$-directions. While the Heisenberg Uncertainty Principle precludes simultaneous monitoring of $X$- and $Z$-flips on a single qubit, it is possible to encode quantum information in large arrays of entangled qubits that enable accurate monitoring of all errors in the system, provided the error rate is low. Another crucial requirement is that errors cannot be correlated. Here, we characterize a superconducting multiqubit circuit and find that charge fluctuations are highly correlated on a length scale over 600~$\mu$m; moreover, discrete charge jumps are accompanied by a strong transient suppression of qubit energy relaxation time across the millimeter-scale chip. The resulting correlated errors are explained in terms of the charging event and phonon-mediated quasiparticle poisoning associated with absorption of gamma rays and cosmic-ray muons in the qubit substrate. Robust quantum error correction will require the development of mitigation strategies to protect multiqubit arrays from correlated errors due to particle impacts.


title: A family of concordance homomorphisms from Khovanov homology

abstract: By considering a version of Khovanov homology incorporating both the Lee and $E(-1)$ differentials, we construct a $1$-parameter family of concordance homomorphisms similar to the Upsilon invariant from knot Floer homology. This invariant gives lower bounds on the slice genus and can be used to prove that certain infinite families of pretzel knots are linearly independent in the smooth concordance group.


title: Reinforcement Learning Agents for Ubisoft's Roller Champions

abstract: In recent years, Reinforcement Learning (RL) has seen increasing popularity in research and popular culture. However, skepticism still surrounds the practicality of RL in modern video game development. In this paper, we demonstrate by example that RL can be a great tool for Artificial Intelligence (AI) design in modern, non-trivial video games. We present our RL system for Ubisoft's Roller Champions, a 3v3 Competitive Multiplayer Sports Game played on an oval-shaped skating arena. Our system is designed to keep up with agile, fast-paced development, taking 1--4 days to train a new model following gameplay changes. The AIs are adapted for various game modes, including a 2v2 mode, a Training with Bots mode, in addition to the Classic game mode where they replace players who have disconnected. We observe that the AIs develop sophisticated co-ordinated strategies, and can aid in balancing the game as an added bonus. Please see the accompanying video at https://vimeo.com/466780171 (password: rollerRWRL2020) for examples.


title: Tensor product of modules over a Lie conformal algebra

abstract: We find a necessary and sufficient condition for the existence of the tensor product of modules over a Lie conformal algebra. We provide two algebraic constructions of the tensor product. We show the relation between tensor product and conformal linear maps. We prove commutativity of the tensor product.


title: Autocatalytic systems and recombination: a reaction network perspective

abstract: Autocatalytic systems are very often incorporated in the "origin of life" models, a connection that has been analyzed in the context of the classical hypercycles introduced by Manfred Eigen. We investigate the dynamics of certain networks called bimolecular autocatalytic systems. In particular, we consider the dynamics corresponding to the relative populations in these networks, and show that they can be analyzed by studying well-chosen autonomous polynomial dynamical systems. Moreover, we find that one can use results from reaction network theory to prove persistence and permanence of several types of bimolecular autocatalytic systems called autocatalytic recombination networks.


title: Artificial Intelligence & Cooperation

abstract: The rise of Artificial Intelligence (AI) will bring with it an ever-increasing willingness to cede decision-making to machines. But rather than just giving machines the power to make decisions that affect us, we need ways to work cooperatively with AI systems. There is a vital need for research in "AI and Cooperation" that seeks to understand the ways in which systems of AIs and systems of AIs with people can engender cooperative behavior. Trust in AI is also key: trust that is intrinsic and trust that can only be earned over time. Here we use the term "AI" in its broadest sense, as employed by the recent 20-Year Community Roadmap for AI Research (Gil and Selman, 2019), including but certainly not limited to, recent advances in deep learning. With success, cooperation between humans and AIs can build society just as human-human cooperation has. Whether coming from an intrinsic willingness to be helpful, or driven through self-interest, human societies have grown strong and the human species has found success through cooperation. We cooperate "in the small" -- as family units, with neighbors, with co-workers, with strangers -- and "in the large" as a global community that seeks cooperative outcomes around questions of commerce, climate change, and disarmament. Cooperation has evolved in nature also, in cells and among animals. While many cases involving cooperation between humans and AIs will be asymmetric, with the human ultimately in control, AI systems are growing so complex that, even today, it is impossible for the human to fully comprehend their reasoning, recommendations, and actions when functioning simply as passive observers.


title: SensiX: A Platform for Collaborative Machine Learning on the Edge

abstract: The emergence of multiple sensory devices on or near a human body is uncovering new dynamics of extreme edge computing. In this, a powerful and resource-rich edge device such as a smartphone or a Wi-Fi gateway is transformed into a personal edge, collaborating with multiple devices to offer remarkable sensory al eapplications, while harnessing the power of locality, availability, and proximity. Naturally, this transformation pushes us to rethink how to construct accurate, robust, and efficient sensory systems at personal edge. For instance, how do we build a reliable activity tracker with multiple on-body IMU-equipped devices? While the accuracy of sensing models is improving, their runtime performance still suffers, especially under this emerging multi-device, personal edge environments. Two prime caveats that impact their performance are device and data variabilities, contributed by several runtime factors, including device availability, data quality, and device placement. To this end, we present SensiX, a personal edge platform that stays between sensor data and sensing models, and ensures best-effort inference under any condition while coping with device and data variabilities without demanding model engineering. SensiX externalises model execution away from applications, and comprises of two essential functions, a translation operator for principled mapping of device-to-device data and a quality-aware selection operator to systematically choose the right execution path as a function of model accuracy. We report the design and implementation of SensiX and demonstrate its efficacy in developing motion and audio-based multi-device sensing systems. Our evaluation shows that SensiX offers a 7-13% increase in overall accuracy and up to 30% increase across different environment dynamics at the expense of 3mW power overhead.


title: Data-based Discovery of Governing Equations

abstract: Most common mechanistic models are traditionally presented in mathematical forms to explain a given physical phenomenon. Machine learning algorithms, on the other hand, provide a mechanism to map the input data to output without explicitly describing the underlying physical process that generated the data. We propose a Data-based Physics Discovery (DPD) framework for automatic discovery of governing equations from observed data. Without a prior definition of the model structure, first a free-form of the equation is discovered, and then calibrated and validated against the available data. In addition to the observed data, the DPD framework can utilize available prior physical models, and domain expert feedback. When prior models are available, the DPD framework can discover an additive or multiplicative correction term represented symbolically. The correction term can be a function of the existing input variable to the prior model, or a newly introduced variable. In case a prior model is not available, the DPD framework discovers a new data-based standalone model governing the observations. We demonstrate the performance of the proposed framework on a real-world application in the aerospace industry.


title: Photonic realization of the kappa-deformed Dirac equation

abstract: We show an implementation of a kappa-deformed Dirac equation in tight-binding arrays of photonic waveguides. This is done with a special configuration of couplings extending to second nearest neighbors. Geometric manipulations can control these evanescent couplings. A careful study of wave packet propagation is presented, including the effects of deformation parameters on Zitterbewegung or trembling motion. In this way, we demonstrate how to recreate the effects of a flat noncommutative spacetime -i.e., kappa-Minkowski spacetime -in simple experimental setups. We touch upon elastic realizations in the section of Conclusions.


title: Excess Power, Energy and Intensity of Stochastic Fields in Quasi-Static and Dynamic Environments

abstract: The excess power, energy and intensity of a random electromagnetic field above a high threshold level are characterized based on a Slepian--Kac model for upcrossings. For quasi-static fields, the probability distribution of the excess intensity in its regression approximation evolves from $\chi^2_3$ to $\chi^2_2$ when the threshold level increases. The excursion area associated with excess energy exhibits a chi-cubed ($\chi^3_2$) distribution above asymptotically high thresholds, where excursions are parabolic. For dynamic fields, the dependence of the electrical and environmental modulations of the excess power on the hybrid modulation index and threshold level are established. The normalized effective power relative to the quasi-static power increases non-monotonically when this index increases. The mean and standard deviation of the dynamic excess power are obtained in closed form and validated by Monte Carlo simulation.

categories: physics.app-ph physics.class-ph physics.data-an stat.AP


title: Anisotropies in Compressible MHD Turbulence: Probing Magnetic Fields and Measuring Magnetization

abstract: Probing magnetic fields in the interstellar medium (ISM) is notoriously challenging. Motivated by the modern theories of magnetohydrodynamic (MHD) turbulence and turbulence anisotropy, we introduce the Structure-Function Analysis (SFA) as a new approach to measure the magnetic field orientation and estimate the magnetization. We analyze the statistics of turbulent velocities in three-dimensional compressible MHD simulations through the second-order structure functions in both local and global reference frames. In the sub-Alfvenic turbulence with the magnetic energy larger than the turbulent energy, the SFA of turbulent velocities measured in the directions perpendicular and parallel to the magnetic field can be significantly different. Their ratio has a power-law dependence on the Alfven Mach number $M_A$, which is inversely proportional to the magnetic field strength. We demonstrate that the anisotropic structure functions of turbulent velocities can be used to estimate both the orientation and strength of magnetic fields. With turbulent velocities measured using different tracers, our approach can be generally applied to probing the magnetic fields in the multi-phase interstellar medium.


title: Data-Driven System Identification of Linear Quantum Systems Coupled to Time-Varying Coherent Inputs

abstract: In this paper, we develop a system identification algorithm to identify a model for unknown linear quantum systems driven by time-varying coherent states, based on empirical single-shot continuous homodyne measurement data of the system's output. The proposed algorithm identifies a model that satisfies the physical realizability conditions for linear quantum systems, challenging constraints not encountered in classical (non-quantum) linear system identification. Numerical examples on a multiple-input multiple-output optical cavity model are presented to illustrate an application of the identification algorithm.


title: Recursive Generation of The Semi-Classical Expansion in Arbitrary Dimension

abstract: We present a recursive procedure, which is based on the small time expansion of the propagator, in order to generate a semi-classical expansion of the \textit{quantum action} for a quantum mechanical potential in arbitrary dimensions. In the method we use the spectral information emerges from the singularities of the propagator on the complex $t$ plane, which are handled by the $i\ve$ prescription and basic complex analysis. This feature allows for generalization to higher dimensions. We illustrate the procedure by providing simple examples in non-relativistic quantum mechanics.


title: NEO Population, Velocity Bias, and Impact Risk from an ATLAS Analysis

abstract: We estimate the total population of near-Earth objects (NEOs) in the Solar System, using an extensive, `Solar System to pixels' fake-asteroid simulation to debias detections of real NEOs by the ATLAS survey. Down to absolute magnitudes $H=25$ and 27.6 (diameters of $\sim 34$ and 10 meters, respectively, for 15% albedo), we find total populations of $(3.72 \pm 0.49) \times 10^5$ and $(1.59 \pm 0.45) \times 10^7$ NEOs, respectively. Most plausible sources of error tend toward underestimation, so the true populations are likely larger. We find the distribution of $H$ magnitudes steepens for NEOs fainter than $H \sim 22.5$, making small asteroids more common than extrapolation from brighter $H$ mags would predict. Our simulation indicates a strong bias against detecting small but dangerous asteroids that encounter Earth with high relative velocities -- i.e., asteroids in highly inclined and/or eccentric orbits. Worldwide NEO discovery statistics indicate this bias affects global NEO detection capability, to the point that an observational census of small asteroids in such orbits is probably not currently feasible. Prompt and aggressive followup of NEO candidates, combined with closer collaborations between segments of the global NEO community, can increase detection rates for these dangerous objects.


title: Provable Defense against Privacy Leakage in Federated Learning from Representation Perspective

abstract: Federated learning (FL) is a popular distributed learning framework that can reduce privacy risks by not explicitly sharing private data. However, recent works demonstrated that sharing model updates makes FL vulnerable to inference attacks. In this work, we show our key observation that the data representation leakage from gradients is the essential cause of privacy leakage in FL. We also provide an analysis of this observation to explain how the data presentation is leaked. Based on this observation, we propose a defense against model inversion attack in FL. The key idea of our defense is learning to perturb data representation such that the quality of the reconstructed data is severely degraded, while FL performance is maintained. In addition, we derive certified robustness guarantee to FL and convergence guarantee to FedAvg, after applying our defense. To evaluate our defense, we conduct experiments on MNIST and CIFAR10 for defending against the DLG attack and GS attack. Without sacrificing accuracy, the results demonstrate that our proposed defense can increase the mean squared error between the reconstructed data and the raw data by as much as more than 160X for both DLG attack and GS attack, compared with baseline defense methods. The privacy of the FL system is significantly improved.


title: Mesoscopic photogrammetry with an unstabilized phone camera

abstract: We present a feature-free photogrammetric technique that enables quantitative 3D mesoscopic (mm-scale height variation) imaging with tens-of-micron accuracy from sequences of images acquired by a smartphone at close range (several cm) under freehand motion without additional hardware. Our end-to-end, pixel-intensity-based approach jointly registers and stitches all the images by estimating a coaligned height map, which acts as a pixel-wise radial deformation field that orthorectifies each camera image to allow homographic registration. The height maps themselves are reparameterized as the output of an untrained encoder-decoder convolutional neural network (CNN) with the raw camera images as the input, which effectively removes many reconstruction artifacts. Our method also jointly estimates both the camera's dynamic 6D pose and its distortion using a nonparametric model, the latter of which is especially important in mesoscopic applications when using cameras not designed for imaging at short working distances, such as smartphone cameras. We also propose strategies for reducing computation time and memory, applicable to other multi-frame registration problems. Finally, we demonstrate our method using sequences of multi-megapixel images captured by an unstabilized smartphone on a variety of samples (e.g., painting brushstrokes, circuit board, seeds).


abstract: A bijection $(l,r)$ between $M^2$ and $M$ is said to be a pairing function with no cycles, if any composition of its coordinate functions has no fixed point. We compute here the Grothendieck ring of the pairing function without cycles to be isomorphic to $\mathbb{Z}^2\simeq \mathbb{Z}[X]/(X-X^2)$. More generally, for any $n\in \mathbb{N}^*$ and any bijetion without cycles betwen $M$ and $M^n$, the exact same method proves that $K_0(M)=\mathbb{Z}[X]/(X-X^n)$.


title: Interactive Weak Supervision: Learning Useful Heuristics for Data Labeling

abstract: Obtaining large annotated datasets is critical for training successful machine learning models and it is often a bottleneck in practice. Weak supervision offers a promising alternative for producing labeled datasets without ground truth annotations by generating probabilistic labels using multiple noisy heuristics. This process can scale to large datasets and has demonstrated state of the art performance in diverse domains such as healthcare and e-commerce. One practical issue with learning from user-generated heuristics is that their creation requires creativity, foresight, and domain expertise from those who hand-craft them, a process which can be tedious and subjective. We develop the first framework for interactive weak supervision in which a method proposes heuristics and learns from user feedback given on each proposed heuristic. Our experiments demonstrate that only a small number of feedback iterations are needed to train models that achieve highly competitive test set performance without access to ground truth training labels. We conduct user studies, which show that users are able to effectively provide feedback on heuristics and that test set results track the performance of simulated oracles.


title: KNN Classification with One-step Computation

abstract: KNN classification is an improvisational learning mode, in which they are carried out only when a test data is predicted that set a suitable K value and search the K nearest neighbors from the whole training sample space, referred them to the lazy part of KNN classification. This lazy part has been the bottleneck problem of applying KNN classification due to the complete search of K nearest neighbors. In this paper, a one-step computation is proposed to replace the lazy part of KNN classification. The one-step computation actually transforms the lazy part to a matrix computation as follows. Given a test data, training samples are first applied to fit the test data with the least squares loss function. And then, a relationship matrix is generated by weighting all training samples according to their influence on the test data. Finally, a group lasso is employed to perform sparse learning of the relationship matrix. In this way, setting K value and searching K nearest neighbors are both integrated to a unified computation. In addition, a new classification rule is proposed for improving the performance of one-step KNN classification. The proposed approach is experimentally evaluated, and demonstrated that the one-step KNN classification is efficient and promising


title: Reinforced Multi-Teacher Selection for Knowledge Distillation

abstract: In natural language processing (NLP) tasks, slow inference speed and huge footprints in GPU usage remain the bottleneck of applying pre-trained deep models in production. As a popular method for model compression, knowledge distillation transfers knowledge from one or multiple large (teacher) models to a small (student) model. When multiple teacher models are available in distillation, the state-of-the-art methods assign a fixed weight to a teacher model in the whole distillation. Furthermore, most of the existing methods allocate an equal weight to every teacher model. In this paper, we observe that, due to the complexity of training examples and the differences in student model capability, learning differentially from teacher models can lead to better performance of student models distilled. We systematically develop a reinforced method to dynamically assign weights to teacher models for different training instances and optimize the performance of student model. Our extensive experimental results on several NLP tasks clearly verify the feasibility and effectiveness of our approach.


title: The Rise of AI-Driven Simulators: Building a New Crystal Ball

abstract: The use of computational simulation is by now so pervasive in society that it is no exaggeration to say that continued U.S. and international prosperity, security, and health depend in part on continued improvements in simulation capabilities. What if we could predict weather two weeks out, guide the design of new drugs for new viral diseases, or manage new manufacturing processes that cut production costs and times by an order of magnitude? What if we could predict collective human behavior, for example, response to an evacuation request during a natural disaster, or labor response to fiscal stimulus? (See also the companion CCC Quad Paper on Pandemic Informatics, which discusses features that would be essential to solving large-scale problems like preparation for, and response to, the inevitable next pandemic.) The past decade has brought remarkable advances in complementary areas: in sensors, which can now capture enormous amounts of data about the world, and in AI methods capable of learning to extract predictive patterns from those data. These advances may lead to a new era in computational simulation, in which sensors of many kinds are used to produce vast quantities of data, AI methods identify patterns in those data, and new AI-driven simulators combine machine-learned and mathematical rules to make accurate and actionable predictions. At the same time, there are new challenges -- computers in some important regards are no longer getting faster, and in some areas we are reaching the limits of mathematical understanding, or at least of our ability to translate mathematical understanding into efficient simulation. In this paper, we lay out some themes that we envision forming part of a cohesive, multi-disciplinary, and application-inspired research agenda on AI-driven simulators.


title: Simultaneous gain profile design and noise figure prediction for Raman amplifiers using machine learning

abstract: A machine learning framework predicting pump powers and noise figure profile for a target distributed Raman amplifier gain profile is experimentally demonstrated. We employ a single-layer neural network to learn the mapping from the gain profiles to the pump powers and noise figures. The obtained results show highly-accurate gain profile designs and noise figure predictions, with a maximum error on average of ~0.3dB. This framework provides the comprehensive characterization of the Raman amplifier and thus is a valuable tool for predicting the performance of the next-generation optical communication systems, expected to employ Raman amplification.


title: Mapping the Space of Chemical Reactions Using Attention-Based Neural Networks

abstract: Organic reactions are usually assigned to classes containing reactions with similar reagents and mechanisms. Reaction classes facilitate the communication of complex concepts and efficient navigation through chemical reaction space. However, the classification process is a tedious task. It requires the identification of the corresponding reaction class template via annotation of the number of molecules in the reactions, the reaction center, and the distinction between reactants and reagents. This work shows that transformer-based models can infer reaction classes from non-annotated, simple text-based representations of chemical reactions. Our best model reaches a classification accuracy of 98.2%. We also show that the learned representations can be used as reaction fingerprints that capture fine-grained differences between reaction classes better than traditional reaction fingerprints. The insights into chemical reaction space enabled by our learned fingerprints are illustrated by an interactive reaction atlas providing visual clustering and similarity searching.


title: Interactive Search Based on Deep Reinforcement Learning

abstract: With the continuous development of machine learning technology, major e-commerce platforms have launched recommendation systems based on it to serve a large number of customers with different needs more efficiently. Compared with traditional supervised learning, reinforcement learning can better capture the user's state transition in the decision-making process, and consider a series of user actions, not just the static characteristics of the user at a certain moment. In theory, it will have a long-term perspective, producing a more effective recommendation. The special requirements of reinforcement learning for data make it need to rely on an offline virtual system for training. Our project mainly establishes a virtual user environment for offline training. At the same time, we tried to improve a reinforcement learning algorithm based on bi-clustering to expand the action space and recommended path space of the recommendation agent.


title: The semiclassical theory for spin dynamics in a disordered system

abstract: We investigate the Drude model of spin dynamics in two-dimensional spin-orbit coupled systems. In the absence of an applied electric field, the spin aligns with the k-dependent effective magnetic field. The influence of disorder (the momentum relaxation time $\tau$) on the system is considered. In the presence of an electric field, the change in momentum causes a change in the effective magnetic field. The in-plane spin can precess around the successive change in the orientation of the effective magnetic field. We find that up to the linear order of the electric field, the spin orientation undergoes Larmor-like and non-Larmor like precession. Furthermore, we find that the non-Larmor motion over a very short time ($t\ll\tau$) exactly equals the result obtained from the Kubo formula. This means that the Kubo formula only captures the system's response over a very short evolution. The spin-Hall conductivity for Larmor and non-Larmor precession in the Rashba system is analytically calculated. We find that the intrinsic spin-Hall current is not a universal constant and correctly drops to zero when the Rashba spin-orbit coupling drops to zero. We also calculate the time-averaged Larmor and non-Larmor spin-Hall conductivities (SHCs) for the k-cubic Rashba system and compare them to experimental values. The time-averaged Larmor SHC vanishes and the non-Larmor SHC is given by $2.1(q/8\pi)$ which is very close to the experimental value $2.2(q/8\pi)$ by Wunderlich et al. [Phys. Rev. Lett. {\bf 94}, 047204 (2005)].


title: Analyzing the Performance of Smart Industry 4.0 Applications on Cloud Computing Systems

abstract: Cloud-based Deep Neural Network (DNN) applications that make latency-sensitive inference are becoming an indispensable part of Industry 4.0. Due to the multi-tenancy and resource heterogeneity, both inherent to the cloud computing environments, the inference time of DNN-based applications are stochastic. Such stochasticity, if not captured, can potentially lead to low Quality of Service (QoS) or even a disaster in critical sectors, such as Oil and Gas industry. To make Industry 4.0 robust, solution architects and researchers need to understand the behavior of DNN-based applications and capture the stochasticity exists in their inference times. Accordingly, in this study, we provide a descriptive analysis of the inference time from two perspectives. First, we perform an application-centric analysis and statistically model the execution time of four categorically different DNN applications on both Amazon and Chameleon clouds. Second, we take a resource-centric approach and analyze a rate-based metric in form of Million Instruction Per Second (MIPS) for heterogeneous machines in the cloud. This non-parametric modeling, achieved via Jackknife and Bootstrap re-sampling methods, provides the confidence interval of MIPS for heterogeneous cloud machines. The findings of this research can be helpful for researchers and cloud solution architects to develop solutions that are robust against the stochastic nature of the inference time of DNN applications in the cloud and can offer a higher QoS to their users and avoid unintended outcomes.


title: SPH modeling of natural convection around a heated horizontal cylinder: A comparison with experiments

abstract: An experimental and numerical smoothed particle hydrodynamics (SPH) analysis was performed for the convective flow arising from a horizontal, thin cylindrical heat source enclosed in a glycerin-filled, slender enclosure at low Rayleigh numbers ($1.18\leq {\rm Ra}\leq 242$). Both the experiments and the SPH calculations were performed for positive ($0.1\leq\Delta T\leq 10$ K) and negative ($-10\leq\Delta T\leq -0.1$ K) temperature differences between the source and the surrounding fluid. In all cases a pair of steady, counter-rotating vortices is formed, accompanied by a plume of vertically ascending flow just above the source for $\Delta T>0$ and a vertically descending flow just below the source for $\Delta T<0$. The maximum flow velocities always occur within the ascending/descending plumes. The SPH predictions are found to match the experimental observations acceptably well with root-mean-square errors in the velocity profiles of the order of $\sim 10^{-5}$ m s$^{-1}$. The fact that the SPH method is able to reveal the detailed features of the flow phenomenon demonstrates the correctness of the approach.


title: Predicting orientation-dependent plastic susceptibility from static structure in amorphous solids via deep learning

abstract: It has been a long-standing materials science challenge to establish structure-property relations in amorphous solids. Here we introduce a rotation-variant local structure representation that enables different predictions for different loading orientations, which is found essential for high-fidelity prediction of the propensity for stress-driven shear transformations. This novel structure representation, when combined with convolutional neural network (CNN), a powerful deep learning algorithm, leads to unprecedented accuracy for identifying atoms with high propensity for shear transformations (i.e., plastic susceptibility), solely from the static structure - the spatial atomic positions - in both two- and three-dimensional model glasses. The data-driven models trained on samples at one composition and a given processing history are found transferrable to glass samples with different processing histories or at different compositions in the same alloy system. Our analysis of the new structure representation also provides valuable insight into key atomic packing features that influence the local mechanical response and its anisotropy in glasses.


title: Interdisciplinary Approaches to Understanding Artificial Intelligence's Impact on Society

abstract: Innovations in AI have focused primarily on the questions of "what" and "how"-algorithms for finding patterns in web searches, for instance-without adequate attention to the possible harms (such as privacy, bias, or manipulation) and without adequate consideration of the societal context in which these systems operate. In part, this is driven by incentives and forces in the tech industry, where a more product-driven focus tends to drown out broader reflective concerns about potential harms and misframings. But this focus on what and how is largely a reflection of the engineering and mathematics-focused training in computer science, which emphasizes the building of tools and development of computational concepts. As a result of this tight technical focus, and the rapid, worldwide explosion in its use, AI has come with a storm of unanticipated socio-technical problems, ranging from algorithms that act in racially or gender-biased ways, get caught in feedback loops that perpetuate inequalities, or enable unprecedented behavioral monitoring surveillance that challenges the fundamental values of free, democratic societies. Given that AI is no longer solely the domain of technologists but rather of society as a whole, we need tighter coupling of computer science and those disciplines that study society and societal values.


title: Next Wave Artificial Intelligence: Robust, Explainable, Adaptable, Ethical, and Accountable

abstract: The history of AI has included several "waves" of ideas. The first wave, from the mid-1950s to the 1980s, focused on logic and symbolic hand-encoded representations of knowledge, the foundations of so-called "expert systems". The second wave, starting in the 1990s, focused on statistics and machine learning, in which, instead of hand-programming rules for behavior, programmers constructed "statistical learning algorithms" that could be trained on large datasets. In the most recent wave research in AI has largely focused on deep (i.e., many-layered) neural networks, which are loosely inspired by the brain and trained by "deep learning" methods. However, while deep neural networks have led to many successes and new capabilities in computer vision, speech recognition, language processing, game-playing, and robotics, their potential for broad application remains limited by several factors. A concerning limitation is that even the most successful of today's AI systems suffer from brittleness-they can fail in unexpected ways when faced with situations that differ sufficiently from ones they have been trained on. This lack of robustness also appears in the vulnerability of AI systems to adversarial attacks, in which an adversary can subtly manipulate data in a way to guarantee a specific wrong answer or action from an AI system. AI systems also can absorb biases-based on gender, race, or other factors-from their training data and further magnify these biases in their subsequent decision-making. Taken together, these various limitations have prevented AI systems such as automatic medical diagnosis or autonomous vehicles from being sufficiently trustworthy for wide deployment. The massive proliferation of AI across society will require radically new ideas to yield technology that will not sacrifice our productivity, our quality of life, or our values.


title: Ultrasensitive Field-Effect Biosensors Enabled by the Unique Electronic Properties of Graphene

abstract: This review provides a critical overview of current developments on nanoelectronic biochemical sensors based on graphene. Composed of a single layer of conjugated carbon atoms, graphene has outstanding high carrier mobility and low intrinsic electrical noise, but a chemically inert surface. Surface functionalization is therefore crucial to unravel graphene sensitivity and selectivity for the detection of targeted analytes. To achieve optimal performance of graphene transistors for biochemical sensing, the tuning of the graphene surface properties via surface functionalization and passivation is highlighted, as well as the tuning of its electrical operation by utilizing multifrequency ambipolar configuration and a high frequency measurement scheme to overcome the Debye screening to achieve low noise and highly sensitive detection. Potential applications and prospectives of ultrasensitive graphene electronic biochemical sensors ranging from environmental monitoring and food safety, healthcare and medical diagnosis, to life science research, are presented as well.


title: Spatially Conditioned Graphs for Detecting Human-Object Interactions

abstract: We address the problem of detecting human-object interactions in images using graphical neural networks. Unlike conventional methods, where nodes send scaled but otherwise identical messages to each of their neighbours, we propose to condition messages between pairs of nodes on their spatial relationships, resulting in different messages going to neighbours of the same node. To this end, we explore various ways of applying spatial conditioning under a multi-branch structure. Through extensive experimentation we demonstrate the advantages of spatial conditioning for the computation of the adjacency structure, messages and the refined graph features. In particular, we empirically show that as the quality of the bounding boxes increases, their coarse appearance features contribute relatively less to the disambiguation of interactions compared to the spatial information. Our method achieves an mAP of 31.33% on HICO-DET and 54.2% on V-COCO, significantly outperforming state-of-the-art on fine-tuned detections.


title: Facile and Ultraclean Graphene-on-Glass Nanopores by Controlled Electrochemical Etching

abstract: A wide range of approaches have been explored to meet the challenges of graphene nanostructure fabrication, all requiring complex and high-end nanofabrication platform and suffering from surface contaminations, potentially giving electrical noise and increasing the thickness of the atomically thin graphene membrane. Here, with the use of an electrical pulse on a low capacitance graphene-on-glass (GOG) membrane, we fabricated clean graphene nanopores on commercially available glass substrates with exceptionally low electrical noise. In situ liquid AFM studies and electrochemical measurements revealed that both graphene nanopore nucleation and growth stem from the electrochemical attack on carbon atoms at defect sites, ensuring the creation of a graphene nanopore. Strikingly, compared to conventional TEM drilled graphene nanopores on SiN supporting membranes, GOG nanopores featured an order-of-magnitude reduced broadband noise, which we ascribed to the electrochemical refreshing of graphene nanopore on mechanically stable glass chips with negligible parasitic capacitance (1 pF). Further experiments on double-stranded DNA translocations demonstrated a greatly reduced current noise, and also confirmed the activation of single nanopores. Therefore, the exceptionally low noise and ease of fabrication will facilitate the understanding of the fundamental property and the application of such atomically thin nanopore sensors.


title: Faster Deterministic Modular Subset Sum

abstract: We consider the Modular Subset Sum problem: given a multiset $X$ of integers from $\mathbb{Z}_m$ and a target integer $t$, decide if there exists a subset of $X$ with a sum equal to $t \pmod{m}$. Recent independent works by Cardinal and Iacono (SOSA'21), and Axiotis et al. (SOSA'21) provided simple and near-linear algorithms for this problem. Cardinal and Iacono gave a randomized algorithm that runs in $O(m \log m)$ time, while Axiotis et al. gave a deterministic algorithm that runs in $O(m \text{ polylog } m)$ time. Both results work by reduction to a text problem, which is solved using a dynamic strings data structure. In this work, we develop a simple data structure, designed specifically to handle the text problem that arises in the algorithms for Modular Subset Sum. Our data structure, which we call the shift-tree, is a simple variant of a segment tree. We provide both a hashing-based and a deterministic variant of the shift-trees. We then apply our data structure to the Modular Subset Sum problem and obtain two algorithms. The first algorithm is Monte-Carlo randomized and matches the $O(m \log m)$ runtime of the Las-Vegas algorithm by Cardinal and Iacono. The second algorithm is fully deterministic and runs in $O(m \log m \cdot \alpha(m))$ time, where $\alpha$ is the inverse Ackermann function.


title: Deep Learning Approach for Matrix Completion Using Manifold Learning

abstract: Matrix completion has received vast amount of attention and research due to its wide applications in various study fields. Existing methods of matrix completion consider only nonlinear (or linear) relations among entries in a data matrix and ignore linear (or nonlinear) relationships latent. This paper introduces a new latent variables model for data matrix which is a combination of linear and nonlinear models and designs a novel deep-neural-network-based matrix completion algorithm to address both linear and nonlinear relations among entries of data matrix. The proposed method consists of two branches. The first branch learns the latent representations of columns and reconstructs the columns of the partially observed matrix through a series of hidden neural network layers. The second branch does the same for the rows. In addition, based on multi-task learning principles, we enforce these two branches work together and introduce a new regularization technique to reduce over-fitting. More specifically, the missing entries of data are recovered as a main task and manifold learning is performed as an auxiliary task. The auxiliary task constrains the weights of the network so it can be considered as a regularizer, improving the main task and reducing over-fitting. Experimental results obtained on the synthetic data and several real-world data verify the effectiveness of the proposed method compared with state-of-the-art matrix completion methods.


title: Three-dimensional black holes via Noether symmetries

abstract: We investigate the Noether symmetries of the Lagrangian for the stationary rotating BTZ-type three-dimensional spacetimes in $f(R)$ theory of gravity. A detailed analysis of Noether symmetries of (2+1)-dimensional rotating BTZ-type black hole spacetime model is presented. Applying the Noether symmetry approach, the first integrals (constants of motion) for each of Noether symmetries are obtained to look for the exact solutions. After solving the first integral equations depending on the form of the function $f(R)$, we derived some new (2+1)-dimensional rotating BTZ-type black hole solutions. We discussed the physical implications of the derived exact solutions. The thermodynamical properties of the obtained BTZ-type black hole solutions are analyzed by making use of the mass $M$ and the angular momentum $J$ in terms of $r_{\pm}$, where $r_+$ is the event horizon and $r_{-}$ is the inner horizon. Further, it is shown that thermodynamic quantities obey the first law, and the Smarr-like formulas of the solutions we found are obtained.


title: Coded sparse matrix computation schemes that leverage partial stragglers

abstract: Distributed matrix computations over large clusters can suffer from the problem of slow or failed worker nodes (called stragglers) which can dominate the overall job execution time. Coded computation utilizes concepts from erasure coding to mitigate the effect of stragglers by running 'coded' copies of tasks comprising a job; stragglers are typically treated as erasures. While this is useful, there are issues with applying, e.g., MDS codes in a straightforward manner. Several practical matrix computation scenarios involve sparse matrices. MDS codes typically require dense linear combinations of submatrices of the original matrices which destroy their inherent sparsity. This is problematic as it results in significantly higher worker computation times. Moreover, treating slow nodes as erasures ignores the potentially useful partial computations performed by them. Furthermore, some MDS techniques also suffer from significant numerical stability issues. In this work we present schemes that allow us to leverage partial computation by stragglers while imposing constraints on the level of coding that is required in generating the encoded submatrices. This significantly reduces the worker computation time as compared to previous approaches and results in improved numerical stability in the decoding process. Exhaustive numerical experiments on Amazon Web Services (AWS) clusters support our findings.


title: Maximizing the number of maximal independent sets of a fixed size

abstract: For a fixed graph G, a maximal independent set is an independent set that is not a proper subset of any other independent set. P. Erd\"os, and independently, J. W. Moon and L. Moser, and R. E. Miller and D. E. Muller, determined the maximum number of maximal independent sets in a graph on n vertices, as well as the extremal graphs. In this paper we maximize the number of maximal independent sets of a fixed size for all graphs of order n and determine the extremal graphs. Our result generalizes the classical result.


title: Gluing and lifting exact model structures for the recollement of exact categories

abstract: In this paper, we first provide an explicit procedure to glue together hereditary exact model structures for the recollement of exact categories. To that end, we use the notion of cotorsion pairs and we investigate the gluing of complete hereditary cotorsion pairs along the recollement of exact categories. Moreover, we study liftings of recollements of hereditary exact model structures to recollements of their associated homotopy categories. This leads to a new method to produce recollements of triangulated categories. Applications are given to contraderived categories, projective stable derived categories and stable categories of Gorenstein injective modules over an upper triangular matrix ring.


title: Two-dimensional adaptive membranes with programmable water and ionic channels

abstract: Membranes are ubiquitous in nature with primary functions that include adaptive filtering and selective transport of chemical and molecular species. Being critical to cellular functions, they are also fundamental in many areas of science and technology. Of particular importance are the adaptive and programmable membranes that can change their permeability or selectivity depending on the environment. Here, we explore implementation of such biological functions in artificial membranes and demonstrate two dimensional self assembled heterostructures of graphene oxide and polyamine macromolecules, forming a network of ionic channels that exhibit regulated permeability of water and monovalent ions. This permeability can be tuned by a change of pH or the presence of certain ions. Unlike traditional membranes, the regulation mechanism reported here relies on specific interactions between the membranes internal components and ions. This allows fabrication of membranes with programmable, predetermined permeability and selectivity, governed by the choice of components, their conformation and their charging state.


title: Power System Dynamic State Estimation Using Extended and Unscented Kalman Filters

abstract: Accurate estimation of power system dynamics is very important for the enhancement of power system reliability, resilience, security, and stability of power system. With the increasing integration of inverter-based distributed energy resources, the knowledge of power system dynamics has become more necessary and critical than ever before for proper control and operation of the power system. Although recent advancement of measurement devices and the transmission technologies have reduced the measurement and transmission error significantly, these measurements are still not completely free from the measurement noises. Therefore, the noisy measurements need to be filtered to obtain the accurate power system operating dynamics. In this work, the power system dynamic states are estimated using extended Kalman filter (EKF) and unscented Kalman filter (UKF). We have performed case studies on Western Electricity Coordinating Council (WECC)'s $3$-machine $9$-bus system and New England $10$-machine $39$-bus. The results show that the UKF and EKF can accurately estimate the power system dynamics. The comparative performance of EKF and UKF for the tested case is also provided. Other Kalman filtering techniques alongwith the machine learning-based estimator will be updated inthis report soon.All the sources code including Newton Raphsonpower flow, admittance matrix calculation, EKF calculation, andUKF calculation are publicly available in Github


title: Adaptive Submodular Meta-Learning

abstract: Meta-Learning has gained increasing attention in the machine learning and artificial intelligence communities. In this paper, we introduce and study an adaptive submodular meta-learning problem. The input of our problem is a set of items, where each item has a random state which is initially unknown. The only way to observe an item's state is to select that item. Our objective is to adaptively select a group of items that achieve the best performance over a set of tasks, where each task is represented as an adaptive submodular function that maps sets of items and their states to a real number. To reduce the computational cost while maintaining a personalized solution for each future task, we first select an initial solution set based on previously observed tasks, then adaptively add the remaining items to the initial solution set when a new task arrives. As compared to the solution where a brand new solution is computed for each new task, our meta-learning based approach leads to lower computational overhead at test time since the initial solution set is pre-computed in the training stage. To solve this problem, we propose a two-phase greedy policy and show that it achieves a $1/2$ approximation ratio for the monotone case. For the non-monotone case, we develop a two-phase randomized greedy policy that achieves a $1/32$ approximation ratio.


title: Anomalous Andreev Reflection on a Torus-Shaped Fermi Surface

abstract: Andreev reflection (AR) refers to the electron-hole conversion at the normal metal-superconductor interface. In a three-dimensional metal with spherical Fermi surface, retro (specular) AR can occur with the sign reversal of all three (a single) components of particle velocity. Here, we predict a novel type of AR with the inversion of two velocity components, dubbed "anomalous-trajectory Andreev reflection" (AAR), which can be realized in a class of materials with torus-shaped Fermi surface, such as doped nodal line semimetals. For its toroidal circle perpendicular to the interface, the Fermi torus doubles the AR channels and generates multiple AR processes. In particular, the AAR and retro AR are found to dominate electron transport in the light and heavy doping regimes, respectively. We show that the AAR visibly manifests as a ridge structure in the spatially resolved nonlocal conductance, in contrast to the peak structure for the retro AR. Our work opens a new avenue for the AR spectroscopy and offers a clear transport signature of torus-shaped Fermi surface.


title: The Genealogy of '$\lor$'

abstract: The use of the symbol $\lor$ for disjunction in formal logic is ubiquitous. Where did it come from? The paper details the evolution of the symbol $\lor$ in its historical and logical context. Some sources say that disjunction in its use as connecting propositions or formulas was introduced by Peano; others suggest that it originated as an abbreviation of the Latin word for "or", vel. We show that the origin of the symbol $\lor$ for disjunction can be traced to Whitehead and Russell's pre-Principia work in formal logic. Because of Principia's influence, its notation was widely adopted by philosophers working in logic (the logical empiricists in the 1920s and 1930s, especially Carnap and early Quine). Hilbert's adoption of $\lor$ in his Grundz\"uge der theoretischen Logic guaranteed its widespread use by mathematical logicians. The origins of other logical symbols are also discussed.


title: Windowed space-time least-squares Petrov-Galerkin method for nonlinear model order reduction

abstract: This work presents the windowed space-time least-squares Petrov-Galerkin method (WST-LSPG) for model reduction of nonlinear parameterized dynamical systems. WST-LSPG is a generalization of the space-time least-squares Petrov-Galerkin method (ST-LSPG). The main drawback of ST-LSPG is that it requires solving a dense space-time system with a space-time basis that is calculated over the entire global time domain, which can be unfeasible for large-scale applications. Instead of using a temporally-global space-time trial subspace and minimizing the discrete-in-time full-order model (FOM) residual over an entire time domain, the proposed WST-LSPG approach addresses this weakness by (1) dividing the time simulation into time windows, (2) devising a unique low-dimensional space-time trial subspace for each window, and (3) minimizing the discrete-in-time space-time residual of the dynamical system over each window. This formulation yields a problem with coupling confined within each window, but sequential across the windows. To enable high-fidelity trial subspaces characterized by a relatively minimal number of basis vectors, this work proposes constructing space-time bases using tensor decompositions for each window. WST-LSPG is equipped with hyper-reduction techniques to further reduce the computational cost. Numerical experiments for the one-dimensional Burgers' equation and the two-dimensional compressible Navier-Stokes equations for flow over a NACA 0012 airfoil demonstrate that WST-LSPG is superior to ST-LSPG in terms of accuracy and computational gain.


title: 6-qubit Optimal Clifford Circuits

abstract: Clifford group lies at the core of quantum computation -- it underlies quantum error correction, its elements can be used to perform magic state distillation and they form randomized benchmarking protocols, Clifford group is used to study quantum entanglement, and more. The ability to utilize Clifford group elements in practice relies heavily on the efficiency of their circuit-level implementation. Finding short circuits is a hard problem; despite Clifford group being finite, its size grows quickly with the number of qubits $n$, limiting known optimal implementations to $n{=}4$ qubits. For $n{=}6$, the number of Clifford group elements is about $2.1{\cdot}10^{23}$. In this paper, we report a set of algorithms, along with their C/C++ implementation, that implicitly synthesize optimal circuits for all 6-qubit Clifford group elements by storing a subset of the latter in a database of size 2.1TB (1KB=1024B). We demonstrate how to extract arbitrary optimal 6-qubit Clifford circuit in $0.0009358$ and $0.0006274$ seconds using consumer- and enterprise-grade computers (hardware) respectively, while relying on this database.


title: An algorithm for onset detection of linguistic segments in continuous electroencephalogram signals

abstract: A Brain Computer Interface based on imagined words can decode the word a subject is thinking on through brain signals to control an external device. In order to build a fully asynchronous Brain Computer Interface based on imagined words in electroencephalogram signals as source, we need to solve the problem of detecting the onset of the imagined words. Although there has been some research in this field, the problem has not been fully solved. In this paper we present an approach to solve this problem by using values from statistics, information theory and chaos theory as features to correctly identify the onset of imagined words in a continuous signal. On detecting the onsets of imagined words, the highest True Positive Rate achieved by our approach was obtained using features based on the generalized Hurst exponent, this True Positive Rate was 0.69 and 0.77 with a timing error tolerance region of 3 and 4 seconds respectively.


title: Smooth Bandit Optimization: Generalization to H\"older Space

abstract: We consider bandit optimization of a smooth reward function, where the goal is cumulative regret minimization. This problem has been studied for $\alpha$-H\"older continuous (including Lipschitz) functions with $0<\alpha\leq 1$. Our main result is in generalization of the reward function to H\"older space with exponent $\alpha>1$ to bridge the gap between Lipschitz bandits and infinitely-differentiable models such as linear bandits. For H\"older continuous functions, approaches based on random sampling in bins of a discretized domain suffices as optimal. In contrast, we propose a class of two-layer algorithms that deploy misspecified linear/polynomial bandit algorithms in bins. We demonstrate that the proposed algorithm can exploit higher-order smoothness of the function by deriving a regret upper bound of $\tilde{O}(T^\frac{d+\alpha}{d+2\alpha})$ for when $\alpha>1$, which matches existing lower bound. We also study adaptation to unknown function smoothness over a continuous scale of H\"older spaces indexed by $\alpha$, with a bandit model selection approach applied with our proposed two-layer algorithms. We show that it achieves regret rate that matches the existing lower bound for adaptation within the $\alpha\leq 1$ subset.


title: Casting Multiple Shadows: High-Dimensional Interactive Data Visualisation with Tours and Embeddings

abstract: Non-linear dimensionality reduction (NLDR) methods such as t-distributed stochastic neighbour embedding (t-SNE) are ubiquitous in the natural sciences, however, the appropriate use of these methods is difficult because of their complex parameterisations; analysts must make trade-offs in order to identify structure in the visualisation of an NLDR technique. We present visual diagnostics for the pragmatic usage of NLDR methods by combining them with a technique called the tour. A tour is a sequence of interpolated linear projections of multivariate data onto a lower dimensional space. The sequence is displayed as a dynamic visualisation, allowing a user to see the shadows the high-dimensional data casts in a lower dimensional view. By linking the tour to an NLDR view, we can preserve global structure and through user interactions like linked brushing observe where the NLDR view may be misleading. We display several case studies from both simulations and single cell transcriptomics, that shows our approach is useful for cluster orientation tasks.


title: Spatial distributions of the fields in guided normal modes of two coupled parallel optical nanofibers

abstract: We study the cross-sectional profiles and spatial distributions of the fields in guided normal modes of two coupled parallel optical nanofibers. We show that the distributions of the components of the field in a guided normal mode of two identical nanofibers are either symmetric or antisymmetric with respect to the radial principal axis and the tangential principal axis in the cross-sectional plane of the fibers. The symmetry of the magnetic field components with respect to the principal axes is opposite to that of the electric field components. We show that, in the case of even $\mathcal{E}_z$-cosine modes, the electric intensity distribution is dominant in the area between the fibers, with a saddle point at the two-fiber center. Meanwhile, in the case of odd $\mathcal{E}_z$-sine modes, the electric intensity distribution at the two-fiber center attains a local minimum of exactly zero. We find that the differences between the results of the coupled mode theory and the exact mode theory are large when the fiber separation distance is small and either the fiber radius is small or the light wavelength is large. We show that, in the case where the two nanofibers are not identical, the intensity distribution is symmetric about the radial principal axis and asymmetric about the tangential principal axis.


title: Random-gate-voltage induced Al'tshuler-Aronov-Spivak effect in topological edge states

abstract: Helical edge states are the hallmark of the quantum spin Hall insulator. Recently, several experiments have observed transport signatures contributed by trivial edge states, making it difficult to distinguish between the topologically trivial and nontrivial phases. Here, we show that helical edge states can be identified by the random-gate-voltage induced $\Phi_0/2$-period oscillation of the averaged electron return probability in the interferometer constructed by the edge states. The random gate voltage can highlight the $\Phi_0/2$-period Al'tshuler-Aronov-Spivak oscillation proportional to $\sin^2(2\pi\Phi/\Phi_0)$ by quenching the $\Phi_0$-period Aharonov-Bohm oscillation. It is found that the helical spin texture induced $\pi$ Berry phase is key to such weak antilocalization behavior with zero return probability at $\Phi=0$. In contrast, the oscillation for the trivial edge states may exhibit either weak localization or antilocalization depending on the strength of the spin-orbit coupling, which have finite return probability at $\Phi=0$. Our results provide an effective way for the identification of the helical edge states. The predicted signature is stabilized by the time-reversal symmetry so that it is robust against disorder and does not require any fine adjustment of system.


title: Hybrid microwave-optical scanning probe for addressing solid-state spins in nanophotonic cavities

abstract: Spin-photon interfaces based on solid-state atomic defects have enabled a variety of key applications in quantum information processing. To maximize the light-matter coupling strength, defects are often placed inside nanoscale devices. Efficiently coupling light and microwave radiation into these structures is an experimental challenge, especially in cryogenic or high vacuum environments with limited sample access. In this work, we demonstrate a fiber-based scanning probe that simultaneously couples light into a planar photonic circuit and delivers high power microwaves for driving electron spin transitions. The optical portion achieves 46% one-way coupling efficiency, while the microwave portion supplies an AC magnetic field with strength up to 9 Gauss. The entire probe can be scanned across a large number of devices inside a $^3$He cryostat without free-space optical access. We demonstrate this technique with silicon nanophotonic circuits coupled to single Er$^{3+}$ ions.


title: Deep Neural Networks Are Effective At Learning High-Dimensional Hilbert-Valued Functions From Limited Data

abstract: Accurate approximation of scalar-valued functions from sample points is a key task in computational science. Recently, machine learning with Deep Neural Networks (DNNs) has emerged as a promising tool for scientific computing, with impressive results achieved on problems where the dimension of the data or problem domain is large. This work broadens this perspective, focusing on approximating functions that are Hilbert-valued, i.e. take values in a separable, but typically infinite-dimensional, Hilbert space. This arises in science and engineering problems, in particular those involving solution of parametric Partial Differential Equations (PDEs). Such problems are challenging: 1) pointwise samples are expensive to acquire, 2) the function domain is high dimensional, and 3) the range lies in a Hilbert space. Our contributions are twofold. First, we present a novel result on DNN training for holomorphic functions with so-called hidden anisotropy. This result introduces a DNN training procedure and full theoretical analysis with explicit guarantees on error and sample complexity. The error bound is explicit in three key errors occurring in the approximation procedure: the best approximation, measurement, and physical discretization errors. Our result shows that there exists a procedure (albeit non-standard) for learning Hilbert-valued functions via DNNs that performs as well as, but no better than current best-in-class schemes. It gives a benchmark lower bound for how well DNNs can perform on such problems. Second, we examine whether better performance can be achieved in practice through different types of architectures and training. We provide preliminary numerical results illustrating practical performance of DNNs on parametric PDEs. We consider different parameters, modifying the DNN architecture to achieve better and competitive results, comparing these to current best-in-class schemes.


title: Effects of a thermal inversion experiment on STEM students learning and application of damped harmonic motion

abstract: There are diverse teaching methodologies to promote both collaborative and individual work in undergraduate physics courses. However, few educational studies seek to understand how students learn and apply new knowledge through open-ended activities that require mathematical modeling and experimentation focused on environmental problems. In this work, we propose a novel home experiment to simulate the dynamics of a particulate under temperature inversion and model it as damped harmonic motion. Twenty six first year students enrolled in STEM majors answered six qualitative questions after designing and developing the experiment. These questions helped analyze the students epistemological beliefs about their learning process of physics topics and its applications. Results showed that this type of open-ended experiments could facilitate the students understanding of physics phenomena. In addition, this experiment showed that it could help physics professors to promote students epistemological development by giving their students the opportunity to search for different sources of knowledge and becoming self-learners instead of looking at the professor as the epistemological authority. At the end, students described this activity as a positive experience that helped them realize alternative ways to apply physics topics in different contexts of their environment.


title: Rainbow Perfect and Near-Perfect Matchings in Complete Graphs with Edges Colored by Circular Distance

abstract: Given an edge-colored complete graph $K_n$ on $n$ vertices, a perfect (respectively, near-perfect) matching $M$ in $K_n$ with an even (respectively, odd) number of vertices is rainbow if all edges have distinct colors. In this paper, we consider an edge coloring of $K_n$ by circular distance, and we denote the resulting complete graph by $K^{\bullet}_n$. We show that when $K^{\bullet}_n$ has an even number of vertices, it contains a rainbow perfect matching if and only if $n=8k$ or $n=8k+2$, where $k$ is a nonnegative integer. In the case of an odd number of vertices, Kirkman matching is known to be a rainbow near-perfect matching in $K^{\bullet}_n$. However, real-world applications sometimes require multiple rainbow near-perfect matchings. We propose a method for using a recursive algorithm to generate multiple rainbow near-perfect matchings in $K^{\bullet}_n$.


title: Charged anisotropic compact objects obeying Karmarkar condition

abstract: This research develops a well-established analytical solution of the Einstein-Maxwell field equations. We analyze the behavior of a spherically symmetric and static interior driven by a charged anisotropic matter distribution. The class I methodology is used to close the system of equations and a suitable relation between the anisotropy factor and the electric field is imposed. The inner geometry of this toy model is described using an ansatz for the radial metric potential corresponding to the well-known isotropic Buchdahl space-time. The main properties are explored in order to determine if the obtained model is appropriate to represent a real compact body such as neutron or quark star. {We have fixed the mass and radii using the data of the compact objects} SMC X-1 and LMC X-4. It was found that the electric field and electric charge have magnitudes of the order of $\sim 10^{21}\ [V/cm]$ and $\sim 10^{20}\ [C]$, respectively. The magnitude of the electric field and electric charge depends on the dimensionless parameter $\chi$. To observe these effects on the total mass, mass-radius ratio and surface gravitational red-shift, we computed numerical data for different values of $\chi$.


title: Sizeable suppression of thermal Hall effect upon isotopic substitution in strontium titanate

abstract: We report measurements of the thermal Hall effect in single crystals of both pristine and isotopically substituted strontium titanate. We discovered a two orders of magnitude difference in the thermal Hall conductivity between $SrTi^{16}O_3$ and $^{18}O$-enriched $SrTi^{18}O_3$ samples. In most temperature ranges, the magnitude of thermal Hall conductivity ($\kappa_{xy}$) in $SrTi^{18}O_3$ is proportional to the magnitude of the longitudinal thermal conductivity ($\kappa_{xx}$), which suggests a phonon-mediated thermal Hall effect. However, they deviate in the temperature of their maxima, and the thermal Hall angle ratio ($|\kappa_{xy}/\kappa_{xx}|$) shows anomalously decreasing behavior below the ferroelectric Curie temperature $T_c$ ~$25 K$. This observation suggests a new underlying mechanism, as the conventional scenario cannot explain such differences within the slight change in phonon spectrum. Notably, the difference in magnitude of thermal Hall conductivity and rapidly decreasing thermal Hall angle ratio in $SrTi^{18}O_3$ is correlated with the strength of quantum critical fluctuations in this displacive ferroelectric. This relation points to a link between the quantum critical physics of strontium titanate and its thermal Hall effect, a possible clue to explain this example of an exotic phenomenon in non-magnetic insulating systems.


title: WITCHER : Detecting Crash Consistency Bugs in Non-volatile Memory Programs

abstract: The advent of non-volatile main memory (NVM) enables the development of crash-consistent software without paying storage stack overhead. However, building a correct crash-consistent program remains very challenging in the presence of a volatile cache. This paper presents WITCHER, a crash consistency bug detector for NVM software, that is (1) scalable -- does not suffer from test space explosion, (2) automatic -- does not require manual source code annotations, and (3) precise -- does not produce false positives. WITCHER first infers a set of "likely invariants" that are believed to be true to be crash consistent by analyzing source codes and NVM access traces. WITCHER automatically composes NVM images that simulate those potentially inconsistent (crashing) states violating the likely invariants. Then WITCHER performs "output equivalence checking" by comparing the output of program executions with and without a simulated crash. It validates if a likely invariant violation under test is a true crash consistency bug. Evaluation with ten persistent data structures, two real-world servers, and five example codes in Intel's PMDK library shows that WITCHER outperforms state-of-the-art tools. WITCHER discovers 37 (32 new) crash consistency bugs, which were all confirmed.


title: Monocular Real-time Full Body Capture with Inter-part Correlations

abstract: We present the first method for real-time full body capture that estimates shape and motion of body and hands together with a dynamic 3D face model from a single color image. Our approach uses a new neural network architecture that exploits correlations between body and hands at high computational efficiency. Unlike previous works, our approach is jointly trained on multiple datasets focusing on hand, body or face separately, without requiring data where all the parts are annotated at the same time, which is much more difficult to create at sufficient variety. The possibility of such multi-dataset training enables superior generalization ability. In contrast to earlier monocular full body methods, our approach captures more expressive 3D face geometry and color by estimating the shape, expression, albedo and illumination parameters of a statistical face model. Our method achieves competitive accuracy on public benchmarks, while being significantly faster and providing more complete face reconstructions.


title: Anomalous Hall resistivity and possible topological Hall effect in the EuAl$_4$ antiferromagnet

abstract: We report the observation of anomalous Hall resistivity in single crystals of EuAl$_4$, a centrosymmetric tetragonal compound, which exhibits coexisting antiferromagnetic (AFM) and charge-density-wave (CDW) orders with onset at $T_\mathrm{N} \sim 15.6$ K and $T_\mathrm{CDW} \sim 140$ K, respectively. In the AFM state, when the magnetic field is applied along the $c$-axis direction, EuAl$_4$ undergoes a series of metamagnetic transitions. Within this field range, we observe a clear hump-like anomaly in the Hall resistivity, representing part of the anomalous Hall resistivity. By considering different scenarios, we conclude that such a hump-like feature is most likely a manifestation of the topological Hall effect, normally occurring in noncentrosymmetric materials known to host nontrivial topological spin textures. In view of this, EuAl$_4$ would represent a rare case where the topological Hall effect not only arises in a centrosymmetric structure, but it also coexists with CDW order.


title: Collisional energy loss and the Chiral Magnetic Effect

abstract: Collisional energy loss of a fast particle in a medium is mostly due to the medium polarization by the electromagnetic fields of the particle. A small fraction of energy is carried away by the Cherenkov radiation. In chiral medium there is an additional contribution to the energy loss due to induction of the anomalous current proportional to the magnetic field. It causes the particle to lose energy mostly in the form of the \emph{chiral} Cherenkov radiation. We employ classical electrodynamics, adequate in a wide range of particle energies, to compute the collisional energy loss by a fast particle in a homogenous chiral plasma and apply the results to Quark-Gluon Plasma and a Weyl semimetal. In the later case photon spectrum is strongly enhanced in the ultraviolet and X-ray regions which makes it amenable to experimental investigation. Our main observation is that while the collisional energy loss in a non-chiral medium is a slow, at most logarithmic, function of energy $E$, the chiral Cherenkov radiation is proportional to $E^2$ when the recoil is neglected and to $E$ when it is taken into account.


title: Fast binomial-code holonomic quantum computation with ultrastrong light-matter coupling

abstract: We propose a protocol for bosonic binomial-code nonadiabatic holonomic quantum computation in a system composed of an artificial atom ultrastrongly coupled to a cavity resonator. In our protocol, the binomial codes, formed by superpositions of Fock states, can greatly save physical resources to correct errors in quantum computation. We apply to the system strong driving fields designed by shortcuts-to-adiabatic methods. This reduces the gate time to tens of nanoseconds. Decoherence of the system accumulated over time can be effectively reduced by such a fast evolution. Noise induced by control imperfections can be suppressed by a systematic-error-sensitivity nullification method, thus the protocol is mostly insensitive to such noises. As a result, this protocol can rapidly ($\sim 35$ ns) generate fault-tolerant and high-fidelity ($\gtrsim 98\%$ with experimentally realistic parameters) quantum gates.


title: Single-pixel Tracking and Imaging under Weak Illumination

abstract: Under weak illumination, tracking and imaging moving object turns out to be hard. By spatially collecting the signal, single pixel imaging schemes promise the capability of image reconstruction from low photon flux. However, due to the requirement on large number of samplings, how to clearly image moving objects is an essential problem for such schemes. Here we present a principle of single pixel tracking and imaging method. Velocity vector of the object is obtained from temporal correlation of the bucket signals in a typical computational ghost imaging system. Then the illumination beam is steered accordingly. Taking the velocity into account, both trajectory and clear image of the object are achieved during its evolution. Since tracking is achieved with bucket signals independently, this scheme is valid for capturing moving object as fast as its displacement within the interval of every sampling keeps larger than the resolution of the optical system. Experimentally, our method works well with the average number of detected photons down to 1.88 photons/speckle.


title: Ultra-bright multiplexed energy-time entangled photon generation from lithium niobate on insulator chip

abstract: High-flux entangled photon source is the key resource for quantum optical study and application. Here it is realized in a lithium niobate on isolator (LNOI) chip, with 2.79*10^11 Hz/mW photon pair rate and 1.53*10^9 Hz/nm/mW spectral brightness. These data are boosted by over two orders of magnitude compared to existing technologies. A 130-nm broad bandwidth is engineered for 8-channel multiplexed energy-time entanglement. Harnessed by high-extinction frequency correlation and Franson interferences up to 99.17% visibility, such energy-time entanglement multiplexing further enhances high-flux data rate, and warrants broad applications in quantum information processing on a chip.


title: A flexible sensitivity analysis approach for unmeasured confounding with multiple treatments and a binary outcome with application to SEER-Medicare lung cancer data

abstract: In the absence of a randomized experiment, a key assumption for drawing causal inference about treatment effects is the ignorable treatment assignment. Violations of the ignorability assumption may lead to biased treatment effect estimates. Sensitivity analysis helps gauge how causal conclusions will be altered in response to the potential magnitude of departure from the ignorability assumption. However, sensitivity analysis approaches for unmeasured confounding in the context of multiple treatments and binary outcomes are scarce. We propose a flexible Monte Carlo sensitivity analysis approach for causal inference in such settings. We first derive the general form of the bias introduced by unmeasured confounding, with emphasis on theoretical properties uniquely relevant to multiple treatments. We then propose methods to encode the impact of unmeasured confounding on potential outcomes and adjust the estimates of causal effects in which the presumed unmeasured confounding is removed. Our proposed methods embed nested multiple imputation within the Bayesian framework, which allow for seamless integration of the uncertainty about the values of the sensitivity parameters and the sampling variability, as well as use of the Bayesian Additive Regression Trees for modeling flexibility. Expansive simulations validate our methods and gain insight into sensitivity analysis with multiple treatments. We use the SEER-Medicare data to demonstrate sensitivity analysis using three treatments for early stage non-small cell lung cancer. The methods developed in this work are readily available in the R package SAMTx.


title: Generative Learning With Euler Particle Transport

abstract: We propose an Euler particle transport (EPT) approach for generative learning. The proposed approach is motivated by the problem of finding an optimal transport map from a reference distribution to a target distribution characterized by the Monge-Ampere equation. Interpreting the infinitesimal linearization of the Monge-Ampere equation from the perspective of gradient flows in measure spaces leads to a stochastic McKean-Vlasov equation. We use the forward Euler method to solve this equation. The resulting forward Euler map pushes forward a reference distribution to the target. This map is the composition of a sequence of simple residual maps, which are computationally stable and easy to train. The key task in training is the estimation of the density ratios or differences that determine the residual maps. We estimate the density ratios (differences) based on the Bregman divergence with a gradient penalty using deep density-ratio (difference) fitting. We show that the proposed density-ratio (difference) estimators do not suffer from the "curse of dimensionality" if data is supported on a lower-dimensional manifold. Numerical experiments with multi-mode synthetic datasets and comparisons with the existing methods on real benchmark datasets support our theoretical results and demonstrate the effectiveness of the proposed method.


title: High efficiency uniform wakefield acceleration of a positron beam using stable asymmetric mode in a hollow channel plasma

abstract: Plasma wakefield acceleration in the blowout regime is particularly promising for high-energy acceleration of electron beams because of its potential to simultaneously provide large acceleration gradients and high energy transfer efficiency while maintaining excellent beam quality. However, no equivalent regime for positron acceleration in plasma wakes has been discovered to-date. We show that after a short propagation distance, an asymmetric electron beam drives a stable wakefield in a hollow plasma channel that can be both accelerating and focusing for a positron beam. A high charge positron bunch placed at a suitable distance behind the drive bunch can beam-load or flatten the longitudinal wakefield and enhance the transverse focusing force, leading to high-efficiency and narrow energy spread acceleration of the positrons. Three-dimensional quasi-static particle-in-cell (PIC) simulations show that over 30% energy extraction efficiency from the wake to the positrons and 1% level energy spread can be simultaneously obtained, and further optimization is feasible.


title: Trace Impurity Transport in Multi-Species Plasmas with Large Particle Fluxes

abstract: A burning fusion plasma has a large inwards flux of fuel and outwards flux of ash. Existing impurity transport theories do not account for a steady-state environment of such large fluxes. In this paper, we extend classical transport theory to account for a background of such fluxes. In a mainly two-ion species magnetized plasma, with oppositely directed density gradients maintained by oppositely directed sources of the two species, trace impurity ions are shown to align with the larger mass ions. This unexpected behavior is derived analytically and simulated numerically using the MITNS (Multiple-Ion Transport Numerical Solver) code. The result suggests a potential advantage of edge-fueled p-B11 fusion in extracting high-Z impurities, and possible applications in plasma-based separation methods.


title: Embedding classical dynamics in a quantum computer

abstract: We develop a framework for simulating measure-preserving, ergodic dynamical systems on a quantum computer. Our approach provides a new operator-theoretic representation of classical dynamics by combining ergodic theory with quantum information science. The resulting quantum embedding of classical dynamics (QECD) enables efficient simulation of spaces of classical observables with exponentially large dimension using a quadratic number of quantum gates. The QECD framework is based on a quantum feature map for representing classical states by density operators on a reproducing kernel Hilbert space, $\mathcal H $, and an embedding of classical observables into self-adjoint operators on $\mathcal H$. In this scheme, quantum states and observables evolve unitarily under the lifted action of Koopman evolution operators of the classical system. Moreover, by virtue of the reproducing property of $\mathcal H$, the quantum system is pointwise-consistent with the underlying classical dynamics. To achieve an exponential quantum computational advantage, we project the state of the quantum system to a density matrix on a $2^n$-dimensional tensor product Hilbert space associated with $n$ qubits. By employing discrete Fourier-Walsh transforms, the evolution operator of the finite-dimensional quantum system is factorized into tensor product form, enabling implementation through a quantum circuit of size $O(n)$. Furthermore, the circuit features a state preparation stage, also of size $O(n)$, and a quantum Fourier transform stage of size $O(n^2)$, which makes predictions of observables possible by measurement in the standard computational basis. We prove theoretical convergence results for these predictions as $n\to\infty$. We present simulated quantum circuit experiments in Qiskit Aer, as well as actual experiments on the IBM Quantum System One.