title: The stratified disk wind of MCG-03-58-007

abstract: Past Suzaku, XMM and NuSTAR observations of the nearby (z=0.0323) bright Seyfert 2 galaxy MCG-03-58-007 revealed the presence of two deep and blue-shifted Fe K-shell absorption line profiles. These could be explained with the presence of two phases of a highly ionized, high column density accretion disk wind outflowing with $v_{out1}\sim -0.1c$ and $v_{out2}\sim -0.2c$. Here we present two new observations of MCG-03-58-007: one was carried out in 2016 with Chandra and one in 2018 with Swift. Both caught MCG-03-58-007 in a brighter state ($F_{\mathrm{2-10\,keV}}\sim 4\times 10^{-12}$ erg cm$^{-2}$ s$^{-1}$) confirming the presence of the fast disk wind. The multi-epoch observations of MCG-03-58-007 covering the period from 2010 to 2018 were then analysed. These data show that the lower velocity component outflowing with $v_{out1}\sim -0.072\pm 0.002c$ is persistent and detected in all the observations, although it is variable in column density in the range $N_\rm{H}\sim 3-8 \times 10^{23}$cm$^{-2}$. In the 2016 Swift observation we detected again the second faster component outflowing with $v_{out2}\sim -0.2c$, with a column density ($N_\rm{H}=7.0^{+5.6}_{-4.1}\times 10^{23}$cm$^{-2}$), similar to that seen during the Suzaku observation. However during the Chandra observation two years earlier, this zone was not present ($N_\rm{H}<1.5\times 10^{23}$cm$^{-2}$), suggesting that this faster zone is intermittent. Overall the multi-epochs observations show that the disk wind in MCG-03-58-007 is not only powerful, but also extremely variable, hence placing MCG-03-58-007 among unique disk winds such as the one seen in the famous QSO PDS456. One of the main results of this investigation is the consideration that these winds could be extremely variable, sometime appearing and sometime disappearing; thus to reach solid and firm conclusions about their energetics multiple observations are mandatory.


title: Reconstruction of Voxels with Position- and Angle-Dependent Weightings

abstract: The reconstruction problem of voxels with individual weightings can be modeled a position- and angle- dependent function in the forward-projection. This changes the system matrix and prohibits to use standard filtered backprojection. In this work we first formulate this reconstruction problem in terms of a system matrix and weighting part. We compute the pseudoinverse and show that the solution is rank-deficient and hence very ill posed. This is a fundamental limitation for reconstruction. We then derive an iterative solution and experimentally show its uperiority to any closed-form solution.


title: Motion Compensated Whole-Heart Coronary Magnetic Resonance Angiography using Focused Navigation (fNAV)

abstract: Background: RSN whole-heart CMRA is a technique that estimates and corrects for respiratory motion. However, RSN has been limited to a 1D rigid correction which is often insufficient for patients with complex respiratory patterns. The goal of this work is therefore to improve the robustness and quality of 3D radial CMRA by incorporating both 3D motion information and nonrigid intra-acquisition correction of the data into a framework called focused navigation (fNAV). Methods: We applied fNAV to 500 data sets from a numerical simulation, 22 healthy volunteers, and 549 cardiac patients. We compared fNAV to RSN and respiratory resolved XD-GRASP reconstructions of the same data and recorded reconstruction times. Motion accuracy was measured as the correlation between fNAV and ground truth for simulations, and fNAV and image registration for in vivo data. Vessel sharpness was measured using Soap-Bubble. Finally, image quality analysis was performed by a blinded expert reviewer who chose the best image for each data set. Results The reconstruction time for fNAV images was significantly higher than RSN (6.1 +/- 2.1 minutes vs 1.4 +/- 0.3, minutes, p<0.025) but significantly lower than XD-GRASP (25.6 +/- 7.1, minutes, p<0.025). There is high correlation between the fNAV, and reference displacement estimates across all data sets (0.73 +/- 0.29). For all data, fNAV lead to significantly sharper vessels than all other reconstructions (p < 0.01). Finally, a blinded reviewer chose fNAV as the best image in 239 out of 571 cases (p = 10-5). Conclusion: fNAV is a promising technique for improving free-breathing 3D radial whole-heart CMRA. This novel approach to respiratory self-navigation can derive 3D nonrigid motion estimations from an acquired 1D signal yielding statistically significant improvement in image sharpness relative to 1D translational correction as well as XD-GRASP reconstructions.


title: Predicting Gas-Particle Partitioning Coefficients of Atmospheric Molecules with Machine Learning

abstract: The formation, properties and lifetime of secondary organic aerosols in the atmosphere are largely determined by gas-particle partitioning coefficients of the participating organic vapours. Since these coefficients are often difficult to measure or compute, we developed a machine learning (ML) model to predict them given molecular structure as input. Our data-driven approach is based on the dataset by Wang et al. (Atmos. Chem. Phys., 17, 7529 (2017)), who computed the partitioning coefficients and saturation vapour pressures of 3414 atmospheric oxidation products from the master chemical mechanism using the COSMOtherm program. We train a kernel ridge regression (KRR) ML model on the saturation vapour pressure ($P_{sat}$), and on two equilibrium partitioning coefficients: between a water-insoluble organic matter phase and the gas phase ($K_{WIOM/G}$), and between an infinitely dilute solution with pure water and the gas phase ($K_{W/G}$). For the input representation of the atomic structure of each organic molecule to the machine, we test different descriptors. Our best ML model predicts $P_{sat}$ and $K_{WIOM/G}$ to within 0.3 and $K_{W/G}$ to within 0.4 logarithmic units of the original COSMOtherm calculations. This is equal or better than the typical accuracy of COSMOtherm predictions compared to experimental data. We then apply our ML model to a dataset of 35,383 molecules that we generated based on a carbon 10 backbone and functionalized with 0 to 6 carboxyl, carbonyl or hydroxyl groups to evaluate its performance for polyfunctional compounds with potentially low $P_{sat}$. The resulting $P_{sat}$ and partitioning coefficient distributions were physico-chemically reasonable, and the volatility predictions for the most highly oxidized compounds were in qualitative agreement with experimentally inferred volatilities of atmospheric oxidation products with similar elemental composition.


title: Spiking Neural Networks -- Part I: Detecting Spatial Patterns

abstract: Spiking Neural Networks (SNNs) are biologically inspired machine learning models that build on dynamic neuronal models processing binary and sparse spiking signals in an event-driven, online, fashion. SNNs can be implemented on neuromorphic computing platforms that are emerging as energy-efficient co-processors for learning and inference. This is the first of a series of three papers that introduce SNNs to an audience of engineers by focusing on models, algorithms, and applications. In this first paper, we first cover neural models used for conventional Artificial Neural Networks (ANNs) and SNNs. Then, we review learning algorithms and applications for SNNs that aim at mimicking the functionality of ANNs by detecting or generating spatial patterns in rate-encoded spiking signals. We specifically discuss ANN-to-SNN conversion and neural sampling. Finally, we validate the capabilities of SNNs for detecting and generating spatial patterns through experiments.


title: Containment strategies after the first wave of COVID-19 using mobility data

abstract: In their response to the COVID-19 outbreak, governments face the dilemma to balance public health and economy. Mobility plays a central role in this dilemma because the movement of people enables both economic activity and virus spread. We use mobility data in the form of counts of travelers between regions, to extend the often-used SEIR models to include mobility between regions. We quantify the trade-off between mobility and infection spread in terms of a single parameter, to be chosen by policy makers, and propose strategies for restricting mobility so that the restrictions are minimal while the infection spread is effectively limited. We consider restrictions where the country is divided into regions, and study scenarios where mobility is allowed within these regions, and disallowed between them. We propose heuristic methods to approximate optimal choices for these regions. We evaluate the obtained restrictions based on our trade-off. The results show that our methods are especially effective when the infections are highly concentrated, e.g., around a few municipalities, as resulting from superspreading events that play an important role in the spread of COVID-19. We demonstrate our method in the example of the Netherlands. The results apply more broadly when mobility data is available.


title: On Erd\H{o}s Chains in the Plane

abstract: Let $P$ be a finite point set in $\mathbb{R}^2$ with the set of distance $n$-chains defined as $$ \Delta_n(P)=\{(|p_1-p_2|,|p_2-p_3|,\ldots,|p_n-p_{n+1}|):p_i \in P\}.$$ We show that for $2\leq n=O_{|P|}(1)$ we have $$|\Delta_n(P)|\gtrsim \frac{|P|^{n}}{\log^{\frac{13}{2}(n-1)}|P|}.$$ Our argument uses the energy construction of Elekes and a general version of Rudnev's rich-line bound implicit in Rudnev's recent hinge paper which allows one to iterate efficiently on highly intersecting nested subsets of Guth-Katz lines. Let $G$ is a simple connected graph on $m=O(1)$ vertices with $m\geq 2$. Define the graph-distance set $\Delta_G(P)$ as $$ \Delta_G(P) = \{ (|p_{i}-p_{j}|)_{\{i,j\}\in E(G)} : p_i,p_j \in P\}.$$ Combining with results of Guth and Katz and Rudnev with the above, if $G$ has a Hamiltonian path we have $$ |\Delta_G(P)| \gtrsim \frac{|P|^{m-1}}{\text{polylog}|P|}. $$ \end{abstract}


title: Dataset: LoED: The LoRaWAN at the Edge Dataset

abstract: This paper presents the LoRaWAN at the Edge Dataset (LoED), an open LoRaWAN packet dataset collected at gateways. Real-world LoRaWAN datasets are important for repeatable sensor-network and communications research and evaluation as, if carefully collected, they provide realistic working assumptions. LoED data is collected from nine gateways over a four month period in a dense urban environment. The dataset contains packet header information and all physical layer properties reported by gateways such as the CRC, RSSI, SNR and spreading factor. Files are provided to analyse the data and get aggregated statistics. The dataset is available at: doi.org/10.5281/zenodo.4121430


title: Renovating Requirements Engineering: First Thoughts to Shape Requirements Engineering as a Profession

