title: GraphBGS: Background Subtraction via Recovery of Graph Signals

abstract: Background subtraction is a fundamental pre-processing task in computer vision. This task becomes challenging in real scenarios due to variations in the background for both static and moving camera sequences. Several deep learning methods for background subtraction have been proposed in the literature with competitive performances. However, these models show performance degradation when tested on unseen videos; and they require huge amount of data to avoid overfitting. Recently, graph-based algorithms have been successful approaching unsupervised and semi-supervised learning problems. Furthermore, the theory of graph signal processing and semi-supervised learning have been combined leading to new insights in the field of machine learning. In this paper, concepts of recovery of graph signals are introduced in the problem of background subtraction. We propose a new algorithm called Graph BackGround Subtraction (GraphBGS), which is composed of: instance segmentation, background initialization, graph construction, graph sampling, and a semi-supervised algorithm inspired from the theory of recovery of graph signals. Our algorithm has the advantage of requiring less labeled data than deep learning methods while having competitive results on both: static and moving camera videos. GraphBGS outperforms unsupervised and supervised methods in several challenging conditions on the publicly available Change Detection (CDNet2014), and UCSD background subtraction databases.


title: RR charge calculation for D brane in B field

abstract: The paper is connected with searches for the Ramond-Ramond charge of D branes in the presence of B field. The consideration of B field inclusion is an important physical and mathematical unsolved problem, which is connected with K group calculations of twisted bundles. Considered two cases of vector bundles, Azumaya and Rosenberg algebras and analyzed their K group realization.


title: Dynamics of the mean-field interacting quantum kicked rotor

abstract: We study the dynamics of the many-body atomic kicked rotor with interactions at the mean-field level, governed by the Gross-Pitaevskii equation. We show that dynamical localization is destroyed by the interaction, and replaced by a subdiffusive behavior. In contrast to results previously obtained from a simplified version of the Gross-Pitaevskii equation, the subdiffusive exponent does not appear to be universal. By studying the phase of the mean-field wave function, we propose a new approximation that describes correctly the dynamics at experimentally relevant times close to the start of subdiffusion, while preserving the reduced computational cost of the former approximation.


title: Counting difficult tree pairs with respect to the rotation distance problem

abstract: Rotation distance between rooted binary trees is the minimum number of simple rotations needed to transform one tree into the other. Computing the rotation distance between a pair of rooted trees can be quickly reduced in cases where there is a common edge between the trees, or where a single rotation introduces a common edge. Tree pairs which do not have such a reduction are difficult tree pairs, where there is no generally known first step. Here, we describe efforts to count and estimate the number of such difficult tree pairs, and find that the fraction decreases exponentially fast toward zero. We also describe how knowing the number of distinct instances of the rotation distance problem is a helpful factor in making the computations more feasible.


title: A Technical Critique of Some Parts of the Free Energy Principle

abstract: We summarize the original formulation of the free energy principle, and highlight some technical issues. We discuss how these issues affect related results involving generalised coordinates and, where appropriate, mention consequences for and reveal, up to now unacknowledged, differences to newer formulations of the free energy principle. In particular, we reveal that various definitions of the "Markov blanket" proposed in different works are not equivalent. We show that crucial steps in the free energy argument which involve rewriting the equations of motion of systems with Markov blankets, are not generally correct without additional (previously unstated) assumptions. We prove by counterexample that the original free energy lemma, when taken at face value, is wrong. We show further that this free energy lemma, when it does hold, implies equality of variational density and ergodic conditional density. The interpretation in terms of Bayesian inference hinges on this point, and we hence conclude that it is not sufficiently justified. Additionally, we highlight that the variational densities presented in newer formulations of the free energy principle and lemma are parameterised by different variables than in older works, leading to a substantially different interpretation of the theory. Note that we only highlight some specific problems in the discussed publications. These problems do not rule out conclusively that the general ideas behind the free energy principle are worth pursuing.


title: Subjective Annotation for a Frame Interpolation Benchmark using Artefact Amplification

abstract: Current benchmarks for optical flow algorithms evaluate the estimation either directly by comparing the predicted flow fields with the ground truth or indirectly by using the predicted flow fields for frame interpolation and then comparing the interpolated frames with the actual frames. In the latter case, objective quality measures such as the mean squared error are typically employed. However, it is well known that for image quality assessment, the actual quality experienced by the user cannot be fully deduced from such simple measures. Hence, we conducted a subjective quality assessment crowdscouring study for the interpolated frames provided by one of the optical flow benchmarks, the Middlebury benchmark. We collected forced-choice paired comparisons between interpolated images and corresponding ground truth. To increase the sensitivity of observers when judging minute difference in paired comparisons we introduced a new method to the field of full-reference quality assessment, called artefact amplification. From the crowdsourcing data, we reconstructed absolute quality scale values according to Thurstone's model. As a result, we obtained a re-ranking of the 155 participating algorithms w.r.t. the visual quality of the interpolated frames. This re-ranking not only shows the necessity of visual quality assessment as another evaluation metric for optical flow and frame interpolation benchmarks, the results also provide the ground truth for designing novel image quality assessment (IQA) methods dedicated to perceptual quality of interpolated images. As a first step, we proposed such a new full-reference method, called WAE-IQA. By weighing the local differences between an interpolated image and its ground truth WAE-IQA performed slightly better than the currently best FR-IQA approach from the literature.


title: A symmetric formula of transformed elasticity tensor in PML domain for elastic wave problem

abstract: The perfectly matched layer is very important for the elastic wave problem in the frequency domain. Generally, the formulas of the elasticity tensor in the perfectly matched layers are derived from the transformed momentum equation. In this note, we proved that the transformed elasticity tensor derived in this way lost its symmetry. Therefore, these formulas are inconsistency in theory and it's hard to explain its numerical performance. We present a new symmetrical formula of elasticity tensor from the weak form. So the theory of elasticity is still applicable in the perfectly matched layers.


title: On the asymmetry of stars at infinity

abstract: Given a bordified space, Karlsson defines an incidence geometry of stars at infinity. These stars and their incidence are closely related to well-understood objects when the space is hyperbolic, CAT(0), or a bounded convex domain with the Hilbert metric. A question stemming from Karlsson's original paper was whether or not the relation of one boundary point being included in a star of another boundary point is symmetric. This paper provides an example demonstrating that this relation in the star boundary of the three-tree Diestel-Leader graph $DL_3(q)$ is not symmetric. In doing so, some interesting bounds on distance in Diestel-Leader graphs are utilized.


title: The sub-fractional CEV model

abstract: The sub-fractional Brownian motion (sfBm) is a stochastic process, characterized by non-stationarity in their increments and long-range dependency, considered as an intermediate step between the standard Brownian motion (Bm) and the fractional Brownian motion (fBm). The mixed process, a linear combination between a Bm and an independent sfBm, called mixed sub-fractional Brownian motion (msfBm), keeps the features of the sfBm adding the semi-martingale property for H>3/4, is a suitable candidate to use in price fluctuation modeling, in particular for option pricing. In this note, we arrive at the European Call price under the Constant Elasticity of Variance (CEV) model driven by a mixed sub-fractional Brownian motion. Empirical tests show the capacity of the proposed model to capture the temporal structure of option prices across different maturities.


title: From constant to variable density inverse extended Born modelling

abstract: For quantitative seismic imaging, iterative least-squares reverse time migration is the recommended approach. The existence of an inverse of the forward modelling operator would considerably reduce the number of required iterations. In the context of the extended model, such a pseudo-inverse exists, built as a weighted version of the adjoint and accounts for the deconvolution, geometrical spreading and uneven illumination. The application of the pseudo-inverse Born modelling is based on constant density acoustic media, which is a limiting factor for practical applications. To consider density perturbation, we propose and investigate two approaches. The first one is a generalization of a recent study proposing to recover acoustic perturbations from angle-dependent response of the pseudo-inverse Born modelling operator. The new version is based on weighted least-squares objective function. The method not only provides more robust results, but also offers the flexibility to include constrains in the objective function in order to reduce the parameters cross-talk. We also propose an alternative approach based on Taylor expansion that does not require any Radon transform. Numerical examples based on simple and the Marmousi2 models using correct and incorrect background models for the variable density Born modelling, verify the effectiveness of the weighted least-squares method when compared with the other two approaches. The Taylor expansion approach appears to contain too many artifacts for a successful applicability.


title: Practical data-driven flood forecasting based on dynamical systems theory: Case studies from Japan

abstract: Data-driven flood forecasting methods are useful, especially for the rivers that lack hydrological information to build physical models. Although these former methods can forecast river stages using only past water levels and rainfall data, they cannot handle previously unexperienced water levels easily, and require a large amount of data to build accurate models. Here, we focus on phase-space reconstruction approaches, and develop a practical data-driven forecasting method to overcome the existing problems. The proposed method can handle the unexperienced water levels and provide forecasts using only a small number of water rise events. We apply the method to data from actual rivers, and it achieved the best forecast performance among existing methods, including a physical runoff model, a data-driven multi-layer perceptron, and a conventional method based on phase-space reconstruction. In addition, the proposed method also forecasted the exceedance of the evacuation warning level 6 h earlier for steep rivers. Given its performance and maintainability, the proposed method can be applied to many actual rivers for early evacuation.


title: Free-ranging dogs do not distinguish between barks without context

