title: Deuteron and Antideuteron Coalescence in Heavy-Ion Collisions: Energy Dependence of the Formation Geometry

abstract: We investigate the collision energy dependence of deuteron and antideuteron emission in the RHIC-BES low- to mid-energy range $\sqrt{s_{NN}} = 4.6-200$ GeV where the formation rate of antinuclei compared to nuclei is strongly suppressed. In the coalescence picture, this can be understood as bulk emission for nuclei in contrast to surface emission for antinuclei. By comparison with experimental data on the coalescence parameter $B_2$, we are able to extract the respective source geometries. This interpretation is further supported by results from the UrQMD transport model, and establishes the following picture: At low energies, nucleons freeze out over the total fireball volume, while antinucleons are annihilated inside the nucleon-rich fireball and can only freeze out on its surface. Towards higher energies, this annihilation effect becomes irrelevant (due to the decreasing baryochemical potential) and the system's freeze-out is driven by the mesons. Thus, the nucleon and antinucleon freeze-out distributions become similar with increasing energy.


title: Existence and Uniqueness of a Fractional Fokker-Planck Equation

abstract: Stochastic differential equations with Levy motion arise the mathematical models for various phenomenon in geophysical and biochemical sciences. The Fokker Planck equation for such a stochastic differential equations is a nonlocal partial differential equations. We prove the existence and uniqueness of the weak solution for this equation.


title: On the flow of states under $T\overline{T}$

abstract: We study the $T\overline{T}$ deformation of two dimensional quantum field theories from a Hamiltonian point of view, focusing on aspects of the theory in Lorentzian signature. Our starting point is a simple rewriting of the spatial integral of the $T\overline{T}$ operator, which directly implies the deformed energy spectrum of the theory. Using this rewriting, we then derive flow equations for various quantities in the deformed theory, such as energy eigenstates, operators, and correlation functions. On the plane, we find that the deformation merely has the effect of implementing successive canonical/Bogoliubov transformations along the flow. This leads us to define a class of non-local, 'dressed' operators (including a dressed stress tensor) which satisfy the same commutation relations as in the undeformed theory. This further implies that on the plane, the deformed theory retains its symmetry algebra, including conformal symmetry, if the original theory is a CFT. On the cylinder the $T\overline{T}$ deformation is much more non-trivial, but even so, correlation functions of certain dressed operators are integral transforms of the original ones. Finally, we propose a tensor network interpretation of our results in the context of AdS/CFT.


title: Collinear expansion for color singlet cross sections

abstract: We demonstrate how to efficiently expand cross sections for color-singlet production at hadron colliders around the kinematic limit of all final state radiation being collinear to one of the incoming hadrons. This expansion is systematically improvable and applicable to a large class of physical observables. We demonstrate the viability of this technique by obtaining the first two terms in the collinear expansion of the rapidity distribution of the gluon fusion Higgs boson production cross section at next-to-next-to leading order (NNLO) in QCD perturbation theory. Furthermore, we illustrate how this technique is used to extract universal building blocks of scattering cross section like the N-jettiness and transverse momentum beam function at NNLO.


title: N-jettiness beam functions at N3LO

abstract: We present the first complete calculation for the quark and gluon $N$-jettiness ($\Tau_N$) beam functions at next-to-next-to-next-to-leading order (N$^3$LO) in perturbative QCD. Our calculation is based on an expansion of the differential Higgs boson and Drell-Yan production cross sections about their collinear limit. This method allows us to employ cutting edge techniques for the computation of cross sections to extract the universal building blocks in question. The class of functions appearing in the matching coefficents for all channels includes iterated integrals with non-rational kernels, thus going beyond the one of harmonic polylogarithms. Our results are a key step in extending the $\Tau_N$ subtraction methods to N$^3$LO, and to resum $\Tau_N$ distributions at N$^3$LL$^\prime$ accuracy both for quark as well as for gluon initiated processes.


title: Distribution and kinematics of 26Al in the Galactic disc

abstract: 26Al is a short-lived radioactive isotope thought to be injected into the interstellar medium (ISM) by massive stellar winds and supernovae. However, all-sky maps of 26Al emission show a distribution with a much larger scale height and faster rotation speed than either massive stars or the cold ISM. We investigate the origin of this discrepancy using an N-body+hydrodynamics simulation of a Milky-Way-like galaxy, self-consistently including self-gravity, star formation, stellar feedback, and 26Al production. We find no evidence that the Milky Way's spiral structure explains the 26Al anomaly. Stars and the 26Al bubbles they produce form along spiral arms, but, because our simulation produces material arms that arise spontaneously rather than propagating arms forced by an external potential, star formation occurs at arm centres rather than leading edges. As a result, we find a scale height and rotation speed for 26Al similar to that of the cold ISM. However, we also show that a synthetic 26Al emission map produced for a possible Solar position at the edge of a large 26Al bubble recovers many of the major qualitative features of the observed 26Al sky. This suggests that the observed anomalous 26Al distribution is the product of foreground emission from the 26Al produced by a nearby, recent supernova.


title: Double Cover of Modular $S_4$ for Flavour Model Building

abstract: We develop the formalism of the finite modular group $\Gamma'_4 \equiv S'_4$, a double cover of the modular permutation group $\Gamma_4 \simeq S_4$, for theories of flavour. The integer weight $k>0$ of the level 4 modular forms indispensable for the formalism can be even or odd. We explicitly construct the lowest-weight ($k=1$) modular forms in terms of two Jacobi theta constants, denoted as $\varepsilon(\tau)$ and $\theta(\tau)$, $\tau$ being the modulus. We show that these forms furnish a 3D representation of $S'_4$ not present for $S_4$. Having derived the $S'_4$ multiplication rules and Clebsch-Gordan coefficients, we construct multiplets of modular forms of weights up to $k=10$. These are expressed as polynomials in $\varepsilon$ and $\theta$, bypassing the need to search for non-linear constraints. We further show that within $S'_4$ there are two options to define the (generalised) CP transformation and we discuss the possible residual symmetries in theories based on modular and CP invariance. Finally, we provide two examples of application of our results, constructing phenomenologically viable lepton flavour models.


title: Eclectic flavor scheme from ten-dimensional string theory -- I. Basic results

abstract: In a consistent top-down approach based on orbifold compactifications, modular and traditional flavor symmetries combine nontrivially to the so-called eclectic flavor symmetry. We extend this scheme from two extra dimensions, discussed previously, to the six extra dimensions of string theory. By doing so, new insights on the nature of CP and its spontaneous breaking emerge. Moreover, we identify a new interpretation of R-symmetries as unbroken remnants from modular symmetries that are associated with geometrically stabilized complex structure moduli. Hence, all symmetries (i.e. modular, traditional flavor, CP and R) share a common origin in string theory: on a technical level, they are given by outer automorphisms of the Narain space group. The eclectic top-down approach leads to a very restrictive scheme with high predictive power. It remains a challenge to connect this with existing bottom-up constructions of modular flavor symmetry.


title: Optimal filters for the moving lens effect

abstract: We assess the prospects for detecting the moving lens effect using cosmological surveys. The bulk motion of cosmological structure induces a small-scale dipolar temperature anisotropy of the cosmic microwave radiation (CMB), centered around halos and oriented along the transverse velocity field. We introduce a set of optimal filters for this signal, and forecast that a high significance detection can be made with upcoming experiments. We discuss the prospects for reconstructing the bulk transverse velocity field on large scales using matched filters, finding good agreement with previous work using quadratic estimators.


title: Optical control of valley Zeeman effect through many-exciton interactions

abstract: Charge carriers in two-dimensional transition metal dichalcogenides (TMDs), such as WSe$_2$, have their spin and valley-pseudospin locked into an optically-addressable index that is proposed as a basis for future information processing. The manipulation of this spin-valley index requires tuning its energy, typically through external magnetic field (B), which is cumbersome. Thus, other efficient routes like all-optical control of spin-valley index are desirable. Here, we show that many-body interactions amongst interlayer excitons in WSe$_2$/MoSe$_2$ heterobilayer induce a steady-state valley Zeeman splitting corresponding to B $\sim$ 6 Tesla. This anomalous splitting, present at incident powers as low as $\mu$Ws, increases with power and enhances, suppresses or even flips the sign of a B-induced splitting. Moreover, the $\it{g}$-factor of valley Zeeman splitting can be tuned by $\sim$ 30 $\%$ with incident power. In addition to valleytronics, our results are relevant for achieving optical non-reciprocity using two-dimensional materials.


title: The GALAH survey: Characterization of emission-line stars with spectral modelling using autoencoders

abstract: We present a neural network autoencoder structure that is able to extract essential latent spectral features from observed spectra and then reconstruct a spectrum from those features. Because of the training with a set of unpeculiar spectra, the network is able to reproduce a spectrum of high signal-to-noise ratio that does not show any spectral peculiarities even if they are present in an observed spectrum. Spectra generated in this manner were used to identify various emission features among spectra acquired by multiple surveys using the HERMES spectrograph at the Anglo-Australian telescope. Emission features were identified by a direct comparison of the observed and generated spectra. Using the described comparison procedure, we discovered 10,364 candidate spectra with a varying degree of H$\alpha$/H$\beta$ emission component produced by different physical mechanisms. A fraction of those spectra belongs to the repeated observation that shows temporal variability in their emission profile. Among emission spectra, we find objects that feature contributions of a nearby rarefied gas (identified through the emission of [NII] and [SII] lines) that was identified in 4004 spectra, which were not all identified as having H$\alpha$ emission. Positions of identified emission-line objects coincide with multiple known regions that harbour young stars. Similarly, detected nebular emission spectra coincide with visually-prominent nebular clouds observable in the red all-sky photographic composites.

categories: astro-ph.SR astro-ph.GA astro-ph.IM physics.data-an


title: Bipolar molecular outflow of the very low-mass star Par-Lup3-4