abstract: Legacy software systems typically include vital data for organizations that use them and should thus to be regularly maintained. Ideally, organizations should rely on Requirements Engineers to understand and manage changes of stakeholder needs and system constraints. However, due to time and cost pressure, and with a heavy focus on implementation, organizations often choose to forgo Requirements Engineers and rather focus on ad-hoc bug fixing and maintenance. This position paper discusses what Requirements Engineers could possibly learn from other similar roles to become crucial for the evolution of legacy systems. Particularly, we compare the roles of Requirements Engineers (according to IREB), Building Architects (according to the German regulations), and Product Owners (according to "The Scrum-Guide"). We discuss overlaps along four dimensions: liability, self-portrayal, core activities, and artifacts. Finally we draw insights from these related fields to foster the concept of a Requirements Engineer as a distinguished profession.


title: How planets grow by pebble accretion. III. Emergence of an interior composition gradient

abstract: During their formation, planets form large, hot atmospheres due to the ongoing accretion of solids. It has been customary to assume that all solids end up at the center constituting a "core" of refractory materials, whereas the envelope remains metal-free. Recent work, as well as observations by the JUNO mission, indicate however that the distinction may not be so clear cut. Indeed, small silicate, pebble-sized particles will sublimate in the atmosphere when they hit the sublimation temperature (T ~ 2,000 K). In this paper we extend previous analytical work to compute the properties of planets under such a pebble accretion scenario. We conduct 1D numerical calculations of the atmosphere of an accreting planet, solving the stellar structure equations, augmented by a non-ideal equation of state that describes a hydrogen/helium-silicate vapor mixture. Calculations terminate at the point where the total mass in metal equals that of the H/He gas, which we numerically confirm as the onset of runaway gas accretion. When pebbles sublimate before reaching the core, insufficient (accretion) energy is available to mix dense, vapor-rich lower layers with the higher layers of lower metallicity. A gradual structure in which Z decreases with radius is therefore a natural outcome of planet formation by pebble accretion. We highlight, furthermore, that (small) pebbles can act as the dominant source of opacity, preventing rapid cooling and presenting a channel for (mini-)Neptunes to survive in gas-rich disks. Nevertheless, once pebble accretion subsides, the atmosphere rapidly clears followed by runaway gas accretion. We consider atmospheric recycling to be the more probable mechanisms that have stalled the growth of these planets' envelopes.


title: Index theorem on $T^2/\mathbb{Z}_N$ orbifolds

abstract: We investigate chiral zero modes and winding numbers at fixed points on $T^2/\mathbb{Z}_N$ orbifolds. It is shown that the Atiyah-Singer index theorem for the chiral zero modes leads to a formula $n_+-n_-=(-V_++V_-)/2N$, where $n_{\pm}$ are the numbers of the $\pm$ chiral zero modes and $V_{\pm}$ are the sums of the winding numbers at the fixed points on $T^2/\mathbb{Z}_N$. This formula is complementary to our zero-mode counting formula on the magnetized orbifolds with non-zero flux background $M \neq 0$, consistently with substituting $M = 0$ for the counting formula $n_+ - n_- = (2M - V_+ + V_-)/2N$.


title: Mpemba-like effect in driven binary mixtures

abstract: The Mpemba effect occurs when two samples at different initial temperatures evolve in such a way that the temperatures cross each other during the relaxation towards equilibrium. In this paper we show the emergence of a Mpemba-like effect in a molecular binary mixture in contact with a thermal reservoir (bath). The interaction between the gaseous particles of the mixture and the thermal reservoir is modeled via a viscous drag force plus a stochastic Langevin-like term. The presence of the external bath couples the time evolution of the total and partial temperatures of each component allowing the appearance of the Mpemba phenomenon, even when the initial temperature differences are of the same order of the temperatures themselves. Analytical results are obtained by considering multitemperature Maxwellian approximations for the velocity distribution functions of each component. The theoretical analysis is carried out for initial states close to and far away (large Mpemba-like effect) from equilibrium. The former situation allows us to develop a simple theory where the time evolution equation for the temperature is linearized around its asymptotic equilibrium solution. This linear theory provides an expression for the crossover time. We also provide a qualitative description of the large Mpemba effect. Our theoretical results agree very well with computer simulations obtained by numerically solving the Enskog kinetic equation by means of the direct simulation Monte Carlo method and by performing molecular dynamics simulations. Finally, preliminary results for driven granular mixtures also show the occurrence of a Mpemba-like effect for inelastic collisions.


title: Coupling Light to Higher Order Transverse Modes of a Near-Concentric Optical Cavity

abstract: Optical cavities in the near-concentric regime have near-degenerate transverse modes; the tight focusing transverse modes in this regime enable strong coupling with atoms. These features provide an interesting platform to explore multi-mode interaction between atoms and light. Here, we use a spatial light modulator (SLM) to shape the phase of an incoming light beam to match several Laguerre-Gaussian (LG) modes of a near-concentric optical cavity. We demonstrate coupling efficiency close to the theoretical prediction for single LG modes and well-defined combinations of them, limited mainly by imperfections in the cavity alignment.


title: Spiking Neural Networks -- Part II: Detecting Spatio-Temporal Patterns

abstract: Inspired by the operation of biological brains, Spiking Neural Networks (SNNs) have the unique ability to detect information encoded in spatio-temporal patterns of spiking signals. Examples of data types requiring spatio-temporal processing include logs of time stamps, e.g., of tweets, and outputs of neural prostheses and neuromorphic sensors. In this paper, the second of a series of three review papers on SNNs, we first review models and training algorithms for the dominant approach that considers SNNs as a Recurrent Neural Network (RNN) and adapt learning rules based on backpropagation through time to the requirements of SNNs. In order to tackle the non-differentiability of the spiking mechanism, state-of-the-art solutions use surrogate gradients that approximate the threshold activation function with a differentiable function. Then, we describe an alternative approach that relies on probabilistic models for spiking neurons, allowing the derivation of local learning rules via stochastic estimates of the gradient. Finally, experiments are provided for neuromorphic data sets, yielding insights on accuracy and convergence under different SNN models.


title: Search for electric dipole moment

abstract: The outstanding progress has been made in reducing the upper bounds on EDM of several particles. Even if significant challenges must be overcome to further improve these limits, it is still one of the best chances to detect new type of interactions beyond the standard model. Analyzing several examples, we highlight a common thread that is visible in different set-ups used for the EDM detection. The electric dipole moment is one of the clear consequences of CP- or T-violating interactions, however it is not the only one. These symmetry-violating interactions enable extra phenomena that unavoidably accompany the EDM-induced spin precession, and they must be taken into account in planning and executing sensitive experiments. After reviewing three typical cases, we suggest conditions for improving the sensitivity of detecting the intrinsic EDM.


title: Freezing of a soft-core fluid in a one-dimensional potential: Appearance of a locked smectic phase

abstract: We investigate the phase behaviour of a two-dimensional colloidal model system of ultra-soft particles on a substrate which varies periodically along one spatial direction. Our calculations are based on mean-field density functional theory for a system of particles interacting via an ultra-soft potential, that is, the generalized exponential model with index four (\mbox{GEM-4}). For suitable substrate periodicities (with commensurability parameter $p=2$), we find a succession of phase transitions from a modulated liquid to a locked smectic and then to a locked floating solid phase. The appearance of a locked smectic phase is consistent with earlier theoretical predictions and experiments for freezing of more repulsive systems on structured surfaces (with $p=2$). However, the present ultra-soft system does not display re-entrant melting. We here investigate the details of the density distributions of the different phases, thereby supplementing earlier work on GEM-4 systems with $p=1$ [Phys. Rev. E \textbf{101}, 012609 (2020)]. Interestingly, the observed succession of phase transitions can be triggered through different paths along which physical control parameters are changed.


title: Spiking Neural Networks -- Part III: Neuromorphic Communications

abstract: Synergies between wireless communications and artificial intelligence are increasingly motivating research at the intersection of the two fields. On the one hand, the presence of more and more wirelessly connected devices, each with its own data, is driving efforts to export advances in machine learning (ML) from high performance computing facilities, where information is stored and processed in a single location, to distributed, privacy-minded, processing at the end user. On the other hand, ML can address algorithm and model deficits in the optimization of communication protocols. However, implementing ML models for learning and inference on battery-powered devices that are connected via bandwidth-constrained channels remains challenging. This paper explores two ways in which Spiking Neural Networks (SNNs) can help address these open problems. First, we discuss federated learning for the distributed training of SNNs, and then describe the integration of neuromorphic sensing, SNNs, and impulse radio technologies for low-power remote inference.


title: A Derived Lagrangian Fibration on the Derived Critical Locus

abstract: We study the symplectic geometry of derived intersections of Lagrangian morphisms. In particular, we show that for a functional $f : X \rightarrow \mathbb{A}_k^1$, the derived critical locus has a natural Lagrangian fibration $\textbf{Crit}(f) \rightarrow X$. In the case where $f$ is non-degenerate and the strict critical locus is smooth, we show that the Lagrangian fibration on the derived critical locus is determined by the Hessian quadratic form.


title: Roton instabilities in the superfluid outer core of neutron stars

abstract: We study the superfluid dynamics of the outer core of neutron stars by means of a hydrodynamic model made of a neutronic superfluid and a protonic superconductor, coupled by both the dynamic entrainment and the Skyrme SLy4 nucleon-nucleon interactions. The resulting nonlinear equations of motion are probed in the search for dynamical instabilities triggered by the relative motion of the superfluids that could be related to observed timing anomalies in pulsars. Through linear analysis, the origin and expected growth of the instabilities is explored for varying nuclear-matter density. Differently from previous findings, the dispersion of linear excitations in our model shows rotonic structures below the pair-breaking energy threshold, which lies at the origin of the dynamical instabilities, and could eventually lead to emergent vorticity along with modulations of the superfluid density.


title: UAV-Assisted and Intelligent Reflecting Surfaces-Supported Terahertz Communications

