title: Asymptotical study of two-layered discrete waveguide with a weak coupling

abstract: A thin two-layered waveguide is considered. The governing equations for this waveguide is a matrix Klein--Gordon equation of dimension~2. A formal solution of this system in the form of a double integral can be obtained by using Fourier transformation. Then, the double integral can be reduced to a single integral with the help of residue integration with respect to the time frequency. However, such an integral can be difficult to estimate since it involves branching and oscillating functions. This integral is studied asymptotically. A zone diagram technique is proposed to represent the set of possible asymptotic formulae. The zone diagram generalizes the concept of far-field and near-field zones.


title: LEKID sensitivity for space applications between 80 and 600 GHz

abstract: We report the design, fabrication and testing of Lumped Element Kinetic Inductance Detectors (LEKID) showing performance in line with the requirements of the next generation space telescopes operating in the spectral range from 80 to 600 GHz. This range is of particular interest for Cosmic Microwave Background (CMB) studies. For this purpose we have designed and fabricated 100-pixel arrays covering five distinct bands. These wafers have been measured via multiplexing, where a full array is read out using a single pair of lines. We adopted a custom cold black-body installed in front of the detectors and regulated at temperatures between 1 K and 20 K. We will describe in the present paper the main design considerations, the fabrication processes, the testing and the data analysis.


title: Precise Fabrication of Uniform Molecular Gaps for Active Nanoscale Devices

abstract: Molecules with versatile functionalities and well-defined structures, can serve as building blocks for extreme nanoscale devices. This requires their precise integration into functional heterojunctions, most commonly in the form of metal-molecule-metal architectures. Structural damage and nonuniformities caused by current fabrication techniques, however, limit their effective incorporation. Here, we present a hybrid fabrication approach enabling uniform molecular gaps. Template-stripped lithographically-patterned gold electrodes with sub-nanometer roughness are used as the bottom contacts upon which the molecular layer is formed through self-assembly. The top contacts are assembled using dielectrophoretic trapping of colloidal gold nanorods, resulting in uniform sub-5 nm junctions. In these electrically-active designs, we further explore the possibility of mechanical tunability. The presence of molecules may help control sub-nanometer mechanical modulation which is conventionally difficult to achieve due to instabilities caused by surface adhesive forces. Our approach is versatile, providing a platform to develop and study active molecular gaps towards functional nanodevices.


title: Near-infrared time-series photometry in the field of Cygnus OB2 association II. Mapping the variability of candidate members

abstract: We present the results of a J, H, and K photometric variability survey of the central 0.78 square degrees of the young OB association Cygnus OB2. We used data observed with the Wide-Field CAMera at the United Kingdom Infrared Telescope in 2007 (spanning 217 days) to investigate the light curves of 5083 low mass candidate members in the association and explore the occurrence and main characteristics of their near-infrared variability. We identified 2529 stars ($\sim$50$\%$ of the sample) with significant variability with time-scales ranging from days to months. We classified the variable stars into the following three groups according to their light curve morphology: periodic variability (1697 stars), occultation variability (124 stars), and other types of variability (726 stars). We verified that the disk-bearing stars in our sample are significantly more variable in the near-infrared than diskless stars, with a steep increase in the disk-fraction among stars with higher variability amplitude. We investigated the trajectories described by variable stars in the color-space and measured slopes for 335 stars describing linear trajectories. Based on the trajectories in the color-space, we inferred that the sample analyzed is composed of a mix of young stars presenting variability due to hot and cold spots, extinction by circumstellar material, and changes in the disk emission in the near-infrared. We contemplated using the use of near-infrared variability to identify disk-bearing stars and verified that 53.4$\%$ of the known disk-bearing stars in our sample could have been identified as such based solely on their variability. We present 18 newly identified disk-bearing stars and 14 eclipsing binary candidates among CygOB2 lower-mass members.


title: ProVe -- Self-supervised pipeline for automated product replacement and cold-starting based on neural language models

abstract: In retail vertical industries, businesses are dealing with human limitation of quickly understanding and adapting to new purchasing behaviors. Moreover, retail businesses need to overcome the human limitation of properly managing a massive selection of products/brands/categories. These limitations lead to deficiencies from both commercial (e.g. loss of sales, decrease in customer satisfaction) and operational perspective (e.g. out-of-stock, over-stock). In this paper, we propose a pipeline approach based on Natural Language Understanding, for recommending the most suitable replacements for products that are out-of-stock. Moreover, we will propose a solution for managing products that were newly introduced in a retailer's portfolio with almost no transactional history. This solution will help businesses: automatically assign the new products to the right category; recommend complementary products for cross-sell from day 1; perform sales predictions even with almost no transactional history. Finally, the vector space model resulted by applying the pipeline presented in this paper is directly used as semantic information in deep learning-based demand forecasting solutions, leading to more accurate predictions. The whole research and experimentation process have been done using real-life private transactional data, however the source code is available on https://github.com/Lummetry/ProVe


title: Anomalous Cyclotron Motion in Graphene Superlattice Cavities

abstract: We consider graphene superlattice miniband fermions probed by electronic interferometry in magneto-transport experiments. By decoding the observed Fabry-P\'erot interference patterns together with our corresponding quantum transport simulations, we find that the Dirac quasiparticles originating from the superlattice minibands do not undergo conventional cyclotron motion but follow more subtle trajectories. In particular, dynamics at low magnetic fields is characterized by peculiar, straight trajectory segments. Our results provide new insights into superlattice miniband fermions and open up novel possibilities to use periodic potentials in electron optics experiments.


title: Pullbacks of $\kappa$ classes on $\overline{\mathcal{M}}_{0,n}$

abstract: The moduli space $\overline{\mathcal{M}}_{0,n}$ carries a codimension-$d$ cycle class $\kappa_{d}$. We consider the subspace $\mathcal{K}^{d}_{n}$ of $A^d(\overline{\mathcal{M}}_{0,n},\mathbb{Q})$ spanned by pullbacks of $\kappa_d$ via forgetful maps. We find a permutation basis for $\mathcal{K}^{d}_{n}$, and describe its annihilator under the intersection pairing in terms of $d$-dimensional boundary strata. As an application, we give a new permutation basis of the divisor class group of $\overline{\mathcal{M}}_{0,n}$.


title: The Initial Mass Function in the ELT era

abstract: The initial mass function (IMF) is an important, yet enigmatic aspect of the star formation process. The two major open questions regarding the IMF are: is the IMF constant regardless of environment? Is the IMF a universal property of star formation? The next generation of extremely large telescopes will allow us to observe further, fainter and more compact stellar clusters than is possible with current facilities. In these proceeding we present our study looking at just how much will these future observatories improve our knowledge of the IMF.


title: Endogenous Treatment Effect Estimation with some Invalid and Irrelevant Instruments

abstract: Instrumental variables (IV) regression is a popular method for the estimation of the endogenous treatment effects. Conventional IV methods require all the instruments are relevant and valid. However, this is impractical especially in high-dimensional models when we consider a large set of candidate IVs. In this paper, we propose an IV estimator robust to the existence of both the invalid and irrelevant instruments (called R2IVE) for the estimation of endogenous treatment effects. This paper extends the scope of Kang et al. (2016) by considering a true high-dimensional IV model and a nonparametric reduced form equation. It is shown that our procedure can select the relevant and valid instruments consistently and the proposed R2IVE is root-n consistent and asymptotically normal. Monte Carlo simulations demonstrate that the R2IVE performs favorably compared to the existing high-dimensional IV estimators (such as, NAIVE (Fan and Zhong, 2018) and sisVIVE (Kang et al., 2016)) when invalid instruments exist. In the empirical study, we revisit the classic question of trade and growth (Frankel and Romer, 1999).


title: On the convergence of the Metropolis algorithm with fixed-order updates for multivariate binary probability distributions

abstract: The Metropolis algorithm is arguably the most fundamental Markov chain Monte Carlo (MCMC) method. But the algorithm is not guaranteed to converge to the desired distribution in the case of multivariate binary distributions (e.g., Ising models or stochastic neural networks such as Boltzmann machines) if the variables (sites or neurons) are updated in a fixed order, a setting commonly used in practice. The reason is that the corresponding Markov chain may not be irreducible. We propose a modified Metropolis transition operator that behaves almost always identically to the standard Metropolis operator and prove that it ensures irreducibility and convergence to the limiting distribution in the multivariate binary case with fixed-order updates. The result provides an explanation for the behaviour of Metropolis MCMC in that setting and closes a long-standing theoretical gap. We experimentally studied the standard and modified Metropolis operator for models were they actually behave differently. If the standard algorithm also converges, the modified operator exhibits similar (if not better) performance in terms of convergence speed.


title: A Hennessy-Milner Theorem for ATL with Imperfect Information

abstract: We show that a history-based variant of alternating bisimulation with imperfect information allows it to be related to a variant of Alternating-time Temporal Logic (ATL) with imperfect information by a full Hennessy-Milner theorem. The variant of ATL we consider has a common knowledge semantics, which requires that the uniform strategy available for a coalition to accomplish some goal must be common knowledge inside the coalition, while other semantic variants of ATL with imperfect information do not accommodate a Hennessy-Milner theorem. We also show that the existence of a history-based alternating bisimulation between two finite Concurrent Game Structures with imperfect information (iCGS) is undecidable.


title: Effects of Low Pressure Condition on Partial Discharges in WBG Power Electronics Modules

