title: Randomized compiling for scalable quantum computing on a noisy superconducting quantum processor

abstract: The successful implementation of algorithms on quantum processors relies on the accurate control of quantum bits (qubits) to perform logic gate operations. In this era of noisy intermediate-scale quantum (NISQ) computing, systematic miscalibrations, drift, and crosstalk in the control of qubits can lead to a coherent form of error which has no classical analog. Coherent errors severely limit the performance of quantum algorithms in an unpredictable manner, and mitigating their impact is necessary for realizing reliable quantum computations. Moreover, the average error rates measured by randomized benchmarking and related protocols are not sensitive to the full impact of coherent errors, and therefore do not reliably predict the global performance of quantum algorithms, leaving us unprepared to validate the accuracy of future large-scale quantum computations. Randomized compiling is a protocol designed to overcome these performance limitations by converting coherent errors into stochastic noise, dramatically reducing unpredictable errors in quantum algorithms and enabling accurate predictions of algorithmic performance from error rates measured via cycle benchmarking. In this work, we demonstrate significant performance gains under randomized compiling for the four-qubit quantum Fourier transform algorithm and for random circuits of variable depth on a superconducting quantum processor. Additionally, we accurately predict algorithm performance using experimentally-measured error rates. Our results demonstrate that randomized compiling can be utilized to leverage and predict the capabilities of modern-day noisy quantum processors, paving the way forward for scalable quantum computing.


title: Conditional probability and interferences in generalized measurements with or without definite causal order

abstract: In the context of generalized measurement theory, the Gleason-Busch theorem assures the unique form of the associated probability function. Recently, in Flatt et al. Phys. Rev. A 96, 062125 (2017), the case of subsequent measurements has been treated, with the derivation of the L\"uders rule and its generalization (Kraus update rule). Here we investigate the special case of subsequent measurements where an intermediate measurement is a composition of two measurements (a or b) and the case where the causal order is not defined (a and b or b and a). In both cases interference effects can arise. We show that the associated probability cannot be written univocally, and the distributive property on its arguments cannot be taken for granted. The two probability expressions correspond to the Born rule and the classical probability; they are related to the intrinsic possibility of obtaining definite results for the intermediate measurement. For indefinite causal order, a causal inequality is also deduced. The frontier between the two cases is investigated in the framework of generalized measurements with a toy model, a Mach-Zehnder interferometer with a movable beam splitter.


title: CoVer: Collaborative Light-Node-Only Verification and Data Availability for Blockchains

abstract: Validating a blockchain incurs heavy computation, communication, and storage costs. As a result, clients with limited resources, called light nodes, cannot verify transactions independently and must trust full nodes, making them vulnerable to security attacks. Motivated by this problem, we ask a fundamental question: can light nodes securely validate without any full nodes? We answer affirmatively by proposing CoVer, a decentralized protocol that allows a group of light nodes to collaboratively verify blocks even under a dishonest majority, achieving the same level of security for block validation as full nodes while only requiring a fraction of the work. In particular, work per node scales down proportionally with the number of participants (up to a log factor), resulting in computation, communication, and storage requirements that are sublinear in block size. Our main contributions are light-node-only protocols for fraud proofs and data availability.


title: Phonon-assisted phase separation in strongly correlated systems

abstract: We relate the phase separation observed in many crystals with pronounced electron correlations to the regions of negative electron compressibility. They were found in several models describing strong electron correlations. At low temperatures, these regions arise near chemical potentials corresponding to the change of the ground state in the site Hamiltonian. The negative electron compressibility leads to the separation of the system into electron-rich and electron-poor domains. The energy released in the course of this separation is absorbed by phonons. Another role of phonons is to give a definite form -- stripes or checkerboards -- to lattice distortions and domains of different electron concentrations. The shape, direction, and periodicity of such textures are determined by wave vectors of lattice distortions, which most strongly scatter electrons.


title: 4D Dyck triangle and its projections

abstract: The classic Dyck triangle, the Catalan triangle, and the Catalan convolution matrix are plane projections of the multidimensional Dyck triangle. In the Dyck path, each node is uniquely determined by two of four interrelated parameters: (i) the position of the current parenthesis, (ii) the current unbalance of the parentheses, (iii) the number of viewed left parentheses, and (iv) the same for right parentheses. The last two parameters can be redefined, respectively, as the index of the current Catalan number and the index of the summand in the decomposition of the Catalan number into the sum of squares (Dyck squares). For the 4D Dyck triangle, we consider six 2D projections (some of them are not yet in demand) and four 3D projections.


title: Quantum Annealing Approaches to the Phase-Unwrapping Problem in Synthetic-Aperture Radar Imaging

abstract: The focus of this work is to explore the use of quantum annealing solvers for the problem of phase unwrapping of synthetic aperture radar (SAR) images. Although solutions to this problem exist based on network programming, these techniques do not scale well to larger-sized images. Our approach involves formulating the problem as a quadratic unconstrained binary optimization (QUBO) problem, which can be solved using a quantum annealer. Given that present embodiments of quantum annealers remain limited in the number of qubits they possess, we decompose the problem into a set of subproblems that can be solved individually. These individual solutions are close to optimal up to an integer constant, with one constant per sub-image. In a second phase, these integer constants are determined as a solution to yet another QUBO problem. We test our approach with a variety of software-based QUBO solvers and on a variety of images, both synthetic and real. Additionally, we experiment using D-Wave Systems's quantum annealer, the D-Wave 2000Q. The software-based solvers obtain high-quality solutions comparable to state-of-the-art phase-unwrapping solvers. We are currently working on optimally mapping the problem onto the restricted topology of the quantum annealer to improve the quality of the solution.


title: Enhanced nonlinear optical response of single metal-dielectric nanocavities resonating in the near-infrared

abstract: Harmonic generation mechanisms are of great interest in nanoscience and nanotechnology, since they allow generating visible light by using near-infrared radiation, which is particularly suitable for its endless applications in bio-nanophotonics and opto-electronics. In this context, multilayer metal-dielectric nanocavities are widely used for light confinement and waveguiding at the nanoscale. They exhibit intense and localized resonances that can be conveniently tuned in the near-infrared and are therefore ideal for enhancing nonlinear effects in this spectral range. In this work, we experimentally investigate the nonlinear optical response of multilayer metal-dielectric nanocavities. By engineering their absorption efficiency and exploiting their intrinsic interface-induced symmetry breaking, we achieve one order of magnitude higher second-harmonic generation efficiency compared to gold nanostructures featuring the same geometry and resonant behavior. In particular, while the third order nonlinear susceptibility is comparable with that of bulk Au, we estimate a second order nonlinear susceptibility of the order of 1 pm/V, which is comparable with that of typical nonlinear crystals. We envision that our system, which combines the advantages of both plasmonic and dielectric materials, might enable the realization of composite and multi-functional nano-systems for an efficient manipulation of nonlinear optical processes at the nanoscale.


title: Proportional reinsurance for fractional Brownian risk model

abstract: This paper investigates ruin probabilities for a two-dimensional fractional Brownian risk model with a proportional reinsurance scheme. We focus on joint and simultaneous ruin probabilities in a finite-time horizon. The risk processes of both insurance and reinsurance companies are composed of a large number of i.i.d. sub-risk processes, representing independent businesses. We derive the asymptotics as the initial capital tends to infinity.


title: An Enhanced Energy Management System Including a Real-Time Load-Redistribution Threat Analysis Tool and Cyber-Physical SCED

abstract: It is possible to launch undetectable load-redistribution (LR) attacks against power systems, even in systems with protection schemes. Therefore, detecting LR attacks in power systems and establishing a corrective action to provide secured operating points are imperative. In this paper, we develop a systematic real-time LR threat analysis (RTLRTA) tool, which can flag LR attacks and identify all affected transmission assets. Since attackers might use random deviations to create LR attacks, we introduce an optimization model to generate random LR attacks. Hence, we can determine accurate thresholds for our detection index and test the tool's functionality when there are random LR attacks. Additionally, based on an estimation for the actual loads in the post-attack stage, we design a set of physical line flow security constraints (PLFSCs) and add it to the security-constrained economic dispatch (SCED) model. We call the new model cyber-physical SCED (CPSCED), which can appropriately respond to the identified LR attacks and provide secured dispatch points. We generate multiple scenarios of random LR attacks and noise errors for different target lines in the $2383$-bus Polish test system to validate our proposed methods' accuracy and functionality in detecting LR attacks and responding to them.


title: Ubiquitous topological states of phonons in solids: Silicon as a model material

abstract: Research on topological physics of phonons has attracted enormous interest but demands appropriate model materials. Our {\it ab initio} calculations identify silicon as an ideal candidate material containing extraordinarily rich topological phonon states. In silicon, we identify various topological nodal lines protected by glide mirror or mirror symmetries and characterized by quantized Berry phase $\pi$, which gives drumhead surface states observable from any surface orientations. Remarkably, a novel type of topological nexus phonon is discovered, which is featured by double Fermi-arc-like surface states and distinguished from Weyl phonons by requiring neither inversion nor time-reversal symmetry breaking. Versatile topological states can be created from the nexus phonons, such as Hopf nodal link by strain. Furthermore, we generalize the symmetry analysis to other centrosymmetric systems and find numerous candidate materials, demonstrating the ubiquitous existence of topological phonons in solids. These findings open up new opportunities for studying topological phonons in realistic materials and their influence on surface physics.


title: Robust edge states in magnetic domain-wall racetrack