abstract: In this paper, unmanned aerial vehicles (UAVs) and intelligent reflective surface (IRS) are utilized to support terahertz (THz) communications. To this end, the joint optimization of UAV's trajectory, the phase shift of IRS, the allocation of THz sub-bands, and the power control is investigated to maximize the minimum average achievable rate of all the users. An iteration algorithm based on successive Convex Approximation with the Rate constraint penalty (CAR) is developed to obtain UAV's trajectory, and the IRS phase shift is formulated as a closed-form expression with introduced pricing factors. Simulation results show that the proposed scheme significantly enhances the rate performance of the whole system.


title: Distortion Representations of Multivariate Distributions

abstract: The univariate distorted distribution were introduced in risk theory to represent changes (distortions) in the expected distributions of some risks. Later they were also applied to represent distributions of order statistics, coherent systems, proportional hazard rate (PHR) and proportional reversed hazard rate (PRHR) models, etc. In this paper we extend this concept to the multivariate setup. We show that, in some cases, they are a valid alternative to the copula representations especially when the marginal distributions may not be easily handled. Several relevant examples illustrate the applications of such representations in statistical modeling. They include the study of paired (dependent) ordered data, joint residual lifetimes, order statistics and coherent systems.


title: Fully developed relativistic turbulence

abstract: We use a simple model consisting of energy-momentum tensor conservation and a Maxwell-Cattaneo equation for its viscous part to study nonlinear phenomena in a real relativistic fluid. We focus on new types of behavior without nonrelativistic equivalents, such as an entropy cascade driven by fluctuations in the tensor degrees of freedom of the theory. We write down the von K\'arm\'an-Howarth equations for this kind of turbulence, and consider the correlations corresponding to fully developed turbulence.


title: Range-Net: A High Precision Streaming SVD for Big Data Applications

abstract: In a Big Data setting computing the dominant SVD factors is restrictive due to the main memory requirements. Recently introduced streaming Randomized SVD schemes work under the restrictive assumption that the singular value spectrum of the data has exponential decay. This is seldom true for any practical data. Although these methods are claimed to be applicable to scientific computations due to associated tail-energy error bounds, the approximation errors in the singular vectors and values are high when the aforementioned assumption does not hold. Furthermore from a practical perspective, oversampling can still be memory intensive or worse can exceed the feature dimension of the data. To address these issues, we present Range-Net as an alternative to randomized SVD that satisfies the tail-energy lower bound given by Eckart-Young-Mirsky (EYM) theorem. Range-Net is a deterministic two-stage neural optimization approach with random initialization, where the main memory requirement depends explicitly on the feature dimension and desired rank, independent of the sample dimension. The data samples are read in a streaming setting with the network minimization problem converging to the desired rank-r approximation. Range-Net is fully interpretable where all the network outputs and weights have a specific meaning. We provide theoretical guarantees that Range-Net extracted SVD factors satisfy EYM tail-energy lower bound at machine precision. Our numerical experiments on real data at various scales confirms this bound. A comparison against the state of the art streaming Randomized SVD shows that Range-Net accuracy is better by six orders of magnitude while being memory efficient.


title: Efficient, Simple and Automated Negative Sampling for Knowledge Graph Embedding

abstract: Negative sampling, which samples negative triplets from non-observed ones in knowledge graph (KG), is an essential step in KG embedding. Recently, generative adversarial network (GAN), has been introduced in negative sampling. By sampling negative triplets with large gradients, these methods avoid the problem of vanishing gradient and thus obtain better performance. However, they make the original model more complex and harder to train. In this paper, motivated by the observation that negative triplets with large gradients are important but rare, we propose to directly keep track of them with the cache. In this way, our method acts as a "distilled" version of previous GAN-based methods, which does not waste training time on additional parameters to fit the full distribution of negative triplets. However, how to sample from and update the cache are two critical questions. We propose to solve these issues by automated machine learning techniques. The automated version also covers GAN-based methods as special cases. Theoretical explanation of NSCaching is also provided, justifying the superior over fixed sampling scheme. Besides, we further extend NSCaching with skip-gram model for graph embedding. Finally, extensive experiments show that our method can gain significant improvements on various KG embedding models and the skip-gram model, and outperforms the state-of-the-art negative sampling methods.


title: New interfaces for musical expression

abstract: The rapid evolution of electronics, digital media, advanced materials, and other areas of technology, is opening up unprecedented opportunities for musical interface inventors and designers. The possibilities afforded by these new technologies carry with them the challenges of a complex and often confusing array of choices for musical composers and performers. New musical technologies are at least partly responsible for the current explosion of new musical forms, some of which are controversial and challenge traditional definitions of music. Alternative musical controllers, currently the leading edge of the ongoing dialogue between technology and musical culture, involve many of the issues covered at past CHI meetings. This workshop brings together interface experts interested in musical controllers and musicians and composers involved in the development of new musical interfaces.


title: Infection spreading and recovery in a square lattice

abstract: We investigate spreading and recovery of disease in a square lattice, and in particular, emphasize the role of the initial distribution of infected patches in the network, on the progression of an endemic and initiation of a recovery process, if any, due to migration of both the susceptible and infected hosts. The disease starts in the lattice with three possible initial distribution patterns of infected and infection-free sites, infected core patches (ICP), infected peripheral patches (IPP) and randomly distributed infected patches (RDIP). Our results show that infection spreads monotonically in the lattice with increasing migration without showing any sign of recovery in the ICP case. In the IPP case, it follows a similar monotonic progression with increasing migration, however, a self-organized healing process starts for higher migration, leading the lattice to full recovery at a critical rate of migration. Encouragingly, for the initial RDIP arrangement, chances of recovery are much higher with a lower rate of critical migration. An eigenvalue based semi-analytical study is made to determine the critical migration rate for realizing a stable infection-free lattice. The initial fraction of infected patches and the force of infection play significant roles in the self-organized recovery. They follow an exponential law, for the RDIP case, that governs the recovery process. For the frustrating case of ICP arrangement, we propose a random rewiring of links in the lattice allowing long-distance migratory paths that effectively initiate a recovery process. Global prevalence of infection thereby declines and progressively improves with the rewiring probability that follows a power law with the critical migration and leads to the birth of emergent infection-free networks.


title: A Comparison of Discrete Latent Variable Models for Speech Representation Learning

abstract: Neural latent variable models enable the discovery of interesting structure in speech audio data. This paper presents a comparison of two different approaches which are broadly based on predicting future time-steps or auto-encoding the input signal. Our study compares the representations learned by vq-vae and vq-wav2vec in terms of sub-word unit discovery and phoneme recognition performance. Results show that future time-step prediction with vq-wav2vec achieves better performance. The best system achieves an error rate of 13.22 on the ZeroSpeech 2019 ABX phoneme discrimination challenge


title: Virtual Alignment Method and its application to the dental prostheses and diagnosis

abstract: The recent proposal of a new alignment solution for X-ray tomography, Virtual alignment method (VAM) allowed a more accurate method to remove the possible errors that limit the resolution and clarity of the reconstructed image. In the field of dentistry, the movement of patients during the scanning poses as one of the major factors hindering the final reconstructed image quality. Here, the patient's movement was artificially given to the projection image set and the newly proposed algorithm using the sinogram and the fixed point was applied to the tooth sample to compare the reconstruction image to the actual projection image set. The new alignment method showed promising results by reducing the margin of errors down to a few micrometer, which will allow the production of high-quality dental prostheses with accuracy and precision. We hope that the newly proposed alignment method can be further investigated to be applied more readily in the filed of dentistry ot provide better quality images of patients to make a more accurate diagnosis and prostheses.


title: Another estimating the absolute value of Mertens function

abstract: Through an inversion approach, we suggest a possible estimation for the absolute value of Mertens function $\vert M(x) \vert$ that $ \left\vert M(x) \right\vert \sim \left[\frac{1}{\pi \sqrt{\varepsilon}(x+\varepsilon)}\right]\sqrt{x}$ (where $x$ is an appropriately large real number, and $\varepsilon$ ($0<\varepsilon<1$) is a small real number which makes $2x+\varepsilon$ to be an integer). For any large $x$, we can always find an $\varepsilon$, so that $\vert M(x) \vert < \left[\frac{1}{\pi \sqrt{\varepsilon}(x+\varepsilon)}\right]\sqrt{x}$.


title: Align-Refine: Non-Autoregressive Speech Recognition via Iterative Realignment

abstract: Non-autoregressive models greatly improve decoding speed over typical sequence-to-sequence models, but suffer from degraded performance. Infilling and iterative refinement models make up some of this gap by editing the outputs of a non-autoregressive model, but are constrained in the edits that they can make. We propose iterative realignment, where refinements occur over latent alignments rather than output sequence space. We demonstrate this in speech recognition with Align-Refine, an end-to-end Transformer-based model which refines connectionist temporal classification (CTC) alignments to allow length-changing insertions and deletions. Align-Refine outperforms Imputer and Mask-CTC, matching an autoregressive baseline on WSJ at 1/14th the real-time factor and attaining a LibriSpeech test-other WER of 9.0% without an LM. Our model is strong even in one iteration with a shallower decoder.


title: Global Sentiment Analysis Of COVID-19 Tweets Over Time

abstract: The Coronavirus pandemic has affected the normal course of life. People around the world have taken to social media to express their opinions and general emotions regarding this phenomenon that has taken over the world by storm. The social networking site, Twitter showed an unprecedented increase in tweets related to the novel Coronavirus in a very short span of time. This paper presents the global sentiment analysis of tweets related to Coronavirus and how the sentiment of people in different countries has changed over time. Furthermore, to determine the impact of Coronavirus on daily aspects of life, tweets related to Work From Home (WFH) and Online Learning were scraped and the change in sentiment over time was observed. In addition, various Machine Learning models such as Long Short Term Memory (LSTM) and Artificial Neural Networks (ANN) were implemented for sentiment classification and their accuracies were determined. Exploratory data analysis was also performed for a dataset providing information about the number of confirmed cases on a per-day basis in a few of the worst-hit countries to provide a comparison between the change in sentiment with the change in cases since the start of this pandemic till June 2020.