abstract: The aviation industry aims to reduce CO2 emission by reducing energy consumption and benefitting from more electrical systems than those based on fossil fuels. The more electric aircraft (MEA) can take advantage of the drives based on wide bandgap (WBG) based power modules that are lighter and can bear higher voltages and currents. However, the fast-rise and repetitive voltage pulses generated by WBG-based systems can endanger the insulating property of dielectric materials due to the partial discharges. PDs cause the local breakdown of insulation materials in regions with high electric field magnitude. In the case of power electronic modules, silicone gel, which is a predominant type of encapsulation material, is susceptible to these PDs. In this study, it is aimed to investigate the impact of rise time on various PD characteristics including PD true charge magnitude, inception and extinction fields, duration, and so forth. Besides, those systems that are anticipated to operate under harsh environmental conditions such as naval or aviation industries, may expect additional threat. Therefore, this paper puts forth the combination of low pressure conditions and fast rise, high frequency square wave pulses. The results demonstrate that more intense and more intense ones are going to be imposed at lower pressures. COMSOL Multiphysics interfaced with MATLAB is used to simulate the PD detection process based on the experimental data found in the literature.


title: Ab initio determination of magnetic ground state of pyrochlore Y$_2$Mn$_2$O$_7$

abstract: There are two discrepant experimental results on the magnetic ground state of Y$_{2}$Mn$_{2}$O$_{7}$, one study proposes a spin glass state, while another introduces the material as a ferromagnet. In this study, we attempt to resolve this issue by employing density functional theory and Monte Carlo simulations. We derive different spin models by varying the Hubbard $U$ parameter in ab initio GGA+$U$ calculations. For the most range of Hubbard $U$, We obtain that the leading terms in the spin Hamiltonian are bi-quadratic and the nearest neighbor Heisenberg exchange interactions. By comparing Monte Carlo simulations of these models with the experiments, we find a ferromagnetic ground state for Y$_{2}$Mn$_{2}$O$_{7}$ as the most compatible with experiments. We also consider Y$_{2}$Mo$_{2}$O$_{7}$ as a prototype of the defect-free pyrochlore system with spin-glass behavior and compare it with Y$_{2}$Mn$_{2}$O$_{7}$. The orbital degrees of freedom are considered as a leading factor in converting a defect-free pyrochlore such as Y$_{2}$Mn$_{2}$O$_{7}$ to a spin glass system. By changing the $d$ orbital occupations of Mo atoms, our GGA+$U$ calculations for Y$_{2}$Mo$_{2}$O$_{7}$ indicate many nearly degenerate states with different $d$ orbital orientations which reveals $d$ orbital degrees of freedom in this material. While for Y$_{2}$Mn$_{2}$O$_{7}$, we find a single ground state with a fixed orbital orientation. Consequently, all of our ab initio approaches confirm Y$_{2}$Mn$_{2}$O$_{7}$ as a ferromagnetic system.


title: Truncated moment sequences and a solution to the channel separability problem

abstract: We consider the problem of separability of quantum channels via the Choi matrix representation given by the Choi-Jamio{\l}kowski isomorphism. We explore three classes of separability across different cuts between systems and ancillae and we provide a solution based on the mapping of the coordinates of the Choi state (in a fixed basis) to a truncated moment sequence (tms) $y$. This results in an algorithm which gives a separability certificate using semidefinite programming. The computational complexity and the performance of it depend on the number of variables $n$ in the tms and on the size of the moment matrix $M_t(y)$ of order $t$. We exploit the algorithm to numerically investigate separability of families of 2-qubit and single-qutrit channels; in the latter case we can provide an answer for examples explored earlier through the criterion based on the negativity $N$, a criterion which remains inconclusive for Choi matrices with $N=0$.


title: Thermal luminosities of cooling neutron stars

abstract: Ages and thermal luminosities of neutron stars, inferred from observations, can be interpreted with the aid of the neutron star cooling theory to gain information on the properties of superdense matter in neutron-star interiors. We present a survey of estimated ages, surface temperatures and thermal luminosities of middle-aged neutron stars with relatively weak or moderately strong magnetic fields, which can be useful for these purposes. The catalogue includes results selected from the literature, supplemented with new results of spectral analysis of a few cooling neutron stars. The data are compared with the theory. We show that overall agreement of theoretical cooling curves with observations improves substantially for models where neutron superfluidity in stellar core is weak.


title: Resource Allocation via Graph Neural Networks in Free Space Optical Fronthaul Networks

abstract: This paper investigates the optimal resource allocation in free space optical (FSO) fronthaul networks. The optimal allocation maximizes an average weighted sum-capacity subject to power limitation and data congestion constraints. Both adaptive power assignment and node selection are considered based on the instantaneous channel state information (CSI) of the links. By parameterizing the resource allocation policy, we formulate the problem as an unsupervised statistical learning problem. We consider the graph neural network (GNN) for the policy parameterization to exploit the FSO network structure with small-scale training parameters. The GNN is shown to retain the permutation equivariance that matches with the permutation equivariance of resource allocation policy in networks. The primal-dual learning algorithm is developed to train the GNN in a model-free manner, where the knowledge of system models is not required. Numerical simulations present the strong performance of the GNN relative to a baseline policy with equal power assignment and random node selection.


title: Efficient Cysteine Conformer Search with Bayesian Optimization

abstract: Finding low-energy molecular conformers is challenging due to the high dimensionality of the search space and the computational cost of accurate quantum chemical methods for determining conformer structures and energies. Here, we combine active-learning Bayesian optimization (BO) algorithms with quantum chemistry methods to address this challenge. Using cysteine as an example, we show that our procedure is both efficient and accurate. After only one thousand single-point calculations and approximately thirty structure relaxations, which is less than 10% computational cost of the current fastest method, we have found the low-energy conformers in good agreement with experimental measurements and reference calculations.


title: Database Reconstruction from Noisy Volumes: A Cache Side-Channel Attack on SQLite

abstract: We demonstrate the feasibility of database reconstruction under a cache side-channel attack on SQLite. Specifically, we present a Flush+Reload attack on SQLite that obtains approximate (or "noisy") volumes of range queries made to a private database. We then present several algorithms that, taken together, reconstruct nearly the exact database in varied experimental conditions, given these approximate volumes. Our reconstruction algorithms employ novel techniques for the approximate/noisy setting, including a noise-tolerant clique-finding algorithm, a "Match & Extend" algorithm for extrapolating volumes that are omitted from the clique, and a "Noise Reduction Step" that makes use of a closest vector problem (CVP) solver to improve the overall accuracy of the reconstructed database. The time complexity of our attacks grows quickly with the size of the range of the queried attribute, but scales well to large databases. Experimental results show that we can reconstruct databases of size 100,000 and ranges of size 12 with error percentage of 0.11 % in under 12 hours on a personal laptop.


title: The Nonuniversality of Wealth Distribution Tails Near Wealth Condensation Criticality

abstract: In this work, we modify the affine wealth model of wealth distributions to examine the effects of nonconstant redistribution on the very wealthy. Previous studies of this model, restricted to flat redistribution schemes, have demonstrated the presence of a phase transition to a partially wealth-condensed state, or "partial oligarchy," at the critical value of an order parameter. These studies have also indicated the presence of an exponential tail in wealth distribution precisely at criticality. Away from criticality, the tail was observed to be Gaussian. In this work, we generalize the flat redistribution within the affine wealth model to allow for an essentially arbitrary redistribution policy. We show that the exponential tail observed near criticality in prior work is, in fact, a special case of a much broader class of critical, slower-than-Gaussian decays that depend sensitively on the corresponding asymptotic behavior of the progressive redistribution model used. We thereby demonstrate that the functional form of the tail of the wealth distribution in a near-critical society is not universal in nature but rather entirely determined by the specifics of public policy decisions. This is significant because most major economies today are observed to be near-critical.


title: A Unifying Framework for Reinforcement Learning and Planning

abstract: Sequential decision making, commonly formalized as optimization of a Markov Decision Process, is a key challenge in artificial intelligence. Two successful approaches to MDP optimization are reinforcement learning and planning, which both largely have their own research communities. However, if both research fields solve the same problem, then we might be able to disentangle the common factors in their solution approaches. Therefore, this paper presents a unifying algorithmic framework for reinforcement learning and planning (FRAP), which identifies underlying dimensions on which MDP planning and learning algorithms have to decide. At the end of the paper, we compare a variety of well-known planning, model-free and model-based RL algorithms along these dimensions. Altogether, the framework may help provide deeper insight in the algorithmic design space of planning and reinforcement learning.


title: Effect of population migration and punctuated lockdown on the spread of infectious diseases

abstract: One of the critical measures to control infectious diseases is a lockdown. Once past the lockdown stage in many parts of the world, the crucial question now concerns the effects of relaxing the lockdown and finding the best ways to implement further lockdown(s), if required, to control the spread. With the relaxation of lockdown, people migrate to different cities and enhance the spread of the disease. This work presents the population migration model for n-cities and applies the model for migration between two and three cities. The reproduction number is calculated, and the effect of the migration rate is analyzed. A punctuated lockdown is implemented to simulate a protocol of repeated lockdowns that limits the resurgence of infections. A damped oscillatory behavior is observed with multiple peaks over a period.


title: Non-Local Thermodynamic Equilibrium Transmission Spectrum Modelling of HD209458b

abstract: Context - Exoplanetary upper atmospheres are low density environments where radiative processes can compete with collisional ones and introduce non-local thermodynamic equilibrium (NLTE) effects into transmission spectra. Aims - We develop a NLTE radiative transfer framework capable of modelling exoplanetary transmission spectra over a wide range of planetary properties. Methods - We adapt the NLTE spectral synthesis code Cloudy to produce an atmospheric structure and atomic transmission spectrum in both NLTE and local thermodynamic equilibrium (LTE) for the hot Jupiter HD209458b, given a published T-P profile and assuming solar metallicity. Selected spectral features, including H$\alpha$, Na I D, He I $\lambda$10830, Fe I & II ultra-violet (UV) bands, and C, O and Si UV lines, are compared with literature observations and models where available. The strength of NLTE effects are measured for individual spectral lines to identify which features are most strongly affected. Results - The developed modelling framework computing NLTE synthetic spectra reproduces literature results for the He I $\lambda$10830 triplet, the Na I D lines, and the forest of Fe I lines in the optical. Individual spectral lines in the NLTE spectrum exhibit up to 40 % stronger absorption relative to the LTE spectrum.