abstract: Controllable artificial pinning is indispensable in numerous domain-wall (DW) devices, such as memory, sensor, logic gate, and neuromorphic computing hardware. The high-accuracy determination of the effective spring constant of the pinning potential, however, remains challenging, because the extrinsic pinning is often mixed up with intrinsic ones caused by materials defects and randomness. Here, we study the collective dynamics of interacting DWs in a racetrack with pinning sites of alternate distances. By mapping the governing equations of DW motion to the Su-Schrieffer-Heeger model and evaluating the quantized Zak phase, we predict two topologically distinct phases in the racetrack. Robust edge state emerges at either one or both ends depending on the parity of the DW number and the ratio of alternating intersite lengths. We show that the in-gap DW oscillation frequency has a fixed value which depends only on the geometrical shape of the pinning notch, and is insensitive to device imperfections and inhomogeneities. We propose to accurately quantify the spring coefficient that equals the square of the robust DW frequency multiplied by its constant mass. Our findings suggest as well that the DW racetrack is an ideal platform to study the topological phase transition.


title: Spectral asymptotics for the semiclassical Bochner Laplacian of a line bundle with constant rank curvature

abstract: The goal of this paper is manyfold. Firstly, we want to give a short introduction to the Bochner Laplacian and explain why it acts locally as a magnetic Laplacian. Secondly, given a confining magnetic field, we use Agmon-like estimates to reduce its spectral study to magnetic Laplacians, in the semiclassical limit. Finally, we use this to translate already-known spectral asymptotics for the magnetic Laplacian to the Bochner Laplacian.


title: On evolution equations for Lie groupoids

abstract: Using the calculus of Fourier integral operators on Lie groupoids developped in [18], we study the fundamental solution of the evolution equation ($\partial$ $\partial$t + iP)u = 0 where P is a self adjoint elliptic order one G-pseudodifferential operator on the Lie groupoid G. Along the way, we continue the study of distributions on Lie groupoids done in [17] by adding the reduced C *-algebra of G in the picture and we investigate the local nature of the regularizing operators of [32].


title: A Method of Using Quantum Kinematic Gravitational Effect to Describe the Anomalous Acceleration of 1I/2017 U1'Oumuamua

abstract: The Newtonian gravity considers gravity is a force that action at a distance, however, if gravity is quantized, then it should not be an instantaneous force, but propagates toward a direction at the speed of light by gravitational quanta. The quantum kinematic gravitational effect may be able to adapt the Newton's gravitational equation to quantum gravity theory, and may improve its accuracy in some dynamic scenes. Here is a possible method of using quantum kinematic gravitational effect to describe the anomalous acceleration of 1I/2017 U1'Oumuamua, the margin of error between the calculated result and the astronomical observations during Oct 19, 2017 - May 3, 2018 is approximately 8% - 10%. The accumulation of extra velocity caused by quantum kinematic gravitational effect may provide extra kinetic energy, which may be one of the reasons the Moon orbital expansion.


title: $3$-setwise intersecting families of the symmetric group

abstract: Given two positive integers $n\geq 3$ and $t\leq n$, the permutations $\sigma,\pi \in \operatorname{Sym}(n)$ are $t$-setwise intersecting if they agree (setwise) on a $t$-subset of $\{1,2,\ldots,n\}$. A family $\mathcal{F} \subset \operatorname{Sym}(n)$ is $t$-setwise intersecting if any two permutations of $\mathcal{F}$ are $t$-setwise intersecting. Ellis [Journal of Combinatorial Theory, Series A, 119(4), 825--849, 2012] conjectured that if $t\leq n$ and $\mathcal{F} \subset \operatorname{Sym}(n)$ is a $t$-setwise intersecting family, then $|\mathcal{F}|\leq t!(n-t)!$ and equality holds only if $\mathcal{F}$ is a coset of a setwise stablizer of a $t$-subset of $\{1,2,\ldots,n\}$. In this paper, we prove that if $n\geq 11$ and $\mathcal{F}$ is $3$-setwise intersecting, then $|\mathcal{F}|\leq 6(n-3)!$. Moreover, we prove that the characteristic vector of a $3$-setwise intersecting family of maximum size lies in the sum of the eigenspaces induced by the permutation module of $\operatorname{Sym}(n)$ acting on the $3$-subsets of $\{1,2,\ldots,n\}$.


title: A perfect X-ray beam splitter and its applications to time-domain interferometry and quantum optics exploiting free-electron lasers

abstract: X-ray free-electron lasers (FEL) deliver ultrabright X-ray pulses, but not the sequences of phase-coherent pulses required for time-domain interferometry and control of quantum states. For conventional split-and-delay schemes to produce such sequences the challenge stems from extreme stability requirements when splitting Angstrom wavelength beams where tiniest path length differences introduce phase jitter. We describe an FEL mode based on selective electron bunch degradation and transverse beam shaping in the accelerator, combined with a self-seeded photon emission scheme. Instead of splitting the photon pulses after their generation by the FEL, we split the electron bunch in the accelerator, prior to photon generation, to obtain phase-locked X-ray pulses with sub-femtosecond duration. Time-domain interferometry becomes possible, enabling the concomitant program of classical and quantum optics experiments with X-rays. The scheme leads to new scientific benefits of cutting-edge FELs with attosecond and/or high-repetition rate capabilities, ranging from the X-ray analog of Fourier transform infrared spectroscopy to damage-free measurements.

categories: physics.acc-ph cond-mat.mtrl-sci cond-mat.str-el physics.atm-clus quant-ph


title: Deep Group-wise Variational Diffeomorphic Image Registration

abstract: Deep neural networks are increasingly used for pair-wise image registration. We propose to extend current learning-based image registration to allow simultaneous registration of multiple images. To achieve this, we build upon the pair-wise variational and diffeomorphic VoxelMorph approach and present a general mathematical framework that enables both registration of multiple images to their geodesic average and registration in which any of the available images can be used as a fixed image. In addition, we provide a likelihood based on normalized mutual information, a well-known image similarity metric in registration, between multiple images, and a prior that allows for explicit control over the viscous fluid energy to effectively regularize deformations. We trained and evaluated our approach using intra-patient registration of breast MRI and Thoracic 4DCT exams acquired over multiple time points. Comparison with Elastix and VoxelMorph demonstrates competitive quantitative performance of the proposed method in terms of image similarity and reference landmark distances at significantly faster registration.


title: Analysis of the weighted kappa and its maximum with Markov moves

abstract: In this paper the notion of Markov move from Algebraic Statistics is used to analyze the weighted kappa indices in rater agreement problems. In particular, the problem of the maximum kappa and its dependence on the choice of the weighting schemes are discussed. The Markov moves are also used in a simulated annealing algorithm to actually find the configuration of maximum agreement.


title: Lupus DANCe. Census of stars and 6D structure with Gaia-DR2 data

abstract: Context: Lupus is recognised as one of the closest star-forming regions, but the lack of trigonometric parallaxes in the pre-Gaia era hampered many studies on the kinematic properties of this region and led to incomplete censuses of its stellar population. Aims: We use the second data release of the Gaia space mission combined with published ancillary radial velocity data to revise the census of stars and investigate the 6D structure of the Lupus complex. Methods: We performed a new membership analysis of the Lupus association based on astrometric and photometric data over a field of 160 deg2 around the main molecular clouds of the complex and compared the properties of the various subgroups in this region. Results: We identified 137 high-probability members of the Lupus association of young stars, including 47 stars that had never been reported as members before. Many of the historically known stars associated with the Lupus region identified in previous studies are more likely to be field stars or members of the adjacent Scorpius-Centaurus association. Our new sample of members covers the magnitude and mass range from G=8 to G=18 mag and from 0.03 to 2.4Msun, respectively. We compared the kinematic properties of the stars projected towards the molecular clouds Lupus 1 to 6 and showed that these subgroups are located at roughly the same distance (about 160~pc) and move with the same spatial velocity. Our age estimates inferred from stellar models show that the Lupus subgroups are coeval (with median ages ranging from about 1 to 3 Myr). The Lupus association appears to be younger than the population of young stars in the Corona-Australis star-forming region recently investigated by our team using a similar methodology. The initial mass function of the Lupus association inferred from the distribution of spectral types shows little variation compared to other star-forming regions.


title: Prediction of two-dimensional ferromagnetic ferroelectric VOF$_2$ monolayer

abstract: Nowadays, designing and searching for materials with multiple functional characteristics are the key to achieving high-performance electronic devices. Among many candidates, two-dimensional multiferroic materials have great potential to be applied in highly-integrated magnetoelectric devices, such as high-density non-volatile memories. Here, we predict a two-dimensional material VOF2 monolayer to possess intrinsic ferroelectric and ferromagnetic properties. The VOF$_2$ monolayer own the largest in-plane ferroelectric polarization ($332$ pC/m) in the family of VO$X_2$ ($X$: halogen) oxyhalides. Different from other VO$X_2$ monolayers whose magnetic ground states are antiferromagnetic or noncollinear spiral textures, the VOF$_2$ monolayer owns a robust ferromagnetic ground state, which is rare but highly desirable. Our theoretical prediction provides a good candidate and starting point for the further pursuit of more two-dimensional multiferroic materials with high-performance magnetoelectricity.


title: The nature of the Milky Way's stellar halo revealed by the three integrals of motion