title: Multi-XScience: A Large-scale Dataset for Extreme Multi-document Summarization of Scientific Articles

abstract: Multi-document summarization is a challenging task for which there exists little large-scale datasets. We propose Multi-XScience, a large-scale multi-document summarization dataset created from scientific articles. Multi-XScience introduces a challenging multi-document summarization task: writing the related-work section of a paper based on its abstract and the articles it references. Our work is inspired by extreme summarization, a dataset construction protocol that favours abstractive modeling approaches. Descriptive statistics and empirical results---using several state-of-the-art models trained on the Multi-XScience dataset---reveal that Multi-XScience is well suited for abstractive models.


title: Scientific intuition inspired by machine learning generated hypotheses

abstract: Machine learning with application to questions in the physical sciences has become a widely used tool, successfully applied to classification, regression and optimization tasks in many areas. Research focus mostly lies in improving the accuracy of the machine learning models in numerical predictions, while scientific understanding is still almost exclusively generated by human researchers analysing numerical results and drawing conclusions. In this work, we shift the focus on the insights and the knowledge obtained by the machine learning models themselves. In particular, we study how it can be extracted and used to inspire human scientists to increase their intuitions and understanding of natural systems. We apply gradient boosting in decision trees to extract human interpretable insights from big data sets from chemistry and physics. In chemistry, we not only rediscover widely know rules of thumb but also find new interesting motifs that tell us how to control solubility and energy levels of organic molecules. At the same time, in quantum physics, we gain new understanding on experiments for quantum entanglement. The ability to go beyond numerics and to enter the realm of scientific insight and hypothesis generation opens the door to use machine learning to accelerate the discovery of conceptual understanding in some of the most challenging domains of science.

categories: cs.LG cs.AI cs.CE physics.chem-ph quant-ph


title: Space-dependent diffusion with stochastic resetting: A first-passage study

abstract: We explore the effect of stochastic resetting on the first-passage properties of space-dependent diffusion in presence of a constant bias. In our analytically tractable model system, a particle diffusing in a linear potential $U(x)\propto\mu |x|$ with a spatially varying diffusion coefficient $D(x)=D_0|x|$ undergoes stochastic resetting, i.e., returns to its initial position $x_0$ at random intervals of time, with a constant rate $r$. Considering an absorbing boundary placed at $x_a<x_0$, we first derive an exact expression of the survival probability of the diffusing particle in the Laplace space and then explore its first-passage to the origin as a limiting case of that general result. In the limit $x_a\to0$, we derive an exact analytic expression for the first-passage time distribution of the underlying process. Once resetting is introduced, the system is observed to exhibit a series of dynamical transitions in terms of a sole parameter, $\nu=(1+\mu D_0^{-1})$, that captures the interplay of the drift and the diffusion. Constructing a full phase diagram in terms of $\nu$, we show that for $\nu<0$, i.e., when the potential is strongly repulsive, the particle can never reach the origin. In contrast, for weakly repulsive or attractive potential ($\nu>0$), it eventually reaches the origin. Resetting accelerates such first-passage when $\nu<3$, but hinders its completion for $\nu>3$. A resetting transition is therefore observed at $\nu=3$, and we provide a comprehensive analysis of the same. The present study paves the way for an array of theoretical and experimental works that combine stochastic resetting with inhomogeneous diffusion in a conservative force-field.


title: An inverse spectral problem for second-order functional-differential pencils with two delays

abstract: We consider a second order functional-differential pencil with two constant delays of the argument and study the inverse problem of recovering its coefficients from the spectra of two boundary value problems with one common boundary condition. The uniqueness theorem is proved and a constructive procedure for solving this inverse problem along with necessary and sufficient conditions for its solvability is obtained. Moreover, we give a survey on the contemporary state of the inverse spectral theory for operators with delay. The pencil under consideration generalizes Sturm-Liouville-type operators with delay, which allows us to illustrate essential results in this direction, including recently solved open questions.


title: Structural and electronic properties of realistic two-dimensional amorphous topological insulators

abstract: We investigate the structure and electronic spectra properties of two-dimensional amorphous bismuthene structures and show that these systems are topological insulators. We employ realistic modeling of amorphous geometries together with density functional theory for electronic structure calculations. We investigate the system topological properties throughout the amorphization process and find that the robustness of the topological phase is associated with the spin-orbit coupling strength and size of the pristine topological gap. Using recursive non-equilibrium Green's function, we study the electronic transport properties of nanoribbons devices with lengths comparable to experimentally synthesized materials. We find a $2e^2/h$ conductance plateau within the topological gap and an onset of Anderson localization at the trivial insulator phase.


title: GPS Measurements of Precipitable Water Vapor Can Improve Survey Calibration: A Demonstration from KPNO and the Mayall z-band Legacy Survey

abstract: We here show that dual-band GPS measurements of precipitable water vapor (PWV) at KPNO predict the overall per-image sensitivity of the Mayall z-band Legacy Survey (MzLS). The per-image variation in the brightness of individual stars is strongly correlated with the measured PWV and the color of the star. We use synthetic stellar spectra and TAPAS transmission models to predict the expected PWV-induced photometric errors and find good agreement with the observations. We also find that PWV absorption can be well-approximated by a linear relationship with PWV_eff and present an update on the traditional treatment in the literature. Within the range of reasonable observing conditions, the MzLS zero point varies with a standard deviation of 127 mmag. This variation is dominated by a gray secular trend with time, consistent with a gradual accumulation of contamination on optical surfaces that accounts for ~114 mmag of variation. Correcting for PWV based on a suite of stellar spectra and detailed PWV absorption models accounts for another 47 mmag of zero-point variation. The MzLS per-image sensitivity is decreased by ~40 mmag per effective mm of PWV. The difference between blue (r-z < 0.5 mag) and red (1.2 mag < r-z) stars increases by 3.25 mmag per effective mm of PWV. These results show the need for high-precision photometric surveys to simultaneously monitor PWV. We find that this GPS system provides more precise PWV measurements than using differential measurements of stars of different colors and recommend that observatories install dual-band GPS as a low-maintenance, relatively low cost, auxiliary calibration system. We extend our results of the need for well-calibrated PWV measurements by presenting calculations of the PWV photometric impact on three science cases of interest: stellar photometry, supernova cosmology, and quasar identification and variability.


title: Discrete Systems in Thermal Physics and Engineering -- A Glance from Non-Equilibrium Thermodynamics

abstract: Non-equilibrium processes in Schottky systems generate by projection onto the equilibrium subspace reversible accompanying processes for which the non-equilibrium variables are functions of the equilibrium ones. The embedding theorem which guarantees the compatibility of the accompanying processes with the non-equilibrium entropy is proved. The non-equilibrium entropy is defined as a state function on the non-equilibrium state space containing the contact temperature as a non-equilibrium variable. If the entropy production does not depend on the internal energy, the contact temperature changes into the thermostatic temperature also in non-equilibrium, a fact which allows to use temperature as a primitive concept in non-equilibrium. The dissipation inequality is revisited, and an efficiency of generalized cyclic processes beyond the Carnot process is achieved.


title: Deep generative factorization for speech signal

abstract: Various information factors are blended in speech signals, which forms the primary difficulty for most speech information processing tasks. An intuitive idea is to factorize speech signal into individual information factors (e.g., phonetic content and speaker trait), though it turns out to be highly challenging. This paper presents a speech factorization approach based on a novel factorial discriminative normalization flow model (factorial DNF). Experiments conducted on a two-factor case that involves phonetic content and speaker trait demonstrates that the proposed factorial DNF has powerful capability to factorize speech signals and outperforms several comparative models in terms of information representation and manipulation.


title: Squeezing value of cross-domain labels: a decoupled scoring approach for speaker verification

abstract: Domain mismatch often occurs in real applications and causes serious performance reduction on speaker verification systems. The common wisdom is to collect cross-domain data and train a multi-domain PLDA model, with the hope to learn a domain-independent speaker subspace. In this paper, we firstly present an empirical study to show that simply adding cross-domain data does not help performance in conditions with enrollment-test mismatch. Careful analysis shows that this striking result is caused by the incoherent statistics between the enrollment and test conditions. Based on this analysis, we present a decoupled scoring approach that can maximally squeeze the value of cross-domain labels and obtain optimal verification scores when the enrollment and test are mismatched. When the statistics are coherent, the new formulation falls back to the conventional PLDA. Experimental results on cross-channel test show that the proposed approach is highly effective and is a principle solution to domain mismatch.


title: GMACO-P: GPU assisted Preemptive MACO algorithm for enabling Smart Transportation

abstract: Vehicular Ad-hoc NETworks (VANETs) are developing at a very fast pace to enable smart transportation in urban cities, by designing some mechanisms for decreasing travel time for commuters by reducing congestion. Inefficient Traffic signals and routing mechanisms are the major factors that contribute to the increase of road congestion. For smoother traffic movement and reducing congestion on the roads, the waiting time at intersections must be reduced and an optimal path should be chosen simultaneously. In this paper, A GPU assisted Preemptive MACO (GMACO-P) algorithm has been proposed to minimize the total travel time of the commuters. GMACO-P is an improvement of MACO-P algorithm that uses the harnessing the power of the GPU to provide faster computations for further minimizing the travel time. The MACO-P algorithm is based on an existing MACO algorithm that avoid the path with the congestion. The MACO-P algorithm reduces the average queue length at intersections by incorporating preemption that ensures less waiting time. In this paper, GMACO-P algorithm is proposed harnessing the power of GPU to improve MACO-P to further reduce the travel time. The GMACO-P algorithm is executed with CUDA toolkit 7.5 using C language and the obtained results were compared with existing Dijkstra, ACO, MACO, MACO-P, parallel implementation of the Dijkstra, ACO and MACO algorithms. Obtained results show the significant reduction in the travel time after using the proposed GMACO-P algorithm.


