title: Well-resolved velocity fields using discontinuous Galerkin shallow water solutions

abstract: Computational models based on the depth-averaged shallow water equations (SWE) offer an efficient choice to analyse velocity fields around hydraulic structures. Second-order finite volume (FV2) solvers have often been used for this purpose subject to adding an eddy viscosity term at sub-meter resolution, but have been shown to fall short of capturing small-scale field transients emerging from wave-structure interactions. The second-order discontinuous Galerkin (DG2) alternative is significantly more resistant to the growth of numerical diffusion and leads to faster convergence rates. These properties make the DG2 solver a promising modelling tool for detailed velocity field predictions. This paper focuses on exploring this DG2 capability with reference to an FV2 counterpart for a selection of test cases that require well-resolved velocity field predictions. The findings of this work lead to identifying a particular setting for the DG2 solver that allows for obtaining more accurate and efficient depth-averaged velocity fields incorporating small-scale transients.


title: Hadronization studies at the LHC with ALICE

abstract: Studies of the production of heavy-flavour baryons are of prominent importance to investigate hadronization mechanisms at the LHC, in particular through the study of the evolution of the baryon-over-meson production ratio. Measurements performed in pp and p--Pb collisions at the LHC have revealed unexpected features, qualitatively similar to what was observed in heavy-ion collisions and, in the charm sector, not in line with the expectations based on previous measurements from $\rm e^+e^-$ colliders and in ep collisions. These results suggest that charmed baryon formation might not be universal and that the baryon-over-meson ratio depends on the collision system or multiplicity. A review of ALICE measurements of charmed baryons, including $\rm \Lambda_c^+/D^0$ as a function of charged-particle multiplicity in pp, p--Pb and Pb--Pb collisions, $\rm \Sigma_c^{0, +, ++}/D^0$ and $\rm \Xi_c^{0, +}/D^0$ as a function of $p_{\rm T}$ in pp collisions and $\rm \Gamma(\Xi_c^0\rightarrow\Xi^-e^+\nu_e)/\Gamma(\Xi_c^0\rightarrow\Xi^-\pi^+)$, will be presented. Comparison to phenomenological models will be also discussed. Emphasis will be given to the discussion of the impact of these studies on the understanding of hadronization processes.


title: A quiver invariant theoretic approach to Radial Isotropy and Paulsen's Problem for matrix frames

abstract: In this paper, we view matrix frames as representations of quivers and study them within the general framework of quiver invariant theory. We are thus led to consider the large class of semi-stable matrix frames. Within this class, we are particularly interested in radial isotropic and Parseval matrix frames. Using methods from quiver invariant theory [CD19], we first prove a far reaching generalization of Barthe's Radial Isotropy Theorem [Bar98] to matrix frames (see Theorems 1(3) and 29). With this tool at our disposal, we provide a quiver invariant theoretic approach to the Paulsen problem for matrix frames. We show in Theorem 2 that for any given $\varepsilon$-nearly equal-norm Parseval frame $\mathcal{F}$ of $n$ matrices with $d$ rows there exists an equal-norm Parseval frame $\mathcal{W}$ of $n$ matrices with $d$ rows such that $\mathsf{dist}^2(\mathcal{F},\mathcal{W})\leq 46 \epsilon d^2$.


title: A QUBO Algorithm to Compute Eigenvectors of Symmetric Matrices

abstract: We describe an algorithm to compute the extremal eigenvalues and corresponding eigenvectors of a symmetric matrix by solving a sequence of Quadratic Binary Optimization problems. This algorithm is robust across many different classes of symmetric matrices, can compute the eigenvector/eigenvalue pair to essentially arbitrary precision, and with minor modifications can also solve the generalized eigenvalue problem. Performance is analyzed on small random matrices and selected larger matrices from practical applications.


title: Building Load Control using Distributionally Robust Chance-Constrained Programs with Right-Hand Side Uncertainty and the Risk-Adjustable Variants

abstract: Aggregation of heating, ventilation, and air conditioning (HVAC) loads can provide reserves to absorb volatile renewable energy, especially solar photo-voltaic (PV) generation. In this paper, we decide HVAC control schedules under uncertain PV generation, using a distributionally robust chance-constrained (DRCC) building load control model under two typical ambiguity sets: the moment-based and Wasserstein ambiguity sets. We derive mixed-integer linear programming (MILP) reformulations for DRCC problems under both sets. Especially, for the Wasserstein ambiguity set, we utilize the right-hand side (RHS) uncertainty to derive a more compact MILP reformulation than the commonly known MILP reformulations with big-M constants. All the results also apply to general individual chance constraints with RHS uncertainty. Furthermore, we propose an adjustable chance-constrained variant to achieve a trade-off between the operational risk and costs. We derive MILP reformulations under the Wasserstein ambiguity set and second-order conic programming (SOCP) reformulations under the moment-based set. Using real-world data, we conduct computational studies to demonstrate the efficiency of the solution approaches and the effectiveness of the solutions.


title: Resonant Electro-optic Imaging for Microscopy at Nanosecond Resolution

abstract: We demonstrate an electro-optic wide-field method to enable fluorescence lifetime microscopy (FLIM) with high throughput and single-molecule sensitivity. Resonantly driven Pockels cells are used to efficiently gate images at 39 MHz, allowing fluorescence lifetime to be captured on standard camera sensors. Lifetime imaging of single molecules is enabled in wide-field with exposure times of less than 100 milliseconds. This capability allows combination of wide-field FLIM with single-molecule super-resolution localization microscopy. Fast single-molecule dynamics such as FRET and molecular binding events are captured from wide-field images without prior spatial knowledge. A lifetime sensitivity of 1.9 times the photon shot-noise limit is achieved, and high throughput is shown by acquiring wide-field FLIM images with millisecond exposure and $>10^8$ photons per frame. Resonant electro-optic FLIM allows lifetime contrast in any wide-field microscopy method.


title: Ordered intricacy of Shilnikov saddle-focus homoclinics in symmetric systems

abstract: Using the technique of Poincar\'{e} return maps, we disclose an intricate order of the subsequent homoclinics near the primary homoclinic bifurcation of the Shilnikov saddle-focus in systems with reflection symmetry. We also reveal the admissible shapes of the corresponding bifurcation curves in a parameter plane of such systems. The scalability ratio of geometry and organization is proven to be universal for such homoclinic bifurcations of higher orders. Two applications with similar dynamics due to the Shilnikov saddle-foci, a smooth adaptation of the Chua circuit and a 3D normal form, are used to illustrate the theory


title: SPECTRE: Defending Against Backdoor Attacks Using Robust Statistics

abstract: Modern machine learning increasingly requires training on a large collection of data from multiple sources, not all of which can be trusted. A particularly concerning scenario is when a small fraction of poisoned data changes the behavior of the trained model when triggered by an attacker-specified watermark. Such a compromised model will be deployed unnoticed as the model is accurate otherwise. There have been promising attempts to use the intermediate representations of such a model to separate corrupted examples from clean ones. However, these defenses work only when a certain spectral signature of the poisoned examples is large enough for detection. There is a wide range of attacks that cannot be protected against by the existing defenses. We propose a novel defense algorithm using robust covariance estimation to amplify the spectral signature of corrupted data. This defense provides a clean model, completely removing the backdoor, even in regimes where previous methods have no hope of detecting the poisoned examples. Code and pre-trained models are available at https://github.com/SewoongLab/spectre-defense .


title: Real-time NLOS/LOS Identification for Smartphone-based Indoor Positioning System using WiFi RTT and RSS

abstract: The accuracy of smartphone-based positioning methods using WiFi usually suffers from ranging errors caused by non-line-of-sight (NLOS) conditions. Previous research usually exploits several statistical features from a long time series (hundreds of samples) of WiFi received signal strength (RSS) or WiFi round-trip time (RTT) to achieve a high identification accuracy. However, the long time series or large sample size attributes to high power and time consumption in data collection for both training and testing. This will also undoubtedly be detrimental to user experience as the waiting time of getting enough samples is quite long. Therefore, this paper proposes a new real-time NLOS/LOS identification method for smartphone-based indoor positioning system using WiFi RTT and RSS. Based on our extensive analysis of RSS and RTT features, a machine learning-based method using random forest was chosen and developed to separate the samples for NLOS/LOS conditions. Experiments in different environments show that our method achieves a discrimination accuracy of about 94% with a sample size of 10. Considering the theoretically shortest WiFi ranging interval of 100ms of the RTT-enabled smartphones, our algorithm is able to provide the shortest latency of 1s to get the testing result among all of the state-of-art methods.


title: Improving Hierarchy Storage for Video Streaming in Cloud

abstract: Frequently accessed video streams are pre-transcoded into several formats to satisfy the characteristics of all display devices. Storing several video stream formats imposes a high cost on video stream providers using the old classical way. Alternatively, cloud providers offer a high flexibility of using their services and at a low cost relatively. Therefore, video stream companies adopted cloud technology to store their video streams. Generally, having all video streams stored in one type of cloud storage, the cost rises gradually. More importantly, the variation of the access pattern to frequently accessed video streams impacts negatively the storage cost and increases it significantly. To optimize storage usage and lower its cost, we propose a method that manages the cloud hierarchy storage. Particularly, we develop an algorithm that operates on parts of different videos that are frequently accessed and stores them in their suitable storage type cloud. Experiments came up with promising results on reducing the cost of using cloud storage by 18.75 %.


title: Improved Gamma-Ray Point Source Quantification in Three Dimensions by Modeling Attenuation in the Scene