abstract: We developed a new selection method of halo stars in the phase-space distribution defined by the three integrals of motion in an axisymmetric Galactic potential, ($E$, $L_z$, $I_3$), where $I_3$ is the third integral of motion. The method is used to explore the general chemo-dynamical structure of the halo based on stellar samples from SDSS-SEGUE DR7 and DR16-APOGEE, matched with Gaia-DR2. We found, (a) halo stars can be separated from disk stars by selecting over (1) $0 < L_z < 1500$ \kpckms, $(2I_3)^{1/2} > 1000$ \kpckms (orbital angle $\theta_{\rm orb}$ $>$ 15-20 deg), and $E < -1.5 \times 10^5$ km$^2$ s$^{-2}$, and (2) $L_z < 0$ \kpckms. These selection criteria are free from kinematical biases introduced by the simple high-velocity cuts adopted in recent literature; (b) the averaged, or {\it coarse-grained}, halo phased-space distribution shows a monotonic exponential decrease with increasing $E$ and $I_3$ like the Michie-Bodenheimer models; (c) the inner stellar halo described in \citet{carollo2007,carollo2010} is found to comprise a combination of Gaia Enceladus debris (GE; \citealt{helmi2018}), lowest-$E$ stars (likely in-situ stars), as well as metal-poor prograde stars missed by the high velocity cuts selection; (d) the very metal poor outer halo, ([Fe/H] $< -$2.2), exhibits both retrograde and prograde rotation, with an asymmetric $L_z$ distribution towards high retrograde motions, and larger $\theta_{\rm orb}$ than those possessed by the GE dominated inner halo; (e) the Sgr dSph galaxy could induce a long-range dynamical effect on local halo stars. Implication for the formation of the stellar halo are also discussed.


title: A Niching Indicator-Based Multi-modal Many-objective Optimizer

abstract: Multi-modal multi-objective optimization is to locate (almost) equivalent Pareto optimal solutions as many as possible. Some evolutionary algorithms for multi-modal multi-objective optimization have been proposed in the literature. However, there is no efficient method for multi-modal many-objective optimization, where the number of objectives is more than three. To address this issue, this paper proposes a niching indicator-based multi-modal multi- and many-objective optimization algorithm. In the proposed method, the fitness calculation is performed among a child and its closest individuals in the solution space to maintain the diversity. The performance of the proposed method is evaluated on multi-modal multi-objective test problems with up to 15 objectives. Results show that the proposed method can handle a large number of objectives and find a good approximation of multiple equivalent Pareto optimal solutions. The results also show that the proposed method performs significantly better than eight multi-objective evolutionary algorithms.


title: Performance studies of Time-of-Flight detectors at LHC

abstract: We present results of a toy model study of performance of the Time-of-Flight detectors integrated into forward proton detectors. The goal of the ToF device is to suppress effects of additional soft processes (so called pile-up) accompanying the hard-scale central diffractive event, characterized by two tagged leading protons, one on each side from the interaction point. The method of mitigation of the pile-up effects exemplified in this study is based on measuring a difference between arrival times of these leading protons at the forward proton detectors and hence estimate the z-coordinate of the production vertex. We evaluate effects of the the pile-up background by studying in detail its components, and estimate the performance of the ToF method as a function of the time and spatial resolution of the ToF device and of the number of pile-up interactions per bunch crossing. We also propose a new observable with a potential to efficiently separate central diffractive signal from the harsh pile-up environment.


title: Multi-grained Semantics-aware Graph Neural Networks

abstract: Graph Neural Networks (GNNs) are powerful techniques in representation learning for graphs and have been increasingly deployed in a multitude of different applications that involve node- and graph-wise tasks. Most existing studies solve either the node-wise task or the graph-wise task independently while they are inherently correlated. This work proposes a unified model, AdamGNN, to interactively learn node and graph representations in a mutual-optimisation manner. Compared with existing GNN models and graph pooling methods, AdamGNN enhances the node representation with the learned multi-grained semantics and avoids losing node features and graph structure information during pooling. Specifically, a differentiable pooling operator is proposed to adaptively generate a multi-grained structure that involves meso- and macro-level semantic information in the graph. We also devise the unpooling operator and the flyback aggregator in AdamGNN to better leverage the multi-grained semantics to enhance node representations. The updated node representations can further adjust the graph representation in the next iteration. Experiments on 14 real-world graph datasets show that AdamGNN can significantly outperform 17 competing models on both node- and graph-wise tasks. The ablation studies confirm the effectiveness of AdamGNN's components, and the last empirical analysis further reveals the ingenious ability of AdamGNN in capturing long-range interactions.


title: Optimal Task Assignment to Heterogeneous Federated Learning Devices

abstract: Federated Learning provides new opportunities for training machine learning models while respecting data privacy. This technique is based on heterogeneous devices that work together to iteratively train a model while never sharing their own data. Given the synchronous nature of this training, the performance of Federated Learning systems is dictated by the slowest devices, also known as stragglers. In this paper, we investigate the problem of minimizing the duration of Federated Learning rounds by controlling how much data each device uses for training. We formulate this problem as a makespan minimization problem with identical, independent, and atomic tasks that have to be assigned to heterogeneous resources with non-decreasing cost functions while respecting lower and upper limits of tasks per resource. Based on this formulation, we propose a polynomial-time algorithm named OLAR and prove that it provides optimal schedules. We evaluate OLAR in an extensive experimental evaluation using simulation that includes comparisons to other algorithms from the state of the art and new extensions to them. Our results indicate that OLAR provides optimal solutions with a small execution time. They also show that the presence of lower and upper limits of tasks per resource erase any benefits that suboptimal heuristics could provide in terms of algorithm execution time.


title: Asymptotic decomposition of solutions to random parabolic operators with oscillating coefficients

abstract: We consider Cauchy problem for a divergence form second order parabolic operator with rapidly oscillating coefficients that are periodic in spatial variable and random stationary ergodic in time. As was proved in [25] and [13] in this case the homogenized operator is deterministic. We obtain the leading terms of the asymptotic expansion of the solution , these terms being deterministic functions, and show that a properly renormalized difference between the solution and the said leading terms converges to a solution of some SPDE.


title: Spinor wave equation, relativistic condition, and nonlocality of photon spin

abstract: The purpose of this paper is to derive the photon spin and to deduce its properties from a pair of quantum equations for the photon. To this end, Darwin's equations are reinterpreted so as to meet the need of the quantum mechanics of the photon. It is found that the photon wavefunction transforms under Lorentz transformation as a spinor. The relativistic nature of the photon is expressed through a constraint equation on the wavefunction in such a way that the wave equation, which takes on the form of the Schr\"{o}dinger equation, is not Lorentz covariant unless the constraint equation is taken into account. The wave equation predicts the existence of a kind of spin, an intrinsic degree of freedom. But the constraint equation makes the spin nonlocal in the sense that no unique local density exists for the spin in position space. The nonlocality of the photon spin is a reflection of the nonlocality of the photon itself.


title: Influence of fast emissions and statistical de-excitation on the isospin transport ratio

abstract: Isospin transport ratio is a powerful method to estimate the neutron-proton (n-p) equilibration in heavy-ion collisions, and extensively used to obtain information on the asy-stiffness of the nuclear Equation of State. In fact such a ratio is expected to bypass any perturbations introducing a linear transformation of the chosen observable. In particular, it is supposed to overcome contributions due to emission, either of dynamical or statistical nature, from the primary fragments formed during the collisions. In this paper we explore the validity of this assumption, looking at the quasi-projectile n-p ratio ($N/Z$) in peripheral and semi-peripheral events for Ca+Ca reactions at 35\amev{}, simulated via the Antisymmetrized Molecular Dynamics transport model, coupled to different statistical decay codes. The statistical de-excitation of the primary fragments introduces a linear transformation at relatively high excitation energies (above 2\amev{}) when the residue approaches the Evaporation Attractor Line, while some effect is produced at lower excitation energies due to the occurrence of some non-linearities. As for fast emissions after the end of the projectile-target interaction it is shown that they introduce a non-linear transformation too.


title: MLRSNet: A Multi-label High Spatial Resolution Remote Sensing Dataset for Semantic Scene Understanding

abstract: To better understand scene images in the field of remote sensing, multi-label annotation of scene images is necessary. Moreover, to enhance the performance of deep learning models for dealing with semantic scene understanding tasks, it is vital to train them on large-scale annotated data. However, most existing datasets are annotated by a single label, which cannot describe the complex remote sensing images well because scene images might have multiple land cover classes. Few multi-label high spatial resolution remote sensing datasets have been developed to train deep learning models for multi-label based tasks, such as scene classification and image retrieval. To address this issue, in this paper, we construct a multi-label high spatial resolution remote sensing dataset named MLRSNet for semantic scene understanding with deep learning from the overhead perspective. It is composed of high-resolution optical satellite or aerial images. MLRSNet contains a total of 109,161 samples within 46 scene categories, and each image has at least one of 60 predefined labels. We have designed visual recognition tasks, including multi-label based image classification and image retrieval, in which a wide variety of deep learning approaches are evaluated with MLRSNet. The experimental results demonstrate that MLRSNet is a significant benchmark for future research, and it complements the current widely used datasets such as ImageNet, which fills gaps in multi-label image research. Furthermore, we will continue to expand the MLRSNet. MLRSNet and all related materials have been made publicly available at https://data.mendeley.com/datasets/7j9bv9vwsx/2 and https://github.com/cugbrs/MLRSNet.git.


title: A photoionized accretion disk around a young high-mass star

abstract: We present high spatial resolution ($52$ au) observations of the high-mass young stellar object (HMYSO) G345.4938+01.4677 made using the Atacama Large Millimeter/sub-millimeter Array (ALMA). This O-type HMYSO is located at 2.38 kpc and it is associated with a luminosity of $1.5\times10^5 L_\odot$. We detect circumstellar emission from the H$38\beta$ hydrogen recombination line showing a compact structure rotating perpendicularly to the previously detected radio jet. We interpret this emission as tracing a photo-ionized accretion disk around the HMYSO. While this disk-like structure seems currently too small to sustain continued accretion, the data present direct observational evidence of how disks can effectively survive the photo-ionization feedback from young high-mass stars. We also report the detection of a low-mass young stellar object in the vicinity of the HMYSO and suggest that it forms a high-mass and low-mass star binary system.


title: Lattices and the Geometry of Numbers

