title: ALMA observations and modeling of the rotating outflow in Orion Source I

abstract: We present $^{29}$SiO(J=8--7) $\nu$=0, SiS (J=19--18) $\nu$=0, and $^{28}$SiO (J=8--7) $\nu$=1 molecular line archive observations made with the Atacama Large Millimeter/Submillimeter Array (ALMA) of the molecular outflow associated with Orion Source I. The observations show velocity asymmetries about the flow axis which are interpreted as outflow rotation. We find that the rotation velocity ($\sim$4--8 km s$^{-1}$) decreases with the vertical distance to the disk. In contrast, the cylindrical radius ($\sim$100--300 au), the expansion velocity ($\sim$2--15 km s$^{-1}$), and the axial velocity $v_{\rm z}$ ($\sim$-1--10 km s$^{-1}$) increase with the vertical distance. The mass estimated of the molecular outflow $\mathrm{M}_{\rm outflow}\sim$0.66--1.3 M$_\odot$. Given a kinematic time $\sim$130 yr, this implies a mass loss rate $\dot{\mathrm{M}}_{\rm outflow} \sim 5.1-10 \times 10^{-3}$ M$_\odot$ yr$^{-1}$. This massive outflow sets important contraints on disk wind models. We compare the observations with a model of a shell produced by the interaction between an anisotropic stellar wind and an Ulrich accretion flow that corresponds to a rotating molecular envelope in collapse. We find that the model cylindrical radii are consistent with the $^{29}$SiO(J=8--7) $\nu$=0 data. The expansion velocities and the axial velocities of the model are similar the observed values, except close to the disk ($z\sim\pm$150 au) for the expansion velocity. Nevertheless, the rotation velocities of the model are a factor $\sim$3--10 lower than the observed values. We conclude that the Ulrich flow alone cannot explain the rotation observed and other possibilities should be explored, like the inclusion of the angular momentum of a disk wind.


title: CrowS-Pairs: A Challenge Dataset for Measuring Social Biases in Masked Language Models

abstract: Pretrained language models, especially masked language models (MLMs) have seen success across many NLP tasks. However, there is ample evidence that they use the cultural biases that are undoubtedly present in the corpora they are trained on, implicitly creating harm with biased representations. To measure some forms of social bias in language models against protected demographic groups in the US, we introduce the Crowdsourced Stereotype Pairs benchmark (CrowS-Pairs). CrowS-Pairs has 1508 examples that cover stereotypes dealing with nine types of bias, like race, religion, and age. In CrowS-Pairs a model is presented with two sentences: one that is more stereotyping and another that is less stereotyping. The data focuses on stereotypes about historically disadvantaged groups and contrasts them with advantaged groups. We find that all three of the widely-used MLMs we evaluate substantially favor sentences that express stereotypes in every category in CrowS-Pairs. As work on building less biased models advances, this dataset can be used as a benchmark to evaluate progress.


title: Machine Learning in Airline Crew Pairing to Construct Initial Clusters for Dynamic Constraint Aggregation

abstract: The crew pairing problem (CPP) is generally modelled as a set partitioning problem where the flights have to be partitioned in pairings. A pairing is a sequence of flight legs separated by connection time and rest periods that starts and ends at the same base. Because of the extensive list of complex rules and regulations, determining whether a sequence of flights constitutes a feasible pairing can be quite difficult by itself, making CPP one of the hardest of the airline planning problems. In this paper, we first propose to improve the prototype Baseline solver of Desaulniers et al. (2020) by adding dynamic control strategies to obtain an efficient solver for large-scale CPPs: Commercial-GENCOL-DCA. These solvers are designed to aggregate the flights covering constraints to reduce the size of the problem. Then, we use machine learning (ML) to produce clusters of flights having a high probability of being performed consecutively by the same crew. The solver combines several advanced Operations Research techniques to assemble and modify these clusters, when necessary, to produce a good solution. We show, on monthly CPPs with up to 50 000 flights, that Commercial-GENCOL-DCA with clusters produced by ML-based heuristics outperforms Baseline fed by initial clusters that are pairings of a solution obtained by rolling horizon with GENCOL. The reduction of solution cost averages between 6.8% and 8.52%, which is mainly due to the reduction in the cost of global constraints between 69.79% and 78.11%.


title: Blaschke-Santalo inequality for many functions and geodesic barycenters of measures

abstract: Motivated by the geodesic barycenter problem from optimal transportation theory, we prove a natural generalization of the Blaschke-Santalo inequality and the affine isoperimetric inequalities for many sets and many functions. We derive from it an entropy bound for the total Kantorovich cost appearing in the barycenter problem. We also establish a "pointwise Prekopa-Leindler inequality" and show a monotonicity property of the multimarginal Blaschke-Santalo functional.


title: Switching induced by spin Hall effect in an in-plane magnetized ferromagnet with the easy axis parallel to the current

abstract: Magnetization switching in a fine-structured ferromagnet of nanoscale by the spin-transfer torque excited via the spin Hall effect has attracted much attention because it enables us to manipulate the magnetization without directly applying current to the ferromagnet. However, the switching mechanism is still unclear in regard to the ferromagnet having an in-plane easy axis parallel to the current. Here, we develop an analytical theory of the magnetization switching in this type of ferromagnet, and reveal the threshold current formulas for a deterministic switching. It is clarified that the current should be in between a certain range determined by two threshold currents because the spin-transfer torque due to a large current outside the range brings the magnetization in an energetically unstable state, and causes magnetization precession around the hard axis.


title: Efficient sampling from the Bingham distribution

abstract: We give a algorithm for exact sampling from the Bingham distribution $p(x)\propto \exp(x^\top A x)$ on the sphere $\mathcal S^{d-1}$ with expected runtime of $\operatorname{poly}(d, \lambda_{\max}(A)-\lambda_{\min}(A))$. The algorithm is based on rejection sampling, where the proposal distribution is a polynomial approximation of the pdf, and can be sampled from by explicitly evaluating integrals of polynomials over the sphere. Our algorithm gives exact samples, assuming exact computation of an inverse function of a polynomial. This is in contrast with Markov Chain Monte Carlo algorithms, which are not known to enjoy rapid mixing on this problem, and only give approximate samples. As a direct application, we use this to sample from the posterior distribution of a rank-1 matrix inference problem in polynomial time.


title: Approximation of the Double Travelling Salesman Problem with Multiple Stacks

abstract: The Double Travelling Salesman Problem with Multiple Stacks, DTSPMS, deals with the collect and delivery of n commodities in two distinct cities, where the pickup and the delivery tours are related by LIFO constraints. During the pickup tour, commodities are loaded into a container of k rows, or stacks, with capacity c. This paper focuses on computational aspects of the DTSPMS, which is NP-hard. We first review the complexity of two critical subproblems: deciding whether a given pair of pickup and delivery tours is feasible and, given a loading plan, finding an optimal pair of pickup and delivery tours, are both polynomial under some conditions on k and c. We then prove a (3k)/2 standard approximation for the MinMetrickDTSPMS, where k is a universal constant, and other approximation results for various versions of the problem. We finally present a matching-based heuristic for the 2DTSPMS, which is a special case with k=2 rows, when the distances are symmetric. This yields a 1/2-o(1), 3/4-o(1) and 3/2+o(1) standard approximation for respectively Max2DTSPMS, its restriction Max2DTSPMS-(1,2) with distances 1 and 2, and Min2DTSPMS-(1,2), and a 1/2-o(1) differential approximation for Min2DTSPMS and Max2DTSPMS.


title: Controlled Exciton-Plasmon Coupling in a Mixture of Ultrathin Periodically Aligned Single-Wall Carbon Nanotube Arrays

abstract: We study theoretically the in-plane electromagnetic response and the exciton-plasmon interactions for an experimentally feasible carbon nanotube (CN) film systems composed of parallel aligned periodic semiconducting CN arrays embedded in an ultrathin finite-thickness dielectric. For homogeneous single-CN films, the intertube coupling and thermal broadening bring the exciton and interband plasmon resonances closer together. They can even overlap due to the inhomogeneous broadening for films composed of array mixtures with a slight CN diameter distribution. In such systems the real part of the response function is negative for a broad range of energies (negative refraction band), and the CN film behaves as a hyperbolic metamaterial. We also show that for a properly fabricated two-component CN film, by varying the relative weights of the two constituent CN array components one can tune the optical absorption profile to make the film transmit or absorb light in the neighborhood of an exciton absorption resonance on-demand.


title: Event-Independent Network for Polyphonic Sound Event Localization and Detection

abstract: Polyphonic sound event localization and detection is not only detecting what sound events are happening but localizing corresponding sound sources. This series of tasks was first introduced in DCASE 2019 Task 3. In 2020, the sound event localization and detection task introduces additional challenges in moving sound sources and overlapping-event cases, which include two events of the same type with two different direction-of-arrival (DoA) angles. In this paper, a novel event-independent network for polyphonic sound event localization and detection is proposed. Unlike the two-stage method we proposed in DCASE 2019 Task 3, this new network is fully end-to-end. Inputs to the network are first-order Ambisonics (FOA) time-domain signals, which are then fed into a 1-D convolutional layer to extract acoustic features. The network is then split into two parallel branches. The first branch is for sound event detection (SED), and the second branch is for DoA estimation. There are three types of predictions from the network, SED predictions, DoA predictions, and event activity detection (EAD) predictions that are used to combine the SED and DoA features for on-set and off-set estimation. All of these predictions have the format of two tracks indicating that there are at most two overlapping events. Within each track, there could be at most one event happening. This architecture introduces a problem of track permutation. To address this problem, a frame-level permutation invariant training method is used. Experimental results show that the proposed method can detect polyphonic sound events and their corresponding DoAs. Its performance on the Task 3 dataset is greatly increased as compared with that of the baseline method.


title: Avoidance, Adjacency, and Association in Distributed Systems Design

abstract: Patterns of avoidance, adjacency, and association in complex systems design emerge from the system's underlying logical architecture (functional relationships among components) and physical architecture (component physical properties and spatial location). Understanding the physical--logical architecture interplay that gives rise to patterns of arrangement requires a quantitative approach that bridges both descriptions. Here, we show that statistical physics reveals patterns of avoidance, adjacency, and association across sets of complex, distributed system design solutions. Using an example arrangement problem and tensor network methods, we identify several phenomena in complex systems design, including placement symmetry breaking, propagating correlation, and emergent localization. Our approach generalizes straightforwardly to a broad range of complex systems design settings where it can provide a platform for investigating basic design phenomena.


