title: Identification of AKARI infrared sources by Deep HSC Optical Survey: Construction of New Band-Merged Catalogue in the NEP-Wide field

abstract: The north ecliptic pole (NEP) field is a natural deep field location for many satellite observations. It has been targeted manytimes since it was surveyed by the AKARI space telescope with its unique wavelength coverage from the near- to mid-infrared(mid-IR). Many follow-up observations have been carried out and made this field one of the most frequently observed areas witha variety of facilities, accumulating abundant panchromatic data from X-ray to radio wavelength range. Recently, a deep opticalsurvey with the Hyper Suprime-Cam (HSC) at the Subaru telescope covered the NEP-Wide (NEPW) field, which enabled us toidentify faint sources in the near- and mid-IR bands, and to improve the photometric redshift (photo-z) estimation. In this work,we present newly identified AKARI sources by the HSC survey, along with multi-band photometry for 91,861 AKARI sourcesobserved over the NEPW field. We release a new band-merged catalogue combining various photometric data from GALEXUV to the submillimetre (sub-mm) bands (e.g., Herschel/SPIRE, JCMT/SCUBA-2). About 20,000 AKARI sources are newlymatched to the HSC data, most of which seem to be faint galaxies in the near- to mid-infrared AKARI bands. This cataloguemotivates a variety of current research, and will be increasingly useful as recently launched (eROSITA/ART-XC) and futurespace missions (such as JWST, Euclid, and SPHEREx) plan to take deep observations in the NEP field.


title: Classification of three generation models by orbifolding magnetized $T^2 \times T^2$

abstract: We study orbifolding by the $\mathbb{Z}_2^{\rm (per)}$ permutaion of $T^2_1 \times T^2_2$ with magnetic fluxes and its twisted orbifolds. We classify the possible three generation models which lead to non-vanishing Yukawa couplings on the magnetized $T^2_1 \times T^2_2$ and orbifolds including the $\mathbb{Z}_2^{\rm (per)}$ permutation and $\mathbb{Z}_2^{\rm (t)}$ twist. We also study the modular symmetry on such orbifold models. As an illustrating model, we examine the realization of quark masses and mixing angles.


title: Forecasting Black Sigatoka Infection Risks with Latent Neural ODEs

abstract: Black Sigatoka disease severely decreases global banana production, and climate change aggravates the problem by altering fungal species distributions. Due to the heavy financial burden of managing this infectious disease, farmers in developing countries face significant banana crop losses. Though scientists have produced mathematical models of infectious diseases, adapting these models to incorporate climate effects is difficult. We present MR. NODE (Multiple predictoR Neural ODE), a neural network that models the dynamics of black Sigatoka infection learnt directly from data via Neural Ordinary Differential Equations. Our method encodes external predictor factors into the latent space in addition to the variable that we infer, and it can also predict the infection risk at an arbitrary point in time. Empirically, we demonstrate on historical climate data that our method has superior generalization performance on time points up to one month in the future and unseen irregularities. We believe that our method can be a useful tool to control the spread of black Sigatoka.


title: Jos\'e Monteiro da Rocha (1734-1819) and his work of 1782 on the determination of comet orbits

abstract: In 1782 Jos\'e Monteiro da Rocha, astronomer and professor of the University of Coimbra, presented in a public session of the Royal Academy of Sciences of Lisbon a memoir on the problem of the determination of the comets' orbits. Only in 1799, the "Determina\c{c}\~ao das Orbitas dos Cometas" (Determination of the orbits of comets) would be published in the Academy's memoires. In that work, Monteiro da Rocha presents a method for solving the problem of the determination of the parabolic orbit of a comet from three observations. Monteiro da Rocha's method is essentially the same method proposed by Olbers and published under von Zach's sponsorship two years before, in 1797. To have been written and published in Portuguese was certainly a hindrance for its dissemination among the international astronomical community. In this article, we intend to present Monteiro da Rocha's method and trying to see to what extent Gomes Teixeira's assertion (Teixeira 1934) that Monteiro da Rocha and Olbers must figure together in the history of astronomy, as the first inventors of a practical and easy method for the determination of parabolic orbits of comets, is justified.


title: Drinfeld Hecke algebras for symmetric groups in positive characteristic

abstract: We investigate deformations of skew group algebras arising from the action of the symmetric group on polynomial rings over fields of arbitrary characteristic. Over the real or complex numbers, Lusztig's graded affine Hecke algebra and analogs are all isomorphic to Drinfeld Hecke algebras, which include the symplectic reflection algebras and rational Cherednik algebras. Over fields of prime characteristic, new deformations arise that capture both a disruption of the group action and also a disruption of the commutativity relations defining the polynomial ring. We classify deformations for the symmetric group acting in its natural (reducible) reflection representation.


title: Generalized solutions to bounded-confidence models

abstract: Bounded-confidence models in social dynamics describe multi-agent systems, where each individual interacts only locally with others. Several models are written as systems of ordinary differential equations with discontinuous right-hand side: this is a direct consequence of restricting interactions to a bounded region with non-vanishing strength at the boundary. Various works in the literature analyzed properties of solutions, such as barycenter invariance and clustering. On the other side, the problem of giving a precise definition of solution, from an analytical point of view, was often overlooked. However, a rich literature proposing different concepts of solution to discontinuous differential equations is available. Using several concepts of solution, we show how existence is granted under general assumptions, while uniqueness may fail even in dimension one, but holds for almost every initial conditions. Consequently, various properties of solutions depend on the used definition and initial conditions.


title: Lectures on Hausdorff and Gromov-Hausdorff Distance Geometry

abstract: The course was given at Peking University, Fall 2019. We discuss the following subjects: (1) Introduction to general topology, hyperspaces, metric and pseudometric spaces, graph theory. (2) Graphs in metric spaces, minimum spanning tree, Steiner minimal tree, Gromov minimal filling. (3) Hausdorff distance, Vietoris topology, Limits theory, inheritance of completeness, total boundedness, compactness by hyperspaces. (4) Gromov-Hausdorff distance, triangle inequality, positive definiteness for isometry classes of compact spaces, counterexample for boundedly compact spaces. (5) Gromov-Hausdorff distance for separable spaces in terms of their isometric images in \ell_\infty, correspondences, Gromov-Hausdorff distance in terms of correspondences. (6) Epsilon-isometries and Gromov-Hausdorff distance. (7) Irreducible correspondences and Gromov-Hausdorff distance. (8) Gromov-Hausdorff convergence, inheritance of metric and topological properties while Gromov-Hausdorff convergence. (9) Gromov-Hausdorff space (GH-space), optimal correspondences, existence of closed optimal correspondences for compact metric spaces, GH-space is geodesic. (10) Cover number, packing number, total boundedness, completeness, and separability of GH-space. (11) mst-spectrum in terms of GH-distances to simplexes, Steiner problem in GH-space. (12) GH-distance to simplexes with more points, GH-distance to simplexes with at most the same number of points. (13) Generalized Borsuk problem, solution of Generalized Borsuk problem in terms of GH-distances, clique covering number and chromatic number of simple graphs, their dualities, calculating these numbers in terms of GH-distances.


title: Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

abstract: PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.


title: Maximally-hard spectra from diffusive shock-acceleration in pulsar-wind nebulae

abstract: The processes leading to the exceptionally hard radio spectra of pulsar-wind nebulae (PWNe) are not yet understood. Radio photon spectral indices among $29$ PWNe from the literature show an approximately normal, $\alpha=0.2\pm0.2$ distribution. We present $\sim 3\sigma$ evidence for a distinct sub-population of PWNe, with a hard spectrum $\alpha=0.01\pm0.06$ near the termination shock and significantly softer elsewhere, possibly due to a recent evacuation of the shock surroundings. Such spectra, especially in the hard sub-population, suggest a Fermi process, such as diffusive shock acceleration, at its extreme, $\alpha=0$ limit. We show that this limit is approached for sufficiently anisotropic small-angle scattering, enhanced on either side of the shock for particles approaching the shock front. In the upstream, the spectral hardening is mostly associated with an enhanced energy gain, possibly driven by the same beamed particles crossing the shock. Downstream, the main effect is a diminished escape probability, but this lowers the acceleration efficiency to $\lesssim25\%$ for $\alpha=0.3$ and $\lesssim1\%$ for $\alpha=0.03$.


title: MaX-DeepLab: End-to-End Panoptic Segmentation with Mask Transformers

abstract: We present MaX-DeepLab, the first end-to-end model for panoptic segmentation. Our approach simplifies the current pipeline that depends heavily on surrogate sub-tasks and hand-designed components, such as box detection, non-maximum suppression, thing-stuff merging, etc. Although these sub-tasks are tackled by area experts, they fail to comprehensively solve the target task. By contrast, our MaX-DeepLab directly predicts class-labeled masks with a mask transformer, and is trained with a panoptic quality inspired loss via bipartite matching. Our mask transformer employs a dual-path architecture that introduces a global memory path in addition to a CNN path, allowing direct communication with any CNN layers. As a result, MaX-DeepLab shows a significant 7.1% PQ gain in the box-free regime on the challenging COCO dataset, closing the gap between box-based and box-free methods for the first time. A small variant of MaX-DeepLab improves 3.0% PQ over DETR with similar parameters and M-Adds. Furthermore, MaX-DeepLab, without test time augmentation, achieves new state-of-the-art 51.3% PQ on COCO test-dev set. Code is available at https://github.com/google-research/deeplab2.


title: Design of an optomagnonic crystal: towards optimal magnon-photon mode matching at the microscale