abstract: In this paper we discuss about properties of lattices and its application in theoretical and algorithmic number theory. This result of Minkowski regarding the lattices initiated the subject of Geometry of Numbers, which uses geometry to study the properties of algebraic numbers. It has application on various other fields of mathematics especially the study of Diophantine equations, analysis of functional analysis etc. This paper will review all the major developments that have occurred in the field of geometry of numbers. In this paper we shall first give a broad overview of the concept of lattice and then discuss about the geometrical properties it has and its applications.


title: CariMe: Unpaired Caricature Generation with Multiple Exaggerations

abstract: Caricature generation aims to translate real photos into caricatures with artistic styles and shape exaggerations while maintaining the identity of the subject. Different from the generic image-to-image translation, drawing a caricature automatically is a more challenging task due to the existence of various spacial deformations. Previous caricature generation methods are obsessed with predicting definite image warping from a given photo while ignoring the intrinsic representation and distribution for exaggerations in caricatures. This limits their ability on diverse exaggeration generation. In this paper, we generalize the caricature generation problem from instance-level warping prediction to distribution-level deformation modeling. Based on this assumption, we present the first exploration for unpaired CARIcature generation with Multiple Exaggerations (CariMe). Technically, we propose a Multi-exaggeration Warper network to learn the distribution-level mapping from photo to facial exaggerations. This makes it possible to generate diverse and reasonable exaggerations from randomly sampled warp codes given one input photo. To better represent the facial exaggeration and produce fine-grained warping, a deformation-field-based warping method is also proposed, which helps us to capture more detailed exaggerations than other point-based warping methods. Experiments and two perceptual studies prove the superiority of our method comparing with other state-of-the-art methods, showing the improvement of our work on caricature generation.


title: WeChat Neural Machine Translation Systems for WMT20

abstract: We participate in the WMT 2020 shared news translation task on Chinese to English. Our system is based on the Transformer (Vaswani et al., 2017a) with effective variants and the DTMT (Meng and Zhang, 2019) architecture. In our experiments, we employ data selection, several synthetic data generation approaches (i.e., back-translation, knowledge distillation, and iterative in-domain knowledge transfer), advanced finetuning approaches and self-bleu based model ensemble. Our constrained Chinese to English system achieves 36.9 case-sensitive BLEU score, which is the highest among all submissions.


title: Cosmological Vlasov-Poisson Simulations of Structure Formation with Relic Neutrinos: Nonlinear Clustering and the Neutrino Mass

abstract: We present the results of cosmological simulations of large-scale structure formation with massive neutrinos. The phase-space distribution of the cosmic relic neutrinos is followed, for the first time, by directly integrating the six-dimensional Vlasov-Poisson equations. Our novel approach allows us to represent free streaming and clustering of neutrinos, and their gravitational interaction with cold dark matter accurately. We thus obtain solutions for the collisionless dynamics independent of conventional N-body methods. We perform a suite of hybrid N-body/Vlasov simulations with varying the neutrino mass, and systematically examine the dynamical effects of massive neutrinos on large-scale structure formation. Our simulations show characteristic large-scale clustering of the neutrinos and their coherent streaming motions relative to dark matter. The effective local neutrino "temperature" around massive galaxy clusters varies by several percent with respect to the cosmic mean; the neutrinos in clusters can be hotter or colder depending on the neutrino mass. We study a number of statistics of the large-scale structure and of dark matter halos in comparison with those obtained by N-body simulations and/or by perturbation theory. Our simulations mark an important milestone in numerical cosmology, and pave a new way to study cosmic structure formation with massive neutrinos.


title: Green function formalism for resonant interaction of x-rays with nuclei in structured media

abstract: The resonant interaction between x-ray photons and nuclei is one of the most exciting subjects of the burgeoning field of x-ray quantum optics. A resourceful platform used so far are thin-film x-ray cavities with embedded layers or M\"ossbauer nuclei such as $^{57}\mathrm{Fe}$. A new quantum optical model based on the classical electromagnetic Green's function is developed to investigate theoretically the nuclear response inside the x-ray cavity. The model is versatile and provides an intuitive picture about the influence of the cavity structure on the resulting spectra. We test its predictive powers with the help of the semiclassical coherent scattering formalism simulations and discuss our results for increasing complexity of layer structures.


title: Muckenhoupt-type conditions on weighted Morrey spaces

abstract: We define a Muckenhoup-type condition on weighted Morrey spaces using the K\"othe dual of the space. We show that the condition is necessary and sufficient for the boundedness of the maximal operator defined with balls centered at the origin on weighted Morrey spaces. A modified condition characterizes the weighted inequalities for the Calder\'on operator. We also show that the Muckenhoup-type condition is necessary and sufficient for the boundedness on weighted local Morrey spaces of the usual Hardy-Littlewood maximal operator, simplifying the previous characterization of Nakamura-Sawano-Tanaka. For the same operator, in the case of global Morrey spaces the condition is necessary and for the sufficiency we add a local $A_p$ condition. We can extrapolate from Lebesgue $A_p$-weighted inequalities to weighted global and local Morrey spaces in a very general setting, with applications to many operators.


title: Magnetic quadri-dipolar stars rotating in vacuum

abstract: Main sequence stars and compact objects like white dwarfs and neutron stars are usually embedded in magnetic fields that strongly deviate from a pure dipole located right at the stellar centre. An off-centred dipole can sometimes better adjust existing data and offer a simple geometric picture to include multipolar fields. However, such configurations are usually to restrictive, limiting multipolar components to strength less than the underlying dipole. In this paper, we consider the most general lowest order multipolar combination given by a dipole and a quadrupole magnetic field association in vacuum. Following the general formalism for multipolar field computations, we derive the full electromagnetic field outside a rotating quadridipole. Exact analytical expressions for the Poynting flux and the electromagnetic kick are given. Such geometry is useful to study the magnetosphere of neutron stars for which more and more compelling observations reveals hints for at least quadridipolar fields. We also show that for sufficiently high quadrupole components at the stellar surface, the electromagnetic kick imprinted to a neutron star can reach thousands of km/s for a millisecond period at birth.


title: Microbunching Instability Characterisation via Temporally Modulated Laser Pulses

abstract: High-brightness electron bunches, such as those generated and accelerated in free-electron lasers (FELs), can develop small-scale structure in the longitudinal phase space. This causes variations in the slice energy spread and current profile of the bunch which then undergo amplification, in an effect known as the microbunching instability. By imposing energy spread modulations on the bunch in the low-energy section of an accelerator, using an undulator and a modulated laser pulse in the centre of a dispersive chicane, it is possible tomanipulate the bunch longitudinal phase space. This allows for the control and study of the instability in unprecedented detail. We report measurements and analysis of such modulated electron bunches in the 2Dspectro-temporal domain at the FERMI FEL, for three different bunch compression schemes. We also perform corresponding simulations of these experiments and show that the codes are indeed able to reproduce the measurements across a wide spectral range. This detailed experimental verification of the ability of codes to capture the essential beam dynamics of the microbunching instability will benefit the design and performance of future FELs.


title: An Efficient Numerical Method for Forward-Backward Stochastic Differential Equations Driven by $G$-Brownian motion

abstract: In this paper, we study the numerical method for solving forward-backward stochastic differential equations driven by $G$-Brownian motion ($G$-FBSDEs) which correspond to fully nonlinear partial differential equations (PDEs). First, we give an approximate conditional $G$-expectation and obtain feasible methods to calculate the distribution of $G$-Brownian motion. On this basis, some efficient numerical schemes for $G$-FBSDEs are then proposed. We rigorously analyze errors of the proposed schemes and prove the convergence results. Finally, several numerical experiments are given to demonstrate the accuracy of our method.


title: Filling the Black Hole Mass Gap: Avoiding Pair Instability in Massive Stars through Addition of Non-Nuclear Energy

abstract: In standard stellar evolution, stars with masses ranging from approximately $150$ to $240 M_\odot$ are expected to evolve to a pair instability supernova with no black hole (BH) remnant. This evolutionary behavior leads to a predicted gap in the black hole mass function from approximately $50$ to $140 M_\odot$. Yet the LIGO and Virgo Collaborations[1] recently discovered black holes of masses $66 M_\odot$ and $85 M_\odot$ in the gravitational wave event GW190521. We propose a new method to populate the BH mass gap. If an energy source is added throughout the star in addition to nuclear fusion, it is possible for the altered evolution to avoid the complete destruction of a pair instability supernova, and instead a BH remnant is left behind. An example of an extra energy source is dark matter annihilation within the star, but our results hold more generally. We show this phenomenon by exploring the effect of adding an energy source independent of temperature and density to a $180 M_\odot$ star, using the MESA one-dimensional stellar evolution software. If $\sim50$\% of the star's energy is due to this new source, the star is capable of avoiding the pair instability entirely and evolving towards a core-collapse supernova and ultimately a BH remnant with mass $\sim 120 M_\odot$.


title: Efficient Construction of a Control Modular Adder on a Carry-Lookahead Adder Using Relative-phase Toffoli Gates

abstract: Control modular addition is a core arithmetic function, and we must consider the computational cost for actual quantum computers to realize efficient implementation. To achieve a low computational cost in a control modular adder, we focus on minimizing KQ, defined by the product of the number of qubits and the depth of the circuit. In this paper, we construct an efficient control modular adder with small KQ by using relative-phase Toffoli gates in two major types of quantum computers: Fault-Tolerant Quantum Computers (FTQ) on the Logical layer and Noisy Intermediate-Scale Quantum Computers (NISQ). We give a more efficient construction compared to Van Meter and Itoh's, based on a carry-lookahead adder. In FTQ, $T$ gates incur heavy cost due to distillation, which fabricates ancilla for running $T$ gates with high accuracy but consumes a lot of specially prepared ancilla qubits and a lot of time. Thus, we must reduce the number of $T$ gates. We propose a new control modular adder that uses only 20% of the number of $T$ gates of the original. Moreover, when we take distillation into consideration, we find that we minimize $\text{KQ}_{T}$ (the product of the number of qubits and $T$-depth) by running $\Theta\left(n / \sqrt{\log n} \right)$ $T$ gates simultaneously. In NISQ, CNOT gates are the major error source. We propose a new control modular adder that uses only 35% of the number of CNOT gates of the original. Moreover, we show that the $\text{KQ}_{\text{CX}}$ (the product of the number of qubits and CNOT-depth) of our circuit is 38% of the original. Thus, we realize an efficient control modular adder, improving prospects for the efficient execution of arithmetic in quantum computers.