abstract: Very low-mass stars are known to have jets and outflows, which is indicative of a scaled-down version of low-mass star formation. However, only very few outflows in very low-mass sources are well characterized. We characterize the bipolar molecular outflow of the very low-mass star Par-Lup3-4, a 0.12 M$_{\odot}$ object known to power an optical jet. We observed Par-Lup3-4 with ALMA in Bands 6 and 7, detecting both the continuum and CO molecular gas. In particular, we studied three main emission lines: CO(2-1), CO(3-2), and $^{13}$CO(3-2). Our observations reveal for the first time the base of a bipolar molecular outflow in a very low-mass star, as well as a stream of material moving perpendicular to the primary outflow of this source. The primary outflow morphology is consistent with the previously determined jet orientation and disk inclination. The outflow mass is $9.5\times10^{-7}\mathrm{M}_{\odot}$ , with an outflow rate of $4.3\times10^{-9}\mathrm{M}_{\odot}\mathrm{yr}^{-1}$ A new fitting to the spectral energy distribution suggests that Par-Lup3-4 may be a binary system. We have characterized Par-Lup3-4 in detail, and its properties are consistent with those reported in other very low-mass sources. This source provides further evidence that very low-mass sources form as a scaled-down version of low-mass stars.


title: Inferring the lensing rate of LIGO-Virgo sources from the stochastic gravitational wave background

abstract: Strong lensing of gravitational waves is more likely for distant sources but predicted event rates are highly uncertain with many astrophysical origins proposed. Here we open a new avenue to estimate the event rate of strongly lensed systems by exploring the amplitude of the stochastic gravitational wave background (SGWB). This method can provide a direct upper bound on the high redshift binary coalescing rates, which can be translated into an upper bound on the expected rate of strongly lensed systems. We show that from the ongoing analysis of the Laser Interferometer Gravitational-wave Observatory (LIGO)-Virgo and in the future from the LIGO-Virgo design sensitivity stringent bounds on the lensing event rate can be imposed using the SGWB signal. Combining measurements of loud gravitational wave events with an unresolved stochastic background detection will improve estimates of the numbers of lensed events at high redshift. The proposed method is going to play a crucial in understanding the population of lensed and unlensed systems from gravitational wave observations.


title: Massive black hole merger rates: the effect of kpc separation wandering and supernova feedback

abstract: We revisit the predictions for the merger rate of massive black hole binaries detectable by the Laser Interferometer Space Antenna (LISA) and their background signal for pulsar-timing arrays. We focus on the effect of the delays between the merger of galaxies and the final coalescence of black hole binaries, and on supernova feedback on the black hole growth. By utilizing a semi-analytic galaxy formation model, not only do we account for the driving the evolution of binaries at separations $\lesssim 1$ pc (gas-driven migration, stellar hardening and triple/quadruple massive black hole systems), but we also improve on previous studies by accounting for the time spent by black hole pairs from kpc down to pc separation. We also include the effect of supernova feedback, which may eject gas from the nuclear region of low-mass galaxies, thus hampering the growth of black holes via accretion and suppressing their orbital migration in circumbinary disks. Despite including these novel physical effects, we predict that the LISA detection rate should still be $\gtrsim 2 \mbox{yr}^{-1}$, irrespective of the model for the black hole seeds at high redshifts. Scenarios where black holes form from $\sim100 M_\odot$ seeds are more significantly impacted by supernova feedback. We also find that for detectable events, the merging black holes typically have mass ratios between $\sim 0.1$ and $1$. Predictions for the stochastic background in the band of pulsar-timing array experiments are instead rather robust, and show only a mild dependence on the model.


title: Modeling Photoionized Turbulent Material in the Circumgalactic Medium II: Effect of Turbulence within a Stratified Medium

abstract: The circumgalactic medium (CGM) of nearby star-forming galaxies shows clear indications of OVI absorption accompanied by little to no detectable NV absorption. This unusual spectral signature, accompanied by highly non-uniform absorption from lower ionization state species, indicates that the CGM must be viewed as a dynamic, multiphase medium, such as occurs in the presence of turbulence. Motivated by previous isotropic turbulent simulations, we carry out chemodynamical simulations of stratified media in a Navarro-Frenk-White (NFW) gravitational potential with a total mass of $10^{12}$ solar masses and turbulence that decreases radially. The simulations assume a metallicity of 0.3 solar, a redshift zero metagalatic UV background, and they track ionizations, recombinations, and species-by-species radiative cooling using the MAIHEM package. We compare a suite of ionic column densities with the COS-Halos sample of low-redshift star-forming galaxies. Turbulence with an average one-dimensional velocity dispersion approximately 40 km/s, corresponding to an energy injection rate of approximately $10^{49}$ erg/yr, produces a CGM that matches many of the observed ionic column densities and ratios. In this simulation, the NVI to OVI ratio is suppressed from its equilibrium value due to a combination of radiative cooling and cooling from turbulent mixing. This level of turbulence is consistent with expectations from observations of better constrained, higher-mass systems, and could be sustained by energy input from supernovae, gas inflows, and dynamical friction from dark matter subhalos. We also conduct a higher resolution run which yields smaller-scale structures, but remains in agreement with observations.


title: Exoplanet detection and its dependence on stochastic sampling of the stellar Initial Mass Function

abstract: Young Moving Groups (YMGs) are close (<100pc), coherent collections of young (<100Myr) stars that appear to have formed in the same star-forming molecular cloud. As such we would expect their individual initial mass functions (IMFs) to be similar to other star-forming regions, and by extension the Galactic field. Their close proximity to the Sun and their young ages means that YMGs are promising locations to search for young forming exoplanets. However, due to their low numbers of stars, stochastic sampling of the IMF means their stellar populations could vary significantly. We determine the range of planet-hosting stars (spectral types A, G and M) possible from sampling the IMF multiple times, and find that some YMGs appear deficient in M-dwarfs. We then use these data to show that the expected probability of detecting terrestrial magma ocean planets is highly dependent on the exact numbers of stars produced through stochastic sampling of the IMF.


title: Fault-tolerant quantum error correction using error weight parities

abstract: In quantum error correction using imperfect primitives, errors of high weight arising from a few faults are major concerns since they might not be correctable by the quantum error correcting code. Fortunately, some errors of different weights are logically equivalent and the same correction procedure is applicable to all equivalent errors, thus correcting high-weight errors is sometimes possible. In this work, we introduce a technique called weight parity error correction (WPEC) which can correct Pauli error of any weight in some stabilizer codes provided that the parity of the weight of the error is known. We show that the technique is applicable to concatenated codes constructed from the [[7,1,3]] Steane code or the [[23,1,7]] Golay code. We also provide a fault-tolerant error correction protocol using WPEC for the [[49,1,9]] concatenated Steane code which can correct up to 3 faults and requires only 2 ancillas.


title: Semi-device-dependent blind quantum tomography

abstract: Extracting tomographic information about quantum states is a crucial task in the quest towards devising high-precision quantum devices. Current schemes typically require measurement devices for tomography that are a priori calibrated to high precision. Ironically, the accuracy of the measurement calibration is fundamentally limited by the accuracy of state preparation, establishing a vicious cycle. Here, we prove that this cycle can be broken and the dependence on the measurement device's calibration significantly relaxed. We show that exploiting the natural low-rank structure of quantum states of interest suffices to arrive at a highly scalable `blind' tomography scheme with a classically efficient post-processing algorithm. We further improve the efficiency of our scheme by making use of the sparse structure of the calibrations. This is achieved by relaxing the blind quantum tomography problem to the de-mixing of a sparse sum of low-rank matrices. We prove that the proposed algorithm recovers a low-rank quantum state and the calibration provided that the measurement model exhibits a restricted isometry property. For generic measurements, we show that it requires a close-to-optimal number of measurement settings. Complementing these conceptual and mathematical insights, we numerically demonstrate that robust blind quantum tomography is possible in a practical setting inspired by an implementation of trapped ions.


title: Quantum computer-aided design: digital quantum simulation of quantum processors

abstract: With the increasing size of quantum processors, sub-modules that constitute the processor hardware will become too large to accurately simulate on a classical computer. Therefore, one would soon have to fabricate and test each new design primitive and parameter choice in time-consuming coordination between design, fabrication, and experimental validation. Here we show how one can design and test the performance of next-generation quantum hardware -- by using existing quantum computers. Focusing on superconducting transmon processors as a prominent hardware platform, we compute the static and dynamic properties of individual and coupled transmons. We show how the energy spectra of transmons can be obtained by variational hybrid quantum-classical algorithms that are well-suited for near-term noisy quantum computers. In addition, single- and two-qubit gate simulations are demonstrated via Suzuki-Trotter decomposition. Our methods pave a promising way towards designing candidate quantum processors when the demands of calculating sub-module properties exceed the capabilities of classical computing resources.


title: Entangling logical qubits with lattice surgery

abstract: Future quantum computers will require quantum error correction for faithful operation. The correction capabilities come with an overhead for performing fault-tolerant logical operations on the encoded qubits. One of the most resource efficient ways to implement logical operations is lattice surgery, where groups of physical qubits, arranged on lattices, can be merged and split to realize entangling gates and teleport logical information. Here, we report on the experimental realization of lattice surgery between two topologically encoded qubits in a 10-qubit ion trap quantum information processor. In particular, we demonstrate entanglement between two logical qubits and we implement logical state teleportation.


title: Circumnuclear Molecular Gas in Low-redshift Quasars and Matched Star-forming Galaxies

abstract: A series of gravitational instabilities in a circumnuclear gas disk (CND) are required to trigger gas transport to a central supermassive black hole (SMBH) and ignite Active Galactic Nuclei (AGNs). A test of this scenario is to investigate whether an enhanced molecular gas mass surface density ($\Sigma_{\rm mol}$) is found in the CND-scale of quasars relative to a comparison sample of inactive galaxies. Here we performed sub-kpc resolution CO(2-1) observations with ALMA of four low-redshift ($z \sim 0.06$), luminous ($\sim 10^{45}$ erg s$^{-1}$) quasars with each matched to a different star-forming galaxy, having similar redshift, stellar mass, and star-formation rate. We detected CO(2-1) emission from all quasars, which show diverse morphologies. Contrary to expectations, $\Sigma_{\rm mol}$ of the quasar sample, computed from the CO(2-1) luminosity, tends to be smaller than the comparison sample at $r < 500$ pc; there is no systematic enhancement of $\Sigma_{\rm mol}$ in our quasars. We discuss four possible scenarios that would explain the lower molecular gas content (or CO(2-1) luminosity as an actual observable) at the CND-scale of quasars, i.e., AGN-driven outflows, gas-rich minor mergers, time-delay between the onsets of a starburst-phase and a quasar-phase, and X-ray-dominated region (XDR) effects on the gas chemical abundance and excitation. While not extensively discussed in the literature, XDR effects can have an impact on molecular mass measurements particularly in the vicinity of luminous quasar nuclei; therefore higher resolution molecular gas observations, which are now viable using ALMA, need to be considered.