abstract: Canids display a vast diversity of social organizations, from solitary-living to pairs to packs. Domestic dogs have descended from pack-living gray wolf-like ancestors. Unlike their group living ancestors, free-ranging dogs are facultatively social, preferring to forage solitarily. They are scavengers by nature, mostly dependent on human garbage and generosity for their sustenance. Free-ranging dogs are highly territorial, often defending their territories using vocalizations. Vocal communication plays a critical role between inter and intraspecies and group interaction and maintaining their social dynamics. Barking is the most common among the different types of vocalizations of dogs. Dogs have a broad hearing range and can respond to sounds over long distances. Domestic dogs have been shown to have the ability to distinguish between barking in different contexts. Since free-ranging dogs regularly engage in various kinds of interactions with each other, it is interesting to know whether they are capable of distinguishing between vocalizations of their own and other groups. In this study, a playback experiment was used to test if dogs can distinguish between barking of their own group member from a non-group member. Though dogs respond to barking from other groups in territorial exchanges, they did not respond differently to the self and other group barking in the playback experiments. This suggests a role of context in the interactions between dogs and opens up possibilities for future studies on the comparison of the responses of dogs in playback experiments with their natural behavior through long-term observations.


title: SeismiQB -- a novel framework for deep learning with seismic data

abstract: In recent years, Deep Neural Networks were successfully adopted in numerous domains to solve various image-related tasks, ranging from simple classification to fine borders annotation. Naturally, many researches proposed to use it to solve geological problems. Unfortunately, many of the seismic processing tools were developed years before the era of machine learning, including the most popular SEG-Y data format for storing seismic cubes. Its slow loading speed heavily hampers experimentation speed, which is essential for getting acceptable results. Worse yet, there is no widely-used format for storing surfaces inside the volume (for example, seismic horizons). To address these problems, we've developed an open-sourced Python framework with emphasis on working with neural networks, that provides convenient tools for (i) fast loading seismic cubes in multiple data formats and converting between them, (ii) generating crops of desired shape and augmenting them with various transformations, and (iii) pairing cube data with labeled horizons or other types of geobodies.


title: Experimental passive state preparation for continuous variable quantum communications

abstract: In the Gaussian-modulated coherent state quantum key distribution (QKD) protocol, the sender first generates Gaussian distributed random numbers and then encodes them on weak laser pulses actively by performing amplitude and phase modulations. Recently, an equivalent passive QKD scheme was proposed by exploring the intrinsic field fluctuations of a thermal source [B. Qi, P. G. Evans, and W. P. Grice, Phys. Rev. A 97, 012317 (2018)]. This passive QKD scheme is especially appealing for chip-scale implementation since no active modulations are required. In this paper, we conduct an experimental study of the passively encoded QKD scheme using an off-the-shelf amplified spontaneous emission source operated in continuous-wave mode. Our results show that the excess noise introduced by the passive state preparation scheme can be effectively suppressed by applying optical attenuation and secure key could be generated over metro-area distances.


title: Chiral effective Lagrangian for heavy-light mesons from QCD

abstract: We derive the chiral effective Lagrangian for heavy-light mesons in the heavy quark limit from QCD under proper approximations. The low energy constants in the effective Lagrangian are expressed in terms of the light quark self-energy. With typical forms of the running coupling constant of QCD and the quark self-energy obtained from Dyson-Schwinger equations as well as lattice QCD, we estimate the low energy constants in the model and the strong decay widths. A comparison with data and some discussions of the numerical results are presented.


title: Beam energy dependence of net-$\Lambda$ fluctuations measured by the STAR experiment at RHIC

abstract: The measurements of particle multiplicity distributions have generated considerable interest in understanding the fluctuations of conserved quantum numbers in the Quantum Chromodynamics (QCD) hadronization regime, in particular near a possible critical point and near the chemical freeze-out. We report the measurement of efficiency and centrality bin width corrected cumulant ratios ($C_{2}/C_{1}$, $C_{3}/C_{2}$) of net-$\Lambda$ distributions, in the context of both strangeness and baryon number conservation, as a function of collision energy, centrality and rapidity. The results are for Au + Au collisions at five beam energies ($\sqrt{s_{NN}}$ = 19.6, 27, 39, 62.4 and 200 GeV) recorded with the Solenoidal Tracker at RHIC (STAR). We compare our results to the Poisson and negative binomial (NBD) expectations, as well as to Ultra-relativistic Quantum Molecular Dynamics (UrQMD) and Hadron Resonance Gas (HRG) model predictions. Both NBD and Poisson baselines agree with data within the statistical and systematic uncertainties. The ratios of the measured cumulants show no features of critical fluctuations. The chemical freeze-out temperatures extracted from a recent HRG calculation, which was successfully used to describe the net-proton, net-kaon and net-charge data, indicate $\Lambda$ freeze-out conditions similar to those of kaons. However, large deviations are found when comparing to temperatures obtained from net-proton fluctuations. The net-$\Lambda$ cumulants show a weak, but finite, dependence on the rapidity coverage in the acceptance of the detector, which can be attributed to quantum number conservation.


title: IPPO: A Privacy-Aware Architecture for Decentralized Data-sharing

abstract: Online trackers personalize ads campaigns, exponentially increasing their efficacy compared to traditional channels. The downside of this is that thousands of mostly unknown systems own our profiles and violate our privacy without our awareness. IPPO turns the table and re-empower users of their data, through anonymised data publishing via a Blockchain-based Decentralized Data Marketplace. We also propose a service based on machine learning and big data analytics to automatically identify web trackers and build Privacy Labels (PLs), based on the nutrition labels concept. This paper describes the motivation, the vision, the architecture and the research challenges related to IPPO.


title: Electronic bridge excitation in highly charged Th-229 ions

abstract: The excitation of the 8 eV $^{229m}$Th isomer through the electronic bridge mechanism in highly charged ions is investigated theoretically. By exploiting the rich level scheme of open $4f$ orbitals and the robustness of highly charged ions against photoionization, a pulsed high-intensity optical laser can be used to efficiently drive the nuclear transition by coupling it to the electronic shell. We show how to implement a promising electronic bridge scheme in an electron beam ion trap starting from a metastable electronic state. This setup would avoid the need for a tunable vacuum ultraviolet laser. Based on our theoretical predictions, determining the isomer energy with an uncertainty of $10^{-5}$ eV could be achieved in one day of measurement time using realistic laser parameters.


title: An efficient sampling algorithm for difficult tree pairs

abstract: It is an open question whether there exists a polynomial-time algorithm for computing the rotation distances between pairs of extended ordered binary trees. The problem of computing the rotation distance between an arbitrary pair of trees, (S, T), can be efficiently reduced to the problem of computing the rotation distance between a difficult pair of trees (S', T'), where there is no known first step which is guaranteed to be the beginning of a minimal length path. Of interest, therefore, is how to sample such difficult pairs of trees of a fixed size. We show that it is possible to do so efficiently, and present such an algorithm that runs in time $O(n^4)$.


title: InChorus: Designing Consistent Multimodal Interactions for Data Visualization on Tablet Devices

abstract: While tablet devices are a promising platform for data visualization, supporting consistent interactions across different types of visualizations on tablets remains an open challenge. In this paper, we present multimodal interactions that function consistently across different visualizations, supporting common operations during visual data analysis. By considering standard interface elements (e.g., axes, marks) and grounding our design in a set of core concepts including operations, parameters, targets, and instruments, we systematically develop interactions applicable to different visualization types. To exemplify how the proposed interactions collectively facilitate data exploration, we employ them in a tablet-based system, InChorus that supports pen, touch, and speech input. Based on a study with 12 participants performing replication and fact-checking tasks with InChorus, we discuss how participants adapted to using multimodal input and highlight considerations for future multimodal visualization systems.


title: Solution of the unconditional extremum problem for a liner-fractional integral functional on a set of probability measures and its application in the theory of optimal control of semi-Markov processes

abstract: In this paper, a new method for solving the problem of optimal control of semi-Markov processes with finitely many states is considered. A new form of the assertion on an extremum of a liner-fractional integral functional given on a set of probability measures is formulated and proved. This form underlies the theorem of optimal control strategy for semi-Markov processes. It is proved that the solution of the optimal control problem for a semi-Markov process with finitely many states is completely determined by the extremum properties of the so-called test function of the liner-fractional integral functional which is the control quality index. At the same time, an explicit analytic representation was obtained for this test function in terms of the initial probability characteristics of the semi-Markov model.


title: Derivation of a refined 6-parameter shell model: Descent from the three-dimensional Cosserat elasticity using a method of classical shell theory

abstract: Starting from the three-dimensional Cosserat elasticity, we derive a two-dimensional model for isotropic elastic shells. For the dimensional reduction, we employ a derivation method similar to that used in classical shell theory, as presented systematically by Steigmann in [J. Elast. \textbf{111}: 91-107, 2013]. As a result, we obtain a geometrically nonlinear Cosserat shell model with a specific form of the strain-energy density, which has a simple expression with coefficients depending on the initial curvature tensor and on three-dimensional material constants. The explicit forms of the stress-strain relations and the local equilibrium equations are also recorded. Finally, we compare our results with other 6-parameter shell models and discuss the relation to the classical Koiter shell model.


title: Eye of the Mind: Image Processing for Social Coding

abstract: Developers are increasingly sharing images in social coding environments alongside the growth in visual interactions within social networks. The analysis of the ratio between the textual and visual content of Mozilla's change requests and in Q/As of StackOverflow programming revealed a steady increase in sharing images over the past five years. Developers' shared images are meaningful and are providing complementary information compared to their associated text. Often, the shared images are essential in understanding the change requests, questions, or the responses submitted. Relying on these observations, we delve into the potential of automatic completion of textual software artifacts with visual content.


title: TailorGAN: Making User-Defined Fashion Designs

abstract: Attribute editing has become an important and emerging topic of computer vision. In this paper, we consider a task: given a reference garment image A and another image B with target attribute (collar/sleeve), generate a photo-realistic image which combines the texture from reference A and the new attribute from reference B. The highly convoluted attributes and the lack of paired data are the main challenges to the task. To overcome those limitations, we propose a novel self-supervised model to synthesize garment images with disentangled attributes (e.g., collar and sleeves) without paired data. Our method consists of a reconstruction learning step and an adversarial learning step. The model learns texture and location information through reconstruction learning. And, the model's capability is generalized to achieve single-attribute manipulation by adversarial learning. Meanwhile, we compose a new dataset, named GarmentSet, with annotation of landmarks of collars and sleeves on clean garment images. Extensive experiments on this dataset and real-world samples demonstrate that our method can synthesize much better results than the state-of-the-art methods in both quantitative and qualitative comparisons.


