title: Automatic tumour segmentation in H&E-stained whole-slide images of the pancreas

abstract: Pancreatic cancer will soon be the second leading cause of cancer-related death in Western society. Imaging techniques such as CT, MRI and ultrasound typically help providing the initial diagnosis, but histopathological assessment is still the gold standard for final confirmation of disease presence and prognosis. In recent years machine learning approaches and pathomics pipelines have shown potential in improving diagnostics and prognostics in other cancerous entities, such as breast and prostate cancer. A crucial first step in these pipelines is typically identification and segmentation of the tumour area. Ideally this step is done automatically to prevent time consuming manual annotation. We propose a multi-task convolutional neural network to balance disease detection and segmentation accuracy. We validated our approach on a dataset of 29 patients (for a total of 58 slides) at different resolutions. The best single task segmentation network achieved a median Dice of 0.885 (0.122) IQR at a resolution of 15.56 $\mu$m. Our multi-task network improved on that with a median Dice score of 0.934 (0.077) IQR.


title: Learning to automate cryo-electron microscopy data collection with Ptolemy

abstract: Over the past decade, cryogenic electron microscopy (cryo-EM) has emerged as a primary method for determining near-native, near-atomic resolution 3D structures of biological macromolecules. In order to meet increasing demand for cryo-EM, automated methods to improve throughput and efficiency while lowering costs are needed. Currently, all high-magnification cryo-EM data collection softwares require human input and manual tuning of parameters. Expert operators must navigate low- and medium-magnification images to find good high-magnification collection locations. Automating this is non-trivial: the images suffer from low signal-to-noise ratio and are affected by a range of experimental parameters that can differ for each collection session. Here, we use various computer vision algorithms, including mixture models, convolutional neural networks, and U-Nets to develop the first pipeline to automate low- and medium-magnification targeting. Learned models in this pipeline are trained on a large internal dataset of images from real world cryo-EM data collection sessions, labeled with locations that were selected by operators. Using these models, we show that we can effectively detect and classify regions of interest in low- and medium-magnification images, and can generalize to unseen sessions, as well as to images captured using different microscopes from external facilities. We expect our open-source pipeline, Ptolemy, will be both immediately useful as a tool for automation of cryo-EM data collection, and serve as a foundation for future advanced methods for efficient and automated cryo-EM microscopy.


title: Robust End-to-End Focal Liver Lesion Detection using Unregistered Multiphase Computed Tomography Images

abstract: The computer-aided diagnosis of focal liver lesions (FLLs) can help improve workflow and enable correct diagnoses; FLL detection is the first step in such a computer-aided diagnosis. Despite the recent success of deep-learning-based approaches in detecting FLLs, current methods are not sufficiently robust for assessing misaligned multiphase data. By introducing an attention-guided multiphase alignment in feature space, this study presents a fully automated, end-to-end learning framework for detecting FLLs from multiphase computed tomography (CT) images. Our method is robust to misaligned multiphase images owing to its complete learning-based approach, which reduces the sensitivity of the model's performance to the quality of registration and enables a standalone deployment of the model in clinical practice. Evaluation on a large-scale dataset with 280 patients confirmed that our method outperformed previous state-of-the-art methods and significantly reduced the performance degradation for detecting FLLs using misaligned multiphase CT images. The robustness of the proposed method can enhance the clinical adoption of the deep-learning-based computer-aided detection system.


title: Greybody Factor for a Rotating Bardeen Black Hole by Perfect Fluid Dark Matter

abstract: In this paper, the greybody factor is studied analytically for a rotating regular Bardeen black hole surrounded by perfect fluid dark matter. Firstly, we examine the behavior of effective potential by using the radial equation of motion developed from the Klein-Gordon equation. We then consider tortoise coordinate to convert the radial equation into Schr\"odinger form equation. We solve the radial equation of motion and obtain two different asymptotic solutions in terms of hypergeometric function measured at distinct regimes so called near and far-field horizons. These solutions are smoothly matched over the whole radial coordinate in an intermediate regime to check their viability. Finally, we measure the absorption probability for massless scalar field and examine the effect of perfect fluid dark matter. It is concluded that both the effective potential and greybody factor increase with perfect fluid dark matter.


title: Improving mathematical questioning in teacher training

abstract: High-fidelity, AI-based simulated classroom systems enable teachers to rehearse effective teaching strategies. However, dialogue-oriented open-ended conversations such as teaching a student about scale factors can be difficult to model. This paper builds a text-based interactive conversational agent to help teachers practice mathematical questioning skills based on the well-known Instructional Quality Assessment. We take a human-centered approach to designing our system, relying on advances in deep learning, uncertainty quantification, and natural language processing while acknowledging the limitations of conversational agents for specific pedagogical needs. Using experts' input directly during the simulation, we demonstrate how conversation success rate and high user satisfaction can be achieved.


title: Testing Reactive Systems Using Behavioural Programming, a Model Centric Approach

abstract: Testing is a significant aspect of software development. As systems become complex and their use becomes critical to the security and the function of society, the need for testing methodologies that ensure reliability and detect faults as early as possible becomes critical. The most promising approach is the model-based approach where a model is developed that defines how the system is expected to behave and how it is meant to react. The tests are derived from the model and an analysis of the test results is conducted based on it. We will investigate the prospects of using the Behavioral Programming (BP) for a model-based testing (MBT) approach that we will develop. We will develop a natural language for representing the requirements. The model will be fed to algorithms that we will develop. This includes algorithms for the automatic creation of minimal sets of test cases that cover all of the system's requirements, analysing the results of the tests, and other tools that support the testing process. The focus of our methodology will be to find faults caused by the interaction between different requirements in ways that are difficult for the testers to detect. Specifically, we will focus our attention to concurrency issues such as deadlocks and logical race condition. We will use a variety of methods that are made possible by BP, such as non-deterministic execution of scenarios and use of in-code model-checking for building test scenarios and for finding minimal coverage of the test scenarios for the system requirements using Combinatorial Test Design (CTD) methodologies. We will develop a proof-of-concept tool kit which will allow us to demonstrate and evaluate the above mentioned capabilities. We will compare the performance of our tools with the performance of manual testers and of other model-based tools using comparison criteria that we will define and develop.


title: Sh 2-301: a blistered H II region undergoing star formation

abstract: We present a multiwavelength study of the H II region Sh 2-301 (S301) using deep optical data, near-infrared data, radio continuum data and other archival data at longer wavelengths. A cluster of young stellar objects (YSOs) is identified in the north-east (NE) direction of S301. The H{\alpha} and radio continuum images trace the distribution of the ionized gas surrounding a massive star ALS 207, and the S301 H II region is bounded by an arc-like structure of gas and dust emission in the south-eastern direction. The north-western part of S301 seems to be devoid of gas and dust emission, while the presence of molecular material between the NE cluster and the central massive star ALS 207 is found. The distribution of warm dust emission, ionized gas, and neutral hydrogen together suggests a blistered morphology of the S301 H II region powered by ALS 207, which appears to be located near the edge of the cloud. The location of the NE cluster embedded in the cold molecular cloud is found opposite to the blistered morphology. There is a noticeable age difference investigated between the massive star and the NE cluster. This age difference, pressure calculation, photodissociation regions (PDRs), and the distribution of YSOs favour the positive feedback of the massive star ALS 207 in S301. On a wider scale of S301, the H II region and the young stellar cluster are depicted toward the central region of a hub-filamentary system, which is evident in the infrared images.


title: Galactic tides and the Crater II dwarf spheroidal: a challenge to LCDM?

abstract: The unusually low velocity dispersion and large size of Crater II pose a challenge to our understanding of dwarf galaxies in the Lambda Cold Dark Matter (LCDM) cosmogony. The low velocity dispersion suggests either a dark halo mass much lower than the minimum expected from hydrogen cooling limit arguments, or one that is in the late stages of extreme tidal stripping. The tidal interpretation has been favoured in recent work and is supported by the small pericentric distances consistent with available kinematic estimates. We use N-body simulations to examine this interpretation in detail, assuming a Navarro-Frenk-White (NFW) profile for Crater II's progenitor halo. Our main finding is that, although the low velocity dispersion can indeed result from the effect of tides, the large size of Crater II is inconsistent with this hypothesis. This is because galaxies stripped to match the observed velocity dispersion are also reduced to sizes much smaller than the observed half-light radius of Crater II. Unless its size has been substantially overestimated, reconciling this system with LCDM requires that either (i) it is not bound and near equilibrium (unlikely, given its crossing time is shorter than the time elapsed since pericentre), or that (ii) its progenitor halo deviates from the assumed NFW profile. The latter alternative may signal that baryons can affect the inner halo cusp even in extremely faint dwarfs or, more intriguingly, may signal effects associated with the intimate nature of the dark matter, such as finite self-interactions, or other such deviations from the canonical LCDM paradigm.


title: DeepZipper: A Novel Deep Learning Architecture for Lensed Supernovae Identification

abstract: Large-scale astronomical surveys have the potential to capture data on large numbers of strongly gravitationally lensed supernovae (LSNe). To facilitate timely analysis and spectroscopic follow-up before the supernova fades, an LSN needs to be identified soon after it begins. To quickly identify LSNe in optical survey datasets, we designed ZipperNet, a multi-branch deep neural network that combines convolutional layers (traditionally used for images) with long short-term memory (LSTM) layers (traditionally used for time series). We tested ZipperNet on the task of classifying objects from four categories -- no lens, galaxy-galaxy lens, lensed type Ia supernova, lensed core-collapse supernova -- within high-fidelity simulations of three cosmic survey data sets -- the Dark Energy Survey (DES), Rubin Observatory's Legacy Survey of Space and Time (LSST), and a Dark Energy Spectroscopic Instrument (DESI) imaging survey. Among our results, we find that for the LSST-like dataset, ZipperNet classifies LSNe with a receiver operating characteristic area under the curve of 0.97, predicts the spectroscopic type of the lensed supernovae with 79\% accuracy, and demonstrates similarly high performance for LSNe 1-2 epochs after first detection. We anticipate that a model like ZipperNet, which simultaneously incorporates spatial and temporal information, can play a significant role in the rapid identification of lensed transient systems in cosmic survey experiments.


