title: Forward doubly-virtual Compton scattering off the nucleon in chiral perturbation theory: II. Spin polarizabilities and moments of polarized structure functions

abstract: We examine the polarized doubly-virtual Compton scattering (VVCS) off the nucleon using chiral perturbation theory ($\chi$PT). The polarized VVCS contains a wealth of information on the spin structure of the nucleon which is relevant to the calculation of the two-photon-exchange effects in atomic spectroscopy and electron scattering. We report on a complete next-to-leading-order (NLO) calculation of the polarized VVCS amplitudes $S_1(\nu, Q^2)$ and $S_2(\nu, Q^2)$, and the corresponding polarized spin structure functions $g_1(x, Q^2)$ and $g_2(x,Q^2)$. Our results for the moments of polarized structure functions, partially related to different spin polarizabilities, are compared to other theoretical predictions and "data-driven" evaluations, as well as to the recent Jefferson Lab measurements. By expanding the results in powers of the inverse nucleon mass, we reproduce the known "heavy-baryon" expressions. This serves as a check of our calculation, as well as demonstrates the differences between the manifestly Lorentz-invariant baryon $\chi$PT (B$\chi$PT) and heavy-baryon (HB$\chi$PT) frameworks.

categories: hep-ph hep-lat nucl-ex nucl-th physics.atom-ph


title: A New Multi-Wavelength Census of Blazars

abstract: Context:Blazars are the rarest and most powerful active galactic nuclei, playing a crucial and growing role in today multi-frequency and multi-messenger astrophysics. Current blazar catalogs, however, are incomplete and particularly depleted at low Galactic latitudes. Aims: We aim at augmenting the current blazar census to build a catalog of blazar candidates with homogeneous sky coverage that can provide candidate counterparts to unassociated gamma-ray sources, sources of high-energy neutrino emission, and ultra-high energy cosmic rays. Methods: Starting from the ALMA Calibrator Catalog we built a catalog of 1580 blazar candidates (ALMA Blazar Candidates, ABC) for which we collect multi-wavelength information. We also compared ABC sources with existing blazar catalogs. Results: The ABC catalogue fills the lack of low Galactic latitude sources in current blazar catalogues. ABC sources are significantly dimmer than known blazars in Gaia g band, and they appear bluer in SDSS and WISE colors. The majority of ABC sources (~ 90%) have optical spectra that classify them as QSO, while the remaining sources resulted galactic objects. ABC sources are similar in X-rays to known blazar, while in gamma-rays they are on average dimmer and softer, indicating a significant contribution of FSRQ sources. Making use of WISE colours, we classified 715 ABC sources as candidate gamma-ray blazar of different classes. Conclusions: We built a new catalogue of 1580 candidate blazars with a rich multi-wavelength data-set, filling the lack of low Galactic latitude sources in current blazar catalogues. This will be particularly important to identify the source population of high energy neutrinos or ultra-high energy cosmic rays. The data collected by the upcoming LSST surveys will provide a key tool to investigate the possible blazar nature of these sources.


title: Peeling off foregrounds with the constrained moment ILC method to unveil primordial CMB $B$-modes

abstract: Galactic foregrounds are the main obstacle to observations of the cosmic microwave background (CMB) $B$-mode polarization. In addition to obscuring the inflationary $B$-mode signal by several orders of magnitude, Galactic foregrounds have non-trivial spectral signatures that are partially unknown and distorted by averaging effects along the line-of-sight, within the pixel/beam window, and by various analysis choices (e.g., spherical harmonic transforms and filters). Statistical moment expansion methods provide a powerful tool for modeling the effective Galactic foreground emission resulting from these averaging effects in CMB observations, while blind component separation treatments can handle unknown foregrounds. In this work, we combine these two approaches to develop a new semi-blind component separation method at the intersection of parametric and blind methods, called constrained moment ILC (cMILC). This method adds several constraints to the standard ILC method to de-project the main statistical moments of the Galactic foreground emission. Applications to maps are performed in needlet space and when compared to the NILC method, this helps significantly reducing residual foreground contamination (bias, variance, and skewness) in the reconstructed CMB $B$-mode map, power spectrum, and tensor-to-scalar ratio. We consider sky-simulations for experimental settings similar to those of LiteBIRD and PICO, illustrating which trade-offs between residual foreground biases and degradation of the constraint on $r$ can be expected within the new cMILC framework. We also outline several directions that require more work in preparation for the coming analysis challenges.


title: Anomalies, Black strings and the charged Cardy formula

abstract: We derive the general anomaly polynomial for a class of two-dimensional CFTs arising as twisted compactifications of a higher-dimensional theory on compact manifolds $\mathcal{M}_d$, including the contribution of the isometries of $\mathcal{M}_d$. We then use the result to perform a counting of microstates for electrically charged and rotating supersymmetric black strings in AdS$_5\times S^5$ and AdS$_7\times S^4$ with horizon topology BTZ$ \ltimes S^2$ and BTZ$ \ltimes S^2 \times \Sigma_\mathfrak{g}$, respectively, where $\Sigma_\mathfrak{g}$ is a Riemann surface. We explicitly construct the latter class of solutions by uplifting a class of four-dimensional rotating black holes. We provide a microscopic explanation of the entropy of such black holes by using a charged version of the Cardy formula.


title: The Newman-Penrose Map and the Classical Double Copy

abstract: Gauge-gravity duality is arguably our best hope for understanding quantum gravity. Considerable progress has been made in relating scattering amplitudes in certain gravity theories to those in gauge theories---a correspondence dubbed the "double copy". Recently, double copies have also been realized in a classical setting, as maps between exact solutions of gauge theories and gravity. We present here a novel map between a certain class of real, exact solutions of Einstein's equations and self-dual solutions of the flat-space vacuum Maxwell equations. This map, which we call the "Newman-Penrose map", is well-defined even for non-vacuum, non-stationary spacetimes, providing a systematic framework for exploring gravity solutions in the context of the double copy that have not been previously studied in this setting. To illustrate this, we present here the Newman-Penrose map for the Schwarzschild and Kerr black holes, and Kinnersley's photon rocket.


title: Collisional picture of quantum optics with giant emitters

abstract: The effective description of the weak interaction between an emitter and a bosonic field as a sequence of two-body collisions provides a simple intuitive picture compared to traditional quantum optics methods as well as an effective calculation tool of the joint emitter-field dynamics. Here, this collisional approach is extended to many emitters (atoms or resonators), each generally interacting with the field at many coupling points ("giant" emitter). In the regime of negligible delays, the unitary describing each collision in particular features a contribution of a chiral origin resulting in an effective Hamiltonian. The picture is applied to derive a Lindblad master equation (ME) of a set of giant atoms coupled to a (generally chiral) waveguide field in an arbitrary white-noise Gaussian state, which condenses into a single equation and extends a variety of quantum optics and waveguide-QED MEs. The effective Hamiltonian and jump operators corresponding to a selected photodetection scheme are also worked out.


title: The Mass-Size Relation and the Constancy of GMC Surface Densities in the Milky Way

abstract: We use two existing molecular cloud catalogs derived from the same CO survey and two catalogs derived from local dust extinction surveys to investigate the nature of the GMC mass-size relation in the Galaxy. We find that the four surveys are well described by $M_{GMC} \sim R^2$ implying a constant mean surface density, $\Sigma_{GMC}$, for the cataloged clouds. However, the scaling coefficients and scatter differ significantly between the CO and extinction derived relations. We find that the additional scatter seen in the CO relations is due to a systematic variation in $\Sigma_{GMC}$ with Galactic radius that is unobservable in the local extinction data. We decompose this radial variation of $\Sigma_{GMC}$ into two components, a linear negative gradient with Galactic radius and a broad peak coincident with the molecular ring and superposed on the linear gradient. We show that the former may be due to a radial dependence of X$_{CO}$ on metallicity while the latter likely results from a combination of increased surface densities of individual GMCs and a systematic upward bias in the measurements of $\Sigma_{GMC}$ due to cloud blending in the molecular ring. We attribute the difference in scaling coefficients between the CO and extinction data to an underestimate of X$_{CO}$. We recalibrate the CO observations of nearby GMCs using extinction measurements to find that locally X$_{CO}$ $=$ 3.6$\pm$0.3 $\times$ 10$^{20}$ cm$^{-2}$ (K-km/s)$^{-1}$. We conclude that outside the molecular ring the GMC population of the Galaxy can be described to relatively good precision by a constant $\Sigma_{GMC}$ of 35 M$_\odot$ pc$^{-2}$.


title: Surface Operators in Superspace

abstract: We generalize the geometrical formulation of Wilson loops recently introduced in arXiv:2003.01729v2 to the description of Wilson Surfaces. For N=(2,0) theory in six dimensions, we provide an explicit derivation of BPS Wilson Surfaces with non-trivial coupling to scalars, together with their manifestly supersymmetric version. We derive explicit conditions which allow to classify these operators in terms of the number of preserved supercharges. We also discuss kappa-symmetry and prove that BPS conditions in six dimensions arise from kappa-symmetry invariance in eleven dimensions. Finally, we discuss super-Wilson Surfaces - and higher dimensional operators - as objects charged under global $p$-form (super)symmetries generated by tensorial supercurrents. To this end, the construction of conserved supercurrents in supermanifolds and of the corresponding conserved charges is developed in details.


title: The intrinsic SFRF and sSFRF of galaxies: comparing SDSS observation with IllustrisTNG simulation