abstract: We put forward the concept of an optomagnonic crystal: a periodically patterned structure at the microscale based on a magnetic dielectric, which can co-localize magnon and photon modes. The co-localization in small volumes can result in large values of the photon-magnon coupling at the single quanta level, which opens perspectives for quantum information processing and quantum conversion schemes with these systems. We study theoretically a simple geometry consisting of a one-dimensional array of holes with an abrupt defect, considering the ferrimagnet Yttrium Iron Garnet (YIG) as the basis material. We show that both magnon and photon modes can be localized at the defect, and use symmetry arguments to select an optimal pair of modes in order to maximize the coupling. We show that an optomagnonic coupling in the kHz range is achievable in this geometry, and discuss possible optimization routes in order to improve both coupling strengths and optical losses.


title: AGN Triality of Triple Mergers: Detection of Faint X-ray Point Sources

abstract: We present results from our X-ray analysis of the first systematic search for triple AGN in nearby (z<0.077) triple galaxy mergers. We analyze archival Chandra observations of 7 triple galaxy mergers with BAYMAX (Bayesian Analysis of Multiple AGN in X-rays), fitting each observation with single, dual, and triple X-ray point source models. In doing so, we conclude that 1 triple merger has one X-ray point source (SDSS J0858+1822, although it's unlikely to be an AGN); 5 triple mergers are likely composed of two X-ray point sources (NGC 3341, SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356$-$1016); and one system is composed of three X-ray point sources (SDSS J0849+1114). By fitting the individual X-ray spectra of each point source, we analyze the 2-7 keV luminosities as well as the levels of obscuration associated with each potential AGN. We find that 4/5 dual X-ray point source systems have primary and secondary point sources with bright X-ray luminosities (L_2-7 kev >10^40 erg s^-1), possibly associated with 4 new undetected dual AGN. The dual and triple point source systems are found to have physical separations between 3-9 kpc and flux ratios between 2x10^-3 - 0.84. A multi-wavelength analysis to determine the origin of the X-ray point sources discovered in this work is presented in our companion paper (Foord et al. 2020c).


title: Anomalous Floquet Chiral Topological Superconductivity in a Topological Insulator Sandwich Structure

abstract: We show that Floquet chiral topological superconductivity arises naturally in Josephson junctions made of magnetic topological insulator-superconductor sandwich structures. The Josephson phase modulation associated with an applied bias voltage across the junction drives the system into the anomalous Floquet chiral topological superconductor hosting chiral Majorana edge modes in the quasienergy spectrum, with the bulk Floquet bands carrying zero Chern numbers. The bias voltage acts as a tuning parameter enabling novel dynamical topological quantum phase transitions driving the system into a myriad of exotic Majorana-carrying Floquet topological superconducting phases. Our theory establishes a new paradigm for realizing Floquet chiral topological superconductivity in solid-state systems, which should be experimentally directly accessible.

categories: cond-mat.mes-hall cond-mat.mtrl-sci cond-mat.str-el cond-mat.supr-con quant-ph


title: Solvable Theory of a Strange Metal at the Breakdown of a Heavy Fermi Liquid

abstract: We introduce an effective theory for quantum critical points (QCPs) in heavy fermion systems, involving a change in carrier density without symmetry breaking. Our new theory captures a strongly coupled metallic QCP, leading to robust marginal Fermi liquid transport phenomenology, and associated linear in temperature ($T$) "strange metal" resistivity, all within a controlled large $N$ limit. In the parameter regime of strong damping of emergent bosonic excitations, the QCP also displays a near-universal "Planckian" transport lifetime, $\tau_{\mathrm{tr}}\sim\hbar/(k_BT)$. This is contrasted with the conventional so-called "slave boson" theory of the Kondo breakdown, where the large $N$ limit describes a weak coupling fixed point and non-trivial transport behavior may only be obtained through uncontrolled $1/N$ corrections. We also compute the weak-field Hall coefficient within the effective model as the system is tuned across the transition. We further find that between the two plateaus, reflecting the different carrier densities in the two Fermi liquid phases, the Hall coefficient can develop a peak in the critical crossover regime, like in recent experimental findings, in the parameter regime of weak boson damping.


title: The Exomoon Corridor: Half of all exomoons exhibit TTV frequencies within a narrow window due to aliasing

abstract: Exomoons are expected to produce potentially detectable transit timing variations (TTVs) upon their parent planet. Unfortunately, distinguishing moon-induced TTVs from other sources, in particular planet-planet interactions, has severely impeded its usefulness as a tool for identifying exomoon candidates. A key feature of exomoon TTVs is that they will always be undersampled, due to the simple fact that we can only observe the TTVs once per transit/planetary period. We show that it is possible to analytically express the aliased TTV periodicity as a function of planet and moon period. Further, we show that inverting an aliased TTV period back to a true moon period is fraught with hundreds of harmonic modes. However, a unique aspect of these TTV aliases is that they are predicted to occur at consistently short periods, irrespective of what model one assumes for the underlying moon population. Specifically, 50% of all exomoons are expected to induce TTVs with a period between 2 to 4 cycles, a range that planet-planet interactions rarely manifest at. This provides an exciting and simple tool for quickly identifying exomoons candidates and brings the TTV method back to the fore as an exomoon hunting strategy. Applying this method to the candidate, Kepler-1625b i, reveals that its TTV periodicity centers around the median period expected for exomoons.


title: A Proposed Network to Detect Axion Quark Nugget Dark Matter

abstract: A network of synchronized detectors can increase the likelihood of discovering the QCD axion, within the Axion Quark Nugget (AQN) dark matter model. A similar network can also discriminate the X-rays emitted by the AQNs from the background signal. These networks can provide information on the directionality of the dark matter flux (if any), as well as its velocity distribution, and can therefore test the Standard Halo Model. We show that the optimal configuration to detect AQN-induced axions is a triangular network of stations 100 km apart. For X-rays, the optimal network is an array of tetrahedral units.


title: Holomorphic Anomalies, Fourfolds and Fluxes

abstract: We investigate holomorphic anomalies of partition functions underlying string compactifications on Calabi-Yau fourfolds with background fluxes. For elliptic fourfolds the partition functions have an alternative interpretation as elliptic genera of N=1 supersymmetric string theories in four dimensions, or as generating functions for relative Gromov-Witten invariants of fourfolds with fluxes. We derive the holomorphic anomaly equations by starting from the BCOV formalism of topological strings, and translating them into geometrical terms. The result can be recast into modular and elliptic anomaly equations. As a new feature, as compared to threefolds, we find an extra contribution which is given by a gravitational descendant invariant. This leads to linear terms in the anomaly equations, which support an algebra of derivatives mapping between partition functions of the various flux sectors. These geometric features are mirrored by certain properties of quasi-Jacobi forms. We also offer an interpretation of the physics from the viewpoint of the worldsheet theory.


title: From hydrodynamics to N-body simulations of star clusters: mergers and rotation

abstract: We present a new method to obtain more realistic initial conditions for N-body simulations of young star clusters. We start from the outputs of hydrodynamical simulations of molecular cloud collapse, in which star formation is modelled with sink particles. In our approach, we instantaneously remove gas from these hydrodynamical simulation outputs to mock the end of the gas-embedded phase, induced by stellar feedback. We then enforce a realistic initial mass function by splitting or joining the sink particles based on their mass and position. Such initial conditions contain more consistent information on the spatial distribution and the kinematical and dynamical states of young star clusters, which are fundamental to properly study these systems. For example, by applying our method to a set of previously run hydrodynamical simulations, we found that the early evolution of young star clusters is affected by gas removal and by the early dry merging of sub-structures. This early evolution can either quickly erase the rotation acquired by our (sub-)clusters in their embedded phase or "fuel" it by feeding of angular momentum by sub-structure mergers, before two-body relaxation acts on longer timescales.


title: Infrared interferometric imaging of the compact dust disk around the AGB star HR3126 with the bipolar Toby Jug Nebula

abstract: The asymptotic giant branch (AGB) star HR3126, associated with the arcminute-scale bipolar Toby Jug Nebula, provides a rare opportunity to study the emergence of bipolar structures at the end of the AGB phase. We carried out long-baseline interferometric observations with AMBER and GRAVITY (2--2.45 micron) at the Very Large Telescope Interferometer, speckle interferometric observations with VLT/NACO (2.24 micron), and imaging with SPHERE-ZIMPOL (0.55 micron) and VISIR (7.9--19.5 micron). The images reconstructed in the continuum at 2.1--2.29 micron from the AMBER+GRAVITY data reveal the central star surrounded by an elliptical ring-like structure with a semimajor and semiminor axis of 5.3 and 3.5 mas, respectively. The ring is interpreted as the inner rim of an equatorial dust disk viewed from an inclination angle of ~50 degrees, and its axis is approximately aligned with the bipolar nebula. The disk is surprisingly compact, with an inner radius of a mere 3.5 Rstar (2 au). Our 2-D radiative transfer modeling shows that an optically thick flared disk with silicate grains as large as ~4 micron can reproduce the observed continuum images and the spectral energy distribution. The images obtained in the CO first overtone bands reveal elongated extended emission around the central star, suggesting the oblateness of the star's atmosphere or the presence of a CO gas disk inside the dust cavity. The object is unresolved with SPHERE-ZIMPOL, NACO, and VISIR. If the disk formed together with the bipolar nebula, the grain growth from sub-micron to a few microns should have taken place over the nebula's dynamical age of ~3900 yrs. The non-detection of a companion in the reconstructed images implies that either its 2.2 micron brightness is more than ~30 times lower than that of the red giant or it might have been shredded due to binary interaction.


title: AGN Triality of Triple Mergers: Multi-wavelength Classifications

abstract: We present results from a multi-wavelength analysis searching for multiple AGN systems in nearby (z<0.077) triple galaxy mergers. Combining archival Chandra, SDSS, WISE, and VLA observations, we quantify the rate of nearby triple AGN, as well as investigate possible connections between SMBH accretion and merger environments. Analyzing the multi-wavelength observations of 7 triple galaxy mergers, we find that 1 triple merger has a single AGN (NGC 3341); we discover, for the first time, 4 likely dual AGN (SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356-1016); we confirm one triple AGN system, SDSS J0849+1114; and 1 triple merger in our sample remains ambiguous (SDSS J0858+1822). Analyzing the WISE data, we find a trend of increasing N_H (associated with the primary AGN) as a function of increasing W1-W2 color, reflecting that the motions of gas and dust are coupled in merging environments, where large amount of both can be funneled into the active central region during mergers. Additionally, we find that the one triple AGN system in our sample has the highest levels of N_H and W1-W2 color, while the dual AGN candidates all have lower levels; these results are consistent with theoretical merger simulations that suggest higher levels of nuclear gas are more likely to activate AGN in mergers.


title: Most of the cool CGM of star-forming galaxies is not produced by supernova feedback

abstract: The characterization of the large amount of gas residing in the galaxy halos, the so called circumgalactic medium (CGM), is crucial to understand galaxy evolution across cosmic time. We focus here on the the cool ($T\sim10^4$ K) phase of this medium around star-forming galaxies in the local universe, whose properties and dynamics are poorly understood. We developed semi-analytical parametric models to describe the cool CGM as an outflow of gas clouds from the central galaxy, as a result of supernova explosions in the disc (galactic wind). The cloud motion is driven by the galaxy gravitational pull and by the interactions with the hot ($T\sim10^6$ K) coronal gas. Through a bayesian analysis, we compare the predictions of our models with the data of the COS-Halos and COS-GASS surveys, which provide accurate kinematic information of the cool CGM around more than 40 low-redshift star-forming galaxies, probing distances up to the galaxy virial radii. Our findings clearly show that a supernova-driven outflow model is not suitable to describe the dynamics of the cool circumgalactic gas. Indeed, to reproduce the data, we need extreme scenarios, with initial outflow velocities and mass loading factors that would lead to unphysically high energy coupling from the supernovae to the gas and with supernova efficiencies largely exceeding unity. This strongly suggests that, since the outflows cannot reproduce most of the cool gas absorbers, the latter are likely the result of cosmological inflow in the outer galaxy halos, in analogy to what we have previously found for early-type galaxies.


title: Strategies to Detect Dark-Matter Decays with Line-Intensity Mapping

abstract: The nature of dark matter is a longstanding mystery in cosmology, which can be studied with laboratory or collider experiments, as well as astrophysical and cosmological observations. In this work, we propose realistic and efficient strategies to detect radiative products from dark-matter decays with line-intensity mapping (LIM) experiments. This radiation will behave as a line interloper for the atomic and molecular spectral lines targeted by LIM surveys. The most distinctive signatures of the contribution from dark-matter radiative decays are an extra anisotropy on the LIM power spectrum due to projection effects, as well as a narrowing and a shift towards higher intensities of the voxel intensity distribution. We forecast the minimum rate of decays into two photons that LIM surveys will be sensitive to as function of the dark-matter mass in the range $\sim 10^{-6}-10$ eV, and discuss how to reinterpret such results for dark matter that decays into a photon and another particle. We find that both the power spectrum and the voxel intensity distribution are expected to be very sensitive to the dark-matter contribution, with the voxel intensity distribution being more promising for most experiments considered. Interpreting our results in terms of the axion, we show that LIM surveys will be extremely competitive to detect its decay products, improving several orders of magnitudes (depending on the mass) the sensitivity of laboratory and astrophysical searches, especially in the mass range $\sim 1-10$ eV.


title: Probing Large-scale UV Background Inhomogeneity Associated with Quasars using Metal Absorption

abstract: We study large-scale UV background inhomogeneity in three-dimensions associated with the observed quasar populations at high redshift. We do this by measuring metal absorption through quasar absorption spectrum stacking as a function distance to closest quasar in SDSS-IV/eBOSS on 10s of comoving megaparsec scales. We study both intergalactic medium absorbers and mixed circumgalactic medium absorbers and probe absorption in O VI, C IV, and Si IV, and Si III. Overall stronger high ionization species absorption is seen closer to quasars at $2.4<z<3.1$. O VI absorption shows a particularly strong change, with effects in C IV evident in some cases, and more marginal effects in Si III and Si IV. We further study $2.05<z<2.4$ (with weak signs of increasing homogeneity with time) and explore the study of metal absorption as a function of integrated SDSS-r band flux quasar flux (yielding consistent but less significant results). While the metal absorption does show sensitivity to large-scale 3D quasar proximity, the current incomplete quasar samples limit detailed interpretation. This work does, however, demonstrate that UV background inhomogeneities exist on scales of several 10s of comoving megaparsecs associated with quasars and that they can be measured with precision by examining metal absorption in the intergalactic medium.


title: Higgs alignment and novel $CP$-Violating observables in two-Higgs-doublet models

abstract: Null results from searches for new physics at the Large Hadron Collider (LHC) tend to enforce the belief that new particles must be much heavier than the weak scale. We undertake a systematic study of the interplay between Higgs alignment and $CP$-violation in complex two-Higgs-doublet models, which enables us to construct a $CP$-violating scenario where new Higgs bosons are close to the weak scale after including stringent constraints from the electric dipole moment and measurements at the LHC. In addition, we propose a novel potential signal of $CP$-violation in the Higgs-to-Higgs decays, $h_3\to h_2 h_1$, where $h_3, h_2$, and $h_1$ are the heaviest, second heaviest and the Standard Model-like neutral Higgs bosons, respectively. The decay could manifest itself in triple boson final states in $h_1h_1h_1$ and $h_1h_1Z$, which are quite distinct and provide unique venues for new measurements at the LHC.


title: Electroweak Restoration at the LHC and Beyond: The $Vh$ Channel

abstract: The LHC is exploring electroweak (EW) physics at the scale EW symmetry is broken. As the LHC and new high energy colliders push our understanding of the Standard Model to ever-higher energies, it will be possible to probe not only the breaking of but also the restoration of EW symmetry. We propose to observe EW restoration in double EW boson production via the convergence of the Goldstone boson equivalence theorem. This convergence is most easily measured in the vector boson plus Higgs production, $Vh$, which is dominated by the longitudinal polarizations. We define EW restoration by carefully taking the limit of zero Higgs vacuum expectation value (vev). EW restoration is then measured through the ratio of the $p_T^h$ distributions between $Vh$ production in the Standard Model and Goldstone boson plus Higgs production in the zero vev theory, where $p_T^h$ is the Higgs transverse momentum. As EW symmetry is restored, this ratio converges to one at high energy. We present a method to extract this ratio from collider data. With a full signal and background analysis, we demonstrate that the 14 TeV HL-LHC can confirm that this ratio converges to one to 40% precision while at the 27 TeV HE-LHC the precision will be 6%. We also investigate statistical tests to quantify the convergence at high energies. Our analysis provides a roadmap for how to stress test the Goldstone boson equivalence theorem and our understanding of spontaneously broken symmetries, in addition to confirming the restoration of EW symmetry.


title: Morphological evolution of supermassive black hole merger hosts and multimessenger signatures

abstract: With projects such as Laser Interferometer Space Antenna (LISA) and Pulsar Timing Arrays expected to detect gravitational waves from supermassive black hole mergers in the near future, it is key that we understand what we expect those detections to be, and maximize what we can learn from them. To address this, we study the mergers of supermassive black holes in the Illustris simulation, the overall rate of mergers, and the correlation between merging black holes and their host galaxies. We find that these mergers occur in typical galaxies along the $M_{\rm{BH}}-M_*$ relation, and that between LISA and PTAs we expect to probe the full range of galaxy masses. As galaxy mergers can trigger increased star formation, we find that galaxies hosting low-mass black hole mergers tend to show a slight increase in star formation rates compared to a mass-matched sample. However, high-mass merger hosts have typical star formation rates, due to a combination of low gas fractions and powerful AGN feedback. Although minor black hole mergers do not correlate with disturbed morphologies, major mergers (especially at high-masses) tend to show morphological evidence of recent galaxy mergers which survives for ~500 Myr. This is on the same scale as the infall/hardening time of the merging black holes, suggesting that electromagnetic followups to gravitational wave signals may not be able to observe this correlation. We further find that incorporating a realistic timescale delay for the black hole mergers could shift the distribution of merger masses toward higher-masses, decreasing the rate of LISA detections while increasing the rate of PTA detections.


title: Hot Stars With Kepler Planets Have High Obliquities

abstract: It has been known for a decade that hot stars with hot Jupiters tend to have high obliquities. Less is known about the degree of spin-orbit alignment for hot stars with other kinds of planets. Here, we re-assess the obliquities of hot Kepler stars with transiting planets smaller than Neptune, based on spectroscopic measurements of their projected rotation velocities (vsini). The basis of the method is that a lower obliquity -- all other things being equal -- causes sini to be closer to unity and increases the value of vsini. We sought evidence for this effect using a sample of 150 Kepler stars with effective temperatures between 5950 and 6550K and a control sample of 101 stars with matching spectroscopic properties and random orientations. The planet hosts have systematically higher values of vsini than the control stars, but not by enough to be compatible with perfect spin-orbit alignment. The mean value of sini is 0.856 +/- 0.036, which is 4-sigma away from unity (perfect alignment), and 2-sigma away from pi/4 (random orientations). There is also evidence that the hottest stars have a broader obliquity distribution: when modeled separately, the stars cooler than 6250K have <sini> = 0.928 +/- 0.042, while the hotter stars are consistent with random orientations. This is similar to the pattern previously noted for stars with hot Jupiters. Based on these results, obliquity excitation for early-G and late-F stars appears to be a general outcome of star and planet formation, rather than being exclusively linked to hot Jupiter formation.


title: Characterizing many-body localization via exact disorder-averaged quantum noise

abstract: Many-body localized (MBL) phases of disordered quantum many-particle systems have a number of unique properties, including failure to act as a thermal bath and protection of quantum coherence. Studying MBL is complicated by the effects of rare ergodic regions, necessitating large system sizes and averaging over many disorder configurations. Here, building on the Feynman-Vernon theory of quantum baths, we characterize the quantum noise that a disordered spin system exerts on its parts via an influence matrix (IM). In this approach, disorder averaging is implemented exactly, and the thermodynamic-limit IM obeys a self-consistency equation. Viewed as a wavefunction in the space of trajectories of an individual spin, the IM exhibits slow scaling of temporal entanglement in the MBL phase. This enables efficient matrix product states computations to obtain temporal correlations, providing a benchmark for quantum simulations of non-equilibrium matter. The IM quantum noise formulation provides an alternative starting point for novel rigorous studies of MBL.


title: ARMADA I: Triple Companions Detected in B-Type Binaries alpha Del and nu Gem

abstract: Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrument at the CHARA array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of micro-arcseconds level for <0.2 arcsecond binaries. We test our instrument performance on a known triple system kappa Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30-day companion to alpha Del, and a 50-day companion to nu Gem. We also collected radial velocity data for alpha Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published RV orbit for the inner pair of nu Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history.


title: Dynamic Feature Pyramid Networks for Object Detection

abstract: Feature pyramid network (FPN) is a critical component in modern object detection frameworks. The performance gain in most of the existing FPN variants is mainly attributed to the increase of computational burden. An attempt to enhance the FPN is enriching the spatial information by expanding the receptive fields, which is promising to largely improve the detection accuracy. In this paper, we first investigate how expanding the receptive fields affect the accuracy and computational costs of FPN. We explore a baseline model called inception FPN in which each lateral connection contains convolution filters with different kernel sizes. Moreover, we point out that not all objects need such a complicated calculation and propose a new dynamic FPN (DyFPN). The output features of DyFPN will be calculated by using the adaptively selected branch according to a dynamic gating operation. Therefore, the proposed method can provide a more efficient dynamic inference for achieving a better trade-off between accuracy and computational cost. Extensive experiments conducted on MS-COCO benchmark demonstrate that the proposed DyFPN significantly improves performance with the optimal allocation of computation resources. For instance, replacing inception FPN with DyFPN reduces about 40% of its FLOPs while maintaining similar high performance.


title: Refining Deep Generative Models via Discriminator Gradient Flow

abstract: Deep generative modeling has seen impressive advances in recent years, to the point where it is now commonplace to see simulated samples (e.g., images) that closely resemble real-world data. However, generation quality is generally inconsistent for any given model and can vary dramatically between samples. We introduce Discriminator Gradient flow (DGflow), a new technique that improves generated samples via the gradient flow of entropy-regularized f-divergences between the real and the generated data distributions. The gradient flow takes the form of a non-linear Fokker-Plank equation, which can be easily simulated by sampling from the equivalent McKean-Vlasov process. By refining inferior samples, our technique avoids wasteful sample rejection used by previous methods (DRS & MH-GAN). Compared to existing works that focus on specific GAN variants, we show our refinement approach can be applied to GANs with vector-valued critics and even other deep generative models such as VAEs and Normalizing Flows. Empirical results on multiple synthetic, image, and text datasets demonstrate that DGflow leads to significant improvement in the quality of generated samples for a variety of generative models, outperforming the state-of-the-art Discriminator Optimal Transport (DOT) and Discriminator Driven Latent Sampling (DDLS) methods.


title: Pose-based Sign Language Recognition using GCN and BERT

abstract: Sign language recognition (SLR) plays a crucial role in bridging the communication gap between the hearing and vocally impaired community and the rest of the society. Word-level sign language recognition (WSLR) is the first important step towards understanding and interpreting sign language. However, recognizing signs from videos is a challenging task as the meaning of a word depends on a combination of subtle body motions, hand configurations, and other movements. Recent pose-based architectures for WSLR either model both the spatial and temporal dependencies among the poses in different frames simultaneously or only model the temporal information without fully utilizing the spatial information. We tackle the problem of WSLR using a novel pose-based approach, which captures spatial and temporal information separately and performs late fusion. Our proposed architecture explicitly captures the spatial interactions in the video using a Graph Convolutional Network (GCN). The temporal dependencies between the frames are captured using Bidirectional Encoder Representations from Transformers (BERT). Experimental results on WLASL, a standard word-level sign language recognition dataset show that our model significantly outperforms the state-of-the-art on pose-based methods by achieving an improvement in the prediction accuracy by up to 5%.


title: Analyticity of critical exponents of the $O(N)$ models from nonperturbative renormalization

abstract: We employ the functional renormalization group framework at the second order in the derivative expansion to study the $O(N)$ models continuously varying the number of field components $N$ and the spatial dimensionality $d$. We in particular address the Cardy-Hamber prediction concerning nonanalytical behavior of the critical exponents $\nu$ and $\eta$ across a line in the $(d,N)$ plane, which passes through the point $(2,2)$. By direct numerical evaluation of $\eta(d,N)$ and $\nu^{-1}(d,N)$ as well as analysis of the functional fixed-point profiles, we find clear indications of this line in the form of a crossover between two regimes in the $(d,N)$ plane, however no evidence of discontinuous or singular first and second derivatives of these functions for $d>2$. The computed derivatives of $\eta(d,N)$ and $\nu^{-1}(d,N)$ become increasingly large for $d\to 2$ and $N\to 2$ and it is only in this limit that $\eta(d,N)$ and $\nu^{-1}(d,N)$ as obtained by us are evidently nonanalytical. By scanning the dependence of the subleading eigenvalue of the RG transformation on $N$ for $d>2$ we find no indication of its vanishing as anticipated by the Cardy-Hamber scenario. For dimensionality $d$ approaching 3 there are no signatures of the Cardy-Hamber line even as a crossover and its existence in the form of a nonanalyticity of the anticipated form is excluded.


title: Interpretable Phase Detection and Classification with Persistent Homology

abstract: We apply persistent homology to the task of discovering and characterizing phase transitions, using lattice spin models from statistical physics for working examples. Persistence images provide a useful representation of the homological data for conducting statistical tasks. To identify the phase transitions, a simple logistic regression on these images is sufficient for the models we consider, and interpretable order parameters are then read from the weights of the regression. Magnetization, frustration and vortex-antivortex structure are identified as relevant features for characterizing phase transitions.


title: Comparison of cell- and vertex-centered finite-volume schemes for flow in fractured porous media

abstract: Flow in fractured porous media is of high relevance in a variety of geotechnical applications, given the fact that they ubiquitously occur in nature and that they can have a substantial impact on the hydraulic properties of rock. As a response to this, an active field of research has developed over the past decades, focussing on the development of mathematical models and numerical methods to describe and simulate the flow processes in fractured rock. In this work, we present a comparison of different finite volume-based numerical schemes for the simulation of flow in fractured porous media by means of numerical experiments. Two novel vertex-centered approaches are presented and compared to well-established numerical schemes in terms of convergence behaviour and their performance on benchmark studies taken from the literature. The new schemes show to produce results that are in good agreement with those of established methods while being computationally less expensive on unstructured simplex grids.


title: Clustering and ordering in cell assemblies with generic asymmetric aligning interactions

abstract: Collective cell migration plays an essential role in various biological processes, such as development or cancer proliferation. While cell-cell interactions are clearly key determinants of collective cell migration -- in addition to individual cells self-propulsion -- the physical mechanisms that control the emergence of cell clustering and collective cell migration are still poorly understood. In particular, observations have shown that binary cell-cell collisions generally lead to anti-alignement of cell polarities and separation of pairs -- a process called contact inhibition of locomotion (CIL), which is expected to disfavor the formation of large scale cell clusters with coherent motion. Here, we adopt a joint experimental and theoretical approach to determine the large scale dynamics of cell assemblies from elementary pairwise cell-cell interaction rules. We quantify experimentally binary cell-cell interactions and show that they can be captured by a minimal equilibrium-like pairwise asymmetric aligning interaction potential that reproduces the CIL phenomenology. We identify its symmetry class, build the corresponding active hydrodynamic theory and show on general grounds that such asymmetric aligning interaction destroys large scale clustering and ordering, leading instead to a liquid-like microphase of cell clusters of finite size and short lived polarity, or to a fully dispersed isotropic phase. Finally, this shows that CIL-like asymmetric interactions in cellular systems -- or general active systems -- control cluster sizes and polarity, and can prevent large scale coarsening and long range polarity, except in the singular regime of dense confluent systems.


title: Search for flares and associated CMEs on late-type main-sequence stars in optical SDSS spectra

abstract: This work aims to detect and classify stellar flares and potential stellar coronal mass ejection (CME) signatures in optical spectra provided by the Sloan Digital Sky Survey (SDSS) data release 14. The sample is constrained to all F, G, K, and M main-sequence type stars, resulting in more than 630,000 stars. This work makes use of the individual spectral exposures provided by the SDSS. An automatic flare search was performed by detecting significant amplitude changes in the $H\alpha$ and $H\beta$ spectral lines after a Gaussian profile was fit to the line core. CMEs were searched for by identifying asymmetries in the Balmer lines caused by the Doppler effect of plasma motions in the line of sight. We identified 281 flares on late-type stars (spectral types K3 to M9). We identified six possible CME candidates showing excess flux in Balmer line wings. Flare energies in $H\alpha$ were calculated and masses of the CME candidates were estimated. The derived $H\alpha$ flare energies range from $3 \times 10^{28} - 2 \times 10^{33}$ erg. The $H\alpha$ flare energy increases with earlier types, while the fraction of flaring times increases with later types. Mass estimates for the CME candidates are in the range of $6 \times 10^{16} - 6 \times 10^{18}$ g, and the highest projected velocities are $\sim300 - 700\:$ km/s. The low detection rate of CMEs we obtained agrees with previous studies, suggesting that for late-type main-sequence stars the CME occurrence rate that can be detected with optical spectroscopy is low.


title: Testable Implications of Multiple Equilibria in Discrete Games with Correlated Types

abstract: We study testable implications of multiple equilibria in discrete games with incomplete information. Unlike de Paula and Tang (2012), we allow the players' private signals to be correlated. In static games, we leverage independence of private types across games whose equilibrium selection is correlated. In dynamic games with serially correlated discrete unobserved heterogeneity, our testable implication builds on the fact that the distribution of a sequence of choices and states are mixtures over equilibria and unobserved heterogeneity. The number of mixture components is a known function of the length of the sequence as well as the cardinality of equilibria and unobserved heterogeneity support. In both static and dynamic cases, these testable implications are implementable using existing statistical tools.


title: Combining Hybrid Input-Output Automaton and Game Theory for Security Modeling of Cyber-Physical Systems

abstract: We consider a security setting in which the Cyber-Physical System (CPS) is composed of subnetworks where each subnetwork is under ownership of one defender. Such CPS can be represented by an attack graph where the defenders are required to invest (subject to a budget constraint) on the graph's edges in order to protect their critical assets (where each defender's critical asset has a certain value to the defender if compromised). We model such CPS using Hybrid Input-Output Automaton (HIOA) where each subnetwork is represented by a HIOA module. We first establish the building blocks needed in our setting. We then present our model that characterizes the continuous time evolution of the investments and discrete transitions between different states (where each state represents different condition and/or perturbation) within the system. Finally, we provide a real-world CPS example to validate our modeling.