abstract: Using a series of detector measurements taken at different locations to localize a source of radiation is a well-studied problem. The source of radiation is sometimes constrained to a single point-like source, in which case the location of the point source can be found using techniques such as maximum likelihood. Recent advancements have shown the ability to locate point sources in 2D and even 3D, but few have studied the effect of intervening material on the problem. In this work we examine gamma-ray data taken from a freely moving system and develop voxelized 3-D models of the scene using data from the onboard LiDAR. Ray casting is used to compute the distance each gamma ray travels through the scene material, which is then used to calculate attenuation assuming a single attenuation coefficient for solids within the geometry. Parameter estimation using maximum likelihood is performed to simultaneously find the attenuation coefficient, source activity, and source position that best match the data. Using a simulation, we validate the ability of this method to reconstruct the true location and activity of a source, along with the true attenuation coefficient of the structure it is inside, and then we apply the method to measured data with sources and find good agreement.


title: Galaxy-galaxy lensing with the DES-CMASS catalogue: measurement and constraints on the galaxy-matter cross-correlation

abstract: The DMASS sample is a photometric sample from the DES Year 1 data set designed to replicate the properties of the CMASS sample from BOSS, in support of a joint analysis of DES and BOSS beyond the small overlapping area. In this paper, we present the measurement of galaxy-galaxy lensing using the DMASS sample as gravitational lenses in the DES Y1 imaging data. We test a number of potential systematics that can bias the galaxy-galaxy lensing signal, including those from shear estimation, photometric redshifts, and observing conditions. After careful systematic tests, we obtain a highly significant detection of the galaxy-galaxy lensing signal, with total $S/N=25.7$. With the measured signal, we assess the feasibility of using DMASS as gravitational lenses equivalent to CMASS, by estimating the galaxy-matter cross-correlation coefficient $r_{\rm cc}$. By jointly fitting the galaxy-galaxy lensing measurement with the galaxy clustering measurement from CMASS, we obtain $r_{\rm cc}=1.09^{+0.12}_{-0.11}$ for the scale cut of $4~h^{-1}{\rm Mpc}$ and $r_{\rm cc}=1.06^{+0.13}_{-0.12}$ for $12~h^{-1}{\rm Mpc}$ in fixed cosmology. By adding the angular galaxy clustering of DMASS, we obtain $r_{\rm cc}=1.06\pm 0.10$ for the scale cut of $4~h^{-1}{\rm Mpc}$ and $r_{\rm cc}=1.03\pm 0.11$ for $12~h^{-1}{\rm Mpc}$. The resulting values of $r_{\rm cc}$ indicate that the lensing signal of DMASS is statistically consistent with the one that would have been measured if CMASS had populated the DES region within the given statistical uncertainty. The measurement of galaxy-galaxy lensing presented in this paper will serve as part of the data vector for the forthcoming cosmology analysis in preparation.


title: Federated Double Deep Q-learning for Joint Delay and Energy Minimization in IoT networks

abstract: In this paper, we propose a federated deep reinforcement learning framework to solve a multi-objective optimization problem, where we consider minimizing the expected long-term task completion delay and energy consumption of IoT devices. This is done by optimizing offloading decisions, computation resource allocation, and transmit power allocation. Since the formulated problem is a mixed-integer non-linear programming (MINLP), we first cast our problem as a multi-agent distributed deep reinforcement learning (DRL) problem and address it using double deep Q-network (DDQN), where the actions are offloading decisions. The immediate cost of each agent is calculated through solving either the transmit power optimization or local computation resource optimization, based on the selected offloading decisions (actions). Then, to enhance the learning speed of IoT devices (agents), we incorporate federated learning (FDL) at the end of each episode. FDL enhances the scalability of the proposed DRL framework, creates a context for cooperation between agents, and minimizes their privacy concerns. Our numerical results demonstrate the efficacy of our proposed federated DDQN framework in terms of learning speed compared to federated deep Q network (DQN) and non-federated DDQN algorithms. In addition, we investigate the impact of batch size, network layers, DDQN target network update frequency on the learning speed of the FDL.


title: JRAF: A Julia Package for Computation of the Relativistic Molecular Auxiliary Functions

abstract: Evaluation of relativistic molecular integrals over exponential-type spinor orbitals require using the relativistic auxiliary functions in prolate spheroidal coordinates. They have derived recently by the author [Physical Review E 91, 023303 (2015)]. They are used in solution of the molecular Dirac equation for electrons moving around Coulomb potential. A series of papers on a method for fully analytical evaluation of relativistic auxiliary functions in following, published [2, 3, 4]. From the computational physics point of view, these works also demonstrate how to deal with the integrals involve product of power functions with non-integer exponents and incomplete gamma functions. The computer program package to calculate these auxiliary functions in high accuracy is presented. It is designed in Julia programming language. It is capable of yielding highly accurate results for molecular integrals over a wide range of orbital parameters and quantum numbers. Additionally, the program package provides numerous tools such as efficient calculation for the angular momentum coefficients arising in the product of two normalized associated Legendre functions centered on different atomic positions, the rotation angular functions used for both complex and real spherical harmonics. Sample calculations are performed for two-center one-electron integrals over non-integer Slater-type orbitals. The results prove that the package is robust.


title: Analytical solution of the cylindrical torsion problem for the relaxed micromorphic continuum and other generalized continua (including full derivations)

abstract: We solve the St.Venant torsion problem for an infinite cylindrical rod whose behaviour is described by a family of isotropic generalized continua, including the relaxed micromorphic and classical micromorphic model. The results can be used to determine the material parameters of these models. Special attention is given to the possible nonphysical stiffness singularity for a vanishing rod diameter, since slender specimens are in general described as stiffer.


title: Matlab code for Lyapunov exponents of fractional-order systems, Part II: The non-commensurate case

abstract: In this paper the Benettin-Wolf algorithm to determine all Lyapunov exponents adapted to a class of non-commensurate fractional-order systems modeled by Caputo's derivative and the corresponding Matlab code are presented. The paper continues the work started in [Danca & Kuznetsov, 2018], where the Matlab code of commensurate fractional-order systems is given. To integrate the extended system, the Adams-Bashforth-Moulton for fractional differential equations is utilized. As the the Matlab program for commensurate-order systems, the program presented in this paper prints and plots all Lyapunov exponents as function of time. The program can be simply adapted to plot the evolution of the Lyapunov exponents as function of orders, or a function of a bifurcation parameter. A special attention is payed to the periodicity of fractional-order systems and its influences. The case of the Lorenz system is considered.


title: Isolating Functions at the Hardware Limit with Virtines

abstract: An important class of applications, including programs that leverage third-party libraries, programs that use user-defined functions in databases, and serverless applications, benefit from isolating the execution of untrusted code at the granularity of individual functions or function invocations. However, existing isolation mechanisms were not designed for this use case; rather, they have been adapted to it. We introduce \textit{virtines}, a new abstraction designed specifically for function granularity isolation, and describe how we build virtines from the ground up by pushing hardware virtualization to its limits. Virtines give developers fine-grained control in deciding which functions should run in isolated environments, and which should not. The virtine abstraction is a general one, and we demonstrate a prototype that adds extensions to the C language. We present a detailed analysis of the overheads of running individual functions in isolated VMs, and guided by those findings, we present Wasp, an embeddable hypervisor that allows programmers to easily use virtines. We describe several representative scenarios that employ individual function isolation, and demonstrate that virtines can be applied in these scenarios with only a few lines of changes to existing codebases and with acceptable slowdowns.


title: Statistical properties of the localization measure of chaotic eigenstates and the spectral statistics in a mixed-type billiard

abstract: We study the quantum localization in the chaotic eigenstates of a billiard with mixed-type phase space, after separating the regular and chaotic eigenstates, in the regime of slightly distorted circle billiard where the classical transport time in the momentum space is still large enough, although the diffusion is not normal. In quantum systems with discrete energy spectrum the Heisenberg time $t_H =2\pi \hbar/\Delta E$, where $\Delta E$ is the mean level spacing, is an important time scale.The classical transport time scale $t_T$ in relation to the Heisenberg time scale $t_H$ (their ratio is the parameter $\alpha=t_H/t_T$) determines the degree of localization of the chaotic eigenstates, whose measure $A$ is based on the information entropy. We show that $A$ is linearly related to normalized inverse participation ratio. The localization of chaotic eigenstates is reflected also in the fractional power-law repulsion between the nearest energy levels in the sense that the probability density to find successive levels on a distance $S$ goes like $\propto S^\beta$ for small $S$, where $0\leq\beta\leq1$, and $\beta=1$ corresponds to completely extended states. We show that the level repulsion exponent $\beta$ is empirically a rational function of $\alpha$, and the mean $\langle A \rangle$ as a function of $\alpha$ is also well approximated by a rational function. In both cases there is some scattering of the empirical data around the mean curve, which is due to the fact that $A$ actually has a distribution, typically with quite complex structure, but in the limit $\alpha\rightarrow \infty$ well described by the beta distribution. Like in other systems, $\beta$ goes from $0$ to $1$ when $\alpha$ goes from $0$ to $\infty$. $\beta$ is a function of $\langle A \rangle$, similar to the quantum kicked rotator and the stadium billiard.


title: High-accuracy liquid-sample $\beta$-NMR setup at ISOLDE

abstract: Recently there has been an increased interest to apply the sensitive $\beta$-decay asymmetry detected nuclear magnetic resonance ($\beta$-NMR) technique to biological studies. A liquid-sample $\beta$-NMR setup was build at ISOLDE to allow such investigations and to use the resolution gain of liquid-state NMR in nuclear physics. As part of this setup a magnetic field locking system, a set of printed circuit board shimming coils, a sample exchange system, a set of compact $\beta$-detectors and a custom experimental vacuum chamber were developed. The main magnetic field was stabilized down to the ppm level by the locking system while allowing the direct determination of the absolute magnetic field. The homogeneity of the magnetic field was improved to $\leq$~5~ppm over the sample volume by the shimming coils. Time spent on changing samples was reduced by a factor of five by the liquid sample exchange system. During experiments it was possible to continuously observe the liquid sample thanks to the custom chamber and compact $\beta$-detectors. The absolute field determination allows for a novel way to reference $\beta$-NMR measurements, removing the need for time consuming reference measurements. The improved accuracy and resolution resulting from these innovations allows the study of the distribution of nuclear magnetization and (bio)chemicals using high-accuracy liquid-sample $\beta$-NMR.


