title: Fourier Series-Based Approximation of Time-Varying Parameters in Ordinary Differential Equations

abstract: Many real-world systems modeled using differential equations involve unknown or uncertain parameters. Standard approaches to address parameter estimation inverse problems in this setting typically focus on estimating constants; yet some unobservable system parameters may vary with time without known evolution models. In this work, we propose a novel approximation method inspired by the Fourier series to estimate time-varying parameters in deterministic dynamical systems modeled with ordinary differential equations. Using ensemble Kalman filtering in conjunction with Fourier series-based approximation models, we detail two possible implementation schemes for sequentially updating the time-varying parameter estimates given noisy observations of the system states. We demonstrate the capabilities of the proposed approach in estimating periodic parameters, both when the period is known and unknown, as well as non-periodic time-varying parameters of different forms with several computed examples using a forced harmonic oscillator. Results emphasize the importance of the frequencies and number of approximation model terms on the time-varying parameter estimates and corresponding dynamical system predictions.


title: Compositional Construction of Abstractions for Infinite Networks of Discrete-Time Switched Systems

abstract: In this paper, we develop a compositional scheme for the construction of continuous approximations for interconnections of infinitely many discrete-time switched systems. An approximation (also known as abstraction) is itself a continuous-space system, which can be used as a replacement of the original (also known as concrete) system in a controller design process. Having designed a controller for the abstract system, it is refined to a more detailed one for the concrete system. We use the notion of so-called simulation functions to quantify the mismatch between the original system and its approximation. In particular, each subsystem in the concrete network and its corresponding one in the abstract network are related through a notion of local simulation functions. We show that if the local simulation functions satisfy certain small-gain type conditions developed for a network containing infinitely many subsystems, then the aggregation of the individual simulation functions provides an overall simulation function quantifying the error between the overall abstraction network and the concrete one. In addition, we show that our methodology results in a scale-free compositional approach for any finite-but-arbitrarily large networks obtained from truncation of an infinite network. We provide a systematic approach to construct local abstractions and simulation functions for networks of linear switched systems. The required conditions are expressed in terms of linear matrix inequalities that can be efficiently computed. We illustrate the effectiveness of our approach through an application to AC islanded microgirds.


title: 5G New Radio for Automotive, Rail, and Air Transport

abstract: The recent and upcoming releases of the 3rd Generation Partnership Project's 5G New Radio specifications include features that are motivated by providing connectivity services to a broad set of verticals, including the automotive, rail, and air transport industries. Currently, several radio access network features are being further enhanced or newly introduced in NR to improve 5G's capability to provide fast, reliable, and non-limiting connectivity for transport applications. In this article, we review the most important characteristics and requirements of a wide range of services that are driven by the desire to help the transport sector to become more sustainable, economically viable, safe, and secure. These requirements will be supported by the evolving and entirely new features of 5G NR systems, including accurate positioning, reference signal design to enable multi-transmission and reception points, service-specific scheduling configuration, and service quality prediction.


title: AdS Euclidean wormholes

abstract: We explore the construction and stability of asymptotically anti-de Sitter Euclidean wormholes in a variety of models. In simple ad hoc low-energy models, it is not hard to construct two-boundary Euclidean wormholes that dominate over disconnected solutions and which are stable (lacking negative modes) in the usual sense of Euclidean quantum gravity. Indeed, the structure of such solutions turns out to strongly resemble that of the Hawking-Page phase transition for AdS-Schwarzschild black holes, in that for boundary sources above some threshold we find both a `large' and a `small' branch of wormhole solutions with the latter being stable and dominating over the disconnected solution for large enough sources. We are also able to construct two-boundary Euclidean wormholes in a variety of string compactifications that dominate over the disconnected solutions we find and that are stable with respect to field-theoretic perturbations. However, as in classic examples investigated by Maldacena and Maoz, the wormholes in these UV-complete settings always suffer from brane-nucleation instabilities (even when sources that one might hope would stabilize such instabilities are tuned to large values). This indicates the existence of additional disconnected solutions with lower action. We discuss the significance of such results for the factorization problem of AdS/CFT.


title: Effect of substituting non-polar chains with polar chains on the structural dynamics of small organic molecule and polymer semiconductors

abstract: The processability and optoelectronic properties of organic semiconductors can be tuned and manipulated via chemical design. The substitution of the alkyl side chains by oligoethers has recently been successful for applications such as bioelectronic sensors and photocatalytic water-splitting. The carbon-oxygen bond in oligoethers is likely to render the system softer and more prone to dynamical disorder that can be detrimental to charge transport for example. We use neutron spectroscopy, X-Ray diffraction (XRD), differential scanning calorimetry (DSC) and polarized optical microscopy to study the effect of the substitution of n-hexyl (Hex) by triethylene glycol (TEG) on the structural dynamics of two organic semiconductors: a phenylene-bithiophene-phenylene (PTTP) molecule and a fluorene-co-dibenzothiophene (FS) polymer. Counterintuitively, inelastic neutron scattering (INS) reveals a softening of the modes of PTTP and FS with Hex side chains, pointing towards an increased dynamical disorder in these systems. However, T-dependent X-Ray and neutron diffraction, INS and DSC evidence an extra reversible transition close to room temperature (RT) for PTTP with TEG side chains. The observed transition, not accompanied by a change in birefringence, can also be observed by quasi-elastic neutron scattering. A fastening of the TEG side chains dynamics is observed in the case of PTTP and not FS. We therefore assign this transition to the melt of the TEG side chains which are promoting dynamical order at RT, but if crystallising, may introduce an extra reversible structural transition above RT leading to thermal instabilities. A deeper understanding of side chain polarity and structural dynamics can help guide materials design and navigate the intricate balance between electronic charge transport and aqueous swelling, sought for a number of emerging organic electronic and bioelectronic applications.


title: Virtual Memory Partitioning for Enhancing Application Performance in Mobile Platforms

abstract: Recently, the amount of running software on smart mobile devices is gradually increasing due to the introduction of application stores. The application store is a type of digital distribution platform for application software, which is provided as a component of an operating system on a smartphone or tablet. Mobile devices have limited memory capacity and, unlike server and desktop systems, due to their mobility they do not have a memory slot that can expand the memory capacity. Low memory killer (LMK) and out-of-memory killer (OOMK) are widely used memory management solutions in mobile systems. They forcibly terminate applications when the available physical memory becomes insufficient. In addition, before the forced termination, the memory shortage incurs thrashing and fragmentation, thus slowing down application performance. Although the existing page reclamation mechanism is designed to secure available memory, it could seriously degrade user responsiveness due to the thrashing. Memory management is therefore still important especially in mobile devices with small memory capacity. This paper presents a new memory partitioning technique that resolves the deterioration of the existing application life cycle induced by LMK and OOMK. It provides a completely isolated virtual memory node at the operating system level. Evaluation results demonstrate that the proposed method improves application execution time under memory shortage, compared with methods in previous studies.


title: Efficient MPI-based Communication for GPU-Accelerated Dask Applications

abstract: Dask is a popular parallel and distributed computing framework, which rivals Apache Spark to enable task-based scalable processing of big data. The Dask Distributed library forms the basis of this computing engine and provides support for adding new communication devices. It currently has two communication devices: one for TCP and the other for high-speed networks using UCX-Py -- a Cython wrapper to UCX. This paper presents the design and implementation of a new communication backend for Dask -- called MPI4Dask -- that is targeted for modern HPC clusters built with GPUs. MPI4Dask exploits mpi4py over MVAPICH2-GDR, which is a GPU-aware implementation of the Message Passing Interface (MPI) standard. MPI4Dask provides point-to-point asynchronous I/O communication coroutines, which are non-blocking concurrent operations defined using the async/await keywords from the Python's asyncio framework. Our latency and throughput comparisons suggest that MPI4Dask outperforms UCX by 6x for 1 Byte message and 4x for large messages (2 MBytes and beyond) respectively. We also conduct comparative performance evaluation of MPI4Dask with UCX using two benchmark applications: 1) sum of cuPy array with its transpose, and 2) cuDF merge. MPI4Dask speeds up the overall execution time of the two applications by an average of 3.47x and 3.11x respectively on an in-house cluster built with NVIDIA Tesla V100 GPUs for 1-6 Dask workers. We also perform scalability analysis of MPI4Dask against UCX for these applications on TACC's Frontera (GPU) system with upto 32 Dask workers on 32 NVIDIA Quadro RTX 5000 GPUs and 256 CPU cores. MPI4Dask speeds up the execution time for cuPy and cuDF applications by an average of 1.71x and 2.91x respectively for 1-32 Dask workers on the Frontera (GPU) system.


title: Privacy-Preserving and Efficient Verification of the Outcome in Genome-Wide Association Studies

abstract: Providing provenance in scientific workflows is essential for reproducibility and auditability purposes. Workflow systems model and record provenance describing the steps performed to obtain the final results of a computation. In this work, we propose a framework that verifies the correctness of the statistical test results that are conducted by a researcher while protecting individuals' privacy in the researcher's dataset. The researcher publishes the workflow of the conducted study, its output, and associated metadata. They keep the research dataset private while providing, as part of the metadata, a partial noisy dataset (that achieves local differential privacy). To check the correctness of the workflow output, a verifier makes use of the workflow, its metadata, and results of another statistical study (using publicly available datasets) to distinguish between correct statistics and incorrect ones. We use case the proposed framework in the genome-wide association studies (GWAS), in which the goal is to identify highly associated point mutations (variants) with a given phenotype. For evaluation, we use real genomic data and show that the correctness of the workflow output can be verified with high accuracy even when the aggregate statistics of a small number of variants are provided. We also quantify the privacy leakage due to the provided workflow and its associated metadata in the GWAS use-case and show that the additional privacy risk due to the provided metadata does not increase the existing privacy risk due to sharing of the research results. Thus, our results show that the workflow output (i.e., research results) can be verified with high confidence in a privacy-preserving way. We believe that this work will be a valuable step towards providing provenance in a privacy-preserving way while providing guarantees to the users about the correctness of the results.