title: Separating Mesoscale and Submesoscale Flows from Clustered Drifter Trajectories

abstract: Drifters deployed in close proximity collectively provide a unique observational data set with which to separate mesoscale and submesoscale flows. In this paper we provide a principled approach for doing so by fitting observed velocities to a local Taylor expansion of the velocity flow field. We demonstrate how to estimate mesoscale and submesoscale quantities that evolve slowly over time, as well as their associated statistical uncertainty. We show that in practice the mesoscale component of our model can explain much first and second-moment variability in drifter velocities, especially at low frequencies. This results in much lower and more meaningful measures of submesoscale diffusivity, which would otherwise be contaminated by unresolved mesoscale flow. We quantify these effects theoretically via computing Lagrangian frequency spectra, and demonstrate the usefulness of our methodology through simulations as well as with real observations from the LatMix deployment of drifters. The outcome of this method is a full Lagrangian decomposition of each drifter trajectory into three components that represent the background, mesoscale, and submesoscale flow.

categories: physics.ao-ph eess.SP physics.data-an physics.flu-dyn stat.AP


title: Wrapped Fivebranes Redux

abstract: The near-horizon region of Neveu-Schwarz fivebranes provides interesting examples of gauge/gravity duality. We revisit the structure of wrapped and/or intersecting fivebranes using the tools of null-gauged WZW models in worldsheet string theory, revealing the effective geometry of the fivebrane throat in a variety of examples. Variant gaugings yield linear dilaton fivebrane throats with AdS3 caps, providing a wealth of information about the near-BPS structure of the corresponding spacetime CFT duals.