title: The lifetimes of planetary debris discs around white dwarfs

abstract: The lifetime of a planetary disc which orbits a white dwarf represents a crucial input parameter into evolutionary models of that system. Here we apply a purely analytical formalism to estimate lifetimes of the debris phase of these discs, before they are ground down into dust or are subject to sublimation from the white dwarf. We compute maximum lifetimes for three different types of white dwarf discs, formed from (i) radiative YORP breakup of exo-asteroids along the giant branch phases at 2-100 au, (ii) radiation-less spin-up disruption of these minor planets at $\sim 1.5-4.5R_{\odot}$, and (iii) tidal disruption of minor or major planets within about $1.3R_{\odot}$. We display these maximum lifetimes as a function of disc mass and extent, constituent planetesimal properties, and representative orbital excitations of eccentricity and inclination. We find that YORP discs with masses up to $10^{24}$ kg live long enough to provide a reservoir of surviving cm-sized pebbles and m- to km-sized boulders that can be perturbed intact to white dwarfs with cooling ages of up to 10 Gyr. Debris discs formed from the spin or tidal disruption of these minor planets or major planets can survive in a steady state for up to respectively 1 Myr or 0.01 Myr, although most tidal discs would leave a steady state within about 1 yr. Our results illustrate that dust-less planetesimal transit detections are plausible, and would provide particularly robust evolutionary constraints. Our formalism can easily be adapted to individual systems and future discoveries.


title: A NICER look at the state transitions of the black hole candidate MAXI J1535-571 during its reflares

abstract: The black hole candidate and X-ray binary MAXI J1535-571 was discovered in September 2017. During the decay of its discovery outburst, and before returning to quiescence, the source underwent at least four reflaring events, with peak luminosities of $\sim$10$^{35-36}$ erg s$^{-1}$ (d/4.1 kpc)$^2$. To investigate the nature of these flares, we analysed a sample of NICER observations taken with almost daily cadence. In this work we present the detailed spectral and timing analysis of the evolution of the four reflares. The higher sensitivity of NICER at lower energies, in comparison with other X-ray detectors, allowed us to constrain the disc component of the spectrum at $\sim$0.5 keV. We found that during each reflare the source appears to trace out a q-shaped track in the hardness-intensity diagram similar to those observed in black hole binaries during full outbursts. MAXI J1535-571 transits between the hard state (valleys) and softer states (peaks) during these flares. Moreover, the Comptonised component is undetected at the peak of the first reflare, while the disc component is undetected during the valleys. Assuming the most likely distance of 4.1 kpc, we find that the hard-to-soft transitions take place at the lowest luminosities ever observed in a black hole transient, while the soft-to-hard transitions occur at some of the lowest luminosities ever reported for such systems.


title: Quantum Computer-Aided design of Quantum Optics Hardware

abstract: The parameters of a quantum system grow exponentially with the number of involved quantum particles. Hence, the associated memory requirement goes well beyond the limit of best classic computers for quantum systems composed of a few dozen particles leading to huge challenges in their numerical simulation. This implied that verification, let alone, design of new quantum devices and experiments, is fundamentally limited to small system size. It is not clear how the full potential of large quantum systems can be exploited. Here, we present the concept of quantum computer designed quantum hardware and apply it to the field of quantum optics. Specifically, we map complex experimental hardware for high-dimensional, many-body entangled photons into a gate-based quantum circuit. We show explicitly how digital quantum simulation of Boson Sampling experiments can be realized. Then we illustrate how to design quantum-optical setups for complex entangled photon systems, such as high-dimensional Greenberger-Horne-Zeilinger states and their derivatives. Since photonic hardware is already on the edge of quantum supremacy (the limit beyond which systems can no longer be calculated classically) and the development of gate-based quantum computers is rapidly advancing, our approach promises to be an useful tool for the future of quantum device design.


title: Geometric Flow Equations for Schwarzschild-AdS Space-time and Hawking-Page Phase Transition

abstract: Following the recent observation that the Ricci flow and the infinite distance swampland conjecture are closely related to each other, we will investigate in this paper geometric flow equations for AdS space-time geometries. First, we consider the so called Yamabe and Ricci-Bourguignon flows and we show that these two flows - in contrast to the Ricci flow - lead to infinite distance fixed points for product spaces like $AdS_d\times S^p$, where $AdS_d$ denotes d-dimensional AdS space and $S^p$ corresponds to a p-dimensional sphere. Second, we consider black hole geometries in AdS space time geometries and their behaviour under the Yamabe and Ricci-Bourguignon flows. Specifically we will examine if and how the AdS black holes will undergo a Hawking-Page phase transition under the Ricci flow, the Yamabe flow and under the general Ricci-Bourguignon flow.


title: Planet Migration in Self-Gravitating Discs: Survival of Planets

abstract: We carry out three-dimensional SPH simulations to study whether planets can survive in self-gravitating protoplanetary discs. The discs modelled here use a cooling prescription that mimics a real disc which is only gravitationally unstable in the outer regions. We do this by modelling the cooling using a simplified method such that the cooling time in the outer parts of the disc is shorter than in the inner regions, as expected in real discs. We find that both giant (> M_Sat) and low mass (< M_Nep) planets initially migrate inwards very rapidly, but are able to slow down in the inner gravitationally stable regions of the disc without needing to open up a gap. This is in contrast to previous studies where the cooling was modelled in a more simplified manner where regardless of mass, the planets were unable to slow down their inward migration. This shows the important effect the thermodynamics has on planet migration. In a broader context, these results show that planets that form in the early stages of the discs' evolution, when they are still quite massive and self-gravitating, can survive.


title: The harmonic polytope

abstract: We study the harmonic polytope, which arose in Ardila, Denham, and Huh's work on the Lagrangian geometry of matroids. We describe its combinatorial structure, showing that it is a $(2n-2)$-dimensional polytope with $(n!)^2(1+\frac12+\cdots+\frac1n)$ vertices and $3^n-3$ facets. We also give a formula for its volume: it is a weighted sum of the degrees of the projective varieties of all the toric ideals of connected bipartite graphs with $n$ edges; or equivalently, a weighted sum of the lattice point counts of all the corresponding trimmed generalized permutahedra.


title: The Importance of Open-Endedness (for the Sake of Open-Endedness)

abstract: A paper in the recent Artificial Life journal special issue on open-ended evolution (OEE) presents a simple evolving computational system that, it is claimed, satisfies all proposed requirements for OEE (Hintze, 2019). Analysis and discussion of the system are used to support the further claims that complexity and diversity are the crucial features of open-endedness, and that we should concentrate on providing proper definitions for those terms rather than engaging in "the quest for open-endedness for the sake of open-endedness" (Hintze, 2019, p. 205). While I wholeheartedly support the pursuit of precise definitions of complexity and diversity in relation to OEE research, I emphatically reject the suggestion that OEE is not a worthy research topic in its own right. In the same issue of the journal, I presented a "high-level conceptual framework to help orient the discussion and implementation of open-endedness in evolutionary systems" (Taylor, 2019). In the current brief contribution I apply my framework to Hinzte's model to understand its limitations. In so doing, I demonstrate the importance of studying open-endedness for the sake of open-endedness.


title: A Census of Star Formation in the Outer Galaxy II: The GLIMPSE360 Field

abstract: We have conducted a study of star formation in the outer Galaxy from 65\degr$< l <$265\degr~in the region observed by the GLIMPSE360 program. This {\it Spitzer} warm mission program mapped the plane of the outer Milky Way with IRAC at 3.6 and 4.5~$\mu$m. We combine the IRAC, {\it WISE}, and 2MASS catalogs and our previous results from another outer Galaxy survey and identify a total of 47,338 Young Stellar Objects (YSOs) across the field spanning $>$180\degr\ in Galactic longitude. Using the $DBSCAN$ method on the combined catalog, we identify 618 clusters or aggregations of YSOs having 5 or more members. We identify 10,476 Class I, 29,604 Class II, and 7,325 anemic Class II/Class III YSOs. The ratio of YSOs identified as members of clusters was 25,528/47,338, or 54\%. We found 100 of the clusters identified have previously measured distances in the {\it WISE} \ion{H}{2} survey. We used these distances in our spectral energy distribution (SED) fitting of the YSOs in these clusters, of which 96 had YSOs with $<3\sigma$ fits. We used the derived masses from the SED model fits to estimate the initial mass function (IMF) in the inner and outer Galaxy clusters: dividing the clusters by Galactocentric distances, the slopes were $\Gamma = 1.87 \pm 0.31$ above 3~M$_{\odot}$ for $R_{Gal} < 11.5$~kpc and $\Gamma = 1.15 \pm 0.24$ above 3~M$_{\odot}$ for $R_{Gal} > 11.5$~kpc. The slope of the combined IMF was found to be $\Gamma = 1.92 \pm 0.42$ above 3~M$_{\odot}$. These values are consistent with each other within the uncertainties, and with literature values in the inner Galaxy high-mass star formation regions. The slopes are likely also consistent with a universal Salpeter IMF.


title: The synchrotron maser emission from relativistic magnetized shocks: Dependence on the pre-shock temperature