title: Controller Synthesis for Hyperproperties

abstract: We investigate the problem of controller synthesis for hyperproperties specified in the temporal logic HyperLTL. Hyperproperties are system properties that relate multiple execution traces. Hyperproperties can elegantly express information-flow policies like noninterference and observational determinism. The controller synthesis problem is to automatically design a controller for a plant that ensures satisfaction of a given specification in the presence of the environment or adversarial actions. We show that the controller synthesis problem is decidable for HyperLTL specifications and finite-state plants. We provide a rigorous complexity analysis for different fragments of HyperLTL and different system types: tree-shaped, acyclic, and general graphs.


title: ({\alpha}, {\beta})-Modules in Graphs

abstract: Modular Decomposition focuses on repeatedly identifying a module M (a collection of vertices that shares exactly the same neighbourhood outside of M) and collapsing it into a single vertex. This notion of exactitude of neighbourhood is very strict, especially when dealing with real world graphs. We study new ways to relax this exactitude condition. However, generalizing modular decomposition is far from obvious. Most of the previous proposals lose algebraic properties of modules and thus most of the nice algorithmic consequences. We introduce the notion of an ({\alpha}, {\beta})-module, a relaxation that allows a bounded number of errors in each node and maintains some of the algebraic structure. It leads to a new combinatorial decomposition with interesting properties. Among the main results in this work, we show that minimal ({\alpha}, {\beta})-modules can be computed in polynomial time, and that every graph admits an ({\alpha},{\beta})-modular decomposition tree, thus generalizing Gallai's Theorem (which corresponds to the case for {\alpha} = {\beta} = 0). Unfortunately we give evidence that computing such a decomposition tree can be difficult.


title: Reheating in small-field inflation on the brane: The Swampland Criteria and observational constraints in light of the PLANCK 2018 results

abstract: We study cosmological inflation and its dynamics in the framework of the Randall-Sundrum II brane model. In particular, we analyze in detail four representative small-field inflationary potentials, namely Natural inflation, Hilltop inflation, Higgs-like inflation, and Exponential SUSY inflation, each characterized by two mass scales. We constrain the parameters for which a viable inflationary Universe emerges using the latest PLANCK results. Furthermore, we investigate whether or not those models in brane cosmology are consistent with the recently proposed Swampland Criteria, and give predictions for the duration of reheating as well as for the reheating temperature after inflation. Our results show that (i) the distance conjecture is satisfied, (ii) the de Sitter conjecture and its refined version may be avoided, and (iii) the allowed range for the five-dimensional Planck mass, $M_5$, is found to be between $10^5~\textrm{TeV}$ and $10^{12}~\textrm{TeV}$. Our main findings indicate that non-thermal leptogenesis cannot work within the framework of RS-II brane cosmology, at least for the inflationary potentials considered here.


title: Control of Plasmons in Doped Topological Insulators via Basis Atoms

abstract: Collective excitations in topologically non-trivial systems have attracted considerable attention in recent years. Here we study plasmons in the Su-Schrieffer-Heeger model whose low-energy electronic band is only partially filled, such that the system is metallic. Using the random phase approximation, we calculate the intra- and inter-band polarization functions and determine the bulk plasmonic dispersion from the dielectric function within the random phase approximation. We find that the sub-lattice basis states strongly affect the polarization functions and therefore control the system's plasmonic excitations. By varying the real-space separation of these local orbitals, one can thus selectively enhance or suppress the plasmonic energies via a tunable trade-off between intra-band and inter-band screening processes. Specifically, this mechanism can be used to stabilize undamped high energy plasmons that have already been reported in related models. We propose scenarios on how to control and observe these effects in experiments.


title: Direct Spatial Implementation of Sparse Matrix Multipliers for Reservoir Computing

abstract: Reservoir computing systems rely on the recurrent multiplication of a very large, sparse, fixed matrix. We argue that direct spatial implementation of these fixed matrices minimizes the work performed in the computation, and allows for significant reduction in latency and power through constant propagation and logic minimization. Bit-serial arithmetic enables massive static matrices to be implemented. We present the structure of our bit-serial matrix multiplier, and evaluate using canonical signed digit representation to further reduce logic utilization. We have implemented these matrices on a large FPGA and provide a cost model that is simple and extensible. These FPGA implementations, on average, reduce latency by 50x up to 86x versus GPU libraries. Comparing against a recent sparse DNN accelerator, we measure a 4.1x to 47x reduction in latency depending on matrix dimension and sparsity. Throughput of the FPGA solution is also competitive for a wide range of matrix dimensions and batch sizes. Finally, we discuss ways these techniques could be deployed in ASICs, making them applicable for dynamic sparse matrix computations.


title: User-Aware Power Management for Mobile Devices

abstract: The power management techniques to extend battery lifespan is becoming increasingly important due to longer user applications' running time in mobile devices. Even when users do not use any applications, battery lifespan decreases continually. It occurs because of service daemons of mobile platform and network-based data synchronization operations. In this paper, we propose a new power management system that recognizes the idle time of the device to reduce the battery consumption of mobile devices.


title: A co-Design approach to develop a smart cooking appliance. Applying a Domain Specific Language for a community supported appliance

abstract: Our environment, whether at work, in public spaces, or at home, is becoming more connected, and increasingly responsive. Meal preparation even when it involves simply heating ready-made food can be perceived as a complex process for people with disabilities. This research aimed to prototype, using a co-Design approach a Community Supported Appliance (CSA) by developing a Domain Specific Language (DSL), precisely created for a semi-automated cooking process. The DSL was shaped and expressed in the idiom of the users and allowed the CSA to support independence for users while performing daily cooking activities.


title: Distributed Compilation System for High-Speed Software Build Processes

abstract: The idle time of personal computers has increased steadily due to the generalization of computer usage and cloud computing. Clustering research aims at utilizing idle computer resources for processing a variable workload on a large number of computers. The workload is processed continually despite the volatile status of the individual computer resources. This paper proposes a distributed compilation system for improving the processing speed of CPU-intensive software compilations. This significantly reduces the compilation time of mass sources by using the idle resources. We expect gains of up to 65% compared to non-distributed compilation systems.


title: Analyzing Epistemic and Aleatoric Uncertainty for Drusen Segmentation in Optical Coherence Tomography Images

abstract: Age-related macular degeneration (AMD) is one of the leading causes of permanent vision loss in people aged over 60 years. Accurate segmentation of biomarkers such as drusen that points to the early stages of AMD is crucial in preventing further vision impairment. However, segmenting drusen is extremely challenging due to their varied sizes and appearances, low contrast and noise resemblance. Most existing literature, therefore, have focused on size estimation of drusen using classification, leaving the challenge of accurate segmentation less tackled. Additionally, obtaining the pixel-wise annotations is extremely costly and such labels can often be noisy, suffering from inter-observer and intra-observer variability. Quantification of uncertainty associated with segmentation tasks offers principled measures to inspect the segmentation output. Realizing its utility in identifying erroneous segmentation and the potential applications in clinical decision making, here we develop a U-Net based drusen segmentation model and quantify the segmentation uncertainty. We investigate epistemic and aleatoric uncertainty capturing model confidence and data uncertainty respectively. We present segmentation results and show how uncertainty can help formulate robust evaluation strategies. We visually inspect the pixel-wise uncertainty and segmentation results on test images. We finally analyze the correlation between segmentation uncertainty and accuracy. Our results demonstrate the utility of leveraging uncertainties in developing and explaining segmentation models for medical image analysis.


title: TAOS-CI: Lightweight & Modular Continuous Integration System for Edge Computing

abstract: With the proliferation of IoT and edge devices, we are observing a lot of consumer electronics becoming yet another IoT and edge devices. Unlike traditional smart devices, such as smart phones, consumer electronics, in general, have significant diversities with fewer number of devices per product model. With such high diversities, the proliferation of edge devices requires frequent and seamless updates of consumer electronics, which makes the manufacturers prone to regressions because the manufacturers have less resource per an instance of software release; i.e., they need to repeat releases by the number of product models times the number of updates. Continuous Integration (CI) systems can help prevent regression bugs from actively developing software packages including the frequently updated device software platforms. The proposed CI system provides a portable and modular software platform automatically inspecting potential issues of incoming changes with the enabled modules: code format and style, performance regressions, static checks on the source code, build and packaging tests, and dynamic checks with the built binary before deploying a platform image on the IoT and edge devices. Besides, our proposed approach is lightweight enough to be hosted in normal desktop computers even for dozens of developers. As a result, it can be easily applied to a lot of various source code repositories. Evaluation results demonstrate that the proposed method drastically improves plug-ins execution time and memory consumption, compared with methods in previous studies.


title: Distilling Large Language Models into Tiny and Effective Students using pQRNN

abstract: Large pre-trained multilingual models like mBERT, XLM-R achieve state of the art results on language understanding tasks. However, they are not well suited for latency critical applications on both servers and edge devices. It's important to reduce the memory and compute resources required by these models. To this end, we propose pQRNN, a projection-based embedding-free neural encoder that is tiny and effective for natural language processing tasks. Without pre-training, pQRNNs significantly outperform LSTM models with pre-trained embeddings despite being 140x smaller. With the same number of parameters, they outperform transformer baselines thereby showcasing their parameter efficiency. Additionally, we show that pQRNNs are effective student architectures for distilling large pre-trained language models. We perform careful ablations which study the effect of pQRNN parameters, data augmentation, and distillation settings. On MTOP, a challenging multilingual semantic parsing dataset, pQRNN students achieve 95.9\% of the performance of an mBERT teacher while being 350x smaller. On mATIS, a popular parsing task, pQRNN students on average are able to get to 97.1\% of the teacher while again being 350x smaller. Our strong results suggest that our approach is great for latency-sensitive applications while being able to leverage large mBERT-like models.