title: Isoptics of log-aesthetic curves

abstract: The log-aesthetic curve has a significant factor in the field of aesthetic design to meet the high industrial requirements. It has much exceptional property and a large number of research papers are published, since its introduction. It can be parametrized by tangential angle, and it is autoevolute, meaning that its evolute is also a log-aesthetic curve. In this paper, we will examine the isoptics of the log-aesthetic curve and determine if it is autoisoptic as well. We will show that the log-aesthetic curve only coincides its isoptic in some specific cases. We will also provide a new formula to calculate the isoptic curve of any tangential angle parametrized curves.


title: LED based photoacoustic NO$_2$ sensor with sub-ppb detection limit

abstract: A high-sensitivity LED based photoacoustic NO$_2$ sensor is demonstrated. Sensitive photoacoustic gas sensors based on incoherent light sources are typically limited by background noise and drifts due to a strong signal generated by light absorbed at the photoacoustic cell walls. Here we reach a sub-ppb detection limit and excellent stability by using cantilever-enhanced photoacoustic detection and performing a two-channel relative measurement. A white-light LED is used as a light source and the spectrum is divided into two wavelength channels with a dichroic filter. The photoacoustic signals generated by the two wavelength channels are measured simultaneously and used to solve the NO$_2$ concentration. The background signal is highly correlated between the two channels and its variations are suppressed in the relative measurement. A noise level below 1 ppb is reached with an averaging time of 70 s. This is, to the best of our knowledge, the first time a sub-ppb detection limit is demonstrated with an LED based photoacoustic NO$_2$ sensor. As LEDs are available at wide selection of emission wavelengths, the results show great potential for development of cost-effective and sensitive detectors for a variety of other trace gasses as well.


title: What are higher-order networks?

abstract: Network-based modeling of complex systems and data using the language of graphs has become an essential topic across a range of different disciplines. Arguably, this graph-based perspective derives its success from the relative simplicity of graphs: A graph consists of nothing more than a set of vertices and a set of edges, describing relationships between pairs of such vertices. This simple combinatorial structure makes graphs interpretable and flexible modeling tools. The simplicity of graphs as system models, however, has been scrutinized in the literature recently. Specifically, it has been argued from a variety of different angles that there is a need for higher-order networks, which go beyond the paradigm of modeling pairwise relationships, as encapsulated by graphs. In this survey article we take stock of these recent developments. Our goals are to clarify (i) what higher-order networks are, (ii) why these are interesting objects of study, and (iii) how they can be used in applications.

categories: cs.SI math.CO math.DS math.ST nlin.AO stat.TH


title: Higher Convexity and Iterated Second Moment Estimates

abstract: We prove bounds for the number of solutions to $$a_1 + \dots + a_k = a_1' + \dots + a_k'$$ over $N$-element sets of reals, which are sufficiently convex or near-convex. A near-convex set will be the image of a set with small additive doubling under a convex function with sufficiently many strictly monotone derivatives. We show, roughly, that every time the number of terms in the equation is doubled, an additional saving of $1$ in the exponent of the trivial bound $N^{2k-1}$ is made, starting from the trivial case $k=1$. In the context of near-convex sets we also provide explicit dependencies on the additive doubling parameters. Higher convexity is necessary for such bounds to hold, as evinced by sets of perfect powers of consecutive integers. We exploit these stronger assumptions using an idea of Garaev, rather than the ubiquitous Szemer\'edi-Trotter theorem, which has not been adapted in earlier results to embrace higher convexity. As an application we prove small improvements for the best known bounds for sumsets of convex sets under additional convexity assumptions.


title: AFDM: A Full Diversity Next Generation Waveform for High Mobility Communications

abstract: We present Affine Frequency Division Multiplexing (AFDM), a new chirp-based multicarrier transceiver scheme for high mobility communications in next-generation wireless systems. AFDM is based on discrete affine Fourier transform (DAFT), a generalization of discrete Fourier transform characterized with two parameters that can be adapted to better cope with doubly dispersive channels. Based on the derived input-output relation, the DAFT parameters underlying AFDM are set in such a way to avoid that time domain channel paths with distinct delays or Doppler frequency shifts overlap in the DAFT domain. The resulting DAFT domain impulse response thus conveys a full delay-Doppler representation of the channel. We show that AFDM can achieve the full diversity of linear time-varying (LTV) channels. Our analytical results are validated through numerical simulations, which evince that AFDM outperforms state-of-the-art multicarrier schemes in terms of bit error rate (BER) in doubly dispersive channels.


title: Backup Control Barrier Functions: Formulation and Comparative Study

abstract: The backup control barrier function (CBF) was recently proposed as a tractable formulation that guarantees the feasibility of the CBF quadratic programming (QP) via an implicitly defined control invariant set. The control invariant set is based on a fixed backup policy and evaluated online by forward integrating the dynamics under the backup policy. This paper is intended as a tutorial of the backup CBF approach and a comparative study to some benchmarks. First, the backup CBF approach is presented step by step with the underlying math explained in detail. Second, we prove that the backup CBF always has a relative degree 1 under mild assumptions. Third, the backup CBF approach is compared with benchmarks such as Hamilton Jacobi PDE and Sum-of-Squares on the computation of control invariant sets, which shows that one can obtain a control invariant set close to the maximum control invariant set under a good backup policy for many practical problems.


title: The Quantum Double Slit Experiment With Local Elements of Reality

abstract: We present a fully local treatment of the double slit experiment in the formalism of quantum field theory. Our exposition is predominantly pedagogical in nature and exemplifies the fact that there is an entirely local description of the quantum double slit interference that does not suffer from any supposed paradoxes usually related to the wave-particle duality. The wave-particle duality indeed vanishes in favour of the field picture in which particles should not be regarded as the primary elements of reality and only represent excitations of some specific field configurations. Our treatment is general and can be applied to any other phenomenon involving quantum interference of any bosonic or fermionic field, both spatially and temporally. For completeness, we present the full treatment of single qubit interference in the same spirit.


title: Manipulating topology in tailored artificial graphene nanoribbons

abstract: Topological phases of matter give rise to exotic physics that can be leveraged for next generation quantum computation and spintronic devices. Thus, the search for topological phases and the quantum states that they exhibit have become the subject of a massive research effort in condensed matter physics. Topologically protected states have been produced in a variety of systems, including artificial lattices, graphene nanoribbons (GNRs) and bismuth bilayers. Despite these advances, the real-time manipulation of individual topological states and their relative coupling, a necessary feature for the realization of topological qubits, remains elusive. Guided by first-principles calculations, we spatially manipulate robust, zero-dimensional topological states by altering the topological invariants of quasi-one-dimensional artificial graphene nanostructures. This is achieved by positioning carbon monoxide molecules on a copper surface to confine its surface state electrons into artificial atoms positioned to emulate the low-energy electronic structure of graphene derivatives. Ultimately, we demonstrate control over the coupling between adjacent topological states that are finely engineered and simulate complex Hamiltonians. Our atomic synthesis gives access to an infinite range of nanoribbon geometries, including those beyond the current reach of synthetic chemistry, and thus provides an ideal platform for the design and study of novel topological and quantum states of matter.


title: Z-Sequence: Photometric redshift predictions for galaxy clusters with sequential random k-nearest neighbours

abstract: We introduce Z-Sequence, a novel empirical model that utilises photometric measurements of observed galaxies within a specified search radius to estimate the photometric redshift of galaxy clusters. Z-Sequence itself is composed of a machine learning ensemble based on the k-nearest neighbours algorithm. We implement an automated feature selection strategy that iteratively determines appropriate combinations of filters and colours to minimise photometric redshift prediction error. We intend for Z-Sequence to be a standalone technique but it can be combined with cluster finders that do not intrinsically predict redshift, such as our own DEEP-CEE. In this proof-of-concept study we train, fine-tune and test Z-Sequence on publicly available cluster catalogues derived from the Sloan Digital Sky Survey. We determine the photometric redshift prediction error of Z-Sequence via the median value of $|\Delta z|/(1+z)$ (across a photometric redshift range of $0.05 \le \textit{z} \le 0.6$) to be $\sim0.01$ when applying a small search radius. The photometric redshift prediction error for test samples increases by 30-50 per cent when the search radius is enlarged, likely due to line-of-sight interloping galaxies. Eventually, we aim to apply Z-Sequence to upcoming imaging surveys such as the Legacy Survey of Space and Time to provide photometric redshift estimates for large samples of as yet undiscovered and distant clusters.


title: The derivative expansion in asymptotically safe quantum gravity: general setup and quartic order

abstract: We present a general framework to systematically study the derivative expansion of asymptotically safe quantum gravity. It is based on an exact decoupling and cancellation of different modes in the Landau limit, and implements a correct mode count as well as a regularisation based on geometrical considerations. It is applicable independent of the truncation order. To illustrate the power of the framework, we discuss the quartic order of the derivative expansion and its fixed point structure as well as physical implications.


title: Efficient discretization and preconditioning of the singularly perturbed Reaction Diffusion problem

abstract: We consider the reaction diffusion problem and present efficient ways to discretize and precondition in the singular perturbed case when the reaction term dominates the equation. Using the concepts of optimal test norm and saddle point reformulation, we provide efficient discretization processes for uniform and non-uniform meshes. We present a preconditioning strategy that works for a large range of the perturbation parameter. Numerical examples to illustrate the efficiency of the method are included for a problem on the unit square.


title: Spin-orbit resonances of high-eccentricity asteroids: regular, switching, and jumping