abstract: Electromagnetic precursor waves generated by the synchrotron maser instability at relativistic magnetized shocks have been recently invoked to explain the coherent radio emission of Fast Radio Bursts. By means of two-dimensional particle-in-cell simulations, we explore the properties of the precursor waves in relativistic electron-positron perpendicular shocks as a function of the pre-shock magnetization $\sigma\gtrsim 1$ (i.e., the ratio of incoming Poynting flux to particle energy flux) and thermal spread $\Delta \gamma \equiv kT/mc^2=10^{-5}-10^{-1}$. We measure the fraction $f_\xi$ of total incoming energy that is converted into precursor waves, as computed in the post-shock frame. At fixed magnetization, we find that $f_{\xi}$ is nearly independent of temperature as long as $\Delta \gamma \lesssim 10^{-1.5}$ (with only a modest decrease of a factor of three from $\Delta\gamma =10^{-5}$ to $\Delta\gamma =10^{-1.5}$), but it drops by nearly two orders of magnitude for $\Delta \gamma \gtrsim 10^{-1}$. At fixed temperature, the scaling with magnetization $f_\xi\sim 10^{-3}\,\sigma^{-1}$ is consistent with our earlier one-dimensional results. For our reference $\sigma=1$, the power spectrum of precursor waves is relatively broad (fractional width $\sim 1-3$) for cold temperatures, whereas it shows pronounced line-like features with fractional width $\sim 0.2$ for $10^{-3} \lesssim \Delta \gamma \lesssim 10^{-1.5} $. For $\sigma\gtrsim 1$, the precursor waves are beamed within an angle $\simeq \sigma^{-1/2}$ from the shock normal (as measured in the post-shock frame), as required so they can outrun the shock. Our results can provide physically-grounded inputs for FRB emission models based on maser emission from relativistic shocks.


title: Chern-Simons Invariants and Heterotic Superpotentials

abstract: The superpotential in four-dimensional heterotic effective theories contains terms arising from holomorphic Chern-Simons invariants associated to the gauge and tangent bundles of the compactification geometry. These effects are crucial for a number of key features of the theory, including vacuum stability and moduli stabilization. Despite their importance, few tools exist in the literature to compute such effects in a given heterotic vacuum. In this work we present new techniques to explicitly determine holomorphic Chern-Simons invariants in heterotic string compactifications. The key technical ingredient in our computations are real bundle morphisms between the gauge and tangent bundles. We find that there are large classes of examples, beyond the standard embedding, where the Chern-Simons superpotential vanishes. We also provide explicit examples for non-flat bundles where it is non-vanishing and fractionally quantized, generalizing previous results for Wilson lines.


title: Asymptotic behavior of a network of neurons with random linear interactions

abstract: We study the asymptotic behavior for asymmetric neuronal dynamics in a network of linear Hopfield neurons. The interaction between the neurons is modeled by random couplings which are centered i.i.d. random variables with finite moments of all orders. We prove that if the initial condition of the network is a set of i.i.d. random variables and independent of the synaptic weights, each component of the limit system is described as the sum of the corresponding coordinate of the initial condition with a centered Gaussian process whose covariance function can be described in terms of a modified Bessel function. This process is not Markovian. The convergence is in law almost surely with respect to the random weights. Our method is essentially based on the method of moments to obtain a Central Limit Theorem.


title: Tunable quantum interference using a topological source of indistinguishable photon pairs

abstract: Sources of quantum light, in particular correlated photon pairs that are indistinguishable in all degrees of freedom, are the fundamental resource that enables continuous-variable quantum computation and paradigms such as Gaussian boson sampling. Nanophotonic systems offer a scalable platform for implementing sources of indistinguishable correlated photon pairs. However, such sources have so far relied on the use of a single component, such as a single waveguide or a ring resonator, which offers limited ability to tune the spectral and temporal correlations between photons. Here, we demonstrate the use of a topological photonic system comprising a two-dimensional array of ring resonators to generate indistinguishable photon pairs with dynamically tunable spectral and temporal correlations. Specifically, we realize dual-pump spontaneous four-wave mixing in this array of silicon ring resonators that exhibits topological edge states. We show that the linear dispersion of the edge states over a broad bandwidth allows us to tune the correlations, and therefore, quantum interference between photons by simply tuning the two pump frequencies in the edge band. Furthermore, we demonstrate energy-time entanglement between generated photons. We also show that our topological source is inherently protected against fabrication disorders. Our results pave the way for scalable and tunable sources of squeezed light that are indispensable for quantum information processing using continuous variables.


title: Hierarchical hyperbolicity of graph products

abstract: We show that any graph product of finitely generated groups is hierarchically hyperbolic relative to its vertex groups. We apply this result to answer two questions of Behrstock, Hagen, and Sisto: we show that the syllable metric on any graph product forms a hierarchically hyperbolic space, and that graph products of hierarchically hyperbolic groups are themselves hierarchically hyperbolic groups. This last result is a strengthening of a result of Berlai and Robbio by removing the need for extra hypotheses on the vertex groups. We also answer two questions of Genevois about the geometry of the electrification of a graph product of finite groups.


title: Augmented fidelities for single qubit gates

abstract: An average gate fidelity is a standard performance metric to quantify deviation between an ideal unitary gate transformation and its realistic experimental implementation. The average is taken with respect to states uniformly distributed over the full Hilbert space. We analytically (single-qubit) and numerically (two-qubit) show how this average changes if the uniform distribution condition is relaxed, replaced by parametrized distributions - polar cap and von Mises-Fisher distributions - and how the resulting fidelities can differentiate certain noise models. In particular, we demonstrate that Pauli channels with different noise rates along the three axes can be faithfully distinguished using these augmented fidelities.


title: Fermionic systems for quantum information people

abstract: The operator algebra of fermionic modes is isomorphic to that of qubits, the difference between them is twofold: the embedding of subalgebras corresponding to mode subsets and multiqubit subsystems on the one hand, and the parity superselection in the fermionic case on the other. We discuss these two fundamental differences extensively, and illustrate these through the Jordan--Wigner representation in a coherent, self-contained, pedagogical way, from the point of view of quantum information theory. Our perspective leads us to develop useful new tools for the treatment of fermionic systems, such as the fermionic (quasi-)tensor product, fermionic canonical embedding, fermionic partial trace, fermionic products of maps and fermionic embeddings of maps. We formulate these by direct, easily applicable formulas, without mode permutations, for arbitrary partitionings of the modes. It is also shown that fermionic reduced states can be calculated by the fermionic partial trace, containing the proper phase factors. We also consider variants of the notions of fermionic mode correlation and entanglement, which can be endowed with the usual, local operation based motivation, if the parity superselection rule is imposed. We also elucidate some other fundamental points, related to joint map extensions, which make the parity superselection inevitable in the description of fermionic systems.


title: A Real-Time Closed-Form Model for Nonlinearity Modeling in Ultra-Wide-Band Optical Fiber Links Accounting for Inter-channel Stimulated Raman Scattering and Co-Propagating Raman Amplification

abstract: In this paper, we present a novel closed-form model (CFM) for accurate and fast evaluation of nonlinear interference in modern ultrawideband coherent optical fiber communication systems. Starting from the Gaussian noise model (GN model), using reasonable approximations and machine-learning optimized improvements, we achieve an accurate CFM capable of handling ultrawide band (C+L or wider) optical fiber systems in the presence of Inter-channel Stimulated Raman Scattering (ISRS) and forward-pumped Raman amplification.


title: Towards Understanding Fast Adversarial Training

abstract: Current neural-network-based classifiers are susceptible to adversarial examples. The most empirically successful approach to defending against such adversarial examples is adversarial training, which incorporates a strong self-attack during training to enhance its robustness. This approach, however, is computationally expensive and hence is hard to scale up. A recent work, called fast adversarial training, has shown that it is possible to markedly reduce computation time without sacrificing significant performance. This approach incorporates simple self-attacks, yet it can only run for a limited number of training epochs, resulting in sub-optimal performance. In this paper, we conduct experiments to understand the behavior of fast adversarial training and show the key to its success is the ability to recover from overfitting to weak attacks. We then extend our findings to improve fast adversarial training, demonstrating superior robust accuracy to strong adversarial training, with much-reduced training time.


title: Motion by mean curvature in interacting particle systems

abstract: There are a number of situations in which rescaled interacting particle systems have been shown to converge to a reaction diffusion equation (RDE) with a bistable reaction term. These RDEs have traveling wave solutions. When the speed of the wave is nonzero, block constructions have been used to prove the existence or nonexistence of nontrivial stationary distributions. Here, we follow the approach in a paper by Etheridge, Freeman, and Pennington to show that in a wide variety of examples when the RDE limit has a bistable reaction term and traveling waves have speed 0, one can run time faster and further rescale space to obtain convergence to motion by mean curvature. This opens up the possibility of proving that the sexual reproduction model with fast stirring has a discontinuous phase transition, and that in Region 2 of the phase diagram for the nonlinear voter model studied by Molofsky et al there were two nontrivial stationary distributions.


title: The FRET-based structural dynamics challenge -- community contributions to consistent and open science practices

abstract: Single-molecule F\"{o}rster resonance energy transfer (smFRET) has become a mainstream technique for probing biomolecular structural dynamics. The rapid and wide adoption of the technique by an ever-increasing number of groups has generated many improvements and variations in the technique itself, in methods for sample preparation and characterization, in analysis of the data from such experiments, and in analysis codes and algorithms. Recently, several labs that employ smFRET have joined forces to try to bring the smFRET community together in adopting a consensus on how to perform experiments and analyze results for achieving quantitative structural information. These recent efforts include multi-lab blind-tests to assess the accuracy and precision of smFRET between different labs using different procedures, the formal assembly of the FRET community and development of smFRET procedures to be considered for entries in the wwPDB. Here we delve into the different approaches and viewpoints in the field. This position paper describes the current "state-of-the field", points to unresolved methodological issues for quantitative structural studies, provides a set of 'soft recommendations' about which an emerging consensus exists, and a list of resources that are openly available. To make further progress, we strongly encourage 'open science' practices. We hope that this position paper will provide a roadmap for newcomers to the field, as well as a reference for seasoned practitioners.


title: Generation of inflationary perturbations in the continuous spontaneous localization model: The second order power spectrum