title: Cloud-Based Content Cooperation System to Assist Collaborative Learning Environment

abstract: Online educational systems running on smart devices have the advantage of allowing users to learn online regardless of the location of the users. In particular, data synchronization enables users to cooperate on contents in real time anywhere by sharing their files via cloud storage. However, users cannot collaborate by simultaneously modifying files that are shared with each other. In this paper, we propose a content collaboration method and a history management technique that are based on distributed system structure and can synchronize data shared in the cloud for multiple users and multiple devices.


title: Fano lineshapes and Rabi splittings: Can they be artificially generated or obscured by the numerical aperture?

abstract: Fano resonances and Rabi splittings are routinely reported in the scientific literature. Asymmetric resonance lineshapes are usually associated with Fano resonances, and two split peaks in the spectrum are often attributed to a Rabi splitting. True Fano resonances and Rabi splittings are unequivocal signatures of coherent coupling between subsystems. However, can the same spectral lineshapes characterizing Fano resonances and Rabi splittings arise from a purely incoherent sum of intensities? Here we answer this question through experiments with a tunable Fabry-P\'erot cavity containing a CsPbBr\textsubscript{3} perovskite crystal. By measuring the transmission and photoluminescence of this system using microscope objectives with different numerical aperture ($NA$), we find that even a modest $NA=0.4$ can artificially generate Fano resonances and Rabi splittings. We furthermore show that this modest $NA$ can obscure the anti-crossing of a bona fide strongly coupled light-matter system. Through transfer matrix calculations we confirm that these spectral artefacts are due to the incoherent sum of transmitted intensities at different angles captured by the $NA$. Our results are relevant to the wide nanophotonics community, characterizing dispersive optical systems with high numerical aperture microscope objectives. We conclude with general guidelines to avoid pitfalls in the characterization of such optical systems.


title: Common kernel-smoothed proper orthogonal decomposition (CKSPOD): An efficient reduced-order model for emulation of spatiotemporally evolving flow dynamics

abstract: In the present study, we propose a new surrogate model, called common kernel-smoothed proper orthogonal decomposition (CKSPOD), to efficiently emulate the spatiotemporal evolution of fluid flow dynamics. The proposed surrogate model integrates and extends recent developments in Gaussian process learning, high-fidelity simulations, projection-based model reduction, uncertainty quantification, and experimental design, rendering a systematic, multidisciplinary framework. The novelty of the CKSPOD emulation lies in the construction of a common Gram matrix, which results from the Hadamard product of Gram matrices of all observed design settings. The Gram matrix is a spatially averaged temporal correlation matrix and contains the temporal dynamics of the corresponding sampling point. The common Gram matrix synthesizes the temporal dynamics by transferring POD modes into spatial functions at each observed design setting, which remedies the phase-difference issue encountered in the kernel-smoothed POD (KSPOD) emulation, a recent fluid flow emulator proposed in Chang et al. (2020). The CKSPOD methodology is demonstrated through a model study of flow dynamics of swirl injectors with three design parameters. A total of 30 training design settings and 8 validation design settings are included. Both qualitative and quantitative results show that the CKSPOD emulation outperforms the KSPOD emulation for all validation cases, and is capable of capturing small-scale wave structures on the liquid-film surface faithfully. The turbulent kinetic energy prediction using CKSPOD reveals lower predictive uncertainty than KSPOD, thereby allowing for more accurate and precise flow predictions. The turnaround time of the CKSPOD emulation is about 5 orders of magnitude faster than the corresponding high-fidelity simulation, which enables an efficient and scalable framework for design exploration and optimization.


title: Generative Replay-based Continual Zero-Shot Learning

abstract: Zero-shot learning is a new paradigm to classify objects from classes that are not available at training time. Zero-shot learning (ZSL) methods have attracted considerable attention in recent years because of their ability to classify unseen/novel class examples. Most of the existing approaches on ZSL works when all the samples from seen classes are available to train the model, which does not suit real life. In this paper, we tackle this hindrance by developing a generative replay-based continual ZSL (GRCZSL). The proposed method endows traditional ZSL to learn from streaming data and acquire new knowledge without forgetting the previous tasks' gained experience. We handle catastrophic forgetting in GRCZSL by replaying the synthetic samples of seen classes, which have appeared in the earlier tasks. These synthetic samples are synthesized using the trained conditional variational autoencoder (VAE) over the immediate past task. Moreover, we only require the current and immediate previous VAE at any time for training and testing. The proposed GRZSL method is developed for a single-head setting of continual learning, simulating a real-world problem setting. In this setting, task identity is given during training but unavailable during testing. GRCZSL performance is evaluated on five benchmark datasets for the generalized setup of ZSL with fixed and dynamic (incremental class) settings of continual learning. The existing class setting presented recently in the literature is not suitable for a class-incremental setting. Therefore, this paper proposes a new setting to address this issue. Experimental results show that the proposed method significantly outperforms the baseline and the state-of-the-art method and makes it more suitable for real-world applications.


title: Iterative Optimisation with an Innovation CNN for Pose Refinement

abstract: Object pose estimation from a single RGB image is a challenging problem due to variable lighting conditions and viewpoint changes. The most accurate pose estimation networks implement pose refinement via reprojection of a known, textured 3D model, however, such methods cannot be applied without high quality 3D models of the observed objects. In this work we propose an approach, namely an Innovation CNN, to object pose estimation refinement that overcomes the requirement for reprojecting a textured 3D model. Our approach improves initial pose estimation progressively by applying the Innovation CNN iteratively in a stochastic gradient descent (SGD) framework. We evaluate our method on the popular LINEMOD and Occlusion LINEMOD datasets and obtain state-of-the-art performance on both datasets.


title: A Self-Updating K-Contingency List for Smart Grid System

abstract: A reliable decision making by the operator in a smart grid is contingent upon correct analysis of intra-and-interdependencies between its entities and also on accurate identification of the most critical entities at a given point of time. A measurement based self-updating contingency list can provide real-time information to the operator about current system condition which can help the operator to take the required action. In this paper, the underlying intra-and-interdependencies between entities for a given power-communication network is captured using a dependency model called Modified Implicative Interdependency Model (MIIM) [1]. Given an integer K, the event-driven self-updating contingency list problem gives the list of K-most critical entities, failure of which maximizes the network damage at the current time. Owing to the problem being NP complete, a fast heuristic method to generate a real-time contingency list using system measurements is provided here. The validation of the work is done by comparing the contingency list obtained for different K values using the MIIM model on a smart grid of IEEE 14-Bus system with that obtained by simulating the smart grid using a co-simulation system formed by MATPOWER and Java Network Simulator (JNS). The results also indicate that the network damage predicted by both the ILP based solution [2] and the proposed heuristic solution using MIIM are more realistic compared to that obtained using another dependency model called Implicative Interdependency Model (IIM) [3].


title: Meshless Fragile Points Methods Based on Petrov-Galerkin Weak-Forms for Transient Heat Conduction Problems in Complex Anisotropic Nonhomogeneous Media

abstract: Three kinds of Fragile Points Methods based on Petrov-Galerkin weak-forms (PG-FPMs) are proposed for analyzing heat conduction problems in nonhomogeneous anisotropic media. This is a follow-up of the previous study on the original FPM based on a symmetric Galerkin weak-form. The trial function is piecewise-continuous, written as local Taylor expansions at the Fragile Points. A modified Radial Basis Function-based Differential Quadrature (RBF-DQ) method is employed for establishing the local approximation. The Dirac delta function, Heaviside step function, and the local fundamental solution of the governing equation are alternatively used as test functions. Vanishing or pure contour integral formulation in subdomains or on local boundaries can be obtained. Extensive numerical examples in 2D and 3D are provided as validations. The collocation method (PG-FPM-1) is superior in transient analysis with arbitrary point distribution and domain partition. The finite volume method (PG-FPM-2) shows the best efficiency, saving 25% to 50% computational time comparing with the Galerkin FPM. The singular solution method (PG-FPM-3) is highly efficient in steady-state analysis. The anisotropy and nonhomogeneity give rise to no difficulties in all the methods. The proposed PG-FPM approaches represent an improvement to the original Galerkin FPM, as well as to other meshless methods in earlier literature.


title: Consecutive primes which are widely digitally delicate

abstract: We show that for every positive integer $k$, there exist $k$ consecutive primes having the property that if any digit of any one of the primes, including any of the infinitely many leading zero digits, is changed, then that prime becomes composite.


title: Parabolic velocity profile causes drift of inertial prolate spheroids -- but gravity is stronger

abstract: Motion of elongated particles in shear is studied. In applications where particles are much heavier than the carrying fluid, e.g. aerosols, the influence of particle inertia dominates the particle dynamics. Assuming that the particle only experiences a local linear velocity profile, its rotational and translational motion are independent. However, we show that quadratic terms of the local velocity profile combined with particle inertia cause a lateral drift of prolate spheroidal particles. We find that this drift is maximal when particle inertial forces are of the same order of magnitude as viscous forces, and that both extremely light and extremely heavy particles have negligible drift. In the non-inertial case, the particle rotates according to the local linear velocity profile, with each instantaneous orientation corresponding to a velocity that gives zero force on the particle. This results in a translational motion in the flow direction with periodic velocity fluctuations. With added particle inertia, the particle is slow to react to the surrounding fluid motion and the particle will switch between slower and faster rotation compared to the zero-force solution. The final motion that gives zero integrated force over a rotational period, is a motion with a lateral drift. We show that this drift is purely an effect of the non-sphericity of the particle and its translational inertia, while rotational inertia is negligible. Finally, although this inertial drift will contribute to the lateral motion of heavy elongated particles in channel flow, sedimentation due to gravity will dominate in any practical application on earth.


title: Energy estimates and convergence of weak solutions of the porous medium equation