title: Basics of observational entropy

abstract: These notes provide a brief primer on the basic aspects of "observational entropy" (also known as "quantum coarse-grained entropy"), a general framework for applying the concept of coarse-graining to quantum systems. We review the basic formalism, survey applications to thermodynamics, make a connection to quantum correlations and entanglement entropy, compare to the corresponding classical theory, and discuss a generalization based on POVM measurements.


title: Opinion dynamics on tie-decay networks

abstract: In social networks, interaction patterns typically change over time. We study opinion dynamics on tie-decay networks in which tie strength increases instantaneously when there is an interaction and decays exponentially between interactions. Specifically, we formulate continuous-time Laplacian dynamics and a discrete-time DeGroot model of opinion dynamics on these tie-decay networks, and we carry out numerical computations for the continuous-time Laplacian dynamics. We examine the speed of convergence by studying the spectral gaps of combinatorial Laplacian matrices of tie-decay networks. First, we compare the spectral gaps of the Laplacian matrices of tie-decay networks that we construct from empirical data with the spectral gaps for corresponding randomized and aggregate networks. We find that the spectral gaps for the empirical networks tend to be smaller than those for the randomized and aggregate networks. Second, we study the spectral gap as a function of the tie-decay rate and time. Intuitively, we expect small tie-decay rates to lead to fast convergence because the influence of each interaction between two nodes lasts longer for smaller decay rates. Moreover, as time progresses and more interactions occur, we expect eventual convergence. However, we demonstrate that the spectral gap need not decrease monotonically with respect to the decay rate or increase monotonically with respect to time. Our results highlight the importance of the interplay between the times that edges strengthen and decay in temporal networks.


title: Quantum hydrodynamics of spin winding

abstract: An easy-plane spin winding in a quantum spin chain can be treated as a transport quantity, which propagates along the chain but has a finite lifetime due to phase slips. In a hydrodynamic formulation for the winding dynamics, the quantum continuity equation acquires a source term due to the transverse vorticity flow. The latter reflects the phase slips and generally compromises the global conservation law. A linear-response formalism for the nonlocal winding transport then reduces to a Kubo response for the winding flow along the spin chain, in conjunction with the parasitic vorticity flow transverse to it. One-dimensional topological hydrodynamics can be recovered when the vorticity flow is asymptotically small. Starting with a microscopic spin-chain formulation, we focus on the asymptotic behavior of the winding transport based on the renormalized sine-Gordon equation, incorporating phase slips as well as Gilbert damping. A generic electrical device is proposed to manifest this physics. We thus suggest winding conductivity as a tangible concept that can characterize low-energy dynamics in a broad class of quantum magnets.


title: Entropy Regularization for Mean Field Games with Learning

abstract: Entropy regularization has been extensively adopted to improve the efficiency, the stability, and the convergence of algorithms in reinforcement learning. This paper analyzes both quantitatively and qualitatively the impact of entropy regularization for Mean Field Game (MFG) with learning in a finite time horizon. Our study provides a theoretical justification that entropy regularization yields time-dependent policies and, furthermore, helps stabilizing and accelerating convergence to the game equilibrium. In addition, this study leads to a policy-gradient algorithm for exploration in MFG. Under this algorithm, agents are able to learn the optimal exploration scheduling, with stable and fast convergence to the game equilibrium.


title: Beyond Equilibrium Temperature: How the Atmosphere/Interior Connection Affects the Onset of Methane, Ammonia, and Clouds in Warm Transiting Giant Planets

abstract: The atmospheric pressure-temperature profiles for transiting giant planets cross a range of chemical transitions. Here we show that the particular shape of these irradiated profiles for warm giant planets below 1300 K lead to striking differences in the behavior of non-equilibrium chemistry compared to brown dwarfs of similar temperatures. Our particular focus is H$_2$O, CO, CH$_4$, CO$_2$, and NH$_3$ in Jupiter- and Neptune-class planets. We show the cooling history of a planet, which depends most significantly on planetary mass and age, can have a dominant effect on abundances in the visible atmosphere, often swamping trends one might expect based on Teq alone. The onset of detectable CH$_4$ in spectra can be delayed to lower Teq for some planets compared to equilibrium, or pushed to higher Teq. The detectability of NH$_3$ is typically enhanced compared to equilibrium expectations, which is opposite to the brown dwarf case. We find that both CH$_4$ and NH$_3$ can become detectable at around the same Teq (at Teq values that vary with mass and metallicity) whereas these "onset" temperatures are widely spaced for brown dwarfs. We suggest observational strategies to search for atmospheric trends and stress that non-equilibrium chemistry and clouds can serve as probes of atmospheric physics. As examples of atmospheric complexity, we assess three Neptune-class planets GJ 436b, GJ 3470b, and WASP-107, all around Teq=700 K. Tidal heating due to eccentricity damping in all three planets heats the deep atmosphere by thousands of degrees, and may explain the absence of CH$_4$ in these cool atmospheres. Atmospheric abundances must be interpreted in the context of physical characteristics of the planet.


title: Extended chromatic symmetric functions and equality of ribbon Schur functions

abstract: We prove a general inclusion-exclusion relation for the extended chromatic symmetric function of a weighted graph, which specializes to (extended) $k$-deletion, and we give two methods to obtain numerous new bases from weighted graphs for the algebra of symmetric functions. Moreover, we classify when two weighted paths have equal extended chromatic symmetric functions by proving this is equivalent to the classification of equal ribbon Schur functions. This latter classification is dependent on the operation composition of compositions, which we generalize to composition of graphs. We then apply our generalization to obtain infinitely many families of weighted graphs whose members have equal extended chromatic symmetric functions.


title: DEEPMIR: A DEEP neural network for differential detection of cerebral Microbleeds and IRon deposits in MRI

abstract: Lobar cerebral microbleeds (CMBs) and localized non-hemorrhage iron deposits in the basal ganglia have been associated with brain aging, vascular disease and neurodegenerative disorders. Particularly, CMBs are small lesions and require multiple neuroimaging modalities for accurate detection. Quantitative susceptibility mapping (QSM) derived from in vivo magnetic resonance imaging (MRI) is necessary to differentiate between iron content and mineralization. We set out to develop a deep learning-based segmentation method suitable for segmenting both CMBs and iron deposits. We included a convenience sample of 24 participants from the MESA cohort and used T2-weighted images, susceptibility weighted imaging (SWI), and QSM to segment the two types of lesions. We developed a protocol for simultaneous manual annotation of CMBs and non-hemorrhage iron deposits in the basal ganglia. This manual annotation was then used to train a deep convolution neural network (CNN). Specifically, we adapted the U-Net model with a higher number of resolution layers to be able to detect small lesions such as CMBs from standard resolution MRI. We tested different combinations of the three modalities to determine the most informative data sources for the detection tasks. In the detection of CMBs using single class and multiclass models, we achieved an average sensitivity and precision of between 0.84-0.88 and 0.40-0.59, respectively. The same framework detected non-hemorrhage iron deposits with an average sensitivity and precision of about 0.75-0.81 and 0.62-0.75, respectively. Our results showed that deep learning could automate the detection of small vessel disease lesions and including multimodal MR data (particularly QSM) can improve the detection of CMB and non-hemorrhage iron deposits with sensitivity and precision that is compatible with use in large-scale research studies.


title: Regarding the Euler-Plateau Problem with Elastic Modulus

abstract: We study equilibrium configurations for the Euler-Plateau energy with elastic modulus, which couples an energy functional of Euler-Plateau type with a total curvature term often present in models for the free energy of biomembranes. It is shown that the potential minimizers of this energy are highly dependent on the choice of physical rigidity parameters, and that the area of critical surfaces can be computed entirely from their boundary data. When the elastic modulus does not vanish, it is shown that axially symmetric critical immersions and critical immersions of disk type are necessarily planar domains bounded by area-constrained elasticae. The cases of topological genus zero with multiple boundary components and unrestricted genus with control on the geodesic torsion are also discussed, and sufficient conditions are given which establish the same conclusion in these cases.


title: Learning from Mistakes: Combining Ontologies via Self-Training for Dialogue Generation

abstract: Natural language generators (NLGs) for task-oriented dialogue typically take a meaning representation (MR) as input. They are trained end-to-end with a corpus of MR/utterance pairs, where the MRs cover a specific set of dialogue acts and domain attributes. Creation of such datasets is labor-intensive and time-consuming. Therefore, dialogue systems for new domain ontologies would benefit from using data for pre-existing ontologies. Here we explore, for the first time, whether it is possible to train an NLG for a new larger ontology using existing training sets for the restaurant domain, where each set is based on a different ontology. We create a new, larger combined ontology, and then train an NLG to produce utterances covering it. For example, if one dataset has attributes for family-friendly and rating information, and the other has attributes for decor and service, our aim is an NLG for the combined ontology that can produce utterances that realize values for family-friendly, rating, decor and service. Initial experiments with a baseline neural sequence-to-sequence model show that this task is surprisingly challenging. We then develop a novel self-training method that identifies (errorful) model outputs, automatically constructs a corrected MR input to form a new (MR, utterance) training pair, and then repeatedly adds these new instances back into the training data. We then test the resulting model on a new test set. The result is a self-trained model whose performance is an absolute 75.4% improvement over the baseline model. We also report a human qualitative evaluation of the final model showing that it achieves high naturalness, semantic coherence and grammaticality


title: Shadow cast by a rotating and nonlinear magnetic-charged black hole in perfect fluid dark matter

