title: The dependence of evanescent wave polarization on the losses of guided optical modes

abstract: Spin-momentum locking of evanescent waves describes the relationship between the propagation constant of an evanescent mode and the polarization of its electromagnetic field, giving rise to applications in light nano-routing and polarimetry among many others. The use of complex numbers in physics is a powerful representation in areas such as quantum mechanics or electromagnetism; it is well known that a lossy waveguide can be modelled with the addition of an imaginary part to the propagation constant. Here we explore how these losses are entangled with the polarization of the associated evanescent tails for the waveguide, revealing a well-defined mapping between waveguide losses and the Poincar\'e sphere of polarizations, in what could be understood as a "polarization-loss locking" of evanescent waves. We analyse the implications for near-field directional coupling of sources to waveguides, as optimized dipoles must take into account the losses for a perfectly unidirectional excitation. We also reveal the potential advantage of calculating the angular spectrum of a source defined in a complex, rather than the traditionally purely real, transverse wavevector space formalism.


title: Low precision logarithmic number systems: Beyond base-2

abstract: Logarithmic number systems (LNS) are used to represent real numbers in many applications using a constant base raised to a fixed-point exponent making its distribution exponential. This greatly simplifies hardware multiply, divide and square root. LNS with base-2 is most common, but in this paper we show that for low-precision LNS the choice of base has a significant impact. We make four main contributions. First, LNS is not closed under addition and subtraction, so the result is approximate. We show that choosing a suitable base can manipulate the distribution to reduce the average error. Second, we show that low-precision LNS addition and subtraction can be implemented efficiently in logic rather than commonly used ROM lookup tables, the complexity of which can be reduced by an appropriate choice of base. A similar effect is shown where the result of arithmetic has greater precision than the input. Third, where input data from external sources is not expected to be in LNS, we can reduce the conversion error by selecting a LNS base to match the expected distribution of the input. Thus, there is no one base which gives the global optimum, and base selection is a trade-off between different factors. Fourth, we show that circuits realized in LNS require lower area and power consumption for short word lengths.


title: Two-flux tunable Aharonov-Bohm caging in a photonic lattice

abstract: We study the Aharonov-Bohm caging effect in a one-dimensional lattice of theta-shaped units defining a chain of interconnected plaquettes, each one threaded by two synthetic flux lines. In the proposed system, light trapping results from the destructive interference of waves propagating along three arms, this implies that the caging effect is tunable and it can be controlled by changing the tunnel couplings $J$. These features reflect on the diffraction pattern allowing to establish a clear connection between the lattice topology and the resulting AB interference.


title: On the stickiness of CO$_{2}$ and H$_{2}$O ice particles

abstract: Laboratory experiments revealed that CO$_{2}$ ice particles stick less efficiently than H$_{2}$O ice particles, and there is an order of magnitude difference in the threshold velocity for sticking. However, the surface energies and elastic moduli of CO$_{2}$ and H$_{2}$O ices are comparable, and the reason why CO$_{2}$ ice particles were poorly sticky compared to H$_{2}$O ice particles was unclear. Here we investigate the effects of viscoelastic dissipation on the threshold velocity for sticking of ice particles using the viscoelastic contact model derived by Krijt et al. We find that the threshold velocity for sticking of CO$_{2}$ ice particles reported in experimental studies is comparable to that predicted for perfectly elastic spheres. In contrast, the threshold velocity for sticking of H$_{2}$O ice particles is an order of magnitude higher than that predicted for perfectly elastic spheres. Therefore, we conclude that the large difference in stickiness between CO$_{2}$ and H$_{2}$O ice particles would mainly originate from the difference in the strength of viscoelastic dissipation, which is controlled by the viscoelastic relaxation time.


title: DeepGLEAM: A hybrid mechanistic and deep learning model for COVID-19 forecasting

abstract: We introduce DeepGLEAM, a hybrid model for COVID-19 forecasting. DeepGLEAM combines a mechanistic stochastic simulation model GLEAM with deep learning. It uses deep learning to learn the correction terms from GLEAM, which leads to improved performance. We further integrate various uncertainty quantification methods to generate confidence intervals. We demonstrate DeepGLEAM on real-world COVID-19 mortality forecasting tasks.


title: Learning Depth via Leveraging Semantics: Self-supervised Monocular Depth Estimation with Both Implicit and Explicit Semantic Guidance

abstract: Self-supervised depth estimation has made a great success in learning depth from unlabeled image sequences. While the mappings between image and pixel-wise depth are well-studied in current methods, the correlation between image, depth and scene semantics, however, is less considered. This hinders the network to better understand the real geometry of the scene, since the contextual clues, contribute not only the latent representations of scene depth, but also the straight constraints for depth map. In this paper, we leverage the two benefits by proposing the implicit and explicit semantic guidance for accurate self-supervised depth estimation. We propose a Semantic-aware Spatial Feature Alignment (SSFA) scheme to effectively align implicit semantic features with depth features for scene-aware depth estimation. We also propose a semantic-guided ranking loss to explicitly constrain the estimated depth maps to be consistent with real scene contextual properties. Both semantic label noise and prediction uncertainty is considered to yield reliable depth supervisions. Extensive experimental results show that our method produces high quality depth maps which are consistently superior either on complex scenes or diverse semantic categories, and outperforms the state-of-the-art methods by a significant margin.


title: Theoretical Study of Isotope Production in The Peripheral Heavy-ion Collision 136Xe + Pb at 1 GeV/nucleon

abstract: We have studied the fragment yields emitted from the fragmentation of excited projectile nuclei in peripheral collisions of 136Xe+Pb at 1 GeV/nucleon, and measured with the high-resolution magnetic spectrometer, the Fragment Separator (FRS) of GSI. The mass, charge and isotope distributions of nuclear fragments formed in the reactions were calculated within a statistical ensemble approach and compared to the experimental data. The ensemble of excited projectilelike source nuclei were created in the framework of a previous analysis of similar reactions performed at 600 MeV/nucleon (ALADIN-experiments, GSI). In addition, the ensemble of the low-excited compound nuclei is involved in the analysis. The overall agreement between theory and experiment was very satisfactory in reproducing the experimental data of isotope yields measured in the heavy-ion collisions. It is confirmed that a broad range of elements can be reproduced within a universal statistical disintegration of the excited projectile residues.


title: Destination similarity based on implicit user interest

abstract: With the digitization of travel industry, it is more and more important to understand users from their online behaviors. However, online travel industry data are more challenging to analyze due to extra sparseness, dispersed user history actions, fast change of user interest and lack of direct or indirect feedbacks. In this work, a new similarity method is proposed to measure the destination similarity in terms of implicit user interest. By comparing the proposed method to several other widely used similarity measures in recommender systems, the proposed method achieves a significant improvement on travel data. Key words: Destination similarity, Travel industry, Recommender System, Implicit user interest


title: The chromatic number of a two families of generalized Kneser graphs related to finite generalized quadrangles and finite projective 3-spaces

abstract: Let $\Gamma$ be the graph whose vertices are the chambers of the finite projective space $PG(3,q)$ with two vertices being adjacent when the corresponding chambers are in general position. It is known that the independence number of this graph is $(q^2+q+1)(q+1)^2$. For $q\ge 43$ we determine the largest independent set of $\Gamma$ and show that every maximal independent set that is not a largest one has at most constant times $q^3$ elements. For $q\ge 47$, this information is then used to show that $\Gamma$ has chromatic number $q^2+q$. Furthermore, for many families of generalized quadrangles we prove similar results for the graph that is built in the same way on the chambers of the generalized quadrangle.


title: Can single photon excitation of two spatially separated modes lead to a violation of Bell inequality via homodyne measurements?