abstract: We study the convergence of the weak solution of the porous medium equation with a type of Robin boundary conditions, by tuning a parameter either to zero or to infinity. The convergence is in the strong sense, with respect to the $L^2$-norm, and the limiting function solves the same equation with Neumann (resp. Dirichlet) boundary conditions when the parameter is taken to zero (resp. infinity). Our approach is to consider an underlying microscopic dynamics whose space-time evolution of the density is ruled by the solution of those equations and from this, we derive sufficiently strong energy estimates which are the keystone to the proof of our convergence result.


title: Event vertex and time reconstruction in large volume liquid scintillator detector

abstract: Large-volume liquid scintillator detectors with ultra-low background levels have been widely used to study neutrino physics and search for dark matter. Event vertex and event time are not only useful for event selection but also essential for the reconstruction of event energy. In this study, four event vertex and event time reconstruction algorithms using charge and time information collected by photomultiplier tubes were analyzed comprehensively. The effects of photomultiplier tube properties were also investigated. The results indicate that the transit time spread is the main effect degrading the vertex reconstruction, while the effect of dark noise is limited. In addition, when the event is close to the detector boundary, the charge information provides better performance for vertex reconstruction than the time information.


title: Length functions on groups and rigidity

abstract: Let $G$ be a group. A function $l:G\rightarrow \lbrack 0,\infty )$ is called a length function if (1) $l(g^{n})=|n|l(g)$ for any $g\in G$ and $n\in \mathbb{Z};$ (2) $l(hgh^{-1})=l(g)$ for any $h,g\in G;$ and (3) $l(ab)\leq l(a)+l(b)$ for commuting elements $a,b.$ Such length functions exist in many branches of mathematics, mainly as stable word lengths, stable norms, smooth measure-theoretic entropy, translation lengths on $\mathrm{CAT}(0)$ spaces and Gromov $\delta $% -hyperbolic spaces, stable norms of quasi-cocycles, rotation numbers of circle homeomorphisms, dynamical degrees of birational maps and so on. We study length functions on Lie groups, Gromov hyperbolic groups, arithmetic subgroups, matrix groups over rings and Cremona groups. As applications, we prove that every group homomorphism from an arithmetic subgroup of a simple algebraic $\mathbb{Q}$-group of $\mathbb{Q}$-rank at least $2,$ or a finite-index subgroup of the elementary group $E_{n}(R)$ $(n\geq 3)$ over an associative ring, or the Cremona group $\mathrm{Bir}(P_{\mathbb{C}}^{2})$ to any group $G$ having a purely positive length function must have its image finite. Here $G$ can be outer automorphism group $\mathrm{Out}(F_{n})$ of free groups, mapping classes group $\mathrm{MCG}(\Sigma_{g})$, $\mathrm{CAT}% (0)$ groups or Gromov hyperbolic groups, or the group $\mathrm{Diff}(\Sigma ,\omega )$ of diffeomorphisms of a hyperbolic closed surface preserving an area form $\omega .$


title: Newcomer OSS-Candidates: Characterizing Contributions of Novice Developers to GitHub

abstract: The ability of an Open Source Software (OSS) project to attract, onboard, and retain any newcomer is vital to its livelihood. Although, evidence suggests an upsurge in novice developers joining social coding platforms (such as GitHub), the extent to which their activities result in a OSS contribution is unknown. Henceforth, we execute the protocols of a registered report to study activities of a "Newcomer OSS-Candidate", who is a novice developer that is new to that social coding platform, and has the intention to later onboard an OSS project. Using GitHub as a case platform, we analyze 171 identified Newcomer OSS-Candidates to characterize their contribution activities. Results show that Newcomer OSS-Candidates are likely to target software based repositories (i.e., 66%), and their first contributions are mainly associated with development (commits) and maintenance (PRs). Newcomer OSS-Candidates are less likely to practice social coding, but eventually end up onboarding (i.e., 30% quantitative, 70% follow-up survey) an OSS project. Furthermore, they cite finding a way to start as the most challenging barrier to contribute. Our work reveals insights on how newcomers to social coding platforms are potential sources of OSS contributions.


title: Applications of artificial intelligence in drug development using real-world data

abstract: The US Food and Drug Administration (FDA) has been actively promoting the use of real-world data (RWD) in drug development. RWD can generate important real-world evidence reflecting the real-world clinical environment where the treatments are used. Meanwhile, artificial intelligence (AI), especially machine- and deep-learning (ML/DL) methods, have been increasingly used across many stages of the drug development process. Advancements in AI have also provided new strategies to analyze large, multidimensional RWD. Thus, we conducted a rapid review of articles from the past 20 years, to provide an overview of the drug development studies that use both AI and RWD. We found that the most popular applications were adverse event detection, trial recruitment, and drug repurposing. Here, we also discuss current research gaps and future opportunities.


title: A symbolic information approach to characterize response-related differences in cortical activity during a Go/No-Go task

abstract: How the brain processes information from external stimuli in order to perceive the world and act on it is one of the greatest questions in neuroscience. To address this question different time series analyzes techniques have been employed to characterize the statistical properties of brain signals during cognitive tasks. Typically response-specific processes are addressed by comparing the time course of average event-related potentials in different trials type. Here we analyze monkey Local Field Potentials data during visual pattern discrimination called Go/No-Go task in the light of information theory quantifiers. We show that the Bandt-Pompe symbolization methodology to calculate entropy and complexity of data is a useful tool to distinguish response-related differences between Go and No-Go trials. We propose to use an asymmetry index to statistically validate trial type differences. Moreover, by using the multi-scale approach and embedding time delays to downsample the data we can estimate the important time scales in which the relevant information is been processed.


title: Aharonov-Bohm-Like Scattering in the Generalized Uncertainty Principle-corrected Quantum Mechanics

abstract: We discuss classical electrodynamics and the Aharonov-Bohm effect in the presence of the minimal length. In the former we derive the classical equation of motion and the corresponding Lagrangian. In the latter we adopt the generalized uncertainty principle (GUP) and compute the scattering cross section up to the first-order of the GUP parameter $\beta$. Even though the minimal length exists, the cross section is invariant under the simultaneous change $\phi \rightarrow -\phi$, $\alpha' \rightarrow -\alpha'$, where $\phi$ and $\alpha'$ are azimuthal angle and magnetic flux parameter. However, unlike the usual Aharonv-Bohm scattering the cross section exhibits discontinuous behavior at every integer $\alpha'$. The symmetries, which the cross section has in the absence of GUP, are shown to be explicitly broken at the level of ${\cal O} (\beta)$.


title: Quasi-solvable lattice models for $\mathrm{Sp}_{2n}$ and $\mathrm{SO}_{2n+1}$ Demazure atoms and characters

abstract: We construct colored lattice models whose partition functions represent symplectic and odd orthogonal Demazure characters and atoms. We show that our lattice models are not solvable, but we are able to show the existence of sufficiently many solutions of the Yang--Baxter equation that allows us to compute functional equations for the corresponding partition functions. From these functional equations, we determine that the partition function of our models are the Demazure atoms and characters for the symplectic and odd orthogonal Lie groups. We coin our lattice models as quasi-solvable. We use the natural bijection of admissible states in our models with Proctor patterns to give a right key algorithm for reverse King tableaux and Sundaram tableaux.


title: Minimizing Age of Incorrect Information for Unreliable Channel with Power Constraint

abstract: Age of Incorrect Information (AoII) is a newly introduced performance metric that considers communication goals. Therefore, comparing with traditional performance metrics and the recently introduced metric - Age of Information (AoI), AoII achieves better performance in many real-life applications. However, the fundamental nature of AoII has been elusive so far. In this paper, we consider the AoII in a system where a transmitter sends updates about a multi-state Markovian source to a remote receiver through an unreliable channel. The communication goal is to minimize AoII subject to a power constraint. We cast the problem into a Constrained Markov Decision Process (CMDP) and prove that the optimal policy is a mixture of two deterministic threshold policies. Afterward, by leveraging the notion of Relative Value Iteration (RVI) and the structural properties of threshold policy, we propose an efficient algorithm to find the threshold policies as well as the mixing coefficient. Lastly, numerical results are laid out to highlight the performance of AoII-optimal policy.


title: Study of Pre-processing Defenses against Adversarial Attacks on State-of-the-art Speaker Recognition Systems

abstract: Adversarial examples to speaker recognition (SR) systems are generated by adding a carefully crafted noise to the speech signal to make the system fail while being imperceptible to humans. Such attacks pose severe security risks, making it vital to deep-dive and understand how much the state-of-the-art SR systems are vulnerable to these attacks. Moreover, it is of greater importance to propose defenses that can protect the systems against these attacks. Addressing these concerns, this paper at first investigates how state-of-the-art x-vector based SR systems are affected by white-box adversarial attacks, i.e., when the adversary has full knowledge of the system. x-Vector based SR systems are evaluated against white-box adversarial attacks common in the literature like fast gradient sign method (FGSM), basic iterative method (BIM)--a.k.a. iterative-FGSM--, projected gradient descent (PGD), and Carlini-Wagner (CW) attack. To mitigate against these attacks, the paper proposes four pre-processing defenses. It evaluates them against powerful adaptive white-box adversarial attacks, i.e., when the adversary has full knowledge of the system, including the defense. The four pre-processing defenses--viz. randomized smoothing, DefenseGAN, variational autoencoder (VAE), and Parallel WaveGAN vocoder (PWG) are compared against the baseline defense of adversarial training. Conclusions indicate that SR systems were extremely vulnerable under BIM, PGD, and CW attacks. Among the proposed pre-processing defenses, PWG combined with randomized smoothing offers the most protection against the attacks, with accuracy averaging 93% compared to 52% in the undefended system and an absolute improvement >90% for BIM attacks with $L_\infty>0.001$ and CW attack.