abstract: Few solar system asteroids and comets are found in high eccentricity orbits ($e > 0.9$) but in the primordial planetesimal disks and in exoplanet systems around dying stars such objects are believed to be common. For 2006 HY51, the main belt asteroid with the highest known eccentricity 0.9684, we investigate the probable rotational states today using our computer-efficient chaotic process simulation method. Starting with random initial conditions, we find that this asteroid is inevitably captured into stable spin-orbit resonances typically within tens to a hundred Myr. The resonances are confirmed by direct integration of the equation of motion in the vicinity of end-points. Most resonances are located at high spin values above 960 times the mean motion (such as 964:1 or 4169:4), corresponding to rotation periods of a few days. We discover three types of resonance in the high-eccentricity regime: 1) regular circulation with weakly librating aphelion velocities and integer-number spin-orbit commensurabilities; 2) switching resonances of higher order with orientation alternating between aligned (0 or $\pi$) and sidewise ($\pi/2$) angles at aphelia and perihelia; 3) jumping resonances with aphelion spin alternating between two quantum states in the absence of spin-orbit commensurability. The islands of equilibrium are numerous at high spin rates but small in parameter space area, so that it takes millions of orbits of chaotic wandering to accidentally entrap in one of them. We discuss the implications of this discovery for the origins and destiny of high-eccentricity objects and the prospects of extending this analysis to the full 3D treatment.


title: Preparation in Bohmian Mechanics

abstract: According to Bohmian mechanics, we see the particle, not the pilot wave. But to make predictions we need to know the wave. How do we learn about the wave to make predictions, if we only see the particle? I show that the puzzle can be solved, but only thanks to decoherence.


title: Design not Lost in Translation: A Case Study of an Intimate-Space Socially Assistive Robot for Emotion Regulation

abstract: We present a Research-through-Design case study of the design and development of an intimate-space tangible device perhaps best understood as a socially assistive robot, aimed at scaffolding children's efforts at emotional regulation. This case study covers the initial research device development, as well as knowledge transfer to a product development company towards translating the research into a workable commercial product that could also serve as a robust research product for field trials. Key contributions to the literature include: 1. sharing of lessons learned from the knowledge transfer process that can be useful to others interested in developing robust products, whether commercial or research, that preserve design values, while allowing for large scale deployment and research; 2. articulation of a design space in HCI/HRI--Human Robot Interaction--of intimate space socially assistive robots, with the current artifact as a central exemplar, contextualized alongside other related HRI artifacts.


title: Thermodynamic stability versus Kinetic Accessibility: Pareto Fronts for Programmable Self-Assembly

abstract: A challenge in designing self-assembling building blocks is to ensure the target state is both thermodynamically stable and kinetically accessible. These two objectives are known to be typically in competition, but it is not known how to simultaneously optimize them. We consider this problem through the lens of multi-objective optimization theory: we develop a genetic algorithm to compute the Pareto fronts characterizing the tradeoff between equilibrium probability and folding rate, for a model system of small polymers of colloids with tunable short-ranged interaction energies. We use a coarse-grained model for the particles' dynamics that allows us to efficiently search over parameters, for systems small enough to be enumerated. For most target states there is a tradeoff when the number of types of particles is small, with medium-weak bonds favouring fast folding, and strong bonds favouring high equilibrium probability. The tradeoff disappears when the number of particle types reaches a value $m_*$, that is usually much less than the total number of particles. This general approach of computing Pareto fronts allows one to identify the minimum number of design parameters to avoid a thermodynamic-kinetic tradeoff. However, we argue, by contrasting our coarse-grained model's predictions with those of Brownian dynamics simulations, that particles with short-ranged isotropic interactions should generically have a tradeoff, and avoiding it in larger systems will require orientation-dependent interactions.


title: A Machine Learning Framework for Real-time Inverse Modeling and Multi-objective Process Optimization of Composites for Active Manufacturing Control

abstract: For manufacturing of aerospace composites, several parts may be processed simultaneously using convective heating in an autoclave. Due to uncertainties including tool placement, convective Boundary Conditions (BCs) vary in each run. As a result, temperature histories in some of the parts may not conform to process specifications due to under-curing or over-heating. Thermochemical analysis using Finite Element (FE) simulations are typically conducted prior to fabrication based on assumed range of BCs. This, however, introduces unnecessary constraints on the design. To monitor the process, thermocouples (TCs) are placed under tools near critical locations. The TC data may be used to back-calculate BCs using trial-and-error FE analysis. However, since the inverse heat transfer problem is ill-posed, many solutions are obtained for given TC data. In this study, a novel machine learning (ML) framework is presented capable of optimizing air temperature cycle in real-time based on TC data from multiple parts, for active control of manufacturing. The framework consists of two recurrent Neural Networks (NN) for inverse modeling of the ill-posed curing problem at the speed of 300 simulations/second, and a classification NN for multi-objective optimization of the air temperature at the speed of 35,000 simulations/second. A virtual demonstration of the framework for process optimization of three composite parts with data from three TCs is presented.


title: Ultimate fate of apparent horizons during a binary black hole merger I: Locating and understanding axisymmetric marginally outer trapped surfaces

abstract: In classical numerical relativity, marginally outer trapped surfaces (MOTSs) are the main tool to locate and characterize black holes. For five decades it has been known that during a binary merger, a new outer horizon forms around the initial apparent horizons of the individual holes once they are sufficiently close together. However the ultimate fate of those initial horizons has remained a subject of speculation. Recent axisymmetric studies have shed new light on this process and this pair of papers essentially completes that line of research: we resolve the key features of the post-swallowing axisymmetric evolution of the initial horizons. This first paper introduces a new shooting-method for finding axisymmetric MOTSs along with a reinterpretation of the stability operator as the analogue of the Jacobi equation for families of MOTSs. Here, these tools are used to study exact solutions and initial data. In the sequel paper they are applied to black hole mergers.


title: Ultimate fate of apparent horizons during a binary black hole merger II: Horizons weaving back and forth in time

abstract: In this second part of a two-part paper, we discuss numerical simulations of a head-on merger of two non-spinning black holes. We resolve the fate of the original two apparent horizons by showing that after intersecting, their world tubes "turn around" and continue backwards in time. Using the method presented in the first paper to locate these surfaces, we resolve several such world tubes evolving and connecting through various bifurcations and annihilations. This also draws a consistent picture of the full merger in terms of apparent horizons, or more generally, marginally outer trapped surfaces (MOTSs). The MOTS stability operator provides a natural mechanism to identify MOTSs which should be thought of as black hole boundaries. These are the two initial ones and the final remnant. All other MOTSs lie in the interior and are neither stable nor inner trapped.


title: Characterizing Outflows in the Cygnus X Region

abstract: In this paper, we perform an analysis of 13 outflows in the Cygnus X star-forming region. We use the James Clerk Maxwell Telescope observations of $^{13}$CO(3-2) and C$^{18}$O(3-2) molecular emission lines combined with archival $^{12}$CO(3-2) data. Using these new observations, we measure the mechanical properties of the outflows, and identify the associated protostars, finding their properties consistent with previous surveys of outflows throughout the Milky Way. Finally, we develop and test a method to measure the same properties using the existing $^{12}$CO(3-2) line data alone, finding the properties agree to within a factor of 2.


title: Global well-posedness for $H^{-1}(\mathbb{R})$ perturbations of KdV with exotic spatial asymptotics

abstract: Given a suitable solution $V(t,x)$ to the Korteweg--de Vries equation on the real line, we prove global well-posedness for initial data $u(0,x) \in V(0,x) + H^{-1}(\mathbb{R})$. Our conditions on $V$ do include regularity but do not impose any assumptions on spatial asymptotics. We show that periodic profiles $V(0,x) \in H^5(\mathbb{R}/\mathbb{Z})$ satisfy our hypotheses. In particular, we can treat localized perturbations of the much-studied periodic traveling wave solutions (cnoidal waves) of KdV. In our companion paper we show that smooth step-like initial data also satisfy our hypotheses. We employ the method of commuting flows introduced by Killip and Vi\c{s}an; in the special case $V\equiv 0$, we recover their sharp $H^{-1}(\mathbb{R})$ result.


title: Restoring degraded speech via a modified diffusion model

abstract: There are many deterministic mathematical operations (e.g. compression, clipping, downsampling) that degrade speech quality considerably. In this paper we introduce a neural network architecture, based on a modification of the DiffWave model, that aims to restore the original speech signal. DiffWave, a recently published diffusion-based vocoder, has shown state-of-the-art synthesized speech quality and relatively shorter waveform generation times, with only a small set of parameters. We replace the mel-spectrum upsampler in DiffWave with a deep CNN upsampler, which is trained to alter the degraded speech mel-spectrum to match that of the original speech. The model is trained using the original speech waveform, but conditioned on the degraded speech mel-spectrum. Post-training, only the degraded mel-spectrum is used as input and the model generates an estimate of the original speech. Our model results in improved speech quality (original DiffWave model as baseline) on several different experiments. These include improving the quality of speech degraded by LPC-10 compression, AMR-NB compression, and signal clipping. Compared to the original DiffWave architecture, our scheme achieves better performance on several objective perceptual metrics and in subjective comparisons. Improvements over baseline are further amplified in a out-of-corpus evaluation setting.


title: Earnings-21: A Practical Benchmark for ASR in the Wild

abstract: Commonly used speech corpora inadequately challenge academic and commercial ASR systems. In particular, speech corpora lack metadata needed for detailed analysis and WER measurement. In response, we present Earnings-21, a 39-hour corpus of earnings calls containing entity-dense speech from nine different financial sectors. This corpus is intended to benchmark ASR systems in the wild with special attention towards named entity recognition. We benchmark four commercial ASR models, two internal models built with open-source tools, and an open-source LibriSpeech model and discuss their differences in performance on Earnings-21. Using our recently released fstalign tool, we provide a candid analysis of each model's recognition capabilities under different partitions. Our analysis finds that ASR accuracy for certain NER categories is poor, presenting a significant impediment to transcript comprehension and usage. Earnings-21 bridges academic and commercial ASR system evaluation and enables further research on entity modeling and WER on real world audio.


title: Scalable Predictive Time-Series Analysis of COVID-19: Cases and Fatalities