title: Extending the SAGA Survey (xSAGA) I: Satellite Radial Profiles as a Function of Host Galaxy Properties

abstract: We present "Extending the Satellites Around Galactic Analogs Survey" (xSAGA), a method for identifying low-$z$ galaxies on the basis of optical imaging, and results on the spatial distributions of xSAGA satellites around host galaxies. Using spectroscopic redshift catalogs from the SAGA Survey as a training data set, we have optimized a convolutional neural network (CNN) to identify $z < 0.03$ galaxies from more distant objects using image cutouts from the DESI Legacy Imaging Surveys. From the sample of $> 100,000$ CNN-selected low-$z$ galaxies, we identify $>20,000$ probable satellites located between 36-300 projected kpc from NASA-Sloan Atlas central galaxies in the stellar mass range $9.5 < \log(M_\star/M_\odot) < 11$. We characterize the incompleteness and contamination for CNN-selected samples, and apply corrections in order to estimate the true number of satellites as a function of projected radial distance from their hosts. Satellite richness depends strongly on host stellar mass, such that more massive host galaxies have more satellites, and on host morphology, such that elliptical hosts have more satellites than disky hosts with comparable stellar masses. We also find a strong inverse correlation between satellite richness and the magnitude gap between a host and its brightest satellite. The normalized satellite radial distribution between 36-300 kpc does not depend strongly on host stellar mass, morphology, or magnitude gap. The satellite abundances and radial distributions we measure are in reasonable agreement with predictions from hydrodynamic simulations. Our results deliver unprecedented statistical power for studying satellite galaxy populations, and highlight the promise of using machine learning for extending galaxy samples of wide-area surveys.


title: Supersymmetric spindles

abstract: In the context of holography, we analyse aspects of supersymmetric geometries based on two-dimensional orbifolds known as spindles. By analysing spin$^c$ spinors on a spindle with an azimuthal rotation symmetry we show that under rather general conditions there are just two possibilities, called the `twist' and the `anti-twist', which are determined by the quantized magnetic flux through the spindle. A special case of the twist is the standard topological twist which is associated with constant and chiral spinors. We construct solutions of $D=5$ and $D=4$ STU gauged supergravity theories that are dual to D3-branes and M2-branes wrapping spindles, respectively, which realize both the anti-twist, as seen before, but also the twist. For the wrapped D3-brane solutions we reproduce the central charge of the gravity solution from the dual field theory by analysing the anomaly polynomial of $\mathcal{N}=4$ SYM theory. We also discuss M5-branes wrapped on spindles both from a gravity and a field theory point of view.


title: Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback

abstract: Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$.


title: Forecasting the potential of weak lensing magnification to enhance LSST large-scale structure analyses

abstract: Recent works have shown that weak lensing magnification must be included in upcoming large-scale structure analyses, such as for the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), to avoid biasing the cosmological results. In this work we investigate whether including magnification has a positive impact on the precision of the cosmological constraints, as well as being necessary to avoid bias. We forecast this using an LSST mock catalog and a halo model to calculate the galaxy power spectra. We find that including magnification has little effect on the precision of the cosmological parameter constraints for an LSST galaxy clustering analysis, where the halo model parameters are additionally constrained by the galaxy luminosity function. In particular, we find that for the LSST gold sample ($i < 25.3$) including weak lensing magnification only improves the galaxy clustering constraint on $\Omega_{\rm{m}}$ by a factor of 1.03, and when using a very deep LSST mock sample ($i<26.5$) by a factor of 1.3. Since magnification predominantly contributes to the clustering measurement and provides similar information to that of cosmic shear, this improvement would be reduced for a combined galaxy clustering and shear analysis. We also confirm that not modelling weak lensing magnification will catastrophically bias the cosmological results from LSST. Magnification must therefore be included in LSST large-scale structure analyses even though it does not significantly enhance the precision of the cosmological constraints.


title: Neutral Gas within 20,000 Schwarzschild radii of Sagittarius A*

abstract: Murchikova et al 2019 discovered a disk of cool ionized gas within 20,000 Schwarzschild radii of the Milky Way's Galactic Center black hole Sagittarius A*. They further demonstrated that the ionizing photon flux in the region is enough to keep the disk ionized, but there is not ample excess of this radiation. This raised the possibility that some neutral gas could also be in the region shielded within the cool ionized clumps. Here we present ALMA observations of a broad 1.3 millimeter hydrogen recombination line H30alpha: n = 31 -> 30, conducted during the flyby of the S0-2 star by Sgr A*. We report that the velocity-integrated H30alpha line flux two month prior to the S0-2 pericenter passage is about 20% larger than it was one month prior to the passage. The S0-2 is a strong source of ionizing radiation moving at several thousand kilometers per second during the approach. Such a source is capable of ionising parcels of neural gas along its trajectory, resulting in variation of the recombination line spectra from epoch to epoch. We conclude that there are at least (6.6 +- 3.3) x 10^{-6} Msun of neutral gas within 20,000 Schwarzschild radii of Sgr A*.


title: Critical Dynamics: multiplicative noise fixed point in two dimensional systems

abstract: We study the critical dynamics of a real scalar field in two dimensions near a continuous phase transition. We have built up and solved Dynamical Renormalization Group equations at one-loop approximation. We have found that, different form the case $d\lesssim 4$, characterized by a Wilson-Fisher fixed point with dynamical critical exponent $z=2+ O(\epsilon^2)$, the critical dynamics is dominated by a novel multiplicative noise fixed point. The zeroes of the beta function depend on the stochastic prescription used to define the Wiener integrals. However, the critical exponents and the anomalous dimension do not depend on the prescription used. Thus, even though each stochastic prescription produces different dynamical evolutions, all of them are in the same universality class.


title: The Linear Template Fit

abstract: The estimation of parameters from data is a common problem in many areas of the physical sciences, and frequently used algorithms rely on sets of simulated data which are fit to data. In this article, an analytic solution for simulation-based parameter estimation problems is presented. The matrix formalism, termed the Linear Template Fit, calculates the best estimators for the parameters of interest. It combines a linear regression with the method of least squares. The algorithm uses only predictions calculated for a few values of the parameters of interest, which have been made available prior to its execution. The Linear Template Fit is particularly suited for performance critical applications and parameter estimation problems with computationally intense simulations, which are otherwise often limited in their usability for statistical inference. Equations for error propagation are discussed in detail and are given in closed analytic form. For the solution of problems with a nonlinear dependence on the parameters of interest, the Quadratic Template Fit is introduced. As an example application, a determination of the strong coupling constant from inclusive jet cross section data at the CERN Large Hadron Collider is studied and compared with previously published results.

categories: physics.data-an hep-ex hep-ph stat.AP stat.ME


title: Application of a Steady-State Accretion Disk Model to Spectrophotometry and High-Resolution Spectra of Two Recent FU Ori Outbursts

abstract: We apply a conventional accretion disk model to the FU Ori-type objects HBC 722 and Gaia 17bpi. Our base model is a steady-state, thin Keplerian disk featuring a modified Shakura-Sunyaev temperature profile, with each annulus radiating as an area-weighted spectrum given by a NextGen atmosphere at the appropriate temperature. We explore departures from the standard model by altering the temperature distribution in the innermost region of the disk to account for "boundary region"-like effects. We consider the overall spectral energy distribution (SED) as well as medium- and high-resolution spectra in evaluating best-fit models to the data. Parameter degeneracies are studied via a Markov Chain Monte Carlo (MCMC) parameter estimation technique. Allowing all parameters to vary, we find accretion rates for HBC 722 of $\dot{M} = 10^{-4.90} M_\odot \textrm{ yr}^{-1}\; {}^{+0.99}_{-0.40} \textrm{ dex}$ and for Gaia 17bpi of $\dot{M} = 10^{-6.70} M_\odot \textrm{ yr}^{-1}\; {}^{+0.46}_{-0.36} \textrm{ dex}$; the corresponding maximum disk temperatures are $7100_{-500}^{+300}$ K and $7900_{-400}^{+900}$ K, respectively. While the accretion rate of HBC 722 is on the same order as other FU Ori-type objects, Gaia 17bpi has a lower rate than previously reported as typical, commensurate with its lower luminosity. Alternate models that fix some disk or stellar parameters are also presented, with tighter confidence intervals on the remaining fitted parameters. In order to improve upon the somewhat large credible intervals for the $\dot{M}$ values, and make progress on boundary layer characterization, flux-calibrated ultraviolet spectroscopy is needed.


title: The effect of spiral arms on the S\'ersic photometry of galaxies

abstract: Context. The S\'ersic profile is a widely-used model to describe the surface brightness distribution of galaxies. Spiral galaxies, however, are qualitatively different from a S\'ersic model. Aims. The goal of this study is to assess how accurately the total flux and half-light radius of a galaxy with spiral arms can be recovered when fitted with a S\'ersic profile. Methods. I selected a sample of bulge-dominated galaxies with spiral arms. Using photometric data from the Hyper Suprime-Cam survey, I estimated the contribution of the spiral arms to their total flux. Then I generated simulated images of galaxies with similar characteristics, fitted them with a S\'ersic model, and quantified the error on the determination of the total flux and half-light radius. Results. Spiral arms can introduce biases on the photometry of galaxies in a way that depends on the underlying smooth surface brightness profile, the location of the arms, and the depth of the photometric data. A set of spiral arms accounting for 10% of the flux of a bulge-dominated galaxy typically causes the total flux and the half-light radius to be overestimated by 15% and 30%, respectively. This bias, however, is much smaller if the galaxy is disk-dominated. Conclusions. Galaxies with a prominent bulge and a non-zero contribution from spiral arms are the most susceptible to biases in the total flux and half-light radius, when fitted with a S\'ersic profile. If photometric measurements with high accuracy are required, then measurements over finite apertures are to be preferred over global estimates of the flux.