title: Single Neuron Segmentation using Graph-based Global Reasoning with Auxiliary Skeleton Loss from 3D Optical Microscope Images

abstract: One of the critical steps in improving accurate single neuron reconstruction from three-dimensional (3D) optical microscope images is the neuronal structure segmentation. However, they are always hard to segment due to the lack in quality. Despite a series of attempts to apply convolutional neural networks (CNNs) on this task, noise and disconnected gaps are still challenging to alleviate with the neglect of the non-local features of graph-like tubular neural structures. Hence, we present an end-to-end segmentation network by jointly considering the local appearance and the global geometry traits through graph reasoning and a skeleton-based auxiliary loss. The evaluation results on the Janelia dataset from the BigNeuron project demonstrate that our proposed method exceeds the counterpart algorithms in performance.


title: Construction of linearly independent non-orthogonal AGP states

abstract: We show how to construct a linearly independent set of antisymmetrized geminal power (AGP) states, which allows us to rewrite our recently introduced geminal replacement models as linear combinations of non-orthogonal AGPs. This greatly simplifies the evaluation of matrix elements and permits us to introduce an AGP-based selective configuration interaction method, which can reach arbitrary excitation levels relative to a reference AGP, balancing accuracy and cost as we see fit.


title: A Two-Stage Deep Learning Detection Classifier for the ATLAS Asteroid Survey