abstract: We derived an exact solution of the spherically symmetric Hayward black hole surrounded by perfect fluid dark matter (PFDM). By applying the Newman-Janis algorithm, we generalized it to the corresponding rotating black hole. Then, we studied the shadows of rotating Hayward black hole in PFDM. The apparent shape of the shadow depends upon the black hole spin $a$, the magnetic charge $Q$ and the PFDM intensity parameter $k$ ($k<0$). The shadow is a perfect circle in the non-rotating case ($a=0$) and a deformed one in the rotating case ($a\neq{0}$). For a fixed value of $a$, the size of the shadow increases with the increasing $\vert{k}\vert$, but decreases with the increasing $Q$. We further investigated the black hole emission rate. We found that the emission rate decreases with the increasing $\vert{k}\vert$ (or $Q$) and the peak of the emission shifts to lower frequency. Finally, we discussed the observational prospects corresponding to the supermassive black hole $\mathrm{Sgr\ A^{*}}$ at the center of the Milky Way.


title: Efficient sorting, duplicate removal, grouping, and aggregation

abstract: Database query processing requires algorithms for duplicate removal, grouping, and aggregation. Three algorithms exist: in-stream aggregation is most efficient by far but requires sorted input; sort-based aggregation relies on external merge sort; and hash aggregation relies on an in-memory hash table plus hash partitioning to temporary storage. Cost-based query optimization chooses which algorithm to use based on several factors including the sort order of the input, input and output sizes, and the need for sorted output. For example, hash-based aggregation is ideal for output smaller than the available memory (e.g., TPC-H Query 1), whereas sorting the entire input and aggregating after sorting are preferable when both aggregation input and output are large and the output needs to be sorted for a subsequent operation such as a merge join. Unfortunately, the size information required for a sound choice is often inaccurate or unavailable during query optimization, leading to sub-optimal algorithm choices. In response, this paper introduces a new algorithm for sort-based duplicate removal, grouping, and aggregation. The new algorithm always performs at least as well as both traditional hash-based and traditional sort-based algorithms. It can serve as a system's only aggregation algorithm for unsorted inputs, thus preventing erroneous algorithm choices. Furthermore, the new algorithm produces sorted output that can speed up subsequent operations. Google's F1 Query uses the new algorithm in production workloads that aggregate petabytes of data every day.


title: Examining the rhetorical capacities of neural language models

abstract: Recently, neural language models (LMs) have demonstrated impressive abilities in generating high-quality discourse. While many recent papers have analyzed the syntactic aspects encoded in LMs, there has been no analysis to date of the inter-sentential, rhetorical knowledge. In this paper, we propose a method that quantitatively evaluates the rhetorical capacities of neural LMs. We examine the capacities of neural LMs understanding the rhetoric of discourse by evaluating their abilities to encode a set of linguistic features derived from Rhetorical Structure Theory (RST). Our experiments show that BERT-based LMs outperform other Transformer LMs, revealing the richer discourse knowledge in their intermediate layer representations. In addition, GPT-2 and XLNet apparently encode less rhetorical knowledge, and we suggest an explanation drawing from linguistic philosophy. Our method shows an avenue towards quantifying the rhetorical capacities of neural LMs.


title: Deformable Kernel Convolutional Network for Video Extreme Super-Resolution

abstract: Video super-resolution, which attempts to reconstruct high-resolution video frames from their corresponding low-resolution versions, has received increasingly more attention in recent years. Most existing approaches opt to use deformable convolution to temporally align neighboring frames and apply traditional spatial attention mechanism (convolution based) to enhance reconstructed features. However, such spatial-only strategies cannot fully utilize temporal dependency among video frames. In this paper, we propose a novel deep learning based VSR algorithm, named Deformable Kernel Spatial Attention Network (DKSAN). Thanks to newly designed Deformable Kernel Convolution Alignment (DKC_Align) and Deformable Kernel Spatial Attention (DKSA) modules, DKSAN can better exploit both spatial and temporal redundancies to facilitate the information propagation across different layers. We have tested DKSAN on AIM2020 Video Extreme Super-Resolution Challenge to super-resolve videos with a scale factor as large as 16. Experimental results demonstrate that our proposed DKSAN can achieve both better subjective and objective performance compared with the existing state-of-the-art EDVR on Vid3oC and IntVID datasets.


title: Statistics of a Family of Piecewise Linear Maps

abstract: We study statistical properties of the truncated flat spot map $f_t(x)$. In particular, we investigate whether for large $n$, the deviations $\sum_{i=0}^{n-1} \left(f_t^i(x_0)-\frac 12\right)$ upon rescaling satisfy a $Q$-Gaussian distribution if $x_0$ and $t$ are both independently and uniformly distributed on the unit circle. This was motivated by the fact that if $f_t$ is the rotation by $t$, then it has been shown that in this case the rescaled deviations are distributed as a $Q$-Gaussian with $Q=2$ (a Cauchy distribution). This is the only case where a non-trivial (i.e. $Q\neq 1$) $Q$-Gaussian has been analytically established in a conservative dynamical system. In this note, however, we prove that for the family considered here, $\lim_n S_n/n$ converges to a random variable with a curious distribution which is clearly not a $Q$-Gaussian or any other standard smooth distribution.


title: GeoD: Consensus-based Geodesic Distributed Pose Graph Optimization

abstract: We present a consensus-based distributed pose graph optimization algorithm for obtaining an estimate of the 3D translation and rotation of each pose in a pose graph, given noisy relative measurements between poses. The algorithm, called GeoD, implements a continuous time distributed consensus protocol to minimize the geodesic pose graph error. GeoD is distributed over the pose graph itself, with a separate computation thread for each node in the graph, and messages are passed only between neighboring nodes in the graph. We leverage tools from Lyapunov theory and multi-agent consensus to prove the convergence of the algorithm. We identify two new consistency conditions sufficient for convergence: pairwise consistency of relative rotation measurements, and minimal consistency of relative translation measurements. GeoD incorporates a simple one step distributed initialization to satisfy both conditions. We demonstrate GeoD on simulated and real world SLAM datasets. We compare to a centralized pose graph optimizer with an optimality certificate (SE-Sync) and a Distributed Gauss-Seidel (DGS) method. On average, GeoD converges 20 times more quickly than DGS to a value with 3.4 times less error when compared to the global minimum provided by SE-Sync. GeoD scales more favorably with graph size than DGS, converging over 100 times faster on graphs larger than 1000 poses. Lastly, we test GeoD on a multi-UAV vision-based SLAM scenario, where the UAVs estimate their pose trajectories in a distributed manner using the relative poses extracted from their on board camera images. We show qualitative performance that is better than either the centralized SE-Sync or the distributed DGS methods.


title: Universal Effectiveness of High-Depth Circuits in Variational Eigenproblems

abstract: We explore the effectiveness of variational quantum circuits in simulating the ground states of quantum many-body Hamiltonians. We show that generic high-depth circuits, performing a sequence of layer unitaries of the same form, can accurately approximate the desired states. We demonstrate their universal success by using two Hamiltonian systems with very different properties: the transverse field Ising model and the Sachdev-Ye-Kitaev model. The energy landscape of the high-depth circuits has a proper structure for the gradient-based optimization, i.e. the presence of local extrema -- near any random initial points -- reaching the ground level energy. We further test the circuit's capability of replicating random quantum states by minimizing the Euclidean distance.


title: Anisotropic spin distribution and perpendicular magnetic anisotropy in the layered ferromagnetic semiconductor (Ba,K)(Zn,Mn)$_{2}$As$_{2}$

abstract: Perpendicular magnetic anisotropy of the new ferromagnetic semiconductor (Ba,K)(Zn,Mn)$_{2}$As$_{2}$ is studied by angle-dependent x-ray magnetic circular dichroism measurements. The large magnetic anisotropy with the anisotropy field of 0.85 T is deduced by fitting the Stoner-Wohlfarth model to the magnetic-field-angle dependence of the projected magnetic moment. Transverse XMCD spectra highlights the anisotropic distribution of Mn 3$d$ electrons, where the $d_{xz}$ and $d_{yz}$ orbitals are less populated than the $d_{xy}$ state because of the $D_{2d}$ splitting arising from the elongated MnAs$_{4}$ tetrahedra. It is suggested that the magnetic anisotropy originates from the degeneracy lifting of $p$-$d_{xz}$, $d_{yz}$ hybridized states at the Fermi level and resulting energy gain due to spin-orbit coupling when spins are aligned along the $z$ direction.


title: Direct Observation of Ion Micromotion in a Linear Paul Trap

abstract: In this paper, direct observation of micromotion for multiple ions in a laser-cooled trapped ion crystal is discussed along with a novel measurement technique for micromotion amplitude. Micromotion is directly observed using a time-resolving, single-photon sensitive camera that provides both fluorescence and position data for each ion on the nanosecond time scale. Micromotion amplitude and phase for each ion in the crystal are measured, allowing this method to be sensitive to tilts and shifts of the ion chain from the null of the radiofrequency quadrupole potential in the linear trap. Spatial resolution makes this micromotion detection technique suitable for complex ion configurations, including two-dimensional geometries. It does not require any additional equipment or laser beams, and the modulation of the cooling lasers or trap voltages is not necessary for detection, as it is in other methods.


title: High-temperature superconductivity on the verge of a structural instability in lanthanum superhydride

abstract: A possibility of high, room-temperature superconductivity was predicted for metallic hydrogen in the 1960s. However, metallization and superconductivity of hydrogen are yet to be unambiguously demonstrated in the laboratory and may require pressures as high as 5 million atmospheres. Rare earth based "superhydrides" such as LaH10 can be considered a close approximation of metallic hydrogen even though they form at moderately lower pressures. In superhydrides the predominance of H-H metallic bonds and high superconducting transition temperatures bear the hallmarks of metallic hydrogen. Still, experimental studies revealing the key factors controlling their superconductivity are scarce. Here, we report on the pressure and magnetic field response of the superconducting order observed in LaH10. For LaH10 we find a correlation between superconductivity and a structural instability, strongly affecting the lattice vibrations responsible for the superconductivity.


title: Unknown Delay for Adversarial Bandit Setting with Multiple Play

abstract: This paper addresses the problem of unknown delays in adversarial multi-armed bandit (MAB) with multiple play. Existing work on similar game setting focused on only the case where the learner selects an arm in each round. However, there are lots of applications in robotics where a learner needs to select more than one arm per round. It is therefore worthwhile to investigate the effect of delay when multiple arms are chosen. The multiple arms chosen per round in this setting are such that they experience the same amount of delay. There can be an aggregation of feedback losses from different combinations of arms selected at different rounds, and the learner is faced with the challenge of associating the feedback losses to the arms producing them. To address this problem, this paper proposes a delayed exponential, exploitation and exploration for multiple play (DEXP3.M) algorithm. The regret bound is only slightly worse than the regret of DEXP3 already proposed for the single play setting with unknown delay.