title: Analytic Classification of Germs of Parabolic Antoholomorphic Diffeomorphisms of Codimension k

abstract: We investigate the local dynamics of antiholomorphic diffeomorphisms around a parabolic fixed point. We first give a normal form. Then we give a complete classification including a modulus space for antiholomorphic germs with a parabolic fixed point under analytic conjugacy. We then study some geometric applications: existence of real analytic invariant curve, existence of holomorphic and antiholomorphic roots of holomorphic and antiholomorphic parabolic germs, commuting holomorphic and antiholomorphic parabolic germs.


title: A Brillouin Laser Optical Atomic Clock

abstract: Over the last decade, optical atomic clocks have surpassed their microwave counterparts and now offer the ability to measure time with an increase in precision of two orders of magnitude or more. This performance increase is compelling not only for enabling new science, such as geodetic measurements of the earth, searches for dark matter, and investigations into possible long-term variations of fundamental physics constants but also for revolutionizing existing technology, such as the global positioning system (GPS). A significant remaining challenge is to transition these optical clocks to non-laboratory environments, which requires the ruggedization and miniaturization of the atomic reference and clock laser along with their supporting lasers and electronics. Here, using a compact stimulated Brillouin scattering (SBS) laser to interrogate a $^8$$^8$Sr$^+$ ion, we demonstrate a promising component of a portable optical atomic clock architecture. In order to bring the stability of the SBS laser to a level suitable for clock operation, we utilize a self-referencing technique to compensate for temperature drift of the laser to within $170$ nK. Our SBS optical clock achieves a short-term stability of $3.9 \times 10^{-14}$ at $1$ s---an order of magnitude improvement over state-of-the-art microwave clocks. Based on this technology, a future GPS employing portable SBS clocks offers the potential for distance measurements with a 100-fold increase in resolution.


title: Detection of Fe I in the atmosphere of the ultra-hot Jupiter WASP-121b, and a new likelihood-based approach for Doppler-resolved spectroscopy

abstract: High-resolution Doppler-resolved spectroscopy has opened up a new window into the atmospheres of both transiting and non-transiting exoplanets. Here, we present VLT/UVES observations of a transit of WASP-121b, an 'ultra-hot' Jupiter previously found to exhibit a temperature inversion and detections of multiple species at optical wavelengths. We present initial results using the blue arm of UVES ($\approx$3700-5000A), recovering a clear signal of neutral Fe in the planet's atmosphere at >8$\sigma$, which could contribute to (or even fully explain) the temperature inversion in the stratosphere. However, using standard cross-correlation methods, it is difficult to extract physical parameters such as temperature and abundances. Recent pioneering efforts have sought to develop likelihood `mappings' that can be used to directly fit models to high-resolution datasets. We introduce a new framework that directly computes the likelihood of the model fit to the data, and can be used to explore the posterior distribution of parameterised model atmospheres via MCMC techniques. Our method also recovers the physical extent of the atmosphere, as well as account for time- and wavelength-dependent uncertainties. We measure a temperature of $3710^{+490}_{-510}$K, indicating a higher temperature in the upper atmosphere when compared to low-resolution observations. We also show that the Fe I signal is physically separated from the exospheric Fe II. However, the temperature measurements are highly degenerate with aerosol properties; detection of additional species, using more sophisticated atmospheric models, or combining these methods with low-resolution spectra should help break these degeneracies.


title: Are the Carriers of Diffuse Interstellar Bands and Extended Red Emission the same?

abstract: We report the first spectroscopic observations of a background star seen through the region between the ionization front and the dissociation front of the nebula IC 63. This photodissociation region (PDR) exhibits intense extended red emission (ERE) attributed to fluorescence by large molecules/ions. We detected strong diffuse interstellar bands (DIB) in the stellar spectrum, including an exceptionally strong and broad DIB at $\lambda$4428. The detection of strong DIBs in association with ERE could be consistent with the suggestion that the carriers of DIBs and ERE are identical. The likely ERE process is recurrent fluorescence, enabled by inverse internal conversions from highly excited vibrational levels of the ground state to low-lying electronic states with subsequent transitions to ground. This provides a path to rapid radiative cooling for molecules/molecular ions, greatly enhancing their ability to survive in a strongly irradiated environment. The ratio of the equivalent widths (EW) of DIBs $\lambda$5797 and $\lambda$5780 in IC 63 is the same as that observed in the low-density interstellar medium with UV interstellar radiation fields (ISRF) weaker by at least two orders of magnitude. This falsifies suggestions that the ratio of these two DIBs can serve as a measure of the UV strength of the ISRF. Observations of the nebular spectrum of the PDR of IC 63 at locations immediately adjacent to where DIBs were detected failed to reveal any presence of sharp emission features seen in the spectrum of the Red Rectangle nebula. This casts doubts upon proposals that the carriers of these features are the same as those of DIBs seen at slightly shorter wavelengths.


title: Universality classes of spin transport in one-dimensional isotropic magnets: the onset of logarithmic anomalies

abstract: We report a systematic study of finite-temperature spin transport in quantum and classical one-dimensional magnets with isotropic spin interactions, including both integrable and non-integrable models. Employing a phenomenological framework based on a generalized Burgers' equation in a time-dependent stochastic environment, we identify four different universality classes of spin fluctuations. These comprise, aside from normal spin diffusion, three types of superdiffusive transport: the KPZ universality class and two distinct types of anomalous diffusion with multiplicative logarithmic corrections. Our predictions are supported by extensive numerical simulations on various examples of quantum and classical chains. Contrary to common belief, we demonstrate that even non-integrable spin chains can display a diverging spin diffusion constant at finite temperatures.


title: Domain wall induced spin-polarized flat bands in antiferromagnetic topological insulators

abstract: A flat band in fermionic system is a dispersionless single-particle state with a diverging effective mass and nearly zero group velocity. These flat bands are expected to support exotic properties in the ground state, which might be important for a wide range of promising physical phenomena. For many applications it is highly desirable to have such states in Dirac materials, but so far they have been reported only in non-magnetic Dirac systems. In this work we propose a realization of topologically protected spin-polarized flat bands generated by domain walls in planar magnetic topological insulators. Using first-principles material design we suggest a family of intrinsic antiferromagnetic topological insulators with an in-plane sublattice magnetization and a high N\'eel temperature. Such systems can host domain walls in a natural manner. For these materials, we demonstrate the existence of spin-polarized flat bands in the vicinity of the Fermi level and discuss their properties and potential applications.


title: Sinogram super-resolution and denoising convolutional neural network (SRCN) for limited data photoacoustic tomography

abstract: The quality of the reconstructed photoacoustic image largely depends on the amount of photoacoustic (PA) boundary data available, which in turn is proportional to the number of detectors employed. In case of limited data (owing to less number of detectors due to cost/instrumentation constraints), the reconstructed PA images suffer from artifacts and are often noisy. In this work, for the first time, a deep learning based model was developed to super resolve and denoise the photoacoustic sinogram data. The proposed method was compared with existing nearest neighbor interpolation and wavelet based denoising techniques and was shown to outperform them both in numerical and in-vivo cases. The improvement obtained in Root Mean Square Error (RMSE) and Peak Signal to Noise Ratio (PSNR) for the reconstructed PA image using the sinogram data that was super-resolved and denoised using proposed neural network based method was as high as 41.70 % and 6.93 dB respectively compared to utilizing limited sinogram data.


title: Addendum and Erratum: Mapping cones for morphisms involving a band complex in the bounded derived category of a gentle algebra

abstract: In this note we correct two oversights in [Mapping cones in the bounded derived category of a gentle algebra, J. Algebra 530 (2019), 163--194, also arXiv:1609.09688] which only occur when a band complex is involved. As a consequence we see that the mapping cone of a morphism between two band complexes can decompose into arbitrarily many indecomposable direct summands.


title: On Visibility Graphs of Convex Fans and Terrains

abstract: For two points in the closure of a simple polygon $P$, we say that they see each other in $P$ if the line segment uniting them does not intersect the exterior of $P$. The visibility graph of $P$ is the graph whose vertex set is the vertex set of $P$ and two vertices are joined by an edge if they see each other in $P$. The characterization of visibility graphs has been an open problem for decades, and a significant effort has been made to characterize visibility graphs of restricted polygon classes. Among them is the convex fan: a simple polygon with a convex vertex (i.e. whose internal angle in less than 180 degrees) that sees every other point in the closure of the polygon (often called kernel vertex). We show that visibility graphs of convex fans are equivalent to visibility graphs of terrains with the addition of a universal vertex, that is, a vertex that is adjacent to every other vertex.


title: A Novel Methodology for designing Policies in Mobile Crowdsensing Systems

abstract: Mobile crowdsensing is a people-centric sensing system based on users' contributions and incentive mechanisms aim at stimulating them. In our work, we have rethought the design of incentive mechanisms through a game-theoretic methodology. Thus, we have introduced a multi-layer social sensing framework, where humans as social sensors interact on multiple social layers and various services. We have proposed to weigh these dynamic interactions by including the concept of homophily and we have modelled the evolutionary dynamics of sensing behaviours by defining a mathematical framework based on multiplex EGT, quantifying the impact of homophily, network heterogeneity and various social dilemmas. We have detected the configurations of social dilemmas and network structures that lead to the emergence and sustainability of human cooperation. Moreover, we have defined and evaluated local and global Nash equilibrium points by including the concepts of homophily and heterogeneity. We have analytically defined and measured novel statistical measures of social honesty, QoI and users' behavioural reputation scores based on the evolutionary dynamics. We have defined the Decision Support System and a novel incentive mechanism by operating on the policies in terms of users' reputation scores, that also incorporate users' behaviours other than quality and quantity of contributions. Experimentally, we have considered the Waze dataset on vehicular traffic monitoring application and derived the disbursement of incentives comparing our method with baselines. Results demonstrate that our methodology, which also includes the local (microscopic) spatio-temporal distribution of behaviours, is able to better discriminate users' behaviours. This multi-scale characterisation of users represents a novel research direction and paves the way for novel policies on mobile crowdsensing systems.