abstract: In this paper we present a two-step neural network model to separate detections of solar system objects from optical and electronic artifacts in data obtained with the "Asteroid Terrestrial-impact Last Alert System" (ATLAS), a near-Earth asteroid sky survey system [arXiv:1802.00879]. A convolutional neural network [arXiv:1807.10912] is used to classify small "postage-stamp" images of candidate detections of astronomical sources into eight classes, followed by a multi-layered perceptron that provides a probability that a temporal sequence of four candidate detections represents a real astronomical source. The goal of this work is to reduce the time delay between Near-Earth Object (NEO) detections and submission to the Minor Planet Center. Due to the rare and hazardous nature of NEOs [Harris and D'Abramo, 2015], a low false negative rate is a priority for the model. We show that the model reaches 99.6\% accuracy on real asteroids in ATLAS data with a 0.4\% false negative rate. Deployment of this model on ATLAS has reduced the amount of NEO candidates that astronomers must screen by 90%, thereby bringing ATLAS one step closer to full autonomy.


title: Implicit shock tracking for unsteady flows by the method of lines

abstract: A recently developed high-order implicit shock tracking (HOIST) framework for resolving discontinuous solutions of inviscid, steady conservation laws [41, 43] is extended to the unsteady case. Central to the framework is an optimization problem which simultaneously computes a discontinuity-aligned mesh and the corresponding high-order approximation to the flow, which provides nonlinear stabilization and a high-order approximation to the solution. This work extends the implicit shock tracking framework to the case of unsteady conservation laws using a method of lines discretization via a diagonally implicit Runge-Kutta method by "solving a steady problem at each timestep". We formulate and solve an optimization problem that produces a feature-aligned mesh and solution at each Runge-Kutta stage of each timestep, and advance this solution in time by standard Runge-Kutta update formulas. A Rankine-Hugoniot based prediction of the shock location together with a high-order, untangling mesh smoothing procedure provides a high-quality initial guess for the optimization problem at each time, which results in rapid convergence of the sequential quadratic programing (SQP) optimization solver. This method is shown to deliver highly accurate solutions on coarse, high-order discretizations without nonlinear stabilization and recover the design accuracy of the Runge-Kutta scheme. We demonstrate this framework on a series of inviscid, unsteady conservation laws in both one- and two- dimensions. We also verify that our method is able to recover the design order of accuracy of our time integrator in the presence of a strong discontinuity.


title: A Contextualist Decision Theory

abstract: Decision theorists propose a normative theory of rational choice. Traditionally, they assume that they should provide some constant and invariant principles as criteria for rational decisions, and indirectly, for agents. They seek a decision theory that invaribably works for all agents all the time. They believe that a rational agent should follow a certain principle, perhaps the principle of maximizing expected utility everywhere, all the time. As a result of the given context, these principles are considered, in this sense, context-independent. Furthermore, decision theorists usually assume that the relevant agents at work are ideal agents, and they believe that non-ideal agents should follow them so that their decisions qualify as rational. These principles are universal rules. I will refer to this context-independent and universal approach in traditional decision theory as Invariantism. This approach is, implicitly or explicitly, adopted by theories which are proposed on the basis of these two assumptions.


title: The convergence rate of of multivariate operators on simplex in Orlicz space

abstract: The approximation of functions in Orlicz space by multivariate operators on simplex is considered. The convergence rate is given by using modulus of smoothness.


title: Hadronic decays of Higgs boson at NNLO matched with parton shower

abstract: We present predictions for hadronic decays of the Higgs boson at next-to-next-to-leading order (NNLO) in QCD matched with parton shower based on the POWHEG framework. Those include decays into bottom quarks with full bottom-quark mass dependence, light quarks, and gluons in the heavy top quark effective theory. Our calculations describe exclusive decays of the Higgs boson with leading logarithmic accuracy in the Sudakov region and next-to-leading order (NLO) accuracy matched with parton shower in the three-jet region, with normalizations fixed to the partial width at NNLO. We estimated remaining perturbative uncertainties taking typical event shape variables as an example and demonstrated the need of future improvements on both parton shower and matrix element calculations. The calculations can be used immediately in evaluations of the physics performances of detector designs for future Higgs factories.


title: SGA: A Robust Algorithm for Partial Recovery of Tree-Structured Graphical Models with Noisy Samples

abstract: We consider learning Ising tree models when the observations from the nodes are corrupted by independent but non-identically distributed noise with unknown statistics. Katiyar et al. (2020) showed that although the exact tree structure cannot be recovered, one can recover a partial tree structure; that is, a structure belonging to the equivalence class containing the true tree. This paper presents a systematic improvement of Katiyar et al. (2020). First, we present a novel impossibility result by deriving a bound on the necessary number of samples for partial recovery. Second, we derive a significantly improved sample complexity result in which the dependence on the minimum correlation $\rho_{\min}$ is $\rho_{\min}^{-8}$ instead of $\rho_{\min}^{-24}$. Finally, we propose Symmetrized Geometric Averaging (SGA), a more statistically robust algorithm for partial tree recovery. We provide error exponent analyses and extensive numerical results on a variety of trees to show that the sample complexity of SGA is significantly better than the algorithm of Katiyar et al. (2020). SGA can be readily extended to Gaussian models and is shown via numerical experiments to be similarly superior.


title: Performance Analysis for Cache-enabled Cellular Networks with Cooperative Transmission

abstract: The large amount of deployed smart devices put tremendous traffic pressure on networks. Caching at the edge has been widely studied as a promising technique to solve this problem. To further improve the successful transmission probability (STP) of cache-enabled cellular networks (CEN), we combine the cooperative transmission technique with CEN and propose a novel transmission scheme. Local channel state information (CSI) is introduced at each cooperative base station (BS) to enhance the strength of the signal received by the user. A tight approximation for the STP of this scheme is derived using tools from stochastic geometry. The optimal content placement strategy of this scheme is obtained using a numerical method to maximize the STP. Simulation results demonstrate the optimal strategy achieves significant gains in STP over several comparative baselines with the proposed scheme.


title: Understanding the Tradeoffs in Client-side Privacy for Downstream Speech Tasks

abstract: As users increasingly rely on cloud-based computing services, it is important to ensure that uploaded speech data remains private. Existing solutions rely either on server-side methods or focus on hiding speaker identity. While these approaches reduce certain security concerns, they do not give users client-side control over whether their biometric information is sent to the server. In this paper, we formally define client-side privacy and discuss its three unique technical challenges: (1) direct manipulation of raw data on client devices, (2) adaptability with a broad range of server-side processing models, and (3) low time and space complexity for compatibility with limited-bandwidth devices. Solving these challenges requires new models that achieve high-fidelity reconstruction, privacy preservation of sensitive personal attributes, and efficiency during training and inference. As a step towards client-side privacy for speech recognition, we investigate three techniques spanning signal processing, disentangled representation learning, and adversarial training. Through a series of gender and accent masking tasks, we observe that each method has its unique strengths, but none manage to effectively balance the trade-offs between performance, privacy, and complexity. These insights call for more research in client-side privacy to ensure a safer deployment of cloud-based speech processing services.


title: High efficient multipartite entanglement purification using hyperentanglement

abstract: Multipartite entanglement plays an important role in controlled quantum teleportation, quantum secret sharing, quantum metrology and some other important quantum information branches. However, the maximally multipartite entangled state will degrade into the mixed state because of the noise. We present an efficient multipartite entanglement purification protocol (EPP) which can distill the high quality entangled states from low quality entangled states for N-photon systems in a Greenberger-Horne-Zeilinger (GHZ) state in only linear optics. After performing the protocol, the spatial-mode entanglement is used to purify the polarization entanglement and one pair of high quality polarization entangled state will be obtained. This EPP has several advantages. Firstly, with the same purification success probability, this EPP only requires one pair of multipartite GHZ state, while existing EPPs usually require two pairs of multipartite GHZ state. Secondly, if consider the practical transmission and detector efficiency, this EPP may be extremely useful for the ratio of purification efficiency is increased rapidly with both the number of photons and the transmission distance. Thirdly, this protocol requires linear optics and does not add additional measurement operations, so that it is feasible for experiment. All these advantages will make this protocol have potential application for future quantum information processing.


title: Long time asymptotics of large data in the Kadomtsev-Petviashvili models

abstract: We consider the Kadomtsev-Petviashvili (KP) equations posed on $\mathbb{R}^2$. For both equations, we provide sequential in time asymptotic descriptions of solutions, of arbitrarily large data, inside regions not containing lumps or line solitons, and under minimal regularity assumptions. The proof involves the introduction of two new virial identities adapted to the KP dynamics, showing decay in large regions of space, especially in the KP-I case, where no monotonicity property was previously known. Our results do not require the use of the integrability of KP and are adaptable to well-posed perturbations of KP.


title: How does COVID-19 change insurance and vaccine demand? Evidence from short-panel data in Japan

abstract: In this study, we explored how the coronavirus disease (COVID-19) affected the demand for insurance and vaccines in Japan from mid-March to mid-April 2020. Through independent internet surveys, respondents were asked hypothetical questions concerning the demand for insurance and vaccines for protection against COVID-19. Using the collected short-panel data, after controlling for individual characteristics using the fixed effects model, the key findings, within the context of the pandemic, were as follows: (1) Contrary to extant studies, the demand for insurance by females was smaller than that by their male counterparts; (2) The gap in demand for insurance between genders increased as the pandemic prevailed; (3) The demand for a vaccine by females was higher than that for males; and (4) As COVID-19 spread throughout Japan, demand for insurance decreased, whereas the demand for a vaccine increased.


title: Snapshot Hyperspectral Imaging Based on Weighted High-order Singular Value Regularization

abstract: Snapshot hyperspectral imaging can capture the 3D hyperspectral image (HSI) with a single 2D measurement and has attracted increasing attention recently. Recovering the underlying HSI from the compressive measurement is an ill-posed problem and exploiting the image prior is essential for solving this ill-posed problem. However, existing reconstruction methods always start from modeling image prior with the 1D vector or 2D matrix and cannot fully exploit the structurally spectral-spatial nature in 3D HSI, thus leading to a poor fidelity. In this paper, we propose an effective high-order tensor optimization based method to boost the reconstruction fidelity for snapshot hyperspectral imaging. We first build high-order tensors by exploiting the spatial-spectral correlation in HSI. Then, we propose a weight high-order singular value regularization (WHOSVR) based low-rank tensor recovery model to characterize the structure prior of HSI. By integrating the structure prior in WHOSVR with the system imaging process, we develop an optimization framework for HSI reconstruction, which is finally solved via the alternating minimization algorithm. Extensive experiments implemented on two representative systems demonstrate that our method outperforms state-of-the-art methods.


title: Numerical analysis of a deep learning formulation of elastic full waveform inversion with high order total variation regularization in different parameterization

abstract: We have formulated elastic seismic full waveform inversion (FWI) within a deep learning environment. In our formulation, a recurrent neural network is set up with rules enforcing elastic wave propagation, with the wavefield projected onto a measurement surface acting as the synthetic data to be compared with observed seismic data. Gradients for iterative updating of an elastic model, with a variety of parameterizations and misfit functionals, can be efficiently constructed within the network through the automatic differential method. With this method, the inversion based on complex misfits can be calculated. We survey the impact of different complex misfits (based on the l2, l1 ) with high order total variation (TV) regulations on multiparameter elastic FWI recovery of models within velocity/density, modulus/density, and stiffness parameter/density parameterizations. We analyze parameter cross-talk. Inversion results on simple and complex models show that the RNN elastic FWI with high order TV regulation using l1 norm can help mitigate cross-talk issues with gradient-based optimization methods.


title: Differentially Private SGD with Non-Smooth Losses

abstract: In this paper, we are concerned with differentially private {stochastic gradient descent (SGD)} algorithms in the setting of stochastic convex optimization (SCO). Most of the existing work requires the loss to be Lipschitz continuous and strongly smooth, and the model parameter to be uniformly bounded. However, these assumptions are restrictive as many popular losses violate these conditions including the hinge loss for SVM, the absolute loss in robust regression, and even the least square loss in an unbounded domain. We significantly relax these restrictive assumptions and establish privacy and generalization (utility) guarantees for private SGD algorithms using output and gradient perturbations associated with non-smooth convex losses. Specifically, the loss function is relaxed to have an $\alpha$-H\"{o}lder continuous gradient (referred to as $\alpha$-H\"{o}lder smoothness) which instantiates the Lipschitz continuity ($\alpha=0$) and the strong smoothness ($\alpha=1$). We prove that noisy SGD with $\alpha$-H\"older smooth losses using gradient perturbation can guarantee $(\epsilon,\delta)$-differential privacy (DP) and attain optimal excess population risk $\mathcal{O}\Big(\frac{\sqrt{d\log(1/\delta)}}{n\epsilon}+\frac{1}{\sqrt{n}}\Big)$, up to logarithmic terms, with the gradient complexity $ \mathcal{O}( n^{2-\alpha\over 1+\alpha}+ n).$ This shows an important trade-off between $\alpha$-H\"older smoothness of the loss and the computational complexity for private SGD with statistically optimal performance. In particular, our results indicate that $\alpha$-H\"older smoothness with $\alpha\ge {1/2}$ is sufficient to guarantee $(\epsilon,\delta)$-DP of noisy SGD algorithms while achieving optimal excess risk with the linear gradient complexity $\mathcal{O}(n).$


title: A Two-stream Neural Network for Pose-based Hand Gesture Recognition

abstract: Pose based hand gesture recognition has been widely studied in the recent years. Compared with full body action recognition, hand gesture involves joints that are more spatially closely distributed with stronger collaboration. This nature requires a different approach from action recognition to capturing the complex spatial features. Many gesture categories, such as "Grab" and "Pinch", have very similar motion or temporal patterns posing a challenge on temporal processing. To address these challenges, this paper proposes a two-stream neural network with one stream being a self-attention based graph convolutional network (SAGCN) extracting the short-term temporal information and hierarchical spatial information, and the other being a residual-connection enhanced bidirectional Independently Recurrent Neural Network (RBi-IndRNN) for extracting long-term temporal information. The self-attention based graph convolutional network has a dynamic self-attention mechanism to adaptively exploit the relationships of all hand joints in addition to the fixed topology and local feature extraction in the GCN. On the other hand, the residual-connection enhanced Bi-IndRNN extends an IndRNN with the capability of bidirectional processing for temporal modelling. The two streams are fused together for recognition. The Dynamic Hand Gesture dataset and First-Person Hand Action dataset are used to validate its effectiveness, and our method achieves state-of-the-art performance.


title: Single-shot Dispersion Sampling for Optical Pulse Reconstruction

abstract: We present a novel approach to single-shot characterization of the spectral phase of broadband laser pulses. Our method is inexpensive, insensitive to alignment and combines the simplicity and robustness of the dispersion scan technique, that does not require spatio-temporal pulse overlap, with the advantages of single-shot pulse characterization methods such as single-shot frequency-resolved optical gating at a real-time reconstruction rate of several Hz.


title: Machine Learning Percolation Model

abstract: Recent advances in machine learning have become increasingly popular in the applications of phase transitions and critical phenomena. By machine learning approaches, we try to identify the physical characteristics in the two-dimensional percolation model. To achieve this, we adopt Monte Carlo simulation to generate dataset at first, and then we employ several approaches to analyze the dataset. Four kinds of convolutional neural networks (CNNs), one variational autoencoder (VAE), one convolutional VAE (cVAE), one principal component analysis (PCA), and one $k$-means are used for identifying order parameter, the permeability, and the critical transition point. The former three kinds of CNNs can simulate the two order parameters and the permeability with high accuracy, and good extrapolating performance. The former two kinds of CNNs have high anti-noise ability. To validate the robustness of the former three kinds of CNNs, we also use the VAE and the cVAE to generate new percolating configurations to add perturbations into the raw configurations. We find that there is no difference by using the raw or the perturbed configurations to identify the physical characteristics, under the prerequisite of corresponding labels. In the case of lacking labels, we use unsupervised learning to detect the physical characteristics. The PCA, a classical unsupervised learning, performs well when identifying the permeability but fails to deduce order parameter. Hence, we apply the fourth kinds of CNNs with different preset thresholds, and identify a new order parameter and the critical transition point. Our findings indicate that the effectiveness of machine learning still needs to be evaluated in the applications of phase transitions and critical phenomena.


title: REPOSE: Distributed Top-k Trajectory Similarity Search with Local Reference Point Tries

abstract: Trajectory similarity computation is a fundamental component in a variety of real-world applications, such as ridesharing, road planning, and transportation optimization. Recent advances in mobile devices have enabled an unprecedented increase in the amount of available trajectory data such that efficient query processing can no longer be supported by a single machine. As a result, means of performing distributed in-memory trajectory similarity search are called for. However, existing distributed proposals suffer from either computing resource waste or are unable to support the range of similarity measures that are being used. We propose a distributed in-memory management framework called REPOSE for processing top-k trajectory similarity queries on Spark. We develop a reference point trie (RP-Trie) index to organize trajectory data for local search. In addition, we design a novel heterogeneous global partitioning strategy to eliminate load imbalance in distributed settings. We report on extensive experiments with real-world data that offer insight into the performance of the solution, and show that the solution is capable of outperforming the state-of-the-art proposals.


title: A practical method for recovering Sturm-Liouville problems from the Weyl function

abstract: In the paper we propose a direct method for recovering the Sturm-Liouville potential from the Weyl-Titchmarsh $m$-function given on a countable set of points. We show that using the Fourier-Legendre series expansion of the transmutation operator integral kernel the problem reduces to an infinite linear system of equations, which is uniquely solvable if so is the original problem. The solution of this linear system allows one to reconstruct the characteristic determinant and hence to obtain the eigenvalues as its zeros and to compute the corresponding norming constants. As a result, the original inverse problem is transformed to an inverse problem with a given spectral density function, for which the direct method of solution from arXiv:2010.15275 is applied. The proposed method leads to an efficient numerical algorithm for solving a variety of inverse problems. In particular, the problems in which two spectra or some parts of three or more spectra are given, the problems in which the eigenvalues depend on a variable boundary parameter (including spectral parameter dependent boundary conditions), problems with a partially known potential and partial inverse problems on quantum graphs.


title: Rationality of even-dimensional intersections of two real quadrics

abstract: We study rationality constructions for smooth complete intersections of two quadrics over nonclosed fields. Over the real numbers, we establish a criterion for rationality in dimension four.


title: Sobolev Training for Physics Informed Neural Networks

abstract: Physics Informed Neural Networks (PINNs) is a promising application of deep learning. The smooth architecture of a fully connected neural network is appropriate for finding the solutions of PDEs; the corresponding loss function can also be intuitively designed and guarantees the convergence for various kinds of PDEs. However, the rate of convergence has been considered as a weakness of this approach. This paper proposes Sobolev-PINNs, a novel loss function for the training of PINNs, making the training substantially efficient. Inspired by the recent studies that incorporate derivative information for the training of neural networks, we develop a loss function that guides a neural network to reduce the error in the corresponding Sobolev space. Surprisingly, a simple modification of the loss function can make the training process similar to \textit{Sobolev Training} although PINNs is not a fully supervised learning task. We provide several theoretical justifications that the proposed loss functions upper bound the error in the corresponding Sobolev spaces for the viscous Burgers equation and the kinetic Fokker--Planck equation. We also present several simulation results, which show that compared with the traditional $L^2$ loss function, the proposed loss function guides the neural network to a significantly faster convergence. Moreover, we provide the empirical evidence that shows that the proposed loss function, together with the iterative sampling techniques, performs better in solving high dimensional PDEs.


title: Inferring solar differential rotation through normal-mode coupling using Bayesian statistics

abstract: Normal-mode helioseismic data analysis uses observed solar oscillation spectra to infer perturbations in the solar interior due to global and local-scale flows and structural asphericity. Differential rotation, the dominant global-scale axisymmetric perturbation, has been tightly constrained primarily using measurements of frequency splittings via "$a$-coefficients". However, the frequency-splitting formalism invokes the approximation that multiplets are isolated. This assumption is inaccurate for modes at high angular degrees. Analysing eigenfunction corrections, which respect cross coupling of modes across multiplets, is a more accurate approach. However, applying standard inversion techniques using these cross-spectral measurements yields $a$-coefficients with a significantly wider spread than the well-constrained results from frequency splittings. In this study, we apply Bayesian statistics to infer $a$-coefficients due to differential rotation from cross spectra for both $f$-modes and $p$-modes. We demonstrate that this technique works reasonably well for modes with angular degrees $\ell=50-291$. The inferred $a_3-$coefficients are found to be within $1$ nHz of the frequency splitting values for $\ell > 200$. We also show that the technique fails at $\ell < 50$ owing to the insensitivity of the measurement to the perturbation. These results serve to further establish mode coupling as an important helioseismic technique with which to infer internal structure and dynamics, both axisymmetric (e.g., meridional circulation) and non-axisymmetric perturbations.


title: AS-Net: Fast Photoacoustic Reconstruction with Multi-feature Fusion from Sparse Data

abstract: Photoacoustic (PA) imaging is a biomedical imaging modality capable of acquiring high-contrast images of optical absorption at depths much greater than traditional optical imaging techniques. However, practical instrumentation and geometry limit the number of available acoustic sensors surrounding the imaging target, which results in the sparsity of sensor data. Conventional PA image reconstruction methods give severe artifacts when they are applied directly to the sparse PA data. In this paper, we firstly propose to employ a novel signal processing method to make sparse PA raw data more suitable for the neural network, concurrently speeding up image reconstruction. Then we propose Attention Steered Network (AS-Net) for PA reconstruction with multi-feature fusion. AS-Net is validated on different datasets, including simulated photoacoustic data from fundus vasculature phantoms and experimental data from in vivo fish and mice. Notably, the method is also able to eliminate some artifacts present in the ground truth for in vivo data. Results demonstrated that our method provides superior reconstructions at a faster speed.


title: Direct and inverse scattering problems for the first-order discrete system associated with the derivative NLS system

abstract: The direct and inverse scattering problems are analyzed for a first-order discrete system associated with the semi-discrete version of the derivative NLS system. The Jost solutions, the scattering coefficients, the bound-state dependency and norming constants are investigated and related to the corresponding quantities for two particular discrete linear systems associated with the semi-discrete version of the NLS system. The bound-state data set with any multiplicities is described in an elegant manner in terms of a pair of constant matrix triplets. Several methods are presented to the solve the inverse problem. One of these methods involves a discrete Marchenko system using as input the scattering data set consisting of the scattering coefficients and the bound-state information, and this method is presented in a way generalizable to other first-order systems both in the discrete and continuous cases for which a Marchenko system is not yet available. Finally, using the time-evolved scattering data set, the inverse scattering transform is applied on the corresponding semi-discrete derivative NLS system, and in the reflectionless case certain explicit solution formulas are presented in closed form expressed in terms of the two matrix triplets.


title: Numerical methods for backward stochastic differential equations: A survey

abstract: Backward Stochastic Differential Equations (BSDEs) have been widely employed in various areas of social and natural sciences, such as the pricing and hedging of financial derivatives, stochastic optimal control problems, optimal stopping problems and gene expression. Most BSDEs cannot be solved analytically and thus numerical methods must be applied to approximate their solutions. There have been a variety of numerical methods proposed over the past few decades as well as many more currently being developed. For the most part, they exist in a complex and scattered manner with each requiring a variety of assumptions and conditions. The aim of the present work is thus to systematically survey various numerical methods for BSDEs, and in particular, compare and categorize them, for further developments and improvements. To achieve this goal, we focus primarily on the core features of each method based on an extensive collection of 333 references: the main assumptions, the numerical algorithm itself, key convergence properties and advantages and disadvantages, to provide an up-to-date coverage of numerical methods for BSDEs, with insightful summaries of each and a useful comparison and categorization.


title: Prior Preference Learning from Experts:Designing a Reward with Active Inference

abstract: Active inference may be defined as Bayesian modeling of a brain with a biologically plausible model of the agent. Its primary idea relies on the free energy principle and the prior preference of the agent. An agent will choose an action that leads to its prior preference for a future observation. In this paper, we claim that active inference can be interpreted using reinforcement learning (RL) algorithms and find a theoretical connection between them. We extend the concept of expected free energy (EFE), which is a core quantity in active inference, and claim that EFE can be treated as a negative value function. Motivated by the concept of prior preference and a theoretical connection, we propose a simple but novel method for learning a prior preference from experts. This illustrates that the problem with inverse RL can be approached with a new perspective of active inference. Experimental results of prior preference learning show the possibility of active inference with EFE-based rewards and its application to an inverse RL problem.


title: Spin Hall effect of light under arbitrarily polarized and unpolarized light

abstract: The spin Hall effect of light (SHEL), which refers to a spin-dependent and transverse splitting at refraction and reflection phenomena, inherently depends on the polarization states of the incidence. Most of the previous research have focused on a horizontally or vertically polarized incidence, in which the analytic formula of the shift is well-formulated and SHEL appears symmetrically in both shift and intensity. However, the SHEL under an arbitrarily polarized or unpolarized incidence has remained largely unexplored. Whereas the SHEL under other polarization is sensitive to incident polarization and is asymmetrical, here we demonstrate that the SHEL is independent of the incident polarization and is symmetrical in shift if Fresnel coefficients of the two linear polarization are the same. The independence of the shift with respect to the incident polarization is proved both analytically and numerically. Moreover, we prove that under an unpolarized incidence composed of a large number of completely random polarization states, the reflected beam is split in half into two circularly polarized components that undergo the same amount of splitting but in opposite directions. This result means that under unpolarized incidence, the SHEL occurs exactly the same as under horizontally or vertically polarized incidence. We believe that the incident-polarization-independent SHEL can broaden the applicability of the SHEL to cover optical systems in which polarization is ill-defined.


title: A Rich Type System for Quantum Programs

abstract: We show that Gottesman's semantics (GROUP22, 1998) for Clifford circuits based on the Heisenberg representation can be treated as a type system that can efficiently characterize a common subset of quantum programs. Our applications include (i) certifying whether auxiliary qubits can be safely disposed of, (ii) determining if a system is separable across a given bi-partition, (iii) checking the transversality of a gate with respect to a given stabilizer code, and (iv) typing post-measurement states for computational basis measurements. Further, this type system is extended to accommodate universal quantum computing by deriving types for the $T$-gate, multiply-controlled unitaries such as the Toffoli gate, and some gate injection circuits that use associated magic states. These types allow us to prove a lower bound on the number of $T$ gates necessary to perform a multiply-controlled $Z$ gate.


title: Hessian-Aware Pruning and Optimal Neural Implant

abstract: Pruning is an effective method to reduce the memory footprint and FLOPs associated with neural network models. However, existing structured-pruning methods often result in significant accuracy degradation for moderate pruning levels. To address this problem, we introduce a new Hessian Aware Pruning (HAP) method coupled with a Neural Implant approach that uses second-order sensitivity as a metric for structured pruning. The basic idea is to prune insensitive components and to use a Neural Implant for moderately sensitive components, instead of completely pruning them. For the latter approach, the moderately sensitive components are replaced with with a low rank implant that is smaller and less computationally expensive than the original component. We use the relative Hessian trace to measure sensitivity, as opposed to the magnitude based sensitivity metric commonly used in the literature. We test HAP for both computer vision tasks and natural language tasks, and we achieve new state-of-the-art results. Specifically, HAP achieves less than $0.1\%$/$0.5\%$ degradation on PreResNet29/ResNet50 (CIFAR-10/ImageNet) with more than 70\%/50\% of parameters pruned. Meanwhile, HAP also achieves significantly better performance (up to 0.8\% with 60\% of parameters pruned) as compared to gradient based method for head pruning on transformer-based models. The framework has been open sourced and available online.


title: Optimising optical tweezers experiments for magnetic resonance sensing with nanodiamonds

abstract: In this article we explore the requirements for enabling high quality optically detected magnetic resonance (ODMR) spectroscopy in a conventional gradient force optical tweezers using nanodiamonds containing nitrogen-vacancy (NV$^{-}$) centres. We find that modulation of the infrared (1064 nm) trapping laser during spectroscopic measurements substantially improves the ODMR contrast compared with continuous wave trapping. The work is significant as it allows trapping and quantum sensing protocols to be performed in conditions where signal contrast is substantially reduced. We demonstrate the utility of the technique by resolving NV$^{-}$ spin projections within the ODMR spectrum. Manipulating the orientation of the nanodiamond via the trapping laser polarisation, we observe changes in spectral features. Theoretical modelling then allows us to infer the crystallographic orientation of the NV$^{-}$. This is an essential capability for future magnetic sensing applications of optically trapped nanodiamonds.


title: Enriching Non-Autoregressive Transformer with Syntactic and SemanticStructures for Neural Machine Translation

abstract: The non-autoregressive models have boosted the efficiency of neural machine translation through parallelized decoding at the cost of effectiveness when comparing with the autoregressive counterparts. In this paper, we claim that the syntactic and semantic structures among natural language are critical for non-autoregressive machine translation and can further improve the performance. However, these structures are rarely considered in the existing non-autoregressive models. Inspired by this intuition, we propose to incorporate the explicit syntactic and semantic structures of languages into a non-autoregressive Transformer, for the task of neural machine translation. Moreover, we also consider the intermediate latent alignment within target sentences to better learn the long-term token dependencies. Experimental results on two real-world datasets (i.e., WMT14 En-De and WMT16 En-Ro) show that our model achieves a significantly faster speed, as well as keeps the translation quality when compared with several state-of-the-art non-autoregressive models.


title: Binary Polar Codes Based on Bit Error Probability

abstract: This paper introduces techniques to construct binary polar source/channel codes based on the bit error probability of successive-cancellation decoding. The polarization lemma is reconstructed based on the bit error probability and then techniques to compute the bit error probability are introduced. These techniques can be applied to the construction of polar codes and the computation of lower and upper bounds of the block decoding error probability.


title: Learning to Simulate High Energy Particle Collisions from Unlabeled Data

abstract: In many scientific fields which rely on statistical inference, simulations are often used to map from theoretical models to experimental data, allowing scientists to test model predictions against experimental results. Experimental data is often reconstructed from indirect measurements causing the aggregate transformation from theoretical models to experimental data to be poorly-described analytically. Instead, numerical simulations are used at great computational cost. We introduce Optimal-Transport-based Unfolding and Simulation (OTUS), a fast simulator based on unsupervised machine-learning that is capable of predicting experimental data from theoretical models. Without the aid of current simulation information, OTUS trains a probabilistic autoencoder to transform directly between theoretical models and experimental data. Identifying the probabilistic autoencoder's latent space with the space of theoretical models causes the decoder network to become a fast, predictive simulator with the potential to replace current, computationally-costly simulators. Here, we provide proof-of-principle results on two particle physics examples, $Z$-boson and top-quark decays, but stress that OTUS can be widely applied to other fields.


title: HAWC Search for High-Mass Microquasars

abstract: Microquasars with high-mass companion stars are promising very-high-energy (VHE; 0.1-100 TeV) gamma-ray emitters, but their behaviors above 10 TeV are poorly known. Using the High Altitude Water Cherenkov (HAWC) observatory, we search for excess gamma-ray emission coincident with the positions of known high-mass microquasars (HMMQs). No significant emission is observed for LS 5039, Cygnus X-1, Cygnus X-3, and SS 433 with 1,523 days of HAWC data. We set the most stringent limit above 10 TeV obtained to date on each individual source. Under the assumption that HMMQs produce gamma rays via a common mechanism, we have performed source-stacking searches, considering two different scenarios: I) gamma-ray luminosity is a fraction $\epsilon_\gamma$ of the microquasar jet luminosity, and II) very-high-energy gamma rays are produced by relativistic electrons up-scattering the radiation field of the companion star in a magnetic field $B$. We obtain $\epsilon_\gamma < 5.4\times 10^{-6}$ for scenario I, which tightly constrains models that suggest observable high-energy neutrino emission by HMMQs. In the case of scenario II, the non-detection of VHE gamma rays yields a strong magnetic field, which challenges synchrotron radiation as the dominant mechanism of the microquasar emission between 10 keV and 10 MeV.