title: Apps Against the Spread: Privacy Implications and User Acceptance of COVID-19-Related Smartphone Apps on Three Continents

abstract: The COVID-19 pandemic has fueled the development of smartphone applications to assist disease management. Many "corona apps" require widespread adoption to be effective, which has sparked public debates about the privacy, security, and societal implications of government-backed health applications. We conducted a representative online study in Germany (n = 1,003), the US (n = 1,003), and China (n = 1,019) to investigate user acceptance of corona apps, using a vignette design based on the contextual integrity framework. We explored apps for contact tracing, symptom checks, quarantine enforcement, health certificates, and mere information. Our results provide insights into data processing practices that foster adoption and reveal significant differences between countries, with user acceptance being highest in China and lowest in the US. Chinese participants prefer the collection of personalized data, while German and US participants favor anonymity. Across countries, contact tracing is viewed more positively than quarantine enforcement, and technical malfunctions negatively impact user acceptance.


title: $\mu$NAS: Constrained Neural Architecture Search for Microcontrollers

abstract: IoT devices are powered by microcontroller units (MCUs) which are extremely resource-scarce: a typical MCU may have an underpowered processor and around 64 KB of memory and persistent storage, which is orders of magnitude fewer computational resources than is typically required for deep learning. Designing neural networks for such a platform requires an intricate balance between keeping high predictive performance (accuracy) while achieving low memory and storage usage and inference latency. This is extremely challenging to achieve manually, so in this work, we build a neural architecture search (NAS) system, called $\mu$NAS, to automate the design of such small-yet-powerful MCU-level networks. $\mu$NAS explicitly targets the three primary aspects of resource scarcity of MCUs: the size of RAM, persistent storage and processor speed. $\mu$NAS represents a significant advance in resource-efficient models, especially for "mid-tier" MCUs with memory requirements ranging from 0.5 KB to 64 KB. We show that on a variety of image classification datasets $\mu$NAS is able to (a) improve top-1 classification accuracy by up to 4.8%, or (b) reduce memory footprint by 4--13x, or (c) reduce the number of multiply-accumulate operations by at least 2x, compared to existing MCU specialist literature and resource-efficient models.


title: A simple but efficient concept of blended teaching of mathematics for engineering students during the COVID-19 pandemic

abstract: We present a simple but efficient concept for the realization of blended teaching of mathematics and its applications in theoretical mechanics that was conceived, tested and implemented at the University of Trento, Italy, during the COVID-19 pandemic. The concept foresees traditional blackboard lectures with a reduced number of students present in the lecture hall, while the same lectures are simultaneously made available to the remaining students via high quality low-bandwidth online streaming. Based on our first assumption that traditional blackboard lectures, including the gestures and the facial expressions of the professor, are still a very efficient and highly appreciated means of teaching mathematics, this paper deliberately does not want to propose a novel pedagogical concept of how to teach mathematics, but rather presents a technical concept how to preserve the quality of traditional blackboard lectures even during the pandemic and how to make them available to the students at home via online streaming with adequate audio and video quality at low internet bandwidth. The second assumption is that the teaching of mathematics is a dynamic creative process that requires the physical presence of students in the lecture hall as audience so that the professor can instantaneously fine-tune the evolution of the lecture according to his/her perception of the level of attention and the facial expressions of the students. The third assumption of this paper is that students need to have the possibility to interact with each other personally. We report on the necessary hardware, software and logistics, and on the perception of the proposed blended lectures by students from civil and environmental engineering at the University of Trento, compared to traditional lectures and also compared to the pure online lectures that were needed as emergency measure at the beginning of the pandemic.


title: Human towers or castells modelling. 158th European Study Groups with Industry (long report)

abstract: Human towers or castells are human structures played in festivals mainly in Catalonia. These unique cultural and traditional displays have become very popular in the last years, but they date from the XVIII century. On 2010 they became part of the Unesco Representative List of the Intangible Cultural Heritage of Humanity. Safety is very important in the performance of castells. To this end, it is crucial to understand the mechanisms that allow a castell to be built and, more importantly, the factors that may cause its collapse. This work is focused on the mechanical aspects that make a pilar (the simplest structure in the castells) to stand. We suggest three different but complementary approaches for the running stage of a pilar (stage where it has been built and has not yet collapsed): the $N$-link pendulum as a first dynamical model, the response of the castellers as a control problem, and a static analysis to capture the feasibility of a given configuration. We include some preliminary simulations to better understand the previous approaches, which seem to match with qualitative perceptions that castellers have of a castell. Possible future developments are also discussed. To our knowledge, this work represents the first one to study the castells from a mechanical point of view.


title: Eulers Graph World - Purity, Regularity and Evenness -Law of Nature? - Constructions and Examples

abstract: We propose a Law of Nature? Viz., Pure Regularity Occurs at Na\"ive Levels and Regularity has Affinity with Evenness. In a series of three papers, it was established that regular Euler graphs with only one type of (pure) cycles are nonexistent; Regular Euler graphs with only two types of cycles are possible in one of the six cases, viz., regular bipartite Euler graphs of degree >2; Evenness plays role in unveiling regularity; Lastly, K5 is a regular Euler graph with three types of cycles (0,1,3); This is the only known graph with the property; It is conjectured that regular Euler graphs of order >5 with only three cycle types are nonexistent and this is proved true in part cases in each of the four cases. Some constructions and examples are given for the Euler graphs under (mod 4) satisfying intersection (combined cycle) rules. These infinite classes of Euler graphs serve as candidates for gracefulness. Infinite families of graceful graphs are presented in Case-0.


title: Geometry of Arrangements that Determine Shapes

abstract: Shape grammars compute over shapes which are defined in the universe $U^*$. Shapes in the universe $U^*$ are analogous to line drawings that can be physically realized in the plane. Any shape is embedded or contained in an arrangement of points and lines in the plane called, respectively, registration marks and construction lines, that satisfy special incidence laws. In this paper, arrangements that contain shapes are studied as incidence structures and the finite geometries they give rise to as such are characterized. In particular, arrangements that contain shapes are distinguished into those that give rise to finite near-linear and linear spaces, and those that do not give rise to any proper form of geometry (in the strict mathematical sense). Arrangements that constitute finite geometries (near-linear and linear spaces) give an alternative characterization of determinate rules in shape grammars.


title: Quadratic Formula: Revisiting a Proof Through the Lens of Transformations

abstract: In this paper, we derive the quadratic formula as a consequence of constructively proving the existence of standard and factored forms for general form real quadratic functions. Emphasis is put on connections to graphing of corresponding parabolas through function transformation perspectives of scaling and translation.


title: Neural Networked Assisted Tree Search for the Personnel Rostering Problem

abstract: The personnel rostering problem is the problem of finding an optimal way to assign employees to shifts, subject to a set of hard constraints which all valid solutions must follow, and a set of soft constraints which define the relative quality of valid solutions. The problem has received significant attention in the literature and is addressed by a large number of exact and metaheuristic methods. In order to make the complex and costly design of heuristics for the personnel rostering problem automatic, we propose a new method combined Deep Neural Network and Tree Search. By treating schedules as matrices, the neural network can predict the distance between the current solution and the optimal solution. It can select solution strategies by analyzing existing (near-)optimal solutions to personnel rostering problem instances. Combined with branch and bound, the network can give every node a probability which indicates the distance between it and the optimal one, so that a well-informed choice can be made on which branch to choose next and to prune the search tree.


title: Randomized double and triple Kaczmarz for solving extended normal equations

abstract: The randomized Kaczmarz algorithm has received considerable attention recently because of its simplicity, speed, and the ability to approximately solve large-scale linear systems of equations. In this paper we propose randomized double and triple Kaczmarz algorithms to solve extended normal equations of the form $\bf A^\top Ax=A^\top b-c$. The proposed algorithms avoid forming $\bf A^\top A$ explicitly and work for {\it arbitrary} $\mbf A\in\mbbr^{m\times n}$ (full rank or rank deficient, $m\geq n$ or $m<n$). {\it Tight} upper bounds showing exponential convergence in the mean square sense of the proposed algorithms are presented and numerical experiments are given to illustrate the theoretical results.


title: On (non-)monotonicity and phase diagram of finitary random interlacement

abstract: In this paper, we study the evolution of a Finitary Random Interlacement (FRI) with respect to the expected length of each fiber. In contrast to the previously proved phase transition between sufficiently large and small fiber length, we show that for $d=3,4$, FRI is NOT stochastically monotone as fiber length increasing. At the same time, numerical evidences still strongly support the existence of a unique and sharp phase transition on the existence of a unique infinite cluster, while the critical value for phase transition is estimated to be an inversely proportional function with respect to the system intensity.


title: RH-Net: Improving Neural Relation Extraction via Reinforcement Learning and Hierarchical Relational Searching

abstract: Distant supervision (DS) aims to generate large-scale heuristic labeling corpus, which is widely used for neural relation extraction currently. However, it heavily suffers from noisy labeling and long-tail distributions problem. Many advanced approaches usually separately address two problems, which ignore their mutual interactions. In this paper, we propose a novel framework named RH-Net, which utilizes Reinforcement learning and Hierarchical relational searching module to improve relation extraction. We leverage reinforcement learning to instruct the model to select high-quality instances. We then propose the hierarchical relational searching module to share the semantics from correlative instances between data-rich and data-poor classes. During the iterative process, the two modules keep interacting to alleviate the noisy and long-tail problem simultaneously. Extensive experiments on widely used NYT data set clearly show that our method significant improvements over state-of-the-art baselines.


title: Unraveling the origin of higher success probabilities in quantum versus semi-classical annealing