title: D3Net: A Unified Speaker-Listener Architecture for 3D Dense Captioning and Visual Grounding

abstract: Recent studies on dense captioning and visual grounding in 3D have achieved impressive results. Despite developments in both areas, the limited amount of available 3D vision-language data causes overfitting issues for 3D visual grounding and 3D dense captioning methods. Also, how to discriminatively describe objects in complex 3D environments is not fully studied yet. To address these challenges, we present D3Net, an end-to-end neural speaker-listener architecture that can detect, describe and discriminate. Our D3Net unifies dense captioning and visual grounding in 3D in a self-critical manner. This self-critical property of D3Net also introduces discriminability during object caption generation and enables semi-supervised training on ScanNet data with partially annotated descriptions. Our method outperforms SOTA methods in both tasks on the ScanRefer dataset, surpassing the SOTA 3D dense captioning method by a significant margin.


title: Reconstructing Orbits of Galaxies in Extreme Regions (ROGER) II: reliability of projected phase-space in our understanding of galaxy populations

abstract: We connect galaxy properties with their orbital classification by analysing a sample of galaxies with stellar mass $M_{\star} \geq 10^{8.5}h^{-1}M_\odot$ residing in and around massive and isolated galaxy clusters with mass $M_{200} > 10^{15}h^{-1}M_\odot$ at redshift $z=0$. The galaxy population is generated by applying the semi-analytic model of galaxy formation SAG on the cosmological simulation MultiDark Planck 2. We classify galaxies considering their real orbits (3D) and their projected phase-space position using the ROGER code (2D). We define five categories: cluster galaxies, galaxies that have recently fallen into a cluster, backsplash galaxies, infalling galaxies, and interloper galaxies. For each class, we analyse the $g-r$ colour, the specific star formation rate (sSFR), and the stellar age, as a function of the stellar mass. For the 3D classes, we find that cluster galaxies have the lowest sSFR, and are the reddest and the oldest, as expected from environmental effects. Backsplash galaxies have properties intermediate between the cluster and recent infaller galaxies. For each 2D class, we find an important contamination by other classes. We find it necessary to separate the galaxy populations in red and blue to perform a more realistic analysis of the 2D data. For the red population, the 2D results are in good agreement with the 3D predictions. Nevertheless, when the blue population is considered, the 2D analysis only provides reliable results for recent infallers, infalling galaxies and interloper galaxies.


title: Arresting dynamics in hardcore spin models

abstract: We study the dynamics of hardcore spin models on the square and triangular lattice, constructed by analogy to hard spheres, where the translational degrees of freedom of the spheres are replaced by orientational degrees of freedom of spins on a lattice and the packing fraction as a control parameter is replaced by an exclusion angle. In equilibrium, models on both lattices exhibit a Kosterlitz-Thouless transition at an exclusion angle $\Delta_{\rm KT}$. We devise compression protocols for hardcore spins and find that {\it any} protocol that changes the exclusion angle nonadiabatically, if endowed with only local dynamics, fails to compress random initial states beyond an angle $\Delta_{\rm J}> \Delta_{\rm KT}$. This coincides with a doubly algebraic divergence of the relaxation time of compressed states towards equilibrium. We identify a remarkably simple mechanism underpinning this divergent timescale: topological defects involved in the phase ordering kinetics of the system become incompatible with the hardcore spin constraint, leading to a vanishing defect mobility as $\Delta\rightarrow\Delta_{\rm J}$.

categories: cond-mat.stat-mech cond-mat.dis-nn cond-mat.soft cond-mat.str-el


title: Neural Head Avatars from Monocular RGB Videos

abstract: We present Neural Head Avatars, a novel neural representation that explicitly models the surface geometry and appearance of an animatable human avatar that can be used for teleconferencing in AR/VR or other applications in the movie or games industry that rely on a digital human. Our representation can be learned from a monocular RGB portrait video that features a range of different expressions and views. Specifically, we propose a hybrid representation consisting of a morphable model for the coarse shape and expressions of the face, and two feed-forward networks, predicting vertex offsets of the underlying mesh as well as a view- and expression-dependent texture. We demonstrate that this representation is able to accurately extrapolate to unseen poses and view points, and generates natural expressions while providing sharp texture details. Compared to previous works on head avatars, our method provides a disentangled shape and appearance model of the complete human head (including hair) that is compatible with the standard graphics pipeline. Moreover, it quantitatively and qualitatively outperforms current state of the art in terms of reconstruction quality and novel-view synthesis.


title: Is Approximation Universally Defensive Against Adversarial Attacks in Deep Neural Networks?

abstract: Approximate computing is known for its effectiveness in improvising the energy efficiency of deep neural network (DNN) accelerators at the cost of slight accuracy loss. Very recently, the inexact nature of approximate components, such as approximate multipliers have also been reported successful in defending adversarial attacks on DNNs models. Since the approximation errors traverse through the DNN layers as masked or unmasked, this raises a key research question-can approximate computing always offer a defense against adversarial attacks in DNNs, i.e., are they universally defensive? Towards this, we present an extensive adversarial robustness analysis of different approximate DNN accelerators (AxDNNs) using the state-of-the-art approximate multipliers. In particular, we evaluate the impact of ten adversarial attacks on different AxDNNs using the MNIST and CIFAR-10 datasets. Our results demonstrate that adversarial attacks on AxDNNs can cause 53% accuracy loss whereas the same attack may lead to almost no accuracy loss (as low as 0.06%) in the accurate DNN. Thus, approximate computing cannot be referred to as a universal defense strategy against adversarial attacks.


title: The Apollo ATCA Design for the CMS Track Finder and the Pixel Readout at the HL-LHC

abstract: The challenging conditions of the High-Luminosity LHC require tailored hardware designs for the trigger and data acquisition systems. The Apollo platform features a "Service Module" with a powerful system-on-module computer that provides standard ATCA communications and application-specific "Command Module"s with large FPGAs and high-speed optical fiber links. The CMS version of Apollo will be used for the track finder and the pixel readout. It features up to two large FPGAs and more than 100 optical links with speeds up to 25\,Gb/s. We study carefully the design and performance of the board by using customized firmware to test power consumption, heat dissipation, and optical link integrity. This paper presents the results of these performance tests, design updates, and future plans.


title: Online Bayesian Optimization for Beam Alignment in the SECAR Recoil Mass Separator

abstract: The SEparator for CApture Reactions (SECAR) is a next-generation recoil separator system at the Facility for Rare Isotope Beams (FRIB) designed for the direct measurement of capture reactions on unstable nuclei in inverse kinematics. To maximize the performance of the device, careful beam alignment to the central ion optical axis needs to be achieved. This can be difficult to attain through manual tuning by human operators without potentially leaving the system in a sub-optimal and irreproducible state. In this work, we present the first development of online Bayesian optimization with a Gaussian process model to tune an ion beam through a nuclear astrophysics recoil separator. We show that the method achieves small incoming angular deviations (0-1 mrad) in an efficient and reproducible manner that is at least 3 times faster than standard hand-tuning. This method is now routinely used for all separator tuning.


title: Geodesic equation in non-commutative gauge theory of gravity

abstract: In this work, we construct a non-commutative (NC) gauge theory of gravity for any metric with spherical symmetries, where we use a non-diagonal tetrad field. The deformed gauge potentials (tetrad fields) and the components of deformed metric are computed to the second order in the NC parameter $\Theta^{\mu\nu}$, as the application to the Schwarzschild black hole we show that the NC geometry removes the singularity at the origin of the black hole, and increase the event horizon. The non-commutativity correction to the effective potential of the Schwarzschild metric is also computed and we show how this geometry affects the stability condition which it found the NC parameter plays the same role as the mass that can be used to explain the dark matter and we show that the NC Schwarzschild space-time has new stable circular orbits appear near the event horizon that is not allowed by Schwarzschild space-time. The geodesic equations in the NC space and the corrections to the periastron advance in terms of $\Theta$ are obtained. We have also specified the problem of Mercury's perihelion and used the experimental data to estimate the NC parameter $\Theta$, then we show that $\Theta$ of the order $10^{-25}s.kg^{-1}$ gives observable corrections to the movement at a large scale. We show that the NC propriety of the spacetime appears at the High Energy.


title: A QFT for non-semisimple TQFT

abstract: We construct a family of 3d quantum field theories $\mathcal T_{n,k}^A$ that conjecturally provide a physical realization -- and derived generalization -- of non-semisimple mathematical TQFT's based on the modules for the quantum group $U_q(\mathfrak{sl}_n)$ at an even root of unity $q=\text{exp}(i\pi/k)$. The theories $\mathcal T_{n,k}^A$ are defined as topological twists of certain 3d $\mathcal N=4$ Chern-Simons-matter theories, which also admit string/M-theory realizations. They may be thought of as $SU(n)_{k-n}$ Chern-Simons theories, coupled to a twisted $\mathcal N=4$ matter sector (the source of non-semisimplicity). We show that $\mathcal T_{n,k}^A$ admits holomorphic boundary conditions supporting two different logarithmic vertex operator algebras, one of which is an $\mathfrak{sl}_n$-type Feigin-Tipunin algebra; and we conjecture that these two vertex operator algebras are related by a novel logarithmic level-rank duality. (We perform detailed computations to support the conjecture.) We thus relate the category of line operators in $\mathcal T_{n,k}^A$ to the derived category of modules for a boundary Feigin-Tipunin algebra, and -- using a logarithmic Kazhdan-Lusztig-like correspondence that has been established for $n=2$ and expected for general $n$ -- to the derived category of $U_q(\mathfrak{sl}_n)$ modules. We analyze many other key features of $\mathcal T_{n,k}^A$ and match them from quantum-group and VOA perspectives, including deformations by flat $PSL(n,\mathbb C)$ connections, one-form symmetries, and indices of (derived) genus-$g$ state spaces.