abstract: Cosmic inflation, which describes an accelerated expansion of the early Universe, yields the most successful predictions regarding temperature anisotropies in the cosmic microwave background (CMB). Nevertheless, the precise origin of the primordial perturbations and their quantum-to-classical transition is still an open issue. The continuous spontaneous localization model (CSL), in the cosmological context, might be used to provide a solution to the mentioned puzzles by considering an objective reduction of the inflaton wave function. In this work, we calculate the primordial power spectrum at the next leading order in the Hubble flow functions that results from applying the CSL model to slow roll inflation within the semiclassical gravity framework. We employ the method known as uniform approximation along with a second order expansion in the Hubble flow functions. We analyze some features in the CMB temperature and primordial power spectra that could help to distinguish between the standard prediction and our approach.


title: Surface properties of semiconductors from post-illumination photovoltage transient

abstract: Free surfaces of semiconductors respond to light by varying their surface voltage (surface band bending). This surface photovoltage may be easily detected using a Kelvin probe. Modeling the transient temporal behavior of the surface photovoltage after the light is turned off may serve as a means to characterize several key electronic properties of the semiconductor, which are of fundamental importance in numerous electronic device applications, such as transistors and solar cells. In this paper, we develop a model for this temporal behavior and use it to characterize layers and nanowires of several semiconductors. Our results suggest that what has previously been considered to be a logarithmic decay is only approximately so. Due to the known limited frequency bandwidth of the Kelvin probe method, most previous Kelvin-probe-based methods have been limited to slow responding semiconductors. The model we propose extends this range of applicability.


title: An abstract characterization for projections in operator systems

abstract: We show that the set of projections in an operator system can be detected using only the abstract data of the operator system. Specifically, we show that if $p$ is a positive contraction in an operator system $V$ which satisfies certain order-theoretic conditions, then there exists a complete order embedding of $V$ into $B(H)$ mapping $p$ to a projection operator. Moreover, every abstract projection in an operator system $V$ is an honest projection in the C*-envelope of $V$. Using this characterization, we provide an abstract characterization for operator systems spanned by two commuting families of projection-valued measures and discuss applications in quantum information theory.


title: Spin-induced black hole spontaneous scalarization

abstract: We study scalar fields in a black hole background and show that, when the scalar is suitably coupled to curvature, rapid rotation can induce a tachyonic instability. This instability, which is the hallmark of spontaneous scalarization in the linearized regime, is expected to be quenched by nonlinearities and endow the black hole with scalar hair. Hence, our results demonstrate the existence of a broad class of theories that share the same stationary black hole solutions with general relativity at low spins, but which exhibit black hole hair at sufficiently high spins ($a/M\gtrsim 0.5$). This result has clear implications for tests of general relativity and the nature of black holes with gravitational and electromagnetic observations.


title: SOLO: A Corpus of Tweets for Examining the State of Being Alone

abstract: The state of being alone can have a substantial impact on our lives, though experiences with time alone diverge significantly among individuals. Psychologists distinguish between the concept of solitude, a positive state of voluntary aloneness, and the concept of loneliness, a negative state of dissatisfaction with the quality of one's social interactions. Here, for the first time, we conduct a large-scale computational analysis to explore how the terms associated with the state of being alone are used in online language. We present SOLO (State of Being Alone), a corpus of over 4 million tweets collected with query terms 'solitude', 'lonely', and 'loneliness'. We use SOLO to analyze the language and emotions associated with the state of being alone. We show that the term 'solitude' tends to co-occur with more positive, high-dominance words (e.g., enjoy, bliss) while the terms 'lonely' and 'loneliness' frequently co-occur with negative, low-dominance words (e.g., scared, depressed), which confirms the conceptual distinctions made in psychology. We also show that women are more likely to report on negative feelings of being lonely as compared to men, and there are more teenagers among the tweeters that use the word 'lonely' than among the tweeters that use the word 'solitude'.


title: The Jahn-Teller effect and spin--orbit coupling: friends or foes?

abstract: The Jahn-Teller effect is one of the most fundamental phenomena important not only for physics, but also for chemistry and material science. Solving the Jahn-Teller problem and taking into account strong electron correlations we show that quantum entanglement of the spin and orbital degrees of freedom via spin--orbit coupling strongly affects this effect. Depending on the number of $d$ electrons it may quench (electronic configurations $t_{2g}^2$, $t_{2g}^4$, and $t_{2g}^5$), partially suppress ($t_{2g}^1$) or in contrast induce ($t_{2g}^3$) Jahn-Teller distortions. Moreover, in certain situations, interplay between the Jahn-Teller effect and spin--orbit coupling promotes formation of the "Mexican hat" energy surface facilitating various quantum phenomena.


title: Dependence of Earth's Thermal Radiation on Five Most Abundant Greenhouse Gases

abstract: The atmospheric temperatures and concentrations of Earth's five most important, greenhouse gases, H$_2$O, CO$_2$, O$_3$, N$_2$O and CH$_4$ control the cloud-free, thermal radiative flux from the Earth to outer space. Over 1/3 million lines having strengths as low as $10^{-27}$ cm of the HITRAN database were used to evaluate the dependence of the forcing on the gas concentrations. For a hypothetical, optically thin atmosphere, where there is negligible saturation of the absorption bands, or interference of one type of greenhouse gas with others, the per-molecule forcings are of order $10^{-22}$ W for H$_2$O, CO$_2$, O$_3$, N$_2$O and CH$_4$. For current atmospheric concentrations, the per-molecule forcings of the abundant greenhouse gases H$_2$O and CO$_2$ are suppressed by four orders of magnitude. The forcings of the less abundant greenhouse gases, O$_3$, N$_2$O and CH$_4$, are also suppressed, but much less so. For current concentrations, the per-molecule forcings are two to three orders of magnitude greater for O$_3$, N$_2$O and CH$_4$, than those of H$_2$O or CO$_2$. Doubling the current concentrations of CO$_2$, N$_2$O or CH$_4$ increases the forcings by a few per cent. These forcing results are close to previously published values even though the calculations did not utilize either a CO$_2$ or H$_2$O continuum. The change in surface temperature due to CO$_2$ doubling is estimated taking into account radiative-convective equilibrium of the atmosphere as well as water feedback for the cases of fixed absolute and relative humidities as well as the effect of using a pseudoadiabatic lapse rate to model the troposphere temperature. Satellite spectral measurements at various latitudes are in excellent quantitative agreement with modelled intensities.


title: Enhancing the effect of quantum many-body scars on dynamics by minimising the effective dimension

abstract: Quantum many-body scarring is believed to be the mechanism behind long-lived coherent oscillations in interacting Rydberg atom chains. These persistent oscillations are due to the large overlap of the many-body scars with certain initial states. We show that the "effective dimension" is a useful measure for identifying non-thermalising initial states in many-body scarred systems. By minimising the effective dimension we find physically reasonable initial states of the Rydberg chain that lead to more pronounced and longer lived oscillations, accentuating the effect of the many-body scars on the dynamics.


title: Steady gradient K\"ahler-Ricci solitons on crepant resolutions of Calabi-Yau cones

abstract: We show that, up to the flow of the soliton vector field, there exists a unique complete steady gradient K\"ahler-Ricci soliton in every K\"ahler class of an equivariant crepant resolution of a Calabi-Yau cone converging at a polynomial rate to Cao's steady gradient K\"ahler-Ricci soliton on the cone.


title: Statistical mechanical approximations to more efficiently determine polymorph free energy differences for small organic molecules

abstract: Methods to efficiently determine the relative stability of polymorphs of organic crystals are highly desired in crystal structure predictions (CSPs). Current methodologies include use of static lattice phonons, quasi-harmonic approximation (QHA), and computing the full thermodynamic cycle using replica exchange molecular dynamics (REMD). We found that 13 out of the 29 systems minimized from experiment restructured to a lower energy minima when heated using REMD, a phenomena that QHA cannot capture. Here, we present a series of methods that are intermediate in accuracy and expense between QHA and computing the full thermodynamic cycle which can save 42-80% of the computational cost and introduces, on this benchmark, a relatively small (0.16 +/- 0.04 kcal/mol) error relative to the full pseudosupercritical path approach. In particular, a method that Boltzmann weights the harmonic free energy of the trajectory of an REMD replica appears to be an appropriate intermediate between QHA and full thermodynamic cycle using MD when screening crystal polymorph stability.


title: Interactive Music and Synchronous Reactive Programming

abstract: This paper presents Skini, a programming methodology and an execution environment for interactive structured music. With this system, the composer programs his scores in the HipHop.js synchronous reactive language. They are then executed, or played, in live concerts, in interaction with the audience. The system aims at helping composers to find a good balance between the determinism of the compositions and the nondeterminism of the interactions with the public. Each execution of a Skini score yields to a different but aesthetically consistent interpretation. This work raises many questions in the musical fields. How to combine composition and interaction? How to control the musical style when the audience influences what is to play next? What are the possible connections with generative music? These are important questions for the Skini system but they are out of the scope of this paper that focuses exclusively on the computer science aspects of the system. From that perspective, the main questions are how to program the scores and in which language? General purpose languages are inappropriate because their elementary constructs (i.e., variables, functions, loops, etc.) do not match the constructions needed to express music and musical constraints. We show that synchronous programming languages are a much better fit because they rely on temporal constructs that can be directly used to represent musical scores and because their malleability enables composers to experiment easily with artistic variations of their initial scores. The paper mostly focuses on scores programming. It exposes the process a composer should follow from his very first musical intuitions up to the generation of a musical artifact. The paper presents some excerpts of the programming of a classical music composition that it then precisely relates to an actual recording. Examples of techno music and jazz are also presented, with audio artifact, to demonstrate the versatility of the system. Finally, brief presentations of past live concerts are presented as an evidence of viability of the system.


title: SMIE: Weakness is Power!: Auto-indentation with incomplete information

abstract: Automatic indentation of source code is fundamentally a simple matter of parsing the code and then applying language- and style-specific rules about relative indentation of the various constructs. Yet, in practice, full parsing is not always an option, either because of quirks of the language, or because the code is temporarily syntactically incorrect, or because of an incomplete or broken grammar. I present the design of Emacs's Simple-Minded Indentation Engine (SMIE), which gets its power from the weakness of the underlying parsing technique. It makes it possible to perform local parsing, which is hence unaffected by irrelevant surrounding code. This provides a form of graceful degradation in the face of incomplete, erroneous, or just plain problematic information.