title: Generation and verification of 27-qubit Greenberger-Horne-Zeilinger states in a superconducting quantum computer

abstract: Generating and detecting genuine multipartite entanglement (GME) of sizeable quantum states prepared on physical devices is an important benchmark for highlighting the progress of near-term quantum computers. A common approach to certify GME is to prepare a Greenberger-Horne-Zeilinger (GHZ) state and measure a GHZ fidelity of at least 0.5. We measure the fidelities using multiple quantum coherences of GHZ states on 11 to 27 qubits prepared on the IBM Quantum ibmq_montreal device. Combinations of quantum readout error mitigation (QREM) and parity verification error detection are applied to the states. A fidelity of $0.546 \pm 0.017$ was recorded for a 27-qubit GHZ state when QREM was used, demonstrating GME across the full device with a confidence level of 98.6%. We benchmarked the effect of parity verification on GHZ fidelity for two GHZ state preparation embeddings on the heavy-hexagon architecture. The results show that the effect of parity verification, while relatively modest, led to a detectable improvement of GHZ fidelity.


title: On Optimization of 1/2-Approximation Path Cover Algorithm

abstract: In this paper, we propose a deterministic algorithm that approximates the optimal path cover on weighted undirected graphs. Based on the 1/2-Approximation Path Cover Algorithm by Moran et al., we add a procedure to remove the redundant edges as the algorithm progresses. Our optimized algorithm not only significantly reduces the computation time but also maintains the theoretical guarantee of the original 1/2-Approximation Path Cover Algorithm. To test the time complexity, we conduct numerical tests on graphs with various structures and random weights, from structured ring graphs to random graphs, such as Erdos-Renyi graphs. The tests demonstrate the effectiveness of our proposed algorithm on graphs, especially those with high degree nodes, and the advantages expand as the graph gets larger. Moreover, we also launch tests on various graphs/networks derived from a wide range of real-world problems to suggest the effectiveness and applicability of the proposed algorithm.