title: Strong solutions to the Stefan problem with Gibbs-Thomson correction and boundary contact

abstract: We prove existence and uniqueness of strong solutions to the two-phase Stefan problem with Gibbs-Thomson law where the free interface forms a ninety degree contact angle with the fixed boundary. We also discuss existence of global solutions and convergence to equilibria.


title: Log-concavity and log-convexity of moments of averages of i.i.d. random variables

abstract: We show that the sequence of moments of order less than 1 of averages of i.i.d. positive random variables is log-concave. For moments of order at least 1, we conjecture that the sequence is log-convex and show that this holds eventually for integer moments (after neglecting the first $p^2$ terms of the sequence).


title: Combining PRNU and noiseprint for robust and efficient device source identification

abstract: PRNU-based image processing is a key asset in digital multimedia forensics. It allows for reliable device identification and effective detection and localization of image forgeries, in very general conditions. However, performance impairs significantly in challenging conditions involving low quality and quantity of data. These include working on compressed and cropped images, or estimating the camera PRNU pattern based on only a few images. To boost the performance of PRNU-based analyses in such conditions we propose to leverage the image noiseprint, a recently proposed camera-model fingerprint that has proved effective for several forensic tasks. Numerical experiments on datasets widely used for source identification prove that the proposed method ensures a significant performance improvement in a wide range of challenging situations.


title: Many-body Electronic Structure of NdNiO$_2$ and CaCuO$_2$

abstract: The demonstration of superconductivity in nickelate analogues of high $T_c$ cuprates provides new perspectives on the physics of correlated electron materials. The degree to which the nickelate electronic structure is similar to that of cuprates is an important open question. This paper presents results of a comparative study of the many-body electronic structure and theoretical phase diagram of the isostructural materials CaCuO$_2$ and NdNiO$_2$. Important differences include the proximity of the oxygen $2p$ bands to the Fermi level, the bandwidth of the transition metal-derived $3d$ bands, and the presence, in NdNiO$_2$, of both Nd-derived $5d$ states crossing the Fermi level and a van Hove singularity that crosses the Fermi level as the out of plane momentum is varied. The low energy physics of NdNiO$_2$ is found to be that of a single Ni-derived correlated band, with additional accompanying weakly correlated bands of Nd-derived states that dope the Ni-derived band. The effective correlation strength of the Ni-derived $d$-band crossing the Fermi level in NdNiO$_2$ is found to be greater than that of the Cu-derived $d$-band in CaCuO$_2$, but the predicted magnetic transition temperature of NdNiO$_2$ is substantially lower than that of CaCuO$_2$ because of the smaller bandwidth.


title: Globe-hopping