title: Structure and Internal Kinematics of Nine Inner Milky Way Globular Clusters

abstract: This study constitutes part of a larger effort aimed at better characterizing the Galactic globular clusters (GGCs) located towards the inner Milky Way bulge and disk. Here, we focus on internal kinematics of nine GGCs, obtained from space-based imaging over time baselines of $>$9 years. We exploit multiple avenues to assess the dynamical state of the target GGCs, constructing radial profiles of projected stellar density, proper motion dispersion, and anisotropy. We posit that two-thirds (6/9) of our target GGCs are in an advanced dynamical state, and are close to (or have recently undergone) core collapse, supported by at least two lines of evidence: First, we find relatively steep proper motion dispersion profiles, in accord with literature values for core-collapsed GGCs. Second, we find that our sample is, in the mean, isotropic even out to their half-light radii, although one of our target clusters (NGC 6380) is tangentially anisotropic at $>$1$\sigma$ beyond its half-light radius, in accord with theoretical predictions for clusters evolving in strong tidal fields. Our proper motion dispersion and anisotropy profiles are made publicly available.


title: GPI 2.0: Upgrades to the IFS including new spectral modes

abstract: The Gemini Planet Imager (GPI) is a high-contrast imaging instrument designed to directly image and characterize exoplanets. GPI is currently undergoing several upgrades to improve performance. In this paper, we discuss the upgrades to the GPI IFS. This primarily focuses on the design and performance improvements of new prisms and filters. This includes an improved high-resolution prism which will provide more evenly dispersed spectra across y, J, H and K-bands. Additionally, we discuss the design and implementation of a new low-resolution mode and prism which allow for imaging of all four bands (y, J, H and K-bands) simultaneously at R=10. We explore the possibility of using a multiband filter which would block the light between the four spectral bands. We discuss possible performance improvements from the multiband filter, if implemented. Finally we explore the possibility of making small changes to the optical design to improve the IFS's performance near the edge of the field of view.