title: Spin dynamics in NaFeAs and NaFe$_{0.53}$Cu$_{0.47}$As probed by resonant inelastic X-ray scattering

abstract: The parent compounds of iron-based superconductors are magnetically-ordered bad metals, with superconductivity appearing near a putative magnetic quantum critical point. The presence of both Hubbard repulsion and Hund's coupling leads to rich physics in these multiorbital systems, and motivated descriptions of magnetism in terms of itinerant electrons or localized spins. The NaFe$_{1-x}$Cu$_x$As series consists of magnetically-ordered bad metal ($x=0$), superconducting ($x\approx0.02$) and magnetically-ordered semiconducing/insulating ($x\approx0.5$) phases, providing a platform to investigate the connection between superconductivity, magnetism and electronic correlations. Here we use X-ray absorption spectroscopy and resonant inelastic X-ray scattering to study the valence state of Fe and spin dynamics in two NaFe$_{1-x}$Cu$_x$As compounds ($x=0$ and 0.47). We find that magnetism in both compounds arises from Fe$^{2+}$ atoms, and exhibits underdamped dispersive spin waves in their respective ordered states. The dispersion of spin excitations in NaFe$_{0.53}$Cu$_{0.47}$As is consistent with being quasi-one-dimensional. Compared to NaFeAs, the band top of spin waves in NaFe$_{0.53}$Cu$_{0.47}$As is slightly softened with significantly more spectral weight of the spin excitations. Our results indicate the spin dynamics in NaFe$_{0.53}$Cu$_{0.47}$As arise from localized magnetic moments and suggest the iron-based superconductors are proximate to a correlated insulating state with localized iron moments.


title: Strongly coupled Yukawa plasma layer in a harmonic trap

abstract: Observations made in dusty plasma experiments suggest that an ensemble of electrically charged solid particles, confined in an elongated trap, develops structural inhomogeneities. With narrowing the trap the particles tend to form layers oriented parallel with the trap walls. In this work we present theoretical and numerical results on the structure of three-dimensional many-particle systems with screened Coulomb (Yukawa) inter-particle interaction in the strongly coupled liquid phase, confined in one-dimensional harmonic trap, forming quasi-2D configurations. Particle density profiles are calculated by means of the hypernetted chain approximation (HNC), showing clear signs of layer formation. The mechanism behind the formation of layer structure is discussed and a method to predict the number of layers is presented. Molecular dynamics (MD) simulations provide validation of the theoretical results and detailed microscopic insights.


title: Implications of High Polarization Degree for the Surface State of Ryugu

abstract: The asteroid exploration project "Hayabusa2" has successfully returned samples from the asteroid (162173) Ryugu. In this study, we measured the linear polarization degrees of Ryugu using four ground-based telescopes from 2020 September 27 to December 25, covering a wide-phase angle (Sun-target-observer's angle) range from 28$^\circ$ to 104$^\circ$. We found that the polarization degree of Ryugu reached 53$\%$ around a phase angle of 100$^\circ$, the highest value among all asteroids and comets thus far reported. The high polarization degree of Ryugu can be attributed to the scattering properties of its surface layers, in particular the relatively small contribution of multiply-scattered light. Our polarimetric results indicate that Ryugu's surface is covered with large grains. On the basis of a comparison with polarimetric measurements of pulverized meteorites, we can infer the presence of submillimeter-sized grains on the surface layer of Ryugu. We also conjecture that this size boundary represents the grains that compose the aggregate. It is likely that a very brittle structure has been lost in the recovered samples, although they may hold a record of its evolution. Our data will be invaluable for future experiments aimed at reproducing the surface structure of Ryugu.