title: An unsupervised deep learning approach in solving partial integro-differential equations

abstract: We investigate solving partial integro-differential equations (PIDEs) using unsupervised deep learning in this paper. To price options, assuming underlying processes follow Levy processes, we require to solve PIDEs. In supervised deep learning, pre-calculated labels are used to train neural networks to fit the solution of the PIDE. In an unsupervised deep learning, neural networks are employed as the solution, and the derivatives and the integrals in the PIDE are calculated based on the neural network. By matching the PIDE and its boundary conditions, the neural network gives an accurate solution of the PIDE. Once trained, it would be fast for calculating options values as well as option Greeks.


title: Magnetic reconnection and plasmoid formation driven by the two-dimensional resistive kink instability in a current hole configuration

abstract: We investigate the nonlinear evolution of the $m = 1$ resistive internal kink mode in a two-dimensional (2D) configuration containing a central region of negative current density, also known as the "current hole" setup. The finite-element code FINMHD is used to solve a reduced set of incompressible Magnetohydrodynamic (MHD) equations with a current-vorticity formulation. First, the kink instability linearly develops in agreement with the general theory of resistive internal kink mode, and it subsequently leads to the formation of a current sheet. At relatively low Lundquist number $S$, a magnetic reconnection process proceeds with a rate predicted by the Sweet-Parker regime. Conversely, when $S$ exceeds a critical value that is $S_c \simeq 10^4$, the current sheet is disrupted by the formation of plasmoids on a slightly sub-Alfv\'enic time scale. In the latter case, a stochastic reconnection regime exhibiting Petshek-type features enriched by plasmoids is reached. A relatively fast normalized reconnection rate value of order $0.02$ is also measured. Finally, we compare our results with those obtained in similar 2D previous studies using ideal MHD instabilities to initiate the process, and discuss their relevance for the general theory of plasmoid chains formation and associated fast reconnection regime.


title: Microwave diagnostics of magnetic field strengths in solar flaring loops

abstract: We have performed microwave diagnostics of the magnetic field strengths in solar flare loops based on the theory of gyrosynchrotron emission. From Nobeyama Radioheliograph observations of three flare events at 17 and 34 GHz, we obtained the degree of circular polarization and the spectral index of microwave flux density, which were then used to map the magnetic field strengths in post-flare loops. Our results show that the magnetic field strength typically decreases from ~800 G near the loop footpoints to ~100 G at a height of 10--25 Mm. Comparison of our results with magnetic field modeling using a flux rope insertion method is also discussed. Our study demonstrates the potential of microwave imaging observations, even at only two frequencies, in diagnosing the coronal magnetic field of flaring regions.


title: SASO: Joint 3D Semantic-Instance Segmentation via Multi-scale Semantic Association and Salient Point Clustering Optimization

abstract: We propose a novel 3D point cloud segmentation framework named SASO, which jointly performs semantic and instance segmentation tasks. For semantic segmentation task, inspired by the inherent correlation among objects in spatial context, we propose a Multi-scale Semantic Association (MSA) module to explore the constructive effects of the semantic context information. For instance segmentation task, different from previous works that utilize clustering only in inference procedure, we propose a Salient Point Clustering Optimization (SPCO) module to introduce a clustering procedure into the training process and impel the network focusing on points that are difficult to be distinguished. In addition, because of the inherent structures of indoor scenes, the imbalance problem of the category distribution is rarely considered but severely limits the performance of 3D scene perception. To address this issue, we introduce an adaptive Water Filling Sampling (WFS) algorithm to balance the category distribution of training data. Extensive experiments demonstrate that our method outperforms the state-of-the-art methods on benchmark datasets in both semantic segmentation and instance segmentation tasks.


title: Geometry of the X-ray source 1H 0707-495

abstract: We investigate constraints for the size and location of the X-ray source in 1H 0707-495 determined from the shape of the relativistically smeared reflection from the accretion disc. We develop a new code to model an extended X-ray source and we apply it to all archival XMM observations of 1H 0707-495. Contrary to Wilkins et al. we find that the relativistic reflection in this source is not consistent with an extended uniform corona. Instead, we find that the X-ray source must be very compact, with the size of at most a gravitational radius, and located at most at a few gravitational radii from the black hole horizon. A uniform extended corona does indeed produce an emissivity which is like a twice broken power law, but the inner emissivity is fixed by the source geometry rather than being a free parameter. In 1H0707-495, reflection from the inner disc is much stronger than expected for a uniformly extended source. Including the effect of ionised absorption from a wind does not change this conclusion, but including scattered emission (and more complex absorption) from the wind can dramatically change the reflection parameters.


title: Observational Constraints on the Regularized 4D Einstein-Gauss-Bonnet Theory of Gravity

abstract: In this paper we study the observational constraints that can be imposed on the coupling parameter, $\hat \alpha$, of the regularized version of the 4-dimensional Einstein-Gauss-Bonnet theory of gravity. We use the scalar-tensor field equations of this theory to perform a thorough investigation of its slow-motion and weak-field limit, and apply our results to observations of a wide array of physical systems that admit such a description. We find that the LAGEOS satellites are the most constraining, requiring $| \hat \alpha | \lesssim 10^{10} \,{\rm m}^2$. This constraint suggests that the possibility of large deviations from general relativity is small in all systems except the very early universe ($t<10^{-3}\, {\rm s}$), or the immediate vicinity of stellar-mass black holes ($M\lesssim100\, M_{\odot}$). We then consider constraints that can be imposed on this theory from cosmology, black hole systems, and table-top experiments. It is found that early universe inflation prohibits all but the smallest negative values of $\hat \alpha$, while observations of binary black hole systems are likely to offer the tightest constraints on positive values, leading to overall bounds $0 \lesssim \hat \alpha \lesssim 10^8 \, {\rm m}^2$.


title: The clustering dynamics of primordial black boles in $N$-body simulations

abstract: We explore the possibility that Dark Matter (DM) may be explained by a non-uniform background of approximately stellar-mass clusters of Primordial Black Holes (PBHs), by simulating the evolution them from recombination to the present with over 5000 realisations using a Newtonian $ N $-body code. We compute the cluster rate of evaporation, and extract the binary and merged sub-populations along with their parent and merger tree histories, lifetimes and formation rates; the dynamical and orbital parameter profiles, the degree of mass segregation and dynamical friction, and power spectrum of close encounters. Overall, we find that PBHs can constitute a viable DM candidate, and that their clustering presents a rich phenomenology throughout the history of the Universe. We show that binary systems constitute about 9.5\% of all PBHs at present, with mass ratios of $ \bar{q}_{\rm B} = 0.154 $, and total masses of $ \bar{m}_{\rm T,\,B} = 303\,M_\odot$. Merged PBHs are rare, about 0.0023\% of all PBHs at present, with mass ratios of $ \bar{q}_{\rm B}= 0.965 $ with total and chirp masses of $ \bar{m}_{\rm T,\,B}= 1670\,M_\odot$ and $ \bar{m}_{c,{\rm M}} = 642\,M_\odot $ respectively. We find that cluster puffing up and evaporation leads to bubbles of these PBHs of order 1 kpc containing at present times about 36\% of objects and mass, with hundred pc sized cores. We also find that these PBH sub-haloes are distributed in wider PBH haloes of order hundreds of kpc, containing about 63\% of objects and mass, coinciding with the sizes of galactic halos. We find at last high rates of close encounters of massive Black Holes ($ M \sim 1000\,M_\odot$), with $ \Gamma^{\mathrm{S}} = (1.2^{+5.9}_{-0.9}) \times 10^{7} \mathrm{yr^{-1} Gpc^{-3}}$ and mergers with $\Gamma^{\mathrm{M}} = 1337 \pm 41 \mathrm{yr^{-1} Gpc^{-3}} $.


title: Angular momentum alignment-to-orientation conversion in the ground state of Rb atoms at room temperature

abstract: We investigated experimentally and theoretically angular momentum alignment-to-orientation conversion created by the joint interaction of laser radiation and an external magnetic field with atomic rubidium at room temperature. In particular we were interested in alignment-to-orientation conversion in atomic ground state. Experimentally the laser frequency was fixed to the hyperfine transitions of $D_1$ line of rubidium. We used a theoretical model for signal simulations that takes into account all neighboring hyperfine levels, the mixing of magnetic sublevels in an external magnetic field, the coherence properties of the exciting laser radiation, and the Doppler effect. The experiments were carried out by exciting the atoms with linearly polarized laser radiation. Two oppositely circularly polarized laser induced fluorescence (LIF) components were detected and afterwards their difference was taken. The combined LIF signals originating from the hyperfine magnetic sublevel transitions of $^{85}$Rb and $^{87}$Rb rubidium isotopes were included. The alignment-to-orientation conversion can be undoubtedly identified in the difference signals for various laser frequencies as well as change in signal shapes can be observed when the laser power density is increased. We studied the formation and the underlying physical processes of the observed signal of the LIF components and their difference by performing the analysis of the influence of incoherent and coherent effects. We performed simulations of theoretical signals that showed the influence of ground-state coherent effects on the LIF difference signal.


title: Pre-training via Paraphrasing

abstract: We introduce MARGE, a pre-trained sequence-to-sequence model learned with an unsupervised multi-lingual multi-document paraphrasing objective. MARGE provides an alternative to the dominant masked language modeling paradigm, where we self-supervise the reconstruction of target text by retrieving a set of related texts (in many languages) and conditioning on them to maximize the likelihood of generating the original. We show it is possible to jointly learn to do retrieval and reconstruction, given only a random initialization. The objective noisily captures aspects of paraphrase, translation, multi-document summarization, and information retrieval, allowing for strong zero-shot performance on several tasks. For example, with no additional task-specific training we achieve BLEU scores of up to 35.8 for document translation. We further show that fine-tuning gives strong performance on a range of discriminative and generative tasks in many languages, making MARGE the most generally applicable pre-training method to date.