abstract: COVID 19 is an acute disease that started spreading throughout the world, beginning in December 2019. It has spread worldwide and has affected more than 7 million people, and 200 thousand people have died due to this infection as of Oct 2020. In this paper, we have forecasted the number of deaths and the confirmed cases in Los Angeles and New York of the United States using the traditional and Big Data platforms based on the Times Series: ARIMA and ETS. We also implemented a more sophisticated time-series forecast model using Facebook Prophet API. Furthermore, we developed the classification models: Logistic Regression and Random Forest regression to show that the Weather does not affect the number of the confirmed cases. The models are built and run in legacy systems (Azure ML Studio) and Big Data systems (Oracle Cloud and Databricks). Besides, we present the accuracy of the models.


title: Making squeezed-coherent states concrete by determining their wavefunction

abstract: With the successes of the Laser Interferometer Gravitational-wave Observatory, we anticipate increased interest in working with squeezed states in the undergraduate and graduate quantum-mechanics classroom. Because squeezed-coherent states are minimum uncertainty states, their wavefunctions in position and momentum space must be Gaussians. But this result is rarely discussed in treatments of squeezed states in quantum textbooks or quantum optics textbooks. In this work, we show three different ways to construct the wavefunction for squeezed-coherent states: (i) a differential equation-based approach; (ii) an approach that uses an expansion in terms of the simple-harmonic oscillator wavefunctions; and (iii) a fully operator-based approach. We do this to illustrate that the concept of the wavefunction can be introduced no matter what methodology an instructor wishes to use. We hope that working with the wavefunction will help demystify the concept of a squeezed-coherent state.


title: Towards Automated Acceptance testing for industrial robots

abstract: Industrial robots are important machines applied in numerous modern industries that execute repetitive tasks with high accuracy, replacing or supporting dangerous jobs. In this kind of system, with increased complexity in which cost is related to the time the system keeps working, the system must operate with a minimum number of failures. In other words, a quality aspect important in industry is reliability. We hypothesize that Automated Acceptance Testing improves reliability for industrial robot program. We present the research question, the motivation for this study, our hypothesis and future research efforts.


title: On the Analytic Invariants and Semiroots of Plane Branches

abstract: The value semigroup of a $k$-semiroot $C_k$ of a plane branch $C$ allow us to recover part of the value semigroup $\Gamma =\langle v_0,\ldots ,v_g\rangle$ of $C$, that is, it is related to topological invariants of $C$. In this paper we consider the set of values of differentials $\Lambda_k$ of $C_k$, that is an analytical invariant, and we show how it determine part of the set of values of differentials $\Lambda$ of $C$. As a consequence, in a fixed topological class, we relate the Tjurina number $\tau$ of $C$ with the Tjurina number of $C_k$. In particular, we show that $\tau\leq \mu-\frac{3n_g-2}{4}\mu_{g-1}$ where $n_g=gcd(v_0,\ldots ,v_{g-1})$, $\mu$ and $\mu_{g-1}$ denote the Milnor number of $C$ and $C_{g-1}$ respectively. If $n_g=2$, we have that $\tau=\mu-\mu_{g-1}$ for any curve in the topological class determined by $\Gamma$ that is a generalization of a result obtained by Luengo and Pfister.


title: Optimal Cost Design for Model Predictive Control

abstract: Many robotics domains use some form of nonconvex model predictive control (MPC) for planning, which sets a reduced time horizon, performs trajectory optimization, and replans at every step. The actual task typically requires a much longer horizon than is computationally tractable, and is specified via a cost function that cumulates over that full horizon. For instance, an autonomous car may have a cost function that makes a desired trade-off between efficiency, safety, and obeying traffic laws. In this work, we challenge the common assumption that the cost we optimize using MPC should be the same as the ground truth cost for the task (plus a terminal cost). MPC solvers can suffer from short planning horizons, local optima, incorrect dynamics models, and, importantly, fail to account for future replanning ability. Thus, we propose that in many tasks it could be beneficial to purposefully choose a different cost function for MPC to optimize: one that results in the MPC rollout having low ground truth cost, rather than the MPC planned trajectory. We formalize this as an optimal cost design problem, and propose a zeroth-order optimization-based approach that enables us to design optimal costs for an MPC planning robot in continuous MDPs. We test our approach in an autonomous driving domain where we find costs different from the ground truth that implicitly compensate for replanning, short horizon, incorrect dynamics models, and local minima issues. As an example, the learned cost incentivizes MPC to delay its decision until later, implicitly accounting for the fact that it will get more information in the future and be able to make a better decision. Code and videos available at https://sites.google.com/berkeley.edu/ocd-mpc/.


title: The Role of Dust, UV Luminosity and Large-scale Environment on the Escape of Lya Photons: A Case Study of a Protocluster field at z = 3.1

abstract: We present a detailed characterization of the Lya properties for 93 Lya emitters (LAEs) at z~3.1 selected from the D1 field of the Canada-France-Hawaii-Telescope Legacy Survey, including 24 members of a massive protocluster. The median-stacked Lya image shows an extended Lya halo (LAH) surrounding the galaxy with the exponential scale length 4.9+/-0.7kpc, which accounts for roughly half of the total line flux. Accounting for the LAH contribution, the total Lya escape fraction, f_esc, is 40+/-26%. Combining the dataset with existing measurements, we find a dependence of f_esc on the galaxy's UV slope (beta) and UV luminosity (L_UV). The simultaneous use of both parameters allows prediction of f_esc within 0.18dex, a substantial improvement over 0.23dex when only beta is used. The correlation between f_esc and E(B-V) suggests that Lya photons undergo interstellar dust attenuation in a similar manner to continuum photons. Yet, Lya transmission is typically higher than that expected for continuum photons at similar wavelength by a factor, which depends on UV luminosity, up to 2 in the samples we studied. These results hint at complex geometries and physical conditions of the interstellar medium, which affect the Lya transmission or production. Alternatively, the dust law may change with luminosity leading to over-or under-estimation of f_esc. Finally, we report that protocluster member LAEs tend to be bluer and more UV-luminous than their field cousins, resulting in systematically higher f_esc values. We speculate that it may be due to the widespread formation of young low-mass galaxies in dense gas-rich environments.


title: PROFIT: Projection-based Test in Longitudinal Functional Data

abstract: In many modern applications, a dependent functional response is observed for each subject over repeated time, leading to longitudinal functional data. In this paper, we propose a novel statistical procedure to test whether the mean function varies over time. Our approach relies on reducing the dimension of the response using data-driven orthogonal projections and it employs a likelihood-based hypothesis testing. We investigate the methodology theoretically and discuss a computationally efficient implementation. The proposed test maintains the type I error rate, and shows excellent power to detect departures from the null hypothesis in finite sample simulation studies. We apply our method to the longitudinal diffusion tensor imaging study of multiple sclerosis (MS) patients to formally assess whether the brain's health tissue, as summarized by fractional anisotropy (FA) profile, degrades over time during the study period.


title: Localization in Quantum Walks with a Single Lattice Defect: A Comparative Study

abstract: We study how a single lattice defect in a discrete time quantum walk affects the return probability of a quantum particle. This defect at the starting position is modeled by a quantum coin that is distinct from the others over the lattice. This coin has a dependence on $\omega$ which quantifies the intensity of the localization. For some sorts of lattice defects, we show how the localization can have a dependence just on $\omega$, and also, the polar $\alpha$ and azimuth $\beta$ angles of the initial qubit by numerical calculations. We propose a lattice defect whose localization has additional dependence on $\beta+\omega$, leading to extra localization profiles. We compare the quantum walks with our lattice defect to the earlier ones, and we discuss their spreading and survival probability.


title: Energetics and Kinetics Requirements for Organic Solar Cells to 2 Break the 20% Power Conversion Efficiency Barrier

abstract: The thermodynamic limit for the efficiency of solar cells is predominantly defined by the energy bandgap of the used semiconductor. In case of organic solar cells both energetics and kinetics of three different species play role: excitons, charge transfer states and charge separated states. In this work, we clarify the effect of the relative energetics and kinetics of these species on the recombination and generation dynamics. Making use of detailed balance, we develop an analytical framework describing how the intricate interplay between the different species influence the photocurrent generation, the recombination, and the open-circuit voltage in organic solar cells. Furthermore, we clarify the essential requirements for equilibrium between excitons, CT states and charge carriers to occur. Finally, we find that the photovoltaic parameters are not only determined by the relative energy level between the different states but also by the kinetic rate constants. These findings provide vital insights into the operation of state-of-art non-fullerene organic solar cells with low offsets.


title: Joint Mean-Vector and Var-Matrix estimation for Locally Stationary VAR(1) processes

abstract: During the last two decades, locally stationary processes have been widely studied in the time series literature. In this paper we consider the locally-stationary vector-auto-regression model of order one, or LS-VAR(1), and estimate its parameters by weighted least squares. The LS-VAR(1) we consider allows for a smoothly time-varying non-diagonal VAR matrix, as well as for a smoothly time-varying non-zero mean. The weighting scheme is based on kernel smoothers. The time-varying mean and the time-varying VAR matrix are estimated jointly, and the definition of the local-linear weighting matrix is provided in closed-from. The quality of the estimated curves is illustrated through simulation results.


title: Model Checking for Verification of Quantum Circuits

abstract: In this talk, we will describe a framework for assertion-based verification (ABV) of quantum circuits by applying model checking techniques for quantum systems developed in our previous work, in which: (i) Noiseless and noisy quantum circuits are modelled as operator- and super-operator-valued transition systems, respectively, both of which can be further represented by tensor networks. (ii) Quantum assertions are specified by a temporal extension of Birkhoff-von Neumann quantum logic. Their semantics is defined based on the design decision: they will be used in verification of quantum circuits by simulation on classical computers or human reasoning rather than by quantum physics experiments (e.g. testing through measurements); (iii) Algorithms for reachability analysis and model checking of quantum circuits are developed based on contraction of tensor networks. We observe that many optimisation techniques for computing relational products used in BDD-based model checking algorithms can be generalised for contracting tensor networks of quantum circuits.