abstract: Quantum annealing aims at finding optimal solutions to complex optimization problems using a suitable quantum many body Hamiltonian encoding the solution in its ground state. To find the solution one typically evolves the ground state of a soluble initial Hamiltonian adiabatically to the ground state of the designated final Hamiltonian. Here we explore whether and when a full quantum representation of the dynamics leads to higher probability to end up in the desired ground when compared to a classical mean field approximation. As simple, nontrivial example we target the ground state of interacting bosons trapped in a tight binding lattice with small local defect by turning on long range interactions. Already two atoms in four sites interacting via two cavity modes prove complex enough to exhibit significant differences between the full quantum model and a mean field approximation for the cavity fields mediating the interactions. We find a large parameter region of highly successful quantum annealing where the semi-classical approach largely fails. Here we see strong evidence for the importance of entanglement to end close to the optimal solution. The quantum model also reduces the minimal time for a high target occupation probability. In contrast to naive expectations that enlarging the Hilbert space is beneficial, different numerical cut-offs of the Hilbert space reveal an improved performance for lower cut-offs, i.e. an nonphysical reduced Hilbert space, for short simulation times. Hence a less faithful representation of the full quantum dynamics sometimes creates a higher numerical success probability in shorter time. However, a sufficiently high cut-off proves relevant to obtain near perfect fidelity for long simulations times in a single run. Overall our results exhibit a clear improvement based on a quantum model versus simulations based on a classical field approximation.


title: Algorithms for $q$-ary Error-Correcting Codes with Limited Magnitude and Feedback

abstract: Berlekamp and Zigangirov completely determined the capacity error function for binary error correcting codes with noiseless feedback. It is still an unsolved problem if the upper bound for the capacity error function in the non-binary case of Ahlswede, Lebedev, and Deppe is sharp. We consider wraparound channels with limited magnitude and noiseless feedback. We completely determine the capacity error function for all $q$-ary wraparound channels with a magnitude of level $r$. All of our algorithms use partial noiseless feedback. Furthermore, a special case of the problem is equivalent to Shannon's zero-error problem.


title: Physics-Based Deep Learning for Fiber-Optic Communication Systems

abstract: We propose a new machine-learning approach for fiber-optic communication systems whose signal propagation is governed by the nonlinear Schr\"odinger equation (NLSE). Our main observation is that the popular split-step method (SSM) for numerically solving the NLSE has essentially the same functional form as a deep multi-layer neural network; in both cases, one alternates linear steps and pointwise nonlinearities. We exploit this connection by parameterizing the SSM and viewing the linear steps as general linear functions, similar to the weight matrices in a neural network. The resulting physics-based machine-learning model has several advantages over "black-box" function approximators. For example, it allows us to examine and interpret the learned solutions in order to understand why they perform well. As an application, low-complexity nonlinear equalization is considered, where the task is to efficiently invert the NLSE. This is commonly referred to as digital backpropagation (DBP). Rather than employing neural networks, the proposed algorithm, dubbed learned DBP (LDBP), uses the physics-based model with trainable filters in each step and its complexity is reduced by progressively pruning filter taps during gradient descent. Our main finding is that the filters can be pruned to remarkably short lengths-as few as 3 taps/step-without sacrificing performance. As a result, the complexity can be reduced by orders of magnitude in comparison to prior work. By inspecting the filter responses, an additional theoretical justification for the learned parameter configurations is provided. Our work illustrates that combining data-driven optimization with existing domain knowledge can generate new insights into old communications problems.


title: Superconductivity and Fermi Surface Nesting in the Candidate Dirac Semimetal NbC

abstract: We report the synthesis of single-crystal NbC, a transition metal carbide with various unusual properties. Transport, magnetic susceptibility, and specific heat measurements demonstrate that NbC is a conventional superconductor with a superconducting transition temperature ($T_c$) of 11.5 K. Our theoretical calculations show that NbC is a type-II Dirac semimetal with strong Fermi surface nesting, which is supported by our ARPES measurement results. We also observed the superconducting gaps of NbC using angle-resolved photoemission spectroscopy (ARPES) and found some unconventional behaviors. These intriguing superconducting and topological properties, combined with the high corrosion resistance, make NbC an ideal platform for both fundamental research and device applications.


title: Concentric mixtures of Mallows models for top-$k$ rankings: sampling and identifiability

abstract: In this paper, we consider mixtures of two Mallows models for top-$k$ rankings, both with the same location parameter but with different scale parameters, i.e., a mixture of concentric Mallows models. This situation arises when we have a heterogeneous population of voters formed by two homogeneous populations, one of which is a subpopulation of expert voters while the other includes the non-expert voters. We propose efficient sampling algorithms for Mallows top-$k$ rankings. We show the identifiability of both components, and the learnability of their respective parameters in this setting by, first, bounding the sample complexity for the Borda algorithm with top-$k$ rankings and second, proposing polynomial time algorithm for the separation of the rankings in each component. Finally, since the rank aggregation will suffer from a large amount of noise introduced by the non-expert voters, we adapt the Borda algorithm to be able to recover the ground truth consensus ranking which is especially consistent with the expert rankings.


title: A Simple and Efficient Registration of 3D Point Cloud and Image Data for Indoor Mobile Mapping System

abstract: Registration of 3D LiDAR point clouds with optical images is critical in the combination of multi-source data. Geometric misalignment originally exists in the pose data between LiDAR point clouds and optical images. To improve the accuracy of the initial pose and the applicability of the integration of 3D points and image data, we develop a simple but efficient registration method. We firstly extract point features from LiDAR point clouds and images: point features is extracted from single-frame LiDAR and point features from images using classical Canny method. Cost map is subsequently built based on Canny image edge detection. The optimization direction is guided by the cost map where low cost represents the the desired direction, and loss function is also considered to improve the robustness of the the purposed method. Experiments show pleasant results.


title: Above-ground biomass change estimation using national forest inventory data with Sentinel-2 and Landsat 8

abstract: This study aimed at estimating total forest above-ground net change (Delta AGB, Mt) over five years (2014-2019) based on model-assisted estimation utilizing freely available satellite imagery. The study was conducted for a boreal forest area (approx. 1.4 Mill hectares) in Norway where bi-temporal national forest inventory (NFI), Sentinel-2, and Landsat data were available. Biomass change was modelled based on a direct approach. The precision of estimates using only the NFI data in a basic expansion estimator were compared to four different alternative model-assisted estimates using 1) Sentinel-2 or Landsat data, and 2) using bi- or uni-temporal remotely sensed data. We found that the use of remotely sensed data improved the precision of the purely field-based estimates by a factor of up to three. The most precise estimates were found for the model-assisted estimation using bi-temporal Sentinel-2 (standard error; SE= 1.7 Mt). However, the decrease in precision when using Landsat data was small (SE= 1.92 Mt). In addition, we found that Delta AGB could be precisely estimated also when remotely sensed data were available only at the end of the monitoring period. We conclude that satellite optical data can considerably improve Delta AGB estimates, even in those cases where repeated and coincident NFI data are available. The free availability, global coverage, frequent update, and long-term time horizon make data from programs such as Sentinel-2 and Landsat a valuable data source for a consistent and durable monitoring of forest carbon dynamics.


title: Source Enumeration via RMT Estimator Based on Linear Shrinkage Estimation of Noise Eigenvalues Using Relatively Few Samples

abstract: Estimating the number of signals embedded in noise is a fundamental problem in array signal processing. The classic RMT estimator based on random matrix theory (RMT) tends to under-estimate the number of signals as it does not consider the non-negligible bias term among eigenvalues for finite sample size. Moreover, the RMT estimator suffers from uncertainty in noise variance estimation problem. In order to overcome these problems, we firstly derive a more accurate expression for the distribution of the sample eigenvalues and the bias term among eigenvalues by utilizing the linear shrinkage (LS) estimate of noise sample eigenvalues. Then, we analyze the effect of the bias term among eigenvalues on the estimation performance of the RMT estimator, and derive the increased under-estimation probability of the RMT estimator incurred by this bias term. Based on these results, we propose a novel RMT estimator based on LS estimate of noise eigenvalues (termed as LS-RMT estimator) by incorporating the bias term into the decision criterion of the RMT estimator. As the LS-RMT estimator incorporates this bias term among eigenvalues into the decision criterion of the RMT estimator, it can detect signal eigenvalues immersed in this bias term. Therefore, the LS-RMT estimator can overcome the higher under-estimation probability of the RMT estimator incurred by the bias term among eigenvalues, and also avoids the uncertainty in the noise variance estimation suffered by the RMT estimator as the noise variance is estimated under the assumption that the eigenvalue being tested is arising from noise. Finally, extensive simulation results are presented to show that the proposed LS-RMT estimator outperforms the existing estimators.


title: Wild Local Structures of Automorphic Lie Algebras

abstract: We study automorphic Lie algebras using a family of evaluation maps parametrised by the representations of the associative algebra of functions. This provides a descending chain of ideals for the automorphic Lie algebra which is used to prove that it is of wild representation type. We show that the associated quotients of the automorphic Lie algebra are isomorphic to twisted truncated polynomial current algebras. When a simple Lie algebra is used in the construction, this allows us to describe the local Lie structure of the automorphic Lie algebra in terms of affine Kac-Moody algebras.


title: A Weaker Faithfulness Assumption based on Triple Interactions

abstract: One of the core assumptions in causal discovery is the faithfulness assumption, i.e., assuming that independencies found in the data are due to separations in the true causal graph. This assumption can, however, be violated in many ways, including xor connections, deterministic functions or cancelling paths. In this work, we propose a weaker assumption that we call $2$-adjacency faithfulness. In contrast to adjacency faithfulness, which assumes that there is no conditional independence between each pair of variables that are connected in the causal graph, we only require no conditional independence between a node and a subset of its Markov blanket that can contain up to two nodes. Equivalently, we adapt orientation faithfulness to this setting. We further propose a sound orientation rule for causal discovery that applies under weaker assumptions. As a proof of concept, we derive a modified Grow and Shrink algorithm that recovers the Markov blanket of a target node and prove its correctness under strictly weaker assumptions than the standard faithfulness assumption.


title: End-to-end trainable network for degraded license plate detection via vehicle-plate relation mining