title: Portability of Scientific Workflows in NGS Data Analysis: A Case Study

abstract: The analysis of next-generation sequencing (NGS) data requires complex computational workflows consisting of dozens of autonomously developed yet interdependent processing steps. Whenever large amounts of data need to be processed, these workflows must be executed on a parallel and/or distributed systems to ensure reasonable runtime. Porting a workflow developed for a particular system on a particular hardware infrastructure to another system or to another infrastructure is non-trivial, which poses a major impediment to the scientific necessities of workflow reproducibility and workflow reusability. In this work, we describe our efforts to port a state-of-the-art workflow for the detection of specific variants in whole-exome sequencing of mice. The workflow originally was developed in the scientific workflow system snakemake for execution on a high-performance cluster controlled by Sun Grid Engine. In the project, we ported it to the scientific workflow system SaasFee that can execute workflows on (multi-core) stand-alone servers or on clusters of arbitrary sizes using the Hadoop. The purpose of this port was that also owners of low-cost hardware infrastructures, for which Hadoop was made for, become able to use the workflow. Although both the source and the target system are called scientific workflow systems, they differ in numerous aspects, ranging from the workflow languages to the scheduling mechanisms and the file access interfaces. These differences resulted in various problems, some expected and more unexpected, that had to be resolved before the workflow could be run with equal semantics. As a side-effect, we also report cost/runtime ratios for a state-of-the-art NGS workflow on very different hardware platforms: A comparably cheap stand-alone server (80 threads), a mid-cost, mid-sized cluster (552 threads), and a high-end HPC system (3784 threads).


title: Defining Estimands Using a Mix of Strategies to Handle Intercurrent Events in Clinical Trials

abstract: Randomized controlled trials (RCT) are the gold standard for evaluation of the efficacy and safety of investigational interventions. If every patient in an RCT were to adhere to the randomized treatment, one could simply analyze the complete data to infer the treatment effect. However, intercurrent events (ICEs) including the use of concomitant medication for unsatisfactory efficacy, treatment discontinuation due to adverse events, or lack of efficacy, may lead to interventions that deviate from the original treatment assignment. Therefore, defining the appropriate estimand (the appropriate parameter to be estimated) based on the primary objective of the study is critical prior to determining the statistical analysis method and analyzing the data. The International Council for Harmonisation (ICH) E9 (R1), published on November 20, 2019, provided 5 strategies to define the estimand: treatment policy, hypothetical, composite variable, while on treatment and principal stratum. In this article, we propose an estimand using a mix of strategies in handling ICEs. This estimand is an average of the null treatment difference for those with ICEs potentially related to safety and the treatment difference for the other patients if they would complete the assigned treatments. Two examples from clinical trials evaluating anti-diabetes treatments are provided to illustrate the estimation of this proposed estimand and to compare it with the estimates for estimands using hypothetical and treatment policy strategies in handling ICEs.


title: Model-Based Generalization Under Parameter Uncertainty Using Path Integral Control

abstract: This work addresses the problem of robot interaction in complex environments where online control and adaptation is necessary. By expanding the sample space in the free energy formulation of path integral control, we derive a natural extension to the path integral control that embeds uncertainty into action and provides robustness for model-based robot planning. Our algorithm is applied to a diverse set of tasks using different robots and validate our results in simulation and real-world experiments. We further show that our method is capable of running in real-time without loss of performance. Videos of the experiments as well as additional implementation details can be found at https://sites.google.com/view/emppi.


title: Attention and Encoder-Decoder based models for transforming articulatory movements at different speaking rates

abstract: While speaking at different rates, articulators (like tongue, lips) tend to move differently and the enunciations are also of different durations. In the past, affine transformation and DNN have been used to transform articulatory movements from neutral to fast(N2F) and neutral to slow(N2S) speaking rates [1]. In this work, we improve over the existing transformation techniques by modeling rate specific durations and their transformation using AstNet, an encoder-decoder framework with attention. In the current work, we propose an encoder-decoder architecture using LSTMs which generates smoother predicted articulatory trajectories. For modeling duration variations across speaking rates, we deploy attention network, which eliminates the needto align trajectories in different rates using DTW. We performa phoneme specific duration analysis to examine how well duration is transformed using the proposed AstNet. As the range of articulatory motions is correlated with speaking rate, we also analyze amplitude of the transformed articulatory movements at different rates compared to their original counterparts, to examine how well the proposed AstNet predicts the extent of articulatory movements in N2F and N2S. We observe that AstNet could model both duration and extent of articulatory movements better than the existing transformation techniques resulting in more accurate transformed articulatory trajectories.


title: A Linear Algebraic Approach to Model Parallelism in Deep Learning

abstract: Training deep neural networks (DNNs) in large-cluster computing environments is increasingly necessary, as networks grow in size and complexity. Local memory and processing limitations require robust data and model parallelism for crossing compute node boundaries. We propose a linear-algebraic approach to model parallelism in deep learning, which allows parallel distribution of any tensor in the DNN. Rather than rely on automatic differentiation tools, which do not universally support distributed memory parallelism models, we show that parallel data movement operations, e.g., broadcast, sum-reduce, and halo exchange, are linear operators, and by defining the relevant spaces and inner products, we manually develop the adjoint, or backward, operators required for gradient-based training of DNNs. We build distributed DNN layers using these parallel primitives, composed with sequential layer implementations, and demonstrate their application by building and training a distributed DNN using DistDL, a PyTorch and MPI-based distributed deep learning toolkit.


title: Embedding calculus and smooth structures

abstract: We study the dependence of the embedding calculus Taylor tower on the smooth structures of the source and target. We prove that embedding calculus does not distinguish exotic smooth structures in dimension 4, implying a negative answer to a question of Viro. In contrast, we show that embedding calculus does distinguish certain exotic spheres in higher dimensions. As a technical tool of independent interest, we prove an isotopy extension theorem for the limit of the embedding calculus tower, which we use to investigate several further examples.


title: SAveRUNNER: a network-based algorithm for drug repurposing and its application to COVID-19

abstract: The novelty of new human coronavirus COVID-19/SARS-CoV-2 and the lack of effective drugs and vaccines gave rise to a wide variety of strategies employed to fight this worldwide pandemic. Many of these strategies rely on the repositioning of existing drugs that could shorten the time and reduce the cost compared to de novo drug discovery. In this study, we presented a new network-based algorithm for drug repositioning, called SAveRUNNER (Searching off-lAbel dRUg aNd NEtwoRk), which predicts drug-disease associations by quantifying the interplay between the drug targets and the disease-specific proteins in the human interactome via a novel network-based similarity measure that prioritizes associations between drugs and diseases locating in the same network neighborhoods. Specifically, we applied SAveRUNNER on a panel of 14 selected diseases with a consolidated knowledge about their disease-causing genes and that have been found to be related to COVID-19 for genetic similarity, comorbidity, or for their association to drugs tentatively repurposed to treat COVID-19. Focusing specifically on SARS subnetwork, we identified 282 repurposable drugs, including some the most rumored off-label drugs for COVID-19 treatments, as well as a new combination therapy of 5 drugs, actually used in clinical practice. Furthermore, to maximize the efficiency of putative downstream validation experiments, we prioritized 24 potential anti-SARS-CoV repurposable drugs based on their network-based similarity values. These top-ranked drugs include ACE-inhibitors, monoclonal antibodies, and thrombin inhibitors. Finally, our findings were in-silico validated by performing a gene set enrichment analysis, which confirmed that most of the network-predicted repurposable drugs may have a potential treatment effect against human coronavirus infections.


title: The High Inclination Solar Mission

abstract: The High Inclination Solar Mission (HISM) is a concept for an out-of-the-ecliptic mission for observing the Sun and the heliosphere. The mission profile is largely based on the Solar Polar Imager concept: initially spiraling in to a 0.48 AU ecliptic orbit, then increasing the orbital inclination at a rate of $\sim 10$ degrees per year, ultimately reaching a heliographic inclination of $>$75 degrees. The orbital profile is achieved using solar sails derived from the technology currently being developed for the Solar Cruiser mission, currently under development. HISM remote sensing instruments comprise an imaging spectropolarimeter (Doppler imager / magnetograph) and a visible light coronagraph. The in-situ instruments include a Faraday cup, an ion composition spectrometer, and magnetometers. Plasma wave measurements are made with electrical antennas and high speed magnetometers. The $7,000\,\mathrm{m}^2$ sail used in the mission assessment is a direct extension of the 4-quadrant $1,666\,\mathrm{m}^2$ Solar Cruiser design and employs the same type of high strength composite boom, deployment mechanism, and membrane technology. The sail system modelled is spun (~1 rpm) to assure required boom characteristics with margin. The spacecraft bus features a fine-pointing 3-axis stabilized instrument platform that allows full science observations as soon as the spacecraft reaches a solar distance of 0.48 AU.


title: Embedding Directed Graphs in Potential Fields Using FastMap-D

abstract: Embedding undirected graphs in a Euclidean space has many computational benefits. FastMap is an efficient embedding algorithm that facilitates a geometric interpretation of problems posed on undirected graphs. However, Euclidean distances are inherently symmetric and, thus, Euclidean embeddings cannot be used for directed graphs. In this paper, we present FastMap-D, an efficient generalization of FastMap to directed graphs. FastMap-D embeds vertices using a potential field to capture the asymmetry between the pairwise distances in directed graphs. FastMap-D learns a potential function to define the potential field using a machine learning module. In experiments on various kinds of directed graphs, we demonstrate the advantage of FastMap-D over other approaches.


title: The size, shape, density and ring of the dwarf planet Haumea from a stellar occultation