title: The delay time distribution of Type-Ia supernovae in galaxy clusters: the impact of extended star-formation histories

abstract: The delay time distribution (DTD) of Type-Ia supernovae (SNe Ia) is important for understanding chemical evolution, SN Ia progenitors, and SN Ia physics. Past estimates of the DTD in galaxy clusters have been deduced from SN Ia rates measured in cluster samples observed at various redshifts, corresponding to different time intervals after a presumed initial brief burst of star formation. A recent analysis of a cluster sample at $z=1.13-1.75$ confirmed indications from previous studies of lower-redshift clusters, that the DTD has a power-law form, ${\rm DTD}(t)=R_1 (t/{\rm Gyr})^\alpha$, with amplitude $R_1$, at delay $t=1~\rm Gyr$, several times higher than measured in field-galaxy environments. This implied that SNe Ia are somehow produced in larger numbers by the stellar populations in clusters. This conclusion, however, could have been affected by the implicit assumption that the stars were formed in a single brief starburst at high $z$. Here, we re-derive the DTD from the cluster SN Ia data, but relax the single-burst assumption. Instead, we allow for a range of star-formation histories and dust extinctions for each cluster. Via MCMC modeling, we simultaneously fit, using stellar population synthesis models and DTD models, the integrated galaxy-light photometry in several bands, and the SN Ia numbers discovered in each cluster. With these more-realistic assumptions, we find a best-fit DTD with power-law index $\alpha=-1.09_{-0.12}^{+0.15}$, and amplitude $R_1=0.41_{-0.10}^{+0.12}\times 10^{-12}~{\rm yr}^{-1}{\rm M}_\odot^{-1}$. We confirm a cluster-environment DTD with a larger amplitude than the field-galaxy DTD, by a factor $\sim2-3$ (at $3.8\sigma$). Cluster and field DTDs have consistent slopes of $\alpha\approx-1.1$.


title: Using math in physics: 5. Functional dependence

abstract: When students are learning to use math in physics, one of the most important ideas they need to learn is that equations are not just calculational tools; they represent relationships between physical variables that change together (covary). How much a change in one variable or parameter is associated with a change in another depends on how they appear in the equation: their functional dependence. Understanding this sort of relationship is rarely taught in introductory mathematics classes, and students who have not yet learned to blend conceptual ideas with mathematical symbols may not see the relevance and power of this idea. We need to explicitly teach functional dependence as part of our effort to help students to learn to use math productively in science.


title: The Exoplanet Transmission Spectroscopy Imager (ETSI)