abstract: The star formation rate function (SFRF) and specific star formation rate function (sSFRF) from the observation are impacted by the Eddington bias, due to the uncertainties on the estimated SFR. We develop a novel method to correct the Eddington bias and obtained the intrinsic SFRF and sSFRF from the Sloan Digital Sky Survey Data Release 7. The intrinsic SFRF is in good agreement with measurements from previous data in the literature that relied on UV SFRs but its high star-forming end is slightly lower than those IR and radio tracers. We demonstrate that the intrinsic sSFRF from SDSS has a bi-modal form with the one peak found at ${\rm sSFR \sim 10^{-9.7} yr^{-1}}$ representing the star-forming objects while the other peak is found at ${\rm sSFR \sim 10^{-12} yr^{-1}}$ representing the quenched population. Furthermore, we compare our observations with the predictions from the IllustrisTNG and Illustris simulations and show that the ``TNG'' model performs much better than its predecessor. However, we show that the simulated SFRF and cosmic star formation density (CSFRD) of TNG simulations are highly dependent on resolution, reflecting the limitations of the model and today state-of-the-art simulations. We demonstrate that the bi-modal, two peaked sSFRF implied by the SDSS observations does not appear in TNG regardless of the adopted box-size or resolution. This tension reflects the need for inclusion of an additional efficient quenching mechanism to the TNG model.


title: Effect of a strong magnetic field on gravity-mode period spacings in red giant stars

abstract: When a star evolves into a red giant, the enhanced coupling between core-based gravity modes and envelope-based pressure modes forms mixed modes, allowing its deep interior to be probed by asteroseismology. The ability to obtain information about stellar interiors is important for constraining theories of stellar structure and evolution, for which the origin of various discrepancies between prediction and observation are still under debate. Ongoing speculation surrounds the possibility that some red giant stars may harbour strong (dynamically significant) magnetic fields in their cores, but interpretation of the observational data remains controversial. In part, this is tied to shortfalls in our understanding of the effects of strong fields on the seismic properties of gravity modes, which lies beyond the regime of standard perturbative methods. Here we seek to investigate the effect of a strong magnetic field on the asymptotic period spacings of gravity modes. We use a Hamiltonian ray approach to measure the volume of phase space occupied by mode-forming rays, this being roughly proportional to the average density of modes (number of modes per unit frequency interval). A strong field appears to systematically increase this by about 10%, which predicts a ~10% smaller period spacing. Evidence of near integrability in the ray dynamics hints that the gravity-mode spectrum may still exhibit pseudo-regularities under a strong field.


title: Parton distributions and lattice QCD calculations: toward 3D structure

abstract: The strong force which binds hadrons is described by the theory of Quantum Chromodynamics (QCD). Determining the character and manifestations of QCD is one of the most important and challenging outstanding issues necessary for a comprehensive understanding of the structure of hadrons. Within the context of the QCD parton picture, the Parton Distribution Functions (PDFs) have been remarkably successful in describing a wide variety of processes. However, these PDFs have generally been confined to the description of collinear partons within the hadron. New experiments and facilities provide the opportunity to additionally explore the transverse structure of hadrons which is described by Generalized Parton Distributions (GPDs) and Transverse Momentum Dependent Parton Distribution Functions (TMD PDFs). In our previous review, we compared and contrasted the two main approaches used to determine the collinear PDFs: the first based on perturbative QCD factorization theorems, and the second based on lattice QCD calculations. In the present report, we provide an update of recent progress on the collinear PDFs, and also expand the scope to encompass the generalized PDFs (GPDs and TMD PDFs). We review the current state of the various calculations, and consider what new data might be available in the near future. We also examine how a shared effort can foster dialog between the PDF and Lattice QCD communities, and yield improvements for these generalized PDFs.


title: Axion Miniclusters Made Easy

abstract: We use a modified version of the Peak Patch excursion set formalism to compute the mass and size distribution of QCD axion miniclusters from a fully non-Gaussian initial density field obtained from numerical simulations of axion string decay. We find strong agreement with N-Body simulations at a significantly lower computational cost. We employ a spherical collapse model and provide fitting functions for the modified barrier in the radiation era. The halo mass function at $z=629$ has a power-law distribution $M^{-0.6}$ for masses within the range $10^{-15}\lesssim M\lesssim 10^{-10}M_{\odot}$, with all masses scaling as $(m_a/50\mu\mathrm{eV})^{-0.5}$. We construct merger trees to estimate the collapse redshift and concentration mass relation, $C(M)$, which is well described using analytical results from the initial power spectrum and linear growth. Using the calibrated analytic results to extrapolate to $z=0$, our method predicts a mean concentration $C\sim \mathcal{O}(\text{few})\times10^4$. The low computational cost of our method makes future investigation of the statistics of rare, dense miniclusters easy to achieve.


title: The Data Processing Pipeline for the MUSE Instrument

abstract: Processing of raw data from modern astronomical instruments is nowadays often carried out using dedicated software, so-called "pipelines" which are largely run in automated operation. In this paper we describe the data reduction pipeline of the Multi Unit Spectroscopic Explorer (MUSE) integral field spectrograph operated at ESO's Paranal observatory. This spectrograph is a complex machine: it records data of 1152 separate spatial elements on detectors in its 24 integral field units. Efficiently handling such data requires sophisticated software, a high degree of automation and parallelization. We describe the algorithms of all processing steps that operate on calibrations and science data in detail, and explain how the raw science data gets transformed into calibrated datacubes. We finally check the quality of selected procedures and output data products, and demonstrate that the pipeline provides datacubes ready for scientific analysis.


title: Casting a graph net to catch dark showers

abstract: Strongly interacting dark sectors predict novel LHC signatures such as semi-visible jets resulting from dark showers that contain both stable and unstable dark mesons. Distinguishing such semi-visible jets from large QCD backgrounds is difficult and constitutes an exciting challenge for jet classification. In this article we explore the potential of supervised deep neural networks to identify semi-visible jets. We show that dynamic graph convolutional neural networks operating on so-called particle clouds outperform convolutional neural networks analysing jet images as well as other neural networks based on Lorentz vectors. We investigate how the performance depends on the properties of the dark shower and discuss training on mixed samples as a strategy to reduce model dependence. By modifying an existing mono-jet analysis we show that LHC sensitivity to dark sectors can be enhanced by more than an order of magnitude by using the dynamic graph network as a dark shower tagger.


title: Forecasting Chemical Abundance Precision for Extragalactic Stellar Archaeology

abstract: Increasingly powerful and multiplexed spectroscopic facilities promise detailed chemical abundance patterns for millions of resolved stars in galaxies beyond the Milky Way (MW). Here, we employ the Cram\'er-Rao Lower Bound (CRLB) to forecast the precision to which stellar abundances for metal-poor, low-mass stars outside the MW can be measured for 41 current (e.g., Keck, MMT, VLT, DESI) and planned (e.g., MSE, JWST, ELTs) spectrograph configurations. We show that moderate resolution ($R\lesssim5000$) spectroscopy at blue-optical wavelengths ($\lambda\lesssim4500$ \AA) (i) enables the recovery of 2-4 times as many elements as red-optical spectroscopy ($5000\lesssim\lambda\lesssim10000$ \AA) at similar or higher resolutions ($R\sim 10000$) and (ii) can constrain the abundances of several neutron capture elements to $\lesssim$0.3 dex. We further show that high-resolution ($R\gtrsim 20000$), low S/N ($\sim$10 pixel$^{-1}$) spectra contain rich abundance information when modeled with full spectral fitting techniques. We demonstrate that JWST/NIRSpec and ELTs can recover (i) $\sim$10 and 30 elements, respectively, for metal-poor red giants throughout the Local Group and (ii) [Fe/H] and [$\alpha$/Fe] for resolved stars in galaxies out to several Mpc with modest integration times. We show that select literature abundances are within a factor of $\sim$2 (or better) of our CRLBs. We suggest that, like ETCs, CRLBs should be used when planning stellar spectroscopic observations. We include an open source python package, \texttt{Chem-I-Calc}, that allows users to compute CRLBs for spectrographs of their choosing.


title: The Pristine Inner Galaxy Survey (PIGS) II: Uncovering the most metal-poor populations in the inner Milky Way

abstract: Metal-poor stars are important tools for tracing the early history of the Milky Way, and for learning about the first generations of stars. Simulations suggest that the oldest metal-poor stars are to be found in the inner Galaxy. Typical bulge surveys, however, lack low metallicity ([Fe/H] < -1.0) stars because the inner Galaxy is predominantly metal-rich. The aim of the Pristine Inner Galaxy Survey (PIGS) is to study the metal-poor and very metal-poor (VMP, [Fe/H] < -2.0) stars in this region. In PIGS, metal-poor targets for spectroscopic follow-up are selected from metallicity-sensitive CaHK photometry from the CFHT. This work presents the ~250 deg^2 photometric survey as well as intermediate-resolution spectroscopic follow-up observations for ~8000 stars using AAOmega on the AAT. The spectra are analysed using two independent tools: ULySS with an empirical spectral library, and FERRE with a library of synthetic spectra. The comparison between the two methods enables a robust determination of the stellar parameters and their uncertainties. We present a sample of 1300 VMP stars -- the largest sample of VMP stars in the inner Galaxy to date. Additionally, our spectroscopic dataset includes ~1700 horizontal branch stars, which are useful metal-poor standard candles. We furthermore show that PIGS photometry selects VMP stars with unprecedented efficiency: 86%/80% (lower/higher extinction) of the best candidates satisfy [Fe/H] < -2.0, as do 80%/63% of a larger, less strictly selected sample. We discuss future applications of this unique dataset that will further our understanding of the chemical and dynamical evolution of the innermost regions of our Galaxy.


title: Sixfold fermion near the Fermi level in cubic PtBi2

abstract: We show that the cubic compound PtBi2, is a topological semimetal hosting a sixfold band touching point in close proximity to the Fermi level. Using angle-resolved photoemission spectroscopy, we map the bandstructure of the system, which is in good agreement with results from density functional theory. Further, by employing a low energy effective Hamiltonian valid close to the crossing point, we study the effect of a magnetic field on the sixfold fermion. The latter splits into a total of twenty Weyl cones for a Zeeman field oriented in the diagonal, [111] direction. Our results mark cubic PtBi2, as an ideal candidate to study the transport properties of gapless topological systems beyond Dirac and Weyl semimetals.


title: On the training dynamics of deep networks with $L_2$ regularization