abstract: Among the four known transneptunian dwarf planets, Haumea is an exotic, very elongated, and fast rotating body. In contrast to the other dwarf planets, its size, shape, albedo, and density are not well constrained. Here we report results of a multi-chord stellar occultation, observed on 2017 January 21. Secondary events observed around the main body are consistent with the presence of a ring of opacity 0.5, width 70 km, and radius 2,287$_{-45}^{+75}$ km. The Centaur Chariklo was the first body other than a giant planet to show a ring system and the Centaur Chiron was later found to possess something similar to Chariklo's rings. Haumea is the first body outside the Centaur population with a ring. The ring is coplanar with both Haumea's equator and the orbit of its satellite Hi'iaka. Its radius places close to the 3:1 mean motion resonance with Haumea's spin period. The occultation by the main body provides an instantaneous elliptical limb with axes 1,704 $\pm$ 4 km x 1,138 $\pm$ 26 km. Combined with rotational light-curves, it constrains Haumea's 3D orientation and its triaxial shape, which is inconsistent with a homogeneous body in hydrostatic equilibrium. Haumea's largest axis is at least 2,322 $\pm$ 60 km, larger than thought before. This implies an upper limit of 1,885 $\pm$ 80 kg m$^{-3}$ for Haumea's density, smaller and less puzzling than previous estimations, and a geometric albedo of 0.51 $\pm$ 0.02, also smaller than previous estimations. No global N$_2$ or CH$_4$ atmosphere with pressures larger than 15 and 50 nbar (3-$\sigma$ limits), respectively, is detected.


title: Solar Neutrino Detection Sensitivity in DARWIN via Electron Scattering

abstract: We detail the sensitivity of the liquid xenon (LXe) DARWIN observatory to solar neutrinos via elastic electron scattering. We find that DARWIN will have the potential to measure the fluxes of five solar neutrino components: $pp$, $^7$Be, $^{13}$N, $^{15}$O and $pep$. The precision of the $^{13}$N, $^{15}$O and $pep$ components is hindered by the double-beta decay of $^{136}$Xe and, thus, would benefit from a depleted target. A high-statistics observation of $pp$ neutrinos would allow us to infer the values of the weak mixing angle, $\sin^2\theta_w$, and the electron-type neutrino survival probability, $P_e$, in the electron recoil energy region from a few keV up to 200 keV for the first time, with relative precision of 5% and 4%, respectively, at an exposure of 300 ty. An observation of $pp$ and $^7$Be neutrinos would constrain the neutrino-inferred solar luminosity down to 0.2%. A combination of all flux measurements would distinguish between the high (GS98) and low metallicity (AGS09) solar models with 2.1-2.5$\sigma$ significance, independent of external measurements from other experiments or a measurement of $^8$B neutrinos through coherent elastic neutrino-nucleus scattering in DARWIN. Finally, we demonstrate that with a depleted target DARWIN may be sensitive to the neutrino capture process of $^{131}$Xe.

categories: physics.ins-det astro-ph.SR hep-ex hep-ph


title: Elements with north-south dynamics and free subgroups have positive density in Thompson's groups T and V

abstract: We show that the density of elements in Thompson's groups T and V which have north-south dynamics acting on the circle is positive with respect to a stratification in terms of size. We show that the fraction of element pairs which generate a free subgroup of rank two is positive in both groups with respect to a stratification based on diagram size, which shows that the associated subgroup spectra of $T$ and $V$ contain free groups.


title: Policy Learning of MDPs with Mixed Continuous/Discrete Variables: A Case Study on Model-Free Control of Markovian Jump Systems

abstract: Markovian jump linear systems (MJLS) are an important class of dynamical systems that arise in many control applications. In this paper, we introduce the problem of controlling unknown (discrete-time) MJLS as a new benchmark for policy-based reinforcement learning of Markov decision processes (MDPs) with mixed continuous/discrete state variables. Compared with the traditional linear quadratic regulator (LQR), our proposed problem leads to a special hybrid MDP (with mixed continuous and discrete variables) and poses significant new challenges due to the appearance of an underlying Markov jump parameter governing the mode of the system dynamics. Specifically, the state of a MJLS does not form a Markov chain and hence one cannot study the MJLS control problem as a MDP with solely continuous state variable. However, one can augment the state and the jump parameter to obtain a MDP with a mixed continuous/discrete state space. We discuss how control theory sheds light on the policy parameterization of such hybrid MDPs. Then we modify the widely used natural policy gradient method to directly learn the optimal state feedback control policy for MJLS without identifying either the system dynamics or the transition probability of the switching parameter. We implement the (data-driven) natural policy gradient method on different MJLS examples. Our simulation results suggest that the natural gradient method can efficiently learn the optimal controller for MJLS with unknown dynamics.


title: Counting Cards: Exploiting Variance and Data Distributions for Robust Compute In-Memory

abstract: Compute in-memory (CIM) is a promising technique that minimizes data transport, the primary performance bottleneck and energy cost of most data intensive applications. This has found wide-spread adoption in accelerating neural networks for machine learning applications. Utilizing a crossbar architecture with emerging non-volatile memories (eNVM) such as dense resistive random access memory (RRAM) or phase change random access memory (PCRAM), various forms of neural networks can be implemented to greatly reduce power and increase on chip memory capacity. However, compute in-memory faces its own limitations at both the circuit and the device levels. In this work, we explore the impact of device variation and peripheral circuit design constraints. Furthermore, we propose a new algorithm based on device variance and neural network weight distributions to increase both performance and accuracy for compute-in memory based designs. We demonstrate a 27% power improvement and 23% performance improvement for low and high variance eNVM, while satisfying a programmable threshold for a target error tolerance, which depends on the application.


title: Absolute continuity of the harmonic measure on low dimensional rectifiable sets

abstract: We consider a uniformly rectifiable set $\Gamma \subset \mathbb R^n$ of dimension $d<n-1$. By using degenerate elliptic operators on the complement $\Omega = \mathbb R^n \setminus \Gamma$, Guy David, Svitlana Mayboroda, and the author introduced a notion of harmonic measure on $\Gamma$. We prove in the present article that this harmonic measure on $\Gamma$ satisfies the $A^\infty$-property, that is the harmonic measure and the $d$-dimension Hausdorff measure on $\Gamma$ are mutually absolutely continuous in a quantitative and scale invariant way. Thus, we give an alternate proof of a recent theorem of David and Mayboroda, which itself extends a result of Hofmann and Martell to the case where the uniformly rectifiable set $\Gamma$ is not of codimension 1. The proof is surprisingly simple - in particular does not follow the route used by David and Mayboroda, or by Hofmann and Martell - but is specific to the case when $d<n-1$.


title: How individual behaviors drive inequality in online community sizes: an agent-based simulation

abstract: Why are online community sizes so extremely unequal? Most answers to this question have pointed to general mathematical processes drawn from physics like cumulative advantage. These explanations provide little insight into specific social dynamics or decisions that individuals make when joining and leaving communities. In addition, explanations in terms of cumulative advantage do not draw from the enormous body of social computing research that studies individual behavior. Our work bridges this divide by testing whether two influential social mechanisms used to explain community joining can also explain the distribution of community sizes. Using agent-based simulations, we evaluate how well individual-level processes of social exposure and decisions based on individual expected benefits reproduce empirical community size data from Reddit. Our simulations contribute to social computing theory by providing evidence that both processes together---but neither alone---generate realistic distributions of community sizes. Our results also illustrate the potential value of agent-based simulation to online community researchers to both evaluate and bridge individual and group-level theories.


title: Dual-Wavelength ALMA Observations of Dust Rings in Protoplanetary Disks

abstract: We present new Atacama Large Millimeter/submillimeter Array (ALMA) observations for three protoplanetary disks in Taurus at 2.9\,mm and comparisons with previous 1.3\,mm data both at an angular resolution of $\sim0.''1$ (15\,au for the distance of Taurus). In the single-ring disk DS Tau, double-ring disk GO Tau, and multiple-ring disk DL Tau, the same rings are detected at both wavelengths, with radial locations spanning from 50 to 120\,au. To quantify the dust emission morphology, the observed visibilities are modeled with a parametric prescription for the radial intensity profile. The disk outer radii, taken as 95\% of the total flux encircled in the model intensity profiles, are consistent at both wavelengths for the three disks. Dust evolution models show that dust trapping in local pressure maxima in the outer disk could explain the observed patterns. Dust rings are mostly unresolved. The marginally resolved ring in DS Tau shows a tentatively narrower ring at the longer wavelength, an observational feature expected from efficient dust trapping. The spectral index ($\alpha_{\rm mm}$) increases outward and exhibits local minima that correspond to the peaks of dust rings, indicative of the changes in grain properties across the disks. The low optical depths ($\tau\sim$0.1--0.2 at 2.9\,mm and 0.2--0.4 at 1.3\,mm) in the dust rings suggest that grains in the rings may have grown to millimeter sizes. The ubiquitous dust rings in protoplanetary disks modify the overall dynamics and evolution of dust grains, likely paving the way towards the new generation of planet formation.


title: Effects of algorithmic flagging on fairness: quasi-experimental evidence from Wikipedia

abstract: Online community moderators often rely on social signals such as whether or not a user has an account or a profile page as clues that users may cause problems. Reliance on these clues can lead to overprofiling bias when moderators focus on these signals but overlook the misbehavior of others. We propose that algorithmic flagging systems deployed to improve the efficiency of moderation work can also make moderation actions more fair to these users by reducing reliance on social signals and making norm violations by everyone else more visible. We analyze moderator behavior in Wikipedia as mediated by RCFilters, a system which displays social signals and algorithmic flags, and estimate the causal effect of being flagged on moderator actions. We show that algorithmically flagged edits are reverted more often, especially those by established editors with positive social signals, and that flagging decreases the likelihood that moderation actions will be undone. Our results suggest that algorithmic flagging systems can lead to increased fairness in some contexts but that the relationship is complex and contingent.


title: SIDU: Similarity Difference and Uniqueness Method for Explainable AI

abstract: A new brand of technical artificial intelligence ( Explainable AI ) research has focused on trying to open up the 'black box' and provide some explainability. This paper presents a novel visual explanation method for deep learning networks in the form of a saliency map that can effectively localize entire object regions. In contrast to the current state-of-the art methods, the proposed method shows quite promising visual explanations that can gain greater trust of human expert. Both quantitative and qualitative evaluations are carried out on both general and clinical data sets to confirm the effectiveness of the proposed method.


title: Semigroups for dynamical processes on metric graphs