title: Tunable cryogenic THz cavity for strong light-matter coupling in complex materials

abstract: We report here the realization and commissioning of an experiment dedicated to the study of the optical properties of light matter hybrids constituted of crystalline samples embedded in an optical cavity. The experimental assembly developed offers the unique opportunity to study the weak and strong coupling regime between a tunable optical cavity in cryogenic environment and low energy degrees of freedom such as phonons, magnons or charge fluctuations. We describe here the setup developed which allows the positioning of crystalline samples in an optical cavity of different quality factor, the tuning of the cavity length at cryogenic temperatures and its optical characterization with a broadband time domain THz spectrometer (0.2-6 THz). We demonstrate the versatility of the setup by studying the vibrational strong coupling in CuGeO3 single crystal at cryogenic temperatures.


title: Baseline-free scanned-wavelength direct absorption spectroscopy: theoretical analysis and experimental validation

abstract: Applicability and accuracy of traditional scanned-wavelength direct absorption spectroscopy (SDAS) diagnostic method is largely affected by combustion state such as temperature and pressure. To resolve this problem, an innovative baseline-free scanned-wavelength direct absorption spectroscopy (BF-SDAS) method is proposed in this work. This method does not need the zero-absorption regions for baseline fitting and can exclude the influence of baseline fitting error on the final result. It is especially suitable for the cases when the scanning range is too narrow or the absorption linewidth is too large to acquire the baseline. In this paper, the BF-SDAS is first derived and theoretically analyzed based on the Fourier analysis, and then the accuracy of this method is validated by measuring the concentration of H2O molecules in a static gas absorption cell at room temperature.


title: Operator Growth from Global Out-of-time-order Correlators

abstract: In the context of chaotic quantum many-body systems, we show that operator growth, as diagnosed by out-of-time-order correlators of local operators, also leaves a sharp imprint in out-of-time-order correlators of global operators. In particular, the characteristic spacetime shape of growing local operators can be accessed using global measurements without any local control or readout. Building on an earlier conjectured phase diagram for operator growth in chaotic systems with power-law interactions, we show that existing nuclear spin data for out-of-time-order correlators of global operators are well fit by our theory. We also predict super-polynomial operator growth in dipolar systems in 3d and discuss the potential observation of this physics in future experiments with nuclear spins and ultra-cold polar molecules.


title: Lattice models from CFT on surfaces with holes I: Torus partition function via two lattice cells

abstract: We construct a one-parameter family of lattice models starting from a two-dimensional rational conformal field theory on a torus with a regular lattice of holes, each of which is equipped with a conformal boundary condition. The lattice model is obtained by cutting the surface into triangles with clipped-off edges using open channel factorisation. The parameter is given by the hole radius. At finite radius, high energy states are suppressed and the model is effectively finite. In the zero-radius limit, it recovers the CFT amplitude exactly. In the touching hole limit, one obtains a topological field theory. If one chooses a special conformal boundary condition which we call "cloaking boundary condition", then for each value of the radius the fusion category of topological line defects of the CFT is contained in the lattice model. The fact that the full topological symmetry of the initial CFT is realised exactly is a key feature of our lattice models. We provide an explicit recursive procedure to evaluate the interaction vertex on arbitrary states. As an example, we study the lattice model obtained from the Ising CFT on a torus with one hole, decomposed into two lattice cells. We numerically compare the truncated lattice model to the CFT expression obtained from expanding the boundary state in terms of the hole radius and we find good agreement at intermediate values of the radius.

categories: cond-mat.stat-mech cond-mat.str-el hep-th math-ph math.MP


title: Can scalars matter? A vector bilepton phenomenologic example

abstract: Expecting scalar contributions to be less important to a given observable signal than vector ones, in usual scenarios, is a natural intuition. Such assertion should hold for a great part of physically relevant parameter space within most models of interest. In this work, besides stressing that the proposition is not general, we show how the alternative possibility may give light to interesting phenomena.


title: A Generic Graph Sparsification Framework using Deep Reinforcement Learning

abstract: The interconnectedness and interdependence of modern graphs are growing ever more complex, causing enormous resources for processing, storage, communication, and decision-making of these graphs. In this work, we focus on the task graph sparsification: an edge-reduced graph of a similar structure to the original graph is produced while various user-defined graph metrics are largely preserved. Existing graph sparsification methods are mostly sampling-based, which introduce high computation complexity in general and lack of flexibility for a different reduction objective. We present SparRL, the first generic and effective graph sparsification framework enabled by deep reinforcement learning. SparRL can easily adapt to different reduction goals and promise graph-size-independent complexity. Extensive experiments show that SparRL outperforms all prevailing sparsification methods in producing high-quality sparsified graphs concerning a variety of objectives.


title: Theoretical Analysis of an XGBoost Framework for Product Cannibalization

abstract: This paper is an extension of our work where we presented a three-stage XGBoost algorithm for forecasting sales under product cannibalization scenario. Previously we developed the model based on our intuition and provided empirical evidence on its performance. In this study we would briefly go over the algorithm and then provide mathematical reasoning behind its working.


title: Constant Scalar Curvature Sasaki Metrics and Projective Bundles

abstract: In this paper we consider the Boothby-Wang construction over twist 1 stage 3 Bott orbifolds given in terms of the log pair $(S_{\bf n},\Delta_{\bf m})$. We give explicit constant scalar curvature (CSC) Sasaki metrics either directly from CSC K\"ahler orbifold metrics or by using the weighted extremal approach of Apostolov and Calderbank. The Sasaki 7-manifolds (orbifolds) are finitely covered by compact simply connected manifolds (orbifolds) with the rational homology of the 2-fold connected sum of $S^2\times S^5$.


title: Energy Spectra vs. Element Abundances in Solar Energetic Particles and the Roles of Magnetic Reconnection and Shock Acceleration

abstract: We reexamine the relationship between energy spectral indices and element abundance enhancements in solar energetic particle (SEP) events at energies of a few MeV/amu. We find a correlated behavior only in the largest gradual SEP4 events when all ions are accelerated from the ambient coronal plasma by shock waves driven by fast, wide coronal mass ejections (CMEs). This correlated abundance behavior can track complex time variations in the spectral indices during an event. In other (SEP3) events, CME-driven shock waves, days apart, sample seed particles from a single pool of suprathermal impulsive ions contributed earlier. Of the smaller, Fe-rich, impulsive SEP events, previously related to magnetic reconnection in solar jets, over half are subsequently reaccelerated by CME-driven shock waves (SEP2) causing typical ion intensities to have a 64% correlation with shock speed. In these SEP2 events, onset of shock acceleration is signaled by a new component in the abundances, large proton excesses. The remaining SEP1 events lack evidence of shock acceleration. However, for all these events (SEP1 - SEP3) with abundances determined by magnetic reconnection, spectra and abundances are decoupled.


title: Observation of coexisting weak localization and superconducting fluctuations in strained Sn1-xInxTe thin films

abstract: Topological superconductors have attracted tremendous excitement as they are predicted to host Majorana zero modes that can be utilized for topological quantum computing. Candidate topological superconductor Sn1-xInxTe thin films (0<x<0.3) grown by molecular beam epitaxy and strained in the (111) plane are shown to host three coexisting quantum effects: localization, antilocalization and superconducting fluctuations above the critical temperature Tc. An analysis of the normal state magnetoresistance reveals these effects. Weak localization is consistently observed in superconducting samples, indicating that superconductivity originates dominantly from trivial valence band states that may be strongly spin-orbit split. A large enhancement of the conductivity is observed above Tc, indicating that quantum coherent quasiparticle effects coexist with superconducting fluctuations. Our results motivate a re-examination of the debated pairing symmetry of this material when subjected to quantum confinement and lattice strain.


title: Trajectory Clustering Performance Evaluation: If we know the answer, it's not clustering

abstract: Advancements in Intelligent Traffic Systems (ITS) have made huge amounts of traffic data available through automatic data collection. A big part of this data is stored as trajectories of moving vehicles and road users. Automatic analysis of this data with minimal human supervision would both lower the costs and eliminate subjectivity of the analysis. Trajectory clustering is an unsupervised task. In this paper, we perform a comprehensive comparison of similarity measures, clustering algorithms and evaluation measures using trajectory data from seven intersections. We also propose a method to automatically generate trajectory reference clusters based on their origin and destination points to be used for label-based evaluation measures. Therefore, the entire procedure remains unsupervised both in clustering and evaluation levels. Finally, we use a combination of evaluation measures to find the top performing similarity measures and clustering algorithms for each intersection. The results show that there is no single combination of distance and clustering algorithm that is always among the top ten clustering setups.