abstract: License plate detection is the first and essential step of the license plate recognition system and is still challenging in real applications, such as on-road scenarios. In particular, small-sized and oblique license plates, mainly caused by the distant and mobile camera, are difficult to detect. In this work, we propose a novel and applicable method for degraded license plate detection via vehicle-plate relation mining, which localizes the license plate in a coarse-to-fine scheme. First, we propose to estimate the local region around the license plate by using the relationships between the vehicle and the license plate, which can greatly reduce the search area and precisely detect very small-sized license plates. Second, we propose to predict the quadrilateral bounding box in the local region by regressing the four corners of the license plate to robustly detect oblique license plates. Moreover, the whole network can be trained in an end-to-end manner. Extensive experiments verify the effectiveness of our proposed method for small-sized and oblique license plates. Codes are available at https://github.com/chensonglu/LPD-end-to-end.


title: Galaxy Skew-Spectra in Redshift-Space

abstract: Modern galaxy surveys focus on the galaxy power spectrum or 2-point correlation function to test and constrain cosmological models. Additional information comes from higher-order N-point functions, but their analysis is challenging. A simple solution is to compute the cross-power spectrum between the squared galaxy density and the galaxy density. Being simple to measure and to plot, this skew-spectrum shares many of the familiar useful properties of the standard galaxy power spectrum. We show that by computing multiple quadratic fields and correlating them with the density, all contributions to the tree-level redshift-space galaxy bispectrum can be captured with skew-spectra. Using synthetic datasets, we show that our measurement pipeline matches analytical predictions, and that their dependence on galaxy bias parameters and the logarithmic growth rate is as expected theoretically.


title: Random Shifting Intelligent Reflecting Surface for OTP Encrypted Data Transmission

abstract: In this paper, we propose a novel encrypted data transmission scheme using an intelligent reflecting surface (IRS) to generate secret keys in wireless communication networks. We show that perfectly secure one-time pad (OTP) communications can be established by using a simple random phase shifting of the IRS elements. To maximize the secure transmission rate, we design an optimal time slot allocation algorithm for the IRS secret key generation and the encrypted data transmission phases. Moreover, a theoretical expression of the key generation rate is derived based on Poisson point process (PPP) for the practical scenario when eavesdroppers' channel state information (CSI) is unavailable. Simulation results show that employing our IRS-based scheme can significantly improve the encrypted data transmission performance for a wide-range of wireless channel gains and system parameters.


title: Leveraging speaker attribute information using multi task learning for speaker verification and diarization

abstract: Deep speaker embeddings have become the leading method for encoding speaker identity in speaker recognition tasks. The embedding space should ideally capture the variations between all possible speakers, encoding the multiple acoustic aspects that make up a speaker's identity, whilst being robust to non-speaker acoustic variation. Deep speaker embeddings are normally trained discriminatively, predicting speaker identity labels on the training data. We hypothesise that additionally predicting speaker-related auxiliary variables -- such as age and nationality -- may yield representations that are better able to generalise to unseen speakers. We propose a framework for making use of auxiliary label information, even when it is only available for speech corpora mismatched to the target application. On a test set of US Supreme Court recordings, we show that by leveraging two additional forms of speaker attribute information derived respectively from the matched training data, and VoxCeleb corpus, we improve the performance of our deep speaker embeddings for both verification and diarization tasks, achieving a relative improvement of 26.2% in DER and 6.7% in EER compared to baselines using speaker labels only. This improvement is obtained despite the auxiliary labels having been scraped from the web and being potentially noisy.


title: A Method of Generating Measurable Panoramic Image for Indoor Mobile Measurement System

abstract: This paper designs a technique route to generate high-quality panoramic image with depth information, which involves two critical research hotspots: fusion of LiDAR and image data and image stitching. For the fusion of 3D points and image data, since a sparse depth map can be firstly generated by projecting LiDAR point onto the RGB image plane based on our reliable calibrated and synchronized sensors, we adopt a parameter self-adaptive framework to produce 2D dense depth map. For image stitching, optimal seamline for the overlapping area is searched using a graph-cuts-based method to alleviate the geometric influence and image blending based on the pyramid multi-band is utilized to eliminate the photometric effects near the stitching line. Since each pixel is associated with a depth value, we design this depth value as a radius in the spherical projection which can further project the panoramic image to the world coordinate and consequently produces a high-quality measurable panoramic image. The purposed method is tested on the data from our data collection platform and presents a satisfactory application prospects.


title: Cross-lingual Machine Reading Comprehension with Language Branch Knowledge Distillation

abstract: Cross-lingual Machine Reading Comprehension (CLMRC) remains a challenging problem due to the lack of large-scale annotated datasets in low-source languages, such as Arabic, Hindi, and Vietnamese. Many previous approaches use translation data by translating from a rich-source language, such as English, to low-source languages as auxiliary supervision. However, how to effectively leverage translation data and reduce the impact of noise introduced by translation remains onerous. In this paper, we tackle this challenge and enhance the cross-lingual transferring performance by a novel augmentation approach named Language Branch Machine Reading Comprehension (LBMRC). A language branch is a group of passages in one single language paired with questions in all target languages. We train multiple machine reading comprehension (MRC) models proficient in individual language based on LBMRC. Then, we devise a multilingual distillation approach to amalgamate knowledge from multiple language branch models to a single model for all target languages. Combining the LBMRC and multilingual distillation can be more robust to the data noises, therefore, improving the model's cross-lingual ability. Meanwhile, the produced single multilingual model is applicable to all target languages, which saves the cost of training, inference, and maintenance for multiple models. Extensive experiments on two CLMRC benchmarks clearly show the effectiveness of our proposed method.


title: High Resolution Dielectric Characterization of Single Cells and Microparticles Using Integrated Microfluidic Microwave Sensors

abstract: Microwave sensors can probe intrinsic material properties of analytes in a microfluidic channel at physiologically relevant ion concentrations. While microwave sensors have been used to detect single cells and microparticles in earlier studies, the synergistic use and comparative analysis of microwave sensors with optical microscopy for material classification and size tracking applications have been scarcely investigated so far. Here we combined microwave and optical sensing to differentiate microscale objects based on their dielectric properties. We designed and fabricated two types of planar sensor: a Coplanar Waveguide Resonator (CPW) and a Split-Ring Resonator (SRR). Both sensors possessed sensing electrodes with a narrow gap to detect single cells passing through a microfluidic channel integrated on the same chip. We also show that standalone microwave sensors can track the relative changes in cellular size in real-time. In sensing single 20-micron diameter polystyrene particles, Signal-to-Noise ratio values of approximately 100 for CPW and 70 for SRR sensors were obtained. These findings demonstrate that microwave sensing technology can serve as a complementary technique for single-cell biophysical experiments and microscale pollutant screening.


title: $1/2$-conjectures on the domination game and claw-free graphs

abstract: Let $\gamma_g(G)$ be the game domination number of a graph $G$. Rall conjectured that if $G$ is a traceable graph, then $\gamma_g(G) \le \left\lceil \frac{1}{2}n(G)\right\rceil$. Our main result verifies the conjecture over the class of line graphs. Moreover, in this paper we put forward the conjecture that if $\delta(G) \geq 2$, then $\gamma_g(G) \leq \left\lceil \frac{1}{2}n(G) \right\rceil$. We show that both conjectures hold true for claw-free cubic graphs. We further prove the upper bound $\gamma_g(G) \le \left\lceil \frac{11}{20} \, n(G) \right\rceil$ over the class of claw-free graphs of minimum degree at least $2$. Computer experiments supporting the new conjecture and sharpness examples are also presented.


title: Behavior Priors for Efficient Reinforcement Learning

abstract: As we deploy reinforcement learning agents to solve increasingly challenging problems, methods that allow us to inject prior knowledge about the structure of the world and effective solution strategies becomes increasingly important. In this work we consider how information and architectural constraints can be combined with ideas from the probabilistic modeling literature to learn behavior priors that capture the common movement and interaction patterns that are shared across a set of related tasks or contexts. For example the day-to day behavior of humans comprises distinctive locomotion and manipulation patterns that recur across many different situations and goals. We discuss how such behavior patterns can be captured using probabilistic trajectory models and how these can be integrated effectively into reinforcement learning schemes, e.g.\ to facilitate multi-task and transfer learning. We then extend these ideas to latent variable models and consider a formulation to learn hierarchical priors that capture different aspects of the behavior in reusable modules. We discuss how such latent variable formulations connect to related work on hierarchical reinforcement learning (HRL) and mutual information and curiosity based objectives, thereby offering an alternative perspective on existing ideas. We demonstrate the effectiveness of our framework by applying it to a range of simulated continuous control domains.


title: Effects of early intense bombardment on megaregolith evolution and on lunar (and planetary) surface samples

abstract: Impact rates in the first 500 Myr of the solar system are critical to an understanding of lunar geological history, but they have been controversial. The widely accepted, post-Apollo paradigm of early lunar impact cratering (ca. 1975-2014) proposed very low or negligible impact cratering in the period from accretion (>4.4 Ga) to about 4.0 Ga ago, followed by a 170-million-year-long spike of cataclysmic cratering, during which most prominent multi-ring impact basins formed at age of about 3.9 Ga. More recent dynamical models suggest very early intense impact rates, declining throughout the period from accretion until an age of about 3.0 Ga. These models remove the basin-forming spike. This shift has important consequences on megaregolith evolution and properties of rock samples that can be collected on the lunar surface today. We adopt the Morbidelli et al. (2018) "accretion tail" model of early intense bombardment, declining as a function of time. We find effects differing from the previous models: early crater saturation and supersaturation; disturbance of magma ocean solidification; deep early megaregolith; and erosive destruction of the earliest multi-ring basins, their impact melts, and their ejecta blankets. Our results explain observations such as differences in numbers of early lunar impact melts vs. numbers of early igneous crustal rocks, highland breccias containing impact melts as old as 4.35 Ga, absence of a 170 Myr-long spike in impact melt ages at 3.9 Ga among lunar and asteroidal meteorites, and GRAIL observations of lunar crustal structure.


title: Synthetic superfluid chemistry with vortex-trapped quantum impurities