abstract: We present the operator semigroups approach to first- and second-order dynamical systems taking place on metric graphs. We briefly survey the existing results and focus on the well-posedness of the problems with standard vertex conditions. Finally, we show two applications to biological models.


title: Annihilation trajectory of defects in smectic-C films

abstract: In a 2D liquid crystal, each topological defect has a topological charge and a characteristic orientation, and hence can be regarded as an oriented particle. Theories predict that the trajectories of annihilating defects depend on their relative orientation. Recently, these predictions have been tested in experiments on smectic-C films. Those experiments find curved trajectories that are similar to the predictions, but the detailed relationship between the defect orientations and the far-field director is different. To understand this difference, we extend the previous theories by adding the effects of elastic anisotropy, and find that it significantly changes the curved trajectories.


title: Upgrade of IAEA recommended data of selected nuclear reactions for production of PET and SPECT isotopes

abstract: An IAEA research project was dedicated to the compilation, evaluation and recommendation of cross-section data for the accelerator production of $^{11}$C, $^{13}$N, $^{15}$O, $^{18}$F, $^{64}$Cu, and $^{124}$I positron-emitting radionuclides clinically used for PET imaging, and for the accelerator production of two gamma emitters, $^{81}$Rb and $^{123}$I, used in SPECT imaging. Cross sections for 19 charged-particle induced reactions that can be employed for radionuclide accelerator production were evaluated including uncertainties. The resulting reference cross-section data were obtained from Pad\'e fits to selected and corrected experimental data, and integral thick target yields were subsequently deduced. Uncertainties in the fitted results were estimated via a Pad\'e least-squares method with the addition of a 4\% assessed systematic uncertainty to the estimated experimental uncertainty. Experimental data were also compared with theoretical predictions available from the TENDL library. All of the numerical reference cross-section data with their corresponding uncertainties and deduced integral thick target yields are available on-line at the IAEA-NDS medical portal https://www-nds.iaea.org/medportal/ and also at the IAEA-NDS web page https://www-nds.iaea.org/medical/.


title: Exact pointwise estimates for polynomial approximation with Hermite interpolation

abstract: We establish best possible pointwise (up to a constant multiple) estimates for approximation, on a finite interval, by polynomials that satisfy finitely many (Hermite) interpolation conditions, and show that these estimates cannot be improved. In particular, we show that {\bf any} algebraic polynomial of degree $n$ approximating a function $f\in C^r(I)$, $I=[-1,1]$, at the classical pointwise rate $\rho_n^r(x) \omega_k(f^{(r)}, \rho_n(x))$, where $\rho_n(x)=n^{-1}\sqrt{1-x^2}+n^{-2}$, and (Hermite) interpolating $f$ and its derivatives up to the order $r$ at a point $x_0\in I$, has the best possible pointwise rate of (simultaneous) approximation of $f$ near $x_0$. Several applications are given.


title: Kinetics of the surface-nucleated transformation of spherical particles and new model for grain-boundary nucleated transformations

abstract: Equations for the transformed volume fraction of a spherical particle with nucleation on its surface, both nonisothermal and isothermal, are derived in the framework of Kolmogorov method adapted for this problem. Characteristic parameters governing the transformation kinetics are determined; the latter is studied with particular emphasis on the Avrami exponent temporal behavior. It is shown that the surface-nucleated transformation qualitatively differs from the bulk-nucleated one at large values of the characteristic parameters due to the one-dimensional radial growth of the new phase occurring after the complete transformation of the surface itself at the early stage of the process. This effect also manifests itself in the considered ensemble of size-distributed particles and in the grain-boundary nucleated transformations. The logarithmic normal distribution inherent for the particles obtained by grinding is employed for numerical calculations and shown to stretch temporally the volume-fraction and Avrami-exponent dependences for the ensemble of identical particles. A new model for grain-boundary nucleated transformations alternative to the Cahn model of random planes is offered; it is based on the ensemble of size-distributed spherical particles with the possibility for a growing nucleus to cross grain boundaries. The kinetics of this process is shown to be governed by the same characteristic parameter, as for a single particle, and qualitatively differs from the Cahn-model one. In particular, the logarithmic volume-fraction plot at large values of the governing parameter ends by a characteristic bend observed on experimental curves for the crystallization of bulk metallic glasses. This peculiarity together with the form of the plot as a whole directly indicates to the grain (polycluster) structure of metallic glasses and nucleation at intercluster boundaries.


title: Zappa-Sz\'{e}p product of a Fell bundle and a groupoid

abstract: We define the Zappa-Sz\'{e}p product of a Fell bundle by a groupoid, which turns out to be a Fell bundle over the Zappa-Sz\'{e}p product of the underlying groupoids. Under certain assumptions, every Fell bundle over the Zappa-Sz\'{e}p product of groupoids arises in this manner. We then study the representation associated with the Zappa-Sz\'{e}p product Fell bundle and show its relation to covariant representations. Finally, we study the associated universal C*-algebra, which turns out to be a C*-blend, generalizing an earlier result about the Zappa-Sz\'{e}p product of groupoid C*-algebras. In the case of discrete groups, the universal C*-algebra of a Fell bundle embeds injectively inside the universal C*-algebra of the Zappa-Sz\'{e}p product Fell bundle.


title: Negative refraction in the relativistic electron gas

abstract: We show that a gas of relativistic electrons is a left-handed material at low frequencies by computing the effective electric permittivity and effective magnetic permeability that appear in Maxwell's equations in terms of the responses appearing in the constitutive relations, and showing that the former are both negative below the {\it same} frequency, which coincides with the zero-momentum frequency of longitudinal plasmons. We also show, by explicit computation, that the photonic mode of the electromagnetic radiation does not dissipate energy, confirming that it propagates in the gas with the speed of light in vacuum, and that the medium is transparent to it. We then combine those results to show that the gas has a negative effective index of refraction $n_{\rm eff}=-1$. We illustrate the consequences of this fact for Snell's law, and for the reflection and transmission coefficients of the gas.


title: A liquid nitrogen cooled superconducting transition edge sensor with ultra-high responsivity and GHz operation speeds

abstract: Photodetectors based on nano-structured superconducting thin films are currently some of the most sensitive quantum sensors and are key enabling technologies in such broad areas as quantum information, quantum computation and radio-astronomy. However, their broader use is held back by the low operation temperatures which require expensive cryostats. Here, we demonstrate a nitrogen cooled superconducting transition edge sensor, which shows orders of magnitude improved performance characteristics of any superconducting detector operated above 77K, with a responsivity of 9.61x10^4 V/W, noise equivalent power of 15.9 fW/Hz-1/2 and operation speeds up to GHz frequencies. It is based on van der Waals heterostructures of the high temperature superconductor Bi2Sr2CaCu2O8, which are shaped into nano-wires with ultra-small form factor. To highlight the versatility of the detector we demonstrate its fabrication and operation on a telecom grade SiN waveguide chip. Our detector significantly relaxes the demands of practical applications of superconducting detectors and displays its huge potential for photonics based quantum applications.


title: Where are the Dangerous Intersections for Pedestrians and Cyclists: A Colocation-Based Approach

abstract: Pedestrians and cyclists are vulnerable road users. They are at greater risk for being killed in a crash than other road users. The percentage of fatal crashes that involve a pedestrian or cyclist is higher than the overall percentage of total trips taken by both modes. Because of this risk, finding ways to minimize problematic street environments is critical. Understanding traffic safety spatial patterns and identifying dangerous locations with significantly high crash risks for pedestrians and cyclists is essential in order to design possible countermeasures to improve road safety. This research develops two indicators for examining spatial correlation patterns between elements of the built environment (intersections) and crashes (pedestrian- or cyclist-involved). The global colocation quotient detects the overall connection in an area while the local colocation quotient identifies the locations of high-risk intersections. To illustrate our approach, we applied the methods to inspect the colocation patterns between pedestrian- or cyclist-vehicle crashes and intersections in Houston, Texas and we identified among many intersections the ones that significantly attract crashes. We also scrutinized those intersections, discussed possible attributes leading to high colocation of crashes and proposed corresponding countermeasures.


title: Earnings Prediction with Deep Learning

abstract: In the financial sector, a reliable forecast the future financial performance of a company is of great importance for investors' investment decisions. In this paper we compare long-term short-term memory (LSTM) networks to temporal convolution network (TCNs) in the prediction of future earnings per share (EPS). The experimental analysis is based on quarterly financial reporting data and daily stock market returns. For a broad sample of US firms, we find that both LSTMs outperform the naive persistent model with up to 30.0% more accurate predictions, while TCNs achieve and an improvement of 30.8%. Both types of networks are at least as accurate as analysts and exceed them by up to 12.2% (LSTM) and 13.2% (TCN).


title: Data-driven fracture mechanics

abstract: We present a new data-driven paradigm for variational brittle fracture mechanics. The fracture-related material modeling assumptions are removed and the governing equations stemming from variational principles are combined with a set of discrete data points, leading to a model-free data-driven method of solution. The solution at a given load step is identified as the point within the data set that best satisfies either the Kuhn-Tucker conditions stemming from the variational fracture problem or global minimization of a suitable energy functional, leading to data-driven counterparts of both the local and the global minimization approaches of variational fracture mechanics. Both formulations are tested on different test configurations with and without noise and for Griffith and R-curve type fracture behavior.


title: Tensor Completion Made Practical

abstract: Tensor completion is a natural higher-order generalization of matrix completion where the goal is to recover a low-rank tensor from sparse observations of its entries. Existing algorithms are either heuristic without provable guarantees, based on solving large semidefinite programs which are impractical to run, or make strong assumptions such as requiring the factors to be nearly orthogonal. In this paper we introduce a new variant of alternating minimization, which in turn is inspired by understanding how the progress measures that guide convergence of alternating minimization in the matrix setting need to be adapted to the tensor setting. We show strong provable guarantees, including showing that our algorithm converges linearly to the true tensors even when the factors are highly correlated and can be implemented in nearly linear time. Moreover our algorithm is also highly practical and we show that we can complete third order tensors with a thousand dimensions from observing a tiny fraction of its entries. In contrast, and somewhat surprisingly, we show that the standard version of alternating minimization, without our new twist, can converge at a drastically slower rate in practice.