title: Numerical computation of second order vacuum perturbations of Kerr black holes

abstract: Motivated by the desire to understand the leading order nonlinear gravitational wave interactions around arbitrarily rapidly rotating Kerr black holes, we describe a numerical code designed to compute second order vacuum perturbations on such spacetimes. A general discussion of the formalism we use is presented in (arXiv:2008.11770); here we show how we numerically implement that formalism with a particular choice of coordinates and tetrad conditions, and give example results for black holes with dimensionless spin parameters $a=0.7$ and $a=0.998$. We first solve the Teukolsky equation for the linearly perturbed Weyl scalar $\Psi_4^{(1)}$, followed by direct reconstruction of the spacetime metric from $\Psi_4^{(1)}$, and then solve for the dynamics of the second order perturbed Weyl scalar $\Psi_4^{(2)}$. This code is a first step toward a more general purpose second order code, and we outline how our basic approach could be further developed to address current questions of interest, including extending the analysis of ringdown in black hole mergers to before the linear regime, exploring gravitational wave "turbulence" around near-extremal Kerr black holes, and studying the physics of extreme mass ratio inspiral.


title: Bayesian Meta-reinforcement Learning for Traffic Signal Control

abstract: In recent years, there has been increasing amount of interest around meta reinforcement learning methods for traffic signal control, which have achieved better performance compared with traditional control methods. However, previous methods lack robustness in adaptation and stability in training process in complex situations, which largely limits its application in real-world traffic signal control. In this paper, we propose a novel value-based Bayesian meta-reinforcement learning framework BM-DQN to robustly speed up the learning process in new scenarios by utilizing well-trained prior knowledge learned from existing scenarios. This framework is based on our proposed fast-adaptation variation to Gradient-EM Bayesian Meta-learning and the fast-update advantage of DQN, which allows for fast adaptation to new scenarios with continual learning ability and robustness to uncertainty. The experiments on restricted 2D navigation and traffic signal control show that our proposed framework adapts more quickly and robustly in new scenarios than previous methods, and specifically, much better continual learning ability in heterogeneous scenarios.


title: Tuning the topology of $p$-wave superconductivity in an analytically solvable two-band model

abstract: We introduce and solve a two-band model of spinless fermions with $p_x$-wave pairing on a square lattice. The model reduces to the well-known extended Harper-Hofstadter model with half-flux quanta per plaquette and weakly coupled Kitaev chains in two respective limits. We show that its phase diagram contains a topologically nontrivial weak pairing phase as well as a trivial strong pairing phase as the ratio of the pairing amplitude and hopping is tuned. Introducing periodic driving to the model, we observe a cascade of Floquet phases with well defined quasienergy gaps and featuring chiral Majorana edge modes at the zero- or $\pi$-gap, or both. Dynamical topological invariants are obtained to characterize each phase and to explain the emergence of edge modes in the anomalous phase where all the quasienergy bands have zero Chern number. Analytical solution is achieved by exploiting a generalized mirror symmetry of the model, so that the effective Hamiltonian is decomposed into that of spin-$1/2$ in magnetic field, and the loop unitary operator becomes spin rotations. We further show the dynamical invariants manifest as the Hopf linking numbers.


title: Neighbourhood Bootstrap for Respondent-Driven Sampling