abstract: We study the role of $L_2$ regularization in deep learning, and uncover simple relations between the performance of the model, the $L_2$ coefficient, the learning rate, and the number of training steps. These empirical relations hold when the network is overparameterized. They can be used to predict the optimal regularization parameter of a given model. In addition, based on these observations we propose a dynamical schedule for the regularization parameter that improves performance and speeds up training. We test these proposals in modern image classification settings. Finally, we show that these empirical relations can be understood theoretically in the context of infinitely wide networks. We derive the gradient flow dynamics of such networks, and compare the role of $L_2$ regularization in this context with that of linear models.


title: $SL(2,R)$ lattices as information processors

abstract: Black holes past their Page times should act as efficient scramblers and information mirrors. The information of the infalling bits are rapidly encoded by the old black hole in the Hawking quanta, but it should take time that is exponential in the Page time entropy to decode the interior. Motivated by the features of fragmentation instability of near-extremal black holes, we construct a simple phenomenological model of the black hole as a lattice of interacting nearly $AdS_2$ throats with gravitational hair charges propagating over the lattice. We study the microstate solutions and their response to shocks. The energy of the shocks are almost wholly absorbed by the total ADM mass of the $AdS_2$ throats, but the information of their locations and time-ordering come out in the hair oscillations, which decouple from the final microstate to which the full system quickly relaxes. We discuss the Hayden-Preskill protocol of decoding infalling information. We also construct generalizations of our model involving a lattice of $AdS_2$ throats networked via wormholes and their analogues in the form of tensor networks of SYK spin-states.


title: High-redshift Extreme Variability Quasars from Sloan Digital Sky Survey Multi-Epoch Spectroscopy

abstract: We perform a systematic search for high-redshift ($z >$ 1.5) extreme variability quasars (EVQs) using repeat spectra from the Sixteenth Data Release of Sloan Digital Sky Survey, which provides a baseline spanning up to $\sim$18 yrs in the observed frame. We compile a sample of 348 EVQs with a maximum continuum variability at rest frame 1450 Angstrom of more than 100% (i.e., $\delta$V $\equiv$ (Max$-$Min)/Mean $>$1). The EVQs show a range of emission line variability, including 23 where at least one line in our redshift range disappears below detectability, which can then be seen as analogous to low-redshift changing-look quasars (CLQs)". Importantly, spurious CLQs caused by SDSS problematic spectral flux calibration, e.g., fiber drop issue, have been rejected. The similar properties (e.g., continuum/line, difference-composite spectra and Eddington ratio) of normal EVQs and CLQs, implies that they are basically the same physical population with analogous intrinsic variability mechanisms, as a tail of a continuous distribution of normal quasar properties. In addition, we find no reliable evidence ($\lesssim$ 1$\sigma$) to support that the CLQs are a subset of EVQs with less efficient accretion. Finally, we also confirm the anti-breathing of C IV (i.e., line width increases as luminosity increases) in EVQs, and find that in addition to $\sim$ 0.4 dex systematic uncertainty in single-epoch C IV virial black hole mass estimates, an extra scatter of $\sim$ 0.3 dex will be introduced by extreme variability.


title: Quantum trimer models and topological SU(3) spin liquids on the kagome lattice

abstract: We construct and study quantum trimer models and resonating SU(3)-singlet models on the kagome lattice, which generalize quantum dimer models and the Resonating Valence Bond wavefunctions to a trimer and SU(3) setting. We demonstrate that these models carry a Z_3 symmetry which originates in the structure of trimers and the SU(3) representation theory, and which becomes the only symmetry under renormalization. Based on this, we construct simple and exact parent Hamiltonians for the model which exhibit a topological 9-fold degenerate ground space. A combination of analytical reasoning and numerical analysis reveals that the quantum order ultimately displayed by the model depends on the relative weight assigned to different types of trimers -- it can display either Z_3 topological order or form a symmetry-broken trimer crystal, and in addition possesses a point with an enhanced U(1) symmetry and critical behavior. Our results accordingly hold for the SU(3) model, where the two natural choices for trimer weights give rise to either a topological spin liquid or a system with symmetry-broken order, respectively. Our work thus demonstrates the suitability of resonating trimer and SU(3)-singlet ansatzes to model SU(3) topological spin liquids on the kagome lattice.


title: Universal Relation for Life-span Energy Consumption in Living Organisms: Insights for the origin of ageing

abstract: Metabolic energy consumption has long been thought to play a major role in the aging process ({\it 1}). Across species, a gram of tissue on average expends about the same amount of energy during life-span ({\it 2}). Energy restriction has also been shown that increases maximum life-span ({\it 3}) and retards age-associated changes ({\it 4}). However, there are significant exceptions to a universal energy consumption during life-span, mainly coming from the inter-class comparison ({\it 5, 6}). Here we present a unique relation for life-span energy consumption, valid for $\sim$300 species representing all classes of living organisms, from unicellular ones to the largest mammals. The relation has an average scatter of only 0.3 dex, with 95\% ($\rm 2-\sigma$) of the organisms having departures less than a factor of $\pi$ from the relation, despite the $\sim$20 orders of magnitude difference in body mass, reducing any possible inter-class variation in the relation to only a geometrical factor. This result can be interpreted as supporting evidence for the existence of an approximately constant total number $\rm N_r \sim 10^8$ of respiration cycles per lifetime for all organisms, effectively predetermining the extension of life by the basic energetics of respiration, being an incentive for future studies that investigate the relation of such constant $\rm N_r$ cycles per lifetime with the production rates of free radicals and oxidants, which may give definite constraints on the origin of ageing.


title: AdS$_3$ gravity and random CFT

abstract: We compute the path integral of three-dimensional gravity with negative cosmological constant on spaces which are topologically a torus times an interval. These are Euclidean wormholes, which smoothly interpolate between two asymptotically Euclidean AdS$_3$ regions with torus boundary. From our results we obtain the spectral correlations between BTZ black hole microstates near threshold, as well as extract the spectral form factor at fixed momentum, which has linear growth in time with small fluctuations around it. The low-energy limit of these correlations is precisely that of a double-scaled random matrix ensemble with Virasoro symmetry. Our findings suggest that if pure three-dimensional gravity has a holographic dual, then the dual is an ensemble which generalizes random matrix theory.


title: Open giant magnons suspended between dual giant gravitons in ${\cal N}=4$ SYM

abstract: We study classical solutions to the Nambu-Goto string on $AdS_2 \times S^1$ and $AdS_3 \times S^1$ corresponding to strings stretched between wrapped branes extending in $AdS$. The solutions are obtained by analytic continuation of giant magnon solutions, cut at the position of the branes. These solutions carry one or two $SO(2,4)$ charges and a single $SO(6)$ charge. We compute their energies and show their relation to $\frac{1}{2}$-BPS geometries. Their relevance to the $SL(2)$ sector of $\mathcal{N}=4$ SYM is also discussed.


title: Masses of Doubly Heavy Tetraquarks with Error Bars

abstract: In the heavy-quark limit, the two heavy quarks in a doubly heavy baryon or a doubly heavy tetraquark are bound by their color-Coulomb potential into a compact diquark. The doubly heavy hadrons are related by the approximate heavy-quark--diquark symmetry of QCD to the heavy hadrons obtained by replacing the heavy diquark by a heavy antiquark. Effective field theories can be used to expand the masses of singly heavy hadrons and doubly heavy hadrons in inverse powers of the heavy quark masses. The coefficients in the expansions for doubly heavy tetraquarks can be determined from those for heavy mesons, heavy baryons, and doubly heavy baryons using heavy-quark--diquark symmetry. We predict the masses of the ground-state doubly heavy tetraquarks with error bars using as inputs the masses of heavy mesons and heavy baryons measured in experiments and the masses of doubly heavy baryons calculated using lattice QCD. The only doubly heavy tetraquarks predicted to be stable with respect to strong decays are $bb$ tetraquarks with light flavor $\bar u \bar d$, $\bar s \bar u$ and $\bar s \bar d$.


title: Measuring the surface mass density ellipticity of redMaPPer galaxy clusters using weak-lensing

abstract: In this work we study the shape of the projected surface mass density distribution of galaxy clusters using weak-lensing stacking techniques. In particular, we constrain the average aligned component of the projected ellipticity, $\epsilon$, for a sample of redMaPPer clusters ($0.1 \leq z < 0.4$). We consider six different proxies for the cluster orientation and measure $\epsilon$ for three ranges of projected distances from the cluster centres. The mass distribution in the inner region (up to $700\,$kpc) is better traced by the cluster galaxies with a higher membership probability, while the outer region (from $700\,$kpc up to $5\,$Mpc) is better traced by the inclusion of less probable galaxy cluster members. The fitted ellipticity in the inner region is $\epsilon = 0.21 \pm 0.04$, in agreement with previous estimates. We also study the relation between $\epsilon$ and the cluster mean redshift and richness. By splitting the sample in two redshift ranges according to the median redshift, we obtain larger $\epsilon$ values for clusters at higher redshifts, consistent with the expectation from simulations. In addition, we obtain higher ellipticity values in the outer region of clusters at low redshifts. We discuss several systematic effects that might affect the measured lensing ellipticities and their relation to the derived ellipticity of the mass distribution.


title: Multifractality in the general cosmological solution of Einstein's equations

abstract: We demonstrate the scale invariance of the vacuum Bianchi type IX equations and use this to argue for the possibility of multifractal turbulence as a realisation of the suggestion by Belinski that there will be a fragmentation of local regions of inhomogeneous Mixmaster chaos on approach to an initial inhomogeneous cosmological singularity. Differences between the gravitational and hydrodynamical situations are outlined. Various potential obstacles to this picture of gravitational turbulence are discussed.


title: Simulating exotic phases of matter with bond-directed interactions with arrays of Majorana-Cooper pair boxes