abstract: We present the design of a novel instrument tuned to detect transiting exoplanet atmospheres. The instrument, which we call the exoplanet transmission spectroscopy imager (ETSI), makes use of a new technique called common-path multi-band imaging (CMI). ETSI uses a prism and multi-band filter to simultaneously image 15 spectral bandpasses on two detectors from $430-975nm$ (with a average spectral resolution of $R = \lambda/\Delta\lambda = 23$) during exoplanet transits of a bright star. A prototype of the instrument achieved photon-noise limited results which were below the atmospheric amplitude scintillation noise limit. ETSI can detect the presence and composition of an exoplanet atmosphere in a relatively short time on a modest-size telescope. We show the optical design of the instrument. Further, we discuss design trades of the prism and multi-band filter which are driven by the science of the ETSI instrument. We describe the upcoming survey with ETSI that will measure dozens of exoplanet atmosphere spectra in $\sim2$ years on a two meter telescope. Finally, we will discuss how ETSI will be a powerful means for follow up on all gas giant exoplanets that transit bright stars, including a multitude of recently identified TESS (NASA's Transiting Exoplanet Survey Satellite) exoplanets.


title: Wireless Secret Sharing Game between Two Legitimate Users and an Eavesdropper

abstract: Wireless secret sharing is crucial to information security in the era of Internet of Things. One method is to utilize the effect of the randomness of the wireless channel in the data link layer to generate the common secret between two legitimate users Alice and Bob. This paper studies this secret sharing mechanism from the perspective of game theory. In particular, we formulate a non-cooperative zero-sum game between the legitimate users and an eavesdropper Eve. In a symmetrical game where Eve has the same probability of successfully receiving a packet from Alice and Bob when the transmission distance is the same, we show that both pure and mixed strategy Nash equilibria exist. In an asymmetric game where Eve has different probabilities of successfully receiving a packet from Alice and Bob, a pure strategy may not exist; in this case, we show how a mixed strategy Nash equilibrium can be found.


title: Searching for odderon exchange in exclusive $pp \to pp \phi$ and $pp \to pp \phi \phi$ reactions at the LHC

abstract: The possibility to use the exclusive $pp \to pp \phi$ and $pp \to pp \phi \phi$ reactions in identifying the odderon exchange is discussed. So far there is no unambiguous experimental evidence for the odderon, the charge conjugation $C = -1$ counterpart of the $C = +1$ pomeron. Recently proposed tensor-pomeron and vector-odderon model for soft high-energy reactions is applied. For the $p p \to p p \phi$ reaction at high energies the photon-pomeron fusion is the dominant process and the odderon-pomeron fusion is an interesting alternative. Adding odderon exchange improves considerably description of the proton-proton angular correlations measured by the WA102 collaboration. The $p p \to p p \phi \phi$ process via pomeron-pomeron fusion is advantageous one as here the odderon does not couple to protons. The observation of large $M_{\phi \phi}$ and the rapidity difference $Y_{\phi \phi}$ seems well suited to identify odderon exchange. Comparisons with data from the WA102 experiment are made and predictions for LHC experiments are given.


title: Time-Delayed Generalized BSDEs

abstract: We prove the existence and uniqueness of the solution of a BSDE with time-delayed generators in the small delay setting (or equivalently small Lipschitz constant), which employs the Stieltjes integral with respect to an increasing continuous stochastic process. Moreover, we obtain a result of continuity of the solution with regard to the increasing process, assuming only uniform convergence, but not in variation. We also prove the existence in the case of an arbitrary delay by imposing monotonicity and linearity on generators. Lastly, we provide an application of the theoretical framework within an insurance based example.


title: On the Regularity of Optimal Dynamic Blocking Strategies

abstract: The paper studies a dynamic blocking problem, motivated by a model of optimal fire confinement. While the fire can expand with unit speed in all directions, barriers are constructed in real time. An optimal strategy is sought, minimizing the total value of the burned region, plus a construction cost. It is well known that optimal barriers exists. In general, they are a countable union of compact, connected, rectifiable sets. The main result of the present paper shows that optimal barriers are nowhere dense. The proof relies on new estimates on the reachable sets and on optimal trajectories for the fire, solving a minimum time problem in the presence of obstacles.


title: The physical meaning of the magnetic scalar potential, and how to use it to design an electromagnet

abstract: The magnetic scalar potential can be used to design electromagnets accurately and efficiently. I will describe a practical construction algorithm: the prescribed field in a "Target" region (constrained only by Maxwell's equations) specifies boundary conditions which uniquely determine the potential in the surrounding field "Return" region. As required by the physical representation of the magnetic scalar potential, the conducting elements of the resulting coil are directed along equipotentials on the surface of each region, at equal increments of the potential. I give some examples and comments on the limits of precision of the constructed field.


title: Seven steps to reliable cyclic voltammetry measurements for the determination of double layer capacitance

abstract: Discovery of electrocatalytic materials for high-performance energy conversion and storage applications relies on the adequate characterization of their intrinsic activity, which is currently hindered by the dearth of a protocol for consistent and precise determination of double layer capacitance (CDL). Herein, we propose a seven-step method that aims to determine CDL reliably by scan rate-dependent cyclic voltammetry. The method considers three aspects that strongly influence the outcome of the analysis: measurement settings, data collection, and data processing. To illustrate the proposed method, two systems were studied: a resistor-capacitor electric circuit and a glassy carbon disk in an electrochemical cell. With these studies it is demonstrated that when any of the mentioned aspects of the procedure are neglected, substantial deviations of the results are observed with misestimations as large as 61% in the case of the investigated electrochemical system. Moreover, we propose allometric regression as a more suitable model than linear regression for the determination of CDL for both the ideal and the non-ideal systems investigated. We stress the importance of assessing the accuracy of not only highly specialized electrochemical methods, but also of those that are well-known and commonly used as it is the case of the voltammetric methods. The methodology proposed herein is not limited to the determination of CDL, but can be effectively applied to any other voltammetry-based analysis that aims to deliver quantitative results.


title: Adversarial Robustness Across Representation Spaces

abstract: Adversarial robustness corresponds to the susceptibility of deep neural networks to imperceptible perturbations made at test time. In the context of image tasks, many algorithms have been proposed to make neural networks robust to adversarial perturbations made to the input pixels. These perturbations are typically measured in an $\ell_p$ norm. However, robustness often holds only for the specific attack used for training. In this work we extend the above setting to consider the problem of training of deep neural networks that can be made simultaneously robust to perturbations applied in multiple natural representation spaces. For the case of image data, examples include the standard pixel representation as well as the representation in the discrete cosine transform~(DCT) basis. We design a theoretically sound algorithm with formal guarantees for the above problem. Furthermore, our guarantees also hold when the goal is to require robustness with respect to multiple $\ell_p$ norm based attacks. We then derive an efficient practical implementation and demonstrate the effectiveness of our approach on standard datasets for image classification.


title: Generator Parameter Estimation by Q-Learning Based on PMU Measurements

abstract: In this paper, a novel Q-learning based approach is proposed for estimating the parameters of synchronous generators using PMU measurements. Event playback is used to generate model outputs under different parameters for training the agent in Q-learning. We assume that the exact values of some parameters in the model are not known by the agent in Q-learning. Then, an optimal history-dependent policy for the exploration-exploitation trade-off is planned. With given prior knowledge, the parameter vector can be viewed as states with a specific reward, which is a function of the fitting error compared with the measurements. The agent takes an action (either increasing or decreasing the parameter) and the estimated parameter will move to a new state. Based on the reward function, the optimal action policy will move the parameter set to a state with the highest reward. If multiple events are available, they will be used sequentially so that the updated $\mathbfcal{Q}$-value can be utilized to improve the computational efficiency. The effectiveness of the proposed approach is validated by estimating the parameters of the dynamic model of a synchronous generator.


title: Reduced models of cardiomyocytes excitability: comparing Karma and FitzHugh-Nagumo

abstract: Since Noble adapted in 1962 the model of Hodgkin and Huxley to fit Purkinje fibres the refinement of models for cardiomyocytes has continued. Most of these models are high-dimensional systems of coupled equations so that the possible mathematical analysis is quite limited, even numerically. This has inspired the development of reduced, phenomenological models that preserve qualitatively the main featuture of cardiomyocyte's dynamics. In this paper we present a systematic comparison of the dynamics between two notable low-dimensional models, the FitzHugh-Nagumo model \cite{FitzHugh55, FitzHugh60, FitzHugh61} as a prototype of excitable behaviour and a polynomial version of the Karma model \cite{Karma93, Karma94} which is specifically developed to fit cardiomyocyte's behaviour well. We start by introducing the models and considering their pure ODE versions. We analyse the ODEs employing the main ideas and steps used in the setting of geometric singular perturbation theory. Next, we turn to the spatially extended models, where we focus on travelling wave solutions in 1D. Finally, we perform numerical simulations of the 1D PDE Karma model varying model parameters in order to systematically investigate the impact on wave propagation velocity and shape. In summary, our study provides a reference regarding key similarities as well as key differences of the two models.


title: Stochastic Approximation with Markov Noise: Analysis and applications in reinforcement learning

abstract: We present for the first time an asymptotic convergence analysis of two time-scale stochastic approximation driven by "controlled" Markov noise. In particular, the faster and slower recursions have non-additive controlled Markov noise components in addition to martingale difference noise. We analyze the asymptotic behavior of our framework by relating it to limiting differential inclusions in both time scales that are defined in terms of the ergodic occupation measures associated with the controlled Markov processes. Using a special case of our results, we present a solution to the off-policy convergence problem for temporal-difference learning with linear function approximation. We compile several aspects of the dynamics of stochastic approximation algorithms with Markov iterate-dependent noise when the iterates are not known to be stable beforehand. We achieve the same by extending the lock-in probability (i.e. the probability of convergence to a specific attractor of the limiting o.d.e. given that the iterates are in its domain of attraction after a sufficiently large number of iterations (say) n_0) framework to such recursions. We use these results to prove almost sure convergence of the iterates to the specified attractor when the iterates satisfy an "asymptotic tightness" condition. This, in turn, is shown to be useful in analyzing the tracking ability of general "adaptive" algorithms. Finally, we obtain the first informative error bounds on function approximation for the policy evaluation algorithm proposed by Basu et al. when the aim is to find the risk-sensitive cost represented using exponential utility. We show that this happens due to the absence of difference term in the earlier bound which is always present in all our bounds when the state space is large.


title: Some applications of Wagner's weighted subgraph counting polynomial

abstract: We use Wagner's weighted subgraph counting polynomial to prove that the partition function of the anti-ferromagnetic Ising model on line graphs is real rooted and to prove that roots of the edge cover polynomial have length at most $4$. We moreover discuss how our results relate to efficient algorithms for approximately computing evaluations of these polynomials.


title: On the robustness of minimum norm interpolators and regularized empirical risk minimizers

abstract: This article develops a general theory for minimum norm interpolating estimators and regularized empirical risk minimizers (RERM) in linear models in the presence of additive, potentially adversarial, errors. In particular, no conditions on the errors are imposed. A quantitative bound for the prediction error is given, relating it to the Rademacher complexity of the covariates, the norm of the minimum norm interpolator of the errors and the size of the subdifferential around the true parameter. The general theory is illustrated for Gaussian features and several norms: The $\ell_1$, $\ell_2$, group Lasso and nuclear norms. In case of sparsity or low-rank inducing norms, minimum norm interpolators and RERM yield a prediction error of the order of the average noise level, provided that the overparameterization is at least a logarithmic factor larger than the number of samples and that, in case of RERM, the regularization parameter is small enough. Lower bounds that show near optimality of the results complement the analysis.

categories: math.ST cs.IT cs.NA math.IT math.NA stat.ML stat.TH


title: On the Laplacian spectra of token graphs

abstract: We study the Laplacian spectrum of token graphs, also called symmetric powers of graphs. The $k$-token graph $F_k(G)$ of a graph $G$ is the graph whose vertices are the $k$-subsets of vertices from $G$, two of which being adjacent whenever their symmetric difference is a pair of adjacent vertices in $G$. In this paper, we give a relationship between the Laplacian spectra of any two token graphs of a given graph. In particular, we show that, for any integers $h$ and $k$ such that $1\le h\le k\le \frac{n}{2}$, the Laplacian spectrum of $F_h(G)$ is contained in the Laplacian spectrum of $F_k(G)$. We also show that the double odd graphs and doubled Johnson graphs can be obtained as token graphs of the complete graph $K_n$ and the star $S_{n}=K_{1,n-1}$, respectively. Besides, we obtain a relationship between the spectra of the $k$-token graph of $G$ and the $k$-token graph of its complement $\overline{G}$. This generalizes a well-known property for Laplacian eigenvalues of graphs to token graphs. Finally, the double odd graphs and doubled Johnson graphs provide two infinite families, together with some others, in which the algebraic connectivities of the original graph and its token graph coincide. Moreover, we conjecture that this is the case for any graph $G$ and its token graph.


title: TARdYS: Design and Prototype of an Exoplanet Hunter for TAO using a R6 Echelle Grating

abstract: One limitation in characterizing exoplanet candidates is the availability of infrared, high-resolution spectrographs. An important factor in the scarcity of high precision IR spectrographs is the high cost of these instruments. We present a new optical design, which leads to a cost-effective solution. Our instrument is a high-resolution (R=60,000) infrared spectrograph with a R6 Echelle grating and an image slicer. We compare the best possible performance of quasi-Littrow and White Pupil setups, and prefer the latter because it achieves higher image quality. The instrument is proposed for the University of Tokyo Atacama Observatory (TAO) 6.5 m telescope in Chile. The Tao Aiuc high Resolution (d) Y band Spectrograph (TARdYS) covers 0.843-1.117 um. To reduce the cost, we squeeze 42 spectral orders onto a 1K detector with a semi-cryogenic solution. We obtain excellent resolution even when taking realistic manufacturing and alignment tolerances as well as thermal variations into account. In this paper, we present early results from the prototype of this spectrograph at ambient temperature.


title: Population genetic models of dormancy

abstract: In the present article, we investigate the effects of dormancy on an abstract population genetic level. We first provide a short review of seed bank models in population genetics, and the role of dormancy for the interplay of evolutionary forces in general, before we discuss two recent paradigmatic models, referring to spontaneous resp. simultaneous switching of individuals between the active and the dormant state. We show that both mechanisms give rise to non-trivial mathematical objects, namely the (continuous) seed bank diffusion and the seed bank diffusion with jumps, as well as their dual processes, the seed bank coalescent and the seed bank coalescent with simultaneous switching.


title: Three-dimensional simulations of the inhomogeneous Low Solar Wind

abstract: In the near future, Parker Solar Probe will put theories about the dynamics and nature of the transition between the solar corona and the solar wind to stringent tests. The most popular mechanism aimed to explain the dynamics of the nascent solar wind, including its heating and acceleration is magnetohydrodynamic (MHD) turbulence. Most of the previous models focus on nonlinear cascade induced by interactions of outgoing Alfv\'en waves and their reflections, ignoring effects that might be related to perpendicular structuring of the solar coronal plasma, despite overwhelming evidence for it. In this paper, for the first time, we analyse through 3D MHD numerical simulations the dynamics of the perpendicularly structured solar corona and solar wind, from the low corona to 15 R_Sun. We find that background structuring has a strong effect on the evolution of MHD turbulence, on much faster time scales than in the perpendicularly homogeneous case. On time scales shorter than nonlinear times, linear effects related to phase mixing result in a 1/f perpendicular energy spectrum. As the turbulent cascade develops, we observe a perpendicular (parallel) energy spectrum with the power law index of -3/2 or -5/3 (-2), a steeper perpendicular magnetic field than velocity spectrum, and a strong build-up of negative residual energy. We conclude that the turbulence is most probably generated by the self-cascade of the driven transverse kink waves, referred to previously as `uniturbulence', which might represent the dominant nonlinear energy cascade channel in the pristine solar wind.


title: Some game theoretic marketing attribution models

abstract: In this paper, we propose and analyse two game theoretical models useful to design marketing channels attribution mechanisms based on cooperative TU games and bankruptcy problems, respectively. First, we analyse the Sum Game, a coalitional game introduced by Morales (2016). We extend the ideas introduced in Zhao et al. (2018) and Cano-Berlanga et al. (2017) to the case in which the order and the repetition of channels on the paths to conversion are taken into account. In all studied cases, the Shapley value is proposed as the attribution mechanism. Second, a bankruptcy problem approach is proposed, and a similar analysis is developed relying on the Constrained Equal Loss (CEL) and Proportional (PROP) rules as attribution mechanisms. In particular, it is relevant to note that the class of attribution bankruptcy problems is a proper subclass of bankruptcy problems.


title: Subsums of Null Sequences: A Survey and Open Questions

abstract: We give a survey on the different results involving the topological structure of subsums of null sequences.


title: Nonstationary generalized TASEP in KPZ and jamming regimes

abstract: We study the model of the totally asymmetric exclusion process with generalized update, which compared to the usual totally asymmetric exclusion process, has an additional parameter enhancing clustering of particles. We derive the exact multiparticle distributions of distances travelled by particles on the infinite lattice for two types of initial conditions: step and alternating once. Two different scaling limits of the exact formulas are studied. Under the first scaling associated to Kardar-Parisi-Zhang (KPZ) universality class we prove convergence of joint distributions of the scaled particle positions to finite-dimensional distributions of the universal Airy$_2$ and Airy$_1$ processes. Under the second scaling we prove convergence of the same position distributions to finite-dimensional distributions of two new random processes, which describe the transition between the KPZ regime and the deterministic aggregation regime, in which the particles stick together into a single giant cluster moving as one particle. It is shown that the transitional distributions have the Airy processes and fully correlated Gaussian fluctuations as limiting cases. We also give the heuristic arguments explaining how the non-universal scaling constants appearing from the asymptotic analysis in the KPZ regime are related to the properties of translationally invariant stationary states in the infinite system and how the parameters of the model should scale in the transitional regime.

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


title: Subspace method for multiparameter-eigenvalue problems based on tensor-train representations

abstract: In this paper we solve $m$-parameter eigenvalue problems ($m$EPs), with $m$ any natural number by representing the problem using Tensor-Trains (TT) and designing a method based on this format. $m$EPs typically arise when separation of variables is applied to separable boundary value problems. Often, methods for solving $m$EP are restricted to $m = 3$, due to the fact that, to the best of our knowledge, no available solvers exist for $m>3$ and reasonable size of the involved matrices. In this paper, we prove that computing the eigenvalues of a $m$EP can be recast into computing the eigenvalues of TT-operators. We adapted the algorithm in \cite{Dolgov2014a} for symmetric eigenvalue problems in TT-format to an algorithm for solving generic $m$EPs. This leads to a subspace method whose subspace dimension does not depend on $m$, in contrast to other subspace methods for $m$EPS. This allows us to tackle $m$EPs with $m > 3$ and reasonable size of the matrices. We provide theoretical results and report numerical experiments. The MATLAB code is publicly available.


title: The generating function of Kreweras walks with interacting boundaries is not algebraic

abstract: Beaton, Owczarek and Xu (2019) studied generating functions of Kreweras walks and of reverse Kreweras walks in the quarter plane, with interacting boundaries. They proved that for the reverse Kreweras step set, the generating function is always algebraic, and for the Kreweras step set, the generating function is always D-finite. However, apart from the particular case where the interactions are symmetric in $x$ and~$y$, they left open the question of whether the latter one is algebraic. Using computer algebra tools, we confirm their intuition that the generating function of Kreweras walks is not algebraic, apart from the particular case already identified.


title: Game-Theoretic Malware Detection

abstract: Malware attacks are costly. To mitigate against such attacks, organizations deploy malware detection tools that help them detect and eventually resolve those threats. While running only the best available tool does not provide enough coverage of the potential attacks, running all available tools is prohibitively expensive in terms of financial cost and computing resources. Therefore, an organization typically runs a set of tools that maximizes their coverage given a limited budget. However, how should an organization choose that set? Attackers are strategic, and will change their behavior to preferentially exploit the gaps left by a deterministic choice of tools. To avoid leaving such easily-exploited gaps, the defender must choose a random set. In this paper, we present an approach to compute an optimal randomization over size-bounded sets of available security analysis tools by modeling the relationship between attackers and security analysts as a leader-follower Stackelberg security game. We estimate the parameters of our model by combining the information from the VirusTotal dataset with the more detailed reports from the National Vulnerability Database. In an empirical comparison, our approach outperforms a set of natural baselines under a wide range of assumptions.


title: Designing Voice-Controllable APIs

abstract: The main purpose of a voice command system is to process a sentence in natural language and perform the corresponding action. Although there exist many approaches to map sentences to API (application programming interface) calls, this mapping is usually performed after the API is already implemented, possibly by other programmers. In this paper, we describe how the API developer can use patterns to map sentences to API calls by utilizing the similarities between names and types in the sentences and the API. In the cases when the mapping is not straightforward, we suggest the usage of suitable annotations (attribute-oriented programming).


title: Carbonate-Silicate Cycle Predictions of Earth-like Planetary Climates and Testing the Habitable Zone Concept

abstract: In the conventional habitable zone (HZ) concept, a CO$_{2}$-H$_2$O greenhouse maintains surface liquid water. Through the water-mediated carbonate-silicate weathering cycle, atmospheric CO$_{2}$ partial pressure (pCO$_{2}$) responds to changes in surface temperature, stabilizing the climate over geologic timescales. We show that this weathering feedback ought to produce a log-linear relationship between pCO$_{2}$ and incident flux on Earth-like planets in the HZ. However, this trend has scatter because geophysical and physicochemical parameters can vary, such as land area for weathering and CO$_2$ outgassing fluxes. Using a coupled climate and carbonate-silicate weathering model, we quantify the likely scatter in pCO$_2$ with orbital distance throughout the HZ. From this dispersion, we predict a two-dimensional relationship between incident flux and pCO$_2$ in the HZ and show that it could be detected from at least 83 ($2{\sigma}$) Earth-like exoplanet observations. If fewer Earth-like exoplanets are observed, testing the HZ hypothesis from this relationship could be difficult.


title: Tensor decomposition for bosonic and fermionic scattering amplitudes

abstract: In this paper, we elaborate on a method to decompose multiloop multileg scattering amplitudes into Lorentz-invariant form factors, which exploits the simplifications that arise from considering four-dimensional external states. We propose a simple and general approach that applies to both fermionic and bosonic amplitudes and allows us to identify a minimal number of physically relevant form factors, which can be related one-to-one to the independent helicity amplitudes. We discuss explicitly its applicability to various four- and five-point scattering amplitudes relevant for LHC physics.


title: Automated Creation of a High-Performing Algorithmic Trader via Deep Learning on Level-2 Limit Order Book Data

abstract: We present results demonstrating that an appropriately configured deep learning neural network (DLNN) can automatically learn to be a high-performing algorithmic trading system, operating purely from training-data inputs generated by passive observation of an existing successful trader T. That is, we can point our black-box DLNN system at trader T and successfully have it learn from T's trading activity, such that it trades at least as well as T. Our system, called DeepTrader, takes inputs derived from Level-2 market data, i.e. the market's Limit Order Book (LOB) or Ladder for a tradeable asset. Unusually, DeepTrader makes no explicit prediction of future prices. Instead, we train it purely on input-output pairs where in each pair the input is a snapshot S of Level-2 LOB data taken at the time when T issued a quote Q (i.e. a bid or an ask order) to the market; and DeepTrader's desired output is to produce Q when it is shown S. That is, we train our DLNN by showing it the LOB data S that T saw at the time when T issued quote Q, and in doing so our system comes to behave like T, acting as an algorithmic trader issuing specific quotes in response to specific LOB conditions. We train DeepTrader on large numbers of these S/Q snapshot/quote pairs, and then test it in a variety of market scenarios, evaluating it against other algorithmic trading systems in the public-domain literature, including two that have repeatedly been shown to outperform human traders. Our results demonstrate that DeepTrader learns to match or outperform such existing algorithmic trading systems. We analyse the successful DeepTrader network to identify what features it is relying on, and which features can be ignored. We propose that our methods can in principle create an explainable copy of an arbitrary trader T via "black-box" deep learning methods.


title: Open-Ended Multi-Modal Relational Reasoning for Video Question Answering

abstract: In this paper, we introduce a robotic agent specifically designed to analyze external environments and address participants' questions. The primary focus of this agent is to assist individuals using language-based interactions within video-based scenes. Our proposed method integrates video recognition technology and natural language processing models within the robotic agent. We investigate the crucial factors affecting human-robot interactions by examining pertinent issues arising between participants and robot agents. Methodologically, our experimental findings reveal a positive relationship between trust and interaction efficiency. Furthermore, our model demonstrates a 2\% to 3\% performance enhancement in comparison to other benchmark methods.


title: Multiscale modeling of the effective viscoplastic behavior of Mg$_2$SiO$_4$ wadsleyite: Bridging atomic and polycrystal scales

abstract: The viscoplastic behavior of polycrystalline Mg$_2$SiO$_4$ wadsleyite aggregates, a major high pressure phase of the mantle transition zone of the Earth (depth range: 410 -- 520 km), is obtained by properly bridging several scale transition models. At the very fine nanometric scale corresponding to the disloca-tion core structure, the behavior of thermally activated plastic slip is modeled for strain-rates relevant for laboratory experimental conditions, at high pressure and for a wide range of temperatures, based on the Peierls-Nabarro-Galerkin model. Corresponding single slip reference resolved shear stresses and associated constitutive equations are deduced from Orowan's equation in order to describe the average viscoplastic behavior at the grain scale, for the easiest slip systems. These data have been implemented in two grain-polycrystal scale transition models, a mean-field one (the recent Fully-Optimized Second-Order Viscoplastic Self-Consistent scheme of [44]) allowing rapid evaluation of the effective viscosity of polycrystalline aggregates , and a full-field (FFT based [45] [33]) method allowing investigating stress and strain-rate localization in typical microstructures and heterogeneous activation of slip systems within grains. Calculations have been performed at pressure and temperatures relevant for in-situ conditions. Results are in very good agreement with available mechanical tests conducted at strain-rates typical for laboratory experiments.


title: Quantum-Inspired Classical Algorithm for Slow Feature Analysis

abstract: Recently, there has been a surge of interest for quantum computation for its ability to exponentially speed up algorithms, including machine learning algorithms. However, Tang suggested that the exponential speed up can also be done on a classical computer. In this paper, we proposed an algorithm for slow feature analysis, a machine learning algorithm that extracts the slow-varying features, with a run time O(polylog(n)poly(d)). To achieve this, we assumed necessary preprocessing of the input data as well as the existence of a data structure supporting a particular sampling scheme. The analysis of algorithm borrowed results from matrix perturbation theory, which was crucial for the algorithm's correctness. This work demonstrates the possible application and extent for which quantum-inspired computation can be used.


title: A Study of Checkpointing in Large Scale Training of Deep Neural Networks

abstract: Deep learning (DL) applications are increasingly being deployed on HPC systems, to leverage the massive parallelism and computing power of those systems for DL model training. While significant effort has been put to facilitate distributed training by DL frameworks, fault tolerance has been largely ignored. In this work, we evaluate checkpoint-restart, a common fault tolerance technique in HPC workloads. We perform experiments with three state-of-the-art DL frameworks common in HPC Chainer, PyTorch, and TensorFlow). We evaluate the computational cost of checkpointing, file formats and file sizes, the impact of scale, and deterministic checkpointing. Our evaluation shows some critical differences in checkpoint mechanisms and exposes several bottlenecks in existing checkpointing implementations. We provide discussion points that can aid users in selecting a fault-tolerant framework to use in HPC. We also provide takeaway points that framework developers can use to facilitate better checkpointing of DL workloads in HPC.


title: A Chatbot for Information Security

abstract: Advancements in artificial intelligence (AI), speech recognition systems (ASR), and machine learning have enabled the development of intelligent computer programs called chatbots. Many chatbots have been proposed to provide different services in many areas such as customer service, sales and marketing. However, the use of chatbot as advisers in the field of information security is not yet considered. Furthermore, people, especially normal users who have no technical background, are unaware about many of aspects in information security. Therefore, in this paper we proposed a chatbot that acts as an adviser in information security. The proposed adviser uses a knowledge base with json file. Having such chatbot provides many features including raising the awareness in field of information security by offering accurate advice, based on different opinions from information security expertise, for many users on different. Furthermore, this chatbot is currently deployed through Telegram platform, which is one of widely used social network platforms. The deployment of the proposed chatbot over different platforms is considered as the future work.


title: A Three-Stage Self-Training Framework for Semi-Supervised Semantic Segmentation

abstract: Semantic segmentation has been widely investigated in the community, in which the state of the art techniques are based on supervised models. Those models have reported unprecedented performance at the cost of requiring a large set of high quality segmentation masks. To obtain such annotations is highly expensive and time consuming, in particular, in semantic segmentation where pixel-level annotations are required. In this work, we address this problem by proposing a holistic solution framed as a three-stage self-training framework for semi-supervised semantic segmentation. The key idea of our technique is the extraction of the pseudo-masks statistical information to decrease uncertainty in the predicted probability whilst enforcing segmentation consistency in a multi-task fashion. We achieve this through a three-stage solution. Firstly, we train a segmentation network to produce rough pseudo-masks which predicted probability is highly uncertain. Secondly, we then decrease the uncertainty of the pseudo-masks using a multi-task model that enforces consistency whilst exploiting the rich statistical information of the data. We compare our approach with existing methods for semi-supervised semantic segmentation and demonstrate its state-of-the-art performance with extensive experiments.


title: Real-time gravitational replicas: Formalism and a variational principle

abstract: This work is the first step in a two-part investigation of real-time replica wormholes. Here we study the associated real-time gravitational path integral and construct the variational principle that will define its saddle-points. We also describe the general structure of the resulting real-time replica wormhole saddles, setting the stage for construction of explicit examples. These saddles necessarily involve complex metrics, and thus are accessed by deforming the original real contour of integration. However, the construction of these saddles need not rely on analytic continuation, and our formulation can be used even in the presence of non-analytic boundary-sources. Furthermore, at least for replica- and CPT-symmetric saddles we show that the metrics may be taken to be real in regions spacelike separated from a so-called `splitting surface'. This feature is an important hallmark of unitarity in a field theory dual.


title: Customizing Host IDE for Non-programming Users of Pure Embedded DSLs: A Case Study

abstract: Pure embedding as an implementation strategy of domain-specific languages (DSLs) benefits from low implementation costs. On the other hand, it introduces undesired syntactic noise that impedes involvement of non-programming domain experts. Due to this, pure embedded DSLs are generally not intended for, nor used by, non-programmers. In this work, we try to challenge this state by experimenting with inexpensive customizations of the host IDE (Integrated Development Environment) to reduce the negative impact of syntactic noise. We present several techniques and recommendations based on standard IDE features (e.g., file templates, code folding, etc.) that aim to reduce syntactic noise and generally improve the user experience with pure embedded DSLs. The techniques are presented using a NetBeans IDE case study. The goal of the proposed techniques is to improve the user experience with pure embedded DSLs with a focus on the involvement of non-programming domain experts (or non-programmers in general). The proposed techniques were evaluated using a controlled experiment. The experiment compared a group using Ruby and non-modified RubyMine IDE versus a group using Java and NetBeans IDE customized to use the proposed techniques. Experiment results indicate that even inexpensive host IDE customizations can significantly alleviate issues caused by the syntactic noise: Java with its inflexible syntax performed better than Ruby with its concise syntax.