abstract: We explore the effect of using two-dimensional matter-wave vortices to confine an ensemble of bosonic quantum impurities. This is modelled theoretically using a mass-imbalanced homogeneous two component Gross-Pitaevskii equation where each component has independent atom numbers and equal atomic masses. By changing the mass imbalance of our system we find the shape of the vortices are deformed even at modest imbalances, leading to barrel shaped vortices; which we quantify using a multi-component variational approach. The energy of impurity carrying vortex pairs are computed, revealing a mass-dependent energy splitting. We then compute the excited states of the impurity, which we in turn use to construct `covalent bonds' for vortex pairs. Our work opens a new route to simulating synthetic chemical reactions with superfluid systems.


title: Internal contact modeling for finite strain topology optimization

abstract: The present work proposes an extension of the third medium contact method for solving structural topology optimization problems that involve and exploit self-contact. A new regularization of the void region, which acts as the contact medium, makes the method suitable for cases with very large deformations. The proposed contact method is implemented in a second order topology optimization framework, which employs a coupled simultaneous solution of the mechanical, design update, and adjoint problems. All three problems are derived and presented in weak form, and discretized with finite elements of suitable order. The capabilities and accuracy of the developed method are demonstrated in a topology optimization problem for achieving a desired non-linear force-displacement path.


title: Dark matter constraints and the neutralino sector of the scNMSSM

abstract: The neutralino sector of the semi-constrained next-to-minimal supersymmetric standard model is explored under recent experimental constraints, with special attention to dark matter (DM) limits. The effects of the upper and lower bounds of dark matter relic density and recent direct detection constraints on spin-independent and -dependent cross-sections are thoroughly analyzed. Particularly, we show which regions of the parameter space are ruled out due to the different dark matter constraints and the corresponding model-specific parameters: $\lambda, \kappa, A_{\lambda}$, and $A_{\kappa}$. We analyze all annihilation and co-annihilation processes (with heavier neutralinos and charginos) that contribute to the dark matter relic density. The mass components of the dark matter candidate, the lightest neutralino $\tilde{\chi}_1^0$, are studied, and the decays of heavy neutralinos and charginos, especially $\tilde{\chi}_2^0$ and $\tilde{\chi}_1^+$, into the lightest neutralino are examined. We impose semi-universal boundary conditions at the Grand Unified Theory scale and require a moderate range of $\tan{\beta} \lesssim 10$. We find that the allowed parameter space is associated with a heavy mass spectrum in general and that the lightest neutralino is mostly Higgsino with a mass range that resides mostly between 1000 and 1500 GeV. However, smaller mass values can be achieved if the DM candidate is bino-like or singlino-like.


title: Existence of the Metal-Rich Stellar Halo and High-velocity Thick Disk in the Galaxy

abstract: Based on the second Gaia data release (DR2), combined with the LAMOST and APOGEE spectroscopic surveys, we study the kinematics and metallicity distribution of the high-velocity stars that have a relative speed of at least 220 ${\rm km\ s^{-1}}$ with respect to the local standard of rest in the Galaxy. The rotational velocity distribution of the high-velocity stars with [Fe/H] $>-1.0$ dex can be well described by a two-Gaussian model, with peaks at $V_{\phi}\sim +164.2\pm0.7$ and $V_{\phi}\sim +3.0\pm0.3$ ${\rm km\ s^{-1}}$, associated with the thick disk and halo, respectively. This implies that there should exist a high-velocity thick disk (HVTD) and a metal-rich stellar halo (MRSH) in the Galaxy. The HVTD stars have the same position as the halo in the Toomre diagram and but show the same rotational velocity and metallicity as the canonical thick disk. The MRSH stars have basically the same rotational velocity, orbital eccentricity, and position in the Lindblad and Toomre diagram as the canonical halo stars, but they are more metal-rich. Furthermore, the metallicity distribution function (MDF) of our sample stars are well fitted by a four-Gaussian model, associated with the outer-halo, inner-halo, MRSH, and HVTD, respectively. Chemical and kinematic properties and age imply that the MRSH and HVTD stars may form in situ.


title: Two geometric lemmas for $S^{N-1}$-valued maps and an application to the homogenization of spin systems

abstract: We prove two geometric lemmas for $S^{N-1}$-valued functions that allow to modify sequences of lattice spin functions on a small percentage of nodes during a discrete-to-continuum process so as to have a fixed average. This is used to simplify known formulas for the homogenization of spin systems.


title: Two substructures in the nearby stellar halo found in Gaia and RAVE

abstract: We use the second Gaia data release (Gaia DR2), combined with RAVE spectroscopic surveys, to identify the substructures in the nearby stellar halo. We select 3,845 halo stars kinematically and chemically, and determine their density distribution in energy and angular momentum space. To select the substructures from overdensities, we reshuffle the velocities and estimate their significance. Two statistically significant substructures, GR-1 and GR-2, are identified. GR-1 has a high binding energy and small $z$-angular momentum. GR-2 is metal-rich but retrograde. They are both new substructure, and may be accretion debris of dwarf galaxies.


title: Active Learning for Human-in-the-Loop Customs Inspection

abstract: We study the human-in-the-loop customs inspection scenario, where an AI-assisted algorithm supports customs officers by recommending a set of imported goods to be inspected. If the inspected items are fraudulent, the officers can levy extra duties. Th formed logs are then used as additional training data for successive iterations. Choosing to inspect suspicious items first leads to an immediate gain in customs revenue, yet such inspections may not bring new insights for learning dynamic traffic patterns. On the other hand, inspecting uncertain items can help acquire new knowledge, which will be used as a supplementary training resource to update the selection systems. Based on multiyear customs datasets obtained from three countries, we demonstrate that some degree of exploration is necessary to cope with domain shifts in trade data. The results show that a hybrid strategy of selecting likely fraudulent and uncertain items will eventually outperform the exploitation-only strategy.


title: Associahedra, cyclohedra and inversion of power series

abstract: We introduce the Hopf monoid of sets of cycles and paths, which contains the Fa\`a di Bruno Hopf monoid as a submonoid. We give cancellation-free and grouping-free formulas for its antipode, one in terms of tubings and one in terms of \emph{pointed} noncrossing partitions. We provide an explicit description of the group of characters of this Hopf monoid in terms of pairs of power series. Using graph associahedra, we relate paths and cycles to associahedra and cyclohedra, respectively. We give formulas for inversion in the group of characters in terms of the faces of these polytopes.


title: An Analysis of Parton Distribution Functions of the Pion and the Kaon with the Maximum Entropy Input

abstract: We present pion and kaon parton distribution functions from a global QCD analysis of the experimental data within the framework of dynamical parton model. We use the DGLAP equations with parton-parton recombination corrections and the valence input of uniform distribution which maximizes the information entropy. At our input scale $Q_0^2$, there are no sea quark and gluon distributions. All the sea quarks and gluons of the pion and the kaon are completely generated from the parton splitting processes. The mass-dependent parton splitting kernel is applied for the strange quark distribution in the kaon. The obtained valence quark and sea quark distributions at high $Q^{2}$ ($Q^2>5$ GeV$^2$) are compatible with the existed experimental measurements. Furthermore, the asymptotic behaviours of parton distribution functions at small and large $x$ have been studied for both the pion and the kaon. Lastly, the first three moments of parton distributions at high $Q^{2}$ scale are calculated, which are consistent with other theoretical predictions.


title: An Interior Point-Proximal Method of Multipliers for Positive Semi-Definite Programming

abstract: In this paper we generalize the Interior Point-Proximal Method of Multipliers (IP-PMM) presented in [An Interior Point-Proximal Method of Multipliers for Convex Quadratic Programming, Computational Optimization and Applications, 78, 307--351 (2021)] for the solution of linear positive Semi-Definite Programming (SDP) problems, allowing inexactness in the solution of the associated Newton systems. In particular, we combine an infeasible Interior Point Method (IPM) with the Proximal Method of Multipliers (PMM) and interpret the algorithm (IP-PMM) as a primal-dual regularized IPM, suitable for solving SDP problems. We apply some iterations of an IPM to each sub-problem of the PMM until a satisfactory solution is found. We then update the PMM parameters, form a new IPM neighbourhood, and repeat this process. Given this framework, we prove polynomial complexity of the algorithm, under mild assumptions, and without requiring exact computations for the Newton directions. We furthermore provide a necessary condition for lack of strong duality, which can be used as a basis for constructing detection mechanisms for identifying pathological cases within IP-PMM.


title: Limitations of CMB B-mode template delensing

abstract: Efforts to detect a primordial $B$-mode of CMB polarization generated by inflationary gravitational waves ought to mitigate the large variance associated with the $B$-modes produced by gravitational lensing, a process known as delensing. A popular approach to delensing entails building a lensing $B$-mode template by mimicking the lensing operation, either at gradient order or non-perturbatively, using high-resolution $E$-mode observations and some proxy of the lensing potential. By explicitly calculating all contributions to two-loop order in lensing to the power spectrum of $B$-modes delensed with such a template in the noise-free limit, we are able to show that: (i) corrections to the leading-order calculation of the lensing $B$-mode power spectrum only enter at the $O(1)\,\%$ level because of extensive cancellations between large terms at next-to-leading order; (ii) these cancellations would disappear if a gradient-order template were to be built from unlensed or delensed $E$-modes, giving rise to a residual delensing floor of $O(10)\,\%$ of the original power; (iii) new cancellations arise when the lensed $E$-modes are used in the gradient-order template, allowing for the delensing floor to be as low as $O(1)\,\%$ of the original power in practical applications of this method; and (iv) these new cancellations would disappear for a non-perturbative template constructed from the lensed $E$-modes, reintroducing a residual delensing floor of $O(10)\,\%$. We further show that the gradient-order template outperforms the non-perturbative one in realistic scenarios with noisy estimates of the $E$-mode polarization and lensing potential. We therefore recommend that in practical applications of $B$-mode template delensing, where the template is constructed directly from the (filtered) observed $E$-modes, the gradient-order approach should be used rather than a non-perturbative remapping.