title: Normalized multivariate time series causality analysis and causal graph reconstruction

abstract: Causality analysis is an important problem lying at the heart of science, and is of particular importance in data science and machine learning. An endeavor during the past 16 years viewing causality as real physical notion so as to formulate it from first principles, however, seems to go unnoticed. This study introduces to the community this line of work, with a long-due generalization of the information flow-based bivariate time series causal inference to multivariate series, based on the recent advance in theoretical development. The resulting formula is transparent, and can be implemented as a computationally very efficient algorithm for application. It can be normalized, and tested for statistical significance. Different from the previous work along this line where only information flows are estimated, here an algorithm is also implemented to quantify the influence of a unit to itself. While this forms a challenge in some causal inferences, here it comes naturally, and hence the identification of self-loops in a causal graph is fulfilled automatically as the causalities along edges are inferred. To demonstrate the power of the approach, presented here are two applications in extreme situations. The first is a network of multivariate processes buried in heavy noises (with the noise-to-signal ratio exceeding 100), and the second a network with nearly synchronized chaotic oscillators. In both graphs, confounding processes exist. While it seems to be a huge challenge to reconstruct from given series these causal graphs, an easy application of the algorithm immediately reveals the desideratum. Particularly, the confounding processes have been accurately differentiated. Considering the surge of interest in the community, this study is very timely.


title: Relative tilting theory in abelian categories I: Auslander-Buchweitz-Reiten approximations theory in subcategories and cotorsion pairs

abstract: In this paper we introduce a special kind of relative (co)resolutions associated to a pair of classes of objects in an abelian category $\mathcal{C}.$ We will see that, by studying these relative (co)resolutions, we get a possible generalization of a part of the Auslander-Buchweitz approximation theory that is useful for developing $n$-$\mathcal{X}$-tilting theory in [4]. With this goal, new concepts as $\mathcal{X}$-complete and $\mathcal{X}$-hereditary pairs are introduced as a generalization of complete and hereditary cotorsion pairs. These pairs appear in a natural way in the study of the category of representations of a quiver in an abelian category [5]. Our main results will include an existence theorem for relative approximations, among other results related with closure properties of relative (co)resolution classes and relative homological dimensions which are essential in the development of $n$-$\mathcal{X}$-tilting theory in [4].


title: Spontaneous Scalarization in Scalar-Tensor Theories with Conformal Symmetry as an Attractor

abstract: Motivated by constant-G theory, we introduce a one-parameter family of scalar-tensor theories as an extension of constant-G theory in which the conformal symmetry is a cosmological attractor. Since the model has the coupling function of negative curvature, we expect spontaneous scalarization occurs and the parameter is constrained by pulsar-timing measurements. Modeling neutron stars with realistic equation of states, we study the structure of neutron stars and calculate the effective scalar coupling with the neutron star in these theories. We find that within the parameter region where the observational constraints are satisfied, the effective scalar coupling almost coincides with that derived using the quadratic model with the same curvature. This indicates that the constraints obtained by the quadratic model will be used to limit the curvature of the coupling function universally in the future.


title: Nematicity Liquid in a Trimerized-Kagome Antiferromagnet

abstract: We theoretically study low-temperature properties of the antiferromagnetic Heisenberg model with half-integer spin $S$ on the Kagome lattice with large trimerization. We have derived a low-energy effective model for general $S$ in terms of spin and nematicity operators in triangular units, and studied their low-temperature correlations for S=3/2 by classical Monte Carlo simulations. The previous study for the S=1/2 case [M. Ferrero et al., Phys. Rev. B 68, 214431 (2003)] reported a spin liquid state at low temperatures driven by a glassy behavior of isolated dimers and trimers of nematicities. The results for S=3/2 show that nematicity dimers and trimers are connected by weak links to form a defective planar network. At very low temperatures, nematicities show glassy slow dynamics, and cluster spins continue to fluctuate in a nonperiodically frozen background of nematicities. The characteristic time of glassy dynamics scales with temperature following a power law.


title: A field guide to cultivating computational biology

abstract: Biomedical research centers can empower basic discovery and novel therapeutic strategies by leveraging their large-scale datasets from experiments and patients. This data, together with new technologies to create and analyze it, has ushered in an era of data-driven discovery which requires moving beyond the traditional individual, single-discipline investigator research model. This interdisciplinary niche is where computational biology thrives. It has matured over the past three decades and made major contributions to scientific knowledge and human health, yet researchers in the field often languish in career advancement, publication, and grant review. We propose solutions for individual scientists, institutions, journal publishers, funding agencies, and educators.


title: Can You Trust Your Trust Measure?

abstract: Trust in human-robot interactions (HRI) is measured in two main ways: through subjective questionnaires and through behavioral tasks. To optimize measurements of trust through questionnaires, the field of HRI faces two challenges: the development of standardized measures that apply to a variety of robots with different capabilities, and the exploration of social and relational dimensions of trust in robots (e.g., benevolence). In this paper we look at how different trust questionnaires fare given these challenges that pull in different directions (being general vs. being exploratory) by studying whether people think the items in these questionnaires are applicable to different kinds of robots and interactions. In Study 1 we show that after being presented with a robot (non-humanoid) and an interaction scenario (fire evacuation), participants rated multiple questionnaire items such as "This robot is principled" as "Non-applicable to robots in general" or "Non-applicable to this robot". In Study 2 we show that the frequency of these ratings change (indeed, even for items rated as N/A to robots in general) when a new scenario is presented (game playing with a humanoid robot). Finally, while overall trust scores remained robust to N/A ratings, our results revealed potential fallacies in the way these scores are commonly interpreted. We conclude with recommendations for the development, use and results-reporting of trust questionnaires for future studies, as well as theoretical implications for the field of HRI.


title: The Abundancy Index and Feebly Amicable Numbers

abstract: This research explores the sum of divisors - $\sigma(n)$ - and the abundancy index given by the function $\frac{\sigma(n)}{n}$. We give a generalization of amicable pairs - feebly amicable pairs (also known as harmonious pairs), that is $m,n$ such that $\frac{n}{\sigma(n)}+ \frac{m}{\sigma(m)}=1$. We first give some groundwork in introductory number theory, then the goal of the paper is to determine if all numbers are feebly amicable with at least one other number by using known results about the abundancy index. We establish that not all numbers are feebly amicable with at least one other number. We generate data using the R programming language and give some questions and conjectures.


title: Restriction of exponential sums to hypersurfaces

abstract: We prove moment inequalities for exponential sums with respect to singular measures, whose Fourier decay matches those of curved hypersurfaces. Our emphasis will be on proving estimates that are sharp with respect to the scale parameter $N$, apart from $N^\epsilon$ losses. In a few instances, we manage to remove these losses.


title: Ghost factors in Gauss-sum factorization with transmon qubits

abstract: A challenge in the Gauss sums factorization scheme is the presence of ghost factors - non-factors that behave similarly to actual factors of an integer - which might lead to the misidentification of non-factors as factors or vice versa, especially in the presence of noise. We investigate Type II ghost factors, which are the class of ghost factors that cannot be suppressed with techniques previously laid out in the literature. The presence of Type II ghost factors and the coherence time of the qubit set an upper limit for the total experiment time, and hence the largest factorizable number with this scheme. Discernability is a figure of merit introduced to characterize this behavior. We introduce preprocessing as a strategy to increase the discernability of a system, and demonstrate the technique with a transmon qubit. This can bring the total experiment time of the system closer to its decoherence limit, and increase the largest factorizable number.


title: Collisional flavor instabilities of supernova neutrinos

abstract: A lingering mystery in core-collapse supernova theory is how collective neutrino oscillations affect the dynamics. All previously identified flavor instabilities, some of which might make the effects considerable, are essentially collisionless phenomena. Here it is shown that collisional instabilities exist as well. They are associated with asymmetries between the neutrino and antineutrino interaction rates, are possibly prevalent deep inside supernovae, and pose an unusual instance of decoherent interactions with a thermal environment causing the sustained growth of quantum coherence.


title: Analysis and Modeling of Driver Behavior with Integrated Feedback of Visual and Haptic Information Under Shared Control

abstract: The thesis presents contributions made to the evaluation and design of a haptic guidance system on improving driving performance in cases of normal and degraded visual information, which are based on behavior experiments, modeling and numerical simulations. The effect of shared control on driver behavior in cases of normal and degraded visual information has been successfully evaluated experimentally and numerically. The evaluation results indicate that the proposed haptic guidance system is capable of providing reliable haptic information, and is effective on improving lane following performance in the conditions of visual occlusion from road ahead and declined visual attention under fatigue driving. Moreover, the appropriate degree of haptic guidance is highly related to the reliability of visual information perceived by the driver, which suggests that designing the haptic guidance system based on the reliability of visual information would allow for greater driver acceptance. Furthermore, the parameterized driver model, which considers the integrated feedback of visual and haptic information, is capable of predicting driver behavior under shared control, and has the potential of being used for designing and evaluating the haptic guidance system.


title: Nonparametric estimation of marginal distributions for unordered pairs

abstract: In this article, we consider the estimation of the marginal distributions for pairs of data are recorded, with unobserved order in each pair. New estimators are proposed and their asymptotic properties are established, by proving a Glivenko-Cantelli theorem and a functional central limit result. Results from a simulation study are included and we illustrate the applicability of the method on the homologous chromosomes data.


title: Grasp Synthesis for Novel Objects Using Heuristic-based and Data-driven Active Vision Methods