abstract: Readability criteria, such as distance or neighborhood preservation, are often used to optimize node-link representations of graphs to enable the comprehension of the underlying data. With few exceptions, graph drawing algorithms typically optimize one such criterion, usually at the expense of others. We propose a layout approach, Multicriteria Scalable Graph Drawing via Stochastic Gradient Descent, $(SGD)^2$, that can handle multiple readability criteria. $(SGD)^2$ can optimize any criterion that can be described by a differentiable function. Our approach is flexible and can be used to optimize several criteria that have already been considered earlier (e.g., obtaining ideal edge lengths, stress, neighborhood preservation) as well as other criteria which have not yet been explicitly optimized in such fashion (e.g., node resolution, angular resolution, aspect ratio). The approach is scalable and can handle large graphs. A variation of the underlying approach can also be used to optimize many desirable properties in planar graphs, while maintaining planarity. Finally, we provide quantitative and qualitative evidence of the effectiveness of $(SGD)^2$: we analyze the interactions between criteria, measure the quality of layouts generated from $(SGD)^2$ as well as the runtime behavior, and analyze the impact of sample sizes. The source code is available on github and we also provide an interactive demo for small graphs.


title: Stability results assuming tameness, monster model and continuity of nonsplitting

abstract: Assuming the existence of a monster model, tameness and continuity of nonsplitting in an abstract elementary class (AEC), we extend known superstability results: let $\mu>LS({\bf K})$ be a regular stability cardinal and let $\chi$ be the local character of $\mu$-nonsplitting. The following holds: 1. When $\mu$-nonforking is restricted to $(\mu,\geq\chi)$-limit models ordered by universal extensions, it enjoys invariance, monotonicity, uniqueness, existence, extension and continuity. It also has local character $\chi$. This generalizes Vasey's result which assumed $\mu$-superstability to obtain same properties but with local character $\aleph_0$. 2. There is $\lambda\in[\mu,h(\mu))$ such that if ${\bf K}$ is stable in every cardinal between $\mu$ and $\lambda$, then ${\bf K}$ has $\mu$-symmetry while $\mu$-nonforking in (1) has symmetry. In this case (a) ${\bf K}$ has the uniqueness of $(\mu,\geq\chi)$-limit models: if $M_1,M_2$ are both $(\mu,\geq\chi)$-limit over some $M_0\in K_\mu$, then $M_1\cong_{M_0}M_2$; (b) any increasing chain of $\mu^+$-saturated models of length $\geq\chi$ has a $\mu^+$-saturated union. These generalize VanDieren-Vasey's result and remove the symmetry assumption in Boney-VanDieren and Vasey's result. Under $(<\mu)$-tameness, the conclusions of (1), (2)(a)(b) are equivalent to ${\bf K}$ having the $\chi$-local character of $\mu$-nonsplitting. Grossberg and Vasey gave eventual superstability criteria for tame AECs with a monster model. We remove the high cardinal threshold and reduce the cardinal jump between equivalent superstability criteria. We also add two new superstability criteria to the list: a weaker version of solvability and the boundedness of the $U$-rank.


title: FuseDream: Training-Free Text-to-Image Generation with Improved CLIP+GAN Space Optimization