abstract: It is suggested that networks of Majorana-Cooper pair boxes connected by metallic nanowires can simulate various exotic states of matter. In this simulations Majorana-Cooper boxes play the role of effective spins S=1/2 and the metallic connections generate the Kondo screening and the Ruderman-Kittel-Kasuya-Yosida (RKKY) interaction. Depending on what prevails - whether it is the Kondo effect or the RKKY exchange, one will have either an effective spin model or a Kondo lattice. The list of exotic stets includes the famous hexagonal Kitaev model, a generalization of this model for a Kondo lattice and various spin models with three-spin interactions. A special emphasize is made on the discussion of the Kondo lattice scaenario.


title: Triggerflow: Trigger-based Orchestration of Serverless Workflows

abstract: As more applications are being moved to the Cloud thanks to serverless computing, it is increasingly necessary to support native life cycle execution of those applications in the data center. But existing systems either focus on short-running workflows (like IBM Composer or Amazon Express Workflows) or impose considerable overheads for synchronizing massively parallel jobs (Azure Durable Functions, Amazon Step Functions, Google Cloud Composer). None of them are open systems enabling extensible interception and optimization of custom workflows. We present Triggerflow: an extensible Trigger-based Orchestration architecture for serverless workflows built on top of Knative Eventing and Kubernetes technologies. We demonstrate that Triggerflow is a novel serverless building block capable of constructing different reactive schedulers (State Machines, Directed Acyclic Graphs, Workflow as code). We also validate that it can support high-volume event processing workloads, auto-scale on demand and transparently optimize scientific workflows.


title: Variational Bayesian Monte Carlo with Noisy Likelihoods

abstract: Variational Bayesian Monte Carlo (VBMC) is a recently introduced framework that uses Gaussian process surrogates to perform approximate Bayesian inference in models with black-box, non-cheap likelihoods. In this work, we extend VBMC to deal with noisy log-likelihood evaluations, such as those arising from simulation-based models. We introduce new `global' acquisition functions, such as expected information gain (EIG) and variational interquantile range (VIQR), which are robust to noise and can be efficiently evaluated within the VBMC setting. In a novel, challenging, noisy-inference benchmark comprising of a variety of models with real datasets from computational and cognitive neuroscience, VBMC+VIQR achieves state-of-the-art performance in recovering the ground-truth posteriors and model evidence. In particular, our method vastly outperforms `local' acquisition functions and other surrogate-based inference methods while keeping a small algorithmic cost. Our benchmark corroborates VBMC as a general-purpose technique for sample-efficient black-box Bayesian inference also with noisy models.


title: Multiscale Deep Equilibrium Models

abstract: We propose a new class of implicit networks, the multiscale deep equilibrium model (MDEQ), suited to large-scale and highly hierarchical pattern recognition domains. An MDEQ directly solves for and backpropagates through the equilibrium points of multiple feature resolutions simultaneously, using implicit differentiation to avoid storing intermediate states (and thus requiring only $O(1)$ memory consumption). These simultaneously-learned multi-resolution features allow us to train a single model on a diverse set of tasks and loss functions, such as using a single MDEQ to perform both image classification and semantic segmentation. We illustrate the effectiveness of this approach on two large-scale vision tasks: ImageNet classification and semantic segmentation on high-resolution images from the Cityscapes dataset. In both settings, MDEQs are able to match or exceed the performance of recent competitive computer vision models: the first time such performance and scale have been achieved by an implicit deep learning approach. The code and pre-trained models are at https://github.com/locuslab/mdeq .


title: Opacity Driven Convection and Variability in Accretion Disks around Supermassive Black Holes

abstract: We study the structure of accretion disks around supermassive black holes in the radial range $30\sim 100$ gravitational radii, using a three dimensional radiation magneto-hydrodynamic simulation. For typical conditions in this region of Active Galactic Nuclei (AGN), the Rosseland mean opacity is expected to be larger than the electron scattering value. We show that the iron opacity bump causes the disk to be convective unstable. Turbulence generated by convection puffs up the disk due to additional turbulent pressure support and enhances the local angular momentum transport. This also results in strong fluctuations in surface density and heating of the disk. The opacity drops with increasing temperature and convection is suppressed. The disk cools down and the whole process repeats again. This causes strong oscillations of the disk scale height and luminosity variations by more than a factor of $\approx 3-6$ over a few years' timescale. Since the iron opacity bump will move to different locations of the disk for black holes with different masses and accretion rates, we suggest that this is a physical mechanism that can explain the variability of AGN with a wide range of amplitudes over a time scale of years to decades.


title: ESL: Entropy-guided Self-supervised Learning for Domain Adaptation in Semantic Segmentation