abstract: Respondent-Driven Sampling (RDS) is a form of link-tracing sampling, a sampling technique used for `hard-to-reach' populations that aims to leverage individuals' social relationships to reach potential participants. While the methodological focus has been restricted to the estimation of population proportions, there is a growing interest in the estimation of uncertainty for RDS as recent findings suggest that most variance estimators underestimate variability. Recently, Baraff et al. (2016) proposed the \textit{tree bootstrap} method based on resampling the RDS recruitment tree, and empirically showed that this method outperforms current bootstrap methods. However, some findings suggest that the tree bootstrap (severely) overestimates uncertainty. In this paper, we propose the \textit{neighbourhood} bootstrap method for quantifiying uncertainty in RDS. We prove the consistency of our method under some conditions and investigate its finite sample performance, through a simulation study, under realistic RDS sampling assumptions.


title: Uniqueness of shrinking gradient K\"ahler-Ricci solitons on non-compact toric manifolds

abstract: We show that, up to biholomorphism, there is at most one complete $T^n$-invariant shrinking gradient K\"ahler-Ricci soliton on a non-compact toric manifold $M$. We also establish uniqueness without assuming $T^n$-invariance if the Ricci curvature is bounded and if the soliton vector field lies in the Lie algebra $\mathfrak{t}$ of $T^n$. As an application, we show that, up to isometry, the unique complete shrinking gradient K\"ahler-Ricci soliton with bounded scalar curvature on $\mathbb{CP}^{1} \times \mathbb{C}$ is the standard product metric associated to the Fubini-Study metric on $\mathbb{CP}^{1}$ and the Euclidean metric on $\mathbb{C}$.


title: Lebesgue Measure Preserving Thompson's Monoid

abstract: This paper defines Lebesgue measure preserving Thompson's monoid, denoted by $\mathbb{G}$, which is modeled on Thompson's group $\mathbb{F}$ except that the elements of $\mathbb{G}$ are non-invertible. Moreover, it is required that the elements of $\mathbb{G}$ preserve Lebesgue measure. Monoid $\mathbb{G}$ exhibits very different properties from Thompson's group $\mathbb{F}$. The paper studies a number of algebraic (group-theoretic) and dynamical properties of $\mathbb{G}$ including approximation, mixing, periodicity, entropy, decomposition, generators, and topological conjugacy.


title: Multi-Agent Systems based on Contextual Defeasible Logic considering Focus

abstract: In this paper, we extend previous work on distributed reasoning using Contextual Defeasible Logic (CDL), which enables decentralised distributed reasoning based on a distributed knowledge base, such that the knowledge from different knowledge bases may conflict with each other. However, there are many use case scenarios that are not possible to represent in this model. One kind of such scenarios are the ones that require that agents share and reason with relevant knowledge when issuing a query to others. Another kind of scenarios are those in which the bindings among the agents (defined by means of mapping rules) are not static, such as in knowledge-intensive and dynamic environments. This work presents a multi-agent model based on CDL that not only allows agents to reason with their local knowledge bases and mapping rules, but also allows agents to reason about relevant knowledge (focus) -- which are not known by the agents a priori -- in the context of a specific query. We present a use case scenario, some formalisations of the model proposed, and an initial implementation based on the BDI (Belief-Desire-Intention) agent model.


title: Axion Dark Matter eXperiment: Detailed Design and Operations

abstract: Axion Dark Matter eXperiment (ADMX) ultra low noise haloscope technology has enabled the successful completion of two science runs (1A and 1B) that looked for dark matter axions in the $2.66$ to $3.1$ $\mu$eV mass range with Dine-Fischler-Srednicki-Zhitnisky (DFSZ) sensitivity Ref. [1,2]. Therefore, it is the most sensitive axion search experiment to date in this mass range. We discuss the technological advances made in the last several years to achieve this sensitivity, which includes the implementation of components, such as state-of-the-art quantum limited amplifiers and a dilution refrigerator. Furthermore, we demonstrate the use of a frequency tunable Microstrip Superconducting Quantum Interference Device (SQUID) Amplifier (MSA), in Run 1A, and a Josephson Parametric Amplifier (JPA), in Run 1B, along with novel analysis tools that characterize the system noise temperature.

categories: astro-ph.IM hep-ex physics.ins-det quant-ph


title: A Large Multi-Target Dataset of Common Bengali Handwritten Graphemes

abstract: Latin has historically led the state-of-the-art in handwritten optical character recognition (OCR) research. Adapting existing systems from Latin to alpha-syllabary languages is particularly challenging due to a sharp contrast between their orthographies. The segmentation of graphical constituents corresponding to characters becomes significantly hard due to a cursive writing system and frequent use of diacritics in the alpha-syllabary family of languages. We propose a labeling scheme based on graphemes (linguistic segments of word formation) that makes segmentation in-side alpha-syllabary words linear and present the first dataset of Bengali handwritten graphemes that are commonly used in an everyday context. The dataset contains 411k curated samples of 1295 unique commonly used Bengali graphemes. Additionally, the test set contains 900 uncommon Bengali graphemes for out of dictionary performance evaluation. The dataset is open-sourced as a part of a public Handwritten Grapheme Classification Challenge on Kaggle to benchmark vision algorithms for multi-target grapheme classification. The unique graphemes present in this dataset are selected based on commonality in the Google Bengali ASR corpus. From competition proceedings, we see that deep-learning methods can generalize to a large span of out of dictionary graphemes which are absent during training. Dataset and starter codes at www.kaggle.com/c/bengaliai-cv19.


title: Lowering Helstrom Bound with non-standard coherent states

abstract: In quantum information processing, {using a receiver device to differentiate between two nonorthogonal states leads to a quantum error probability. The minimum possible error is} known as the Helstrom bound. In this work we study and compare quantum limits for states which generalize the Glauber-Sudarshan coherent states, like non-linear, Perelomov, Barut-Girardello, and (modified) Susskind-Glogower coherent states. For some of these, we show that the Helstrom bound can be significantly lowered and even vanish in specific regimes.


title: Renormalization of Twisted Ramond Fields in D1-D5 SCFT$_2$

abstract: We explore the Ramond sector of the deformed two-dimensional $\cal N = (4, 4)$ superconformal $(T^4)^N /S_N$ orbifold theory, describing bound states of D1-D5 brane system in type IIB superstring. We derive the large-$N$ limit of the four-point function of two R-charged twisted Ramond fields and two marginal deformation operators at the free orbifold point. Specific short-distance limits of this function provide several structure constants, the OPE fusion rules and the conformal dimensions of a few non-BPS operators. The second order correction (in the deformation parameter) to the two-point function of the Ramond fields, defined as double integrals over this four-point function, turns out to be UV-divergent, requiring an appropriate renormalization of the fields. We calculate the corrections to the conformal dimensions of the twisted Ramond ground states at the large-$N$ limit. The same integral yields the first-order deviation from zero of the structure constant of the three-point function of two Ramond fields and one deformation operator. Similar results concerning the correction to the two-point function of bare twist operators and their renormalization are also obtained.


title: Isothermal spheres in general relativity and Horava-type gravity

abstract: We construct a toy model for isothermal spheres in Horava gravity, which includes Einstein's gravity if a parameter is appropriately chosen. The equations for the isothermal spheres are derived from the partition function of the gravitating particle system. We confirm that the Newtonian limit of the system coincides with the model of the well-known isothermal sphere. The stability of the isothermal sphere is found to be sensitive to the energy density at the center of the sphere.


title: Information Propagation Model in Hybrid Networks

abstract: It is in practice impossible to describe the topology of a real network or its message propagation process using a single dynamic model. To address this issue, we constructed a new hybrid network model based on scale-free (SF), small-world (SW) features that functions as closely as possible to a real network. And the hybrid propagation model is constructed with susceptible-infected-susceptible (SIS), susceptible-infected-recovered (SIR) and susceptible-infected-recovered-susceptible (SIRS) model mixed in arbitrary proportions. The model applies the concepts of blockbuster effect and implicit node edges to reflect explosive spread as a significant characteristic of information propagation. A theoretical analysis and derivation of the new model in which hybrid networks were simulated revealed that the network degree distribution closely follows a power law. Using an improved similarity function to define the degree of closeness to real network cases, the proposed model was shown to be valid and very close to a real network.


title: C-Arm Non-Circular Orbits: Geometric Calibration, Image Quality, and Avoidance of Metal Artifacts

abstract: Metal artifacts present a frequent challenge to cone-beam CT (CBCT) in image-guided surgery, obscuring visualization of metal instruments and adjacent anatomy. Recent advances in mobile C-arm systems have enabled 3D imaging capacity with non-circular orbits. We extend a previously proposed metal artifacts avoidance (MAA) method to reduce the influence of metal artifacts by prospectively defining a non-circular orbit that avoids metal-induced biases in projection domain. Accurate geometric calibration is an important challenge to accurate 3D image reconstruction for such orbits. We investigate the performance of interpolation-based calibration from a library of circular orbits for any non-circular orbit. We apply the method to non-circular scans acquired for MAA, which involves: (i) coarse 3D localization of metal objects via only two scout views using an end-to-end trained neural network; (ii) calculation of the metal-induced x-ray spectral shift for all possible views; and (iii) identification of the non-circular orbit that minimizes the variations in spectral shift. Non-circular orbits with interpolation-based geometric calibration yielded reasonably accurate 3D image reconstruction. The end-to-end neural network accurately localized metal implants with just two scout views even in complex anatomical scenes, improving Dice coefficient by ~42% compared to a more conventional cascade of separately trained U-nets. In a spine phantom with pedicle screw instrumentation, non-circular orbits identified by the MAA method reduced the magnitude of metal "blomming" artifacts (apparent width of the screw shaft) in CBCT reconstructions by ~70%. The proposed imaging and calibration methods present a practical means to improve image quality in mobile C-arm CBCT by identifying non-circular scan protocols that improve sampling and reduce metal-induced biases in the projection data.


title: Towards Self-learning Edge Intelligence in 6G

abstract: Edge intelligence, also called edge-native artificial intelligence (AI), is an emerging technological framework focusing on seamless integration of AI, communication networks, and mobile edge computing. It has been considered to be one of the key missing components in the existing 5G network and is widely recognized to be one of the most sought-after functions for tomorrow's wireless 6G cellular systems. In this article, we identify the key requirements and challenges of edge-native AI in 6G. A self-learning architecture based on self-supervised Generative Adversarial Nets (GANs) is introduced to \blu{demonstrate the potential performance improvement that can be achieved by automatic data learning and synthesizing at the edge of the network}. We evaluate the performance of our proposed self-learning architecture in a university campus shuttle system connected via a 5G network. Our result shows that the proposed architecture has the potential to identify and classify unknown services that emerge in edge computing networks. Future trends and key research problems for self-learning-enabled 6G edge intelligence are also discussed.


title: Combinatorial invariants for nets of conics in $\text{PG}(2,q)$

abstract: The problem of classifying linear systems of conics in projective planes dates back at least to Jordan, who classified pencils (one-dimensional systems) of conics over $\mathbb{C}$ and $\mathbb{R}$ in 1906--1907. The analogous problem for finite fields $\mathbb{F}_q$ with $q$ odd was solved by Dickson in 1908. In 1914, Wilson attempted to classify nets (two-dimensional systems) of conics over finite fields of odd characteristic, but his classification was incomplete and contained some inaccuracies. In a recent article, we completed Wilson's classification of nets of rank one, namely those containing a repeated line. The aim of the present paper is to introduce and calculate certain combinatorial invariants of these nets, which we expect will be of use in various applications. Our approach is geometric in the sense that we view a net of rank one as a plane in $\text{PG}(5,q)$ that meets the quadric Veronesean in at least one point; two such nets are then equivalent if and only if the corresponding planes belong to the same orbit under the induced action of $\text{PGL}(3,q)$ viewed as a subgroup of $\text{PGL}(6,q)$. We have previously determined the orbits of lines in $\text{PG}(5,q)$ under this action, which correspond to the aforementioned pencils of conics in $\text{PG}(2,q)$. The main contribution of this paper is to determine the line-orbit distribution of a plane $\pi$ corresponding to a net of rank one, namely, the number of lines in $\pi$ belonging to each line orbit. It turns out that this list of invariants completely determines the orbit of $\pi$, and we will use this fact in forthcoming work to develop an efficient algorithm for calculating the orbit of a given net of rank one. As a more immediate application, we also determine the stabilisers of nets of rank one in $\text{PGL}(3,q)$, and hence the orbit sizes.


title: Training Data Augmentation for Deep Learning Radio Frequency Systems

abstract: Applications of machine learning are subject to three major components that contribute to the final performance metrics. Within the category of neural networks, and deep learning specifically, the first two are the architecture for the model being trained and the training approach used. This work focuses on the third component, the data used during training. The primary questions that arise are ``what is in the data'' and ``what within the data matters?'' Looking into the Radio Frequency Machine Learning (RFML) field of Automatic Modulation Classification (AMC) as an example of a tool used for situational awareness, the use of synthetic, captured, and augmented data are examined and compared to provide insights about the quantity and quality of the available data necessary to achieve desired performance levels. There are three questions discussed within this work: (1) how useful a synthetically trained system is expected to be when deployed without considering the environment within the synthesis, (2) how can augmentation be leveraged within the RFML domain, and lastly, (3) what impact knowledge of degradations to the signal caused by the transmission channel contributes to the performance of a system. In general, the examined data types each have useful contributions to a final application, but captured data germane to the intended use case will always provide more significant information and enable the greatest performance. Despite the benefit of captured data, the difficulties and costs that arise from live collection often make the quantity of data needed to achieve peak performance impractical. This paper helps quantify the balance between real and synthetic data, offering concrete examples where training data is parametrically varied in size and source.


title: System Design and Analysis for Energy-Efficient Passive UAV Radar Imaging System using Illuminators of Opportunity

abstract: Unmanned aerial vehicle (UAV) can provide superior flexibility and cost-efficiency for modern radar imaging systems, which is an ideal platform for advanced remote sensing applications using synthetic aperture radar (SAR) technology. In this paper, an energy-efficient passive UAV radar imaging system using illuminators of opportunity is first proposed and investigated. Equipped with a SAR receiver, the UAV platform passively reuses the backscattered signal of the target scene from an external illuminator, such as SAR satellite, GNSS or ground-based stationary commercial illuminators, and achieves bi-static SAR imaging and data communication. The system can provide instant accessibility to the radar image of the interested targets with enhanced platform concealment, which is an essential tool for stealth observation and scene monitoring. The mission concept and system block diagram are first presented with justifications on the advantages of the system. Then, a set of mission performance evaluators is established to quantitatively assess the capability of the system in a comprehensive manner, including UAV navigation, passive SAR imaging and communication. Finally, the validity of the proposed performance evaluators are verified by numerical simulations.


title: Quantum-inspired search method for low-energy states of classical Ising Hamiltonians

abstract: We develop a quantum-inspired numerical procedure for searching low-energy states of a classical Hamiltonian composed of two-body fully-connected random Ising interactions and a random local longitudinal magnetic field. In this method, we introduce infinitesimal quantum interactions that do not commute with the original Ising Hamiltonian, and repeatedly generate and truncate direct product states, inspired by the Krylov subspace method, to obtain the low-energy states of the original classical Ising Hamiltonian. The computational cost is controlled by the form of infinitesimal quantum interactions (e.g., one-body or two-body interactions) and the numbers of infinitesimal interaction terms introduced, different initial states considered, and low-energy states kept during the iteration. For a demonstrate of the method, here we introduce as the infinitesimal quantum interactions pair products of Pauli $X$ operators acting on different sites and on-site Pauli $X$ operators into the random Ising Hamiltonian, in which the numerical cost is $O(N^3)$ per iteration with the system size $N$. We consider 120 instances of the random coupling realizations for the random Ising Hamiltonian with $N$ up to 600 and search the 120 lowest-energy states for each instance. We find that the time-to-solution by the quantum-inspired method proposed here, with parallelization in terms of the different initial states, for searching the ground state of the random Ising Hamiltonian scales approximately as $N^5$ for $N$ up to 600. We also examine the basic physical properties such as the ensemble-averaged ground-state and first-excited energies and the ensemble-averaged number of states in the low-energy region of the random Ising Hamiltonian.


title: Estimation in exponential family Regression based on linked data contaminated by mismatch error

abstract: Identification of matching records in multiple files can be a challenging and error-prone task. Linkage error can considerably affect subsequent statistical analysis based on the resulting linked file. Several recent papers have studied post-linkage linear regression analysis with the response variable in one file and the covariates in a second file from the perspective of the "Broken Sample Problem" and "Permuted Data". In this paper, we present an extension of this line of research to exponential family response given the assumption of a small to moderate number of mismatches. A method based on observation-specific offsets to account for potential mismatches and $\ell_1$-penalization is proposed, and its statistical properties are discussed. We also present sufficient conditions for the recovery of the correct correspondence between covariates and responses if the regression parameter is known. The proposed approach is compared to established baselines, namely the methods by Lahiri-Larsen and Chambers, both theoretically and empirically based on synthetic and real data. The results indicate that substantial improvements over those methods can be achieved even if only limited information about the linkage process is available.


title: Dual Attention Model for Citation Recommendation

abstract: Based on an exponentially increasing number of academic articles, discovering and citing comprehensive and appropriate resources has become a non-trivial task. Conventional citation recommender methods suffer from severe information loss. For example, they do not consider the section of the paper that the user is writing and for which they need to find a citation, the relatedness between the words in the local context (the text span that describes a citation), or the importance on each word from the local context. These shortcomings make such methods insufficient for recommending adequate citations to academic manuscripts. In this study, we propose a novel embedding-based neural network called "dual attention model for citation recommendation (DACR)" to recommend citations during manuscript preparation. Our method adapts embedding of three dimensions of semantic information: words in the local context, structural contexts, and the section on which a user is working. A neural network is designed to maximize the similarity between the embedding of the three input (local context words, section and structural contexts) and the target citation appearing in the context. The core of the neural network is composed of self-attention and additive attention, where the former aims to capture the relatedness between the contextual words and structural context, and the latter aims to learn the importance of them. The experiments on real-world datasets demonstrate the effectiveness of the proposed approach.


title: Quantized critical supercurrent in SrTiO$_3$-based quantum point contacts

abstract: Superconductivity in SrTiO$_3$ occurs at remarkably low carrier densities and therefore, unlike conventional superconductors, can be controlled by electrostatic gates. Here we demonstrate nanoscale weak links connecting superconducting leads, all within a single material, SrTiO$_3$. Ionic liquid gating accumulates carriers in the leads, and local electrostatic gates are tuned to open the weak link. These devices behave as superconducting quantum point contacts with a quantized critical supercurrent. This is a milestone towards establishing SrTiO$_3$ as a single-material platform for mesoscopic superconducting transport experiments, that also intrinsically contains the necessary ingredients to engineer topological superconductivity


title: Filtration Games and Potentially Projective Modules

abstract: The notion of a \textbf{$\boldsymbol{\mathcal{C}}$-filtered} object, where $\mathcal{C}$ is some (typically small) collection of objects in a Grothendieck category, has become ubiquitous since the solution of the Flat Cover Conjecture around the year 2000. We introduce the \textbf{$\boldsymbol{\mathcal{C}}$-Filtration Game of length $\boldsymbol{\omega_1}$} on a module, paying particular attention to the case where $\mathcal{C}$ is the collection of all countably presented, projective modules. We prove that Martin's Maximum implies the determinacy of many $\mathcal{C}$-Filtration Games of length $\omega_1$, which in turn imply the determinacy of certain Ehrenfeucht-Fra\"iss\'{e} games of length $\omega_1$; this allows a significant strengthening of a theorem of Mekler-Shelah-Vaananen \cite{MR1191613}. Also, Martin's Maximum implies that if $R$ is a countable hereditary ring, the class of \textbf{$\boldsymbol{\sigma}$-closed potentially projective modules} -- i.e., those modules that are projective in some $\sigma$-closed forcing extension of the universe -- is closed under $<\aleph_2$-directed limits. We also give an example of a (ZFC-definable) class of abelian groups that, under the ordinary subgroup relation, constitutes an Abstract Elementary Class (AEC) with L\"owenheim-Skolem number $\aleph_1$ in some models in set theory, but fails to be an AEC in other models of set theory.


title: Resolved spectral variations of the centimetre-wavelength continuum from the rho Oph W photo-dissociation-region

abstract: Cm-wavelength radio continuum emission in excess of free-free, synchrotron and Rayleigh-Jeans dust emission (excess microwave emission, EME), and often called `anomalous microwave emission', is bright in molecular cloud regions exposed to UV radiation, i.e. in photo-dissociation regions (PDRs). The EME correlates with IR dust emission on degree angular scales. Resolved observations of well-studied PDRs are needed to compare the spectral variations of the cm-continuum with tracers of physical conditions and of the dust grain population. The EME is particularly bright in the regions of the rho Ophiuchi molecular cloud (rho Oph) that surround the earliest type star in the complex, HD 147889, where the peak signal stems from the filament known as the rho Oph-W PDR. Here we report on ATCA observations of rho Oph-W that resolve the width of the filament. We recover extended emission using a variant of non-parametric image synthesis performed in the sky plane. The multi-frequency 17 GHz to 39 GHz mosaics reveal spectral variations in the cm-wavelength continuum. At ~30 arcsec resolutions, the 17-20 GHz intensities follow tightly the mid-IR, Icm propto I(8 um), despite the breakdown of this correlation on larger scales. However, while the 33-39 GHz filament is parallel to IRAC 8 mum, it is offset by 15-20 arcsec towards the UV source. Such morphological differences in frequency reflect spectral variations, which we quantify spectroscopically as a sharp and steepening high-frequency cutoff, interpreted in terms of the spinning dust emission mechanism as a minimum grain size a_cutoff ~ 6 +- 1A that increases deeper into the PDR.