abstract: Generating images from natural language instructions is an intriguing yet highly challenging task. We approach text-to-image generation by combining the power of the retrained CLIP representation with an off-the-shelf image generator (GANs), optimizing in the latent space of GAN to find images that achieve maximum CLIP score with the given input text. Compared to traditional methods that train generative models from text to image starting from scratch, the CLIP+GAN approach is training-free, zero shot and can be easily customized with different generators. However, optimizing CLIP score in the GAN space casts a highly challenging optimization problem and off-the-shelf optimizers such as Adam fail to yield satisfying results. In this work, we propose a FuseDream pipeline, which improves the CLIP+GAN approach with three key techniques: 1) an AugCLIP score which robustifies the CLIP objective by introducing random augmentation on image. 2) a novel initialization and over-parameterization strategy for optimization which allows us to efficiently navigate the non-convex landscape in GAN space. 3) a composed generation technique which, by leveraging a novel bi-level optimization formulation, can compose multiple images to extend the GAN space and overcome the data-bias. When promoted by different input text, FuseDream can generate high-quality images with varying objects, backgrounds, artistic styles, even novel counterfactual concepts that do not appear in the training data of the GAN we use. Quantitatively, the images generated by FuseDream yield top-level Inception score and FID score on MS COCO dataset, without additional architecture design or training. Our code is publicly available at \url{https://github.com/gnobitab/FuseDream}.


title: Dimension-Free Average Treatment Effect Inference with Deep Neural Networks

abstract: This paper investigates the estimation and inference of the average treatment effect (ATE) using deep neural networks (DNNs) in the potential outcomes framework. Under some regularity conditions, the observed response can be formulated as the response of a mean regression problem with both the confounding variables and the treatment indicator as the independent variables. Using such formulation, we investigate two methods for ATE estimation and inference based on the estimated mean regression function via DNN regression using a specific network architecture. We show that both DNN estimates of ATE are consistent with dimension-free consistency rates under some assumptions on the underlying true mean regression model. Our model assumptions accommodate the potentially complicated dependence structure of the observed response on the covariates, including latent factors and nonlinear interactions between the treatment indicator and confounding variables. We also establish the asymptotic normality of our estimators based on the idea of sample splitting, ensuring precise inference and uncertainty quantification. Simulation studies and real data application justify our theoretical findings and support our DNN estimation and inference methods.


title: Towards Interactive Reinforcement Learning with Intrinsic Feedback

abstract: Reinforcement learning (RL) and brain-computer interfaces (BCI) have experienced significant growth over the past decade. With rising interest in human-in-the-loop (HITL), incorporating human input with RL algorithms has given rise to the sub-field of interactive RL. Adjacently, the field of BCI has long been interested in extracting informative brain signals from neural activity for use in human-computer interactions. A key link between these fields lies in the interpretation of neural activity as feedback such that interactive RL approaches can be employed. We denote this new and emerging medium of feedback as intrinsic feedback. Despite intrinsic feedback's ability to be conveyed automatically and even unconsciously, proper exploration surrounding this key link has largely gone unaddressed by both communities. Thus, to help facilitate a deeper understanding and a more effective utilization, we provide a tutorial-style review covering the motivations, approaches, and open problems of intrinsic feedback and its foundational concepts.


title: Scheduling to Learn In An Unsupervised Online Streaming Model

abstract: An unsupervised online streaming model is considered where samples arrive in an online fashion over $T$ slots. There are $M$ classifiers, whose confusion matrices are unknown a priori. In each slot, at most one sample can be labeled by any classifier. The accuracy of a sample is a function of the set of labels obtained for it from various classifiers. The utility of a sample is a scalar multiple of its accuracy minus the response time (difference of the departure slot and the arrival slot), where the departure slot is also decided by the algorithm. Since each classifier can label at most one sample per slot, there is a tradeoff between obtaining a larger set of labels for a particular sample to improve its accuracy, and its response time. The problem of maximizing the sum of the utilities of all samples is considered, where learning the confusion matrices, sample-classifier matching assignment, and sample departure slot decisions depend on each other. The proposed algorithm first learns the confusion matrices, and then uses a greedy algorithm for sample-classifier matching. A sample departs once its incremental utility turns non-positive. We show that the competitive ratio of the proposed algorithm is $\frac{1}{2}-{\mathcal O}\left(\frac{\log T}{T}\right)$.


title: Artificial Intelligence-driven Image Analysis of Bacterial Cells and Biofilms

abstract: The current study explores an artificial intelligence framework for measuring the structural features from microscopy images of the bacterial biofilms. Desulfovibrio alaskensis G20 (DA-G20) grown on mild steel surfaces is used as a model for sulfate reducing bacteria that are implicated in microbiologically influenced corrosion problems. Our goal is to automate the process of extracting the geometrical properties of the DA-G20 cells from the scanning electron microscopy (SEM) images, which is otherwise a laborious and costly process. These geometric properties are a biofilm phenotype that allow us to understand how the biofilm structurally adapts to the surface properties of the underlying metals, which can lead to better corrosion prevention solutions. We adapt two deep learning models: (a) a deep convolutional neural network (DCNN) model to achieve semantic segmentation of the cells, (d) a mask region-convolutional neural network (Mask R-CNN) model to achieve instance segmentation of the cells. These models are then integrated with moment invariants approach to measure the geometric characteristics of the segmented cells. Our numerical studies confirm that the Mask-RCNN and DCNN methods are 227x and 70x faster respectively, compared to the traditional method of manual identification and measurement of the cell geometric properties by the domain experts.


title: Invariant Priors for Bayesian Quadrature

abstract: Bayesian quadrature (BQ) is a model-based numerical integration method that is able to increase sample efficiency by encoding and leveraging known structure of the integration task at hand. In this paper, we explore priors that encode invariance of the integrand under a set of bijective transformations in the input domain, in particular some unitary transformations, such as rotations, axis-flips, or point symmetries. We show initial results on superior performance in comparison to standard Bayesian quadrature on several synthetic and one real world application.


title: Fast Neural Representations for Direct Volume Rendering

abstract: Despite the potential of neural scene representations to effectively compress 3D scalar fields at high reconstruction quality, the computational complexity of the training and data reconstruction step using scene representation networks limits their use in practical applications. In this paper, we analyze whether scene representation networks can be modified to reduce these limitations and whether such architectures can also be used for temporal reconstruction tasks. We propose a novel design of scene representation networks using GPU tensor cores to integrate the reconstruction seamlessly into on-chip raytracing kernels, and compare the quality and performance of this network to alternative network- and non-network-based compression schemes. The results indicate competitive quality of our design at high compression rates, and significantly faster decoding times and lower memory consumption during data reconstruction. We investigate how density gradients can be computed using the network and show an extension where density, gradient and curvature are predicted jointly. As an alternative to spatial super-resolution approaches for time-varying fields, we propose a solution that builds upon latent-space interpolation to enable random access reconstruction at arbitrary granularity. We summarize our findings in the form of an assessment of the strengths and limitations of scene representation networks \changed{for compression domain volume rendering, and outline future research directions.


title: Challenging the disk instability model: I -- The case of YZ LMi

abstract: Observations of YZ LMi show enhanced emission along the stream trajectory beyond impact at disk rim during outbursts as well as when the quiescent disk is large. We investigated whether these features can be explained in terms of either gas stream overflow or penetration within the frameworks of the disk-instability (DIM) and the mass-transfer instability (MTIM) models of outbursting disks. Gas stream overflow is not possible because the vertical scaleheight of the stream is significantly lower than that of the outer disk and because there is no combination of parameters which enables stream overflow on a larger disk while preventing it on a smaller disk. Stream penetration requires the gas stream to be denser than the outer disk regions. This requirement cannot be met by a low-viscosity DIM disk because its density is significantly larger than that of the gas stream over the whole range of mass transfer rates where the thermal-viscous instability occurs. On the other hand, the high-viscosity MTIM disk has much lower densities which decrease with increasing radius, easily allowing for gas stream penetration during outbursts (when mass transfer rate and stream density increase) as well as in large quiescent disks. The observed features are not consistent with DIM, but can be plausibly explained by MTIM. These results suggest that the outbursts of YZ LMi are the response of a high-viscosity disk to bursts of enhanced mass transfer rate. In this case, the outburst decline timescale of (2-3) d implies a viscosity parameter in the range alpha=3-4.


title: On the Documentation of Refactoring Types

abstract: Commit messages are the atomic level of software documentation. They provide a natural language description of the code change and its purpose. Messages are critical for software maintenance and program comprehension. Unlike documenting feature updates and bug fixes, little is known about how developers document their refactoring activities. Developers can perform multiple refactoring operations, including moving methods, extracting classes, for various reasons. Yet, there is no systematic study that analyzes the extent to which the documentation of refactoring accurately describes the refactoring operations performed at the source code level. Therefore, this paper challenges the ability of refactoring documentation to adequately predict the refactoring types, performed at the commit level. Our analysis relies on the text mining of commit messages to extract the corresponding features that better represent each class. The extraction of text patterns, specific to each refactoring allows the design of a model that verifies the consistency of these patterns with their corresponding refactoring. Such verification process can be achieved via automatically predicting the method-level type of refactoring being applied, namely Extract Method, Inline Method, Move Method, Pull-up Method, Push-down Method, and Rename Method. We compared various classifiers, and a baseline keyword-based approach, in terms of their prediction performance, using a dataset of 5,004 commits. Our main findings show that the complexity of refactoring type prediction varies from one type to another. Rename method and Extract method were found to be the best documented refactoring activities, while Pull-up Method and Push-down Method were the hardest to be identified via textual descriptions. Such findings bring the attention of developers to the necessity of paying more attention to the documentation of these types.


title: LeapfrogLayers: A Trainable Framework for Effective Topological Sampling

abstract: We introduce LeapfrogLayers, an invertible neural network architecture that can be trained to efficiently sample the topology of a 2D $U(1)$ lattice gauge theory. We show an improvement in the integrated autocorrelation time of the topological charge when compared with traditional HMC, and look at how different quantities transform under our model. Our implementation is open source, and is publicly available on github at https://github.com/saforem2/l2hmc-qcd.


title: The Representation Jensen-R\'enyi Divergence

abstract: We introduce a divergence measure between data distributions based on operators in reproducing kernel Hilbert spaces defined by kernels. The empirical estimator of the divergence is computed using the eigenvalues of positive definite Gram matrices that are obtained by evaluating the kernel over pairs of data points. The new measure shares similar properties to Jensen-Shannon divergence. Convergence of the proposed estimators follows from concentration results based on the difference between the ordered spectrum of the Gram matrices and the integral operators associated with the population quantities. The proposed measure of divergence avoids the estimation of the probability distribution underlying the data. Numerical experiments involving comparing distributions and applications to sampling unbalanced data for classification show that the proposed divergence can achieve state of the art results.


title: Wearable Affective Memory Augmentation

abstract: Human memory prioritizes the storage and recall of information that is emotionally-arousing and/or important in a process known as value-directed memory. When experiencing a stream of information (e.g. conversation, book, lecture, etc.), the individual makes conscious and subconscious value assessments of the incoming information and uses this as a metric to determine what to remember. In order to improve automatic recall of memory, previous memory augmentation systems have sensed users' physiological state to determine which sensory media should be prioritized. Here, we propose to prioritize memories using the affective state of individuals that the user is interacting with. Thereby, the proposed wearable Affective Memory Augmentation system uses affective information from the user's social companions in order to facilitate value-directed memory.


title: Differentially Private Exploration in Reinforcement Learning with Linear Representation

abstract: This paper studies privacy-preserving exploration in Markov Decision Processes (MDPs) with linear representation. We first consider the setting of linear-mixture MDPs (Ayoub et al., 2020) (a.k.a.\ model-based setting) and provide an unified framework for analyzing joint and local differential private (DP) exploration. Through this framework, we prove a $\widetilde{O}(K^{3/4}/\sqrt{\epsilon})$ regret bound for $(\epsilon,\delta)$-local DP exploration and a $\widetilde{O}(\sqrt{K/\epsilon})$ regret bound for $(\epsilon,\delta)$-joint DP. We further study privacy-preserving exploration in linear MDPs (Jin et al., 2020) (a.k.a.\ model-free setting) where we provide a $\widetilde{O}\left(K^{\frac{3}{5}}/\epsilon^{\frac{2}{5}}\right)$ regret bound for $(\epsilon,\delta)$-joint DP, with a novel algorithm based on low-switching. Finally, we provide insights into the issues of designing local DP algorithms in this model-free setting.


title: HMC with Normalizing Flows

abstract: We propose using Normalizing Flows as a trainable kernel within the molecular dynamics update of Hamiltonian Monte Carlo (HMC). By learning (invertible) transformations that simplify our dynamics, we can outperform traditional methods at generating independent configurations. We show that, using a carefully constructed network architecture, our approach can be easily scaled to large lattice volumes with minimal retraining effort. The source code for our implementation is publicly available online at https://github.com/nftqcd/fthmc.


title: Improving accuracy and uncertainty quantification of deep learning based quantitative MRI using Monte Carlo dropout

abstract: Dropout is conventionally used during the training phase as regularization method and for quantifying uncertainty in deep learning. We propose to use dropout during training as well as inference steps, and average multiple predictions to improve the accuracy, while reducing and quantifying the uncertainty. The results are evaluated for fractional anisotropy (FA) and mean diffusivity (MD) maps which are obtained from only 3 direction scans. With our method, accuracy can be improved significantly compared to network outputs without dropout, especially when the training dataset is small. Moreover, confidence maps are generated which may aid in diagnosis of unseen pathology or artifacts.


title: String Bit Description of Antiperiodic Fermion Worldsheet Fields

abstract: We study a string bit Hamiltonian whose continuum limit describes anti-periodic (AP) anti-commuting worldsheet fields. We calculate the amplitude for transitions between an AP spin chain and a periodic (P) one in the continuum limit, M-> infinity where M is the bit number of either chain. We also numerically evaluate the corresponding amplitudes at increasing finite M to assess the convergence rate to the continuum. We then give the overlap equations for the transition AP+AP->AP, and numerically solve them for increasing $M$ values at a fixed value of x=K/M, where M is the bit number of the large chain and K is the bit number of one of the smaller chains. For this case, in contrast to the situation with an even number of AP chains, there is an obstacle to directly finding the continuum limit analytically. We suggest an indirect analytic approach to this problem: using the AP->P transition followed by a P->AP+AP transition, each of which has a relatively simple analytic continuum limit. We also show how bosonization of the fermion fields enables an analytic recursive evaluation of the AP$+$AP$\to$AP amplitudes.


title: InfoLM: A New Metric to Evaluate Summarization & Data2Text Generation

abstract: Assessing the quality of natural language generation systems through human annotation is very expensive. Additionally, human annotation campaigns are time-consuming and include non-reusable human labour. In practice, researchers rely on automatic metrics as a proxy of quality. In the last decade, many string-based metrics (e.g., BLEU) have been introduced. However, such metrics usually rely on exact matches and thus, do not robustly handle synonyms. In this paper, we introduce InfoLM a family of untrained metrics that can be viewed as a string-based metric that addresses the aforementioned flaws thanks to a pre-trained masked language model. This family of metrics also makes use of information measures allowing the adaptation of InfoLM to various evaluation criteria. Using direct assessment, we demonstrate that InfoLM achieves statistically significant improvement and over $10$ points of correlation gains in many configurations on both summarization and data2text generation.


title: The Art and Practice of Data Science Pipelines: A Comprehensive Study of Data Science Pipelines In Theory, In-The-Small, and In-The-Large

abstract: Increasingly larger number of software systems today are including data science components for descriptive, predictive, and prescriptive analytics. The collection of data science stages from acquisition, to cleaning/curation, to modeling, and so on are referred to as data science pipelines. To facilitate research and practice on data science pipelines, it is essential to understand their nature. What are the typical stages of a data science pipeline? How are they connected? Do the pipelines differ in the theoretical representations and that in the practice? Today we do not fully understand these architectural characteristics of data science pipelines. In this work, we present a three-pronged comprehensive study to answer this for the state-of-the-art, data science in-the-small, and data science in-the-large. Our study analyzes three datasets: a collection of 71 proposals for data science pipelines and related concepts in theory, a collection of over 105 implementations of curated data science pipelines from Kaggle competitions to understand data science in-the-small, and a collection of 21 mature data science projects from GitHub to understand data science in-the-large. Our study has led to three representations of data science pipelines that capture the essence of our subjects in theory, in-the-small, and in-the-large.


title: PLSUM: Generating PT-BR Wikipedia by Summarizing Multiple Websites

abstract: Wikipedia is an important free source of intelligible knowledge. Despite that, Brazilian Portuguese Wikipedia still lacks descriptions for many subjects. In an effort to expand the Brazilian Wikipedia, we contribute PLSum, a framework for generating wiki-like abstractive summaries from multiple descriptive websites. The framework has an extractive stage followed by an abstractive one. In particular, for the abstractive stage, we fine-tune and compare two recent variations of the Transformer neural network, PTT5, and Longformer. To fine-tune and evaluate the model, we created a dataset with thousands of examples, linking reference websites to Wikipedia. Our results show that it is possible to generate meaningful abstractive summaries from Brazilian Portuguese web content.


title: Online Search With Best-Price and Query-Based Predictions

abstract: In the online (time-series) search problem, a player is presented with a sequence of prices which are revealed in an online manner. In the standard definition of the problem, for each revealed price, the player must decide irrevocably whether to accept or reject it, without knowledge of future prices (other than an upper and a lower bound on their extreme values), and the objective is to minimize the competitive ratio, namely the worst-case ratio between the maximum price in the sequence and the one selected by the player. The problem formulates several applications of decision-making in the face of uncertainty on the revealed samples. Previous work on this problem has largely assumed extreme scenarios in which either the player has almost no information about the input, or the player is provided with some powerful, and error-free advice. In this work, we study learning-augmented algorithms, in which there is a potentially erroneous prediction concerning the input. Specifically, we consider two different settings: the setting in which the prediction is related to the maximum price in the sequence, as well as the setting in which the prediction is obtained as a response to a number of binary queries. For both settings, we provide tight, or near-tight upper and lower bounds on the worst-case performance of search algorithms as a function of the prediction error. We also provide experimental results on data obtained from stock exchange markets that confirm the theoretical analysis, and explain how our techniques can be applicable to other learning-augmented applications.


title: Types and Terms Translated: Unrestricted Resources in Encoding Functions as Processes (Extended Version)

abstract: Type-preserving translations are effective rigorous tools in the study of core programming calculi. In this paper, we develop a new typed translation that connects sequential and concurrent calculi; it is governed by type systems that control resource consumption. Our main contribution is the source language, a new resource $\lambda$-calculus with non-determinism and failures, dubbed \ulamf. In \ulamf, resources are split into linear and unrestricted; failures are explicit and arise from this distinction. We define a type system based on intersection types to control resources and fail-prone computation. The target language is \spi, an existing session-typed $\pi$-calculus that results from a Curry-Howard correspondence between linear logic and session types. Our typed translation subsumes our prior work; interestingly, it treats unrestricted resources in \lamrfailunres as client-server session behaviours in \spi.


title: On the Reliability of Multiple Systems Estimation for the Quantification of Modern Slavery

abstract: The quantification of modern slavery has received increased attention recently as organizations have come together to produce global estimates, where multiple systems estimation (MSE) is often used to this end. Echoing a long-standing controversy, disagreements have re-surfaced regarding the underlying MSE assumptions, the robustness of MSE methodology, and the accuracy of MSE estimates in this application. Our goal is to help address and move past these controversies. To do so, we review MSE, its assumptions, and commonly used models for modern slavery applications. We introduce all of the publicly available modern slavery datasets in the literature, providing a reproducible analysis and highlighting current issues. Specifically, we utilize an internal consistency approach that constructs subsets of data for which ground truth is available, allowing us to evaluate the accuracy of MSE estimators. Next, we propose a characterization of the large sample bias of estimators as a function of misspecified assumptions. Then, we propose an alternative to traditional (e.g., bootstrap-based) assessments of reliability, which allows us to visualize trajectories of MSE estimates to illustrate the robustness of estimates. Finally, our complementary analyses are used to provide guidance regarding the application and reliability of MSE methodology.


title: Smooth rigidity for codimension one Anosov flows

abstract: We introduce the matching functions technique in the setting of Anosov flows. Then we observe that simple periodic cycle functionals (also known as temporal distance functions) provide a source of matching functions for conjugate Anosov flows. For conservative codimension one Anosov flows $\varphi^t\colon M\to M$, $\dim M\ge 4$, these simple periodic cycle functionals are $C^1$ regular and, hence, can be used to improve regularity of the conjugacy. Specifically, we prove that a continuous conjugacy must, in fact, be a $C^1$ diffeomorphism for an open and dense set of codimension one conservative Anosov flows.


title: Hopfian and co-Hopfian modules over Artinian rings

abstract: An $R$-module $M$ is Hopfian (co-Hopfian) if any epic (monic) endomorphism of $M$ is an automorphism. If $R$ is commutative Noetherian, we characterize the co-Hopfian injective $R$-modules, and the Hopfian injectives in the case that $R$ is also reduced. For a commutative Artinian principal ideal ring, we show that $M$ is Hopfian (co-Hopfian) if and only if $M$ is finitely generated if and only if its injective envelope $E(M)$ is Hopfian (co-Hopfian) if and only if $E(M)$ is finitely generated. We identify the obstacle to generalizing this result to arbitrary Artinian principal ideal rings.


title: Quantum speed limit and stability of coherent states in quantum gravity

abstract: Utilizing the program of expectation values in coherent states and its recently developed algorithmic tools, this letter investigates the dynamical properties of cosmological coherent states for Loop Quantum Gravity. To this end, the Quantum Speed Limit is adapted to Quantum Gravity, yielding necessary consistency checks for any proposal of stable families of states. To showcase the strength of the developed tools, they are applied to a prominent model: the Euclidean part of the quantum scalar constraint. We report the variance of this constraint evaluated on a family of coherent states showing that, for short times, this family passes the Quantum Speed Limit test, allowing the transition from one coherent state to another one.


title: What Not to Test (for Cyber-Physical Systems)

abstract: For simulation-based systems, finding a set of test cases with the least cost by exploring multiple goals is a complex task. Domain-specific optimization goals (e.g. maximize output variance) are useful for guiding the rapid selection of test cases via mutation. But evaluating the selected test cases via mutation (that can distinguish the current program from the mutated systems) is a different goal to domain-specific optimizations. While the optimization goals can be used to guide the mutation analysis, that guidance should be viewed as a weak indicator since it can hurt the mutation effectiveness goals by focusing too much on the optimization goals. Based on the above, this paper proposes DoLesS (Domination with Least Squares Approximation) that selects the minimal and effective test cases by averaging over a coarse-grained grid of the information gained from multiple optimizations goals. DoLesS applies an inverted least squares approximation approach to find a minimal set of tests that can distinguish better from worse parts of the optimization goals. When tested on multiple simulation-based systems, DoLesS performs as well or even better as the prior state-of-the-art, while running 80-360 times faster on average (seconds instead of hours).


title: Bottomonium spectral widths at nonzero temperature using maximum likelihood

abstract: We present progress results from the Fastsum collaboration's programme to determine the spectrum of the bottomonium system as a function of temperature using a variety of approaches. In these proceedings, the Maximum Likelihood approach is used with an Ansatz comprising of a Gaussian spectral function for the ground state. Fastsum anisotropic lattices with 2+1 dynamical quark flavours were used with temperatures ranging from 47 to 375 MeV.


title: Experimental study of submerged liquid metal jet in transverse magnetic field

abstract: A liquid metal flow in the form of a submerged round jet entering a square duct in the presence of a transverse magnetic field is studied experimentally. A range of high Reynolds and Hartmann numbers is considered. Flow velocity is measured using electric potential difference probes. A detailed study of the flow in the duct's cross-section about seven jet's diameters downstream of the inlet reveals the dynamics, which is unsteady and dominated by high-amplitude fluctuations resulting from instability of the jet. The flow structure and fluctuation properties are largely determined by the value of the Stuart number N. At moderate N, the mean velocity profile retains a central jet with three-dimensional perturbations increasingly suppressed by the magnetic field as N grows. At higher values of N, the flow becomes quasi-two-dimensional and acquires the form of an asymmetric macrovortex, with high-amplitude velocity fluctuations reemerging.


title: Is RobustBench/AutoAttack a suitable Benchmark for Adversarial Robustness?

abstract: Recently, RobustBench (Croce et al. 2020) has become a widely recognized benchmark for the adversarial robustness of image classification networks. In its most commonly reported sub-task, RobustBench evaluates and ranks the adversarial robustness of trained neural networks on CIFAR10 under AutoAttack (Croce and Hein 2020b) with l-inf perturbations limited to eps = 8/255. With leading scores of the currently best performing models of around 60% of the baseline, it is fair to characterize this benchmark to be quite challenging. Despite its general acceptance in recent literature, we aim to foster discussion about the suitability of RobustBench as a key indicator for robustness which could be generalized to practical applications. Our line of argumentation against this is two-fold and supported by excessive experiments presented in this paper: We argue that I) the alternation of data by AutoAttack with l-inf, eps = 8/255 is unrealistically strong, resulting in close to perfect detection rates of adversarial samples even by simple detection algorithms and human observers. We also show that other attack methods are much harder to detect while achieving similar success rates. II) That results on low-resolution data sets like CIFAR10 do not generalize well to higher resolution images as gradient-based attacks appear to become even more detectable with increasing resolutions.