abstract: We reconsider the all-optical homodyne-measurement based experimental schemes that aim to reveal Bell nonclassicality of a single photon, often termed `nonlocality'. We focus on the schemes put forward by Tan, Walls and Collett (TWC, 1991) and Hardy (1994). In the light of our previous work the Tan, Walls and Collett setup can be described by a precise local hidden variable model, hence the claimed nonclassicality of this proposal is apparent, whereas the nonclassicality proof proposed by Hardy is impeccable. In this work we resolve the following problem: which feature of the Hardy's approach is crucial for its successful confirmation of nonclassicality. The scheme of Hardy differs from the Tan, Walls and Collett setup in two aspects. (i) It introduces a superposition of a single photon excitation with vacuum as the initial state of one of the input modes of a 50-50 beamsplitter, which creates the superposition state of two separable (exit) modes under investigation. (ii) In the final measurements Hardy's proposal utilises a varying strengths of the local oscillator fields, whereas in the TWC case they are constant. In fact the local oscillators in Hardy's scheme are either on or off (the local setting is specified by the presence or absence of the local auxiliary field). We show that it is the varying strength of the local oscillators, from setting to setting, which is the crucial feature enabling violation of local realism in the Hardy setup, whereas it is not necessary to use initial superposition of a single photon excitation with vacuum as the initial state of the input mode. Neither one needs to operate in the fully on/off detection scheme. Despite the failure of the Tan, Walls and Collett scheme in proving Bell nonclassicality, we show that their scheme can serve as an entanglement indicator.


title: Rethinking Eye-blink: Assessing Task Difficulty through Physiological Representation of Spontaneous Blinking

abstract: Continuous assessment of task difficulty and mental workload is essential in improving the usability and accessibility of interactive systems. Eye tracking data has often been investigated to achieve this ability, with reports on the limited role of standard blink metrics. Here, we propose a new approach to the analysis of eye-blink responses for automated estimation of task difficulty. The core module is a time-frequency representation of eye-blink, which aims to capture the richness of information reflected on blinking. In our first study, we show that this method significantly improves the sensitivity to task difficulty. We then demonstrate how to form a framework where the represented patterns are analyzed with multi-dimensional Long Short-Term Memory recurrent neural networks for their non-linear mapping onto difficulty-related parameters. This framework outperformed other methods that used hand-engineered features. This approach works with any built-in camera, without requiring specialized devices. We conclude by discussing how Rethinking Eye-blink can benefit real-world applications.


title: Pocket resonances in low-energy antineutrons reactions with nuclei

abstract: Upon investigating whether the variation of the antineutron-nucleus annihilation cross-sections at very low energies satisfy Bethe-Landau's power law of $\sigma_{\rm ann} (p) \propto 1/p^{\alpha}$ behavior as a function of the antineutron momentum $p$, we uncover unexpected regular oscillatory structures in the low antineutron energy region from 0.001 to 10 MeV, with small amplitudes and narrow periodicity in the logarithm of the antineutron energies, for large-$A$ nuclei such as Pb and Ag. Subsequent semiclassical analyses of the $S$ matrices reveal that these oscillations are pocket resonances that arise from quasi-bound states inside the pocket and the interference between the waves reflecting inside the optical potential pockets with those from beyond the potential barriers, implicit in the nuclear Ramsauer effect. They are the continuation of bound states in the continuum. Experimental observations of these pocket resonances will provide vital information on the properties of the optical model potentials and the nature of the antineutron annihilation process.


title: WKB Electron Wave Functions in a Tightly Focused Laser Beam

abstract: Available laser technology is opening the possibility of testing QED experimentally in the so-called strong-field regime. This calls for developing theoretical tools to investigate strong-field QED processes in electromagnetic fields of complex spacetime structure. Here, we propose a scheme to compute electron wave functions in tightly focused laser beams by taking into account exactly the complex spacetime structure of the fields. The scheme is solely based on the validity of the Wentzel-Kramers-Brillouin (WKB) approximation and the resulting wave functions, unlike previously proposed ones [Phys. Rev. Lett. \textbf{113}, 040402 (2014)], do not rely on approximations on the classical electron trajectory. Moreover, a consistent procedure is indicated to take into account higher-order quantum effects within the WKB approach depending on higher-and-higher powers of the Planck constant. In the case of a plane-wave background field the found wave functions exactly reduce to the Volkov states, which are then written in a new and fully quasiclassical form. Finally, by using the leading-order WKB wave functions to compute the probabilities of nonlinear Compton scattering and nonlinear Breit-Wheeler pair production, it is explicitly shown that, if additionally the energies of the charges are sufficiently large that the latter are not significantly deflected by the field, the corresponding Baier's formulas are exactly reproduced for an otherwise arbitrary classical electron/positron trajectory.


title: The Golden Age of the Mathematical Finance

abstract: This paper is devoted to show that the last quarter of the past century can be considered as the golden age of the Mathematical Finance. In this period the collaboration of great economists and the best generation of probabilists, most of them from the Strasbourg's School led by Paul Andr\'e Meyer, gave rise to the foundations of this discipline. They established the two fundamentals theorems of arbitrage theory, close formulas for options, the main modelling a


title: Heat current and entropy production rate in local non-Markovian quantum dynamics of global Markovian evolution

abstract: We examine the elements of the balance equation of entropy in open quantum evolutions and their response as we go from a Markovian to a non-Markovian situation. In particular, we look at the heat current and entropy production rate in the non-Markovian reduced evolution, as well as a Markovian limit of the same, experienced by one of two interacting systems immersed in a Markovian bath. The analysis naturally leads us to define a heat current deficit and an entropy production rate deficit, which are differences between the global and local versions of the corresponding quantities. The investigation leads, in certain cases, to a complementarity of the time-integrated heat current deficit and the relative entropy of entanglement between the two systems.


title: Bias-Free Scalable Gaussian Processes via Randomized Truncations

abstract: Scalable Gaussian Process methods are computationally attractive, yet introduce modeling biases that require rigorous study. This paper analyzes two common techniques: early truncated conjugate gradients (CG) and random Fourier features (RFF). We find that both methods introduce a systematic bias on the learned hyperparameters: CG tends to underfit while RFF tends to overfit. We address these issues using randomized truncation estimators that eliminate bias in exchange for increased variance. In the case of RFF, we show that the bias-to-variance conversion is indeed a trade-off: the additional variance proves detrimental to optimization. However, in the case of CG, our unbiased learning procedure meaningfully outperforms its biased counterpart with minimal additional computation.


title: Efficient Conditional GAN Transfer with Knowledge Propagation across Classes

abstract: Generative adversarial networks (GANs) have shown impressive results in both unconditional and conditional image generation. In recent literature, it is shown that pre-trained GANs, on a different dataset, can be transferred to improve the image generation from a small target data. The same, however, has not been well-studied in the case of conditional GANs (cGANs), which provides new opportunities for knowledge transfer compared to unconditional setup. In particular, the new classes may borrow knowledge from the related old classes, or share knowledge among themselves to improve the training. This motivates us to study the problem of efficient conditional GAN transfer with knowledge propagation across classes. To address this problem, we introduce a new GAN transfer method to explicitly propagate the knowledge from the old classes to the new classes. The key idea is to enforce the popularly used conditional batch normalization (BN) to learn the class-specific information of the new classes from that of the old classes, with implicit knowledge sharing among the new ones. This allows for an efficient knowledge propagation from the old classes to the new ones, with the BN parameters increasing linearly with the number of new classes. The extensive evaluation demonstrates the clear superiority of the proposed method over state-of-the-art competitors for efficient conditional GAN transfer tasks. The code is available at: https://github.com/mshahbazi72/cGANTransfer


title: Matching Point Sets with Quantum Circuit Learning

abstract: In this work, we propose a parameterised quantum circuit learning approach to point set matching problem. In contrast to previous annealing-based methods, we propose a quantum circuit-based framework whose parameters are optimised via descending the gradients w.r.t a kernel-based loss function. We formulate the shape matching problem into a distribution learning task; that is, to learn the distribution of the optimal transformation parameters. We show that this framework is able to find multiple optimal solutions for symmetric shapes and is more accurate, scalable and robust than the previous annealing-based method. Code, data and pre-trained weights are available at the project page: \href{https://hansen7.github.io/qKC}{https://hansen7.github.io/qKC}


title: A dynamic network model of societal complexity and resilience inspired by Tainter's theory of collapse

abstract: In recent years, several global events have severely disrupted economies and social structures, undermining confidence in the resilience of modern societies. While empirical evidence on the dynamics and drivers of past societal collapse is mounting, a process-based understanding of these dynamics is still in its infancy. Here we aim to identify and illustrate the underlying drivers of such societal instability or even collapse. The inspiration for this work is Joseph Tainter's theory of the "collapse of complex societies", which postulates that the complexity of societies increases as they solve problems, leading to diminishing returns on complexity investments, and ultimately to collapse. In this work we have abstracted this theory into a low-dimensional and stylised model of two classes of networked agents, hereafter referred to as "laborers" and "administrators". We numerically modeled the dynamics of societal complexity, measured as the fraction of "administrators", which is assumed to affect the productivity of connected energy-producing "laborers". We show that collapse becomes increasingly likely as the complexity of the model society continuously increases in response to external stresses that emulate Tainter's abstract notion of problems that societies must solve. We also provide an analytical approximation of the system's dominant dynamics, which matches well with the numerical experiments, and use it to study the influence on network link density, social mobility and productivity. Our work advances the understanding of social-ecological collapse and illustrates its potentially direct link to an ever-increasing societal complexity in response to external shocks or stresses via a self-reinforcing feedback.


title: Simutaneously vanishing higher derived limits without large cardinals

abstract: A question dating to Sibe Marde\v{s}i\'{c} and Andrei Prasolov's 1988 work Strong homology is not additive, and motivating a considerable amount of set theoretic work in the ensuing years, is that of whether it is consistent with the ZFC axioms for the higher derived limits $\mathrm{lim}^n$ $(n>0)$ of a certain inverse system $\mathbf{A}$ indexed by ${^\omega}\omega$ to simultaneously vanish. An equivalent formulation of this question is that of whether it is consistent for all $n$-coherent families of functions indexed by ${^\omega}\omega$ to be trivial. In this paper, we prove that, in any forcing extension given by adjoining $\beth_\omega$-many Cohen reals, $\mathrm{lim}^n \mathbf{A}$ vanishes for all $n > 0$. Our proof involves a detailed combinatorial analysis of the forcing extension and repeated applications of higher dimensional $\Delta$-system lemmas. This work removes all large cardinal hypotheses from the main result of arXiv:1907.11744 and substantially reduces the least value of the continuum known to be compatible with the simultaneous vanishing of $\mathrm{lim}^n \mathbf{A}$ for all $n > 0$.


title: On the Paradox of Certified Training

abstract: Certified defenses based on convex relaxations are an established technique for training provably robust models. The key component is the choice of relaxation, varying from simple intervals to tight polyhedra. Counterintuitively, loose interval-based training often leads to higher certified robustness than what can be achieved with tighter relaxations, which is a well-known but poorly understood paradox. While recent works introduced various improvements aiming to circumvent this issue in practice, the fundamental problem of training models with high certified robustness remains unsolved. In this work, we investigate the underlying reasons behind the paradox and identify two key properties of relaxations, beyond tightness, that impact certified training dynamics: continuity and sensitivity. Our extensive experimental evaluation with a number of popular convex relaxations provides strong evidence that these factors can explain the drop in certified robustness observed for tighter relaxations. We also systematically explore modifications of existing relaxations and discover that improving unfavorable properties is challenging, as such attempts often harm other properties, revealing a complex tradeoff. Our findings represent an important first step towards understanding the intricate optimization challenges involved in certified training.


title: Explaining Neural Scaling Laws

abstract: The population loss of trained deep neural networks often follows precise power-law scaling relations with either the size of the training dataset or the number of parameters in the network. We propose a theory that explains the origins of and connects these scaling laws. We identify variance-limited and resolution-limited scaling behavior for both dataset and model size, for a total of four scaling regimes. The variance-limited scaling follows simply from the existence of a well-behaved infinite data or infinite width limit, while the resolution-limited regime can be explained by positing that models are effectively resolving a smooth data manifold. In the large width limit, this can be equivalently obtained from the spectrum of certain kernels, and we present evidence that large width and large dataset resolution-limited scaling exponents are related by a duality. We exhibit all four scaling regimes in the controlled setting of large random feature and pretrained models and test the predictions empirically on a range of standard architectures and datasets. We also observe several empirical relationships between datasets and scaling exponents under modifications of task and architecture aspect ratio. Our work provides a taxonomy for classifying different scaling regimes, underscores that there can be different mechanisms driving improvements in loss, and lends insight into the microscopic origins of and relationships between scaling exponents.


title: Material absorption-based carrier generation model for modeling optoelectronic devices

abstract: The generation rate of photocarriers in optoelectronic materials is commonly calculated using the Poynting vector in the frequency domain. In time-domain approaches where the nonlinear coupling between electromagnetic (EM) waves and photocarriers can be accounted for, the Poynting vector model is no longer applicable. One main reason is that the photocurrent radiates low-frequency EM waves out of the spectrum of the source, e.g., terahertz (THz) waves are generated in THz photoconductive antennas. These frequency components do not contribute to the photocarrier generation since the corresponding photon energy is smaller than the optoelectronic material's bandgap energy. However, the instantaneous Poynting vector does not distinguish the power flux of different frequency components. This work proposes a material absorption-based model capable of calculating the carrier generation rate accurately in the time domain. Using the Lorentz dispersion model with poles reside in the optical frequency region, the instantaneous optical absorption, which corresponds to the power dissipation in the polarization, is calculated and used to calculate the generation rate. The Lorentz model is formulated with an auxiliary differential equation method that updates the polarization current density, from which the absorbed optical power corresponding to each Lorentz pole is directly calculated in the time domain. Examples show that the proposed model is more accurate than the Poynting vector-based model and is stable even when the generated low-frequency component is strong.


title: SDSS-IV MaNGA: Radial Gradients in Stellar Population Properties of Early-Type and Late-Type Galaxies

abstract: We derive ages, metallicities, and individual element abundances of early- and late-type galaxies (ETGs and LTGs) out to 1.5 R$_e$. We study a large sample of 1900 galaxies spanning $8.6 - 11.3 \log M/M_{\odot}$ in stellar mass, through key absorption features in stacked spectra from the SDSS-IV/MaNGA survey. We use mock galaxy spectra with extended star formation histories to validate our method for LTGs and use corrections to convert the derived ages into luminosity- and mass-weighted quantities. We find flat age and negative metallicity gradients for ETGs and negative age and negative metallicity gradients for LTGs. Age gradients in LTGs steepen with increasing galaxy mass, from $-0.05\pm0.11~\log$ Gyr/R$_e$ for the lowest mass galaxies to $-0.82\pm0.08~\log$ Gyr/R$_e$ for the highest mass ones. This strong gradient-mass relation has a slope of $-0.70\pm0.18$. Comparing local age and metallicity gradients with the velocity dispersion $\sigma$ within galaxies against the global relation with $\sigma$ shows that internal processes regulate metallicity in ETGs but not age, and vice versa for LTGs. We further find that metallicity gradients with respect to local $\sigma$ show a much stronger dependence on galaxy mass than radial metallicity gradients. Both galaxy types display flat [C/Fe] and [Mg/Fe], and negative [Na/Fe] gradients, whereas only LTGs display gradients in [Ca/Fe] and [Ti/Fe]. ETGs have increasingly steep [Na/Fe] gradients with local $\sigma$ reaching $6.50\pm0.78$ dex/$\log$ km/s for the highest masses. [Na/Fe] ratios are correlated with metallicity for both galaxy types across the entire mass range in our sample, providing support for metallicity dependent supernova yields.


title: Proximal and Federated Random Reshuffling

abstract: Random Reshuffling (RR), also known as Stochastic Gradient Descent (SGD) without replacement, is a popular and theoretically grounded method for finite-sum minimization. We propose two new algorithms: Proximal and Federated Random Reshuffing (ProxRR and FedRR). The first algorithm, ProxRR, solves composite convex finite-sum minimization problems in which the objective is the sum of a (potentially non-smooth) convex regularizer and an average of $n$ smooth objectives. We obtain the second algorithm, FedRR, as a special case of ProxRR applied to a reformulation of distributed problems with either homogeneous or heterogeneous data. We study the algorithms' convergence properties with constant and decreasing stepsizes, and show that they have considerable advantages over Proximal and Local SGD. In particular, our methods have superior complexities and ProxRR evaluates the proximal operator once per epoch only. When the proximal operator is expensive to compute, this small difference makes ProxRR up to $n$ times faster than algorithms that evaluate the proximal operator in every iteration. We give examples of practical optimization tasks where the proximal operator is difficult to compute and ProxRR has a clear advantage. Finally, we corroborate our results with experiments on real data sets.


title: Numerical Investigation of Flow Structures Around a Cylindrical Afterbody Under Supersonic Condition

abstract: Large Eddy Simulation (LES) with dynamic Smagorinsky model has been applied to numerically investigate the complicated flow structures that evolve in the near wake of a cylindrical after body aligned with a uniform Mach 2.46 flow. Mean flow field properties obtained from numerical simulations, such as axial velocity, pressure on base surface, have been compared with the experimental measurements as well as with other published results. It has been found that standard k-epsilon model fails to predict the flow properties in the recirculation region where better agreement has been observed between the data obtained from LES and experimental measurements. Flow Statistics like turbulent kinetic energy and primary Reynolds stress have also been calculated and compared with the results obtained from experiments in order to quantitatively assess the ability of LES technique to predict the turbulence properties of flow field in the highly compressible shear layer region. The data obtained from LES has been further analyzed to understand the evolution of coherent structures in the flow field. Proper Orthogonal Decomposition (POD) of the data obtained from central plane in the wake region has been performed in order to reveal the most energetic structures present in the flow field.


title: Crystalline phase discriminating neutron tomography using advanced reconstruction methods

abstract: Time-of-flight neutron imaging offers complementary attenuation contrast to X-ray computed tomography (CT), coupled with the ability to extract additional information from the variation in attenuation as a function of neutron energy (time of flight) at every point (voxel) in the image. In particular Bragg edge positions provide crystallographic information and therefore enable the identification of crystalline phases directly. Here we demonstrate Bragg edge tomography with high spatial and spectral resolution. We propose a new iterative tomographic reconstruction method with a tailored regularisation term to achieve high quality reconstruction from low-count data, where conventional filtered back-projection (FBP) fails. The regularisation acts in a separated mode for spatial and spectral dimensions and favours characteristic piece-wise constant and piece-wise smooth behaviour in the respective dimensions. The proposed method is compared against FBP and a state-of-the-art regulariser for multi-channel tomography on a multi-material phantom. The proposed new regulariser which accommodates specific image properties outperforms both conventional and state-of-the-art methods and therefore facilitates Bragg edge fitting at the voxel level. The proposed method requires significantly shorter exposure to retrieve features of interest. This in turn facilitates more efficient usage of expensive neutron beamline time and enables the full utilisation of state-of-the-art high resolution detectors.


title: Mechanical Performance of 3D Printed Interpenetrating Phase Composites with Spinodal Topologies

abstract: The mechanical response of interpenetrating phase composites (IPCs) with stochastic spinodal topologies is investigated experimentally and numerically. Model polymeric systems are fabricated by Polyjet multi-material printing, with the reinforcing phase taking the topology of a spinodal shell, and the remaining volume filled by a softer matrix. We show that spinodal shell IPCs have comparable compressive strength and stiffness to IPCs with two well-established periodic reinforcements, the Schwarz P triply periodic minimal surface (TPMS) and the octet truss-lattice, while exhibiting far less catastrophic failure and greater damage resistance, particularly at high volume fraction of reinforcing phase. The combination of high stiffness and strength and a long flat plateau after yielding makes spinodal shell IPCs a promising candidate for energy absorption and impact protection applications, where the lack of material softening upon large compressive strains can prevent sudden collapse. Importantly, in contrast with all IPCs with periodic reinforcements, spinodal shell IPCs are amenable to scalable manufacturing via self-assembly techniques.


title: A pedagogical note on the computation of relative entropy of two $n$-mode gaussian states

abstract: We present a formula for the relative entropy S(rho||sigma) of two n mode gaussian states rho, sigma in the boson Fock space. It is shown that the relative entropy has a classical and a quantum part: The classical part consists of a weighted linear combination of relative Shannon entropies of n pairs of Bernouli trials arising from the thermal state composition of the gaussian states rho and sigma. The quantum part has a sum of n terms, that are functions of the annihilation means and the covariance matrices of 1-mode marginals of the gaussian state $\rho'$, which is equivalent to \rho under a disentangling unitary gaussian symmetry operation of the state $\sigma$. A generalized formula for the Petz-Renyi relative entropy S_alpha(rho||sigma) for gaussian states \rho, \sigma is also presented. Furthermore it is shown that the Petz-Renyi relative entropy converges to the limit S(rho||sigma) as alpha increases to 1.


title: Information content in mean pairwise velocity and mean relative velocity between pairs in a triplet

abstract: Velocity field provides a complementary avenue to constrain cosmological information, either through the peculiar velocity surveys or the kinetic Sunyaev Zel'dovich effect. One of the commonly used statistics is the mean radial pairwise velocity. Here, we consider the three-point mean relative velocity, i.e. the mean relative velocities between pairs in a triplet. Using halo catalogs from the Quijote suite of N-body simulations, we first showcase how the analytical prediction for the mean relative velocities between pairs in a triplet achieve better than 4-5% accuracy using standard perturbation theory at leading order for triangular configurations with a minimum separation of $r \geq 50\ h^{-1}$Mpc. Furthermore, we present the three-point relative velocity as a novel probe of neutrino mass estimation. We explore the full cosmological information content of the halo mean pairwise velocities, and the mean relative velocities between halo pairs in a triplet. We undertake this through the Fisher-matrix formalism using 22,000 simulations from the Quijote suite, and considering all triangular configurations with a minimum and a maximum separation of $20\ h^{-1}$Mpc and $120\ h^{-1}$Mpc, respectively. We find that the mean relative velocities in a triplet allows a 1$\sigma$ neutrino mass ($M_\nu$) constraint of 0.065 eV, that is roughly 13 times better than the mean pairwise velocity constraint (0.877 eV). This information gain is not limited only to neutrino mass, but extends to other cosmological parameters: $\Omega_{\mathrm{m}}$, $\Omega_{\mathrm{b}}$, $h$, $n_{\mathrm{s}}$ and $\sigma_{8}$ achieving a gain of 8.9, 11.8, 15.5, 20.9 and 10.9 times respectively. These results illustrate the possibility of exploiting the mean three-point relative velocities for constraining the cosmological parameters accurately from future cosmic microwave background experiments and peculiar velocity surveys.


title: Fast infrared variability from the black-hole candidate MAXI J1535$-$571 and tight constraints on the modelling

abstract: We present the results regarding the analysis of the fast X-ray/infrared (IR) variability of the black-hole transient MAXI J1535$-$571. The data studied in this work consist of two strictly simultaneous observations performed with XMM-Newton (X-rays: 0.7$-$10 keV), VLT/HAWK-I ($K_{\rm s}$ band, 2.2 $\mu$m) and VLT/VISIR ($M$ and $PAH2$_$2$ bands, 4.85 and 11.88 $\mu$m respectively). The cross-correlation function between the X-ray and near-IR light curves shows a strong asymmetric anti-correlation dip at positive lags. We detect a near-IR QPO (2.5 $\sigma$) at $2.07\pm0.09$ Hz simultaneously with an X-ray QPO at approximately the same frequency ($f_0=2.25\pm0.05$). From the cross-spectral analysis a lag consistent with zero was measured between the two oscillations. We also measure a significant correlation between the average near-IR and mid-IR fluxes during the second night, but find no correlation on short timescales. We discuss these results in terms of the two main scenarios for fast IR variability (hot inflow and jet powered by internal shocks). In both cases, our preliminary modelling suggests the presence of a misalignment between disk and jet.


title: The Effects of $\Lambda$CDM Dark Matter Substructure on the Orbital Evolution of Star Clusters

abstract: We present a comprehensive study on how perturbations due to a distribution of $\Lambda$CDM dark matter subhalos can lead to star clusters deviating from their orbits. Through a large suite of massless test particle simulations, we find that (1) subhalos with masses less than $10^8 M_{\odot}$ negligibly affect test particle orbits, (2) perturbations lead to orbital deviations only in environments with substructure fractions $f_{sub} \geq 1\%$, (3) perturbations from denser subhalos produce larger orbital deviations, and (4) subhalo perturbations that are strong relative to the background tidal field lead to larger orbital deviations. To predict how the variation in test particle orbital energy $\sigma_e(t)$ increases with time, we test the applicability of theory derived from single-mass subhalo populations to populations where subhalos have a mass spectrum. We find $\sigma_e(t)$ can be predicted for test particle evolution within a mass spectrum of subhalos by assuming subhalos all have masses equal to the mean subhalo mass and by using the local mean subhalo separation to estimate the change in test particle velocities due to subhalo interactions. Furthermore, the orbital distance variation at an orbital distance $r$ can be calculated via $\sigma_r=2.98 \times 10^{-5} \pm 8 \times 10^{-8} (\rm kpc^{-1} km^{-2} s^{2}) \times r \times \sigma_e$ with a dispersion about the line of best fit equalling 0.08 kpc. Finally, we conclude that clusters that orbit within 100 kpc of Milky Way-like galaxies experience a change no greater than $2\%$ in their dissolution times.


title: Minidisk dynamics in accreting, spinning black hole binaries: Simulations in full general relativity

abstract: We perform magnetohydrodynamic simulations of accreting, equal-mass binary black holes in full general relativity focusing on the impact of black hole spin on the dynamical formation and evolution of minidisks. We find that during the late inspiral the sizes of minidisks are primarily determined by the interplay between the tidal field and the effective innermost stable orbit around each black hole. Our calculations support that a minidisk forms when the Hill sphere around each black hole is significantly larger than the black hole's effective innermost stable orbit. As the binary inspirals, the radius of the Hill sphere decreases, and minidisk sconsequently shrink in size. As a result, electromagnetic signatures associated with minidisks may be expected to gradually disappear prior to merger when there are no more stable orbits within the Hill sphere. In particular, a gradual disappearance of a hard electromagnetic component in the spectrum of such systems could provide a characteristic signature of merging black hole binaries. For a binary of given total mass, the timescale to minidisk "evaporation" should therefore depend on the black hole spins and the mass ratio. We also demonstrate that accreting binary black holes with spin have a higher efficiency for converting accretion power to jet luminosity. These results could provide new ways to estimate black hole spins in the future.


title: Deep learning approach for identification of HII regions during reionization in 21-cm observations

abstract: The upcoming Square Kilometre Array (SKA-Low) will map the distribution of neutral hydrogen during reionization, and produce a tremendous amount of 3D tomographic data. These images cubes will be subject to instrumental limitations, such as noise and limited resolution. Here we present SegU-Net, a stable and reliable method for identification of neutral and ionized regions in these images. SegU-Net is a U-Net architecture based convolutional neural network (CNN) for image segmentation. It is capable of segmenting our image data into meaningful features (ionized and neutral regions) with greater accuracy compared to previous methods. We can estimate the true ionization history from our mock observation of SKA with an observation time of 1000 h with more than 87 per cent accuracy. We also show that SegU-Net can be used to recover various topological summary statistics, such as size distributions and Betti numbers, with a relative difference of only a few per cent. These summary statistics characterise the non-Gaussian nature of the reionization process.


title: Near future MeV telescopes can discover asteroid-mass primordial black hole dark matter

abstract: Primordial black holes (PBHs), formed out of large overdensities in the early Universe, are a viable dark matter (DM) candidate over a broad range of masses. Ultra-light, asteroid-mass PBHs with masses around $10^{17}$ g are particularly interesting as current observations allow them to constitute the entire DM density. PBHs in this mass range emit $\sim$ MeV photons via Hawking radiation which can directly be detected by the gamma ray telescopes, such as the upcoming AMEGO. In this work we forecast how well an instrument with the sensitivity of AMEGO will be able to detect, or rule out, PBHs as a DM candidate, by searching for their evaporating signature when marginalizing over the Galactic and extra-Galactic gamma-ray backgrounds. We find that an instrument with the sensitivity of AMEGO could exclude non-rotating PBHs as the only DM component for masses up to $7 \times 10^{17}$ g at 95% confidence level (C.L.) for a monochromatic mass distribution, improving upon current bounds by nearly an order of magnitude. The forecasted constraints are more stringent for PBHs that have rotation, or which follow extended mass distributions.


title: Simulation of three-dimensional quantum systems with projected entangled-pair states

abstract: Tensor network algorithms have proven to be very powerful tools for studying one- and two-dimensional quantum many-body systems. However, their application to three-dimensional (3D) quantum systems has so far been limited, mostly because the efficient contraction of a 3D tensor network is very challenging. In this paper we develop and benchmark two contraction approaches for infinite projected entangled-pair states (iPEPS) in 3D. The first approach is based on a contraction of a finite cluster of tensors including an effective environment to approximate the full 3D network. The second approach performs a full contraction of the network by first iteratively contracting layers of the network with a boundary iPEPS, followed by a contraction of the resulting quasi-2D network using the corner transfer matrix renormalization group. Benchmark data for the Heisenberg and Bose-Hubbard models on the cubic lattice show that the algorithms provide competitive results compared to other approaches, making iPEPS a promising tool to study challenging open problems in 3D.


title: CP Symmetry and Symplectic Modular Invariance

abstract: We analyze CP symmetry in symplectic modular-invariant supersymmetric theories. We show that for genus $g\ge 3$ the definition of CP is unique, while two independent possibilities are allowed when $g\le 2$. We discuss the transformation properties of moduli, matter multiplets and modular forms in the Siegel upper half plane, as well as in invariant subspaces. We identify CP-conserving surfaces in the fundamental domain of moduli space. We make use of all these elements to build a CP and symplectic invariant model of lepton masses and mixing angles, where known data are well reproduced and observable phases are predicted in terms of a minimum number of parameters.


title: An HST Survey of Protostellar Outflow Cavities: Does Feedback Clear Envelopes?

abstract: We study protostellar envelope and outflow evolution using Hubble Space Telescope NICMOS or WFC3 images of 304 protostars in the Orion Molecular clouds. These near-IR images resolve structures in the envelopes delineated by the scattered light of the central protostars with 80 AU resolution and they complement the 1.2-870 micron spectral energy distributions obtained with the Herschel Orion Protostar Survey program (HOPS). Based on their 1.60 micron morphologies, we classify the protostars into five categories: non-detections, point sources without nebulosity, bipolar cavity sources, unipolar cavity sources, and irregulars. We find point sources without associated nebulosity are the most numerous, and show through monochromatic Monte Carlo radiative transfer modeling that this morphology occurs when protostars are observed at low inclinations or have low envelope densities. We also find that the morphology is correlated with the SED-determined evolutionary class with Class 0 protostars more likely to be non-detections, Class I protostars to show cavities and flat-spectrum protostars to be point sources. Using an edge detection algorithm to trace the projected edges of the cavities, we fit power-laws to the resulting cavity shapes, thereby measuring the cavity half-opening angles and power-law exponents. We find no evidence for the growth of outflow cavities as protostars evolve through the Class I protostar phase, in contradiction with previous studies of smaller samples. We conclude that the decline of mass infall with time cannot be explained by the progressive clearing of envelopes by growing outflow cavities. Furthermore, the low star formation efficiency inferred for molecular cores cannot be explained by envelope clearing alone.


title: The Compton Scattering Total Cross Section at Next-to-Leading Order

abstract: An analytic formula is given for the total scattering cross section of an electron and a photon at order $\alpha^3$. This includes both the double-Compton scattering real-emission contribution as well as the virtual Compton scattering part. When combined with the recent analytic result for the pair-production cross section, the complete $\alpha^3$ cross section is now known. Both the next-to-leading order calculation as well as the pair-production cross section are computed using modern multiloop calculation techniques, where cut diagrams are decomposed into a set of master integrals that are then computed using differential equations.


title: Double detonations of sub-Chandrasekhar mass CO white dwarfs: Can different core and He shell masses explain variations of Type Ia supernovae?

abstract: Sub-Chandrasekhar mass carbon-oxygen white dwarfs (CO WDs) with a surface helium (He) shell have been proposed as progenitors of Type Ia supernovae (SNe Ia). If true, the resulting thermonuclear explosions should be able to account for at least some of the range of SNe Ia observables. To study this, we conduct a parameter study based on 3D simulations of double detonations in CO WDs with a He shell, assuming different core and shell masses. An admixture of C to the shell and solar metallicity are included in the models. The hydrodynamic simulations are carried out using the AREPO code. This allows us to follow the He shell detonation with high numerical resolution, and improves the reliability of predicted nucleosynthetic shell detonation yields. The addition of C to the shell leads to a lower production of 56Ni while including solar metallicity increases the production of IMEs. The production of higher mass elements is further shifted to stable isotopes at solar metallicity. Moreover, we find different core detonation ignition mechanisms depending on the core and shell mass configuration. This has an influence on the ejecta structure. We present the bolometric light curves predicted from our explosion simulations using the radiative transfer code ARTIS, and make comparisons with SNe Ia data. The bolometric light curves of our models show a range of brightnesses, able to account for sub-luminous to normal SNe Ia. We show the model bolometric width-luminosity relation compared to data for a range of viewing angles. We find that, on average, our brighter models lie within the observed data. The ejecta asymmetries produce a wide distribution of observables, which might account for outliers in the data. However, the models overestimate the extent of this compared to data. We also find the bolometric decline rate over 40 days appears systematically faster than data. (abridged)


title: An enquiry on the origins of N-rich stars in the inner Galaxy basedon APOGEE chemical compositions

abstract: Recent evidence based on APOGEE data for stars within a few kpc of the Galactic centre suggests that dissolved globular clusters (GCs) contribute significantly to the stellar mass budget of the inner halo. In this paper we enquire into the origins of tracers of GC dissolution, N-rich stars, that are located in the inner 4 kpc of the Milky Way. From an analysis of the chemical compositions of these stars we establish that about 30% of the N-rich stars previously identified in the inner Galaxy may have an accreted origin. This result is confirmed by an analysis of the kinematic properties of our sample. The specific frequency of N-rich stars is quite large in the accreted population, exceeding that of its in situ counterparts by near an order of magnitude, in disagreement with predictions from numerical simulations. We hope that our numbers provide a useful test to models of GC formation and destruction.


title: Quantum information dynamics in a high-dimensional parity-time-symmetric system

abstract: Non-Hermitian systems with parity-time ($\mathcal{PT}$) symmetry give rise to exceptional points (EPs) with exceptional properties that arise due to the coalescence of eigenvectors. Such systems have been extensively explored in the classical domain, where second or higher order EPs have been proposed or realized. In contrast, quantum information studies of $\mathcal{PT}$-symmetric systems have been confined to systems with a two-dimensional Hilbert space. Here by using a single-photon interferometry setup, we simulate quantum dynamics of a four-dimensional $\mathcal{PT}$-symmetric system across a fourth-order exceptional point. By tracking the coherent, non-unitary evolution of the density matrix of the system in $\mathcal{PT}$-symmetry unbroken and broken regions, we observe the entropy dynamics for both the entire system, and the gain and loss subsystems. Our setup is scalable to the higher-dimensional $\mathcal{PT}$-symmetric systems, and our results point towards the rich dynamics and critical properties.


title: The search for low-mass axion dark matter with ABRACADABRA-10cm

abstract: Two of the most pressing questions in physics are the microscopic nature of the dark matter that comprises 84% of the mass in the universe and the absence of a neutron electric dipole moment. These questions would be resolved by the existence of a hypothetical particle known as the quantum chromodynamics (QCD) axion. In this work, we probe the hypothesis that axions constitute dark matter, using the ABRACADABRA-10cm experiment in a broadband configuration, with world-leading sensitivity. We find no significant evidence for axions, and we present 95% upper limits on the axion-photon coupling down to the world-leading level $g_{a\gamma\gamma}<3.2 \times10^{-11}$ GeV$^{-1}$, representing one of the most sensitive searches for axions in the 0.41 - 8.27 neV mass range. Our work paves a direct path for future experiments capable of confirming or excluding the hypothesis that dark matter is a QCD axion in the mass range motivated by String Theory and Grand Unified Theories.


title: A note on twisted group rings and semilinearization

abstract: In this short note, we construct a right adjoint to the functor which associates to a ring $R$ equipped with a group action its twisted group ring. This right adjoint admits an interpretation as semilinearization, in that it sends an $R$-module to the group of semilinear $R$-module automorphisms of the module. As an immediate corollary, we provide a novel proof of the classical observation that modules over a twisted group ring are modules over the base ring together with a semilinear action.


title: Homotopy Mackey functors of equivariant algebraic $K$-theory

abstract: Given a finite group $G$ acting on a ring $R$, Merling constructed an equivariant algebraic $K$-theory $G$-spectrum, and work of Malkiewich and Merling, as well as work of Barwick, provides an interpretation of this construction as a spectral Mackey functor. This construction is powerful, but highly categorical; as a result the Mackey functors comprising the homotopy are not obvious from the construction and have therefore not yet been calculated. In this work, we provide a computation of the homotopy Mackey functors of equivariant algebraic $K$-theory in terms of a purely algebraic construction. In particular, we construct Mackey functors out of the $n$th algebraic $K$-groups of group rings whose multiplication is twisted by the group action. Restrictions and transfers for these functors admit a tractable algebraic description in that they arise from restriction and extension of scalars along module categories of twisted group rings. In the case where the group action is trivial, our construction recovers work of Dress and Kuku from the 1980's which constructs Mackey functors out of the algebraic $K$-theory of group rings. We develop many families of examples of Mackey functors, both new and old, including $K$-theory of endomorphism rings, the $K$-theory of fixed subrings of Galois extensions, and (topological) Hochschild homology of twisted group rings.


title: Neural Network Libraries: A Deep Learning Framework Designed from Engineers' Perspectives

abstract: While there exist a plethora of deep learning tools and frameworks, the fast-growing complexity of the field brings new demands and challenges, such as more flexible network design, speedy computation on distributed setting, and compatibility between different tools. In this paper, we introduce Neural Network Libraries (https://nnabla.org), a deep learning framework designed from engineer's perspective, with emphasis on usability and compatibility as its core design principles. We elaborate on each of our design principles and its merits, and validate our attempts via experiments.


title: SOAR: A Synthesis Approach for Data Science API Refactoring

abstract: With the growth of the open-source data science community, both the number of data science libraries and the number of versions for the same library are increasing rapidly. To match the evolving APIs from those libraries, open-source organizations often have to exert manual effort to refactor the APIs used in the code base. Moreover, due to the abundance of similar open-source libraries, data scientists working on a certain application may have an abundance of libraries to choose, maintain and migrate between. The manual refactoring between APIs is a tedious and error-prone task. Although recent research efforts were made on performing automatic API refactoring between different languages, previous work relies on statistical learning with collected pairwise training data for the API matching and migration. Using large statistical data for refactoring is not ideal because such training data will not be available for a new library or a new version of the same library. We introduce Synthesis for Open-Source API Refactoring (SOAR), a novel technique that requires no training data to achieve API migration and refactoring. SOAR relies only on the documentation that is readily available at the release of the library to learn API representations and mapping between libraries. Using program synthesis, SOAR automatically computes the correct configuration of arguments to the APIs and any glue code required to invoke those APIs. SOAR also uses the interpreter's error messages when running refactored code to generate logical constraints that can be used to prune the search space. Our empirical evaluation shows that SOAR can successfully refactor 80% of our benchmarks corresponding to deep learning models with up to 44 layers with an average run time of 97.23 seconds, and 90% of the data wrangling benchmarks with an average run time of 17.31 seconds.


title: Operational Annotations: A new method for sequential program verification

abstract: I present a new method for specifying and verifying the partial correctness of sequential programs. The key observation is that, in Hoare logic, assertions are used as selectors of states, that is, an assertion specifies the set of program states that satisfy the assertion. Hence, the usual meaning of the partial correctness Hoare triple $\{f\}\{P\}\{g\}$: if execution is started in any of the states that satisfy assertion $f$, then, upon termination, the resulting state will be some state that satisfies assertion $g$. There are of course other ways to specify a set of states. I propose to replace assertions by terminating programs: a program $A$ specifies a set of states as follows: we start $A$ in any state whatsoever, and all the states that $A$ may terminate in constitute the specified set. I call this set the post-states of $A$. I introduce the operational triple $[A]\, P\, [B]$ to mean: if execution of $P$ is started in any post-state of $A$, then upon termination, the resulting state will be some post-state of $B$. Here, $A$ is the pre-program, and plays the role of a pre-condition, and $B$ is the post-program, and plays the role of a post-condition.


title: Continuous Jets and Backflow Models for the Formation of W50/SS433 in Magnetohydrodynamics Simulations

abstract: The formation mechanism of the W50/SS433 complex has long been a mystery. We propose a new scenario in which the SS433 jets themselves form the W50/SS433 system. We carry out magnetohydrodynamics simulations of two-side jet propagation using the public code CANS+. As found in previous jet studies, when the propagating jet is lighter than the surrounding medium, the shocked plasma flows back from the jet tip to the core. We find that the morphology of light jets is spheroidal at early times, and afterward, the shell and wings are developed by the broadening spherical cocoon. The morphology strongly depends on the density ratio of the injected jet to the surrounding medium. Meanwhile, the ratio of the lengths of the two-side jets depends only on the density profile of the surrounding medium. We also find that most of the jet kinetic energy is dissipated at the oblique shock formed by the interaction between the backflow and beam flow, rather than at the jet terminal shock. The position of the oblique shock is spatially consistent with the X-ray and TeV gamma-ray hotspots of W50.


title: A novel method for object detection using deep learning and CAD models

abstract: Object Detection (OD) is an important computer vision problem for industry, which can be used for quality control in the production lines, among other applications. Recently, Deep Learning (DL) methods have enabled practitioners to train OD models performing well on complex real world images. However, the adoption of these models in industry is still limited by the difficulty and the significant cost of collecting high quality training datasets. On the other hand, when applying OD to the context of production lines, CAD models of the objects to be detected are often available. In this paper, we introduce a fully automated method that uses a CAD model of an object and returns a fully trained OD model for detecting this object. To do this, we created a Blender script that generates realistic labeled datasets of images containing the object, which are then used for training the OD model. The method is validated experimentally on two practical examples, showing that this approach can generate OD models performing well on real images, while being trained only on synthetic images. The proposed method has potential to facilitate the adoption of object detection models in industry as it is easy to adapt for new objects and highly flexible. Hence, it can result in significant costs reduction, gains in productivity and improved products quality.


title: Atom-dimer and dimer-dimer scatterings in a spin-orbit coupled Fermi gas

abstract: Using the diagrammatic approach, here we study how spin-orbit coupling (SOC) affects the fermion-dimer and dimer-dimer scattering lengths in the Born approximation, and benchmark their accuracy with the higher-order approximations. We consider both isotropic and Rashba couplings in three dimensions, and show that the Born approximation gives accurate results in the $1/(m \alpha a_s) \ll -1$ limit, where $m$ is the mass of the fermions, $\alpha$ is the strength of the SOC, and $a_s$ is the $s$-wave scattering length between fermions. This is because the higher-loop contributions form a perturbative series in the $1/(m \alpha a_s) < 0$ region that is controlled by the smallness of the residue $Z$ of the dimer propagator. In sharp contrast, since $Z$ grows with the square-root of the binding energy of the dimer in the $1/(m \alpha a_s) > 0$ region, all of the higher-loop contributions are of similar order.


title: Contrastive latent variable modeling with application to case-control sequencing experiments

abstract: High-throughput RNA-sequencing (RNA-seq) technologies are powerful tools for understanding cellular state. Often it is of interest to quantify and summarize changes in cell state that occur between experimental or biological conditions. Differential expression is typically assessed using univariate tests to measure gene-wise shifts in expression. However, these methods largely ignore changes in transcriptional correlation. Furthermore, there is a need to identify the low-dimensional structure of the gene expression shift to identify collections of genes that change between conditions. Here, we propose contrastive latent variable models designed for count data to create a richer portrait of differential expression in sequencing data. These models disentangle the sources of transcriptional variation in different conditions, in the context of an explicit model of variation at baseline. Moreover, we develop a model-based hypothesis testing framework that can test for global and gene subset-specific changes in expression. We test our model through extensive simulations and analyses with count-based gene expression data from perturbation and observational sequencing experiments. We find that our methods can effectively summarize and quantify complex transcriptional changes in case-control experimental sequencing data.


title: Towards Robust Visual Information Extraction in Real World: New Dataset and Novel Solution

abstract: Visual information extraction (VIE) has attracted considerable attention recently owing to its various advanced applications such as document understanding, automatic marking and intelligent education. Most existing works decoupled this problem into several independent sub-tasks of text spotting (text detection and recognition) and information extraction, which completely ignored the high correlation among them during optimization. In this paper, we propose a robust visual information extraction system (VIES) towards real-world scenarios, which is a unified end-to-end trainable framework for simultaneous text detection, recognition and information extraction by taking a single document image as input and outputting the structured information. Specifically, the information extraction branch collects abundant visual and semantic representations from text spotting for multimodal feature fusion and conversely, provides higher-level semantic clues to contribute to the optimization of text spotting. Moreover, regarding the shortage of public benchmarks, we construct a fully-annotated dataset called EPHOIE (https://github.com/HCIILAB/EPHOIE), which is the first Chinese benchmark for both text spotting and visual information extraction. EPHOIE consists of 1,494 images of examination paper head with complex layouts and background, including a total of 15,771 Chinese handwritten or printed text instances. Compared with the state-of-the-art methods, our VIES shows significant superior performance on the EPHOIE dataset and achieves a 9.01% F-score gain on the widely used SROIE dataset under the end-to-end scenario.


title: Revisiting the details when evaluating a visual tracker

abstract: Visual tracking algorithms are naturally adopted in various applications, there have been several benchmarks and many tracking algorithms, more expected to appear in the future. In this report, I focus on single object tracking and revisit the details of tracker evaluation based on widely used OTB\cite{otb} benchmark by introducing a simpler, accurate, and extensible method for tracker evaluation and comparison. Experimental results suggest that there may not be an absolute winner among tracking algorithms. We have to perform detailed analysis to select suitable trackers for use cases.


title: The Teukolsky--Starobinsky constants: facts and fictions

abstract: The Teukolsky Master Equation describes the dynamics of massless fields with spin on a Kerr black hole. Under separation of variables, spin-reversal for this equation is accomplished through the so-called Teukolsky--Starobinsky identities. These identities are associated to the so-called Teukolsky--Starobinsky constants, which are spin-dependent. We collect some properties of the Teukolsky--Starobinsky constants and dispel some myths present in the literature. We show that, contrary to popular belief, these constants can be negative for spin larger than 2. Such fields thus exhibit a novel form of energy amplification which occurs for non-superradiant frequencies.


title: Learning Deep Neural Networks under Agnostic Corrupted Supervision

abstract: Training deep neural models in the presence of corrupted supervision is challenging as the corrupted data points may significantly impact the generalization performance. To alleviate this problem, we present an efficient robust algorithm that achieves strong guarantees without any assumption on the type of corruption and provides a unified framework for both classification and regression problems. Unlike many existing approaches that quantify the quality of the data points (e.g., based on their individual loss values), and filter them accordingly, the proposed algorithm focuses on controlling the collective impact of data points on the average gradient. Even when a corrupted data point failed to be excluded by our algorithm, the data point will have a very limited impact on the overall loss, as compared with state-of-the-art filtering methods based on loss values. Extensive experiments on multiple benchmark datasets have demonstrated the robustness of our algorithm under different types of corruption.


title: Multivariate Max-Stable Processes and Homogeneous Functionals

abstract: Multivariate max-stable processes are important for both theoretical investigations and various statistical applications motivated by the fact that these are limiting processes, for instance of stationary multivariate regularly varying time series, [1]. In this contribution we explore the relation between homogeneous functionals and multivariate max-stable processes and discuss the connections between multivariate max-stable process and zonoid / max-zonoid equivalence. We illustrate our results considering Brown-Resnick and Smith processes.


title: Kronecker-factored Quasi-Newton Methods for Deep Learning

abstract: Second-order methods have the capability of accelerating optimization by using much richer curvature information than first-order methods. However, most are impractical for deep learning, where the number of training parameters is huge. In Goldfarb et al. (2020), practical quasi-Newton methods were proposed that approximate the Hessian of a multilayer perceptron (MLP) model by a layer-wise block diagonal matrix where each layer's block is further approximated by a Kronecker product corresponding to the structure of the Hessian restricted to that layer. Here, we extend these methods to enable them to be applied to convolutional neural networks (CNNs), by analyzing the Kronecker-factored structure of the Hessian matrix of convolutional layers. Several improvements to the methods in Goldfarb et al. (2020) are also proposed that can be applied to both MLPs and CNNs. These new methods have memory requirements comparable to first-order methods and much less per-iteration time complexity than those in Goldfarb et al. (2020). Moreover, convergence results are proved for a variant under relatively mild conditions. Finally, we compared the performance of our new methods against several state-of-the-art (SOTA) methods on MLP autoencoder and CNN problems, and found that they outperformed the first-order SOTA methods and performed comparably to the second-order SOTA methods.


title: Simple production of cellulose nanofibril microcapsules and the rheology of their suspensions

abstract: Microcapsules are commonly used in applications ranging from therapeutics to personal care products due to their ability to deliver encapsulated species through their porous shells. Here, we demonstrate a simple and scalable approach to fabricate microcapsules with porous shells by interfacial complexation of cellulose nanofibrils and oleylamine, and investigate the rheological properties of suspensions of the resulting microcapsules. The suspensions of neat capsules are viscous liquids whose viscosity increases with volume fraction according to a modified Kreiger-Dougherty relation with a maximum packing fraction of 0.73 and an intrinsic viscosity of 4. When polyacrylic acid (PAA) is added to the internal phase of the microcapsule, however, the suspensions become elastic and display yield stresses with power-law dependencies on capsule volume fraction and PAA concentration. The elasticity appears to originate from associative interactions between microcapsules induced by PAA that resides within the microcapsule shells. These results demonstrate that it is possible to tune the rheological properties of microcapsule suspensions by changing only the composition of the internal phase, thereby providing a novel method to tailor complex fluid rheology.


title: "Shaking in 5 seconds!" A Voluntary Smartphone-based Earthquake Early Warning System

abstract: Public earthquake early warning systems have the potential to reduce individual risk by warning people of an incoming tremor but their development has been hampered by costly infrastructure. Furthermore, users' understanding of such a service and their reactions to warnings remains poorly studied. The smartphone app of the Earthquake Network initiative turns users' smartphones into motion detectors and provides the first example of purely smartphone-based earthquake early warnings, without the need for dedicated seismic station infrastructure and operating in multiple countries. We demonstrate here that early warnings have been emitted in multiple countries even for damaging shaking levels and so this offers an alternative in the many regions unlikely to be covered by conventional early warning systems in the foreseeable future. We also show that although warnings are understood and appreciated by users, notably to get psychologically prepared, only a fraction take protective actions such as "drop, cover and hold".


title: Appearance of Random Matrix Theory in Deep Learning

abstract: We investigate the local spectral statistics of the loss surface Hessians of artificial neural networks, where we discover excellent agreement with Gaussian Orthogonal Ensemble statistics across several network architectures and datasets. These results shed new light on the applicability of Random Matrix Theory to modelling neural networks and suggest a previously unrecognised role for it in the study of loss surfaces in deep learning. Inspired by these observations, we propose a novel model for the true loss surfaces of neural networks, consistent with our observations, which allows for Hessian spectral densities with rank degeneracy and outliers, extensively observed in practice, and predicts a growing independence of loss gradients as a function of distance in weight-space. We further investigate the importance of the true loss surface in neural networks and find, in contrast to previous work, that the exponential hardness of locating the global minimum has practical consequences for achieving state of the art performance.


title: Discovery of Options via Meta-Learned Subgoals

abstract: Temporal abstractions in the form of options have been shown to help reinforcement learning (RL) agents learn faster. However, despite prior work on this topic, the problem of discovering options through interaction with an environment remains a challenge. In this paper, we introduce a novel meta-gradient approach for discovering useful options in multi-task RL environments. Our approach is based on a manager-worker decomposition of the RL agent, in which a manager maximises rewards from the environment by learning a task-dependent policy over both a set of task-independent discovered-options and primitive actions. The option-reward and termination functions that define a subgoal for each option are parameterised as neural networks and trained via meta-gradients to maximise their usefulness. Empirical analysis on gridworld and DeepMind Lab tasks show that: (1) our approach can discover meaningful and diverse temporally-extended options in multi-task RL domains, (2) the discovered options are frequently used by the agent while learning to solve the training tasks, and (3) that the discovered options help a randomly initialised manager learn faster in completely new tasks.


title: Approximation Bounds for Sparse Programs

abstract: We show that sparsity constrained optimization problems over low dimensional spaces tend to have a small duality gap. We use the Shapley-Folkman theorem to derive both data-driven bounds on the duality gap, and an efficient primalization procedure to recover feasible points satisfying these bounds. These error bounds are proportional to the rate of growth of the objective with the target cardinality, which means in particular that the relaxation is nearly tight as soon as the target cardinality is large enough so that only uninformative features are added.


title: Edge Minimizing the Student Conflict Graph

abstract: In many schools, courses are given in sections. Prior to timetabling students need to be assigned to individual sections. We give a hybrid approximation sectioning algorithm that minimizes the number of edges (potential conflicts) in the student conflict graph (SCG). We start with a greedy algorithm to obtain a starting solution and then continue with a constraint programming based algorithm (CP-SAT) that reduces the number of edges. We apply the sectioning algorithm to a highly constrained timetabling model which we specify.


title: Hybrid phonetic-neural model for correction in speech recognition systems

abstract: Automatic speech recognition (ASR) is a relevant area in multiple settings because it provides a natural communication mechanism between applications and users. ASRs often fail in environments that use language specific to particular application domains. Some strategies have been explored to reduce errors in closed ASRs through post-processing, particularly automatic spell checking, and deep learning approaches. In this article, we explore using a deep neural network to refine the results of a phonetic correction algorithm applied to a telesales audio database. The results exhibit a reduction in the word error rate (WER), both in the original transcription and in the phonetic correction, which shows the viability of deep learning models together with post-processing correction strategies to reduce errors made by closed ASRs in specific language domains.


title: Thermal conductivity of CaF$_{2}$ at high pressure

abstract: We study the thermal transport properties of three CaF$_{2}$ polymorphs up to a pressure of 30 GPa using first-principle calculations and an interatomic potential based on machine learning. The lattice thermal conductivity $\kappa$ is computed by iteratively solving the linearized Boltzmann transport equation (BTE) and by taking into account three-phonon scattering. Overall, $\kappa$ increases nearly linearly with pressure, and we show that the recently discovered $\delta$-phase with $P\bar{6}2m$ symmetry and the previously known $\gamma$-CaF$_{2}$ high-pressure phase have significantly lower lattice thermal conductivities than the ambient-thermodynamic cubic fluorite ($Fm\bar{3}m$) structure. We argue that the lower $\kappa$ of these two high-pressure phases stems mainly due to a lower contribution of acoustic modes to $\kappa$ as a result of their small group velocities. We further show that the phonon mean free paths are very short for the $P\bar{6}2m$ and $Pnma$ structures at high temperatures, and resort to the Cahill-Pohl model to assess the lower limit of thermal conductivity in these domains.


title: The Importance of Being a Band: Finite-Sample Exact Distribution-Free Prediction Sets for Functional Data

abstract: Functional Data Analysis represents a field of growing interest in statistics. Despite several studies have been proposed leading to fundamental results, the problem of obtaining valid and efficient prediction sets has not been thoroughly covered. Indeed, the great majority of methods currently in the literature rely on strong distributional assumptions (e.g, Gaussianity), dimension reduction techniques and/or asymptotic arguments. In this work, we propose a new nonparametric approach in the field of Conformal Prediction based on a new family of nonconformity measures inducing conformal predictors able to create closed-form finite-sample valid or exact prediction sets under very minimal distributional assumptions. In addition, our proposal ensures that the prediction sets obtained are bands, an essential feature in the functional setting that allows the visualization and interpretation of such sets. The procedure is also fast, scalable, does not rely on functional dimension reduction techniques and allows the user to select different nonconformity measures depending on the problem at hand always obtaining valid bands. Within this family of measures, we propose also a specific measure leading to prediction bands asymptotically no less efficient than those with constant width.


title: Realizable Universal Adversarial Perturbations for Malware

abstract: Machine learning classifiers are vulnerable to adversarial examples -- input-specific perturbations that manipulate models' output. Universal Adversarial Perturbations (UAPs), which identify noisy patterns that generalize across the input space, allow the attacker to greatly scale up the generation of such examples. Although UAPs have been explored in application domains beyond computer vision, little is known about their properties and implications in the specific context of realizable attacks, such as malware, where attackers must satisfy challenging problem-space constraints. In this paper we explore the challenges and strengths of UAPs in the context of malware classification. We generate sequences of problem-space transformations that induce UAPs in the corresponding feature-space embedding and evaluate their effectiveness across different malware domains. Additionally, we propose adversarial training-based mitigations using knowledge derived from the problem-space transformations, and compare against alternative feature-space defenses. Our experiments limit the effectiveness of a white box Android evasion attack to ~20% at the cost of ~3% TPR at 1% FPR. We additionally show how our method can be adapted to more restrictive domains such as Windows malware. We observe that while adversarial training in the feature space must deal with large and often unconstrained regions, UAPs in the problem space identify specific vulnerabilities that allow us to harden a classifier more effectively, shifting the challenges and associated cost of identifying new universal adversarial transformations back to the attacker.


title: Thermodynamic optimality of glycolytic oscillations

abstract: Temporal order in living matters reflects the self-organizing nature of dynamical processes driven out of thermodynamic equilibrium. Because of functional reason, the period of a biochemical oscillation must be tuned to a specific value with precision; however, according to the thermodynamic uncertainty relation (TUR), the precision of oscillatory period is constrained by the thermodynamic cost of generating it. After reviewing the basics of chemical oscillations using Brusselator as a model system, we study the glycolytic oscillation generated by octameric phosphofructokinase (PFK), which is known to display a period of several minutes. By exploring the phase space of glycolytic oscillations, we find that the glycolytic oscillation under the cellular condition is realized in a cost effective manner. Specifically, over the biologically relevant range of parameter values of glycolysis and octameric PFK, the entropy production from the glycolytic oscillation is minimal when the oscillation period is (5 - 10) minutes. Further, the glycolytic oscillation is found at work near the phase boundary of limit cycles, suggesting that a moderate increase of glucose injection rate leads to the loss of oscillatory dynamics, which is reminiscent of the loss of pulsatile insulin release resulting from elevated blood glucose level.


title: Structural Information Preserving for Graph-to-Text Generation

abstract: The task of graph-to-text generation aims at producing sentences that preserve the meaning of input graphs. As a crucial defect, the current state-of-the-art models may mess up or even drop the core structural information of input graphs when generating outputs. We propose to tackle this problem by leveraging richer training signals that can guide our model for preserving input information. In particular, we introduce two types of autoencoding losses, each individually focusing on different aspects (a.k.a. views) of input graphs. The losses are then back-propagated to better calibrate our model via multi-task training. Experiments on two benchmarks for graph-to-text generation show the effectiveness of our approach over a state-of-the-art baseline. Our code is available at \url{http://github.com/Soistesimmer/AMR-multiview}.


title: Do as I mean, not as I say: Sequence Loss Training for Spoken Language Understanding

abstract: Spoken language understanding (SLU) systems extract transcriptions, as well as semantics of intent or named entities from speech, and are essential components of voice activated systems. SLU models, which either directly extract semantics from audio or are composed of pipelined automatic speech recognition (ASR) and natural language understanding (NLU) models, are typically trained via differentiable cross-entropy losses, even when the relevant performance metrics of interest are word or semantic error rates. In this work, we propose non-differentiable sequence losses based on SLU metrics as a proxy for semantic error and use the REINFORCE trick to train ASR and SLU models with this loss. We show that custom sequence loss training is the state-of-the-art on open SLU datasets and leads to 6% relative improvement in both ASR and NLU performance metrics on large proprietary datasets. We also demonstrate how the semantic sequence loss training paradigm can be used to update ASR and SLU models without transcripts, using semantic feedback alone.


title: Oscillations in a Becker-D\"oring model with injection and depletion

abstract: We study the Becker-D\"oring bubblelator, a variant of the Becker-D\"oring coagulation-fragmentation system that models the growth of clusters by gain or loss of monomers. Motivated by models of gas evolution oscillators from physical chemistry, we incorporate injection of monomers and depletion of large clusters. For a wide range of physical rates, the Becker-D\"oring system itself exhibits a dynamic phase transition as mass density increases past a critical value. We connect the Becker-D\"oring bubblelator to a transport equation coupled with an integrodifferential equation for excess monomer density by formal asymptotics in the near-critical regime. For suitable injection/depletion rates, we argue that time-periodic solutions appear via a Hopf bifurcation. Numerics confirm that the generation and removal of large clusters can become desynchronized, leading to temporal oscillations associated with bursts of large-cluster nucleation.


title: A Hybrid Variance-Reduced Method for Decentralized Stochastic Non-Convex Optimization

abstract: This paper considers decentralized stochastic optimization over a network of $n$ nodes, where each node possesses a smooth non-convex local cost function and the goal of the networked nodes is to find an $\epsilon$-accurate first-order stationary point of the sum of the local costs. We focus on an online setting, where each node accesses its local cost only by means of a stochastic first-order oracle that returns a noisy version of the exact gradient. In this context, we propose a novel single-loop decentralized hybrid variance-reduced stochastic gradient method, called GT-HSGD, that outperforms the existing approaches in terms of both the oracle complexity and practical implementation. The GT-HSGD algorithm implements specialized local hybrid stochastic gradient estimators that are fused over the network to track the global gradient. Remarkably, GT-HSGD achieves a network topology-independent oracle complexity of $O(n^{-1}\epsilon^{-3})$ when the required error tolerance $\epsilon$ is small enough, leading to a linear speedup with respect to the centralized optimal online variance-reduced approaches that operate on a single node. Numerical experiments are provided to illustrate our main technical results.


title: Data-Driven Synthesis of Provably Sound Side Channel Analyses

abstract: We propose a data-driven method for synthesizing a static analyzer to detect side-channel information leaks in cryptographic software. Compared to the conventional way of manually crafting such a static analyzer, which can be labor intensive, error prone and suboptimal, our learning-based technique is not only automated but also provably sound. Our analyzer consists of a set of type-inference rules learned from the training data, i.e., example code snippets annotated with ground truth. Internally, we use syntax-guided synthesis (SyGuS) to generate new features and decision tree learning (DTL) to generate type-inference rules based on these features. We guarantee soundness by formally proving each learned rule via a technique called Datalog query containment checking. We have implemented our technique in the LLVM compiler and used it to detect power side channels in C programs. Our results show that, in addition to being automated and provably sound during synthesis, the learned analyzer also has the same empirical accuracy as two state-of-the-art, manually crafted analyzers while being 300X and 900X faster, respectively.


title: Precise transit and radial-velocity characterization of a resonant pair: a warm Jupiter TOI-216c and eccentric warm Neptune TOI-216b

abstract: TOI-216 hosts a pair of warm, large exoplanets discovered by the TESS Mission. These planets were found to be in or near the 2:1 resonance, and both of them exhibit transit timing variations (TTVs). Precise characterization of the planets' masses and radii, orbital properties, and resonant behavior can test theories for the origins of planets orbiting close to their stars. Previous characterization of the system using the first six sectors of TESS data suffered from a degeneracy between planet mass and orbital eccentricity. Radial velocity measurements using HARPS, FEROS, and PFS break that degeneracy, and an expanded TTV baseline from TESS and an ongoing ground-based transit observing campaign increase the precision of the mass and eccentricity measurements. We determine that TOI-216c is a warm Jupiter, TOI-216b is an eccentric warm Neptune, and that they librate in the 2:1 resonance with a moderate libration amplitude of 60 +/- 2 degrees; small but significant free eccentricity of 0.0222 +0.0005/-0.0003 for TOI-216b; and small but significant mutual inclination of 1.2-3.9 degrees (95% confidence interval). The libration amplitude, free eccentricity, and mutual inclination imply a disturbance of TOI-216b before or after resonance capture, perhaps by an undetected third planet.


title: Unpacking Human Teachers' Intentions For Natural Interactive Task Learning

abstract: Interactive Task Learning (ITL) is an emerging research agenda that studies the design of complex intelligent robots that can acquire new knowledge through natural human teacher-robot learner interactions. ITL methods are particularly useful for designing intelligent robots whose behavior can be adapted by humans collaborating with them. Various research communities are contributing methods for ITL and a large subset of this research is \emph{robot-centered} with a focus on developing algorithms that can learn online, quickly. This paper studies the ITL problem from a \emph{human-centered} perspective to provide guidance for robot design so that human teachers can naturally teach ITL robots. In this paper, we present 1) a qualitative bidirectional analysis of an interactive teaching study (N=10) through which we characterize various aspects of actions intended and executed by human teachers when teaching a robot; 2) an in-depth discussion of the teaching approach employed by two participants to understand the need for personal adaptation to individual teaching styles; and 3) requirements for ITL robot design based on our analyses and informed by a computational theory of collaborative interactions, SharedPlans.


title: ICI-Aware Parameter Estimation for MIMO-OFDM Radar via APES Spatial Filtering

abstract: We propose a novel three-stage delay-Doppler-angle estimation algorithm for a MIMO-OFDM radar in the presence of inter-carrier interference (ICI). First, leveraging the observation that spatial covariance matrix is independent of target delays and Dopplers, we perform angle estimation via the MUSIC algorithm. For each estimated angle, we next formulate the radar delay-Doppler estimation as a joint carrier frequency offset (CFO) and channel estimation problem via an APES (amplitude and phase estimation) spatial filtering approach by transforming the delay-Doppler parameterized radar channel into an unstructured form. In the final stage, delay and Doppler of each target can be recovered from target-specific channel estimates over time and frequency. Simulation results illustrate the superior performance of the proposed algorithm in high-mobility scenarios.


title: Multimodal Data Visualization and Denoising with Integrated Diffusion

abstract: We propose a method called integrated diffusion for combining multimodal datasets, or data gathered via several different measurements on the same system, to create a joint data diffusion operator. As real world data suffers from both local and global noise, we introduce mechanisms to optimally calculate a diffusion operator that reflects the combined information from both modalities. We show the utility of this joint operator in data denoising, visualization and clustering, performing better than other methods to integrate and analyze multimodal data. We apply our method to multi-omic data generated from blood cells, measuring both gene expression and chromatin accessibility. Our approach better visualizes the geometry of the joint data, captures known cross-modality associations and identifies known cellular populations. More generally, integrated diffusion is broadly applicable to multimodal datasets generated in many medical and biological systems.


title: On automatic extraction of on-street parking spaces using park-out events data

abstract: This article proposes two different approaches to automatically create a map for valid on-street car parking spaces. For this, we use car sharing park-out events data. The first one uses spatial aggregation and the second a machine learning algorithm. For the former, we chose rasterization and road sectioning; for the latter we chose decision trees. We compare the results of these approaches and discuss their advantages and disadvantages. Furthermore, we show our results for a neighborhood in the city of Berlin and report a classification accuracy of 91.6\% on the original imbalanced data. Finally, we discuss further work; from gathering more data over a longer period of time to fitting spatial Gaussian densities to the data and the usage of apps for manual validation and annotation of parking spaces to improve ground truth data.


title: Stochastic Gradient Langevin Dynamics with Variance Reduction

abstract: Stochastic gradient Langevin dynamics (SGLD) has gained the attention of optimization researchers due to its global optimization properties. This paper proves an improved convergence property to local minimizers of nonconvex objective functions using SGLD accelerated by variance reductions. Moreover, we prove an ergodicity property of the SGLD scheme, which gives insights on its potential to find global minimizers of nonconvex objectives.


title: Charmonium production in pp, p+Pb and Pb+Pb collisions with CMS experiment

abstract: We report the measurements of the prompt and non-prompt J/{\psi} and {\psi}(2S) nuclear modification factors (RAA) using pp and PbPb data and the RpPb using pPb data collected by CMS at the LHC. All measurements were performed via the di-muon decay channel. The RAA analysis is based on PbPb and pp data samples collected at 5.02 TeV in 2015, corresponding to integrated luminosities of 464 {\mu}b-1 and 28 pb-1. These measurements are performed in the di-muon rapidity range of |y| < 2.4 as a function of centrality, rapidity, and transverse momentum (pT) from pT = 3 GeV/c in the most forward region and up to 50 GeV/c. The RpPb analysis based on pPb data, is performed for prompt J/{\psi} and {\psi}(2S) in different kinematic bins of pT and rapidity.


title: MIMIC-IF: Interpretability and Fairness Evaluation of Deep Learning Models on MIMIC-IV Dataset

abstract: The recent release of large-scale healthcare datasets has greatly propelled the research of data-driven deep learning models for healthcare applications. However, due to the nature of such deep black-boxed models, concerns about interpretability, fairness, and biases in healthcare scenarios where human lives are at stake call for a careful and thorough examinations of both datasets and models. In this work, we focus on MIMIC-IV (Medical Information Mart for Intensive Care, version IV), the largest publicly available healthcare dataset, and conduct comprehensive analyses of dataset representation bias as well as interpretability and prediction fairness of deep learning models for in-hospital mortality prediction. In terms of interpretabilty, we observe that (1) the best performing interpretability method successfully identifies critical features for mortality prediction on various prediction models; (2) demographic features are important for prediction. In terms of fairness, we observe that (1) there exists disparate treatment in prescribing mechanical ventilation among patient groups across ethnicity, gender and age; (2) all of the studied mortality predictors are generally fair while the IMV-LSTM (Interpretable Multi-Variable Long Short-Term Memory) model provides the most accurate and unbiased predictions across all protected groups. We further draw concrete connections between interpretability methods and fairness metrics by showing how feature importance from interpretability methods can be beneficial in quantifying potential disparities in mortality predictors.