title: Electronic friction coefficients from the atom-in-jellium model for $Z=1-92$

abstract: The break-down of the Born-Oppenheimer approximation is an important topic in chemical dynamics on metal surfaces. In this context, the most frequently used "work-horse" is electronic friction theory, commonly relying on friction coefficients obtained from density functional theory (DFT) calculations from the early 80s based on the atom-in-jellium model. However, results are only available for a limited set of jellium densities and elements ($Z=1-18$). In this work, these calculations are revisited by investigating the corresponding friction coefficients for the entire periodic table ($Z=1-92$). Furthermore, friction coefficients obtained by including the electron density gradient on the Generalized Gradient Approximation (GGA) level are presented. Finally, we show that spin polarization and relativistic effects can have sizeable effects on these friction coefficients for some elements.

categories: physics.chem-ph cond-mat.mtrl-sci cond-mat.other physics.plasm-ph


title: Scipp: Scientific data handling with labeled multi-dimensional arrays for C++ and Python

abstract: Scipp is heavily inspired by the Python library xarray. It enriches raw NumPy-like multi-dimensional arrays of data by adding named dimensions and associated coordinates. Multiple arrays are combined into datasets. On top of this, scipp introduces (i) implicit handling of physical units, (ii) implicit propagation of uncertainties, (iii) support for histograms, i.e., bin-edge coordinate axes, which exceed the data's dimension extent by one, and (iv) support for event data. In conjunction these features enable a more natural and more concise user experience. The combination of named dimensions, coordinates, and units helps to drastically reduce the risk for programming errors. The core of scipp is written in C++ to open opportunities for performance improvements that a Python-based solution would not allow for. On top of the C++ core, scipp's Python components provide functionality for plotting and content representations, e.g., for use in Jupyter Notebooks. While none of scipp's concepts in isolation is novel per-se, we are not aware of any project combining all of these aspects in a single coherent software package.


title: Predicting the flow field in a U-bend with deep neural networks

abstract: This paper describes a study based on computational fluid dynamics (CFD) and deep neural networks that focusing on predicting the flow field in differently distorted U-shaped pipes. The main motivation of this work was to get an insight about the justification of the deep learning paradigm in hydrodynamic hull optimisation processes that heavily depend on computing turbulent flow fields and that could be accelerated with models like the one presented. The speed-up can be even several orders of magnitude by surrogating the CFD model with a deep convolutional neural network. An automated geometry creation and evaluation process was set up to generate differently shaped two-dimensional U-bends and to carry out CFD simulation on them. This process resulted in a database with different geometries and the corresponding flow fields (2-dimensional velocity distribution), both represented on 128x128 equidistant grids. This database was used to train an encoder-decoder style deep convolutional neural network to predict the velocity distribution from the geometry. The effect of two different representations of the geometry (binary image and signed distance function) on the predictions was examined, both models gave acceptable predictions with a speed-up of two orders of magnitude.


title: Experimental certification of nonclassicality via phase-space inequalities

abstract: In spite of its fundamental importance in quantum science and technology, the experimental certification of nonclassicality is still a challenging task, especially in realistic scenarios where losses and noise imbue the system. Here, we present the first experimental implementation of the recently introduced phase-space inequalities for nonclassicality certification, which conceptually unite phase-space representations with correlation conditions. We demonstrate the practicality and sensitivity of this approach by studying nonclassicality of a family of noisy and lossy quantum states of light. To this end, we experimentally generate single-photon-added thermal states with various thermal mean photon numbers and detect them at different loss levels. Based on the reconstructed Wigner and Husimi Q functions, the inequality conditions detect nonclassicality despite the fact that the involved distributions are nonnegative, which includes cases of high losses (93%) and cases where other established methods do not reveal nonclassicality. We show the advantages of the implemented approach and discuss possible extensions that assure a wide applicability for quantum science and technologies.


title: On a Brownian motion conditioned to stay in an open set

abstract: Distribution of a Brownian motion conditioned to start from the boundary of an open set $G$ and to stay in $G$ for a finite period of time is studied. Characterizations of such distributions in terms of certain singular stochastic differential equations are obtained. Results are applied to the study of boundaries of clusters in some coalescing stochastic flows on $\mathbb{R}.$


title: NodeSig: Binary Node Embeddings via Random Walk Diffusion

abstract: Graph Representation Learning (GRL) has become a key paradigm in network analysis, with a plethora of interdisciplinary applications. As the scale of networks increases, most of the widely used learning-based graph representation models also face computational challenges. While there is a recent effort toward designing algorithms that solely deal with scalability issues, most of them behave poorly in terms of accuracy on downstream tasks. In this paper, we aim to study models that balance the trade-off between efficiency and accuracy. In particular, we propose NodeSig, a scalable model that computes binary node representations. NodeSig exploits random walk diffusion probabilities via stable random projections towards efficiently computing embeddings in the Hamming space. Our extensive experimental evaluation on various networks has demonstrated that the proposed model achieves a good balance between accuracy and efficiency compared to well-known baseline models on the node classification and link prediction tasks.


title: Active Inference or Control as Inference? A Unifying View

abstract: Active inference (AI) is a persuasive theoretical framework from computational neuroscience that seeks to describe action and perception as inference-based computation. However, this framework has yet to provide practical sensorimotor control algorithms that are competitive with alternative approaches. In this work, we frame active inference through the lens of control as inference (CaI), a body of work that presents trajectory optimization as inference. From the wider view of `probabilistic numerics', CaI offers principled, numerically robust optimal control solvers that provide uncertainty quantification, and can scale to nonlinear problems with approximate inference. We show that AI may be framed as partially-observed CaI when the cost function is defined specifically in the observation states.


title: RefVOS: A Closer Look at Referring Expressions for Video Object Segmentation

abstract: The task of video object segmentation with referring expressions (language-guided VOS) is to, given a linguistic phrase and a video, generate binary masks for the object to which the phrase refers. Our work argues that existing benchmarks used for this task are mainly composed of trivial cases, in which referents can be identified with simple phrases. Our analysis relies on a new categorization of the phrases in the DAVIS-2017 and Actor-Action datasets into trivial and non-trivial REs, with the non-trivial REs annotated with seven RE semantic categories. We leverage this data to analyze the results of RefVOS, a novel neural network that obtains competitive results for the task of language-guided image segmentation and state of the art results for language-guided VOS. Our study indicates that the major challenges for the task are related to understanding motion and static actions.


title: Twisted and Singular gravitating vortices

abstract: We introduce the notion of twisted gravitating vortex on a compact Riemann surface. If the genus of the Riemann surface is greater than 1 and the twisting forms have suitable signs, we prove an existence and uniqueness result for suitable range of the coupling constant generalizing the result of arXiv:1510.03810v2 in the non twisted setting. It is proved via solving a continuity path deforming the coupling constant from 0 for which the system decouples as twisted K\"ahler-Einstein metric and twisted vortices. Moreover, specializing to a family of twisting forms smoothing delta distribution terms, we prove the existence of singular gravitating vortices whose K\"ahler metric has conical singularities and Hermitian metric has parabolic singularities. In the Bogomol'nyi phase, we establish an existence result for singular Einstein-Bogomol'nyi equations, which represents cosmic strings with singularities.


title: Review and Analysis of Three Components of Differential Evolution Mutation Operator in MOEA/D-DE

abstract: A decomposition-based multi-objective evolutionary algorithm with a differential evolution variation operator (MOEA/D-DE) shows high performance on challenging multi-objective problems (MOPs). The DE mutation consists of three key components: a mutation strategy, an index selection method for parent individuals, and a bound-handling method. However, the configuration of the DE mutation operator that should be used for MOEA/D-DE has not been thoroughly investigated in the literature. This configuration choice confuses researchers and users of MOEA/D-DE. To address this issue, we present a review of the existing configurations of the DE mutation operator in MOEA/D-DE and systematically examine the influence of each component on the performance of MOEA/D-DE. Our review reveals that the configuration of the DE mutation operator differs depending on the source code of MOEA/D-DE. In our analysis, a total of 30 configurations (three index selection methods, two mutation strategies, and five bound handling methods) are investigated on 16 MOPs with up to five objectives. Results show that each component significantly affects the performance of MOEA/D-DE. We also present the most suitable configuration of the DE mutation operator, which maximizes the effectiveness of MOEA/D-DE.


title: Comparison of the $n$-categorical nerves

abstract: Our aim is to compare three nerve functors for strict $n$-categories: the Street nerve, the cellular nerve and the multi-simplicial nerve. We show that these three functors are equivalent in some appropriate sense. In particular, the classes of $n$-categorical weak equivalences that they define coincide: they are the Thomason equivalences. We give two applications of this result: the first one states that a Dyer-Kan-type equivalence for Thomason equivalences is a Thomason equivalence; the second one, fundamental, is the stability of the class of Thomason equivalences under the dualities of the category of strict $n$-categories.


title: Dumbbell impurities in 2D crystals of repulsive colloidal spheres induce particle-bound dislocations

abstract: Impurity-induced defects play a crucial role for the properties of crystals, but little is known about impurities with anisotropic shape. Here, we study how colloidal dumbbells distort and interact with a hexagonal crystal of charged colloidal spheres at a fluid interface. We find that subtle changes in the dumbbell length induce a transition from a local distortion to a particle-bound dislocation, and determine how the dumbbell moves inside the repulsive hexagonal lattice. Our results provide new routes towards controlling material properties through particle-bound dislocations.


title: Encrypted control for networked systems -- An illustrative introduction and current challenges

abstract: Cloud computing and distributed computing are becoming ubiquitous in many modern control systems such as smart grids, building automation, robot swarms or intelligent transportation systems. Compared to "isolated" control systems, the advantages of cloud-based and distributed control systems are, in particular, resource pooling and outsourcing, rapid scalability, and high performance. However, these capabilities do not come without risks. In fact, the involved communication and processing of sensitive data via public networks and on third-party platforms promote, among other cyberthreats, eavesdropping and manipulation of data. Encrypted control addresses this security gap and provides confidentiality of the processed data in the entire control loop. This paper presents a tutorial-style introduction to this young but emerging field in the framework of secure control for networked dynamical systems.


title: Robin double-phase problems with singular and superlinear terms

abstract: We consider a nonlinear Robin problem driven by the sum of $p$-Laplacian and $q$-Laplacian (i.e. the $(p,q)$-equation). In the reaction there are competing effects of a singular term and a parametric perturbation $\lambda f(z,x)$, which is Carath\'eodory and $(p-1)$-superlinear at $x\in\mathbb{R},$ without satisfying the Ambrosetti-Rabinowitz condition. Using variational tools, together with truncation and comparison techniques, we prove a bifurcation-type result describing the changes in the set of positive solutions as the parameter $\lambda>0$ varies.


title: Local invariants of braiding quantum gates -- associated link polynomials and entangling power

abstract: For a generic $n$-qubit system, local invariants under the action of $SL(2,\mathbb{C})^{\otimes n}$ characterize non-local properties of entanglement. In general, such properties are not immediately apparent and hard to construct. Here we consider certain two-qubit Yang-Baxter operators, which we dub of the `X-type', and show that their eigenvalues completely determine the non-local properties of the system. Moreover, we apply the Turaev procedure to these operators and obtain their associated link/knot polynomials. We also compute their entangling power and compare it with that of a generic two-qubit operator.