title: On the generalizability of artificial neural networks in spin models

abstract: The applicability of artificial neural networks (ANNs) is typically limited to the models they are trained with and little is known about their generalizability, which is a pressing issue in the practical application of trained ANNs to unseen problems. Here, by using the task of identifying phase transitions in spin models, we establish a systematic generalizability such that simple ANNs trained with the two-dimensional ferromagnetic Ising model can be applied to the ferromagnetic $q$-state Potts model in different dimensions for $q \geq 2$. The same scheme can be applied to the highly nontrivial antiferromagnetic $q$-state Potts model. We demonstrate that similar results can be obtained by reducing the exponentially large state space spanned by the training data to one that comprises only three representative configurations artificially constructed through symmetry considerations. We expect our findings to simplify and accelerate the development of machine learning-assisted tasks in spin-model related disciplines in physics and materials science.

categories: cond-mat.dis-nn cond-mat.stat-mech hep-lat hep-th


title: Cosmological homogeneity scale estimates are dressed

abstract: We investigate number count statistics as measures for transition to homogeneity of the matter distribution in the Universe and analyse how such statistics might be `dressed' by the assumed survey selection function. Since the estimated survey selection function -- which ideally accounts for selection bias in the observed distribution -- is partially degenerate with the estimated underlying distribution of galaxies, the ability to identify the correct survey selection function is of importance for obtaining reliable estimates for clustering statistics. Selection functions of existing galaxy catalogues are modelled from data to resemble the redshift distribution and mean density of the observed galaxies. Proposed estimates of the selection function for upcoming surveys in addition use the angular distribution of galaxies to generate the angular selection function instead of using angular completeness estimates. We argue that such modelling of the selection function could potentially underestimate the deviance from homogeneity at scales probed by existing catalogues. We investigate the impact of conventionally applied methods for estimation of the survey selection function on number count in sphere statistics in a toy model setting. The example density distribution is asymptotically homogeneous, while non-linear density fluctuations are present regionally. We find that density oscillations with period comparable to characteristic scales of the survey are suppressed when conventional estimates of the survey selection function are invoked, resulting in number count statistics that are biased towards homogeneity. For our concrete toy model with maximum density contrasts of 1 and period of the density oscillation comparable in size to the survey radius, we find that the homogeneity scale is underestimated by ~40%, however this quantitative result is dependent on the model setup.


title: The impact of distant fly-bys on the rate of binary primordial black hole mergers

abstract: By performing Monte Carlo simulations of the evolution of binary primordial black hole (PBH) systems, we estimate the effect of distant encounters with single PBHs upon the coalescence time and merger rate of binary PBHs. We find that, for models where PBHs compose a large fraction of dark matter, $f_\mathrm{PBH}\sim 1$, the expected fractional change in coalescence time is negligible, of order $10^{-6}$ for most binaries. For models with significantly lower PBH abundances, $f_\mathrm{PBH}\ll 1$, we find that the average change in binary lifetime due to encounters can be as large as $\mathcal{O}(10^{-2})$, with a small number of binaries experiencing an order unity change in lifetime. In the absence of encounters, we also compare the use of an analytic approximation for the coalescence time to numerically evolving the binary system, finding that the analytic approximation results in an order $10\%$ error in the coalescence time. However, when these effects are taken into consideration, there is a negligible change to the calculated merger rate, placing previous constraints on the PBH abundance arising from observed gravitational wave signals from merging binary black holes on a more secure footing.


title: Computing all $s$-$t$ bridges and articulation points simplified

abstract: Given a directed graph $G$ and a pair of nodes $s$ and $t$, an $s$-$t$ bridge of $G$ is an edge whose removal breaks all $s$-$t$ paths of $G$. Similarly, an $s$-$t$ articulation point of $G$ is a node whose removal breaks all $s$-$t$ paths of $G$. Computing the sequence of all $s$-$t$ bridges of $G$ (as well as the $s$-$t$ articulation points) is a basic graph problem, solvable in linear time using the classical min-cut algorithm. When dealing with cuts of unit size ($s$-$t$ bridges) this algorithm can be simplified to a single graph traversal from $s$ to $t$ avoiding an arbitrary $s$-$t$ path, which is interrupted at the $s$-$t$ bridges. Further, the corresponding proof is also simplified making it independent of the theory of network flows.


title: In situ synthesis of polyynes in a polymer matrix by pulsed laser ablation in liquid

abstract: Polyynes are finite chains formed by sp-hybridized carbon atoms with alternating single and triple bonds and displaying intriguing electronic and optical properties. Pulsed laser ablation in liquid (PLAL) is a well assessed technique for the physical synthesis of hydrogen-capped polyynes in solution, however, their limited stability prevents further exploitation in materials for different applications. In this work, polyynes in poly(vinyl alcohol) (PVA) were produced in a single-step PLAL process by ablating graphite directly in aqueous solution of PVA, investigating the role of polymer concentration. The presence of PVA solution, as a participating medium for PLAL, is shown to favour the formation of polyynes. The addition of Ag colloids to the aqueous PVA/polyynes solution allowed surface-enhanced Raman spectroscopy (SERS) measurements, carried out both on liquid samples and on free-standing nanocomposites, obtained after solvent evaporation. We show that polyynes in the nanocomposite remain stable at least for 11 months, whereas the corresponding PVA/Ag/polyynes solution displayed a strong polyyne reduction already after 3 weeks. These results open the view to further characterizations of the properties of polyyne-based films and materials.


title: Topological entanglement properties of disconnected partitions in the Su-Schrieffer-Heeger model

abstract: We study the disconnected entanglement entropy, $S^D$, of the Su-Schrieffer-Heeger model. $S^D$ is a combination of both connected and disconnected bipartite entanglement entropies that removes all area and volume law contributions, and is thus only sensitive to the non-local entanglement stored within the ground state manifold. Using analytical and numerical computations, we show that $S^D$ behaves as a topological invariant, i.e., it is quantized to either $0$ or $2 \log (2)$ in the topologically trivial and non-trivial phases, respectively. These results also hold in the presence of symmetry-preserving disorder. At the second-order phase transition separating the two phases, $S^D$ displays a system-size scaling behavior akin to those of conventional order parameters, that allows us to compute entanglement critical exponents. To corroborate the topological origin of the quantized values of $S^D$, we show how the latter remain quantized after applying unitary time evolution in the form of a quantum quench, a characteristic feature of topological invariants.


title: Deep-learning Autoencoder for Coherent and Nonlinear Optical Communication

abstract: Motivated by the recent success of end-to-end training of communications in the wireless domain, we strive to adapt the end-to-end-learning idea from the wireless case (i.e., linear) to coherent optical fiber links (i.e., nonlinear). Although, at first glance, it sounds like a straightforward extension, it turns out that several pitfalls exist - in terms of theory but also in terms of practical implementation. This paper analyzes the potential of an autoencoder and limitations for the optical fiber under the influence of Kerr-nonlinearity and chromatic dispersion. As there is no exact capacity limit known and, hence, no analytical perfect system solution available, we set great value to the interpretability on the learnings of the autoencoder. Therefore, we design its architecture to be as close as possible to the structure of a classic communication system, knowing that this may limit its degree of freedom and, thus, its performance. Nevertheless, we were able to achieve an unexpected high gain in terms of spectral efficiency compared to a conventional reference system.


title: The low-luminosity type II SN\,2016aqf: A well-monitored spectral evolution of the Ni/Fe abundance ratio

abstract: Low-luminosity type II supernovae (LL SNe~II) make up the low explosion energy end of core-collapse SNe, but their study and physical understanding remain limited. We present SN\,2016aqf, a LL SN~II with extensive spectral and photometric coverage. We measure a $V$-band peak magnitude of $-14.58$\,mag, a plateau duration of $\sim$100\,days, and an inferred $^{56}$Ni mass of $0.008 \pm 0.002$\,\msun. The peak bolometric luminosity, L$_{\rm bol} \approx 10^{41.4}$\,erg\,s$^{-1}$, and its spectral evolution is typical of other SNe in the class. Using our late-time spectra, we measure the [\ion{O}{i}] $\lambda\lambda6300, 6364$ lines, which we compare against SN II spectral synthesis models to constrain the progenitor zero-age main-sequence mass. We find this to be 12 $\pm$ 3\,\msun. Our extensive late-time spectral coverage of the [\ion{Fe}{ii}] $\lambda7155$ and [\ion{Ni}{ii}] $\lambda7378$ lines permits a measurement of the Ni/Fe abundance ratio, a parameter sensitive to the inner progenitor structure and explosion mechanism dynamics. We measure a constant abundance ratio evolution of $0.081^{+0.009}_{-0.010}$, and argue that the best epochs to measure the ratio are at $\sim$200 -- 300\,days after explosion. We place this measurement in the context of a large sample of SNe II and compare against various physical, light-curve and spectral parameters, in search of trends which might allow indirect ways of constraining this ratio. We do not find correlations predicted by theoretical models; however, this may be the result of the exact choice of parameters and explosion mechanism in the models, the simplicity of them and/or primordial contamination in the measured abundance ratio.


title: Computing the elastic mechanical properties of rod-like DNA nanostructures

abstract: To study the elastic properties of rod-like DNA nanostructures, we perform long simulations of these structure using the oxDNA coarse-grained model. By analysing the fluctuations in these trajectories we obtain estimates of the bend and twist persistence lengths, and the underlying bend and twist elastic moduli and couplings between them. Only on length scales beyond those associated with the spacings between the interhelix crossovers do the bending fluctuations behave like those of a worm-like chain. The obtained bending persistence lengths are much larger than that for double-stranded DNA and increase non-linearly with the number of helices, whereas the twist moduli increase approximately linearly. To within the numerical error in our data, the twist-bend coupling constants are of order zero. That the bending persistence lengths we obtain are generally somewhat higher than in experiment probably reflects both that the simulated origami have no assembly defects and that the oxDNA extensional modulus for double-stranded DNA is too large.


title: Deriving information from missing data: implications for mood prediction

abstract: The availability of mobile technologies has enabled the efficient collection prospective longitudinal, ecologically valid self-reported mood data from psychiatric patients. These data streams have potential for improving the efficiency and accuracy of psychiatric diagnosis as well predicting future mood states enabling earlier intervention. However, missing responses are common in such datasets and there is little consensus as to how this should be dealt with in practice. A signature-based method was used to capture different elements of self-reported mood alongside missing data to both classify diagnostic group and predict future mood in patients with bipolar disorder, borderline personality disorder and healthy controls. The missing-response-incorporated signature-based method achieves roughly 66\% correct diagnosis, with f1 scores for three different clinic groups 59\% (bipolar disorder), 75\% (healthy control) and 61\% (borderline personality disorder) respectively. This was significantly more efficient than the naive model which excluded missing data. Accuracies of predicting subsequent mood states and scores were also improved by inclusion of missing responses. The signature method provided an effective approach to the analysis of prospectively collected mood data where missing data was common and should be considered as an approach in other similar datasets.


title: High Resolution Zero-Shot Domain Adaptation of Synthetically Rendered Face Images

abstract: Generating photorealistic images of human faces at scale remains a prohibitively difficult task using computer graphics approaches. This is because these require the simulation of light to be photorealistic, which in turn requires physically accurate modelling of geometry, materials, and light sources, for both the head and the surrounding scene. Non-photorealistic renders however are increasingly easy to produce. In contrast to computer graphics approaches, generative models learned from more readily available 2D image data have been shown to produce samples of human faces that are hard to distinguish from real data. The process of learning usually corresponds to a loss of control over the shape and appearance of the generated images. For instance, even simple disentangling tasks such as modifying the hair independently of the face, which is trivial to accomplish in a computer graphics approach, remains an open research question. In this work, we propose an algorithm that matches a non-photorealistic, synthetically generated image to a latent vector of a pretrained StyleGAN2 model which, in turn, maps the vector to a photorealistic image of a person of the same pose, expression, hair, and lighting. In contrast to most previous work, we require no synthetic training data. To the best of our knowledge, this is the first algorithm of its kind to work at a resolution of 1K and represents a significant leap forward in visual realism.


title: Numerical analysis of a structure-preserving space-discretization for an anisotropic and heterogeneous boundary controlled N-dimensional wave equation as port-Hamiltonian system

abstract: The anisotropic and heterogeneous $N$-dimensional wave equation, controlled and observed at the boundary, is considered as a port-Hamiltonian system. A recent structure-preserving mixed Galerkin method is applied, leading directly to a finite-dimensional port-Hamiltonian system: its numerical analysis is carried out in a general framework. Optimal choices of mixed finite elements are then proved to reach the best trade-off between the convergence rate and the number of degrees of freedom for the state error. Exta compatibility conditions are identified for the Hamiltonian error to be twice that of the state error, and numerical evidence is provided that some combinations of finite element families meet these conditions. Numerical simulations in 2D are performed to illustrate the main theorems among several choices of classical finite element families. Several test cases are provided, including non-convex domain, anisotropic or hetergoneous cases and absorbing boundary conditions.


title: Beltrami fields exhibit knots and chaos almost surely

abstract: In this paper we show that, with probability 1, a random Beltrami field exhibits chaotic regions that coexist with invariant tori of complicated topologies. The motivation to consider this question, which arises in the study of stationary Euler flows in dimension 3, is V.I. Arnold's 1965 conjecture that a typical Beltrami field exhibits the same complexity as the restriction to an energy hypersurface of a generic Hamiltonian system with two degrees of freedom. The proof hinges on the obtention of asymptotic bounds for the number of horseshoes, zeros, and knotted invariant tori and periodic trajectories that a Gaussian random Beltrami field exhibits, which we obtain through a nontrivial extension of the Nazarov--Sodin theory for Gaussian random monochromatic waves and the application of different tools from the theory of dynamical systems, including KAM theory, Melnikov analysis and hyperbolicity. Our results hold both in the case of Beltrami fields on $\mathbf{R}^3$ and of high-frequency Beltrami fields on the 3-torus.


title: Capillary rise in vuggy media

abstract: Carbonates are highly heterogeneous formations with large variations in pore size distribution and pore space topology, which results in complex multiphase flow behavior. Here we investigate the spontaneous imbibition behavior of fluids in vuggy carbonates. Glass beads of 1.0 mm diameter, with dissolvable inclusions, are sintered to form multiple configurations of heterogeneous vuggy core with variations in matrix porosity, vug size, vug spatial location, and number of vugs. The core fabrication process is repeatable and allows the impact of vug textural properties to be investigated in a controlled manner. Capillary rise experiments are conducted in these proxy vuggy carbonate core and compared with the homogeneous non-vuggy core as reference. Continuous optical imaging is performed to track the position of the air-water interface in the cores. To understand the change in capillary height in the presence of a vug, a volume-of-fluid two-phase numerical simulation is performed in a parallel set of connected and disconnected tubes. Finally x-ray tomography scans are performed to identify the shape of the air-water interface in a select few cores. The results can be summarized as follows: disconnected vugs result in higher capillary rise compared to non-vuggy porous media. The vugs act as capillary barriers, diverting fluid flow to the adjacent connected channels, which ultimately results in a higher overall capillary rise. The results of this work highlight that radius of spontaneous invasion of aqueous phases, such as fracture fluid and hazardous wastes, are affected by vug porosity but not their distribution.


title: Incremental inference of collective graphical models

abstract: We consider incremental inference problems from aggregate data for collective dynamics. In particular, we address the problem of estimating the aggregate marginals of a Markov chain from noisy aggregate observations in an incremental (online) fashion. We propose a sliding window Sinkhorn belief propagation (SW-SBP) algorithm that utilizes a sliding window filter of the most recent noisy aggregate observations along with encoded information from discarded observations. Our algorithm is built upon the recently proposed multi-marginal optimal transport based SBP algorithm that leverages standard belief propagation and Sinkhorn algorithm to solve inference problems from aggregate data. We demonstrate the performance of our algorithm on applications such as inferring population flow from aggregate observations.

categories: stat.ML cs.IT cs.LG cs.SY eess.SY math.IT math.OC


title: Denotational recurrence extraction for amortized analysis

abstract: A typical way of analyzing the time complexity of functional programs is to extract a recurrence expressing the running time of the program in terms of the size of its input, and then to solve the recurrence to obtain a big-O bound. For recurrence extraction to be compositional, it is also necessary to extract recurrences for the size of outputs of helper functions. Previous work has developed techniques for using logical relations to state a formal correctness theorem for a general recurrence extraction translation: a program is bounded by a recurrence when the operational cost is bounded by the extracted cost, and the output value is bounded, according to a value bounding relation defined by induction on types, by the extracted size. This previous work supports higher-order functions by viewing recurrences as programs in a lambda-calculus, or as mathematical entities in a denotational semantics thereof. In this paper, we extend these techniques to support amortized analysis, where costs are rearranged from one portion of a program to another to achieve more precise bounds. We give an intermediate language in which programs can be annotated according to the banker's method of amortized analysis; this language has an affine type system to ensure credits are not spent more than once. We give a recurrence extraction translation of this language into a recurrence language, a simply-typed lambda-calculus with a cost type, and state and prove a bounding logical relation expressing the correctness of this translation. The recurrence language has a denotational semantics in preorders, and we use this semantics to solve recurrences, e.g analyzing binary counters and splay trees.


title: SAR2SAR: a semi-supervised despeckling algorithm for SAR images

abstract: Speckle reduction is a key step in many remote sensing applications. By strongly affecting synthetic aperture radar (SAR) images, it makes them difficult to analyse. Due to the difficulty to model the spatial correlation of speckle, a deep learning algorithm with self-supervision is proposed in this paper: SAR2SAR. Multi-temporal time series are leveraged and the neural network learns to restore SAR images by only looking at noisy acquisitions. To this purpose, the recently proposed noise2noise framework has been employed. The strategy to adapt it to SAR despeckling is presented, based on a compensation of temporal changes and a loss function adapted to the statistics of speckle. A study with synthetic speckle noise is presented to compare the performances of the proposed method with other state-of-the-art filters. Then, results on real images are discussed, to show the potential of the proposed algorithm. The code is made available to allow testing and reproducible research in this field.


title: Statistical Isotropy of the CMB E-mode signal

abstract: We test the statistical isotropy (SI) of the $E$-mode polarization of the cosmic microwave background (CMB) radiation observed by the Planck satellite using two statistics, namely, the contour Minkowski Tensor (CMT) and the Directional statistic ($\mathcal{D}$ statistic). The parameter $\alpha$ obtained from the CMT provides information of the alignment of structures and can be used to infer statistical properties such as Gaussianity and SI of random fields. The $\mathcal{D}$ statistic is based on detecting preferred directionality shown by vectors defined by the field. These two tests are complementary to each other in terms of sensitivity at different angular scales. The CMT is sensitive towards small-scale information present in the CMB map while $\mathcal{D}$ statistic is more sensitive at large-scales. We compute $\alpha$ and $\mathcal{D}$ statistic for the observed $E$-mode of CMB polarization, focusing on the SMICA maps, and compare with the values calculated using FFP10 SMICA simulations which contain both CMB and noise. We find good agreement between the observed data and simulations. Further, in order to specifically analyze the CMB signal in the data, we compare the values of the two statistics obtained from the observed Planck data with the values obtained from isotropic simulations having the same power spectrum, and from SMICA noise simulations. We find no statistically significant deviation from SI using the $\alpha$ parameter. From $\mathcal{D}$ statistic we find that the data shows slight deviation from SI at large angular scales.


title: Hyperfine structure and electric quadrupole transitions in the deuterium molecular ion

abstract: Molecular hydrogen ions are of metrological relevance due to the possibility of precise theoretical evaluation of their spectrum and of external-field-induced shifts. In homonuclear molecular ions the electric dipole $E1$ transitions are strongly suppressed, and of primary laser spectroscopy interest is the electric quadrupole ($E2$) transition spectrum. In continuation of previous work on the H$_2^+$ ion, we report here the results of the calculations of the hyperfine structure of the laser-induced electric quadrupole transitions between a large set of ro-vibrational states of D$_2^+$; the inaccuracies of previous evaluations have been corrected. The effects of the laser polarization are studied in detail. We show that the electric quadrupole moment of the deuteron can in principle be determined with low fractional uncertainty $(\simeq1\times10^{-4})$ by comparing the results presented here with future data from precision spectroscopy of D$_2^+$.


title: Long-range ballistic transport of Brown-Zak fermions in graphene superlattices

abstract: In quantizing magnetic fields, graphene superlattices exhibit a complex fractal spectrum often referred to as the Hofstadter butterfly. It can be viewed as a collection of Landau levels that arise from quantization of Brown-Zak minibands recurring at rational ($p/q$) fractions of the magnetic flux quantum per superlattice unit cell. Here we show that, in graphene-on-boron-nitride superlattices, Brown-Zak fermions can exhibit mobilities above 10$^6$ cm$^2$V$^{-1}$s$^{-1}$ and the mean free path exceeding several micrometers. The exceptional quality of our devices allows us to show that Brown-Zak minibands are $4q$ times degenerate and all the degeneracies (spin, valley and mini-valley) can be lifted by exchange interactions below 1K. We also found negative bend resistance at $1/q$ fractions for electrical probes placed as far as several micrometers apart. The latter observation highlights the fact that Brown-Zak fermions are Bloch quasiparticles propagating in high fields along straight trajectories, just like electrons in zero field.


title: ALMA observations of CS in NGC 1068: chemistry and excitation

abstract: We present results from Atacama Large Millimeter/submillimeter Array (ALMA) observations of CS from the nearby galaxy NGC 1068 ($\sim14$ Mpc). This Seyfert 2 barred galaxy possesses a circumnuclear disc (CND, $r\sim200$ pc) and a starburst ring (SB ring, $r\sim1.3$ kpc). These high-resolution maps ($\sim0.5$", $\sim35$ pc) allow us to analyse specific sub-regions in the galaxy and investigate differences in line intensity ratios and physical conditions, particularly those between the CND and SB ring. Local thermodynamic equilibrium (LTE) analysis of the gas is used to calculate CS densities in each sub-region, followed by non-LTE analysis conducted using the radiative transfer code RADEX to fit observations and constrain gas temperature, CS column density and hydrogen density. Finally, the chemical code UCLCHEM is used to reconstruct the gas, allowing an insight into its origin and chemical history. The density of hydrogen in the CND is found to be $\geq10^5$ cm$^{-2}$, although exact values vary, reaching $10^6$ cm$^{-2}$ at the AGN. The conditions in the two arms of the SB ring appear similar to one another, though the density found ($\sim10^4$ cm$^{-2}$) is lower than in the CND. The temperature in the CND increases from east to west, and is also overall greater than found in the SB ring. These modelling methods indicate the requirement for multi-phase gas components in order to fit the observed emission over the galaxy. A larger number of high resolution transitions across the SLED may allow for further constraining of the conditions, particularly in the SB ring.


title: Global endpoint Strichartz estimates for Schr\"odinger equations on the cylinder $\mathbb{R}\times\mathbb{T}$

abstract: We prove a sharp, global-in-time Strichartz estimate for the Schr\"odinger equation on the cylinder $\mathbb{R}\times\mathbb{T}$.


title: Nearest Neighbour Based Estimates of Gradients: Sharp Nonasymptotic Bounds and Applications

abstract: Motivated by a wide variety of applications, ranging from stochastic optimization to dimension reduction through variable selection, the problem of estimating gradients accurately is of crucial importance in statistics and learning theory. We consider here the classic regression setup, where a real valued square integrable r.v. $Y$ is to be predicted upon observing a (possibly high dimensional) random vector $X$ by means of a predictive function $f(X)$ as accurately as possible in the mean-squared sense and study a nearest-neighbour-based pointwise estimate of the gradient of the optimal predictive function, the regression function $m(x)=\mathbb{E}[Y\mid X=x]$. Under classic smoothness conditions combined with the assumption that the tails of $Y-m(X)$ are sub-Gaussian, we prove nonasymptotic bounds improving upon those obtained for alternative estimation methods. Beyond the novel theoretical results established, several illustrative numerical experiments have been carried out. The latter provide strong empirical evidence that the estimation method proposed works very well for various statistical problems involving gradient estimation, namely dimensionality reduction, stochastic gradient descent optimization and quantifying disentanglement.


title: The nature of $X(3872)$ from high-multiplicity $pp$ collisions

abstract: The structure of exotic resonances that do not trivially fit the usual quark model expectations has been a matter of intense scientific debate during the last two decades. A possible way of estimating the size of these states is to study their behavior when immersed in QCD matter. Recently, LHCb has measured the relative abundance of the exotic $X(3872)$ over the ordinary $\psi(2S)$. We use the comover interaction model to study the yield of a compact $X(3872)$. To confirm the reliability of the model in high-multiplicity $pp$ collisions, we describe the suppression of excited over ground $\Upsilon$ states. With this at hand, we show that the size of the compact $X(3872)$ would be slightly larger than that of the $\psi(2S)$. If the $X(3872)$ is instead assumed to be a meson molecule of large size, we argue that its evolution in QCD matter should be described via a coalescence model, as suggested by data on deuteron production. We show that the predictions of this model for the $X(3872)$ are in contrast with data.


title: K2: Background Survey -- the search for undiscovered transients in Kepler/K2 data

abstract: The K2 mission of the Kepler Space Telescope offers a unique possibility to examine sources of both Galactic and Extra-galactic origin with high cadence photometry. Alongside the multitude of supernovae and quasars detected within targeted galaxies, it is likely that Kepler has serendipitously observed many transients throughout K2. Such events will likely have occurred in background pixels, coincidentally surrounding science targets. Analysing the background pixels presents the possibility to conduct a high cadence survey with areas of a few square degrees per campaign. We demonstrate the capacity to independently recover key K2 transients such as KSN 2015K and SN 2018oh. With this survey, we expect to detect numerous transients and determine the first comprehensive rates for transients with lifetimes $\leq1$ day.


title: Uniqueness in determining the orders of time and spatial fractional derivatives

abstract: We prove the uniqueness in determining both orders of fractional time derivatives and spatial derivatives in diffusion equations by pointwise data. The proof relies on the eigenfunction expansion and the asymptotics of the Mittag-Leffler function.


title: VAMANA: Modeling Binary Black Hole Population with Minimal Assumptions

abstract: The population analysis of compact binaries involves the reconstruction of some of the gravitational wave (GW) signal parameters, such as, the mass and the spin distribution, that gave rise to the observed data. This article introduces VAMANA, which reconstructs the binary black hole population using a mixture model and facilitates excellent density measurement as informed by the data. VAMANA uses a mixture of weighted Gaussians to reconstruct the chirp mass distribution. We expect Gaussian mixtures to provide flexibility in modeling complex distributions and enable us in capturing details in the astrophysical chirp mass distribution. Each of the Gaussian in the mixture is combined with another Gaussian and a power-law to simultaneously model the spin component aligned with the orbital angular momentum and the mass ratio distribution, thus also allowing us to capture their variation with the chirp mass. Additionally, we can also introduce broadband smoothing by restricting the Gaussian mixture to lie within a threshold distance of a predefined reference chirp mass distribution. Using simulated data we show the robustness of our method in reconstructing complex populations for a large number of observations. We also apply our method to the publicly available catalog of GW observations made during LIGO's and Virgo's first and second observation runs and present the reconstructed mass, spin distribution, and the estimated merger rate of binary black holes.


title: Arbitrary positive powers of semicirculant and r-circulant matrices

abstract: We provide a novel recursive method, which does not require any assumption, to compute the entries of the kth power of a semicirculant matrix. As an application, a method for computing the entries of the kth power of r-circulant matrices is also presented.


title: Chemical abundances in Sgr A East: evidence for a Type Iax supernova remnant

abstract: Recent observations have shown a remarkable diversity of observational behaviors and explosion mechanisms in thermonuclear supernovae (SNe). An emerging class of peculiar thermonuclear SNe, called Type Iax, show photometric and spectroscopic behaviors distinct from normal Type Ia. Their origin remains highly controversial, but pure turbulent deflagration of white dwarfs (WDs) has been regarded as the leading formation theory. The large population of Type Iax indicates the existence of unidentified Galactic Type Iax supernova remnants (SNRs). We report evidence that SNR Sgr A East in the Galactic center resulted from a pure turbulent deflagration of a Chandrasekhar-mass carbon-oxygen WD, an explosion mechanism used for Type Iax SNe. Our X-ray spectroscopic study of Sgr A East using 3 Ms of Chandra data shows a low ratio of intermediate-mass elements to Fe and large Mn/Fe and Ni/Fe ratios. This abundance pattern does not accord with the core-collapse or normal Type Ia models. Sgr A East is thus the first Galactic SNR for which a likely Type Iax origin has been proposed and the nearest target for studying this peculiar class. We compared Sgr A East with the Fe-rich SNRs 3C 397 and W49B, which also have high Mn and Cr abundances and were claimed to result from deflagration-to-detonation explosions of Chandrasekhar-mass WDs (although with disputes). Our study shows that they have distinct abundance patterns. The X-ray spectroscopic studies of thermonuclear SNRs provide observational evidence for the theories that there are diverse explosion channels and various metal outputs for Chandrasekhar-mass WDs.


title: Optimal non-classical correlations of light with a levitated nano-sphere

abstract: Nonclassical correlations provide a resource for many applications in quantum technology as well as providing strong evidence that a system is indeed operating in the quantum regime. Optomechanical systems can be arranged to generate quantum entanglement between the mechanics and a mode of travelling light. Here we propose automated optimisation of the production of quantum correlations in such a system, beyond what can be achieved through analytical methods, by applying Bayesian optimisation to the control parameters. Two-mode optomechanical squeezing experiment is simulated using a detailed theoretical model of the system, while the Bayesian optimisation process modifies the controllable parameters in order to maximise the non-classical two-mode squeezing and its detection, independently of the inner workings of the model. The Bayesian optimisation treats the simulations or the experiments as a black box. This we refer to as \emph{theory-blind} optimisation, and the optimisation process is designed to be unaware of whether it is working with a simulation or the actual experimental setup. We find that in the experimentally relevant thermal regimes, the ability to vary and optimise a broad array of control parameters provides access to large values of two-mode squeezing that would otherwise be difficult or intractable to discover. In particular we observe that modulation of the driving frequency around the resonant sideband, when added to the set of control parameters, produces strong nonclassical correlations greater on average than the maximum achieved by optimising over the remaining parameters. We also find that using our optimisation approach raises the upper limit to the thermal regime in which squeezing can be achieved. This extends the range of experimental setups in which non-classical correlations could be generated beyond the region of high quantum cooperativity.


title: Arc-Like Variable Bunch Compressors

abstract: Electron bunch compressors formed of achromat arcs have a natural advantage over the more commonly used chicane compressors in that linearisation of the longitudinal phase space is of the correct sign to compensate for the curvature imprinted by rf acceleration. Here we extend the utility of arc compressors to enable variation of the longitudinal compaction within a fixed footprint. We also show that this variability can be achieved independently order-by-order in momentum deviation. The technique we employ consists of additional dipoles, leading to the advantageous property that variability can be achieved without incurring significant penalty in terms of chromatic degradation. We show this by comparison to an alternative system where additional quadrupoles are utilised to enable variation of momentum compaction. Each of these alternative approaches are being considered in the context of an upgrade of the MAX IV linac, Sweden, to enable a soft X-ray free-electron laser (FEL) in addition to its existing functions.


title: Neutrino interaction classification with a convolutional neural network in the DUNE far detector

abstract: The Deep Underground Neutrino Experiment is a next-generation neutrino oscillation experiment that aims to measure $CP$-violation in the neutrino sector as part of a wider physics program. A deep learning approach based on a convolutional neural network has been developed to provide highly efficient and pure selections of electron neutrino and muon neutrino charged-current interactions. The electron neutrino (antineutrino) selection efficiency peaks at 90% (94%) and exceeds 85% (90%) for reconstructed neutrino energies between 2-5 GeV. The muon neutrino (antineutrino) event selection is found to have a maximum efficiency of 96% (97%) and exceeds 90% (95%) efficiency for reconstructed neutrino energies above 2 GeV. When considering all electron neutrino and antineutrino interactions as signal, a selection purity of 90% is achieved. These event selections are critical to maximize the sensitivity of the experiment to $CP$-violating effects.


title: Case of two electrostatics problems: Can providing a diagram adversely impact introductory physics students' problem solving performance?

abstract: Drawing appropriate diagrams is a useful problem solving heuristic that can transform a problem into a representation that is easier to exploit for solving it. One major focus while helping introductory physics students learn effective problem solving is to help them understand that drawing diagrams can facilitate problem solution. We conducted an investigation in which two different interventions were implemented during recitation quizzes in a large enrollment algebra-based introductory physics course. Students were either (i) asked to solve problems in which the diagrams were drawn for them or (ii) explicitly told to draw a diagram. A comparison group was not given any instruction regarding diagrams. We developed rubrics to score the problem solving performance of students in different intervention groups and investigated ten problems. We found that students who were provided diagrams never performed better and actually performed worse than the other students on three problems, one involving standing sound waves in a tube (discussed elsewhere) and two problems in electricity which we focus on here. These two problems were the only problems in electricity that involved considerations of initial and final conditions, which may partly account for why students provided with diagrams performed significantly worse than students who were not provided with diagrams. In order to explore potential reasons for this finding, we conducted interviews with students and found that some students provided with diagrams may have spent less time on the conceptual analysis and planning stage of the problem solving process. In particular, those provided with the diagram were more likely to jump into the implementation stage of problem solving early without fully analyzing and understanding the problem, which can increase the likelihood of mistakes in solutions.


title: Option Pricing Under a Discrete-Time Markov Switching Stochastic Volatility with Co-Jump Model

abstract: We consider option pricing using a discrete-time Markov switching stochastic volatility with co-jump model, which can model volatility clustering and varying mean-reversion speeds of volatility. For pricing European options, we develop a computationally efficient method for obtaining the probability distribution of average integrated variance (AIV), which is key to option pricing under stochastic-volatility-type models. Building upon the efficiency of the European option pricing approach, we are able to price an American-style option, by converting its pricing into the pricing of a portfolio of European options. Our work also provides constructive guidance for analyzing derivatives based on variance, e.g., the variance swap. Numerical results indicate our methods can be implemented very efficiently and accurately.


title: Object-Centric Learning with Slot Attention

abstract: Learning object-centric representations of complex scenes is a promising step towards enabling efficient abstract reasoning from low-level perceptual features. Yet, most deep learning approaches learn distributed representations that do not capture the compositional properties of natural scenes. In this paper, we present the Slot Attention module, an architectural component that interfaces with perceptual representations such as the output of a convolutional neural network and produces a set of task-dependent abstract representations which we call slots. These slots are exchangeable and can bind to any object in the input by specializing through a competitive procedure over multiple rounds of attention. We empirically demonstrate that Slot Attention can extract object-centric representations that enable generalization to unseen compositions when trained on unsupervised object discovery and supervised property prediction tasks.


title: Cross-Supervised Object Detection

abstract: After learning a new object category from image-level annotations (with no object bounding boxes), humans are remarkably good at precisely localizing those objects. However, building good object localizers (i.e., detectors) currently requires expensive instance-level annotations. While some work has been done on learning detectors from weakly labeled samples (with only class labels), these detectors do poorly at localization. In this work, we show how to build better object detectors from weakly labeled images of new categories by leveraging knowledge learned from fully labeled base categories. We call this novel learning paradigm cross-supervised object detection. We propose a unified framework that combines a detection head trained from instance-level annotations and a recognition head learned from image-level annotations, together with a spatial correlation module that bridges the gap between detection and recognition. These contributions enable us to better detect novel objects with image-level annotations in complex multi-object scenes such as the COCO dataset.


title: A Loss Function for Generative Neural Networks Based on Watson's Perceptual Model

abstract: To train Variational Autoencoders (VAEs) to generate realistic imagery requires a loss function that reflects human perception of image similarity. We propose such a loss function based on Watson's perceptual model, which computes a weighted distance in frequency space and accounts for luminance and contrast masking. We extend the model to color images, increase its robustness to translation by using the Fourier Transform, remove artifacts due to splitting the image into blocks, and make it differentiable. In experiments, VAEs trained with the new loss function generated realistic, high-quality image samples. Compared to using the Euclidean distance and the Structural Similarity Index, the images were less blurry; compared to deep neural network based losses, the new approach required less computational resources and generated images with less artifacts.


title: Exploring one aspect of pedagogical content knowledge of teaching assistants using the Conceptual Survey of Electricity and Magnetism

abstract: Instruction is unlikely to be effective if instructors do not know the common alternate conceptions of introductory physics students and explicitly take into account common student difficulties in their instructional design. Here, we discuss research involving the Conceptual Survey of Electricity and Magnetism (CSEM) to evaluate one aspect of the pedagogical content knowledge of teaching assistants (TAs): knowledge of introductory students' alternate conceptions in electricity and magnetism as revealed by the CSEM. For each item on the CSEM, the TAs were asked to identify the most common incorrect answer choice selected by introductory physics students if they did not know the correct answer after traditional instruction. Then, we used introductory student CSEM post-test data to assess the extent to which TAs were able to identify the most common alternate conception of introductory students in each question on the CSEM. We find that the TAs were thoughtful when attempting to identify common student difficulties and they enjoyed learning about student difficulties this way. However, they struggled to identify many common difficulties of introductory students that persist after traditional instruction. We discuss specific alternate conceptions that persist after traditional instruction, the extent to which TAs were able to identify them, and results from think-aloud interviews with TAs which provided valuable information regarding why the TAs sometimes selected certain alternate conceptions as the most common but were instead very rare among introductory students. We also discuss how tasks such as the one used in this study can be used in professional development programs to engender productive discussions about the importance of being knowledgeable about student alternate conceptions in order to help students learn.


title: Computing Light Transport Gradients using the Adjoint Method

abstract: This paper proposes a new equation from continuous adjoint theory to compute the gradient of quantities governed by the Transport Theory of light. Unlike discrete gradients ala autograd, which work at the code level, we first formulate the continuous theory and then discretize it. The key insight of this paper is that computing gradients in Transport Theory is akin to computing the importance, a quantity adjoint to radiance that satisfies an adjoint equation. Importance tells us where to look for light that matters. This is one of the key insights of this paper. In fact, this mathematical journey started from a whimsical thought that these adjoints might be related. Computing gradients is therefore no more complicated than computing the importance field. This insight and the following paper hopefully will shed some light on this complicated problem and ease the implementations of gradient computations in existing path tracers.


title: Conformal cyclic evolution of phantom energy dominated universe

abstract: From the Wheeler-Dewitt solutions, the scale factor of the initial universe is discussed. In this study, scale factors from Wheeler-Dewitt solutions, loop quantum gravity, and phantom energy dominated stages are compared. Certain modifications have been attempted in scale factor and quantum potentials driven by canonical quantum gravity approaches. Their results are discussed in this work. Despite an increment of phantom energy density, avoidance of Big Rip is reported. Scale factors predicted from various models are discussed in this work. The relationship between scale factors and the smooth continuation of Aeon is discussed by the application of conformal cyclic cosmology. Quantum potentials for various models are correlated and a correction parameter is included in the cosmological constant. Phantom energy dominated, final stage non-singular evolution of the universe is found. Eternal increment of phantom energy density without interacting with dark matter is reported for the consequence of the evolution of the future universe. Also, the non-interacting solutions of phantom energy and dark matter are explained. As the evolution continues even after the final singularity is approached, the validity of conformal cyclic cosmology is predicted. Non zero values for the scale factor for the set of eigenvalues are presented. Results are compared with supersymmetric classical cosmology. The non-interacting solutions are compared with SiBI solutions


title: Intrinsic Reward Driven Imitation Learning via Generative Model

abstract: Imitation learning in a high-dimensional environment is challenging. Most inverse reinforcement learning (IRL) methods fail to outperform the demonstrator in such a high-dimensional environment, e.g., Atari domain. To address this challenge, we propose a novel reward learning module to generate intrinsic reward signals via a generative model. Our generative method can perform better forward state transition and backward action encoding, which improves the module's dynamics modeling ability in the environment. Thus, our module provides the imitation agent both the intrinsic intention of the demonstrator and a better exploration ability, which is critical for the agent to outperform the demonstrator. Empirical results show that our method outperforms state-of-the-art IRL methods on multiple Atari games, even with one-life demonstration. Remarkably, our method achieves performance that is up to 5 times the performance of the demonstration.


title: YBCO-based non-volatile ReRAM tested in Low Earth Orbit

abstract: An YBCO-based test structure corresponding to the family of ReRAM devices associated with the valence change mechanism is presented. We have characterized its electrical response previous to its lift-off to a Low Earth Orbit (LEO) using standard electronics and also with the dedicated LabOSat-01 controller. Similar results were obtained in both cases. After about 200 days at LEO on board a small satellite, electrical tests started on the memory device using the LabOSat-01 controller. We discuss the results of the first 150 tests, performed along a 433-day time interval in space. The memory device remained operational despite the hostile conditions that involved launching, lift-off vibrations, permanent thermal cycling and exposure to ionizing radiation, with doses 3 orders of magnitude greater than the usual ones on Earth. The device showed resistive switching and IV characteristics similar to those measured on Earth, although with changes that follow a smooth drift in time. A detailed study of the electrical transport mechanisms, based on previous models that indicate the existence of various conducting mechanisms through the metal-YBCO interface showed that the observed drift can be associated with a local temperature drift at the LabOSat controller, with no clear evidence that allows determining changes in the underlying microscopic factors. These results show the reliability of complex-oxide non-volatile ReRAM-based devices in order to operate under all the hostile conditions encountered in space-borne applications.


title: On homomorphisms involving a hook Weyl module

abstract: Let $k$ be an infinite field of positive characteristic. We determine all homomorphisms between Weyl modules for $GLn(k)$, where one of the partitions is a hook. As a consequence we obtain a nonvanishing result concerning homomorphisms between Weyl modules for algebraic groups of type B, C and D when one of the partitions is a hook.


title: Extraction of Isotropic Electron-Nuclear Hyperfine Coupling Constants of Paramagnetic Point Defects from Near-Zero Field Magnetoresistance Spectra via Least Squares Fitting to Models Developed from the Stochastic Quantum Liouville Equation

abstract: We report on a method by which we can systematically extract spectroscopic information such as isotropic electron-nuclear hyperfine coupling constants from near-zero field magnetoresistance spectra. The method utilizes a least squares fitting of models developed from the stochastic quantum Liouville equation. We applied our fitting algorithm to two distinct material systems: Si/SiO2 MOSFETs, and a-Si:H MIS capacitors. Our fitted results and hyperfine parameters are in reasonable agreement with existing knowledge of the defects present in the systems. Our work indicates that the NZFMR response and fitting of the NZFMR spectrum via models developed from the stochastic quantum Liouville equation could be a relatively simple yet powerful addition to the family of spin-based techniques used to explore the chemical and structural nature of point defects in semiconductor devices and insulators.


title: Gas kinematics of key prebiotic molecules in GV Tau N revealed with an ALMA, PdBI, and Herschel synergy

abstract: A large effort has been made to detect warm gas in the planet formation zone of circumstellar discs using space and ground-based near infrared facilities. GV Tau N, the most obscured component of the GV Tau system, is an outstanding source, being one of the first targets detected in HCN and the only one detected in CH$_4$ so far. Although near infrared observations have shed light on its chemical content, the physical structure and kinematics of the circumstellar matter remained unknown. We use interferometric images of the HCN 3-2 and $^{13}$CO 3-2 lines, and far-IR observations of $^{13}$CO, HCN, CN and H$_2$O transitions to discern the morphology, kinematics, and chemistry of the dense gas close to the star. These observations constitute the first detection of H$_2$O towards GV Tau N. Moreover, ALMA high spatial resolution (~ 7 au) images of the continuum at 1.1 mm and the HCN 3-2 line resolve different gas components towards GV Tau N, a gaseous disc with R~25 au, an ionized jet, and one (or two) molecular outflows. The asymmetric morphology of the gaseous disc shows that it has been eroded by the jet. All observations can be explained if GV Tau N is binary, and the primary component has a highly inclined individual disc relative to the circumbinary disc. We discuss the origin of the water and the other molecules emission according to this scenario. In particular, we propose that the water emission would come from the disrupted gaseous disc and the molecular outflows.


title: Defective nematogenesis: Gauss curvature in programmable shape-responsive sheets with topological defects

abstract: Flat sheets encoded with patterns of contraction/elongation morph into curved surfaces. If the surfaces bear Gauss curvature, the resulting actuation can be strong and powerful. We deploy the Gauss-Bonnet theorem to deduce the Gauss curvature encoded in a pattern of uniform-magnitude contraction/elongation with spatially varying direction, as is commonly implemented in patterned liquid crystal elastomers. This approach reveals two fundamentally distinct contributions: a structural curvature which depends on the precise form of the pattern, and a topological curvature generated by defects in the contractile direction. These curvatures grow as different functions the contraction/elongation magnitude, explaining the apparent contradiction between previous calculations for simple +1 defects, and smooth defect-free patterns. We verify these structural and topological contributions by conducting numerical shell calculations on sheets encoded with simple higher-order contractile defects to reveal their activated morphology. Finally we calculate the Gauss curvature generated by patterns with spatially varying magnitude and direction, which leads to additional magnitude gradient contributions to the structural term. We anticipate this form will be useful whenever magnitude and direction are natural variables, including in describing the contraction of a muscle along its patterned fiber direction, or a tissue growing by elongating its cells.


title: Constrained-optimization Approach Delivers Superior Classical Performance for Graph Partitioning via Quantum-ready Method

abstract: Graph partitioning is one of an important set of well-known compute-intense (NP-hard) graph problems that devolve to discrete constrained optimization. We sampled solutions to the problem via two different quantum-ready methods to understand the pros and cons of each method. First we formulated and sampled the problem as a quadratic unconstrained binary optimization (QUBO) problem, via the best known QUBO formulation, using a best-in-class QUBO sampler running purely classically. Second, we formulated the problem at a higher level, as a set of constraints and an objective function, and sampled it with a constrained-optimization sampler (which internally samples the problem via QUBOs also sampled classically). We find that both approaches often deliver better partitions than the purpose-built classical graph partitioners. Further, we find that the constrained-optimization approach is often able to deliver better partitions in less time than the bespoke-QUBO approach, without knowledge of the graph-partitioning problem. Stepping back from graph partitioning itself, one key controversial question is whether bespoke algorithms or general tools are more likely to deliver the power of QCs to real-world users. These results bear on that question, though they require confirmation on other problems and instances as well as replacement of the low-level sampler by a QC. Still, this early evidence supports the proposition that general tools will contribute significantly to a range of problems, expanding the impact of QCs. This benefit is independent of the low-level sampler employed, whether software or QC, so reinforces the need for more work on high-level optimization. An early version of such software is commercially available in the cloud today, delivering superior classical performance for some problems, enables quantum-forward organizations to migrate to quantum-ready methods now.


title: Parsummable categories as a strictification of symmetric monoidal categories

abstract: We prove that the homotopy theory of parsummable categories (as defined by Schwede) with respect to the underlying equivalences of categories is equivalent to the usual homotopy theory of symmetric monoidal categories. In particular, this yields a model of symmetric monoidal categories in terms of categories equipped with a strictly commutative, associative, and unital (but only partially defined) operation.


title: Machine learning-based clinical prediction modeling -- A practical guide for clinicians

abstract: In the emerging era of big data, larger available clinical datasets and computational advances have sparked a massive interest in machine learning-based approaches. The number of manuscripts related to machine learning or artificial intelligence has exponentially increased over the past years. As analytical machine learning tools become readily available for clinicians to use, the understanding of key concepts and the awareness of analytical pitfalls are increasingly required for clinicians, investigators, reviewers and editors, who even as experts in their clinical field, sometimes find themselves insufficiently equipped to evaluate machine learning methodologies. In the first section, we provide explanations on the general principles of machine learning, as well as analytical steps required for successful machine learning-based predictive modelling - which is the focus of this series. In further sections, we review the importance of resampling, overfitting and model generalizability as well as feature reduction and selection (Part II), strategies for model evaluation, reporting and discussion of common caveats and other points of significance (Part III), as well as offer a practical guide to classification (Part IV) and regression modelling (Part V), with a complete coding pipeline. Methodological rigor and clarity as well as understanding of the underlying reasoning of the internal workings of a machine learning approach are required, otherwise predictive applications despite being strong analytical tools are not well accepted into the clinical routine. Going forward, machine learning and artificial intelligence shape and influence modern medicine across disciplines including the field of neurosurgery.