title: The conservative lock-in range for PLL with lead-lag filter and triangular phase detector characteristic

abstract: In the present work, a second-order PLL with lead-lag loop filter and triangular phase detector characteristic is analysed. An exact value of the conservative lock-in range is obtained for the considered model. The solution is based on analytical integration of the considered model on the linear segments.


title: Neurosymbolic Systems of Perception & Cognition: The Role of Attention

abstract: A cognitive architecture aimed at cumulative learning must provide the necessary information and control structures to allow agents to learn incrementally and autonomously from their experience. This involves managing an agent's goals as well as continuously relating sensory information to these in its perception-cognition information stack. The more varied the environment of a learning agent is, the more general and flexible must be these mechanisms to handle a wider variety of relevant patterns, tasks, and goal structures. While many researchers agree that information at different levels of abstraction likely differs in its makeup and structure and processing mechanisms, agreement on the particulars of such differences is not generally shared in the research community. A binary processing architecture (often referred to as System-1 and System-2) has been proposed as a model of cognitive processing for low- and high-level information, respectively. We posit that cognition is not binary in this way and that knowledge at any level of abstraction involves what we refer to as neurosymbolic information, meaning that data at both high and low levels must contain both symbolic and subsymbolic information. Further, we argue that the main differentiating factor between the processing of high and low levels of data abstraction can be largely attributed to the nature of the involved attention mechanisms. We describe the key arguments behind this view and review relevant evidence from the literature.