title: Kernel Two-Sample and Independence Tests for Non-Stationary Random Processes

abstract: Two-sample and independence tests with the kernel-based MMD and HSIC have shown remarkable results on i.i.d. data and stationary random processes. However, these statistics are not directly applicable to non-stationary random processes, a prevalent form of data in many scientific disciplines. In this work, we extend the application of MMD and HSIC to non-stationary settings by assuming access to independent realisations of the underlying random process. These realisations - in the form of non-stationary time-series measured on the same temporal grid - can then be viewed as i.i.d. samples from a multivariate probability distribution, to which MMD and HSIC can be applied. We further show how to choose suitable kernels over these high-dimensional spaces by maximising the estimated test power with respect to the kernel hyper-parameters. In experiments on synthetic data, we demonstrate superior performance of our proposed approaches in terms of test power when compared to current state-of-the-art functional or multivariate two-sample and independence tests. Finally, we employ our methods on a real socio-economic dataset as an example application.


title: Kalkayotl: A cluster distance inference code

abstract: Context: Stellar clusters are benchmarks for theories of star formation and evolution. The high precision parallax data of the Gaia mission allows significant improvements in the distance determination to stellar clusters and its stars. In order to have accurate and precise distance determinations, systematics like the parallax spatial correlations need to be accounted for, especially for stars in small sky regions. Aims: Provide the astrophysical community with a free and open code designed to simultaneously infer cluster parameters (i.e. distance and size) and the distances to its stars using Gaia parallax measurements. It includes cluster oriented prior families and is specifically designed to deal with the Gaia parallax spatial correlations. Methods: A Bayesian hierarchical model is created to allow the inference of both the cluster parameters and distances to its stars. Results: Using synthetic data that mimics Gaia parallax uncertainties and spatial correlations, we observe that our cluster oriented prior families result in distance estimates with smaller errors than those obtained with an exponentially decreasing space density prior. In addition, the treatment of the parallax spatial correlations minimizes errors in the estimated cluster size and stellar distances and avoids the underestimation of uncertainties. Although neglecting the parallax spatial correlations has no impact on the accuracy of cluster distance determinations, it underestimates the uncertainties and may result in measurements that are incompatible with the true value. Conclusions: The combination of prior knowledge with the treatment of Gaia parallax spatial correlations produces accurate (error <10%) and trustworthy estimates (i.e. true values contained within the 2$\sigma$ uncertainties) of clusters distances for clusters up to ~5 kpc, and cluster sizes for clusters up to ~1 kpc.


title: Complexity and algorithms for constant diameter augmentation problems

abstract: We study the following problem: for given integers $d,k$ and graph $G$, can we obtain a graph with diameter $d$ via at most $k$ edge deletions ? We determine the computational complexity of this and related problems for different values of $d$.


title: A Molecular Dynamics Investigation of Mechanical Properties of Graphene Reinforced Iron Composite and The Effect of Vacancy Defect Distance from the Matrix-Fiber Interface

abstract: Graphene is a material of excellent mechanical properties, which make it an ideal fiber for reinforcing metal. Since iron is the most used metal in the world, reinforcing iron with graphene can reduce the overall requirement of material in any application where strength is demanded. However, the effect of graphene reinforcement on the mechanical properties of iron needs to be known before the industrial application of the composite. In this paper, we have investigated the mechanical properties of graphene-reinforced iron composite by Molecular Dynamics (MD) method for various conditions. The properties were investigated by applying uniaxial tension on a modeled representative volume element (RVE). The effect of temperature on the mechanical property of the composite was also studied because the knowledge is required for manufacturing products with the composite operating at a wide temperature range. MD analysis also revealed that the initiation of fracture is from the matrix-fiber interface. We also investigated how the distance of vacancy defects from the matrix-fiber interface affects the mechanical properties of the composite, which can be used to select a suitable manufacturing process. The results obtained from this study show that vacancy defects lower the strength at a greater extent as it gets closer to the interface.


title: Molecular Dynamics Simulation for the Analysis of Mechanical Properties and Effect of Stone-Wales and BiVacancy Defect on Carbon Nanotube Reinforced Iron Composites

abstract: Carbon nanotube (CNT) reinforced metal matrix composites (MMCs) are gaining the attention of the researchers because of their demand in space and automobile industries for having low weight and high mechanical properties. Iron is the most used metal in all engineering fields. Therefore, reinforcing iron with CNT can reduce its required amount, which might have a positive economic impact due to the reduced cost of production. However, before the industrial application of any material the mechanical properties under different conditions must be known. In this study, the mechanical properties of iron reinforced separately with single, double and triple wall CNTs are investigated by Molecular Dynamics (MD) simulation. The study revealed that the strength and stiffness of pure iron could be enhanced up to 80.4 % and 57.4 %, respectively, by adding CNTs into iron. We also investigated the effect of fiber volume percentage and temperature on the mechanical properties of the composite having single, double and triplewalled carbon nanotubes individually. As the stone-wales and bi-vacancy defects are inherently introduced in CNTs during manufacturing, their effect on mechanical properties are also investigated in the present study


title: Palladium Zero-Mode Waveguides for Optical Single Molecule Detection with Nanopores

abstract: Holes in metal films block any transmitting light if the wavelength is much larger than the hole diameter, establishing such nanopores as so-called Zero Mode Waveguides (ZMWs). Molecules on the other hand, can still passage through these holes. We use this to detect individual fluorophore-labelled molecules as they travel through a ZMW and thereby traverse from the dark region to the illuminated side, upon which they emit fluorescent light. This is beneficial both for background suppression and to prevent premature bleaching. We use palladium as a novel metal-film material for ZMWs, which is advantageous compared to conventionally used metals. We demonstrate that it is possible to simultaneously detect translocations of individual free fluorophores of different colors. Labeled DNA and protein biomolecules can be detected as well at the single-molecule level with a high signal-to-noise ratio and at high bandwidth, which opens the door to a variety of single-molecule biophysics studies.


title: Jordan Algebras of Symmetric Matrices

abstract: We study linear spaces of symmetric matrices whose reciprocal is also a linear space. These are Jordan algebras. We classify such algebras in low dimensions, and we study the associated Jordan loci in the Grassmannian.


title: McStas (ii): An overview of components, their use, and advice for user contributions

abstract: An key element of the success of McStas is the component layer where users and developers alike are contributing to the description of new physical models and features. In McStas, components realise all physical elements of the simulated instrument from source via optics and samples to detector. In this second review paper of the McStas package, we present an overview of the component classes in McStas: sources, monitors, optics, samples, misc, and contrib. Within each component class we give thorough examples of high-quality components, including their algorithms and example use. We present two example instruments, one for a continuous source and one for a time-of-flight source, that together demonstrate the use of the main component classes. Finally, we give tips and instructions that will allow the reader to write good components and elucidate the pathway of contributing new components to McStas.


title: Control of quantum noise: on the role of dilations

abstract: We show that every finite-dimensional quantum system with Markovian (i.e., GKLS-generated) time evolution has an autonomous unitary dilation which can be dynamically decoupled. Since there is also always an autonomous unitary dilation which cannot be dynamically decoupled, this highlights the role of dilations in the control of quantum noise. We construct our dilation via a time-dependent version of Stinespring in combination with Howland's clock Hamiltonian and certain point-localised states, which may be regarded as a C*-algebraic analogue of improper bra-ket position eigenstates and which are hence of independent mathematical and physical interest.


title: On various types of density of numerical radius attaining operators