title: A parallel subgradient projection algorithm for quasiconvex equilibrium problems under the intersection of convex sets

abstract: In this paper, we studied the equilibrium problem where the bi-function may be quasiconvex with respect to the second variable and the feasible set is the intersection of a finite number of convex sets. We propose a projection-algorithm, where the projection can be computed independently onto each component set. The convergence of the algorithm is investigated and numerical examples for a variational inequality problem involving affine fractional operator are provided to demonstrate the behavior of the algorithm.


title: Who Are the `Silent Spreaders'?: Contact Tracing in Spatio-Temporal Memory Models

abstract: The COVID-19 epidemic has swept the world for over a year. However, a large number of infectious asymptomatic COVID-19 cases (\textit{ACC}s) are still making the breaking up of the transmission chains very difficult. Efforts by epidemiological researchers in many countries have thrown light on the clinical features of ACCs, but there is still a lack of practical approaches to detect ACCs so as to help contain the pandemic. To address the issue of ACCs, this paper presents a neural network model called Spatio-Temporal Episodic Memory for COVID-19 (\textit{STEM-COVID}) to identify ACCs from contact tracing data. Based on the fusion Adaptive Resonance Theory (\textit{ART}), the model encodes a collective spatio-temporal episodic memory of individuals and incorporates an effective mechanism of parallel searches for ACCs. Specifically, the episodic traces of the identified positive cases are used to map out the episodic traces of suspected ACCs using a weighted evidence pooling method. To evaluate the efficacy of STEM-COVID, a realistic agent based simulation model for COVID-19 spreading is implemented based on the recent epidemiological findings on ACCs. The experiments based on rigorous simulation scenarios, manifesting the current situation of COVID-19 spread, show that the STEM-COVID model with weighted evidence pooling has a higher level of accuracy and efficiency for identifying ACCs when compared with several baselines. Moreover, the model displays strong robustness against noisy data and different ACC proportions, which partially reflects the effect of breakthrough infections after vaccination on the virus transmission.


title: Tristable flow states and reversal of the large-scale circulation in two-dimensional circular convection cells

abstract: We present a numerical study of the flow states and reversals of the large-scale circulation (LSC) in a two-dimensional circular Rayleigh-B\'enard cell. Long-time direct numerical simulations are carried out in the Rayleigh number ($Ra$) range $10^{7} \le Ra \le 10^{8}$ and Prandtl number ($Pr$) range $2.0 \le Pr \le 20.0$. We found that a new, long-lived, chaotic flow state exists, in addition to the commonly observed circulation states (the LSC in the clockwise and counterclockwise directions). The circulation states consist of one primary roll in the middle and two secondary rolls near the top and bottom circular walls. The primary roll becomes stronger and larger, while the two secondary rolls diminish, with increasing $Ra$. Our results suggest that the reversal of the LSC is accompanied by the secondary rolls growing, breaking the primary roll and then connecting to form a new primary roll with reversed direction. We mapped out the phase diagram of the existence of the LSC and the reversal in the $Ra$-$Pr$ space, which reveals that the flow is in the circulation states when $Ra$ is large and $Pr$ is small. The reversal of the LSC can only occur in a limited $Pr$ range. The phase diagram can be understood in terms of competition between the thermal and viscous diffusions. We also found that the internal flow states manifested themselves into global properties such as Nusselt and Reynolds numbers.


title: Cosmology with the Einstein Telescope: No Slip Gravity Model and Redshift Specifications

abstract: The Einstein Telescope and other third generation interferometric detectors of gravitational waves are projected to be operational post $2030$. The cosmological signatures of gravitational waves would undoubtedly shed light on any departure from the current gravitational framework. We here confront a specific modified gravity model, the No Slip Gravity model, with forecast observations of gravitational waves. We compare the predicted constraints on the dark energy equation of state parameters $w_0^{}-w_a^{}$, between the modified gravity model and that of Einstein gravity. We show that the No Slip Gravity model mimics closely the constraints from the standard gravitational theory, and that the cosmological constraints are very similar. The use of spectroscopic redshifts, especially in the low--redshift regime, lead to significant improvements in the inferred parameter constraints. We test how well such a prospective gravitational wave dataset would function at testing such models, and find that there are significant degeneracies between the modified gravity model parameters, and the cosmological parameters that determine the distance, due to the gravitational wave dimming effect of the modified theory.


title: A Compare Aggregate Transformer for Understanding Document-grounded Dialogue

abstract: Unstructured documents serving as external knowledge of the dialogues help to generate more informative responses. Previous research focused on knowledge selection (KS) in the document with dialogue. However, dialogue history that is not related to the current dialogue may introduce noise in the KS processing. In this paper, we propose a Compare Aggregate Transformer (CAT) to jointly denoise the dialogue context and aggregate the document information for response generation. We designed two different comparison mechanisms to reduce noise (before and during decoding). In addition, we propose two metrics for evaluating document utilization efficiency based on word overlap. Experimental results on the CMUDoG dataset show that the proposed CAT model outperforms the state-of-the-art approach and strong baselines.


title: Constraint Monotonicity, Epistemic Splitting and Foundedness Could in General Be Too Strong in Answer Set Programming

abstract: Recently, the notions of subjective constraint monotonicity, epistemic splitting, and foundedness have been introduced for epistemic logic programs, with the aim to use them as main criteria respectively intuitions to compare different answer set semantics proposed in the literature on how they comply with these intuitions. In this note, we consider these three notions and demonstrate on some examples that they may be too strong in general and may exclude some desired answer sets respectively world views. In conclusion, these properties should not be regarded as mandatory properties that every answer set semantics must satisfy in general.


title: An Inverse Problem on Determining Second Order Symmetric Tensor for Perturbed Biharmonic Operator

abstract: This article offers a study of the Calder\'on type inverse problem of determining up to second order coefficients of the higher order elliptic operator. Here we show that it is possible to determine an anisotropic second order perturbation given by a symmetric matrix, along with a first order perturbation given by a vector field and a zero-th order potential function inside a bounded domain by measuring the Dirichlet to Neumann map of the perturbed biharmonic operator on the boundary of that domain.


title: Impurity Band Conduction in Si-doped \b{eta}-Ga2O3 Films

abstract: By combining temperature-dependent resistivity and Hall effect measurements, we investigate donor state energy in Si-doped \b{eta}-Ga2O3 films grown using metal-organic vapor phase epitaxy (MOVPE). High magnetic field Hall effect measurements (H = +/-90 kOe) showed non-linear Hall resistance for T < 150 K revealing two-band conduction. Further analyses revealed carrier freeze-out characteristics in both bands yielding donor state energies of ~ 33.7 and ~ 45.6 meV. The former is consistent with the donor energy of Si in \b{eta}-Ga2O3 whereas the latter suggests a residual donor state, likely associated with a DX center. This study provides a critical insight into the impurity band conduction and the defect energy states in \b{eta}-Ga2O3 using high-field magnetotransport measurements.


title: Concentration inequality around the thermal equilibrium measure of Coulomb gases

abstract: This article deals with Coulomb gases at an intermediate temperature regime, in which no structure is observed at the microscopic level, but the mass in confined to a compact set. Our main result is a concentration inequality around the thermal equilibrium measure, stating that with probability exponentially close to $1,$ the empirical measure is $\mathcal{O}\left( \frac{1}{N^{\frac{1}{d}}}\right)$ close to the thermal equilibrium measure. We also prove that this concentration inequality is optimal in some sense. The main new tool are functional inequalities that allow us to compare the bounded Lipschitz norm of a measure to its $H^{-1}$ norm in some cases when the measure does not have compact support.


title: BiLiMO: Bit-Limited MIMO Radar via Task-Based Quantization

abstract: Recent years have witnessed growing interest in reduced cost radar systems operating with low power. Multiple-input multiple-output (MIMO) radar technology is known to achieve high performance sensing by probing with multiple orthogonal waveforms. However, implementing a low cost low power MIMO radar is challenging. One of the reasons for this difficulty stems from the increased cost and power consumption required by analog-to-digital convertors (ADCs) in acquiring the multiple waveforms at the radar receiver. In this work we study reduced cost MIMO radar receivers restricted to operate with low resolution ADCs. We design bit-limited MIMO radar (BiLiMO) receivers which are capable of accurately recovering their targets while operating under strict resolution constraints. This is achieved by applying an additional analog filter to the acquired waveforms, and designing the overall hybrid analog-digital system to facilitate target identification using task-based quantization methods. In particular, we exploit the fact that the target parameters can be recovered from a compressed representation of the received waveforms. We thus tune the acquisition system to recover this representation with minimal distortion, from which the targets can be extracted in digital, and characterize the achievable error in identifying the targets. Our numerical results demonstrate that the proposed BiLiMO receiver operating with a bit budget of one bit per sample achieves target recovery performance which approaches that of costly MIMO radars operating with unlimited resolution ADCs, while substantially outperforming MIMO receivers operating only in the digital domain under the same bit limitations.


title: CAPTURE: A continuum imaging pipeline for the uGMRT

abstract: We present the first fully automated pipeline for making images from the interferometric data obtained from the upgraded Giant Metrewave Radio Telescope (uGMRT) called CAsa Pipeline-cum-Toolkit for Upgraded Giant Metrewave Radio Telescope data REduction - CAPTURE. It is a python program that uses tasks from the NRAO Common Astronomy Software Applications (CASA) to perform the steps of flagging of bad data, calibration, imaging and self-calibration. The salient features of the pipeline are: i) a fully automatic mode to go from the raw data to a self-calibrated continuum image, ii) specialized flagging strategies for short and long baselines that ensure minimal loss of extended structure, iii) flagging of persistent narrow band radio frequency interference (RFI), iv) flexibility for the user to configure the pipeline for step-by-step analysis or special cases and v) analysis of data from the legacy GMRT. CAPTURE is available publicly on github (https://github.com/ruta-k/uGMRT-pipeline, release v1.0.0). The primary beam correction for the uGMRT images produced with CAPTURE is made separately available at https://github.com/ruta-k/uGMRTprimarybeam. We show examples of using CAPTURE on uGMRT and legacy GMRT data. In principle, CAPTURE can be tailored for use with radio interferometric data from other telescopes.


title: Bronnikov-like Wormholes in Einstein-Scalar Gravity

abstract: In this paper, we analyse the global structure of the Bronnikov wormhole, which is the most general spherically-symmetric and static solution in Einstein gravity coupled to a free massless phantom scalar. We then introduce a scalar potential and construct a large class of exact solutions that can be viewed as generalizations of the Bronnikov wormhole. We study the global structure and classify the parameters of these new wormholes. For suitable parameters, some are regular black holes with a bouncing de Sitter spacetime inside the event horizon.


title: Improving Vietnamese Named Entity Recognition from Speech Using Word Capitalization and Punctuation Recovery Models

abstract: Studies on the Named Entity Recognition (NER) task have shown outstanding results that reach human parity on input texts with correct text formattings, such as with proper punctuation and capitalization. However, such conditions are not available in applications where the input is speech, because the text is generated from a speech recognition system (ASR), and that the system does not consider the text formatting. In this paper, we (1) presented the first Vietnamese speech dataset for NER task, and (2) the first pre-trained public large-scale monolingual language model for Vietnamese that achieved the new state-of-the-art for the Vietnamese NER task by 1.3% absolute F1 score comparing to the latest study. And finally, (3) we proposed a new pipeline for NER task from speech that overcomes the text formatting problem by introducing a text capitalization and punctuation recovery model (CaPu) into the pipeline. The model takes input text from an ASR system and performs two tasks at the same time, producing proper text formatting that helps to improve NER performance. Experimental results indicated that the CaPu model helps to improve by nearly 4% of F1-score.


title: Fast Uplink Grant-Free NOMA with Sinusoidal Spreading Sequences

abstract: Uplink (UL) dominated sporadic transmission and stringent latency requirement of massive machine type communication (mMTC) forces researchers to abandon complicated grant-acknowledgment based legacy networks. UL grant-free non-orthogonal multiple access (NOMA) provides an array of features which can be harnessed to efficiently solve the problem of massive random connectivity and latency. Because of the inherent sparsity in user activity pattern in mMTC, the trend of existing literature specifically revolves around compressive sensing based multi user detection (CS-MUD) and Bayesian framework paradigm which employs either random or Zadoff-Chu spreading sequences for non-orthogonal multiple access. In this work, we propose sinusoidal code as candidate spreading sequences. We show that, sinusoidal codes allow some non-iterative algorithms to be employed in context of active user detection, channel estimation and data detection in a UL grant-free mMTC system. This relaxes the requirement of several impractical assumptions considered in the state-of-art algorithms with added advantages of performance guarantees and lower computational cost. Extensive simulation results validate the performance potential of sinusoidal codes in realistic mMTC environments.


title: RRF102: Meeting the TREC-COVID Challenge with a 100+ Runs Ensemble

abstract: In this paper, we report the results of our participation in the TREC-COVID challenge. To meet the challenge of building a search engine for rapidly evolving biomedical collection, we propose a simple yet effective weighted hierarchical rank fusion approach, that ensembles together 102 runs from (a) lexical and semantic retrieval systems, (b) pre-trained and fine-tuned BERT rankers, and (c) relevance feedback runs. Our ablation studies demonstrate the contributions of each of these systems to the overall ensemble. The submitted ensemble runs achieved state-of-the-art performance in rounds 4 and 5 of the TREC-COVID challenge.


title: The global symmetry group of first order differential equations and the global rectification theorem

abstract: Symmetry analysis can provide a suitable change of variables, i.e., in geometric terms, a suitable diffeomorphism that simplifies the given direction field, which can help significantly in solving or studying differential equations. Roughly speaking this is the so-called rectification theorem. The local version of this result is a well-known theorem in the field of ordinary differential equations. In this note we prove a global counterpart when the equation fulfils the Lipschitz condition. Then we use this result to determine the global symmetry group of such an ordinary differential equation. It turns out that, assuming the Lipschitz condition, the full symmetry group is a smooth wreath product of two diffeomorphism groups, and does not depend on the form of the equation, at all.


title: Heteroscedastic Bayesian Optimisation for Stochastic Model Predictive Control

abstract: Model predictive control (MPC) has been successful in applications involving the control of complex physical systems. This class of controllers leverages the information provided by an approximate model of the system's dynamics to simulate the effect of control actions. MPC methods also present a few hyper-parameters which may require a relatively expensive tuning process by demanding interactions with the physical system. Therefore, we investigate fine-tuning MPC methods in the context of stochastic MPC, which presents extra challenges due to the randomness of the controller's actions. In these scenarios, performance outcomes present noise, which is not homogeneous across the domain of possible hyper-parameter settings, but which varies in an input-dependent way. To address these issues, we propose a Bayesian optimisation framework that accounts for heteroscedastic noise to tune hyper-parameters in control problems. Empirical results on benchmark continuous control tasks and a physical robot support the proposed framework's suitability relative to baselines, which do not take heteroscedasticity into account.


title: Gamma-ray and conversion-electron spectroscopy of the high-spin isomer in 145Sm

abstract: Background: High-spin isomers at {\approx}9-MeV excitation energies have been reported in several N = 83 isotones near Z = 64. Spin and parity assignments of J{^\pi} = 49/2+ remain tentative for a number of these states in the odd-A nuclei. Purpose: To study the decay of the (49/2 +) high-spin isomer in 145 Sm, make firm spin and parity assignments to the isomer and states populated in its decay, and investigate the structure of the nucleus. Methods: The 145Sm isomer was populated in the 124Sn(26Mg,5n) reaction. Gamma-ray and conversion-electron data were collected using the Solenogam array. Results: A revised lifetime of t1/2 = 3.52(16) mus was measured for the high-spin isomer. Several new states have been added to the level scheme, and a new state at 8815 keV is proposed as the isomer, based on decay-property systematics, transition strengths, and spin and parity assignments. Firm spin and parity assignments have been made to states up to and including the isomer and the new level scheme is interpreted using shell-model calculations performed with the KShell program. Conclusions: The interpretation of the 49/2+ isomer as a deformed excitation of the core neutrons remains unchanged, although there has been a significant revision of the level scheme below the isomer, and hence significant reinterpretations of the lower-lying states.


title: Robust Model-Free Learning and Control without Prior Knowledge

abstract: We present a simple model-free control algorithm that is able to robustly learn and stabilize an unknown discrete-time linear system with full control and state feedback subject to arbitrary bounded disturbance and noise sequences. The controller does not require any prior knowledge of the system dynamics, disturbances, or noise, yet it can guarantee robust stability and provides asymptotic and worst-case bounds on the state and input trajectories. To the best of our knowledge, this is the first model-free algorithm that comes with such robust stability guarantees without the need to make any prior assumptions about the system. We would like to highlight the new convex geometry-based approach taken towards robust stability analysis which served as a key enabler in our results. We will conclude with simulation results that show that despite the generality and simplicity, the controller demonstrates good closed-loop performance.


title: The vacuum boundary problem for the spherically symmetric compressible Euler equations with positive density and unbounded entropy

abstract: Global stability of the spherically symmetric nonisentropic compressible Euler equations with positive density around global-in-time background affine solutions is shown in the presence of free vacuum boundaries. Vacuum is achieved despite a non-vanishing density by considering a negatively unbounded entropy and we use a novel weighted energy method whereby the exponential of the entropy will act as a changing weight to handle the degeneracy of the vacuum boundary. Spherical symmetry introduces a coordinate singularity near the origin for which we adapt a method developed for the Euler-Poisson system by Guo, Had\v{z}i\'c and Jang to our problem.


title: An Anonymous Trust-Marking Scheme on Blockchain Systems

abstract: During the Coincheck incident, which recorded the largest damages in cryptocurrency history in 2018, it was demonstrated that using Mosaic token can have a certain effect. Although it seems attractive to employ tokens as countermeasures for cryptocurrency leakage, Mosaic is a specific token for the New Economy Movement (NEM) cryptocurrency and is not employed for other blockchain systems or cryptocurrencies. Moreover, although some volunteers tracked leaked NEM using Mosaic in the CoinCheck incident, it would be better to verify that the volunteers can be trusted. Simultaneously, if someone (e.g., who stole cryptocurrencies) can identify the volunteers, then that person or organization may be targets of them. In this paper, we propose an anonymous trust-marking scheme on blockchain systems that is universally applicable to any cryptocurrency. In our scheme, entities called token admitters are allowed to generate tokens adding trustworthiness or untrustworthiness to addresses. Anyone can anonymously verify whether these tokens were issued by a token admitter. Simultaneously, only the designated auditor and no one else, including nondesignated auditors, can identify the token admitters. Our scheme is based on accountable ring signatures and commitment, and is implemented on an elliptic curve called Curve25519, and we confirm that both cryptographic tools are efficient. Moreover, we also confirm that our scheme is applicable to Bitcoin, Ethereum, and NEM.


title: Stochastic Optimal Control for Multivariable Dynamical Systems Using Expectation Maximization

abstract: Trajectory optimization is a fundamental stochastic optimal control problem. This paper deals with a trajectory optimization approach for dynamical systems subject to measurement noise that can be fitted into linear time-varying stochastic models. Exact/complete solutions to these kind of control problems have been deemed analytically intractable in literature because they come under the category of Partially Observable Markov Decision Processes (POMDPs). Therefore, effective solutions with reasonable approximations are widely sought for. We propose a reformulation of stochastic control in a reinforcement learning setting. This type of formulation assimilates the benefits of conventional optimal control procedure, with the advantages of maximum likelihood approaches. Finally, an iterative trajectory optimization paradigm called as Stochastic Optimal Control - Expectation Maximization (SOC-EM) is put-forth. This trajectory optimization procedure exhibits better performance in terms of reduction of cumulative cost-to-go which is proved both theoretically and empirically. Furthermore, we also provide novel theoretical work which is related to uniqueness of control parameter estimates. Analysis of the control covariance matrix is presented, which handles stochasticity through efficiently balancing exploration and exploitation.


title: Moment functions on groups

abstract: The main purpose of this work is to prove characterization theorems for generalized moment functions on groups. According one of the main results these are exponential polynomials that can be described with the aid of complete (exponential) Bell polynomials. These characterizations will be immediate consequences of our main result about the characterization of generalized moment functions of higher rank.


title: Spatially resolved study of the Local Group galaxies

abstract: We have been investigating the metallicity dependence of star and planet formation, particularly focusing on the outer Galaxy with a Galactocentric distance ($R_g$) of $\gtrsim$15 kpc, where the metallicity is determined to be as low as $\sim$$-$1 dex. We have obtained near-infrared (NIR) images of young clusters in the outer Galaxy through the 8.2-m Subaru Telescope and have clearly resolved the cluster members with mass detection limits of $\sim$0.1 $M_\odot$. Consequently, we determined that the initial mass function (IMF) in the outer Galaxy is consistent with that in the solar neighborhood with regard to high-mass slope and IMF peak. Meanwhile, we suggested that the lifetime of protoplanetary disks is significantly shorter than that in the solar neighborhood. We also reported a metallicity dependence of the disk lifetime. Future observations with higher spatial resolution and sensitivity by using ELT will allow us to extend the spatially resolved studies on the IMF and protoplanetary disk evolution to Local Group galaxies. With high spatial resolution and sensitivity (i.e. 0.02 arcsec and $K\sim27$ mag with an adaptive optics), stars with mass of down to $\sim$0.1 $M_\odot$ can be detected and also sufficiently resolved. Based on such a study in a wider variety of environments, we can gain new insights related to environmental effects of star and planet formation.


title: Reading the Lattice QCD Form Factors of the $\Lambda_b \rightarrow \Lambda_c$ Transition Using a C-Code

abstract: The $\Lambda_b \rightarrow \Lambda_c \ell^- \bar{\nu}_\ell$ process has been discussed in the literatures as a probe for new physics. This process receives contributions in new physics models from scalar, vector, and tensor hadronic currents. The form factors of these currents can be obtained from the quark model or the lattice QCD. In this work, we present a c-code to read the scalar, vector, and tensor lattice QCD form factors. A Mathematica code was used in Ref. \cite{Datta:2017aue} to read the form factors. Our c-code is much faster than the Mathematica code where the ratio of the wall clock time for the c-code compared to the Mathematica code is 1:64.2 per data point.


title: Tracking and regret bounds for online zeroth-order Euclidean and Riemannian optimisation

abstract: We study numerical optimisation algorithms that use zeroth-order information to minimise time-varying geodesically-convex cost functions on Riemannian manifolds. In the Euclidean setting, zeroth-order algorithms have received a lot of attention in both the time-varying and time-invariant cases. However, the extension to Riemannian manifolds is much less developed. We focus on Hadamard manifolds, which are a special class of Riemannian manifolds with global nonpositive curvature that offer convenient grounds for the generalisation of convexity notions. Specifically, we derive bounds on the expected instantaneous tracking error, and we provide algorithm parameter values that minimise the algorithm's performance. Our results illustrate how the manifold geometry in terms of the sectional curvature affects these bounds. Additionally, we provide dynamic regret bounds for this online optimisation setting. To the best of our knowledge, these are the first regret bounds even for the Euclidean version of the problem. Lastly, via numerical simulations, we demonstrate the applicability of our algorithm on an online Karcher mean problem.


title: How Macroeconomists Lost Control of Stabilization Policy: Towards Dark Ages

abstract: This paper is a study of the history of the transplant of mathematical tools using negative feedback for macroeconomic stabilization policy from 1948 to 1975 and the subsequent break of the use of control for stabilization policy which occurred from 1975 to 1993. New-classical macroeconomists selected a subset of the tools of control that favored their support of rules against discretionary stabilization policy. The Lucas critique and Kydland and Prescott's time-inconsistency were over-statements that led to the "dark ages" of the prevalence of the stabilization-policy-ineffectiveness idea. These over-statements were later revised following the success of the Taylor rule.


title: Matrix product state recursion methods for strongly correlated quantum systems

abstract: We present a method for extrapolation of real-time dynamical correlation functions which can improve the capability of matrix product state methods to compute spectral functions. Unlike the widely used linear prediction method, which ignores the origin of the data being extrapolated, our recursion methods utilize a representation of the wavefunction in terms of an expansion of the same wavefunction and its translations at earlier times. This recursion method is exact for a noninteracting Fermi system. Surprisingly, the recursion method is also more robust than linear prediction at large interaction strength. We test this method on the Hubbard two-leg ladder and present more accurate results for the spectral function than previous studies.


title: ATP-driven separation of liquid phase condensates in bacteria

abstract: Liquid-liquid phase separated (LLPS) states are key to compartmentalise components in the absence of membranes, however it is unclear whether LLPS condensates are actively and specifically organized in the sub-cellular space and by which mechanisms. Here, we address this question by focusing on the ParABS DNA segregation system, composed of a centromeric-like sequence (parS), a DNA-binding protein (ParB) and a motor (ParA). We show that parS-ParB associate to form nanometer-sized, round condensates. ParB molecules diffuse rapidly within the nucleoid volume, but display confined motions when trapped inside ParB condensates. Single ParB molecules are able to rapidly diffuse between different condensates, and nucleation is strongly favoured by parS. Notably, the ParA motor is required to prevent the fusion of ParB condensates. These results describe a novel active mechanism that splits, segregates and localises non-canonical LLPS condensates in the sub-cellular space.