title: The Gardner problem and cycle slipping bifurcation for type 2 phase-locked loops

abstract: In the present work, a second-order type 2 PLL with a piecewise-linear phase detector characteristic is analysed. An exact solution to the Gardner problem on the lock-in range is obtained for the considered model. The solution is based on a study of cycle slipping bifurcation and improves well-known engineering estimates.


title: Optimal linear optical discrimination of Bell-like states

abstract: Quantum information processing using linear optics is challenging due to the limited set of deterministic operations achievable without using complicated resource-intensive methods. While techniques such as the use of ancillary photons can enhance the information processing capabilities of linear optical systems they are technologically demanding. Therefore, determining the constraints posed by linear optics and optimizing linear optical operations for specific tasks under those constraints, without the use of ancillas, can facilitate their potential implementation. Here, we consider the task of unambiguously discriminating between Bell-like states without the use of ancillary photons. This is a basic problem relevant in diverse settings, for example, in the measurement of the output of an entangling quantum circuit or for entanglement swapping at a quantum repeater station. While it is known that exact Bell states of two qubits can be discriminated with an optimal success probability of 50% we find, surprisingly, that for Bell-like states the optimal probability can be only 25%. We analyze a set of Bell-like states in terms of their distinguishability, entanglement as measured by concurrence, and parameters of the beam-splitter network used for unambiguous discrimination. Further, we provide the linear optical configuration comprised of single photon detectors and beam splitters with input state-dependent parameters that achieves optimal discrimination in the Bell-like case.


title: Test of streaming and triggered readout schemes for the TPEX Lead Tungstate Calorimeter

abstract: Tests of a prototype lead tungstate calorimeter were made over two weeks at the end of September, 2019, at the DESY II Test Beam Facility in Hamburg, Germany. The purpose of these tests was to gain experience with the construction, operation, and performance of a simple lead tungstate calorimeter, and also to compare a traditional triggered readout scheme with a streaming readout approach. These tests are important for the proposed Two-Photon Exchange experiment at DESY and for work towards a future electromagnetic calorimeter that could be used in an Electron-Ion Collider detector. Details on the comparison of streaming and triggered readout schemes are presented here.


title: Universal minimal cost of coherent biochemical oscillations

abstract: Biochemical clocks are essential for virtually all living systems. A biochemical clock that is isolated from an external periodic signal and subjected to fluctuations can oscillate coherently only for a finite number of oscillations. Furthermore, such an autonomous clock can oscillate only if it consumes free energy. What is the minimum amount of free energy consumption required for a certain number of coherent oscillations? We conjecture a universal bound that answers this question. A system that oscillates coherently for $\mathcal{N}$ oscillations has a minimal free energy cost per oscillation of $4\pi^2\mathcal{N} k_B T$. Our bound is valid for general finite Markov processes, is conjectured based on extensive numerical evidence, is illustrated with numerical simulations of a known model for a biochemical oscillator, and applies to existing experimental data.


title: Isomeric trees and the order of Runge--Kutta methods

abstract: The conditions for a Runge--Kutta method to be of order $p$ with $p\ge 5$ for a scalar non-autonomous problem are a proper subset of the order conditions for a vector problem. Nevertheless, Runge--Kutta methods that were derived historically only for scalar problems happened to be of the same order for vector problems. We relate the order conditions for scalar problems to factorisations of the Runge--Kutta trees into "atomic stumps" and enumerate those conditions up to $p=20$. Using a special search procedure over unsatisfied order conditions, new Runge--Kutta methods of "ambiguous orders" five and six are constructed. These are used to verify the validity of the results.


title: Probabilistic Tracking with Deep Factors

abstract: In many applications of computer vision it is important to accurately estimate the trajectory of an object over time by fusing data from a number of sources, of which 2D and 3D imagery is only one. In this paper, we show how to use a deep feature encoding in conjunction with generative densities over the features in a factor-graph based, probabilistic tracking framework. We present a likelihood model that combines a learned feature encoder with generative densities over them, both trained in a supervised manner. We also experiment with directly inferring probability through the use of image classification models that feed into the likelihood formulation. These models are used to implement deep factors that are added to the factor graph to complement other factors that represent domain-specific knowledge such as motion models and/or other prior information. Factors are then optimized together in a non-linear least-squares tracking framework that takes the form of an Extended Kalman Smoother with a Gaussian prior. A key feature of our likelihood model is that it leverages the Lie group properties of the tracked target's pose to apply the feature encoding on an image patch, extracted through a differentiable warp function inspired by spatial transformer networks. To illustrate the proposed approach we evaluate it on a challenging social insect behavior dataset, and show that using deep features does outperform these earlier linear appearance models used in this setting.


title: Recovering H\"older smooth functions from noisy modulo samples

abstract: In signal processing, several applications involve the recovery of a function given noisy modulo samples. The setting considered in this paper is that the samples corrupted by an additive Gaussian noise are wrapped due to the modulo operation. Typical examples of this problem arise in phase unwrapping problems or in the context of self-reset analog to digital converters. We consider a fixed design setting where the modulo samples are given on a regular grid. Then, a three stage recovery strategy is proposed to recover the ground truth signal up to a global integer shift. The first stage denoises the modulo samples by using local polynomial estimators. In the second stage, an unwrapping algorithm is applied to the denoised modulo samples on the grid. Finally, a spline based quasi-interpolant operator is used to yield an estimate of the ground truth function up to a global integer shift. For a function in H\"older class, uniform error rates are given for recovery performance with high probability. This extends recent results obtained by Fanuel and Tyagi for Lipschitz smooth functions wherein $k$NN regression was used in the denoising step.


title: Robust changepoint detection in the variability of multivariate functional data

abstract: We consider the problem of robustly detecting changepoints in the variability of a sequence of independent multivariate functions. We develop a novel changepoint procedure, called the functional Kruskal--Wallis for covariance (FKWC) changepoint procedure, based on rank statistics and multivariate functional data depth. The FKWC changepoint procedure allows the user to test for at most one changepoint (AMOC) or an epidemic period, or to estimate the number and locations of an unknown amount of changepoints in the data. We show that when the ``signal-to-noise'' ratio is bounded below, the changepoint estimates produced by the FKWC procedure attain the minimax localization rate for detecting general changes in distribution in the univariate setting (Theorem 1). We also provide the behavior of the proposed test statistics for the AMOC and epidemic setting under the null hypothesis (Theorem 2) and, as a simple consequence of our main result, these tests are consistent (Corollary 1). In simulation, we show that our method is particularly robust when compared to similar changepoint methods. We present an application of the FKWC procedure to intraday asset returns and f-MRI scans. As a by-product of Theorem 1, we provide a concentration result for integrated functional depth functions (Lemma 2), which may be of general interest.


title: Random Majorana Constellations

abstract: Even the most classical states are still governed by quantum theory. A fantastic array of physical systems can be described by their Majorana constellations of points on the surface of a sphere, where concentrated constellations and highly symmetric distributions correspond to the least and most quantum states, respectively. If these points are chosen randomly, how quantum will the resultant state be, on average? We explore this simple conceptual question in detail, investigating the quantum properties of the resulting random states. We find classical states to be far from the norm, even in the large-number-of-particles limit, where classical intuition often replaces quantum properties, making random Majorana constellations peculiar, intriguing, and useful.


title: Single-lens mass measurement in the high-magnification microlensing event Gaia19bld located in the Galactic disc

abstract: We present the photometric analysis of Gaia19bld, a high-magnification ($A\approx60$) microlensing event located in the southern Galactic plane, which exhibited finite source and microlensing parallax effects. Due to a prompt detection by the Gaia satellite and the very high brightness of $I = 9.05~$mag at the peak, it was possible to collect a complete and unique set of multi-channel follow-up observations, which allowed us to determine all parameters vital for the characterisation of the lens and the source in the microlensing event. Gaia19bld was discovered by the Gaia satellite and was subsequently intensively followed up with a network of ground-based observatories and the Spitzer Space Telescope. We collected multiple high-resolution spectra with Very Large Telescope (VLT)/X-Shooter to characterise the source star. The event was also observed with VLT Interferometer (VLTI)/PIONIER during the peak. Here we focus on the photometric observations and model the light curve composed of data from Gaia, Spitzer, and multiple optical, ground-based observatories. We find the best-fitting solution with parallax and finite source effects. We derived the limit on the luminosity of the lens based on the blended light model and spectroscopic distance. We compute the mass of the lens to be $1.13 \pm 0.03~M_{\odot}$ and derive its distance to be $5.52^{+0.35}_{-0.64}~\mathrm{kpc}$. The lens is likely a main sequence star, however its true nature has yet to be verified by future high-resolution observations. Our results are consistent with interferometric measurements of the angular Einstein radius, emphasising that interferometry can be a new channel for determining the masses of objects that would otherwise remain undetectable, including stellar-mass black holes.