abstract: In this paper, we are interested in studying two properties related to the denseness of the operators which attain their numerical radius: the Bishop-Phelps-Bollob\'as point and operator properties for numerical radius (BPBpp-nu and BPBop-nu, respectively). We prove that every Banach space with micro-transitive norm and second numerical index strictly positive satisfy the BPBpp-nu and that, if the numerical index of $X$ is 1, only one-dimensional spaces enjoy it. On the other hand, we show that the BPBop-nu is a very restrictive property: under some general assumptions, it holds only for one-dimensional spaces. We also consider two weaker properties, the local versions of BPBpp-nu and BPBop-nu, where the $\eta$ which appears in their definition does not depend just on $\epsilon > 0$ but also on a state $(x, x^*)$ or on a numerical radius one operator $T$. We address the relation between the local BPBpp-nu and the strong subdifferentiability of the norm of the space $X$. We show that finite dimensional spaces and $c_0$ are examples of Banach spaces satisfying the local BPBpp-nu, and we exhibit an example of a Banach space with strongly subdifferentiable norm failing it. We finish the paper by showing that finite dimensional spaces satisfy the local BPBop-nu and that, if $X$ has strictly positive numerical index and has the approximation property, this property is equivalent to finite dimensionality.


title: Energy conditions in $f(Q,T)$ gravity

abstract: The recently proposed $f(Q, T)$ gravity (Xu et al. Eur. Phys. J. C \textbf{79} (2019) 708) is an extension of the symmetric teleparallel gravity. The gravitational action $L$ is given by an arbitrary function $f$ of the non-metricity $Q$ and the trace of the matter-energy momentum tensor $T$. In this paper, we examined the essence of some well prompted forms of $f(Q,T)$ gravity models i.e. $f(Q,T)= mQ+bT$ and $f(Q,T)= m Q^{n+1}+b T$ where $m$, $b$, and $n$ are model parameters. We have used the proposed deceleration parameter, which predicts both decelerated and accelerated phases of the Universe, with the transition redshift by recent observations and obtains energy density ($\rho$) and pressure ($p$) to study the various energy conditions for cosmological models. The equation of state parameter ($\omega\simeq -1$) in the present model also supports the accelerating behavior of the Universe. In both, the models, the null, weak, and dominant energy conditions are obeyed with violating strong energy conditions as per the present accelerated expansion.


title: Probabilistic Programs with Stochastic Conditioning

abstract: We tackle the problem of conditioning probabilistic programs on distributions of observable variables. Probabilistic programs are usually conditioned on samples from the joint data distribution, which we refer to as deterministic conditioning. However, in many real-life scenarios, the observations are given as marginal distributions, summary statistics, or samplers. Conventional probabilistic programming systems lack adequate means for modeling and inference in such scenarios. We propose a generalization of deterministic conditioning to stochastic conditioning, that is, conditioning on the marginal distribution of a variable taking a particular form. To this end, we first define the formal notion of stochastic conditioning and discuss its key properties. We then show how to perform inference in the presence of stochastic conditioning. We demonstrate potential usage of stochastic conditioning on several case studies which involve various kinds of stochastic conditioning and are difficult to solve otherwise. Although we present stochastic conditioning in the context of probabilistic programming, our formalization is general and applicable to other settings.


title: Modelling the earth's geomagnetic environment on Cray machines using PETSc and SLEPc

abstract: The British Geological Survey's global geomagnetic model, Model of the Earth's Magnetic Environment (MEME), is an important tool for calculating the earth's magnetic field, which is continually in flux. Whilst the ability to collect data from ground based observation sites and satellites has grown, the memory bound nature of the code has proved a limitation in modelling problem sizes required by modern science. In this paper we describe work replacing the bespoke, sequential, eigen-solver with that of the SLEPc package for solving the system of normal equations. This work had a dual purpose, to break through the memory limit of the code, and thus support the modelling of much larger systems, by supporting execution on distributed machines, and to improve performance. But when adopting SLEPc it was not just the solving of the normal equations, but also fundamentally how we build and distribute the data structures. We describe an approach for building symmetric matrices in a way that provides good load balance and avoids the need for close co-ordination between processes or replication of work. We also study the memory bound nature of the code and combine detailed profiling with software cache prefetching to significantly optimise. Performance and scaling are explored on ARCHER, a Cray XC30, where we achieve a speed up for the solver of 294 times by replacing the model's bespoke approach with SLEPc. This work also provided the ability to model much larger system sizes, up to 100,000 model coefficients. Some of the challenges of modelling systems of this large scale are explored, and mitigations including hybrid MPI+OpenMP along with the use of iterative solvers are also considered. The result of this work is a modern MEME model that is not only capable of simulating problem sizes demanded by state of the art geomagnetism but also acts as further evidence to the utility of the SLEPc libary.


title: Bayesian Policy Search for Stochastic Domains

abstract: AI planning can be cast as inference in probabilistic models, and probabilistic programming was shown to be capable of policy search in partially observable domains. Prior work introduces policy search through Markov chain Monte Carlo in deterministic domains, as well as adapts black-box variational inference to stochastic domains, however not in the strictly Bayesian sense. In this work, we cast policy search in stochastic domains as a Bayesian inference problem and provide a scheme for encoding such problems as nested probabilistic programs. We argue that probabilistic programs for policy search in stochastic domains should involve nested conditioning, and provide an adaption of Lightweight Metropolis-Hastings (LMH) for robust inference in such programs. We apply the proposed scheme to stochastic domains and show that policies of similar quality are learned, despite a simpler and more general inference algorithm. We believe that the proposed variant of LMH is novel and applicable to a wider class of probabilistic programs with nested conditioning.


title: Model order reduction of flow based on a modular geometrical approximation of blood vessels

abstract: We are interested in a reduced order method for the efficient simulation of blood flow in arteries. The blood dynamics is modeled by means of the incompressible Navier-Stokes equations. Our algorithm is based on an approximated domain-decomposition of the target geometry into a number of subdomains obtained from the parametrized deformation of geometrical building blocks (e.g. straight tubes and model bifurcations). On each of these building blocks, we build a set of spectral functions by proper orthogonal decomposition of a large number of snapshots of finite element solutions (offline phase). The global solution of the Navier-Stokes equations on a target geometry is then found by coupling linear combinations of these local basis functions by means of spectral Lagrange multipliers (online phase). Being that the number of reduced degrees of freedom is considerably smaller than their finite element counterpart, this approach allows us to significantly decrease the size of the linear system to be solved in each iteration of the Newton-Raphson algorithm. We achieve large speedups with respect to the full order simulation (in our numerical experiments, the gain is at least of one order of magnitude and grows inversely with respect to the reduced basis size), whilst still retaining satisfactory accuracy for most cardiovascular simulations.


title: A package for computations with sparse resultants

abstract: We introduce the Macaulay2 package SparseResultants, which provides general tools for computing sparse resultants, sparse discriminants, and hyperdeterminants. We give some background on the theory and briefly show how the package works.


title: Joint Persian Word Segmentation Correction and Zero-Width Non-Joiner Recognition Using BERT

abstract: Words are properly segmented in the Persian writing system; in practice, however, these writing rules are often neglected, resulting in single words being written disjointedly and multiple words written without any white spaces between them. This paper addresses the problems of word segmentation and zero-width non-joiner (ZWNJ) recognition in Persian, which we approach jointly as a sequence labeling problem. We achieved a macro-averaged F1-score of 92.40% on a carefully collected corpus of 500 sentences with a high level of difficulty.


title: Action Units Recognition by Pairwise Deep Architecture

abstract: In this paper, we propose a new automatic Action Units (AUs) recognition method used in a competition, Affective Behavior Analysis in-the-wild (ABAW). Our method tackles a problem of AUs label inconsistency among subjects by using pairwise deep architecture. While the baseline score is 0.31, our method achieved 0.67 in validation dataset of the competition.


title: Weighing up the new kid on the block: Impressions of using Vitis for HPC software development

abstract: The use of reconfigurable computing, and FPGAs in particular, has strong potential in the field of High Performance Computing (HPC). However the traditionally high barrier to entry when it comes to programming this technology has, until now, precluded widespread adoption. To popularise reconfigurable computing with communities such as HPC, Xilinx have recently released the first version of Vitis, a platform aimed at making the programming of FPGAs much more a question of software development rather than hardware design. However a key question is how well this technology fulfils the aim, and whether the tooling is mature enough such that software developers using FPGAs to accelerate their codes is now a more realistic proposition, or whether it simply increases the convenience for existing experts. To examine this question we use the Himeno benchmark as a vehicle for exploring the Vitis platform for building, executing and optimising HPC codes, describing the different steps and potential pitfalls of the technology. The outcome of this exploration is a demonstration that, whilst Vitis is an excellent step forwards and significantly lowers the barrier to entry in developing codes for FPGAs, it is not a silver bullet and an underlying understanding of dataflow style algorithmic design and appreciation of the architecture is still key to obtaining good performance on reconfigurable architectures.


title: The geometrically m-step solvable Grothendieck conjecture for genus 0 curves over finitely generated fields

abstract: In this paper, we present some partial results for the geometrically m-step solvable Grothendieck conjecture in anabelian geometry. Among other things, we prove the geometrically 3-step solvable Grothendieck conjecture for genus 0 curves over fields finitely generated over the prime field of arbitrary characteristic.


title: Cost-Sensitive Regularization for Diabetic Retinopathy Grading from Eye Fundus Images