abstract: In this work, we present several heuristic-based and data-driven active vision strategies for viewpoint optimization of an arm-mounted depth camera for the purpose of aiding robotic grasping. These strategies aim to efficiently collect data to boost the performance of an underlying grasp synthesis algorithm. We created an open-source benchmarking platform in simulation (https://github.com/galenbr/2021ActiveVision), and provide an extensive study for assessing the performance of the proposed methods as well as comparing them against various baseline strategies. We also provide an experimental study with a real-world setup by utilizing an existing grasping planning benchmark in the literature. With these analyses, we were able to quantitatively demonstrate the versatility of heuristic methods that prioritize certain types of exploration, and qualitatively show their robustness to both novel objects and the transition from simulation to the real world. We identified scenarios in which our methods did not perform well and scenarios that are objectively difficult, and present a discussion on which avenues for future research show promise.


title: Solids in the space of the Veronese surface in even characteristic

abstract: We classify the orbits of solids in the projective space $\text{PG}(5,q)$, $q$ even, under the setwise stabiliser $K \cong \text{PGL}(3,q)$ of the Veronese surface. For each orbit, we provide an explicit representative $S$ and determine two combinatorial invariants: the point-orbit distribution and the hyperplane-orbit distribution. These invariants characterise the orbits except in two specific cases (in which the orbits are distinguished by their line-orbit distributions). In addition, we determine the stabiliser of $S$ in $K$, thereby obtaining the size of each orbit. As a consequence, we obtain a proof of the classification of pencils of conics in $\text{PG}(2,q)$, $q$ even, which to the best of our knowledge has been heretofore missing in the literature.


title: Weyl Semimetal Made Ideal with a Crystal of Raman Light and Atoms

abstract: Optical lattices are known for their flexibility to emulate condensed matter physics and beyond. Based on an early theoretical proposal [Science Bulletin 65, 2080 (2020)], a recent experiment published by Wang et al. [Science 372, 271 (2021)] accomplishes the first experimental realization of topological band structure of the ideal Weyl semimetal in ultracold atomic matter, prompting fundamental interest in the context of gapless topological physics. With a neat design of 3D spin-orbit interaction, the experiment has probed the gapless band topology through spin texture imaging and quantum quench dynamics. This work has far reaching implications to topological effects and quantum anomaly in condensed matter and high energy physics.


title: Decentralized Federated Averaging

abstract: Federated averaging (FedAvg) is a communication efficient algorithm for the distributed training with an enormous number of clients. In FedAvg, clients keep their data locally for privacy protection; a central parameter server is used to communicate between clients. This central server distributes the parameters to each client and collects the updated parameters from clients. FedAvg is mostly studied in centralized fashions, which requires massive communication between server and clients in each communication. Moreover, attacking the central server can break the whole system's privacy. In this paper, we study the decentralized FedAvg with momentum (DFedAvgM), which is implemented on clients that are connected by an undirected graph. In DFedAvgM, all clients perform stochastic gradient descent with momentum and communicate with their neighbors only. To further reduce the communication cost, we also consider the quantized DFedAvgM. We prove convergence of the (quantized) DFedAvgM under trivial assumptions; the convergence rate can be improved when the loss function satisfies the P{\L} property. Finally, we numerically verify the efficacy of DFedAvgM.


title: Charged black-bounce spacetimes

abstract: Given the recent development of rotating black-bounce-Kerr spacetimes, for both theoretical and observational purposes it becomes interesting to see whether it might be possible to construct black-bounce variants of the entire Kerr-Newman family. Specifically, herein we shall consider black-bounce-Reissner-Nordstr\"om and black-bounce-Kerr-Newman spacetimes as particularly simple and clean everywhere-regular black hole "mimickers" that deviate from the Kerr-Newman family in a precisely controlled and minimal manner, and smoothly interpolate between regular black holes and traversable wormholes. While observationally the electric charges on astrophysical black holes are likely to be extremely low, $|Q|/m \ll 1$, introducing any non-zero electric charge has a significant theoretical impact. In particular, we verify the existence of a Killing tensor (and associated Carter-like constant) but without the full Killing tower of principal tensor and Killing-Yano tensor, also we discuss how, assuming general relativity, the black-bounce-Kerr-Newman solution requires an interesting, non-trivial matter/energy content.


title: Rate compatible reconciliation for continuous-variable quantum key distribution using Raptor-like LDPC codes

abstract: In the practical continuous-variable quantum key distribution (CV-QKD) system, the postprocessing process, particularly the error correction part, significantly impacts the system performance. Multi-edge type low-density parity-check (MET-LDPC) codes are suitable for CV-QKD systems because of their Shannon-limit-approaching performance at a low signal-to-noise ratio (SNR). However, the process of designing a low-rate MET-LDPC code with good performance is extremely complicated. Thus, we introduce Raptor-like LDPC (RL-LDPC) codes into the CV-QKD system, exhibiting both the rate compatible property of the Raptor code and capacity-approaching performance of MET-LDPC codes. Moreover, this technique can significantly reduce the cost of constructing a new matrix. We design the RL-LDPC matrix with a code rate of 0.02 and easily and effectively adjust this rate from 0.016 to 0.034. Simulation results show that we can achieve more than 98% reconciliation efficiency in a range of code rate variation using only one RL-LDPC code that can support high-speed decoding with an SNR less than -16.45 dB. This code allows the system to maintain a high key extraction rate under various SNRs, paving the way for practical applications of CV-QKD systems with different transmission distances.


title: Quantum Discord for multiqubit Systems

abstract: We evaluate analytically the quantum discord for a large family of multiqubit states. It is interesting to note that the quantum discord of three-qubits and five-qubits is the same, as is the quantum discord of two-qubits and six-qubits. We discover that the quantum discord of this family states can be concluded into three categories. The level surfaces of the quantum discord in the three categories is shown through images. Furthermore, we investigated the dynamic behavior of quantum discord under decoherence. For the odd partite systems, we prove the frozen phenomenon of quantum discord doesn't exist under the phase flip channel, while it can be found in the even partite systems.


title: Eigenbackground Revisited: Can We Model the Background with Eigenvectors?

abstract: Using dominant eigenvectors for background modeling (usually known as Eigenbackground) is a common technique in the literature. However, its results suffer from noticeable artifacts. Thus have been many attempts to reduce the artifacts by making some improvements/enhancement in the Eigenbackground algorithm. In this paper, we show the main problem of the Eigenbackground is in its own core and in fact, it is not a good idea to use strongest eigenvectors for modeling the background. Instead, we propose an alternative solution by exploiting the weakest eigenvectors (which are usually thrown away and treated as garbage data) for background modeling. MATLAB codes are available at \url{https://github.com/mamintoosi/Eigenbackground-Revisited}


title: Correspondence between momentum dependent relaxation time and field redefinition of relativistic hydrodynamic theory

abstract: In this article a correspondence has been established between the out of equilibrium system dissipation and the thermodynamic field redefinition of the macroscopic variables through the momentum dependent relaxation time approximation (MDRTA) solution of relativistic transport equation. Here, it has been shown that the out of equilibrium thermodynamic fields are not uniquely defined and are subjected to include dissipative effects from the medium. A second order relativistic hydrodynamic theory has been developed including such dissipative effects. The necessary conditions for developing a hydrodynamic theory has been fulfilled, (i) the thermodynamic identities incorporating such redefined fields have been shown to conserve the energy-momentum tensor perfectly under MDRTA, (ii) the non-negativity of entropy production remains unaffected by the inclusion of such dissipative contributions in hydro fields as long as the independent transport coefficients remain positive.


title: Optical trap for an atom around the midpoint between two coupled identical parallel optical nanofibers

abstract: We study the trapping of a ground-state cesium atom in a small region around the midpoint between two coupled identical parallel optical nanofibers. We suggest to use a blue-detuned guided light field in the odd $\mathcal{E}_z$-sine array mode to produce an optical potential with a local minimum of exact zero at the midpoint between the nanofibers. We find that the effects of the van der Waals potential on the total trapping potential around the minimum point are not significant when the fiber separation distance and the power of the guided light field are large. For appropriate realistic parameters, a net trapping potential with a significant depth of about 1 mK, a large coherence time of several seconds, and a large recoil-heating-limited trap lifetime of several hours can be obtained. We investigate the dependencies of the trapping potential on the power of the guided light field, the fiber radius, the wavelength of light, and the fiber separation distance.


title: Fluctuation effects at the onset of $\bf 2k_F$ density wave order with two pairs of hot spots in two-dimensional metals

abstract: We analyze quantum fluctuation effects at the onset of charge or spin density wave order in two-dimensional metals with an incommensurate nesting ($2k_F$) wave vector connecting two pairs of hot spots on the Fermi surface. We first compute the momentum and frequency dependence of the fermion self-energy near the hot spots to leading order in a perturbation expansion (one loop). Non-Fermi liquid behavior with a linear (in energy) quasi-particle decay rate and a logarithmically vanishing quasi-particle weight is obtained. The momentum dependence of the self-energy entails only finite renormalizations of the Fermi velocity and the Fermi surface curvature at the hot spots. The perturbative one-loop result is not self-consistent and casts doubt on the stability of the $2k_F$ quantum critical point. We construct a self-consistent solution of the one-loop equations with self-energy feedback, where the quantum critical point is stabilized rather than being destroyed by fluctuations, while the non-Fermi liquid behavior as found in the perturbative one-loop calculation is confirmed.


title: $\Gamma$-graphic delta-matroids and their applications

abstract: For an abelian group $\Gamma$, a $\Gamma$-labelled graph is a graph whose vertices are labelled by elements of $\Gamma$. We prove that a certain collection of edge sets of a $\Gamma$-labelled graph forms a delta-matroid, which we call a $\Gamma$-graphic delta-matroid, and provide a polynomial-time algorithm to solve the separation problem, which allows us to apply the symmetric greedy algorithm of Bouchet to find a maximum weight feasible set in such a delta-matroid. We present two algorithmic applications on graphs; Maximum Weight Packing of Trees of Order Not Divisible by $k$ and Maximum Weight $S$-Tree Packing. We also discuss various properties of $\Gamma$-graphic delta-matroids.


title: APRF-Net: Attentive Pseudo-Relevance Feedback Network for Query Categorization

abstract: Query categorization is an essential part of query intent understanding in e-commerce search. A common query categorization task is to select the relevant fine-grained product categories in a product taxonomy. For frequent queries, rich customer behavior (e.g., click-through data) can be used to infer the relevant product categories. However, for more rare queries, which cover a large volume of search traffic, relying solely on customer behavior may not suffice due to the lack of this signal. To improve categorization of rare queries, we adapt the Pseudo-Relevance Feedback (PRF) approach to utilize the latent knowledge embedded in semantically or lexically similar product documents to enrich the representation of the more rare queries. To this end, we propose a novel deep neural model named Attentive Pseudo Relevance Feedback Network (APRF-Net) to enhance the representation of rare queries for query categorization. To demonstrate the effectiveness of our approach, we collect search queries from a large commercial search engine, and compare APRF-Net to state-of-the-art deep learning models for text classification. Our results show that the APRF-Net significantly improves query categorization by 5.9% on F1@1 score over the baselines, which increases to 8.2% improvement for the rare (tail) queries. The findings of this paper can be leveraged for further improvements in search query representation and understanding.


title: In-Situ Assessment of Device-Side Compute Work for Dynamic Load Balancing in a GPU-Accelerated PIC Code

abstract: Maintaining computational load balance is important to the performant behavior of codes which operate under a distributed computing model. This is especially true for GPU architectures, which can suffer from memory oversubscription if improperly load balanced. We present enhancements to traditional load balancing approaches and explicitly target GPU architectures, exploring the resulting performance. A key component of our enhancements is the introduction of several GPU-amenable strategies for assessing compute work. These strategies are implemented and benchmarked to find the most optimal data collection methodology for in-situ assessment of GPU compute work. For the fully kinetic particle-in-cell code WarpX, which supports MPI+CUDA parallelism, we investigate the performance of the improved dynamic load balancing via a strong scaling-based performance model and show that, for a laser-ion acceleration test problem run with up to 6144 GPUs on Summit, the enhanced dynamic load balancing achieves from 62%--74% (88% when running on 6 GPUs) of the theoretically predicted maximum speedup; for the 96-GPU case, we find that dynamic load balancing improves performance relative to baselines without load balancing (3.8x speedup) and with static load balancing (1.2x speedup). Our results provide important insights into dynamic load balancing and performance assessment, and are particularly relevant in the context of distributed memory applications ran on GPUs.

categories: cs.DC physics.acc-ph physics.comp-ph physics.plasm-ph


title: Recording Reusable and Guided Analytics From Interaction Histories

abstract: The use of visual analytics tools has gained popularity in various domains, helping users discover meaningful information from complex and large data sets. Users often face difficulty in disseminating the knowledge discovered without clear recall of their exploration paths and analysis processes. We introduce a visual analysis tool that allows analysts to record reusable and guided analytics from their interaction logs. To capture the analysis process, we use a decision tree whose node embeds visualizations and guide to define a visual analysis task. The tool enables analysts to formalize analysis strategies, build best practices, and guide novices through systematic workflows.


title: Rotation of optically bound particle assembly due to scattering induced spin-orbit coupling of light

abstract: The optical binding of many particles has great potential to achieve the wide-area formation of a "crystal" of small materials. Unlike conventional optical binding, where the whole assembly of targeted particles is irradiated with light, if one can indirectly manipulate remote particles using a single trapped particle through optical binding, the degrees of freedom to create ordered structures will be greatly enhanced. In this Letter, we theoretically investigate the dynamics of the assembly of gold nanoparticles that is manipulated using a single particle trapped by a focused laser. As a result, we demonstrate that the spin--orbit coupling and angular momentum generation of light via scattering induce the assembly and rotational motion of particles through indirect optical force. This result opens the possibility of creating ordered structures with a wide area and manipulating them, controlling local properties using scanning laser beams.


title: Streaming Readout of the CLAS12 Forward Tagger Using TriDAS and JANA2

abstract: An effort is underway to develop streaming readout data acquisition system for the CLAS12 detector in Jefferson Lab's experimental Hall-B. Successful beam tests were performed in the spring and summer of 2020 using a 10GeV electron beam from Jefferson Lab's CEBAF accelerator. The prototype system combined elements of the TriDAS and CODA data acquisition systems with the JANA2 analysis/reconstruction framework. This successfully merged components that included an FPGA stream source, a distributed hit processing system, and software plugins that allowed offline analysis written in C++ to be used for online event filtering. Details of the system design and performance are presented.


title: Spectral-temporal features of repeating ( one-off ) FRBs and Axion Star

abstract: The fast radio bursts ( FRBs ) are energetic radio bursts with millisecond duration only observed at radio frequencies. The generation mechanism is still mysterious. We have proposed a generation mechanism of both repeating and one-off FRBs. They arise from the axion star collision with neutron star or magnetized accretion disk of galactic black hole. Once we accept the existence of the axions, we find that the mechanism well explain previously observed spectral-temporal features. In this paper we show that it also explains recently observed phenomena such as downward drifting in the repeating FRBs, etc.. Analysis of the downward drifting based on Doppler effects has been presented in recent papers, in which a superradiance system of molecular or atom has been proposed as a source of FRBs. We apply the analysis to our mechanism and find that it well explains the relation between the downward drifting rate and the duration of the repeating FRBs. The Doppler effects lead to the fact that the duration of radio burst with higher center frequency is shorter than that of radio burst with lower center frequency in the repeating FRBs. Our generation mechanism naturally explain polarization angle swing observed in the repeating FRB180301 and one-off FRBs. We also discuss the association between the FRB200428 and magnetar SGR J1935+2154. The X ray burst observed just after the observation of the FRB could be triggered by the axion star collision with the magnetar. We also explain the consistency of our generation mechanism with observed spectral-temporal differences in the repeating and one-off FRBs, e.g. longer duration ( smaller flux density ) of repeating FRBs than duration ( flux density ) of one-off FRBs.


title: Transfer training from smaller language model

abstract: Large language models have led to state-of-the-art accuracies across a range of tasks. However,training large language model needs massive computing resource, as more and more open source pre-training models are available, it is worthy to study how to take full advantage of available model. We find a method to save training time and resource cost by changing the small well-trained model to large model. We initialize a larger target model from a smaller source model by copy weight values from source model and padding with zeros or small initialization values on it to make the source and target model have approximate outputs, which is valid due to block matrix multiplication and residual connection in transformer structure. We test the target model on several data sets and find it is still comparable with the source model. When we continue training the target model, the training loss can start from a smaller value.


title: X-ray polarimetry of the Crab nebula with PolarLight: polarization recovery after the glitch and a secular position angle variation

abstract: We report follow-up observations of the Crab nebula with the PolarLight X-ray polarimeter, which revealed a possible variation in polarization associated with a pulsar glitch in 2019. The new observations confirm that the polarization has recovered roughly 100 days after the glitch. With the new observations, we find that the polarization angle (PA) measured with PolarLight from the total nebular emission has a difference of 18.0 +- 4.6 (deg) from that measured 42 years ago with OSO-8, indicating a secular evolution of polarization with either the Crab nebula or pulsar. The long-term variation in PA could be a result of multiple glitches in the history, magnetic reconnection or movement of synchrotron emitting structures in the nebula, or secular evolution of the pulsar magnetic geometry.


title: Convexification-based globally convergent numerical method for a 1D coefficient inverse problem with experimental data

abstract: To compute the spatially distributed dielectric constant from the backscattering data, we study a coefficient inverse problem for a 1D hyperbolic equation. To solve the inverse problem, we establish a new version of Carleman estimate and then employ this estimate to construct a cost functional which is strictly convex on a convex bounded set with an arbitrary diameter in a Hilbert space. The strict convexity property is rigorously proved. This result is called the convexification theorem and is considered as the central analytical result of this paper. Minimizing this convex functional by the gradient descent method, we obtain the desired numerical solution to the coefficient inverse problems. We prove that the gradient descent method generates a sequence converging to the minimizer and we also establish a theorem confirming that the minimizer converges to the true solution as the noise in the measured data and the regularization parameter tend to zero. Unlike the methods that are based on optimization, our convexification method converges globally in the sense that it delivers a good approximation of the exact solution without requiring any initial guess. Results of numerical studies of both computationally simulated and experimental data are presented.


title: Age of information distribution under dynamic service preemption

abstract: Age of Information (AoI) has emerged as an important quality-of-service measure for applications that prioritize delivery of the freshest information, e.g., virtual or augmented reality over mobile devices and wireless sensor networks used in the control of cyber-physical systems. We derive the Laplace transform of the stationary AoI for the M/GI/1/2 system with a "dynamic" service preemption and pushout policy depending on the existing service time of the in-service message. Thus, our system generalizes both the static M/GI/1/2 queue-pushout system without service preemption and the M/GI/1/1 bufferless system with service preemption - two systems considered to provide very good AoI performance. Based on our analysis, for a service-time distribution that is a mixture of deterministic and exponential, we numerically show that the dynamic policy has lower mean AoI than that of these two static policies and also that of the well studied M/GI/1/1 blocking system.


title: BERT-CoQAC: BERT-based Conversational Question Answering in Context

abstract: As one promising way to inquire about any particular information through a dialog with the bot, question answering dialog systems have gained increasing research interests recently. Designing interactive QA systems has always been a challenging task in natural language processing and used as a benchmark to evaluate a machine's ability of natural language understanding. However, such systems often struggle when the question answering is carried out in multiple turns by the users to seek more information based on what they have already learned, thus, giving rise to another complicated form called Conversational Question Answering (CQA). CQA systems are often criticized for not understanding or utilizing the previous context of the conversation when answering the questions. To address the research gap, in this paper, we explore how to integrate conversational history into the neural machine comprehension system. On one hand, we introduce a framework based on a publically available pre-trained language model called BERT for incorporating history turns into the system. On the other hand, we propose a history selection mechanism that selects the turns that are relevant and contributes the most to answer the current question. Experimentation results revealed that our framework is comparable in performance with the state-of-the-art models on the QuAC leader board. We also conduct a number of experiments to show the side effects of using entire context information which brings unnecessary information and noise signals resulting in a decline in the model's performance.