abstract: While fully-supervised deep learning yields good models for urban scene semantic segmentation, these models struggle to generalize to new environments with different lighting or weather conditions for instance. In addition, producing the extensive pixel-level annotations that the task requires comes at a great cost. Unsupervised domain adaptation (UDA) is one approach that tries to address these issues in order to make such systems more scalable. In particular, self-supervised learning (SSL) has recently become an effective strategy for UDA in semantic segmentation. At the core of such methods lies `pseudo-labeling', that is, the practice of assigning high-confident class predictions as pseudo-labels, subsequently used as true labels, for target data. To collect pseudo-labels, previous works often rely on the highest softmax score, which we here argue as an unfavorable confidence measurement. In this work, we propose Entropy-guided Self-supervised Learning (ESL), leveraging entropy as the confidence indicator for producing more accurate pseudo-labels. On different UDA benchmarks, ESL consistently outperforms strong SSL baselines and achieves state-of-the-art results.


title: Does it matter how well I know what you're thinking? Opponent Modelling in an RTS game

abstract: Opponent Modelling tries to predict the future actions of opponents, and is required to perform well in multi-player games. There is a deep literature on learning an opponent model, but much less on how accurate such models must be to be useful. We investigate the sensitivity of Monte Carlo Tree Search (MCTS) and a Rolling Horizon Evolutionary Algorithm (RHEA) to the accuracy of their modelling of the opponent in a simple Real-Time Strategy game. We find that in this domain RHEA is much more sensitive to the accuracy of an opponent model than MCTS. MCTS generally does better even with an inaccurate model, while this will degrade RHEA's performance. We show that faced with an unknown opponent and a low computational budget it is better not to use any explicit model with RHEA, and to model the opponent's actions within the tree as part of the MCTS algorithm.


title: Signatures of secondary acceleration in neutrino flares

abstract: High-energy neutrino flares are interesting prospective counterparts to photon flares, as their detection would guarantee the presence of accelerated hadrons within a source, provide precious information about cosmic-ray acceleration and interactions, and thus impact the subsequent modeling of non-thermal emissions in explosive transients. In these sources, photomeson production can be efficient, producing a large amount of secondary particles, such as charged pions and muons, that decay and produce high-energy neutrinos. Before their decay, secondary particles can experience energy losses and acceleration, which can impact high-energy neutrino spectra and thus affect their detectability. In this work, we focus on the impact of secondary acceleration. We consider a one zone model, characterized mainly by a variability timescale $t_{\rm var}$, a luminosity $L_{\rm bol}$, a bulk Lorentz factor $\Gamma$. The mean magnetic field $B$ is deduced from these parameters. The photon field is modeled by a broken power-law. This generic model allows to evaluate systematically the maximum energy of high-energy neutrinos in the parameter space of explosive transients, and shows that it could be strongly affected by secondary acceleration for a large number of source categories. In order to determine the impact of secondary acceleration on the high-energy neutrino spectrum and in particular on its peak energy and flux, we complement these estimates by several case studies. We show that secondary acceleration can increase the maximum neutrino flux, and produce a secondary peak at the maximum energy in the case of efficient acceleration. Secondary acceleration could therefore enhance the detectability of very-high-energy neutrinos, that will be the target of next generation neutrino detectors such as KM3NeT, IceCube-Gen2, POEMMA or GRAND.


title: Predicting Livelihood Indicators from Community-Generated Street-Level Imagery

abstract: Major decisions from governments and other large organizations rely on measurements of the populace's well-being, but making such measurements at a broad scale is expensive and thus infrequent in much of the developing world. We propose an inexpensive, scalable, and interpretable approach to predict key livelihood indicators from public crowd-sourced street-level imagery. Such imagery can be cheaply collected and more frequently updated compared to traditional surveying methods, while containing plausibly relevant information for a range of livelihood indicators. We propose two approaches to learn from the street-level imagery: (1) a method that creates multi-household cluster representations by detecting informative objects and (2) a graph-based approach that captures the relationships between images. By visualizing what features are important to a model and how they are used, we can help end-user organizations understand the models and offer an alternate approach for index estimation that uses cheaply obtained roadway features. By comparing our results against ground data collected in nationally-representative household surveys, we demonstrate the performance of our approach in accurately predicting indicators of poverty, population, and health and its scalability by testing in two different countries, India and Kenya. Our code is available at https://github.com/sustainlab-group/mapillarygcn.


title: Simultaneous multi-telescope observations of FRB 121102

abstract: We present 11 detections of FRB 121102 in ~3 hours of observations during its 'active' period on the 10th of September 2019. The detections were made using the newly deployed MeerTRAP system and single pulse detection pipeline at the MeerKAT radio telescope in South Africa. Fortuitously, the Nancay radio telescope observations on this day overlapped with the last hour of MeerKAT observations and resulted in 4 simultaneous detections. The observations with MeerKAT's wide band receiver, which extends down to relatively low frequencies (900-1670 MHz usable L-band range), have allowed us to get a detailed look at the complex frequency structure, intensity variations and frequency-dependent sub-pulse drifting. The drift rates we measure for the full-band and sub-banded data are consistent with those published between 600-6500 MHz with a slope of -0.147 +/- 0.014 ms^-1. Two of the detected bursts exhibit fainter 'precursors' separated from the brighter main pulse by ~28 ms and ~34 ms. A follow-up multi-telescope campaign on the 6th and 8th October 2019 to better understand these frequency drifts and structures over a wide and continuous band was undertaken. No detections resulted, indicating that the source was 'inactive' over a broad frequency range during this time.


title: Constraints on low-energy effective theories from weak cosmic censorship

abstract: We examine the weak cosmic censorship conjecture (WCCC) for the extremal charged black hole in possible generalizations of Einstein-Maxwell theory due to the higher order corrections, up to fourth-derivative terms. Our derivation is based on Wald's gedanken experiment to destroy an extremal black hole. We find that, provided the null energy condition for the falling matter, the WCCC is preserved for all possible generalizations. Thus, the WCCC cannot serve as a constraint to the higher order effective theories. We also show that up to first order variations of black hole mass and charge, WCCC is preserved for non-rotating extremal black holes in all $n$-dimensional diffeomorphism-covariant theories of gravity and $U(1)$ gauge field.


title: Ergodicity conditions for general Markov chains in terms of invariant finitely additive measures

abstract: We consider general Markov chains with discrete time in an arbitrary measurable (phase) space and homogeneous in time. Markov chains are defined by the classical transition function which within the framework of the operator treatment generates a conjugate pair of linear Markov operators in the Banach space of measurable bounded functions and in the Banach space of bounded finite additive measures. It is proved that the well-known Doeblin condition $ (D) $ of ergodicity (quasi\-compactness) of the Markov chain is equivalent to the condition $ (*) $: all finitely additive invariant measures of the Markov operator are countably additive i.e. there are no invariant purely finitely additive measures. Under some assumptions, it is proved that the conditions $ (D) $ and $ (*) $ are also equivalent to the condition $ (**) $: the set of invariant finitely additive measures of a Markov operator is finite-dimensional. Ergodic theorems are given.


title: Fast estimation of aperture-mass statistics I: aperture mass variance and an application to the CFHTLenS data

abstract: We explore an alternative method to the usual shear correlation function approach for the estimation of aperture mass statistics in weak lensing survey data. Our approach builds on the direct estimator method of Schneider (1998). In this paper, to test and validate the methodology, we focus on the aperture mass dispersion. After computing the signal and noise for a weighted set of measured ellipticites we show how the direct estimator can be made into a linear order algorithm that enables a fast and efficient computation. We then investigate the applicability of the direct estimator approach in the presence of a real survey mask with holes and chip gaps. For this we use a large ensemble of full ray-tracing mock simulations. By using various weighting schemes for combining information from different apertures we find that inverse variance weighting the individual aperture estimates with an aperture completeness greater than 70 per cent coverage yields an answer that is in close agreement with the standard correlation function approach. We then apply this approach to the CFHTLenS as pilot scheme and find that our method recovers to high accuracy the Kilbinger (2013) result for the variance of both the E and B mode signal, after we correct the catalogue for the shear bias in the lensfit algorithm for pairs closer than 9". We then explore the cosmological information content of the direct estimator using the Fisher information approach. We show that there is a only modest loss in cosmological information from the rejection of apertures that are of low completeness. This method unlocks the door to fast and efficient methods for recovering higher order aperture mass statistics in linear order operations.


title: Runtime Adaptation in Wireless Sensor Nodes Using Structured Learning

abstract: Markov Decision Processes (MDPs) provide important capabilities for facilitating the dynamic adaptation and self-optimization of cyber physical systems at runtime. In recent years, this has primarily taken the form of Reinforcement Learning (RL) techniques that eliminate some MDP components for the purpose of reducing computational requirements. In this work, we show that recent advancements in Compact MDP Models (CMMs) provide sufficient cause to question this trend when designing wireless sensor network nodes. In this work, a novel CMM-based approach to designing self-aware wireless sensor nodes is presented and compared to Q-Learning, a popular RL technique. We show that a certain class of CPS nodes is not well served by RL methods, and contrast RL versus CMM methods in this context. Through both simulation and a prototype implementation, we demonstrate that CMM methods can provide significantly better runtime adaptation performance relative to Q-Learning, with comparable resource requirements.


title: The Landscape of the Proximal Point Method for Nonconvex-Nonconcave Minimax Optimization

abstract: Minimax optimization has become a central tool in machine learning with applications in robust optimization, reinforcement learning, GANs, etc. These applications are often nonconvex-nonconcave, but the existing theory is unable to identify and deal with the fundamental difficulties this poses. In this paper, we study the classic proximal point method (PPM) applied to nonconvex-nonconcave minimax problems. We find that a classic generalization of the Moreau envelope by Attouch and Wets provides key insights. Critically, we show this envelope not only smooths the objective but can convexify and concavify it based on the level of interaction present between the minimizing and maximizing variables. From this, we identify three distinct regions of nonconvex-nonconcave problems. When interaction is sufficiently strong, we derive global linear convergence guarantees. Conversely when the interaction is fairly weak, we derive local linear convergence guarantees with a proper initialization. Between these two settings, we show that PPM may diverge or converge to a limit cycle.


title: Algorithmic Aspects of Temporal Betweenness

abstract: The betweenness centrality of a graph vertex measures how often this vertex is visited on shortest paths between other vertices of the graph. In the analysis of many real-world graphs or networks, betweenness centrality of a vertex is used as an indicator for its relative importance in the network. In particular, it is among the most popular tools in social network analysis. In recent years, a growing number of real-world networks is modeled as temporal graphs, where we have a fixed set of vertices and there is a finite discrete set of time steps and every edge might be present only at some time steps. While shortest paths are straightforward to define in static graphs, temporal paths can be considered "optimal" with respect to many different criteria, including length, arrival time, and overall travel time (shortest, foremost, and fastest paths). This leads to different concepts of temporal betweenness centrality and we provide a systematic study of temporal betweenness variants based on various concepts of optimal temporal paths. Computing the betweenness centrality for vertices in a graph is closely related to counting the number of optimal paths between vertex pairs. We show that counting foremost and fastest paths is computationally intractable (#P-hard) and hence the computation of the corresponding temporal betweenness values is intractable as well. For shortest paths and two selected special cases of foremost paths, we devise polynomial-time algorithms for temporal betweenness computation. Moreover, we also explore the distinction between strict (ascending time labels) and non-strict (non-descending time labels) time labels in temporal paths. In our experiments with established real-world temporal networks, we demonstrate the practical effectiveness of our algorithms, compare the various betweenness concepts, and derive recommendations on their practical use.


title: On Adversarial Bias and the Robustness of Fair Machine Learning

abstract: Optimizing prediction accuracy can come at the expense of fairness. Towards minimizing discrimination against a group, fair machine learning algorithms strive to equalize the behavior of a model across different groups, by imposing a fairness constraint on models. However, we show that giving the same importance to groups of different sizes and distributions, to counteract the effect of bias in training data, can be in conflict with robustness. We analyze data poisoning attacks against group-based fair machine learning, with the focus on equalized odds. An adversary who can control sampling or labeling for a fraction of training data, can reduce the test accuracy significantly beyond what he can achieve on unconstrained models. Adversarial sampling and adversarial labeling attacks can also worsen the model's fairness gap on test data, even though the model satisfies the fairness constraint on training data. We analyze the robustness of fair machine learning through an empirical evaluation of attacks on multiple algorithms and benchmark datasets.


title: Measuring the thermodynamic cost of timekeeping

abstract: All clocks, in some form or another, use the evolution of nature towards higher entropy states to quantify the passage of time. Due to the statistical nature of the second law and corresponding entropy flows, fluctuations fundamentally limit the performance of any clock. This suggests a deep relation between the increase in entropy and the quality of clock ticks. Indeed, minimal models for autonomous clocks in the quantum realm revealed that a linear relation can be derived, where for a limited regime every bit of entropy linearly increases the accuracy of quantum clocks. But can such a linear relation persist as we move towards a more classical system? We answer this in the affirmative by presenting the first experimental investigation of this thermodynamic relation in a nanoscale clock. We stochastically drive a nanometer-thick membrane and read out its displacement with a radio-frequency cavity, allowing us to identify the ticks of a clock. We show theoretically that the maximum possible accuracy for this classical clock is proportional to the entropy created per tick, similar to the known limit for a weakly coupled quantum clock but with a different proportionality constant. We measure both the accuracy and the entropy. Once non-thermal noise is accounted for, we find that there is a linear relation between accuracy and entropy and that the clock operates within an order of magnitude of the theoretical bound.


title: To Pretrain or Not to Pretrain: Examining the Benefits of Pretraining on Resource Rich Tasks

abstract: Pretraining NLP models with variants of Masked Language Model (MLM) objectives has recently led to a significant improvements on many tasks. This paper examines the benefits of pretrained models as a function of the number of training samples used in the downstream task. On several text classification tasks, we show that as the number of training examples grow into the millions, the accuracy gap between finetuning BERT-based model and training vanilla LSTM from scratch narrows to within 1%. Our findings indicate that MLM-based models might reach a diminishing return point as the supervised data size increases significantly.


title: A systematic review and taxonomy of explanations in decision support and recommender systems

abstract: With the recent advances in the field of artificial intelligence, an increasing number of decision-making tasks are delegated to software systems. A key requirement for the success and adoption of such systems is that users must trust system choices or even fully automated decisions. To achieve this, explanation facilities have been widely investigated as a means of establishing trust in these systems since the early years of expert systems. With today's increasingly sophisticated machine learning algorithms, new challenges in the context of explanations, accountability, and trust towards such systems constantly arise. In this work, we systematically review the literature on explanations in advice-giving systems. This is a family of systems that includes recommender systems, which is one of the most successful classes of advice-giving software in practice. We investigate the purposes of explanations as well as how they are generated, presented to users, and evaluated. As a result, we derive a novel comprehensive taxonomy of aspects to be considered when designing explanation facilities for current and future decision support systems. The taxonomy includes a variety of different facets, such as explanation objective, responsiveness, content and presentation. Moreover, we identified several challenges that remain unaddressed so far, for example related to fine-grained issues associated with the presentation of explanations and how explanation facilities are evaluated.


title: Dynamics and large deviation transitions of the XOR-Fredrickson-Andersen kinetically constrained model

abstract: We study a one-dimensional classical stochastic kinetically constrained model (KCM) inspired by Rydberg atoms in their "facilitated" regime, where sites can flip only if a single of their nearest neighbours is excited. We call this model "XOR-FA" to distinguish it from the standard Fredrickson-Andersen (FA) model. We describe the dynamics of the XOR-FA model, including its relation to simple exclusion processes in its domain wall representation. The interesting relaxation dynamics of the XOR-FA is related to the prominence of large dynamical fluctuations that lead to phase transitions between active and inactive dynamical phases as in other KCMs. By means of numerical tensor network methods we study in detail such transitions in the dynamical large deviation regime.


title: Singularities and Phenomenological aspects of Asymptotic Safe Gravity

abstract: Asymptotic Safety (AS) Program for quantum gravity keeps the same fields and symmetries with General Relativity and studies the associated gravitational action as a fundamental part of the complete theory at the nonperturbative level with the help of functional renormalization group (RG) techniques. An important phenomenological task that can test the new point of view of AS approach is the discovery of RG improved cosmologies and black holes. In this work, we analyze the properties of recently found non-singular spherically symmetric and non-singular cosmological solutions. Furthermore, we derive a novel consistent set of modified Einstein field equations, in the spirit of AS, which respects the Bianchi identities. This new set of equations extend previously published modified Einstein equations which arise by adding appropriate covariant kinetic terms to the action.


title: Targeted Maximum Likelihood Estimation of Community-based Causal Effect of Community-Level Stochastic Interventions

abstract: Unlike the commonly used parametric regression models such as mixed models, that can easily violate the required statistical assumptions and result in invalid statistical inference, target maximum likelihood estimation allows more realistic data-generative models and provides double-robust, semi-parametric and efficient estimators. Target maximum likelihood estimators (TMLEs) for the causal effect of a community-level static exposure were previously proposed by Balzer et al. In this manuscript, we build on this work and present identifiability results and develop two semi-parametric efficient TMLEs for the estimation of the causal effect of the single time-point community-level stochastic intervention whose assignment mechanism can depend on measured and unmeasured environmental factors and its individual-level covariates. The first community-level TMLE is developed under a general hierarchical non-parametric structural equation model, which can incorporate pooled individual-level regressions for estimating the outcome mechanism. The second individual-level TMLE is developed under a restricted hierarchical model in which the additional assumption of no covariate interference within communities holds. The proposed TMLEs have several crucial advantages. First, both TMLEs can make use of individual level data in the hierarchical setting, and potentially reduce finite sample bias and improve estimator efficiency. Second, the stochastic intervention framework provides a natural way for defining and estimating casual effects where the exposure variables are continuous or discrete with multiple levels, or even cannot be directly intervened on. Also, the positivity assumption needed for our proposed causal parameters can be weaker than the version of positivity required for other casual parameters.


title: Anomaly-free 2HDMs with a gauged abelian symmetry and two generations of right-handed neutrinos

abstract: Despite the popularity of two-Higgs-doublet models (2HDMs) with a gauged abelian flavor symmetry, the allowed Yukawa textures have only been partly explored so far. In this work, we classify and compare every anomaly-free instance, in the case of having two generations of right-handed neutrinos and a type-I seesaw mechanism. We found 16 valid implementations in total, out of which 11 agree well with current experimental bounds. To our knowledge, neither of these models have been considered previously.


title: A commutator lemma for confined subgroups and applications to groups acting on rooted trees

abstract: A subgroup $H$ of a group $G$ is confined if the $G$-orbit of $H$ under conjugation is bounded away from the trivial subgroup in the space $\operatorname{Sub}(G)$ of subgroups of $G$. We prove a commutator lemma for confined subgroups. For groups of homeomorphisms, this provides the exact analogue for confined subgroups (hence in particular for URSs) of the classical commutator lemma for normal subgroups: if $G$ is a group of homeomorphisms of a Hausdorff space $X$ and $H$ is a confined subgroup of $G$, then $H$ contains the derived subgroup of the rigid stabilizer of some open subset of $X$. We apply this commutator lemma in the setting of groups acting on rooted trees. We prove a theorem describing the structure of URSs of weakly branch groups and of their non-topologically free minimal actions. Among the applications of these results, we show: 1) if $G$ is a finitely generated branch group, the $G$-action on $\partial T$ has the smallest possible growth among all faithful $G$-actions; 2) if $G$ is a finitely generated branch group, then every embedding from $G$ into a group of homeomorphisms of strongly bounded type (e.g. a bounded automaton group) must be spatially realized; 3) if $G$ is a finitely generated weakly branch group, then $G$ does not embed into the group IET of interval exchange transformations.


title: Surface Brightness Fluctuations for constraining the chemical enrichment of massive galaxies

abstract: Based on very deep photometry, Surface Brightness Fluctuations (SBF) have been traditionally used to determine galaxy distances. We have recently computed SBF spectra of stellar populations at moderately high resolution, which are fully based on empirical stars. We show that the SBF spectra provide an unprecedented potential for stellar population studies that, so far, have been tackled on the basis of the mean fluxes. We find that the SBFs are able to unveil metal-poor stellar components at the one percent level, which are not possible to disentangle with the standard analysis. As these metal-poor components correspond to the first stages of the chemical enrichment, the SBF analysis provides stringent constrains on the quenching epoch.


title: Feature Space Saturation during Training

abstract: We propose layer saturation - a simple, online-computable method for analyzing the information processing in neural networks. First, we show that a layer's output can be restricted to the eigenspace of its variance matrix without performance loss. We propose a computationally lightweight method for approximating the variance matrix during training. From the dimension of its lossless eigenspace we derive layer saturation - the ratio between the eigenspace dimension and layer width. We show that saturation seems to indicate which layers contribute to network performance. We demonstrate how to alter layer saturation in a neural network by changing network depth, filter sizes and input resolution. Furthermore, we show that well-chosen input resolution increases network performance by distributing the inference process more evenly across the network.


title: Shape Matters: Understanding the Implicit Bias of the Noise Covariance

abstract: The noise in stochastic gradient descent (SGD) provides a crucial implicit regularization effect for training overparameterized models. Prior theoretical work largely focuses on spherical Gaussian noise, whereas empirical studies demonstrate the phenomenon that parameter-dependent noise -- induced by mini-batches or label perturbation -- is far more effective than Gaussian noise. This paper theoretically characterizes this phenomenon on a quadratically-parameterized model introduced by Vaskevicius et el. and Woodworth et el. We show that in an over-parameterized setting, SGD with label noise recovers the sparse ground-truth with an arbitrary initialization, whereas SGD with Gaussian noise or gradient descent overfits to dense solutions with large norms. Our analysis reveals that parameter-dependent noise introduces a bias towards local minima with smaller noise variance, whereas spherical Gaussian noise does not. Code for our project is publicly available.


title: Evidence of a Type Ia Progenitor for Supernova Remnant 3C 397

abstract: The explosive origin of the young supernova remnant (SNR) 3C 397 (G41.1-0.3) is debated. Its elongated morphology and proximity to a molecular cloud are suggestive of a core-collapse (CC) SN origin, yet recent X-ray studies of heavy metals show chemical yields and line centroid energies consistent with a Type Ia SN. In this paper, we analyze the full X-ray spectrum from 0.7-10 keV of 3C 397 observed with Suzaku and compare the line centroid energies, fluxes, and elemental abundances of intermediate-mass and heavy metals (Mg to Ni) to Type Ia and CC hydrodynamical model predictions. Based on the results, we conclude that 3C 397 likely arises from an energetic Type Ia explosion in a high-density ambient medium, and we show that the progenitor was a near Chandrasekhar mass white dwarf.


title: The Importance of Low Latency to Order Book Imbalance Trading Strategies

abstract: There is a pervasive assumption that low latency access to an exchange is a key factor in the profitability of many high-frequency trading strategies. This belief is evidenced by the "arms race" undertaken by certain financial firms to co-locate with exchange servers. To the best of our knowledge, our study is the first to validate and quantify this assumption in a continuous double auction market with a single exchange similar to the New York Stock Exchange. It is not feasible to conduct this exploration with historical data in which trader identity and location are not reported. Accordingly, we investigate the relationship between latency of access to order book information and profitability of trading strategies exploiting that information with an agent-based interactive discrete event simulation in which thousands of agents pursue archetypal trading strategies. We introduce experimental traders pursuing a low-latency order book imbalance (OBI) strategy in a controlled manner across thousands of simulated trading days, and analyze OBI trader profit while varying distance (latency) from the exchange. Our experiments support that latency is inversely related to profit for the OBI traders, but more interestingly show that latency rank, rather than absolute magnitude, is the key factor in allocating returns among agents pursuing a similar strategy.


title: Voltage-tunable superconducting resonators: a platform for random access quantum memory

abstract: In quantum computing architectures, one important factor is the trade-off between the need to couple qubits to each other and to an external drive and the need to isolate them well enough in order to protect the information for an extended period of time. In the case of superconducting circuits, one approach is to utilize fixed frequency qubits coupled to coplanar waveguide resonators such that the system can be kept in a configuration that is relatively insensitive to noise. Here, we propose a scalable voltage-tunable quantum memory (QuMem) design concept compatible with superconducting qubit platforms. Our design builds on the recent progress in fabrication of Josephson field effect transistors (JJ-FETs) which use InAs quantum wells. The JJ-FET is incorporated into a tunable coupler between a transmission line and a high-quality resonator in order to control the overall inductance of the coupler. A full isolation of the high-quality resonator can be achieved by turning off the JJ-FET. This could allow for long coherence times and protection of the quantum information inside the storage cavity. The proposed design would facilitate the implementation of random access memory for storage of quantum information in between computational gate operations.


title: Efficient Model-Based Reinforcement Learning through Optimistic Policy Search and Planning

abstract: Model-based reinforcement learning algorithms with probabilistic dynamical models are amongst the most data-efficient learning methods. This is often attributed to their ability to distinguish between epistemic and aleatoric uncertainty. However, while most algorithms distinguish these two uncertainties for learning the model, they ignore it when optimizing the policy, which leads to greedy and insufficient exploration. At the same time, there are no practical solvers for optimistic exploration algorithms. In this paper, we propose a practical optimistic exploration algorithm (H-UCRL). H-UCRL reparameterizes the set of plausible models and hallucinates control directly on the epistemic uncertainty. By augmenting the input space with the hallucinated inputs, H-UCRL can be solved using standard greedy planners. Furthermore, we analyze H-UCRL and construct a general regret bound for well-calibrated models, which is provably sublinear in the case of Gaussian Process models. Based on this theoretical foundation, we show how optimistic exploration can be easily combined with state-of-the-art reinforcement learning algorithms and different probabilistic models. Our experiments demonstrate that optimistic exploration significantly speeds-up learning when there are penalties on actions, a setting that is notoriously difficult for existing model-based reinforcement learning algorithms.


title: Bonus Properties of States of Low Energy

abstract: States of Low Energy (SLE) are exact Hadamard states defined on arbitrary Friedmann-Lema\^{i}tre spacetimes. They are constructed from a fiducial state by minimizing the Hamiltonian's expectation value after averaging with a temporal window function. We show the SLE to be expressible solely in terms of the (state independent) commutator function. They also admit a convergent series expansion in powers of the spatial momentum, both for massive and for massless theories. In the massless case the leading infrared behavior is found to be Minkowski-like for all scale factors. This provides a new cure for the infrared divergences in Friedmann-Lema\^{i}tre spacetimes with accelerated expansion. In consequence, massless SLE are viable candidates for pre-inflationary vacua and in a soluble model are shown to entail a qualitatively correct primordial power spectrum.


title: Multi-Image Summarization: Textual Summary from a Set of Cohesive Images

abstract: Multi-sentence summarization is a well studied problem in NLP, while generating image descriptions for a single image is a well studied problem in Computer Vision. However, for applications such as image cluster labeling or web page summarization, summarizing a set of images is also a useful and challenging task. This paper proposes the new task of multi-image summarization, which aims to generate a concise and descriptive textual summary given a coherent set of input images. We propose a model that extends the image-captioning Transformer-based architecture for single image to multi-image. A dense average image feature aggregation network allows the model to focus on a coherent subset of attributes across the input images. We explore various input representations to the Transformer network and empirically show that aggregated image features are superior to individual image embeddings. We additionally show that the performance of the model is further improved by pretraining the model parameters on a single-image captioning task, which appears to be particularly effective in eliminating hallucinations in the output.


title: Transport coefficients of multi-component mixtures of noble gases based on ab initio potentials. Viscosity and thermal conductivity

abstract: The viscosity and thermal conductivity of binary, ternary and quaternary mixtures of helium, neon, argon, and krypton at low density are computed for wide ranges of temperature and molar fractions, applying the Chapman-Enskog method. Ab initio interatomic potentials are employed in order to calculate the omega-integrals. The relative numerical errors of the viscosity and thermal conductivity do not exceed 1.e-6 and 1.e-5, respectively. The relative uncertainty related to the interatomic potential is about 0.1%. A comparison of the present data with results reported in other papers available in the literature shows a significant improvement of accuracy of the transport coefficients considered here.


title: Causal intersectionality for fair ranking

abstract: In this paper we propose a causal modeling approach to intersectional fairness, and a flexible, task-specific method for computing intersectionally fair rankings. Rankings are used in many contexts, ranging from Web search results to college admissions, but causal inference for fair rankings has received limited attention. Additionally, the growing literature on causal fairness has directed little attention to intersectionality. By bringing these issues together in a formal causal framework we make the application of intersectionality in fair machine learning explicit, connected to important real world effects and domain knowledge, and transparent about technical limitations. We experimentally evaluate our approach on real and synthetic datasets, exploring its behaviour under different structural assumptions.


title: Identify the Optimal Locations of Dedicated Bus Lanes to Improve the Stability of Bus Line

abstract: Bus bunching and unevenly dispersed buses along a bus line lead to a low service level and deteriorate the operational stability of the bus line. Speed adjustment as a control means has been proposed to solve the above problem. But due to the difficulty of adjusting speed in the mixed traffic and the discontinuously distributed dedicated bus lanes (DBLs) in a bus line, the existing methods which coordinate or adjust bus speeds in all road segments encounter serious obstacles in practice. To overcome the above problem, we first take into account the influence of the deployment of the DBLs along a bus line on the performance of the speed adjustment strategy. With the funding limit and the maximally allowable influence on other road traffic as main constraints, a nonlinear mixed integer stochastic mathematic model is formulated aiming at improving the whole stability of the bus line operation. The model can be effectively solved by a newly designed branch and bound algorithm. The optimal regulating speed in a DBL is determined with a process of looking multiple critical time point ahead. Numerical experiment verified the effectiveness and efficiency of our new method. Three findings were obtained as follows. To determine the optimal regulating speed, the number of the critical time points to be looked ahead should not be too big. Increasing the number of DBLs improves the stability of a bus line. The more evenly DBLs are distributed along the bus line, the greater the stability of the bus line will be.


title: Generalized and Scalable Optimal Sparse Decision Trees

abstract: Decision tree optimization is notoriously difficult from a computational perspective but essential for the field of interpretable machine learning. Despite efforts over the past 40 years, only recently have optimization breakthroughs been made that have allowed practical algorithms to find optimal decision trees. These new techniques have the potential to trigger a paradigm shift where it is possible to construct sparse decision trees to efficiently optimize a variety of objective functions without relying on greedy splitting and pruning heuristics that often lead to suboptimal solutions. The contribution in this work is to provide a general framework for decision tree optimization that addresses the two significant open problems in the area: treatment of imbalanced data and fully optimizing over continuous variables. We present techniques that produce optimal decision trees over a variety of objectives including F-score, AUC, and partial area under the ROC convex hull. We also introduce a scalable algorithm that produces provably optimal results in the presence of continuous variables and speeds up decision tree construction by several orders of magnitude relative to the state-of-the art.


title: Neural mass modeling of slow-fast dynamics of seizure initiation and abortion

abstract: Epilepsy is a dynamic and complex neurological disease affecting about 1% of the worldwide population, among which 30% of the patients are drug-resistant. Epilepsy is characterized by recurrent episodes of paroxysmal neural discharges (the so-called seizures), which manifest themselves through a large-amplitude rhythmic activity observed in depth-EEG recordings, in particular in local field potentials (LFPs). The signature characterizing the transition to seizures involves complex oscillatory patterns, which could serve as a marker to prevent seizure initiation by triggering appropriate therapeutic neurostimulation methods. To investigate such protocols, neurophysiological lumped-parameter models at the mesoscopic scale, namely neural mass models, are powerful tools that not only mimic the LFP signals but also give insights on the neural mechanisms related to different stages of seizures. Here, we analyze the multiple time-scale dynamics of a neural mass model and explain the underlying structure of the complex oscillations observed before seizure initiation. We investigate population-specific effects of the stimulation and the dependence of stimulation parameters on synaptic timescales. In particular, we show that intermediate stimulation frequencies (>20 Hz) can abort seizures if the timescale difference is pronounced. Those results have the potential in the design of therapeutic brain stimulation protocols based on the neurophysiological properties of tissue.


title: On a minimal solution for the indefinite truncated multidimensional moment problem

abstract: We will consider the indefinite truncated multidimensional moment problem. Necessary and sufficient conditions for a given truncated multisequence to have a signed representing measure $\mu$ with ${\rm card}\,{\rm supp}\, \mu$ as small as possible are given by the existence of a rank preserving extension of a multivariate Hankel matrix (built from the given truncated multisequence) such that the corresponding associated polynomial ideal is real radical. This result is a special case of a more general characterisation of truncated multisequences with a minimal complex representing measure whose support is symmetric with respect to complex conjugation (which we will call {\it quasi-complex}). One motivation for our results is the fact that positive semidefinite truncated multisequence need not have a positive representing measure. Thus, our main result gives the potential for computing a signed representing measure $\mu = \mu_+ - \mu_-$, where ${\rm card} \,\mu_-$ is small. We illustrate this point on concrete examples.


title: Curvature vs degrees of freedom: The case of the critical 2+1 Horava theory

abstract: We present the interesting case of the 2+1 nonprojectable Horava theory formulated at the critical point, where it does not posses local degrees of freedom. The critical point is defined by the value of a coupling constant of the theory. We discuss how, in spite of the absence of degrees of freedom, the theory admits solutions with nonflat or nonconstant curvature. We consider the theory without cosmological constant and without terms of higher order derivatives, hence this is an effect that can be seen at the same order of 2+1 general relativity. We present an exact nonflat solution that is not asymptotically flat. The presence of solutions with nontrivial curvature seems to be related to the relaxing of the asymptotically flat condition. We discuss that there is no analogue of Newtonian potential in this theory, and a broad class of asymptotically flat geometries leads to the restriction that the only solutions that can be found among them are the flat ones.


title: DeshuffleGAN: A Self-Supervised GAN to Improve Structure Learning

abstract: Generative Adversarial Networks (GANs) triggered an increased interest in problem of image generation due to their improved output image quality and versatility for expansion towards new methods. Numerous GAN-based works attempt to improve generation by architectural and loss-based extensions. We argue that one of the crucial points to improve the GAN performance in terms of realism and similarity to the original data distribution is to be able to provide the model with a capability to learn the spatial structure in data. To that end, we propose the DeshuffleGAN to enhance the learning of the discriminator and the generator, via a self-supervision approach. Specifically, we introduce a deshuffling task that solves a puzzle of randomly shuffled image tiles, which in turn helps the DeshuffleGAN learn to increase its expressive capacity for spatial structure and realistic appearance. We provide experimental evidence for the performance improvement in generated images, compared to the baseline methods, which is consistently observed over two different datasets.


title: Revisiting the Impact of Dust Production from Carbon-Rich Wolf-Rayet Binaries

abstract: We present a dust spectral energy distribution (SED) and binary stellar population analysis revisiting the dust production rates (DPRs) in the winds of carbon-rich Wolf-Rayet (WC) binaries and their impact on galactic dust budgets. DustEM SED models of 19 Galactic WC ``dustars" reveal DPRs of $\dot{M}_d\sim10^{-10}-10^{-6}$ M$_\odot$ yr$^{-1}$ and carbon dust condensation fractions, $\chi_C$, between $0.002 - 40\%$. A large ($0.1 - 1.0$ $\mu$m) dust grain size composition is favored for efficient dustars where $\chi_C\gtrsim1\%$. Results for dustars with known orbital periods verify a power-law relation between $\chi_C$, orbital period, WC mass-loss rate, and wind velocity consistent with predictions from theoretical models of dust formation in colliding-wind binaries. We incorporated dust production into Binary Population and Spectral Synthesis (BPASS) models to analyze dust production rates from WC dustars, asymptotic giant branch stars (AGBs), red supergiants (RSGs), and core-collapse supernovae (SNe). BPASS models assuming constant star formation (SF) and a co-eval $10^6$ M$_\odot$ stellar population were performed at low, Large Magellanic Cloud (LMC)-like, and solar metallicities (Z = 0.001, 0.008, and 0.020). Both constant SF and co-eval models show that SNe are net dust destroyers at all metallicities. Constant SF models at LMC-like metallicities show that AGB stars slightly outproduce WC binaries and RSGs by factors of $2-3$, whereas at solar metallicites WC binaries are the dominant source of dust for $\sim60$ Myr until the onset of AGBs, which match the dust input of WC binaries. Co-eval population models show that for "bursty" SF, AGB stars dominate dust production at late times ($t\gtrsim 70$ Myr).


title: Unsupervised Domain Adaptation for Semantic Segmentation of NIR Images through Generative Latent Search

abstract: Segmentation of the pixels corresponding to human skin is an essential first step in multiple applications ranging from surveillance to heart-rate estimation from remote-photoplethysmography. However, the existing literature considers the problem only in the visible-range of the EM-spectrum which limits their utility in low or no light settings where the criticality of the application is higher. To alleviate this problem, we consider the problem of skin segmentation from the Near-infrared images. However, Deep learning based state-of-the-art segmentation techniques demands large amounts of labelled data that is unavailable for the current problem. Therefore we cast the skin segmentation problem as that of target-independent Unsupervised Domain Adaptation (UDA) where we use the data from the Red-channel of the visible-range to develop skin segmentation algorithm on NIR images. We propose a method for target-independent segmentation where the 'nearest-clone' of a target image in the source domain is searched and used as a proxy in the segmentation network trained only on the source domain. We prove the existence of 'nearest-clone' and propose a method to find it through an optimization algorithm over the latent space of a Deep generative model based on variational inference. We demonstrate the efficacy of the proposed method for NIR skin segmentation over the state-of-the-art UDA segmentation methods on the two newly created skin segmentation datasets in NIR domain despite not having access to the target NIR data. Additionally, we report state-of-the-art results for adaption from Synthia to Cityscapes which is a popular setting in Unsupervised Domain Adaptation for semantic segmentation. The code and datasets are available at https://github.com/ambekarsameer96/GLSS.


title: Online Change Point Detection in Molecular Dynamics With Optical Random Features

abstract: Proteins are made of atoms constantly fluctuating, but can occasionally undergo large-scale changes. Such transitions are of biological interest, linking the structure of a protein to its function with a cell. Atomic-level simulations, such as Molecular Dynamics (MD), are used to study these events. However, molecular dynamics simulations produce time series with multiple observables, while changes often only affect a few of them. Therefore, detecting conformational changes has proven to be challenging for most change-point detection algorithms. In this work, we focus on the identification of such events given many noisy observables. In particular, we show that the No-prior-Knowledge Exponential Weighted Moving Average (NEWMA) algorithm can be used along optical hardware to successfully identify these changes in real-time. Our method does not need to distinguish between the background of a protein and the protein itself. For larger simulations, it is faster than using traditional silicon hardware and has a lower memory footprint. This technique may enhance the sampling of the conformational space of molecules. It may also be used to detect change-points in other sequential data with a large number of features.


title: Extrapolatable Relational Reasoning With Comparators in Low-Dimensional Manifolds

abstract: While modern deep neural architectures generalise well when test data is sampled from the same distribution as training data, they fail badly for cases when the test data distribution differs from the training distribution even along a few dimensions. This lack of out-of-distribution generalisation is increasingly manifested when the tasks become more abstract and complex, such as in relational reasoning. In this paper we propose a neuroscience-inspired inductive-biased module that can be readily amalgamated with current neural network architectures to improve out-of-distribution (o.o.d) generalisation performance on relational reasoning tasks. This module learns to project high-dimensional object representations to low-dimensional manifolds for more efficient and generalisable relational comparisons. We show that neural nets with this inductive bias achieve considerably better o.o.d generalisation performance for a range of relational reasoning tasks. We finally analyse the proposed inductive bias module to understand the importance of lower dimension projection, and propose an augmentation to the algorithmic alignment theory to better measure algorithmic alignment with generalisation.


title: Current relaxation in the Random Resistor cum Tunneling Network Model through First-Passage route : Regimes and Time-scales

abstract: Numerically we study the bulk current relaxation in percolative Random Resistor cum Tunneling Network (RRTN) model through a first-passage route. The RRTN considers an extra semi-classical barrier-crossing process over a voltage threshold within a framework of classical RRN bond percolation model. We identify the different temporal regimes of relaxation and corresponding phenomenological time-scales, which fix up the extents of different regimes. These time-scales were previously identified in refs. \cite{relax-physicaA, aksubh}. We investigate on the distributions of these time-scales and observe that there exists a perfect correlation among them in the thermodynamic limit. We conclude that there exists a single time-scale which controls the RRTN dynamics. The variation of mean first-passage time .vs. system size seems to be due to sub-diffusive motion of charge carrier through the network.


title: A multi-stage looking-ahead holding strategy to stabilize a high-frequency bus line

abstract: If a bus line becomes unstable, passengers waiting time will be lengthened and buses capacities will be mismatched. To stabilize a high-frequency bus line, many holding strategies have been proposed. Among these strategies, some need to take oversimplified assumptions to simplify the formulation of a real bus line; some may choose a myopic holding time because they take into account only the local and currently available headway information. To overcome the above shortcomings, we proposed an adaptive holding strategy which continuously reduces the deviations of the instantaneous headways from the Dynamic Target Headway -DTH. By using DTH, the new strategy can make use of the global headway information to determine a proper holding time. To fully estimate the influence of a holding time on the operation of a bus line in a relative long time period, we introduced the multi-stage looking-ahead mechanism into our new strategy. A detailed bus line simulation system was used to realize the above mechanism and to describe all kinds of complex components of a bus line. The numerical experiment demonstrated the effectiveness of our new strategy. Some meaningful insights were uncovered as follow: a. The number of stages to be looked ahead should not be too small or too big; b. The bigger the range of the action set, the better the performance of our strategy; c. To obtain the optimal performance, the interval between two neighbor holding times should not be too small or too big; d. the bigger the number of the control points, the better the performance of our strategy.


title: Supervised Visualization for Data Exploration

abstract: Dimensionality reduction is often used as an initial step in data exploration, either as preprocessing for classification or regression or for visualization. Most dimensionality reduction techniques to date are unsupervised; they do not take class labels into account (e.g., PCA, MDS, t-SNE, Isomap). Such methods require large amounts of data and are often sensitive to noise that may obfuscate important patterns in the data. Various attempts at supervised dimensionality reduction methods that take into account auxiliary annotations (e.g., class labels) have been successfully implemented with goals of increased classification accuracy or improved data visualization. Many of these supervised techniques incorporate labels in the loss function in the form of similarity or dissimilarity matrices, thereby creating over-emphasized separation between class clusters, which does not realistically represent the local and global relationships in the data. In addition, these approaches are often sensitive to parameter tuning, which may be difficult to configure without an explicit quantitative notion of visual superiority. In this paper, we describe a novel supervised visualization technique based on random forest proximities and diffusion-based dimensionality reduction. We show, both qualitatively and quantitatively, the advantages of our approach in retaining local and global structures in data, while emphasizing important variables in the low-dimensional embedding. Importantly, our approach is robust to noise and parameter tuning, thus making it simple to use while producing reliable visualizations for data exploration.


title: Application of Machine Learning to Predict the Risk of Alzheimer's Disease: An Accurate and Practical Solution for Early Diagnostics

abstract: Alzheimer's Disease (AD) ravages the cognitive ability of more than 5 million Americans and creates an enormous strain on the health care system. This paper proposes a machine learning predictive model for AD development without medical imaging and with fewer clinical visits and tests, in hopes of earlier and cheaper diagnoses. That earlier diagnoses could be critical in the effectiveness of any drug or medical treatment to cure this disease. Our model is trained and validated using demographic, biomarker and cognitive test data from two prominent research studies: Alzheimer's Disease Neuroimaging Initiative (ADNI) and Australian Imaging, Biomarker Lifestyle Flagship Study of Aging (AIBL). We systematically explore different machine learning models, pre-processing methods and feature selection techniques. The most performant model demonstrates greater than 90% accuracy and recall in predicting AD, and the results generalize across sub-studies of ADNI and to the independent AIBL study. We also demonstrate that these results are robust to reducing the number of clinical visits or tests per visit. Using a metaclassification algorithm and longitudinal data analysis we are able to produce a "lean" diagnostic protocol with only 3 tests and 4 clinical visits that can predict Alzheimer's development with 87% accuracy and 79% recall. This novel work can be adapted into a practical early diagnostic tool for predicting the development of Alzheimer's that maximizes accuracy while minimizing the number of necessary diagnostic tests and clinical visits.


title: Information technology & astronomical data in Brazil: Perspectives and proposals

abstract: The Commission on Science and Information Technology (CTCI) of the Brazilian Astronomical Society (SAB) is tasked with assisting the Society on issues of astronomical data management, from its handling and the management of data centres and networks, to technical aspects of the archiving, storage and dissemination of data. In this paper we present a summary of the results of a survey recently conducted by the Commission to diagnose the status of several data-related issues within the Brazilian astronomical community, as well as some proposals derived therefrom.