abstract: Assessing the degree of disease severity in biomedical images is a task similar to standard classification but constrained by an underlying structure in the label space. Such a structure reflects the monotonic relationship between different disease grades. In this paper, we propose a straightforward approach to enforce this constraint for the task of predicting Diabetic Retinopathy (DR) severity from eye fundus images based on the well-known notion of Cost-Sensitive classification. We expand standard classification losses with an extra term that acts as a regularizer, imposing greater penalties on predicted grades when they are farther away from the true grade associated to a particular image. Furthermore, we show how to adapt our method to the modelling of label noise in each of the sub-problems associated to DR grading, an approach we refer to as Atomic Sub-Task modeling. This yields models that can implicitly take into account the inherent noise present in DR grade annotations. Our experimental analysis on several public datasets reveals that, when a standard Convolutional Neural Network is trained using this simple strategy, improvements of 3-5\% of quadratic-weighted kappa scores can be achieved at a negligible computational cost. Code to reproduce our results is released at https://github.com/agaldran/cost_sensitive_loss_classification.


title: Open-Set Hypothesis Transfer with Semantic Consistency

abstract: Unsupervised open-set domain adaptation (UODA) is a realistic problem where unlabeled target data contain unknown classes. Prior methods rely on the coexistence of both source and target domain data to perform domain alignment, which greatly limits their applications when source domain data are restricted due to privacy concerns. This paper addresses the challenging hypothesis transfer setting for UODA, where data from source domain are no longer available during adaptation on target domain. We introduce a method that focuses on the semantic consistency under transformation of target data, which is rarely appreciated by previous domain adaptation methods. Specifically, our model first discovers confident predictions and performs classification with pseudo-labels. Then we enforce the model to output consistent and definite predictions on semantically similar inputs. As a result, unlabeled data can be classified into discriminative classes coincided with either source classes or unknown classes. Experimental results show that our model outperforms state-of-the-art methods on UODA benchmarks.


title: Local chemical bonding and structural properties in Ti3AlC2 MAX phase and Ti3C2Tx MXene probed by Ti 1s X-ray absorption spectroscopy

abstract: The chemical bonding within the transition-metal carbide materials MAX phase Ti3AlC2 and MXene Ti3C2Tx is investigated by X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) spectroscopies. MAX phases are inherently nanolaminated materials that consist of alternating layers of Mn+1Xn and monolayers of an A-element from the IIIA or IVA group in the periodic table, where M is a transition metal and X is either carbon or nitrogen. Replacing the A-element with surface termination species Tx will separate the Mn+1Xn-layers forming two-dimensional (2D) flakes of Mn+1XnTx. For Ti3C2Tx the Tx corresponds to fluorine (F) and oxygen (O) covering both sides of every single 2D Mn+1Xn-flake. The Ti K-edge (1s) XANES of both Ti3AlC2 and Ti3C2Tx exhibit characteristic pre-edge absorption regions of C 2p - Ti 3d hybridization with clear crystal-field splitting's while the main-edge absorption features originate from the Ti 1s -> 4p excitation, where only the latter shows sensitivity towards the fcc-site occupation of the termination species. The coordination numbers obtained from EXAFS show that Ti3AlC2 and Ti3C2Tx are highly anisotropic with a strong in-plane contribution for Ti and with a dynamic out-of-plane contribution from the Al monolayers and termination species, respectively. As shown in the temperature-dependent measurements, the O contribution shifts to shorter bond length while the F diminishes as the temperature is raised from room temperature up to 750 {\deg}C.


title: Phonemer at WNUT-2020 Task 2: Sequence Classification Using COVID Twitter BERT and Bagging Ensemble Technique based on Plurality Voting

abstract: This paper presents the approach that we employed to tackle the EMNLP WNUT-2020 Shared Task 2 : Identification of informative COVID-19 English Tweets. The task is to develop a system that automatically identifies whether an English Tweet related to the novel coronavirus (COVID-19) is informative or not. We solve the task in three stages. The first stage involves pre-processing the dataset by filtering only relevant information. This is followed by experimenting with multiple deep learning models like CNNs, RNNs and Transformer based models. In the last stage, we propose an ensemble of the best model trained on different subsets of the provided dataset. Our final approach achieved an F1-score of 0.9037 and we were ranked sixth overall with F1-score as the evaluation criteria.


title: Improving spatial domain based image formation through compressed sensing

abstract: In this paper, we improve image reconstruction in a single-pixel scanning system by selecting an detector optimal field of view. Image reconstruction is based on compressed sensing and image quality is compared to interpolated staring arrays. The image quality comparisons use a "dead leaves" data set, Bayesian estimation and the Peak-Signal-to-Noise Ratio (PSNR) measure. Compressed sensing is explored as an interpolation algorithm and shows with high probability an improved performance compared to Lanczos interpolation. Furthermore, multi-level sampling in a single-pixel scanning system is simulated by dynamically altering the detector field of view. It was shown that multi-level sampling improves the distribution of the Peak-Signal-to-Noise Ratio. We further explore the expected sampling level distributions and PSNR distributions for multi-level sampling. The PSNR distribution indicates that there is a small set of levels which will improve image quality over interpolated staring arrays. We further conclude that multi-level sampling will outperform single-level uniform random sampling on average.


title: Erratum to "Frequency Linear-time Temporal Logic"

abstract: We correct our proof of a theorem stating that satisfiability of frequency linear-time temporal logic is undecidable [TASE 2012].


title: Universal time-series forecasting with mixture predictors

abstract: This book is devoted to the problem of sequential probability forecasting, that is, predicting the probabilities of the next outcome of a growing sequence of observations given the past. This problem is considered in a very general setting that unifies commonly used probabilistic and non-probabilistic settings, trying to make as few as possible assumptions on the mechanism generating the observations. A common form that arises in various formulations of this problem is that of mixture predictors, which are formed as a combination of a finite or infinite set of other predictors attempting to combine their predictive powers. The main subject of this book are such mixture predictors, and the main results demonstrate the universality of this method in a very general probabilistic setting, but also show some of its limitations. While the problems considered are motivated by practical applications, involving, for example, financial, biological or behavioural data, this motivation is left implicit and all the results exposed are theoretical. The book targets graduate students and researchers interested in the problem of sequential prediction, and, more generally, in theoretical analysis of problems in machine learning and non-parametric statistics, as well as mathematical and philosophical foundations of these fields. The material in this volume is presented in a way that presumes familiarity with basic concepts of probability and statistics, up to and including probability distributions over spaces of infinite sequences. Familiarity with the literature on learning or stochastic processes is not required.

categories: cs.LG cs.AI cs.IT math.IT math.ST stat.ML stat.TH


title: Onset of synchronization in coupled Mixmaster oscillators

abstract: We consider the problem of asymptotic synchronization of different spatial points coupled to each other in inhomogeneous spacetime and undergoing chaotic Mixmaster oscillations towards the singularity. We demonstrate that for couplings larger than some threshold value, two Mixmaster spatial points $A,B$, with $A$ in the past of $B$, synchronize and thereby proceed in perfect unison towards the initial singularity. We further show that there is a Lyapunov function for the synchronization dynamics that makes different spatial points able to synchronize exponentially fast in the past direction. We provide an elementary proof of how an arbitrary spatial point responds to the mean field created by the oscillators, leading to their direct interaction through spontaneous synchronization. These results ascribe a clear physical meaning of early-time synchronization leading to a resetting effect for the two BKL maps corresponding to two distinct oscillating spatial points, as the two maps converge to each other to become indistinguishable at the end of synchronization. Our results imply that the universe generically organizes itself through simpler, synchronized, states as it approaches the initial singularity. A discussion of further implications of early-time inhomogeneous Mixmaster synchronization is also provided.


title: Frequency-dependent magnetic susceptibility of magnetic nanoparticles in a polymer solution: a simulation study

abstract: Magnetic composite materials i.e. elastomers, polymer gels, or polymer solutions with embedded magnetic nanoparticles are useful for many technical and bio-medical applications. However, the microscopic details of the coupling mechanisms between the magnetic properties of the particles and the mechanical properties of the (visco)elastic polymer matrix remain unresolved. Here we study the response of a single-domain spherical magnetic nanoparticle that is suspended in a polymer solution to alternating magnetic fields. As interactions we consider only excluded volume interactions with the polymers and hydrodynamic interactions mediated through the solvent. The AC susceptibility spectra are calculated using a linear response Green-Kubo approach, and the influences of changing polymer concentration and polymer length are investigated. Our data is compared to recent measurements of the AC susceptibility for a typical magnetic composite system [Roeben et al., Colloid and Polymer Science, 2014, 2013--2023], and demonstrates the importance of hydrodynamic coupling in such systems.


title: OutbreakFlow: Model-based Bayesian inference of disease outbreak dynamics with invertible neural networks and its application to the COVID-19 pandemics in Germany

abstract: Mathematical models in epidemiology are an indispensable tool to determine the dynamics and important characteristics of infectious diseases. Apart from their scientific merit, these models are often used to inform political decisions and intervention measures during an ongoing outbreak. However, reliably inferring the dynamics of ongoing outbreaks by connecting complex models to real data is still hard and requires either laborious manual parameter fitting or expensive optimization methods which have to be repeated from scratch for every application of a given model. In this work, we address this problem with a novel combination of epidemiological modeling with specialized neural networks. Our approach entails two computational phases: In an initial training phase, a mathematical model describing the epidemic is used as a coach for a neural network, which acquires global knowledge about the full range of possible disease dynamics. In the subsequent inference phase, the trained neural network processes the observed data of an actual outbreak and infers the parameters of the model in order to realistically reproduce the observed dynamics and reliably predict future progression. With its flexible framework, our simulation-based approach is applicable to a variety of epidemiological models. Moreover, since our method is fully Bayesian, it is designed to incorporate all available prior knowledge about plausible parameter values and returns complete joint posterior distributions over these parameters. Application of our method to the early Covid-19 outbreak phase in Germany demonstrates that we are able to obtain reliable probabilistic estimates for important disease characteristics, such as generation time, fraction of undetected infections, likelihood of transmission before symptom onset, and reporting delays using a very moderate amount of real-world observations.