abstract: We consider versions of the grasshopper problem (Goulko and Kent, 2017) on the circle and the sphere, which are relevant to Bell inequalities. For a circle of circumference $2\pi$, we show that for unconstrained lawns of any length and arbitrary jump lengths, the supremum of the probability for the grasshopper's jump to stay on the lawn is one. For antipodal lawns, which by definition contain precisely one of each pair of opposite points and have length $\pi$, we show this is true except when the jump length $\phi$ is of the form $\pi\frac{p}{q}$ with $p,q$ coprime and $p$ odd. For these jump lengths we show the optimal probability is $1 - 1/q$ and construct optimal lawns. For a pair of antipodal lawns, we show that the optimal probability of jumping from one onto the other is $1 - 1/q$ for $p,q$ coprime, $p$ odd and $q$ even, and one in all other cases. For an antipodal lawn on the sphere, it is known (Kent and Pital\'ua-Garc\'ia, 2014) that if $\phi = \pi/q$, where $q \in \mathbb N$, then the optimal retention probability of $1-1/q$ for the grasshopper's jump is provided by a hemispherical lawn. We show that in all other cases where $0<\phi < \pi/2$, hemispherical lawns are not optimal, disproving the hemispherical colouring maximality hypotheses (Kent and Pital\'ua-Garc\'ia, 2014). We discuss the implications for Bell experiments and related cryptographic tests.


title: Scaling VANET Security Through Cooperative Message Verification

abstract: VANET security introduces significant processing overhead for resource-constrained On-Board Units (OBUs). Here, we propose a novel scheme that allows secure Vehicular Communication (VC) systems to scale well beyond network densities for which existing optimization approaches could be workable, without compromising security (and privacy).


title: A workflow for seismic imaging with quantified uncertainty

abstract: The interpretation of seismic images faces challenges due to the presence of several uncertainty sources. Uncertainties exist in data measurements, source positioning, and subsurface geophysical properties. Understanding uncertainties' role and how they influence the outcome is an essential part of the decision-making process in the oil and gas industry. Geophysical imaging is time-consuming. When we add uncertainty quantification, it becomes both time and data-intensive. In this work, we propose a workflow for seismic imaging with quantified uncertainty. We build the workflow upon Bayesian tomography, reverse time migration, and image interpretation based on statistical information. The workflow explores an efficient hybrid parallel computational strategy to decrease the reverse time migration execution time. High levels of data compression are applied to reduce data transfer among workflow activities and data storage. We capture and analyze provenance data at runtime to improve workflow execution, monitoring, and debugging with negligible overhead. Numerical experiments on the Marmousi2 Velocity Model Benchmark demonstrate the workflow capabilities. We observe excellent weak and strong scalability, and results suggest that the use of lossy data compression does not hamper the seismic imaging uncertainty quantification.


title: Trading on the Floor after Sweeping the Book

abstract: Informed traders need to trade fast in order to profit from their private information before it becomes public. Fast electronic markets provide such liquidity. Slow markets provide execution in an auction based trading floor. Hybrid markets combine both execution venues. In its main result, the paper shows that to compensate for their slow and risky executions, trading floors need to be at least twice as deep as the sweeping facility. Furthermore, when a stand-alone trading floor is enhanced with the addition of a sweeping facility, overall informed trading will decline because it is easier for informed traders to extract the full value of their private info.


title: Towards geometric integration of rough differential forms

abstract: We provide a draft of a theory of geometric integration of rough differential forms which are generalizations of classical (smooth) differential forms to similar objects with very low regularity, for instance, involving H\"older continuous functions that may be nowhere differentiable. Borrowing ideas from the theory of rough paths, we show that such a geometric integration can be constructed substituting appropriately differentials with more general asymptotic expansions. This can be seen as the basis of geometric integration similar to that used in geometric measure theory, but without any underlying differentiable structure, thus allowing Lipschitz functions and rectifiable sets to be substituted by far less regular objects (e.g. H\"older functions and their images which may be purely unrectifiable). Our construction includes both the one-dimensional Young integral and multidimensional integrals introduced recently by R. Z\"ust. To simplify the exposition, we limit ourselves to integration of forms of dimensions not exceeding two.


title: Crossing estimates from metric graph and discrete GFF

abstract: We compare level-set percolation for Gaussian free fields (GFFs) defined on a rectangular subset of $\delta \mathbb{Z}^2$ to level-set percolation for GFFs defined on the corresponding metric graph as the mesh size $\delta$ goes to 0. In particular, we look at the probability that there is a path that crosses the rectangle in the horizontal direction on which the field is positive. We show this probability is strictly larger in the discrete graph. In the metric graph case, we show that for appropriate boundary conditions the probability that there exists a closed pivotal edge for the horizontal crossing event decays logarithmically in $\delta$. In the discrete graph case, we compute the limit of the probability of a horizontal crossing for appropriate boundary conditions.


title: Training Normalizing Flows with the Information Bottleneck for Competitive Generative Classification

abstract: The Information Bottleneck (IB) objective uses information theory to formulate a task-performance versus robustness trade-off. It has been successfully applied in the standard discriminative classification setting. We pose the question whether the IB can also be used to train generative likelihood models such as normalizing flows. Since normalizing flows use invertible network architectures (INNs), they are information-preserving by construction. This seems contradictory to the idea of a bottleneck. In this work, firstly, we develop the theory and methodology of IB-INNs, a class of conditional normalizing flows where INNs are trained using the IB objective: Introducing a small amount of {\em controlled} information loss allows for an asymptotically exact formulation of the IB, while keeping the INN's generative capabilities intact. Secondly, we investigate the properties of these models experimentally, specifically used as generative classifiers. This model class offers advantages such as improved uncertainty quantification and out-of-distribution detection, but traditional generative classifier solutions suffer considerably in classification accuracy. We find the trade-off parameter in the IB controls a mix of generative capabilities and accuracy close to standard classifiers. Empirically, our uncertainty estimates in this mixed regime compare favourably to conventional generative and discriminative classifiers.


title: On gauge invariance of transverse momentum dependent distributions at small x

abstract: The interplay between the small x limit of QCD amplitudes and QCD factorization at moderate x has been studied extensively in recent years. It was finally shown that semiclassical formulations of small x physics can have the form of an infinite twist framework involving Transverse Momentum Dependent (TMD) distributions in the eikonal limit. In this work, we demonstrate that small x distributions can be formulated in terms of transverse gauge links. This allows in particular for direct and efficient decompositions of observables into subamplitudes involving gauge invariant sub-operators which span parton distributions.


title: Two new outbursts and transient hard X-rays from 1E 1048.1$-$5937

abstract: Since its discovery, 1E 1048.1$-$5937 has been one of the most active magnetars, both in terms of radiative outbursts, and changes to its spin properties. Here we report on a continuing monitoring campaign with the Neil Gehrels Swift Observatory X-ray Telescope in which we observe two new outbursts from this source. The first outburst occurred in 2016 July, and the second in 2017 December, reaching peak 0.5-10 keV absorbed fluxes of $3.2^{+0.2}_{-0.3}\times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ and $2.2^{+0.2}_{-0.2}\times10^{-11}$ erg s$^{-1}$ cm$^{-2}$, respectively, factors of $\sim$5 and $\sim 4$ above the quiescent flux. Both new outbursts were accompanied by spin-up glitches with amplitudes of $\Delta\nu= 4.47(6)\times10^{-7}$ Hz and $\Delta\nu= 4.32(5)\times10^{-7}$ Hz, respectively. Following the 2016 July outburst, we observe, as for past outbursts, a period of delayed torque fluctuations, which reach a peak spin-down of $1.73\pm0.01$ times the quiescent rate, and which dominates the spin evolution compared to the spin-up glitches. We also report an observation near the peak of the first of these outbursts with NuSTAR in which hard X-ray emission is detected from the source. This emission is well characterized by an absorbed blackbody plus a broken power law, with a power-law index above $13.4\pm0.6$ keV of $0.5_{-0.2}^{+0.3}$, similar to those observed in both persistent and transient magnetars. The hard X-ray results are broadly consistent with models of electron/positron cooling in twisted magnetic field bundles in the outer magnetosphere. However the repeated outbursts and associated torque fluctuations in this source remain puzzling.


title: Coarsened mixtures of hierarchical skew normal kernels for flow cytometry analyses

abstract: Flow cytometry (FCM) is the standard multi-parameter assay for measuring single cell phenotype and functionality. It is commonly used for quantifying the relative frequencies of cell subsets in blood and disaggregated tissues. A typical analysis of FCM data involves cell classification---that is, the identification of cell subgroups in the sample---and comparisons of the cell subgroups across samples or conditions. While modern experiments often necessitate the collection and processing of samples in multiple batches, analysis of FCM data across batches is challenging because differences across samples may occur due to either true biological variation or technical reasons such as antibody lot effects or instrument optics across batches. Thus a critical step in comparative analyses of multi-sample FCM data---yet missing in existing automated methods for analyzing such data---is cross-sample calibration, whose goal is to align corresponding cell subsets across multiple samples in the presence of technical variations, so that biological variations can be meaningfully compared. We introduce a Bayesian nonparametric hierarchical modeling approach for accomplishing both calibration and cell classification simultaneously in a unified probabilistic manner. Three important features of our method make it particularly effective for analyzing multi-sample FCM data: a nonparametric mixture avoids prespecifying the number of cell clusters; a hierarchical skew normal kernel that allows flexibility in the shapes of the cell subsets and cross-sample variation in their locations; and finally the "coarsening" strategy makes inference robust to departures from the model such as heavy-tailness not captured by the skew normal kernels. We demonstrate the merits of our approach in simulated examples and carry out a case study in the analysis of two multi-sample FCM data sets.


title: Design and Analysis of Online Fountain Codes for Intermediate Performance

abstract: For the benefit of improved intermediate performance, recently online fountain codes attract much research attention. However, there is a trade-off between the intermediate performance and the full recovery overhead for online fountain codes, which prevents them to be improved simultaneously. We analyze this trade-off, and propose to improve both of these two performance. We first propose a method called Online Fountain Codes without Build-up phase (OFCNB) where the degree-1 coded symbols are transmitted at first and the build-up phase is removed to improve the intermediate performance. Then we analyze the performance of OFCNB theoretically. Motivated by the analysis results, we propose Systematic Online Fountain Codes (SOFC) to further reduce the full recovery overhead. Theoretical analysis shows that SOFC has better intermediate performance, and it also requires lower full recovery overhead when the channel erasure rate is lower than a constant. Simulation results verify the analyses and demonstrate the superior performance of OFCNB and SOFC in comparison to other online fountain codes.


title: Action Complexity for Semi-Classical Black Holes

abstract: We adapt the complexity as action prescription (CA) to a semi-classical model of two-dimensional dilaton gravity and determine the rate of increase of holographic complexity for an evaporating black hole. The results are consistent with our previous numerical results for semi-classical black hole complexity using a volume prescription (CV) in the same model, but the CA calculation is fully analytic and provides a non-trivial positive test for the holographic representation of the black hole interior.


title: A dynamic black hole corona in an active galaxy through X-ray reverberation mapping

abstract: X-ray reverberation echoes are assumed to be produced in the strongly distorted spacetime around accreting supermassive black holes. This signal allows us to spatially map the geometry of the inner accretion flow - a region which cannot yet be spatially resolved by any telescope - and provides a direct measure of the black hole mass and spin. The reverberation timescale is set by the light travel path between the direct emission from a hot X-ray corona and the reprocessed emission from the inner edge of the accretion disc. However, there is an inherent degeneracy in the reverberation signal between black hole mass, inner disc radius and height of the illuminating corona above the disc. Here, we use a long X-ray observation of the highly-variable active galaxy, IRAS 13224-3809, to track the reverberation signal as the system evolves on timescales of a day. With the inclusion of all the relativistic effects, modelling reveals that the height of the X-ray corona increases with increasing luminosity, providing a dynamic view of the inner accretion region. This simultaneous modelling allows us to break the inherent degeneracies and obtain an independent timing-based estimate for the mass and spin of the black hole. The uncertainty on black hole mass is comparable to the leading optical reverberation method, making X-ray reverberation a powerful technique, particularly for sources with low optical variability.


title: A new class of $p$-adic Lipschitz functions and multidimensional Hensel's Lemma

abstract: In this work we study $p$-adic continuous functions in several variables taking values on $\mathbb{Z}_p$. We describe the orthonormal van der Put base of these functions and study various Lipschitz conditions in several variables, generalizing previous work of Anashin. In particular, we introduce a new class of $p$-adic Lipschitz functions and study some of their properties. We also prove a Hensel's lifting lemma for this new class of functions, generalizing previous results of Yurova and Khrennikov.


title: Computational Approaches for Solving Two-Echelon Vehicle and UAV Routing Problems for Post-Disaster Humanitarian Operations

abstract: Humanitarian logistics service providers have two major responsibilities immediately after a disaster: locating trapped people and routing aid to them. These difficult operations are further hindered by failures in the transportation and telecommunications networks, which are often rendered unusable by the disaster at hand. In this work, we propose a two-echelon vehicle routing framework for performing these operations using aerial uncrewed autonomous vehicles (UAVs or drones) to address the issues associated with these failures. In our proposed framework, we assume that ground vehicles cannot reach the trapped population directly, but they can only transport drones from a depot to some intermediate locations. The drones launched from these locations serve to both identify demands for medical and other aids (e.g., epi-pens, medical supplies, dry food, water) and make deliveries to satisfy them. Specifically, we present a decision framework, in which the resulting optimization problem is formulated as a two-echelon vehicle routing problem with trucks as the first echelon vehicles and for the second echelon vehicles, we consider two types of drones. Hotspot drones have the capability of providing a cell phone and internet reception and hence are used to capture demands. Delivery drones are subsequently employed to satisfy the observed demand. To handle demand uncertainty, we decompose the decision problem into two stages: providing telecommunications capabilities in the first stage thereby capturing demand precisely, and satisfying the resulting demands in the second stage. To solve the resulting models, we propose efficient computational approaches by designing a decomposition algorithm with column generation (CG)-based heuristics to identify optimal drone routes.


title: Neglecting Uncertainties Biases House-Elevation Decisions to Manage Riverine Flood Risks

abstract: Homeowners around the world elevate houses to manage flood risks. Deciding how high to elevate a house poses a nontrivial decision problem. The U.S. Federal Emergency Management Agency (FEMA) recommends elevating existing houses to the Base Flood Elevation (the elevation of the 100-yr flood) plus a freeboard. This recommendation neglects many uncertainties. Here we analyze a case-study of riverine flood risk management using a multi-objective robust decision-making framework in the face of deep uncertainties. While the quantitative results are location-specific, the approach and overall insights are generalizable. We find strong interactions between the economic, engineering, and Earth science uncertainties, illustrating the need for expanding on previous integrated analyses to further understand the nature and strength of these connections. Considering deep uncertainties surrounding flood hazards, the discount rate, the house lifetime, and the fragility can increase the economically optimal house elevation to values well above FEMA recommendation.


title: Rational indices for quantum ground state sectors

abstract: We consider charge transport for interacting many-body systems with a gapped ground state subspace which is finitely degenerate and topologically ordered. To any locality-preserving, charge-conserving unitary that preserves the ground state space, we associate an index that is an integer multiple of $1/p$, where $p$ is the ground state degeneracy. We prove that the index is additive under composition of unitaries. This formalism gives rise to several applications: fractional quantum Hall conductance, a fractional Lieb-Schultz-Mattis theorem that generalizes the standard LSM to systems where the translation-invariance is broken, and the interacting generalization of the Avron-Dana-Zak relation between Hall conductance and the filling factor.

categories: math-ph cond-mat.stat-mech math.MP quant-ph


title: Multiple Ray Received Power Modeling for mmWave Indoor and Outdoor Scenarios

abstract: Millimeter-wave (mmWave) frequency bands are expected to be used for future 5G networks due to the availability of large unused spectrum. However, the attenuation at mmWave frequencies is high. To resolve this issue, the utilization of high gain antennas and beamforming mechanisms are widely investigated in the literature. In this work, we considered mmWave end-to-end propagation modeled by individual ray sources, and explored the effects of the number of rays in the model and radiation patterns of the deployed antennas on the received power. It is shown that taking the dominant two rays is sufficient to model the channel for outdoor open areas as opposed to the indoor corridor which needs five dominant rays to have a good fit for the measurement and simulation results. It is observed that the radiation pattern of the antenna affects the slope of the path loss. Multi-path components increase the received power, thus, for indoor corridor scenarios, path loss according to the link distance is smaller for lower gain antennas due to increased reception of reflected components. For an outdoor open area, the slope of the path loss is found to be very close to that of the free space.


title: Diagonal form of the Varchenko matrices for oriented matroids

abstract: The construction of the Varchenko matrix for hyperplane arrangements, first introduced by Alexandre Varchenko, extends naturally to oriented matroids. In this paper, we generalize the theorem of Gao and Zhang by proving that the Varchenko matrix of an oriented matroid has a diagonal form if and only if the pseudohyperplane arrangement corresponding to the oriented matroid is in semigeneral position, i.e. it does not contain a degeneracy. Furthermore, we show that the Varchenko matrix of a pseudoline arrangement has a block diagonal form. This also provides an alternative combinatorial proof for the Varchenko matrix determinant formula in dimension two.


title: Non-Stationary Saturation of Inhomogeneously Broadened EPR Lines

abstract: Non-stationary saturation of inhomogeneously broadened EPR lines is studied when cross-relaxation has the characteristics of spectral diffusion. A system of generalized kinetic equations is solved in quadratures in this approximation. The result is valid not only when the contribution of the spectral diffusion is negligible or dominant, but also in the intermediate case.


title: A Study of Mental Maps in Immersive Network Visualization

abstract: The visualization of a network influences the quality of the mental map that the viewer develops to understand the network. In this study, we investigate the effects of a 3D immersive visualization environment compared to a traditional 2D desktop environment on the comprehension of a network's structure. We compare the two visualization environments using three tasks--interpreting network structure, memorizing a set of nodes, and identifying the structural changes--commonly used for evaluating the quality of a mental map in network visualization. The results show that participants were able to interpret network structure more accurately when viewing the network in an immersive environment, particularly for larger networks. However, we found that 2D visualizations performed better than immersive visualization for tasks that required spatial memory.


title: Plato Dialogue System: A Flexible Conversational AI Research Platform

abstract: As the field of Spoken Dialogue Systems and Conversational AI grows, so does the need for tools and environments that abstract away implementation details in order to expedite the development process, lower the barrier of entry to the field, and offer a common test-bed for new ideas. In this paper, we present Plato, a flexible Conversational AI platform written in Python that supports any kind of conversational agent architecture, from standard architectures to architectures with jointly-trained components, single- or multi-party interactions, and offline or online training of any conversational agent component. Plato has been designed to be easy to understand and debug and is agnostic to the underlying learning frameworks that train each component.


title: Quantum Optimal Control via Magnus Expansion and Non-Commutative Polynomial Optimization

abstract: Quantum optimal control has numerous important applications ranging from pulse shaping in magnetic-resonance imagining to laser control of chemical reactions and quantum computing. Our objective is to address two major challenges that have limited the success of applications of quantum optimal control so far: non-commutativity inherent in quantum systems and non-convexity of quantum optimal control problems involving more than three quantum levels. Methodologically, we address the non-commutativity of the control Hamiltonian at different times by the use of Magnus expansion. To tackle the non-convexity, we employ non-commutative polynomial optimisation and non-commutative geometry. As a result, we present the first globally convergent methods for quantum optimal control.


title: Unit Testing for MCMC and other Monte Carlo Methods

abstract: We propose approaches for testing implementations of Markov Chain Monte Carlo methods as well as of general Monte Carlo methods. Based on statistical hypothesis tests, these approaches can be used in a unit testing framework to, for example, check if individual steps in a Gibbs sampler or a reversible jump MCMC have the desired invariant distribution. Two exact tests for assessing whether a given Markov chain has a specified invariant distribution are discussed. These and other tests of Monte Carlo methods can be embedded into a sequential method that allows low expected effort if the simulation shows the desired behavior and high power if it does not. Moreover, the false rejection probability can be kept arbitrarily low. For general Monte Carlo methods, this allows testing, for example, if a sampler has a specified distribution or if a sampler produces samples with the desired mean. The methods have been implemented in the R-package MCUnit.


title: Low-latency Cloud-based Volumetric Video Streaming Using Head Motion Prediction

abstract: Volumetric video is an emerging key technology for immersive representation of 3D spaces and objects. Rendering volumetric video requires lots of computational power which is challenging especially for mobile devices. To mitigate this, we developed a streaming system that renders a 2D view from the volumetric video at a cloud server and streams a 2D video stream to the client. However, such network-based processing increases the motion-to-photon (M2P) latency due to the additional network and processing delays. In order to compensate the added latency, prediction of the future user pose is necessary. We developed a head motion prediction model and investigated its potential to reduce the M2P latency for different look-ahead times. Our results show that the presented model reduces the rendering errors caused by the M2P latency compared to a baseline system in which no prediction is performed.


title: 1/N Expansion for Horava-Lifshitz like four-fermion models

abstract: We study a class of four-fermion Gross-Neveu like models in four dimensions with critical exponents $z=2$ and $z=3$. The models with $z=2$ are known to be perturbatively nonrenormalizable but are shown to be renormalizable in the context of the $1/N$ expansion. We calculate explicitly the effective potential for these models.


title: Passively parallel regularized stokeslets

abstract: Stokes flow, discussed by G.G. Stokes in 1851, describes many microscopic biological flow phenomena, including cilia-driven transport and flagellar motility; the need to quantify and understand these flows has motivated decades of mathematical and computational research. Regularized stokeslet methods, which have been used and refined over the past twenty years, offer significant advantages in simplicity of implementation, with a recent modification based on nearest-neighbour interpolation providing significant improvements in efficiency and accuracy. Moreover this method can be implemented with the majority of the computation taking place through built-in linear algebra, entailing that state-of-the-art hardware and software developments in the latter, in particular multicore and GPU computing, can be exploited through minimal modifications ('passive parallelism') to existing MATLAB computer code. Hence, and with widely-available GPU hardware, significant improvements in the efficiency of the regularized stokeslet method can be obtained. The approach is demonstrated through computational experiments on three model biological flows: undulatory propulsion of multiple C. Elegans, simulation of progression and transport by multiple sperm in a geometrically confined region, and left-right symmetry breaking particle transport in the ventral node of the mouse embryo. In general an order-of-magnitude improvement in efficiency is observed. This development further widens the complexity of biological flow systems that are accessible without the need for extensive code development or specialist facilities.


title: Application of peaks theory to the abundance of primordial black holes

abstract: We consider the application of peaks theory to the calculation of the number density of peaks relevant for primordial black hole (PBH) formation. For PBHs, the final mass is related to the amplitude and scale of the perturbation from which it forms, where the scale is defined as the scale at which the compaction function peaks. We therefore extend peaks theory to calculate not only the abundance of peaks of a given amplitude, but peaks of a given amplitude and scale. A simple fitting formula is given in the high-peak limit relevant for PBH formation. We also adapt the calculation to use a Gaussian smoothing function, ensuring convergence regardless of the choice of power spectrum.


title: Companion unit lower Hessenberg matrices

abstract: In recent years there has been a growing interest in companion matrices. There is a deep knowledge of sparse companion matrices, in particular it is known that every sparse companion matrix can be transformed into a unit lower Hessenberg matrix of a particularly simple type by any combination of transposition, permutation similarity and diagonal similarity. The latter is not true for the companion matrices that are non-sparse, although it is known that every non-sparse companion matrix is nonderogatory. In this work the non-sparse companion matrices that are unit lower Hessenberg will be described. A natural generalization is also considered.


title: Learning Sparse Classifiers: Continuous and Mixed Integer Optimization Perspectives

abstract: We consider a discrete optimization formulation for learning sparse classifiers, where the outcome depends upon a linear combination of a small subset of features. Recent work has shown that mixed integer programming (MIP) can be used to solve (to optimality) $\ell_0$-regularized regression problems at scales much larger than what was conventionally considered possible. Despite their usefulness, MIP-based global optimization approaches are significantly slower compared to the relatively mature algorithms for $\ell_1$-regularization and heuristics for nonconvex regularized problems. We aim to bridge this gap in computation times by developing new MIP-based algorithms for $\ell_0$-regularized classification. We propose two classes of scalable algorithms: an exact algorithm that can handle $p\approx 50,000$ features in a few minutes, and approximate algorithms that can address instances with $p\approx 10^6$ in times comparable to the fast $\ell_1$-based algorithms. Our exact algorithm is based on the novel idea of \textsl{integrality generation}, which solves the original problem (with $p$ binary variables) via a sequence of mixed integer programs that involve a small number of binary variables. Our approximate algorithms are based on coordinate descent and local combinatorial search. In addition, we present new estimation error bounds for a class of $\ell_0$-regularized estimators. Experiments on real and synthetic data demonstrate that our approach leads to models with considerably improved statistical performance (especially, variable selection) when compared to competing methods.


title: Gradient descent with momentum --- to accelerate or to super-accelerate?

abstract: We consider gradient descent with `momentum', a widely used method for loss function minimization in machine learning. This method is often used with `Nesterov acceleration', meaning that the gradient is evaluated not at the current position in parameter space, but at the estimated position after one step. In this work, we show that the algorithm can be improved by extending this `acceleration' --- by using the gradient at an estimated position several steps ahead rather than just one step ahead. How far one looks ahead in this `super-acceleration' algorithm is determined by a new hyperparameter. Considering a one-parameter quadratic loss function, the optimal value of the super-acceleration can be exactly calculated and analytically estimated. We show explicitly that super-accelerating the momentum algorithm is beneficial, not only for this idealized problem, but also for several synthetic loss landscapes and for the MNIST classification task with neural networks. Super-acceleration is also easy to incorporate into adaptive algorithms like RMSProp or Adam, and is shown to improve these algorithms.


title: Sign Flip Triangulations of the Amplituhedron

abstract: We present new triangulations of the $m=4$ amplituhedron relevant for scattering amplitudes in planar $\mathcal{N}=4$ super-Yang-Mills, obtained directly from the combinatorial definition of the geometry. Using the "sign flip" characterization of the amplituhedron, we reproduce the canonical forms for the all-multiplicity next-to-maximally helicity violating (NMHV) and next-to-next-to-maximally helicity violating ($\text{N}^2$MHV) tree-level as well as the NMHV one-loop cases, without using any input from traditional amplitudes methods. Our results provide strong evidence for the equivalence of the original definition of the amplituhedron \cite{Arkani-Hamed:2013jha} and the topological one \cite{Arkani-Hamed:2017vfh}, and suggest a new path forward for computing higher loop amplitudes geometrically. In particular, we realize the NMHV one-loop amplituhedron as the intersection of two amplituhedra of lower dimensionality, which is reflected in the novel structure of the corresponding canonical form.


title: Scaled Projected-Directions Methods with Application to Transmission Tomography

abstract: Statistical image reconstruction in X-Ray computed tomography yields large-scale regularized linear least-squares problems with nonnegativity bounds, where the memory footprint of the operator is a concern. Discretizing images in cylindrical coordinates results in significant memory savings, and allows parallel operator-vector products without on-the-fly computation of the operator, without necessarily decreasing image quality. However, it deteriorates the conditioning of the operator. We improve the Hessian conditioning by way of a block-circulant scaling operator and we propose a strategy to handle nondiagonal scaling in the context of projected-directions methods for bound-constrained problems. We describe our implementation of the scaling strategy using two algorithms: TRON, a trust-region method with exact second derivatives, and L-BFGS-B, a linesearch method with a limited-memory quasi-Newton Hessian approximation. We compare our approach with one where a change of variable is made in the problem. On two reconstruction problems, our approach converges faster than the change of variable approach, and achieves much tighter accuracy in terms of optimality residual than a first-order method.


title: A back-end, CMOS compatible ferroelectric Field Effect Transistor for synaptic weights

abstract: Neuromorphic computing architectures enable the dense co-location of memory and processing elements within a single circuit. This co-location removes the communication bottleneck of transferring data between separate memory and computing units as in standard von Neuman architectures for data-critical applications including machine learning. The essential building blocks of neuromorphic systems are non-volatile synaptic elements such as memristors. Key memristor properties include a suitable non-volatile resistance range, continuous linear resistance modulation and symmetric switching. In this work, we demonstrate voltage-controlled, symmetric and analog potentiation and depression of a ferroelectric Hf$_{57}$Zr$_{43}$O$_{2}$ (HZO) field effect transistor (FeFET) with good linearity. Our FeFET operates with a low writing energy (fJ) and fast programming time (40 ns). Retention measurements have been done over 4-bits depth with low noise (1%) in the tungsten oxide (WO$_{x}$) read out channel. By adjusting the channel thickness from 15nm to 8nm, the on/off ratio of the FeFET can be engineered from 1% to 200% with an on-resistance ideally >100 kOhm, depending on the channel geometry. The device concept is using earth-abundant materials, and is compatible with a back end of line (BEOL) integration into complementary metal-oxidesemiconductor (CMOS) processes. It has therefore a great potential for the fabrication of high density, large-scale integrated arrays of artificial analog synapses.


title: LASCA: Learning Assisted Side Channel Delay Analysis for Hardware Trojan Detection

abstract: In this paper, we introduce a Learning Assisted Side Channel delay Analysis (LASCA) methodology for Hardware Trojan detection. Our proposed solution, unlike the prior art, does not require a Golden IC. Instead, it trains a Neural Network to act as a process tracking watchdog for correlating the static timing data (produced at design time) to the delay information obtained from clock frequency sweeping (at test time) for the purpose of Trojan detection. Using the LASCA flow, we detect close to 90% of Hardware Trojans in the simulated scenarios.


title: Bayesian Inference for Big Spatial Data Using Non-stationary Spectral Simulation

abstract: It is increasingly understood that the assumption of stationarity is unrealistic for many spatial processes. In this article, we combine dimension expansion with a spectral method to model big non-stationary spatial fields in a computationally efficient manner. Specifically, we use Mejia and Rodriguez-Iturbe (1974)'s spectral simulation approach to simulate a spatial process with a covariogram at locations that have an expanded dimension. We introduce Bayesian hierarchical modelling to dimension expansion, which originally has only been modeled using a method of moments approach. In particular, we simulate from the posterior distribution using a collapsed Gibbs sampler. Our method is both full rank and non-stationary, and can be applied to big spatial data because it does not involve storing and inverting large covariance matrices. Additionally, we have fewer parameters than many other non-stationary spatial models. We demonstrate the wide applicability of our approach using a simulation study, and an application using ozone data obtained from the National Aeronautics and Space Administration (NASA).


title: On the Smith classes, the van Kampen obstruction and embeddability of $[3]*K$

abstract: In this survey-research paper, we first introduce the theory of Smith classes of complexes with fixed-point free, periodic maps on them. These classes, when defined for the deleted product of a simplicial complex $K$, are the same as the embedding classes of $K$. Embedding classes, in turn, are generalizations of the van Kampen obstruction class for embeddability of a $d$-dimensional complex $K$ into the Euclidean $2d$-space. All of these concepts will be introduced in simple terms. Second, we use the theory introduced in the first part to relate the embedding classes (or the special Smith classes) of the the complex $[3]*K$ with the embedding classes of $K$. Here $[3]*K$ is the join of $K$ with a set of three points. Specifically, we prove that if the $m$-th embedding class of $K$ is non-zero, then the $(m+2)$-nd embedding class of $[3]*K$ is non-zero. We also prove some of the consequences of this theorem for the embeddability of $[3]*K$.


title: Unsupervised Learning of Camera Pose with Compositional Re-estimation

abstract: We consider the problem of unsupervised camera pose estimation. Given an input video sequence, our goal is to estimate the camera pose (i.e. the camera motion) between consecutive frames. Traditionally, this problem is tackled by placing strict constraints on the transformation vector or by incorporating optical flow through a complex pipeline. We propose an alternative approach that utilizes a compositional re-estimation process for camera pose estimation. Given an input, we first estimate a depth map. Our method then iteratively estimates the camera motion based on the estimated depth map. Our approach significantly improves the predicted camera motion both quantitatively and visually. Furthermore, the re-estimation resolves the problem of out-of-boundaries pixels in a novel and simple way. Another advantage of our approach is that it is adaptable to other camera pose estimation approaches. Experimental analysis on KITTI benchmark dataset demonstrates that our method outperforms existing state-of-the-art approaches in unsupervised camera ego-motion estimation.


title: Acyclic complexes and Gorenstein rings

abstract: For a given class of modules $\mathcal{A}$, we denote by $\widetilde{\mathcal{A}}$ the class of exact complexes $X$ having all cycles in $\mathcal{A}$, and by $dw(\mathcal{A})$ the class of complexes $Y$ with all components $Y_j$ in $\mathcal{A}$. We use the notations $\mathcal{GI}$ $(\mathcal{GF}, \mathcal{GP})$ for the class of Gorenstein injective (Gorenstein flat, Gorenstein projective respectively) $R$-modules, $\mathcal{DI}$ for Ding injective modules, and $\mathcal{PGF}$ for projectively coresolved Gorenstein flat modules (see section 2 for definitions). We prove that the following are equivalent over any ring $R$: (1) Every exact complex of injective modules is totally acyclic. (2) Every exact complex of Gorenstein injective modules is in $\widetilde{\mathcal{GI}}$. (3) Every complex in $dw(\mathcal{GI})$ is dg-Gorenstein injective. We show that the analogue result for complexes of flat and Gorenstein flat modules also holds over arbitrary rings. if moreover, the ring is $n$-perfect for some integer $n \ge 0$, then the three equivalent statements for flat and Gorenstein flat modules are also equivalent with their counterparts for projective and projectively coresolved Gorenstein flat modules. We also prove the following characterization of Gorenstein rings: Let $R$ be a commutative coherent ring. The following statements are equivalent: (1) every exact complex of FP-injective modules has all its cycles Ding injective modules. (2) every exact complex of injectives has all its cycles Ding injective modules and every $R$-module M such that $M^+$ is Gorenstein flat is Ding injective. If moreover the ring $R$ has finite Krull dimension then statements (1), (2) above are also equivalent to (3) $R$ is a Gorenstein ring (in the sense of Iwanaga).


title: A scaling limit for the length of the longest cycle in a sparse random digraph

abstract: We discuss the length $\vec{L}_{c,n}$ of the longest directed cycle in the sparse random digraph $D_{n,p},p=c/n$, $c$ constant. We show that for large $c$ there exists a function $\vec{f}(c)$ such that $\vec{L}_{c,n}/n\to \vec{f}(c)$ a.s. The function $\vec{f}(c)=1-\sum_{k=1}^\infty p_k(c)e^{-kc}$ where $p_k$ is a polynomial in $c$. We are only able to explicitly give the values $p_1,p_2$, although we could in principle compute any $p_k$.


title: Solution Bounds, Stability and Estimation of Trapping/Stability Regions of Some Nonlinear Time-Varying Systems

abstract: Estimation of solution norms and stability for time-dependent nonlinear systems is ubiquitous in numerous engineering, natural science and control problems. Yet, practically valuable results are rare in this area. This paper develops a novel approach, which bounds the solution norms, derives the corresponding stability criteria, and estimates the trapping/stability regions for some nonautonomous and nonlinear systems, which arise in various application domains. Our inferences rest on deriving a scalar differential inequality for the norms of solutions to the initial systems. Utility of the Lipschitz inequality linearizes the associated auxiliary differential equation and yields both the upper bounds for the norms of solutions and the relevant stability criteria. To refine these inferences, we introduce a nonlinear extension of the Lipschitz inequality, which improves the developed bounds and allows estimation of the stability basins and trapping regions for the corresponding systems. Finally, we confirm the theoretical results in representative simulations.


title: Estimation of Causal Effects of Multiple Treatments in Observational Studies with a Binary Outcome

abstract: There is a dearth of robust methods to estimate the causal effects of multiple treatments when the outcome is binary. This paper uses two unique sets of simulations to propose and evaluate the use of Bayesian Additive Regression Trees (BART) in such settings. First, we compare BART to several approaches that have been proposed for continuous outcomes, including inverse probability of treatment weighting (IPTW), targeted maximum likelihood estimator (TMLE), vector matching and regression adjustment. Results suggest that under conditions of non-linearity and non-additivity of both the treatment assignment and outcome generating mechanisms, BART, TMLE and IPTW using generalized boosted models (GBM) provide better bias reduction and smaller root mean squared error. BART and TMLE provide more consistent 95 per cent CI coverage and better large-sample convergence property. Second, we supply BART with a strategy to identify a common support region for retaining inferential units and for avoiding extrapolating over areas of the covariate space where common support does not exist. BART retains more inferential units than the generalized propensity score based strategy, and shows lower bias, compared to TMLE or GBM, in a variety of scenarios differing by the degree of covariate overlap. A case study examining the effects of three surgical approaches for non-small cell lung cancer demonstrates the methods.


title: An upper bound on the Chebotarev invariant of a finite group

abstract: A subset $\{g_1, \ldots , g_d\}$ of a finite group $G$ invariably generates $G$ if the set $\{g_1^{x_1}, \ldots, g_d^{x_d}\}$ generates $G$ for every choice of $x_i \in G$. The Chebotarev invariant $C(G)$ of $G$ is the expected value of the random variable $n$ that is minimal subject to the requirement that $n$ randomly chosen elements of $G$ invariably generate $G$. The first author recently showed that $C(G)\le \beta\sqrt{|G|}$ for some absolute constant $\beta$. In this paper we show that, when $G$ is soluble, then $\beta$ is at most $5/3$. We also show that this is best possible. Furthermore, we show that, in general, for each $\epsilon>0$ there exists a constant $c_{\epsilon}$ such that $C(G)\le (1+\epsilon)\sqrt{|G|}+c_{\epsilon}$.


title: K-NN active learning under local smoothness assumption

abstract: There is a large body of work on convergence rates either in passive or active learning. Here we first outline some of the main results that have been obtained, more specifically in a nonparametric setting under assumptions about the smoothness of the regression function (or the boundary between classes) and the margin noise. We discuss the relative merits of these underlying assumptions by putting active learning in perspective with recent work on passive learning. We design an active learning algorithm with a rate of convergence better than in passive learning, using a particular smoothness assumption customized for k-nearest neighbors. Unlike previous active learning algorithms, we use a smoothness assumption that provides a dependence on the marginal distribution of the instance space. Additionally, our algorithm avoids the strong density assumption that supposes the existence of the density function of the marginal distribution of the instance space and is therefore more generally applicable.


title: Bounding the Classical Capacity of Multilevel Damping Quantum Channels

abstract: A recent method to certify the classical capacity of quantum communication channels is applied for general damping channels in finite dimension. The method compares the mutual information obtained by coding on the computational and a Fourier basis, which can be obtained by just two local measurement settings and classical optimization. The results for large representative classes of different damping structures are presented.


title: Algorithms in Multi-Agent Systems: A Holistic Perspective from Reinforcement Learning and Game Theory

abstract: Deep reinforcement learning (RL) has achieved outstanding results in recent years, which has led a dramatic increase in the number of methods and applications. Recent works are exploring learning beyond single-agent scenarios and considering multi-agent scenarios. However, they are faced with lots of challenges and are seeking for help from traditional game-theoretic algorithms, which, in turn, show bright application promise combined with modern algorithms and boosting computing power. In this survey, we first introduce basic concepts and algorithms in single agent RL and multi-agent systems; then, we summarize the related algorithms from three aspects. Solution concepts from game theory give inspiration to algorithms which try to evaluate the agents or find better solutions in multi-agent systems. Fictitious self-play becomes popular and has a great impact on the algorithm of multi-agent reinforcement learning. Counterfactual regret minimization is an important tool to solve games with incomplete information, and has shown great strength when combined with deep learning.


title: The transverse momentum spectrum of low mass Drell-Yan production at next-to-leading order in the parton branching method

abstract: It has been observed in the literature that measurements of low-mass Drell-Yan (DY) transverse momentum spectra at low center-of-mass energies $\sqrt{s}$ are not well described by perturbative QCD calculations in collinear factorization in the region where transverse momenta are comparable with the DY mass. We examine this issue from the standpoint of the Parton Branching (PB) method, combining next-to-leading-order (NLO) calculations of the hard process with the evolution of transverse momentum dependent (TMD) parton distributions. We compare our predictions with experimental measurements at low DY mass, and find very good agreement.In addition we use the low mass DY measurements at low $\sqrt{s}$ to determine the width $q_s$ of the intrinsic Gauss distribution of the PB-TMDs at low evolution scales. We find values close to what has earlier been used in applications of PB -TMDs to high-energy processes at the Large Hadron Collider (LHC) and HERA. We find that at low DY mass and low $\sqrt{s}$ even in the region of $p_t/m_{DY} \sim 1$ the contribution of multiple soft gluon emissions (included in the PB-TMDs) is essential to describe themeasurements, while at larger masses ($m_{DY} \sim m_{Z}$) and LHC energies the contribution from soft gluons in the region of $p_t/m_{DY}\sim 1$ is small.


title: Dynamical self-friction: how mass loss slows you down

abstract: We investigate dynamical self-friction, the process by which material that is stripped from a subhalo torques its remaining bound remnant, which causes it to lose orbital angular momentum. By running idealized simulations of a subhalo orbiting within an analytical host halo potential, we isolate the effect of self-friction from traditional dynamical friction due to the host halo. While at some points in a subhalo's orbit the torque of the stripped material can boost the orbital angular momentum of the remnant, the net effect over the long term is orbital decay regardless of the initial orbital parameters or subhalo mass. In order to quantify the strength of self-friction, we run a suite of simulations spanning typical host-to-subhalo mass ratios and orbital parameters. We find that the time-scale for self-friction, defined as the exponential decay time of the subhalo's orbital angular momentum, scales with mass ratio and orbital circularity similar to standard dynamical friction. The decay time due to self-friction is roughly an order of magnitude longer, suggesting that self-friction only contributes at the 10 percent level. However, along more radial orbits, self-friction can occasionally dominate over dynamical friction close to pericentric passage, where mass stripping is intense. This is also the epoch at which the self-friction torque undergoes large and rapid changes in both magnitude and direction, indicating that self-friction is an important process to consider when modeling pericentric passages of subhaloes and their associated satellite galaxies.


title: A trend in the effective spin distribution of LIGO binary black holes with mass

abstract: Binary black holes (BBHs) detected by gravitational wave (GW) observations could be broadly divided into two formation channels: those formed through field binary evolution and those assembled dynamically in dense stellar systems. Each of these formation channels, and their sub-channels, populate a distinct region in the effective spin-mass ($\chi_{\rm eff}-M$) plane. Depending on the branching ratio of different channels, an ensemble of BBHs could show a trend in this plane. Here we fit a mass-dependent distribution for $\chi_{\rm eff}$ to the GWTC-1 BBHs from the first and second observing runs of Advanced LIGO and Advanced Virgo. We find a negative correlation between mass and the mean effective spin ($\bar{\chi}_{\mathrm{eff}}$), and positive correlation with its dispersion ($\sigma_{\chi_\mathrm{eff}}$) at 75\% and 80\% confidence. This trend is robust against the choice of mass variable, but most pronounced when the mass variable is taken to be the chirp mass of the binary. The result is consistent with significant contributions from both dynamically assembled and field binaries in the GWTC-1 catalog. The upcoming LIGO O3a data release will critically test this interpretation.


title: Supermassive Star Formation via Super Competitive Accretion in Slightly Metal-enriched Clouds

abstract: Direct collapse black hole (DCBH) formation with mass $\gtrsim 10^{5}~M_{\odot}$ is a promising scenario for the origin of high-redshift supermassive black holes. It has usually been supposed that the DCBH can only form in the primordial gas since the metal enrichment enhances the cooling ability and causes the fragmentation into smaller pieces. What actually happens in such an environment, however, has not been explored in detail. Here, we study the impact of the metal enrichment on the clouds, conducting hydrodynamical simulations to follow the cloud evolution in cases with different degree of the metal enrichment $Z/Z_{\odot}=10^{-6}-10^{-3}$. Below $Z/Z_{\odot}=10^{-6}$, metallicity has no effect and supermassive stars form along with a small number of low-mass stars. With more metallicity $Z/Z_{\odot} \gtrsim 5 \times 10^{-6}$, although the dust cooling indeed promotes fragmentation of the cloud core and produces about a few thousand low-mass stars, the accreting flow preferentially feeds the gas to the central massive stars, which grows supermassive as in the primordial case. We term this formation mode as the {\it super competitive accretion}, where only the central few stars grow supermassive while a large number of other stars are competing for the gas reservoir. Once the metallicity exceeds $10^{-3}~Z_{\odot}$ and metal-line cooling becomes operative, the central star cannot grow supermassive due to lowered accretion rate. Supermassive star formation by the super competitive accretion opens up a new window for seed BHs, which relaxes the condition on metallicity and enhances the seed BH abundance.