title: On the regularity of small symbolic powers of edge ideals of graphs

abstract: Assume that $G$ is a graph with edge ideal $I(G)$ and let $I(G)^{(s)}$ denote the $s$-th symbolic power of $I(G)$. It is proved that for every integer $s\geq 1$, $${\rm reg}(I(G)^{(s+1)})\leq \max\bigg\{{\rm reg}(I(G))+2s, {\rm reg}\big(I(G)^{(s+1)}+I(G)^s\big)\bigg\}.$$As a consequence, we conclude that ${\rm reg}(I(G)^{(2)})\leq {\rm reg}(I(G))+2$, and ${\rm reg}(I(G)^{(3)})\leq {\rm reg}(I(G))+4$. Moreover, it is shown that if for some integer $k\geq 1$, the graph $G$ has no odd cycle of length at most $2k-1$, then ${\rm reg}(I(G)^{(s)})\leq 2s+{\rm reg}(I(G))-2$, for every integer $s\leq k+1$. Finally, it is proven that ${\rm reg}(I(G)^{(s)})=2s$, for $s\in \{2, 3, 4\}$, provided that the complementary graph $\overline{G}$ is chordal.


title: Quantum Approximate Counting, Simplified

abstract: In 1998, Brassard, Hoyer, Mosca, and Tapp (BHMT) gave a quantum algorithm for approximate counting. Given a list of $N$ items, $K$ of them marked, their algorithm estimates $K$ to within relative error $\varepsilon$ by making only $O\left( \frac{1}{\varepsilon}\sqrt{\frac{N}{K}}\right) $ queries. Although this speedup is of "Grover" type, the BHMT algorithm has the curious feature of relying on the Quantum Fourier Transform (QFT), more commonly associated with Shor's algorithm. Is this necessary? This paper presents a simplified algorithm, which we prove achieves the same query complexity using Grover iterations only. We also generalize this to a QFT-free algorithm for amplitude estimation. Related approaches to approximate counting were sketched previously by Grover, Abrams and Williams, Suzuki et al., and Wie (the latter two as we were writing this paper), but in all cases without rigorous analysis.


title: Vortices and composite order in $\mathrm{SU}(N)$ theories coupled to Abelian gauge field

abstract: We consider $\mathrm{SU}(N)$-symmetric Ginzburg-Landau models coupled to non-compact Abelian gauge field focusing on the case $N > 2$ at finite temperature. We show that, at least for sufficiently large gauge-field coupling constants, these models have two phase transitions. The intermediate phase between the symmetric and low-temperature phases is a state with composite neutral order and no Meissner effect. In this neutral phase the system spontaneously breaks only the symmetry associated with phase differences and density differences between components. For $N > 2$, in contrast to the $SU(2)$ case, the neutral state cannot be mapped onto an $\mathrm{O}(M)$ model. We term this state ${\mathbb{C}{P}}^{N-1}$-neutral phase. We also show that while $\mathrm{SU}(N)$-symmetric Ginzburg-Landau models are not superconductors or superfluids in the usual sense, their state in external field at sufficiently low temperature is a vortex lattice.


title: Koszul homology of quotients by edge ideals

abstract: We show that the Koszul homology algebra of a quotient by the edge ideal of a forest is generated by the lowest linear strand. This provides a large class of Koszul algebras whose Koszul homology algebras satisfy this property. We obtain this result by constructing the minimal graded free resolution of a quotient by such an edge ideal via the so called iterated mapping cone construction and using the explicit bases of Koszul homology given by Herzog and Maleki. Using these methods we also recover a result of Roth and Van Tuyl on the graded Betti numbers of quotients of edge ideals of trees.


title: Na-ion diffusion and electrochemical performance of NaVO$_3$ anode in Li/Na batteries

abstract: We study Na ion diffusion and electrochemical performance of NaVO$_3$ (NVO) as anode material in Li/Na--ion batteries with the specific capacity of $\approx$350 mAhg$^{-1}$ at the current density 11~mAg$^{-1}$ after 300 cycles. Remarkably, the capacity retains $\ge$200~mAhg$^{-1}$ even after 400~cycles at 44~mAg$^{-1}$ with Coulombic efficiency $>$99\%. The deduced diffusion coefficient from galvanostatic intermittent titration technique (GITT), electrochemical impedance spectroscopy (EIS) and cyclic voltammetry (CV) measurements for NVO as anode in Li--ion battery is in the range of 10$^{-10}$--10$^{-12}$~cm$^2$s$^{-1}$. In case of Na-ion batteries, the NVO electrode exhibits initial capacity of 385~mAhg$^{-1}$ at 7~mAg$^{-1}$ current rate, but the capacity degradation is relatively faster in subsequent cycles. We find the diffusion coefficient of NVO--Na cells similar to that of NVO--Li. On the other hand, our charge discharge measurements suggest that the overall performance of NVO anode is better in Li--ion battery than Na-ion. Moreover, we use the density functional theory to simulate the energetics of Na vacancy formation in the bulk of the NVO structure, which is found to be 0.88~eV higher than that of the most stable (100) surface. Thus, the Na ion incorporation at the surface of the electrode material is more facile compared to the bulk.


title: Phase transitions in $\mathrm{U}(1)^N$ lattice London models

abstract: The nature of the phase transitions in multicomponent gauge theories with Abelian gauge field is an outstanding open question. While an early renormalization-group based study indicated a continuous phase transition for $N = 1$ and $N > 183$, recently it was demonstrated the transition is first order at $N = 2$. We quantitatively study the dependence on $N$ of the degree of discontinuity of the phase transitions from a completely ordered to a completely disordered state. We demonstrate that the phase transition is discontinuous at least up to $N = 7$. The dependence on $N$ is nontrivial in that the degree of discontinuity grows with $N$ at least for low $N$. Furthermore we demonstrate that when, at increased coupling strength, the phase transitions of the neutral and changed sectors of the model split, the superfluid phase transition is also discontinuous, at least for $N = 3$ and $N = 4$. Via a duality argument, this indicates that van der Waals-type interaction between directed loops may be responsible for the first-order phase transitions in these models.


title: Transfer Learning from Partial Annotations for Whole Brain Segmentation

abstract: Brain MR image segmentation is a key task in neuroimaging studies. It is commonly conducted using standard computational tools, such as FSL, SPM, multi-atlas segmentation etc, which are often registration-based and suffer from expensive computation cost. Recently, there is an increased interest using deep neural networks for brain image segmentation, which have demonstrated advantages in both speed and performance. However, neural networks-based approaches normally require a large amount of manual annotations for optimising the massive amount of network parameters. For 3D networks used in volumetric image segmentation, this has become a particular challenge, as a 3D network consists of many more parameters compared to its 2D counterpart. Manual annotation of 3D brain images is extremely time-consuming and requires extensive involvement of trained experts. To address the challenge with limited manual annotations, here we propose a novel multi-task learning framework for brain image segmentation, which utilises a large amount of automatically generated partial annotations together with a small set of manually created full annotations for network training. Our method yields a high performance comparable to state-of-the-art methods for whole brain segmentation.


title: A calibration method structured on Bayesian Inference of the HCM speed-flow relationship for freeways and multilane highways and a temporal analysis of traffic behavior

abstract: This paper presents a calibration method for the speed-flow model of the HCM 2016 for freeways and multilane highways allied to temporal analysis of traffic stream. The proposed method was developed using a sample of more than one million observations collected by 23 traffic sensors on four highways in the state of S{\~a}o Paulo. The method is structured on Bayesian inference and provided for each model parameters a probability distribution function. The free-flow speed and capacity presented a probability density function that approximates a Normal distribution. The segment in which the speed of traffic stream remain constant with the increase of the traffic flow is lower than described in HCM 2016, being in some cases close to zero. Along with the proposed calibration method an analysis of temporal variation is performed which shows a significant variation in traffic behavior for different periods. The free-flow speed, capacity and breakpoint distributions obtained through monthly and annual calibration were considered equal by means of Kolmogorov-Smirnov test, different for the model calibration coefficient.


title: Flows on flow-admissible signed graphs

abstract: In 1983, Bouchet proposed a conjecture that every flow-admissible signed graph admits a nowhere-zero $6$-flow. Bouchet himself proved that such signed graphs admit nowhere-zero $216$-flows and Zyka further proved that such signed graphs admit nowhere-zero $30$-flows. In this paper we show that every flow-admissible signed graph admits a nowhere-zero 11-flow.


title: $X$-secure $T$-private Information Retrieval from MDS Coded Storage with Byzantine and Unresponsive Servers

abstract: The problem of $X$-secure $T$-private information retrieval from MDS coded storage is studied in this paper, where the user wishes to privately retrieve one out of $K$ independent messages that are distributed over $N$ servers according to an MDS code. It is guaranteed that any group of up to $X$ colluding servers learn nothing about the messages and that any group of up to $T$ colluding servers learn nothing about the identity of desired message. A lower bound of achievable rates is proved by presenting a novel scheme based on \emph{cross-subspace alignment} and a successive decoding with interference cancellation strategy. For large number of messages $(K\rightarrow\infty)$ the achieved rate, which we conjecture to be optimal, improves upon the best known rates previously reported in the literature by Raviv and Karpuk, and generalizes an achievable rate for MDS-TPIR previously found by Freij-Hollanti et al. that is also conjectured to be asymptotically optimal. The setting is then expanded to allow unresponsive and Byzantine servers. Finally, the scheme is applied to find a new lower convex hull of (download, upload) pairs of secure and private distributed matrix multiplication that generalizes, and in certain asymptotic settings strictly improves upon the best known previous results.


title: Fermionic eigenvector moment flow

abstract: We exhibit new functions of the eigenvectors of the Dyson Brownian motion which follow an equation similar to the Bourgade-Yau eigenvector moment flow. These observables can be seen as a Fermionic counterpart to the original (Bosonic) ones. By analyzing both Fermionic and Bosonic observables, we obtain new correlations between eigenvectors. The fluctuations $\sum_{\alpha\in I}u_k(\alpha)^2-{\vert I\vert}/{N}$ decorrelate for distinct eigenvectors as the dimension $N$ grows and an optimal estimate on the partial inner product $\sum_{\alpha\in I}u_k(\alpha)u_\ell(\alpha)$ between two eigenvectors is given. These static results obtained by integrable dynamics are stated for generalized Wigner matrices and should apply to wide classes of mean field models.


title: Boundedness of massless scalar waves on Kerr interior backgrounds

abstract: We consider solutions of the massless scalar wave equation $\Box_g\psi=0$, without symmetry, on fixed subextremal Kerr backgrounds $(\mathcal M, g)$. It follows from previous analyses in the Kerr exterior that for solutions $\psi$ arising from sufficiently regular data on a two ended Cauchy hypersurface, the solution and its derivatives decay suitably fast along the event horizon $\mathcal H^+$. Using the derived decay rate, we show that $\psi$ is in fact uniformly bounded, $|\psi|\leq C$, in the black hole interior up to and including the bifurcate Cauchy horizon $\mathcal C\mathcal H^+$, to which $\psi$ in fact extends continuously. In analogy to our previous paper, [30], on boundedness of solutions to the massless scalar wave equation on fixed subextremal Reissner--Nordstr\"om backgrounds, the analysis depends on weighted energy estimates, commutation by angular momentum operators and application of Sobolev embedding. In contrast to the Reissner--Nordstr\"om case the commutation leads to additional error terms that have to be controlled.


title: A Cantilever Torque Magnetometry Method for the Measurement of Hall Conductivity of Highly Resistive Samples

abstract: We present the first measurements of Hall conductivity utilizing a new torque magnetometry method designed for insulators. A Corbino disk exhibits a magnetic dipole moment proportional to Hall conductivity when voltage is applied across a test material. This magnetic dipole moment can be measured through torque magnetometry. The symmetry of this contactless technique allows for the measurement of Hall conductivity in previously inaccessible materials. Finally, a low-temperature noise bound, the lack of systematic errors on dummy devices, and a measurement of the Hall conductivity of sputtered indium tin oxide demonstrate the efficacy of the technique.


title: Dark matter and no dark matter: On the halo mass of NGC 1052

abstract: The NGC 1052 group, and in particular the discovery of two ultra diffuse galaxies with very low internal velocity dispersions, has been the subject of much attention recently. Here we present radial velocities for a sample of 77 globular clusters associated with NGC 1052 obtained on the Keck telescope. Their mean velocity and velocity dispersion are consistent with that of the host galaxy. Using a simple tracer mass estimator, we infer the enclosed dynamical mass and dark matter fraction of NGC 1052. Extrapolating our measurements with an NFW mass profile we infer a total halo mass of 6.2 ($\pm$0.2) $\times$ 10$^{12}$ M$_{\odot}$. This mass is fully consistent with that expected from the stellar mass--halo mass relation, suggesting that NGC 1052 has a normal dark matter halo mass (i.e. it is not deficient in dark matter in contrast to two ultra diffuse galaxies in the group). We present a phase space diagram showing the galaxies that lie within the projected virial radius (390 kpc) of NGC 1052. Finally, we briefly discuss the two dark matter deficient galaxies (NGC 1052--DF and DF4) and consider whether MOND can account for their low observed internal velocity dispersions.


title: High-Order Langevin Diffusion Yields an Accelerated MCMC Algorithm

abstract: We propose a Markov chain Monte Carlo (MCMC) algorithm based on third-order Langevin dynamics for sampling from distributions with log-concave and smooth densities. The higher-order dynamics allow for more flexible discretization schemes, and we develop a specific method that combines splitting with more accurate integration. For a broad class of $d$-dimensional distributions arising from generalized linear models, we prove that the resulting third-order algorithm produces samples from a distribution that is at most $\varepsilon > 0$ in Wasserstein distance from the target distribution in $O\left(\frac{d^{1/4}}{ \varepsilon^{1/2}} \right)$ steps. This result requires only Lipschitz conditions on the gradient. For general strongly convex potentials with $\alpha$-th order smoothness, we prove that the mixing time scales as $O \left(\frac{d^{1/4}}{\varepsilon^{1/2}} + \frac{d^{1/2}}{\varepsilon^{1/(\alpha - 1)}} \right)$.


title: Descents of unipotent cuspidal representations of finite classical groups

abstract: Inspired by the Gan-Gross-Prasad conjecture and the descent problem for classical groups, in this paper we study the descents of unipotent cuspidal representations of orthogonal and symplectic groups over finite fields.


title: Absorption of Fermionic Dark Matter by Nuclear Targets

abstract: Absorption of fermionic dark matter leads to a range of distinct and novel signatures at dark matter direct detection and neutrino experiments. We study the possible signals from fermionic absorption by nuclear targets, which we divide into two classes of four Fermi operators: neutral and charged current. In the neutral current signal, dark matter is absorbed by a target nucleus and a neutrino is emitted. This results in a characteristically different nuclear recoil energy spectrum from that of elastic scattering. The charged current channel leads to induced $\beta$ decays in isotopes which are stable in vacuum as well as shifts of the kinematic endpoint of $ \beta$ spectra in unstable isotopes. To confirm the possibility of observing these signals in light of other constraints, we introduce UV completions of example higher dimensional operators that lead to fermionic absorption signals and study their phenomenology. Most prominently, dark matter which exhibits fermionic absorption signals is necessarily unstable leading to stringent bounds from indirect detection searches. Nevertheless, we find a large viable parameter space in which dark matter is sufficiently long lived and detectable in current and future experiments.


title: Gibbs sampling for game-theoretic modeling of private network upgrades with distributed generation

abstract: Renewable energy is increasingly being curtailed, due to oversupply or network constraints. Curtailment can be partially avoided by smart grid management, but the long term solution is network reinforcement. Network upgrades, however, can be costly, so recent interest has focused on incentivising private investors to participate in network investments. In this paper, we study settings where a private renewable investor constructs a power line, but also provides access to other generators that pay a transmission fee. The decisions on optimal (and interdependent) renewable capacities built by investors, affect the resulting curtailment and profitability of projects, and can be formulated as a Stackelberg game. Optimal capacities rely jointly on stochastic variables, such as the renewable resource at project location. In this paper, we show how Markov chain Monte Carlo (MCMC) and Gibbs sampling techniques, can be used to generate observations from historic resource data and simulate multiple future scenarios. Finally, we validate and apply our game-theoretic formulation of the investment decision, to a real network upgrade problem in the UK.


title: Fast multipole method for 3-D Laplace equation in layered media

abstract: In this paper, a fast multipole method (FMM) is proposed for 3-D Laplace equation in layered media. The potential due to charges embedded in layered media is decomposed into a free space component and four types of reaction field components, and the latter can be associated with the potential of a polarization source defined for each type. New multipole expansions (MEs) and local expansions (LEs), as well as the multipole to local (M2L) translation operators are derived for the reaction components, based on which the FMMs for reaction components are then proposed. The resulting FMM for charge interactions in layered media is a combination of using the classic FMM for the free space components and the new FMMs for the reaction field components. With the help of a recurrence formula for the run-time computation of the Sommerfeld-type integrals used in M2L translation operators, pre-computations of a large number of tables are avoided. The new FMMs for the reaction components are found to be much faster than the classic FMM for the free space components due to the separation of equivalent polarization charges and the associated target charges by a material interface. As a result, the FMM for potential in layered media costs almost the same as the classic FMM in the free space case. Numerical results validate the fast convergence of the MEs for the reaction components, and the O(N) complexity of the FMM with a given truncation number p for charge interactions in 3-D layered media.


title: NuSTAR Discovery of Dead Quasar Engine in Arp 187

abstract: Recent active galactic nucleus (AGN) and quasar surveys have revealed a population showing rapid AGN luminosity variability by a factor of $\sim10$. Here we present the most drastic AGN luminosity decline by a factor of $\gtrsim 10^{3}$ constrained by a NuSTAR X-ray observation of the nearby galaxy Arp 187, which is a promising "dead" quasar whose current activity seems quiet but whose past activity of $L_\mathrm{bol} \sim 10^{46}$ erg s$^{-1}$ is still observable at a large scale by its light echo. The obtained upper bound of the X-ray luminosity is $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.2$, corresponding to $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 42.5$, indicating an inactive central engine. Even if a putative torus model with $N_\mathrm{H} \sim 1.5 \times 10^{24}$ cm$^{-2}$ is assumed, the strong upper-bound still holds with $\log (L_{\rm 2-10 keV}/{\rm erg} {\rm s}^{-1}) < 41.8$ or $\log (L_\mathrm{bol}/{\rm erg} {\rm s}^{-1}) < 43.1$. Given the expected size of the narrow line region, this luminosity decrease by a factor of $\gtrsim 10^3$ must have occurred within $\lesssim 10^4$ yr. This extremely rapid luminosity/accretion shutdown is puzzling and it requires one burst-like accretion mechanism producing a clear outer boundary for an accretion disk. We raise two possible scenarios realizing such an accretion mechanism: a mass accretion 1) by the tidal disruption of a molecular cloud and/or 2) by the gas depletion as a result of vigorous nuclear starformation after rapid mass inflow to the central engine.


title: Atmospheric Cherenkov Telescopes as a Potential Veto Array for Neutrino Astronomy

abstract: The IceCube Neutrino Observatory has revealed the existence of sources of high-energy astrophysical neutrinos. However, identification of the sources is challenging because astrophysical neutrinos are difficult to separate from the background of atmospheric neutrinos produced in cosmic-ray-induced particle cascades in the atmosphere. The efficient detection of air showers in coincidence with detected neutrinos can greatly reduce those backgrounds and increase the sensitivity of neutrino telescopes. Imaging Air Cherenkov Telescopes (IACTs) are sensitive to gamma-ray-induced (and cosmic-ray-induced) air showers in the 50 GeV to 50 TeV range, and can therefore be used as background-identifiers for neutrino observatories. This paper describes the feasibility of an array of small scale, wide field-of-view, cost-effective IACTs as an air shower veto for neutrino astronomy. A surface array of 250 to 750 telescopes would significantly improve the performance of a cubic kilometer-scale detector like IceCube, at a cost of a few percent of the original investment. The number of telescopes in the array can be optimized based on astronomical and geometrical considerations.


title: Faint rapid red transients from Neutron star -- CO white-dwarf mergers

abstract: Mergers of neutron stars (NS) and white dwarfs (WD) may give rise to observable explosive transient events. We use 3D hydrodynamical (SPH) simulations, as well as 2D hydrodynamical-thermonuclear simulations (using the FLASH AMR code) to model the disruption of CO-WDs by NSs, which produce faint transient events. We post-process the simulations using a large nuclear network and make use of the SuperNu radiation-transfer code to predict the observational signatures and detailed properties of these transients. We calculate the light-curves (LC) and spectra for five models of NS - CO-WD mergers. The small yields of Ni56 (few x 0.001Msun) results in faint, rapidly-evolving reddened transients (RRTs) with B (R) - peak magnitudes of ~ -12 (-13) to ~ -13 (-15), much shorter and fainter than both regular and faint/peculiar type-Ia SNe. We show that the spectra of RRTs share some similarities with rapidly - evolving transients such as SN2010x, though RRTs are significantly fainter, especially in the I/R bands, and show far stronger Si lines. We estimate that the upcoming Large Synoptic Survey Telescope could detect RRTs at a rate of ~ 10 - 70 yr^-1, through observations in the R/I bands.


title: Voracious vortices in cataclysmic variables: II. Evidence for the expansion of accretion disc material beyond the Roche lobe of the accretor in HT Cassiopeia during its 2017 superoutburst

abstract: In our earlier Paper I we showed that the accretion disc radius of the dwarf nova HT Cas in its quiescent state has not changed significantly during many years of observations. It has remained consistently large, close to the tidal truncation radius. This result is inconsistent with the modern understanding of the evolution of the disc radius through an outburst cycle. Here we present a new set of time-resolved spectra of HT Cas obtained in the middle of its 2017 superoutburst. We used Doppler tomography to map emission structures in the system, which we compared with those detected in quiescence. We used solutions of the restricted three-body problem to re-discuss the location of emission structures and the disc size of HT Cas in quiescence. The superoutburst spectrum is similar in appearance to the quiescent spectra, although the strength of most of the emission lines decreased. However, the high-excitation lines significantly strengthened in comparison with the Balmer lines. Many lines show a mix of broad emission and narrow absorption components. H$\alpha$ in superoutburst was much narrower than in quiescence. Other emission lines also narrowed in outburst, but they had not become as narrow as H$\alpha$. Doppler maps of H$\alpha$ in quiescence, and of the H$\beta$ and HeI lines in outburst are dominated by a bright emission arc in the right side of the tomograms, which is located at and even beyond the theoretical truncation limit. However, the bulk of the H$\alpha$ emission in outburst has significantly lower velocities. We show that the accretion disc radius of HT Cas during its superoutburst has become hot but remained the same size as it was in quiescence. Instead, we detected cool gas beyond the Roche lobe of the white dwarf that may have been expelled from the hot disc during the superoutburst.


title: SimpleBounce : a simple package for the false vacuum decay

abstract: We present SimpleBounce, a C++ package for finding the bounce solution for the false vacuum decay. This package is based on a flow equation which is proposed by the author and solves Coleman-Glaser-Martin's reduced problem: the minimization problem of the kinetic energy while fixing the potential energy. The bounce configuration is obtained by a scale transformation of the solution of this problem. For models with 1--8 scalar field(s), the bounce action can be calculated with O(0.1) % accuracy in O(0.1) s. This package is available at https://github.com/rsato64/SimpleBounce.


title: Harnessing symmetry-protected topological order for quantum memories

abstract: Spin chains with symmetry-protected edge modes are promising candidates to realize intrinsically robust physical qubits that can be used for the storage and processing of quantum information. In any experimental realization of such physical systems, weak perturbations in the form of induced interactions and disorder are unavoidable and can be detrimental to the stored information. At the same time, the latter may in fact be beneficial; for instance by deliberately inducing disorder which causes the system to localize. In this work, we explore the potential of using an $XZX$ cluster Hamiltonian to encode quantum information into the local edge modes and comprehensively investigate the influence of both many-body interactions and disorder on their stability over time, adding substance to the narrative that many-body localization may stabilize quantum information. We recover the edge state at each time step, allowing us to reconstruct the quantum channel that captures the locally constrained out of equilibrium time evolution. With this representation in hand, we analyze how well classical and quantum information are preserved over time as a function of disorder and interactions. We find that the performance of the edge qubits varies dramatically between disorder realizations. Whereas some show a smooth decoherence over time, a sizeable fraction are rapidly rendered unusable as memories. We also find that the stability of the classical information -- a precursor for the usefulness of the chain as a quantum memory -- depends strongly on the direction in which the bit is encoded. When employing the chain as a genuine quantum memory, encoded qubits are most faithfully recovered for low interaction and high disorder.


title: Intergroup Contact in the Wild: Characterizing Language Differences between Intergroup and Single-group Members in NBA-related Discussion Forums

abstract: Intergroup contact has long been considered as an effective strategy to reduce prejudice between groups. However, recent studies suggest that exposure to opposing groups in online platforms can exacerbate polarization. To further understand the behavior of individuals who actively engage in intergroup contact in practice, we provide a large-scale observational study of intragroup behavioral differences between members with and without intergroup contact. We leverage the existing structure of NBA-related discussion forums on Reddit to study the context of professional sports. We identify fans of each NBA team as members of a group and trace whether they have intergroup contact. Our results show that members with intergroup contact use more negative and abusive language in their affiliated group than those without such contact, after controlling for activity levels. We further quantify different levels of intergroup contact and show that there may exist nonlinear mechanisms regarding how intergroup contact relates to intragroup behavior. Our findings provide complementary evidence to experimental studies in a novel context and also shed light on possible reasons for the different outcomes in prior studies.


title: SOAR TESS Survey. I: Sculpting of TESS planetary systems by stellar companions

abstract: TESS is finding transiting planet candidates around bright, nearby stars across the entire sky. The large field-of-view, however, results in low spatial resolution, therefore multiple stars contribute to almost every TESS light curve. High-angular resolution imaging can detect the previously unknown companions to planetary candidate hosts that dilute the transit depths, lead to host star ambiguity, and in some cases are the source of false-positive transit signals. We use speckle imaging on SOAR to search for companions to 542 TESS planet candidate hosts in the Southern sky. We provide correction factors for the 117 systems with resolved companions due to photometric contamination. The contamination in TESS due to close binaries is similar to that found in surveys of Kepler planet candidates. For the solar-type population, we find a deep deficit of close binary systems with projected stellar separations less than 100 AU among planet candidate hosts (44 observed binaries compared to 124 expected based on field binary statistics). The close binary suppression among TESS planet candidate hosts is similar to that seen for the more distant Kepler population. We also find a large surplus of the TESS planet candidates in wide binary systems, detected in both SOAR and Gaia DR2 (119 observed binaries compared to 77 expected). These wide binaries host almost exclusively giant planets, however, suggesting orbital migration, caused by perturbations from the stellar companion, may lead to planet-planet scattering and suppress the population of small planets in wide binaries. Both trends are also apparent in the M-dwarf planet candidate hosts.


title: A Measure for Quantum Paths, Gravity and Spacetime Microstructure

abstract: The number of classical paths of a given length, connecting any two events in a (pseudo) Riemannian spacetime is, of course, infinite. It is, however, possible to define a useful, finite, measure $N(x_2,x_1;\sigma)$ for the effective number of quantum paths [of length $\sigma$ connecting two events $(x_1,x_2)$] in an arbitrary spacetime. When $x_2=x_1$, this reduces to $C(x,\sigma)$ giving the measure for closed quantum loops of length $\sigma$ containing an event $x$. Both $N(x_2,x_1;\sigma)$ and $C(x,\sigma)$ are well-defined and depend only on the geometry of the spacetime. Various other physical quantities like, for e.g., the effective Lagrangian, can be expressed in terms of $N(x_2,x_1;\sigma)$. The corresponding measure for the total path length contributed by the closed loops, in a spacetime region $\mathcal{V}$, is given by the integral of $L(\sigma;x) \equiv\sigma C(\sigma;x)$ over $\mathcal{V}$. Remarkably enough $L(0;x) \propto R(x)$, the Ricci scalar; i.e, the measure for the total length contributed by infinitesimal closed loops in a region of spacetime gives us the Einstein-Hilbert action. Its variation, when we vary the metric, can provide a new route towards induced/emergent gravity descriptions. In the presence of a background electromagnetic field, the corresponding expressions for $N(x_2,x_1;\sigma)$ and $C(x,\sigma)$ can be related to the holonomies of the field. The measure $N(x_2,x_1;\sigma)$ can also be used to evaluate a wide class of path integrals for which the action and the measure are arbitrary functions of the path length. As an example, I compute a modified path integral which incorporates the zero-point-length in the spacetime. I also describe several other properties of $N(x_2,x_1;\sigma)$ and outline a few simple applications.


title: An empirical infrared transit spectrum of Earth: opacity windows and biosignatures

abstract: The Atmospheric Chemistry Experiment's Fourier Transform Spectrometer on the SCISAT satellite has been measuring infrared transmission spectra of Earth during Solar occultations since 2004. We use these data to build an infrared transit spectrum of Earth. Regions of low atmospheric opacity, known as windows, are of particular interest, as they permit observations of the planet's lower atmosphere. Even in the absence of clouds or refraction, imperfect transmittance leads to a minimum effective thickness of $h_{min} \approx 4$ km in the 10--12$\mu$m opacity window at a spectral resolution of $R=10^3$. Nonetheless, at $R=10^5$, the maximum transmittance at the surface is around 70%. In principle, one can probe the troposphere of an Earth-like planet via high-dispersion transit spectroscopy in the mid-infrared; in practice aerosols and/or refraction likely make this impossible. We simulate the transit spectrum of an Earth-like planet in the TRAPPIST-1 system. We find that a long-term near-infrared campaign with JWST could readily detect CO$_2$, establishing the presence of an atmosphere. A mid-IR campaign or longer NIR campaign would be more challenging, but in principle could detect H$_2$O and the biosignatures O$_3$ and CH$_4$.


title: Characterizing the Nature of the Unresolved Point Sources in the Galactic Center

abstract: The Galactic Center Excess (GCE) of GeV gamma rays can be explained as a signal of annihilating dark matter or of emission from unresolved astrophysical sources, such as millisecond pulsars. Evidence for the latter is provided by a statistical procedure---referred to as Non-Poissonian Template Fitting (NPTF)---that distinguishes the smooth distribution of photons expected for dark matter annihilation from a "clumpy" photon distribution expected for point sources. In this paper, we perform an extensive study of the NPTF on simulated data, exploring its ability to recover the flux and luminosity function of unresolved sources at the Galactic Center. When astrophysical background emission is perfectly modeled, we find that the NPTF successfully distinguishes between the dark matter and point source hypotheses when either component makes up the entirety of the GCE. When the GCE is a mixture of dark matter and point sources, the NPTF may fail to reconstruct the correct contribution of each component. We further study the impact of mismodeling the Galactic diffuse backgrounds, finding that while a dark matter signal could be attributed to point sources in some outlying cases for the scenarios we consider, the significance of a true point source signal remains robust. Our work enables us to comment on a recent study by Leane and Slatyer (2019) that questions prior NPTF conclusions because the method does not recover an artificial dark matter signal injected on actual Fermi data. We demonstrate that the failure of the NPTF to extract an artificial dark matter signal can be natural when point sources are present in the data---with the effect further exacerbated by the presence of diffuse mismodeling---and does not on its own invalidate the conclusions of the NPTF analysis in the Inner Galaxy.


title: Dust survival rates in clumps passing through the Cas A reverse shock I: results for a range of clump densities

abstract: The reverse shock in the ejecta of core-collapse supernovae is potentially able to destroy newly formed dust material. In order to determine dust survival rates, we have performed a set of hydrodynamic simulations using the grid-based code AstroBEAR in order to model a shock wave interacting with clumpy supernova ejecta. Dust motions and destruction rates were computed using our newly developed external, post-processing code Paperboats, which includes gas drag, grain charging, sputtering and grain-grain collisions. We have determined dust destruction rates for the oxygen-rich supernova remnant Cassiopeia A as a function of initial grain sizes and clump gas density. We found that up to 30 % of the carbon dust mass is able to survive the passage of the reverse shock if the initial grain size distribution is narrow with radii around ~10 - 50 nm for high gas densities, or with radii around ~0.5 - 1.5 ${\mu}$m for low and medium gas densities. Silicate grains with initial radii around 10 - 30 nm show survival rates of up to 40 % for medium and high density contrasts, while silicate material with micron sized distributions is mostly destroyed. For both materials, the surviving dust mass is rearranged into a new size distribution that can be approximated by two components: a power-law distribution of small grains and a log-normal distribution of grains having the same size range as the initial distribution. Our results show that grain-grain collisions and sputtering are synergistic and that grain-grain collisions can play a crucial role in determining the surviving dust budget in supernova remnants.


title: The Frustration of being Odd: How Boundary Conditions can destroy Local Order

abstract: A central tenant in the classification of phases is that boundary conditions cannot affect the bulk properties of a system. In this work, we show striking, yet puzzling, evidence of a clear violation of this assumption. We use the prototypical example of an XYZ chain with no external field in a ring geometry with an odd number of sites and both ferromagnetic and antiferromagnetic interactions. In such a setting, even at finite sizes, we are able to calculate directly the spontaneous magnetizations that are traditionally used as order parameters to characterize the system's phases. When ferromagnetic interactions dominate, we recover magnetizations that in the thermodynamic limit lose any knowledge about the boundary conditions and are in complete agreement with standard expectations. On the contrary, when the system is governed by antiferromagnetic interactions, the magnetizations decay algebraically to zero with the system size and are not staggered, despite the AFM coupling. We term this behavior {\it ferromagnetic mesoscopic magnetization}. Hence, in the antiferromagnetic regime, our results show an unexpected dependence of a local, one--spin expectation values on the boundary conditions, which is in contrast with predictions from the general theory.

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


title: Two-Magnon Bound States in the Kitaev Model in a $[111]$-Field

abstract: It is now well established that the Kitaev honeycomb model in a magnetic field along the $[111]$-direction harbors an intermediate gapless quantum spin liquid (QSL) phase sandwiched between a gapped non-abelian QSL at low fields $H< H_{c1}$ and a partially polarized phase at high fields $H> H_{c2}$. Here, we analyze the low field and high field phases and phase transitions in terms of single- and two-magnon excitations using exact diagonalization (ED) and density matrix renormalization group (DMRG) methods. We find that the energy to create a bound state of two-magnons $\Delta_p$ becomes lower than the energy to create a single spin flip $\Delta_s$ near $H_{c2}$. In the entire Kitaev spin liquid $\Delta_p<\Delta_s$ and both gaps vanish at $H_{c2}$. We make testable predictions for magnon pairing that could be observable in Raman scattering measurements on Kitaev QSL candidate materials.


title: Axions as a probe of solar metals

abstract: If axions or axion-like particles exist and are detected, they will not only extend the standard model of particle physics but will also open a new way to probe their sources. Axion helioscopes aim to detect axions which are produced in the core of the sun. Their spectrum contains information about the solar interior and could in principle help to solve the conflict between high and low metallicity solar models. Using the planned International Axion Observatory (IAXO) as an example, we show that helioscopes could measure the strength of characteristic emission peaks caused by the presence of heavier elements with good precision. In order to determine unambiguously the elemental abundances from this information, an improved modelling of the states of atoms inside the solar plasma is required.


title: Higgs-mediated optical amplification in a non-equilibrium superconductor

abstract: The quest for new functionalities in quantum materials has recently been extended to non-equilibrium states, which are interesting both because they exhibit new physical phenomena and because of their potential for high-speed device applications. Notable advances have been made in the creation of metastable phases and in Floquet engineering under external periodic driving. In the context of non-equilibrium superconductivity, examples have included the generation of transient superconductivity above the thermodynamic transition temperature, the excitation of coherent Higgs mode oscillations, and the optical control of the interlayer phase in cuprates. Here, we propose theoretically a novel non-equilibrium phenomenon, through which a prompt quench from a metal to a transient superconducting state could induce large oscillations of the order parameter amplitude. We argue that this oscillating mode could act as a source of parametric amplification of the incident radiation. We report experimental results on optically driven K$_3$C$_{60}$ that are consistent with these predictions. The effect is found to disappear when the onset of the excitation becomes slower than the Higgs mode period, consistent with the theory proposed here. These results open new possibilities for the use of collective modes in many-body systems to induce non-linear optical effects.


title: Zero-site DMRG and the optimal low-rank correction

abstract: A zero-site density matrix renormalization algorithm (DMRG0) is proposed to minimize the energy of matrix product states (MPS). Instead of the site tensors themselves, we propose to optimize sequentially the "message" tensors between neighbor sites, which contain the singular values of the bipartition. This leads to a local minimization step that is independent of the physical dimension of the site. Conceptually, it separates the optimization and decimation steps in DMRG. Furthermore, we introduce two new global perturbations based on the optimal low-rank correction to the current state, which are used to avoid local minima. They are determined variationally as the MPS closest to the one-step correction of the Lanczos or Jacobi-Davidson eigensolver, respectively. These perturbations mainly decrease the energy and are free of hand-tuned parameters. Compared to existing single-site enrichment proposals, our approach gives similar convergence ratios per sweep while the computations are cheaper by construction. Our methods may be useful in systems with many physical degrees of freedom per lattice site. We test our approach on the periodic Heisenberg spin chain for various spins, and on free electrons on the lattice.


title: On the relation between Migdal effect and dark matter-electron scattering in isolated atoms and semiconductors

abstract: A key strategy for the direct detection of sub-GeV dark matter is to search for small ionization signals. These can arise from dark matter-electron scattering or when the dark matter-nucleus scattering process is accompanied by a "Migdal" electron. We show that the theoretical descriptions of both processes are closely related, which allows for a principal mapping between dark matter-electron and dark matter-nucleus scattering rates once the dark matter interactions with matter are specified. We explore this parametric relationship for noble-liquid targets and, for the first time, provide an estimate of the "Migdal" ionization rate in semiconductors that is based on evaluating a crystal form factor that accounts for the semiconductor band structure. We also present new dark-matter-nucleus scattering limits down to dark matter masses of 500 keV using published data from XENON10, XENON100, and a SENSEI prototype Skipper-CCD. For a dark photon mediator, the dark matter-electron scattering rates dominate over the Migdal rates for dark matter masses below 100 MeV. We also provide projections for proposed experiments with xenon and silicon targets.


title: Neutrino Emissions from Tidal Disruption Remnants

abstract: We study high-energy neutrino emissions from tidal disruption remnants (TDRs) around supermassive black holes. The neutrinos are produced by the decay of charged pions originating in ultrarelativistic protons that are accelerated there. In the standard theory of tidal disruption events (TDEs), there are four distinct phases from debris circularization of stellar debris to super- and sub-Eddington to radiatively inefficient accretion flows (RIAFs). In addition, we consider the magnetically arrested disk (MAD) state in both the super-Eddington accretion and RIAF phases. We find that there are three promising cases to produce neutrino emissions: the super-Eddington accretion phase of the MAD state and the RIAF phases of both the non-MAD and MAD states. In the super-Eddington MAD state, the enhanced magnetic field makes it possible to accelerate the protons to $E_{p,max}~0.35 PeV (M_bh/10^{7.7}M_\odot)^{41/48}$ with the other given appropriate parameters. The neutrino energy is then $E_{\nu,pk}~67 TeV (M_bh/10^{7.7}M_\odot)^{41/48}$ at the peak of the energy spectrum. For $M_bh\gtrsim10^{7.7} M_\odot$, the neutrino light curve is proportional to $t^{-65/24}$, while it follows the standard $t^{-5/3}$ decay rate for $M_bh<10^{7.7} M_\odot$. In both cases, the large luminosity and characteristic light curves diagnose the super-Eddington MAD state in TDEs. In the RIAF phase of the non-MAD state, we find $E_{p, max}~0.45 PeV (M_bh/10^7M_\odot)^{5/3}$ and $E_{\nu,pk}~0.35 PeV (M_bh/10^7M_\odot)^{5/3}$, and its light curve is proportional to $t^{-10/3}$. This indicates that one can identify whether the existing RIAFs are the TDE origin or not. TDRs are potentially a population of hidden neutrino sources invisible in gamma rays.


title: Hubble Space Telescope Observations of Mira Variables in the Type Ia Supernova Host NGC 1559: An Alternative Candle to Measure the Hubble Constant

abstract: We present year-long, near-infrared Hubble Space Telescope WFC3 observations used to search for Mira variables in NGC 1559, the host galaxy of the Type Ia supernova (SN Ia) 2005df. This is the first dedicated search for Miras, highly-evolved low-mass stars, in a SN Ia host and subsequently the first calibration of the SN Ia luminosity using Miras in a role historically played by Cepheids. We identify a sample of 115 O-rich Miras with P < 400 days based on their light curve properties. We find that the scatter in the Mira Period-Luminosity Relation (PLR) is comparable to Cepheid PLRs seen in SN Ia supernova host galaxies. Using a sample of O-rich Miras discovered in NGC 4258 with HST F160W and its maser distance, we measure a distance modulus for NGC 1559 of mu1559 = 31.41 +/- 0.050 (statistical) +/- 0.060 (systematic) mag. Based on the light curve of the normal, well-observed, low-reddening SN 2005df, we obtain a measurement of the fiducial SN Ia absolute magnitude of MB0 = -19.27 +/- 0.13 mag. With the Hubble diagram of SNe Ia we find H0 = 72.7 +/- 4.6 kms-1 Mpc-1. Combining the calibration from the NGC 4258 megamaser and the Large Magellanic Cloud detached eclipsing binaries gives a best value of H0 = 73.3 +/- 4.0 km s-1 Mpc-1. This result is within 1-sigma of the Hubble constant derived using Cepheids and multiple calibrating SNe Ia. This is the first of four expected calibrations of the SN Ia luminosity from Miras which should reduce the error in H0 via Miras to ~3%. In light of the present Hubble tension and JWST, Miras have utility in the extragalactic distance scale to check Cepheid distances or calibrate nearby SNe in early-type host galaxies that would be unlikely targets for Cepheid searches.


title: Fundamental energy requirement of reversible quantum operations

abstract: Landauer's principle asserts that any computation has an unavoidable energy cost that grows proportionally to its degree of logical irreversibility. But even a logically reversible operation, when run on a physical processor that operates on different energy levels, requires energy. Here we quantify this energy requirement, providing upper and lower bounds that coincide up to a constant factor. We derive these bounds from a general quantum resource-theoretic argument, which implies that the initial resource requirement for implementing a unitary operation within an error~$\epsilon$ grows like $1/\sqrt \epsilon$ times the amount of resource generated by the operation. Applying these results to quantum circuits, we find that their energy requirement can, by an appropriate design, be made independent of their time complexity.


title: Time varying polarized gamma-rays from GRB 160821A: evidence for ordered magnetic fields

abstract: GRB 160821A is the third most energetic gamma ray burst observed by the {\it Fermi} gamma-ray space telescope. Based on the observations made by Cadmium Zinc Telluride Imager (CZTI) on board {\it AstroSat}, here we report the most conclusive evidence to date of (i) high linear polarization ($66^{+26}_{-27} \%$; $5.3 \sigma$ detection), and (ii) variation of its polarization angle with time happening twice during the rise and decay phase of the burst at $3.5 \sigma$ and $3.1 \sigma$ detections respectively. All confidence levels are reported for two parameters of interest. These observations strongly suggest synchrotron radiation produced in magnetic field lines which are highly ordered on angular scales of $1/\Gamma$, where $\Gamma $ is the Lorentz factor of the outflow.


title: Characterization of nuclearity for Beurling-Bj\"{o}rck spaces

abstract: We characterize the nuclearity of the Beurling-Bj\"{o}rck spaces $\mathcal{S}^{(\omega)}_{(\eta)}(\mathbb{R}^d)$ and $\mathcal{S}^{\{\omega\}}_{\{\eta\}}(\mathbb{R}^d)$ in terms of the defining weight functions $\omega$ and $\eta$.


title: An Application of CNNs to Time Sequenced One Dimensional Data in Radiation Detection

abstract: A Convolutional Neural Network architecture was used to classify various isotopes of time-sequenced gamma-ray spectra, a typical output of a radiation detection system of a type commonly fielded for security or environmental measurement purposes. A two-dimensional surface (waterfall plot) in time-energy space is interpreted as a monochromatic image and standard image-based CNN techniques are applied. This allows for the time-sequenced aspects of features in the data to be discovered by the network, as opposed to standard algorithms which arbitrarily time bin the data to satisfy the intuition of a human spectroscopist. The CNN architecture and results are presented along with a comparison to conventional techniques. The results of this novel application of image processing techniques to radiation data will be presented along with a comparison to more conventional adaptive methods.


title: Eliminating Left Recursion without the Epsilon

abstract: The standard algorithm to eliminate indirect left recursion takes a preventative approach, rewriting a grammar's rules so that indirect left recursion is no longer possible, rather than eliminating it only as and when it occurs. This approach results in many of the rules being lost, so that the parse trees that result are often devoid of the detail that the BNF was supposed to capture in the first place. Furthermore, the standard algorithm results in exponential blow-up as the BNF is rewritten, making it wholly unworkable in practice. To avoid these pitfalls, we revise the standard algorithm to eliminate direct left recursion and then take a graph-theoretic approach to eliminating indirect left recursion. We also extend the algorithm to rewrite the resultant parse trees in order to recover the parse trees that would have resulted if left recursion had not had to be eliminated in the first place. Therefore, aside from a couple of caveats, our algorithm works not just in theory but also in practice.


title: Regularity of Minimizers of a Tensor-valued Variational Obstacle Problem in Three Dimensions

abstract: Motivated by Ball and Majumdar's modification of Landau-de Gennes model for nematic liquid crystals, we study energy-minimizer $Q$ of a tensor-valued variational obstacle problem in a bounded 3-D domain with prescribed boundary data. The energy functional is designed to blow up as $Q$ approaches the obstacle. Under certain assumptions, especially on blow-up profile of the singular bulk potential, we prove higher interior regularity of $Q$, and show that the contact set of $Q$ is either empty, or small with characterization of its Hausdorff dimension. We also prove boundary partial regularity of the energy-minimizer.


title: Note on Interacting Langevin Diffusions: Gradient Structure and Ensemble Kalman Sampler by Garbuno-Inigo, Hoffmann, Li and Stuart

abstract: An interacting system of Langevin dynamics driven particles has been proposed for sampling from a given posterior density by Garbuno-Inigo, Hoffmann, Li and Stuart in Interacting Langevin Diffusions: Gradient Structure and Ensemble Kalman Sampler (arXiv:1903:08866v2). The proposed formulation is primarily studied from a formal mean-field limit perspective, while the theoretical behaviour under a finite particle size is left as an open problem. In this note we demonstrate that the particle-based covariance interaction term requires a non-trivial correction. We also show that the corrected dynamics samples exactly from the desired posterior provided that the empirical covariance matrix of the particle system remains non-singular and the posterior log-density satisfies the standard Bakry-Emery criterion.


title: The first astrophysical result of HISAKI: a search for the EUV He I lines a massive cool core cluster at z=0.7

abstract: Molecular cold gas and star formation have been observed at centers of cool-core clusters, albeit at a level much smaller than expected from the classic cooling model. Feedback from the supermassive black hole is likely to have prevented hot gas from cooling. However, the exact cooling and heating processes are poorly understood. The missing key piece is the link between the hot gas ($10^7$\,K) and cold gas ($10^3$\,K). Using the extreme ultraviolet spectrometer onboard {\sl Hisaki}, we explore a distant galaxy cluster, RCS2 J232727.6-020437, one of the most massive cool-core clusters with a cooling rate of $400$\,M$_{\odot}$\,yr$^{-1}$. We aim to detect gas at intermediate temperatures ($3\times10^4$\,K) emitting He I$\alpha$ and He I$\beta$ at rest wavelengths of 58.4 nm and 53.7 nm, respectively. Our target resides at $z=0.6986$, for which these He I lines shift away from the absorption of the Galaxy. Our findings show that the amount of $10^{4-5}$\,K gas at the center of this cluster is smaller than expected if cooling there was uninhibited, which demonstrates that feedback both operates and is efficient for massive clusters at these epochs.


title: Geometry dependence in linear interface growth

abstract: The effect of geometry in the statistics of \textit{nonlinear} universality classes for interface growth has been widely investigated in recent years and it is well known to yield a split of them into subclasses. In this work, we investigate this for the \textit{linear} classes of Edwards-Wilkinson (EW) and of Mullins-Herring (MH) in one- and two-dimensions. From comparison of analytical results with extensive numerical simulations of several discrete models belonging to these classes, as well as numerical integrations of the growth equations on substrates of fixed size (flat geometry) or expanding linearly in time (radial geometry), we verify that the height distributions (HDs), the spatial and the temporal covariances are universal, but geometry-dependent. In fact, the HDs are always Gaussian and, when defined in terms of the so-called "KPZ ansatz" $[h \simeq v_{\infty} t + (\Gamma t)^{\beta} \chi]$, their probability density functions $P(\chi)$ have mean null, so that all their cumulants are null, except by their variances, which assume different values in the flat and radial cases. The shape of the (rescaled) covariance curves is analyzed in detail and compared with some existing analytical results for them. Overall, these results demonstrate that the splitting of such university classes is quite general, being not restricted to the nonlinear ones.


title: Weak localization measurements of electronic scattering rates in Li-doped epitaxial graphene

abstract: Early experiments on alkali-doped graphene demonstrated that the dopant adatoms modify the conductivity of graphene significantly, as extra carriers enhance conductivity while Coulomb scattering off the adatoms suppresses it. However, conductivity probes the overall scattering rate, so a dominant channel associated with long-range Coulomb scattering will mask weaker short-range channels. We present weak localization measurements of epitaxial graphene with lithium adatoms that separately quantify intra- and intervalley scattering rates, then compare the measurements to tight-binding calculations of expected rates for this system. The intravalley rate is strongly enhanced by Li deposition, consistent with Coulomb scattering off the Li adatoms. A simultaneous enhancement of intervalley scattering is partially explained by extra carriers in the graphene interacting with residual disorder. But differences between measured and calculated rates at high Li coverage may indicate adatom-induced modifications to the band structure that go beyond the applied model. Similar adatom-induced modifications of the graphene bands have recently been observed in ARPES, but a full theoretical understanding of these effects is still in development.


title: The Batalin-Vilkovisky Formalism and the Determinant Line Bundle

abstract: Given a smooth family of massless free fermions parametrized by a base manifold $B$, we show that the (mathematically rigorous) Batalin-Vilkovisky quantization of the observables of this family gives rise to the determinant line bundle for the corresponding family of Dirac operators.


title: Symplectic triangle inequality

abstract: We prove a non-squeezing result for Lagrangian embeddings of the real projective plane into blow-ups of the symplectic ball.


title: Analyzing Customer Feedback for Product Fit Prediction

abstract: One of the biggest hurdles for customers when purchasing fashion online, is the difficulty of finding products with the right fit. In order to provide a better online shopping experience, platforms need to find ways to recommend the right product sizes and the best fitting products to their customers. These recommendation systems, however, require customer feedback in order to estimate the most suitable sizing options. Such feedback is rare and often only available as natural text. In this paper, we examine the extraction of product fit feedback from customer reviews using natural language processing techniques. In particular, we compare traditional methods with more recent transfer learning techniques for text classification, and analyze their results. Our evaluation shows, that the transfer learning approach ULMFit is not only comparatively fast to train, but also achieves highest accuracy on this task. The integration of the extracted information with actual size recommendation systems is left for future work.


title: Linear Noise Approximation of Intensity-Driven Signal Transduction Channels

abstract: Biochemical signal transduction, a form of molecular communication, can be modeled using graphical Markov channels with input-modulated transition rates. Such channel models are strongly non-Gaussian. In this paper we use a linear noise approximation to construct a novel class of Gaussian additive white noise channels that capture essential features of fully- and partially-observed intensity-driven signal transduction. When channel state transitions that are sensitive to the input signal are directly observable, high-frequency information is transduced more efficiently than low-frequency information, and the mutual information rate per bandwidth (spectral efficiency) is significantly greater than when sensitive transitions and observable transitions are disjoint. When both observable and hidden transitions are input-sensitive, we observe a superadditive increase in spectral efficiency.


title: A steganographic approach based on the chaotic fractional map and in the DCT domain

abstract: A steganographic method based on the chaotic fractional map and in the DCT domain is proposed. This method embeds a secret message in some high frequency coefficients of the image using a 128-bit private key and a chaotic fractional map which generate a permutation indicating the positions where the secret bits will be embedded. An experimental work on the validation of the proposed method is also presented, showing performance in imperceptibility, quality, similarity and security analysis of the steganographic system. The proposed algorithm improved the level of imperceptibility and Cachin's security of stego-system analyzed through the values of Peak Signal-to-Noise Ratio (PSNR) and the Relative Entropy (RE).


title: Out the Window: A Crowd-Sourced Dataset for Activity Classification in Security Video

abstract: The Out the Window (OTW) dataset is a crowdsourced activity dataset containing 5,668 instances of 17 activities from the NIST Activities in Extended Video (ActEV) challenge. These videos are crowdsourced from workers on the Amazon Mechanical Turk using a novel scenario acting strategy, which collects multiple instances of natural activities per scenario. Turkers are instructed to lean their mobile device against an upper story window overlooking an outdoor space, walk outside to perform a scenario involving people, vehicles and objects, and finally upload the video to us for annotation. Performance evaluation for activity classification on VIRAT Ground 2.0 shows that the OTW dataset provides an 8.3% improvement in mean classification accuracy, and a 12.5% improvement on the most challenging activities involving people with vehicles.


title: Higgs physics beyond the Standard Model

abstract: We give a brief review of beyond-the-Standard Model (BSM) extensions of the Standard Model (SM) Higgs sector. Going from very simple to more complicated models, our survey includes models with additional scalar singlet fields, the Two Higgs Doublet Model (2HDM) and the Higgs sector of the Minimal Supersymmetric Standard Model (MSSM). We discuss the impact of current experimental results from LHC Higgs searches and measurements and the prospective reach of the LHC in the high-luminosity (HL) phase. We furthermore highlight possible new collider signatures within these models that have not been experimentally probed to date.


title: Synergy between Hund-driven correlations and boson-mediated Superconductivity

abstract: Multiorbital systems such as the iron-based superconductors provide a new avenue to attack the longstanding problem of superconductivity in strongly correlated systems. In this work we study the superconductivity driven by a generic bosonic mechanism in a multiorbital model including the full dynamical electronic correlations induced by the Hubbard U and the Hund's coupling. We show that superconductivity survives much more in a Hund's metal than in an ordinary correlated metal with the same degree of correlation. The crucial role of the redistribution of spectral weight in the Hund's metal reflects also in the enhancement of the orbital-selective character of the superconducting gaps, in agreement with experiments in iron-based superconductors.


title: Densities in certain three-way prime number races

abstract: Let $a_1$, $a_2$, and $a_3$ be distinct reduced residues modulo $q$ satisfying the congruences $a_1^2 \equiv a_2^2 \equiv a_3^2 \pmod q$. We conditionally derive an asymptotic formula, with an error term that has a power savings in $q$, for the logarithmic density of the set of real numbers $x$ for which $\pi(x;q,a_1) > \pi(x;q,a_2) > \pi(x;q,a_3)$. The relationship among the $a_i$ allows us to normalize the error terms for the $\pi(x;q,a_i)$ in an atypical way that creates mutual independence among their distributions, and also allows for a proof technique that uses only elementary tools from probability.


title: Compressive Sampling for Array Cameras

abstract: While design of high performance lenses and image sensors has long been the focus of camera development, the size, weight and power of image data processing components is currently the primary barrier to radical improvements in camera resolution. Here we show that Deep-Learning- Aided Compressive Sampling (DLACS) can reduce operating power on camera-head electronics by 20x. Traditional compressive sampling has to date been primarily applied in the physical sensor layer, we show here that with aid from deep learning algorithms, compressive sampling offers unique power management advantages in digital layer compression.


title: Exploring Exoplanet Cloud Assumptions in \textit{JWST} Transmission Spectra

abstract: Clouds are ubiquitous in extrasolar planet atmospheres and are critical to our understanding of planetary climate and chemistry. They also represent one of the greater challenges to overcome when trying to interpret transit transmission spectra of exoplanet atmospheres as their presence can inhibit precise constraints on atmospheric composition and thermal properties. In this work we take a phenomenological approach towards understanding 1) our ability to constrain bulk cloud properties, and 2) the impact of clouds on constraining various atmospheric properties as obtained through transmission spectroscopy with the \textit{James Webb Space Telescope (JWST)}. We do this by exploring retrievals of atmospheric and cloud properties for a generic "hot-Jupiter" as a function of signal-to-noise ratio (SNR), \textit{JWST} observing modes and four different cloud parameterizations. We find that most key atmospheric and cloud inferences can be well constrained in the wavelength range ($\lambda = $ 0.6 - 11 $\mu$m), with NIRCam ($\lambda =$ 2.5 - 5 $\mu$m) being critical in inferring atmospheric properties and NIRISS + MIRI ($\lambda =$ 0.6 - 2.5, 5 - 11 $\mu$m) being necessary for good constraints on cloud parameters. However, constraining the cloud abundance and therefore the total cloud mass requires an observable cloud base in the transit geometry. While higher SNR observations can place tighter constraints on major parameters such as temperature, metallicity and cloud sedimentation, they are unable to eliminate strong degeneracies among cloud parameters. Our investigation of a generic "warm-Neptune" with photochemical haze parameterization also shows promising results in constraining atmospheric and haze properties in the cooler temperature regime.


title: Nanoscratching of metallic glasses -- An atomistic study

abstract: Tribological properties of materials play an important role in engineering applications. Up to now, a number of experimental studies have identified correlations between tribological parameters and the mechanical response. Using molecular dynamics simulations, we study abrasive wear behavior via nanoscratching of a Cu$_{64.5}$Zr$_{35.5}$ metallic glass. The evolution of the normal and transverse forces and hardness values follows the behavior well known for crystalline substrates. In particular, the generation of the frontal pileup weakens the response of the material to the scratching tip and leads to a decrease of the transverse hardness as compared to the normal hardness. However, metallic glasses soften with increasing temperature, particularly above the glass transition temperature thus showing a higher tendency to structurally relax an applied stress. This plastic response is analyzed focusing on local regions of atoms which underwent strong von-Mises strains, since these are the basis of shear-transformation zones and shear bands. The volume occupied by these atoms increases with temperature, but large increases are only observed above the glass transition temperature. We quantify the generation of plasticity by the concept of plastic efficiency, which relates the generation of plastic volume inside the sample with the formation of external damage, viz. the scratch groove. In comparison to nanoindentation, the generation rate of the plastic volume during nanoscratching is significantly temperature dependent making the glass inside more damage-tolerant at lower temperature but more damage-susceptible at elevated temperatures.


title: Sheaves of maximal intersection and multiplicities of stable log maps

abstract: A great number of theoretical results are known about log Gromov-Witten invariants, but few calculations are worked out. In this paper we restrict to surfaces and to genus 0 stable log maps of maximal tangency. We ask how various natural components of the moduli space contribute to the log Gromov-Witten invariants. The first such calculation by Gross-Pandharipande-Siebert deals with multiple covers over rigid curves in the log Calabi-Yau setting. As a natural continuation, in this paper we compute the contributions of non-rigid irreducible curves in the log Calabi-Yau setting and that of the union of two rigid curves in general position. For the former, we construct and study a moduli space of "logarithmic" 1-dimensional sheaves and compare the resulting multiplicity with tropical multiplicity. For the latter, we explicitly describe the components of the moduli space and work out the logarithmic deformation theory in full, which we then compare with the deformation theory of the analogous relative stable maps.


title: DFPENet-geology: A Deep Learning Framework for High Precision Recognition and Segmentation of Co-seismic Landslides

abstract: Automatic recognition and segmentation methods now become the essential requirement in identifying co-seismic landslides, which are fundamental for disaster assessment and mitigation in large-scale earthquakes. This approach used to be carried out through pixel-based or object-oriented methods. However, due to the massive amount of remote sensing data, variations in different earthquake scenarios, and the efficiency requirement for post-earthquake rescue, these methods are difficult to develop into an accurate, rapid, comprehensive, and general (cross-scene) solution for co-seismic landslide recognition. This paper develops a robust model, Dense Feature Pyramid with Encoder-decoder Network (DFPENet), to understand and fuse the multi-scale features of objects in remote sensing images. The proposed method achieves a competitive segmentation accuracy on the public ISPRS 2D Semantic. Furthermore, a comprehensive and widely-used scheme is proposed for co-seismic landslide recognition, which integrates image features extracted from the DFPENet model, geologic features, temporal resolution, landslide spatial analysis, and transfer learning, while only RGB images are used. To corroborate its feasibility and applicability, the proposed scheme is applied to two earthquake-triggered landslides in Jiuzhaigou (China) and Hokkaido (Japan), using available pre- and post-earthquake remote sensing images.


title: Raman scattering study of lattice and magnetic excitations in CrAs

abstract: We report on the lattice and spin dynamics in the monopnictide CrAs using polarized Raman scattering. This system exhibits a first-order magneto-structural phase transition at T_N~ 265 K, below which the magnetic moments of Cr form an incommensurate double helical magnetic structure while the unit cell volume abruptly expands. Across the transition, the frequencies of the fully symmetric Ag phonons strongly renormalize, which along with the results of first-principle calculations suggest the presence of a sizeable coupling of the phonons to the magnetic degrees of freedom in this compound. In addition, we observe two broad modes around 350 and 1700 cm-1 in the magnetic phase, which we associate with two-magnon Raman scattering. This work provides one of the few examples of magnetic Raman scattering from a non-collinear magnet.


title: Interactive Language Learning by Question Answering

abstract: Humans observe and interact with the world to acquire knowledge. However, most existing machine reading comprehension (MRC) tasks miss the interactive, information-seeking component of comprehension. Such tasks present models with static documents that contain all necessary information, usually concentrated in a single short substring. Thus, models can achieve strong performance through simple word- and phrase-based pattern matching. We address this problem by formulating a novel text-based question answering task: Question Answering with Interactive Text (QAit). In QAit, an agent must interact with a partially observable text-based environment to gather information required to answer questions. QAit poses questions about the existence, location, and attributes of objects found in the environment. The data is built using a text-based game generator that defines the underlying dynamics of interaction with the environment. We propose and evaluate a set of baseline models for the QAit task that includes deep reinforcement learning agents. Experiments show that the task presents a major challenge for machine reading systems, while humans solve it with relative ease.


title: Solutions for the Landsberg unicorn problem in Finsler geometry

abstract: It is still a long-standing open problem in Finsler geometry, is there any regular Landsberg metric which is not Berwaldian. However, there are non-regular Landsberg metrics which are not Berwladian. The known examples are established by G. S. Asanov and Z. Shen. In this paper, we use the Maple program to study some explicit examples of non-Berwaldian Landsberg metrics. In fact, such kinds of examples are very tedious and complicated to investigate. Nonetheless, we use the maple program and Finsler packages to simplify the calculations in an elegant way. Depending on these examples, we manage to figure out some geometric properties of the geodesic spray of a non-Berwaldain Landsberg metric. Deforming this spray in a very specific way, using the metrizability tools of the deformed spray, we get new (very simple) non-Berwaldian Landsberg metrics. Moreover, the powerful of this procedure is investigating a very simple and useful formula for the general class obtained by Z. Shen.


title: Loose ends in a strong force 3-body problem

abstract: Up to symmetries, the orbits of three equal masses under an inverse cube force with zero angular momentum and constant moment of inertia can be reparametrized as the geodesics of a complete, negatively curved metric on a pair of pants. The ends of the pants represent binary collisions. Here we will examine the visibility properties of such negatively curved surfaces, allowing a description of orbits beginning or ending in binary collisions of this 3-body problem.


title: Structure and size of the plastic zone formed during nanoindentation of a metallic glass

abstract: Using molecular dynamics simulation, we study the plastic zone created during nanoindentation of a large CuZr glass system. The plastic zone consists of a core region, in which virtually every atom undergoes plastic rearrangement, and a tail, where the density distribution of the plastically active atoms decays to zero. Compared to crystalline substrates, the plastic zone in metallic glasses is significantly smaller than in crystals. The so-called plastic-zone size factor, which relates the radius of the plastic zone to the contact radius of the indenter with the substrate, assumes values around 1, while in crystals -- depending on the crystal structure -- values of 2--3 are common. The small plastic zone in metallic glasses is caused by the essentially homogeneous deformation in the amorphous matrix, while in crystals heterogeneous dislocations prevail, whose growth leads to a marked extension of the plastic zone.


title: Assessment of Multimodal Choice Behavior for Domestic Passengers Using Bayesian Logistic Regressions: A Case Study in China

abstract: This paper investigates the influencing factors in passengers' multimodal traffic choice behaviors and provides a decision-making basis and improvement strategies. By collecting large individual-level data through a comprehensive field survey that was carried out at the major transportation hubs in Xian City in China from the 1st to the 10th of March 2018, we compared 21 variables from the data with four travel modes including air, high-speed rail (HSR), traditional passenger train, and express bus. Among the variables, 12 variables were used as the independent variables after the correlation analysis and the collinearity test, including age, car ownership, ticketing method, travel purpose, travel distance, fare rate, inter-city travel time per hundred kilometers, safety, comfort, punctuality, access time, and departure time. The regression relationships between the travel mode choice and the independent variables were studied using Bayesian binomial logistic regression (BBL) and Bayesian multiple logistic regression (BMNL). The receiver operating characteristic (ROC) curve was applied to analyze the predictive performance of the regression models. The results show that the area under curve (AUC) of the ROC curves for the air travel, high-speed train, passenger train, and express bus using BMNL were 0.9785, 0.9263, 0.8614, and 0.9025, respectively. The AUC value of the ROC curve of the express bus using BBL was 0.8797, which indicated that the passenger travel choice models established by Bayesian logistic regressions had good prediction ability. The application of Bayesian logistic regression in traffic mode choice is beneficial to enriching the assessments of mode-choice behaviors, and the findings of this study can provide a basis for formulating traffic management strategies in practice for shareholders and policymakers.


title: Choosing between incompatible ideals

abstract: Suppose $\mathcal I$ and $\mathcal J$ are proper ideals on some set $X$. We say that $\mathcal I$ and $\mathcal J$ are incompatible if $\mathcal I \cup \mathcal J$ does not generate a proper ideal. Equivalently, $\mathcal I$ and $\mathcal J$ are incompatible if there is some $A \subseteq X$ such that $A \in \mathcal I$ and $X \setminus A \in \mathcal J$. If some $B \subseteq X$ is either in $\mathcal I \setminus \mathcal J$ or in $\mathcal J \setminus \mathcal I$, then we say that $B$ chooses between $\mathcal I$ and $\mathcal J$. We consider the following Ramsey-theoretic problem: Given several pairs $(\mathcal I_1,\mathcal J_1), (\mathcal I_2,\mathcal J_2), \dots, (\mathcal I_k,\mathcal J_k)$ of incompatible ideals on a set $X$, find some $A \subseteq X$ that chooses between as many of these pairs of ideals as possible. The main theorem is that for every $n \in \mathbb N$, there is some $I(n) \in \mathbb N$ such that given at least $I(n)$ pairs of incompatible ideals on any set $X$, there is some $A \subseteq X$ choosing between at least $n$ of them. This theorem is proved in two main steps. The first step is to identify a (purely finitary) problem in extremal combinatorics, and to show that our problem concerning ideals is equivalent to this combinatorial problem. The second step is to analyze the combinatorial problem in order to show that the number $I(n)$ described above exists, and to put bounds on it. We show $\textstyle \frac{1}{2}n \log_2 n - O(n) \,<\, I(n) \,<\, n \ln n + O(n).$ The upper bound is proved by considering a different but closely related combinatorial problem involving hypergraphs, which may be of independent interest. We also investigate some applications of this theorem to a problem concerning conditionally convergent series.


title: A Magnetic White Dwarf with Five H$\alpha$ Components

abstract: G183$-$35 is an unusual white dwarf that shows an H$\alpha$ line split into five components, instead of the usual three components seen in strongly magnetic white dwarfs. Potential explanations for the unusual set of lines includes a double degenerate system containing two magnetic white dwarfs and/or rotational modulation of a complex magnetic field structure. Here we present time-resolved spectroscopy of G183$-$35 obtained at the Gemini Observatory. These data reveal two sets of absorption lines that appear and disappear over a period of about 4 hours. We also detect low-level (0.2%) variability in optical photometry at the same period. We demonstrate that the spectroscopic and photometric variability can be explained by the presence of spots on the surface of the white dwarf and a change in the average field strength from about 4.6 MG to 6.2 MG. The observed variability is clearly due to G183$-$35's relatively short spin period. However, rotational modulation of a complex magnetic field by itself cannot explain the changes seen in the central H$\alpha$ component. An additional source of variability in the line profiles, most likely due to a chemically inhomogeneous surface composition, is also needed. We propose further observations of similar objects to test this scenario.


title: MFGs for partially reversible investment

abstract: This paper analyzes a class of infinite-time-horizon stochastic games with singular controls motivated from the partially reversible problem. It provides an explicit solution for the mean-field game (MFG) and presents sensitivity analysis to compare the solution for the MFG with that for the single-agent control problem. It shows that in the MFG, model parameters not only affect the optimal strategies as in the single-agent case, but also influence the equilibrium price. It then establishes that the solution to the MFG is an $\epsilon$-Nash Equilibrium to the corresponding $N$-player game, with $\epsilon=O\left(\frac{1}{\sqrt N}\right)$.


title: SpatialNLI: A Spatial Domain Natural Language Interface to Databases Using Spatial Comprehension

abstract: A natural language interface (NLI) to databases is an interface that translates a natural language question to a structured query that is executable by database management systems (DBMS). However, an NLI that is trained in the general domain is hard to apply in the spatial domain due to the idiosyncrasy and expressiveness of the spatial questions. Inspired by the machine comprehension model, we propose a spatial comprehension model that is able to recognize the meaning of spatial entities based on the semantics of the context. The spatial semantics learned from the spatial comprehension model is then injected to the natural language question to ease the burden of capturing the spatial-specific semantics. With our spatial comprehension model and information injection, our NLI for the spatial domain, named SpatialNLI, is able to capture the semantic structure of the question and translate it to the corresponding syntax of an executable query accurately. We also experimentally ascertain that SpatialNLI outperforms state-of-the-art methods.


title: Role of the laser chirp parameter in photoionisation

abstract: Nowadays the development of novel particle accelerators is a hot topic for both experimental and theoretical sciences. In the CERN-AWAKE experiment electrons are accelerated in a cold rubidium plasma, generated by a short, intense laser pulse. We describe the corresponding photoionisation process with ab initio quantum mechanical calculations. In this paper we primarily investigate how the frequency chirping of the applied laser pulse influences the photoionisation. Secondarily we also study how the results of the numerical simulations depend on the shape of the envelope of the laser pulse. Based on our results, we give a recommendation to fine-tune the physical parameters of the applied laser pulse in order to improve the homogeneity of the rubidium plasma.


title: Combinatorics of the geometry of Wilson loop diagrams I: equivalence classes via matroids and polytopes

abstract: Wilson loop diagrams are an important tool in studying scattering amplitudes of SYM $N=4$ theory and are known by previous work to be associated to positroids. We characterize the conditions under which two Wilson loop diagrams give the same positroid, prove that an important subclass of subdiagrams (exact subdiagrams) correspond to uniform matroids, and enumerate the number of different Wilson loop diagrams that correspond to each positroid cell. We also give a correspondence between those positroids which can arise from Wilson loop diagrams and directions in associahedra.


title: Deep Learning Theory Review: An Optimal Control and Dynamical Systems Perspective

abstract: Attempts from different disciplines to provide a fundamental understanding of deep learning have advanced rapidly in recent years, yet a unified framework remains relatively limited. In this article, we provide one possible way to align existing branches of deep learning theory through the lens of dynamical system and optimal control. By viewing deep neural networks as discrete-time nonlinear dynamical systems, we can analyze how information propagates through layers using mean field theory. When optimization algorithms are further recast as controllers, the ultimate goal of training processes can be formulated as an optimal control problem. In addition, we can reveal convergence and generalization properties by studying the stochastic dynamics of optimization algorithms. This viewpoint features a wide range of theoretical study from information bottleneck to statistical physics. It also provides a principled way for hyper-parameter tuning when optimal control theory is introduced. Our framework fits nicely with supervised learning and can be extended to other learning problems, such as Bayesian learning, adversarial training, and specific forms of meta learning, without efforts. The review aims to shed lights on the importance of dynamics and optimal control when developing deep learning theory.


title: Dissipation-enabled hydrodynamic conductivity in a tunable bandgap semiconductor

abstract: Electronic transport in the regime where carrier-carrier collisions are the dominant scattering mechanism has taken on new relevance with the advent of ultraclean two-dimensional materials. Here we present a combined theoretical and experimental study of ambipolar hydrodynamic transport in bilayer graphene demonstrating that the conductivity is given by the sum of two Drude-like terms that describe relative motion between electrons and holes, and the collective motion of the electron-hole plasma. As predicted, the measured conductivity of gapless, charge-neutral bilayer graphene is sample- and temperature-independent over a wide range. Away from neutrality, the electron-hole conductivity collapses to a single curve, and a set of just four fitting parameters provides quantitative agreement between theory and experiment at all densities, temperatures, and gaps measured. This work validates recent theories for dissipation-enabled hydrodynamic conductivity and creates a link between semiconductor physics and the emerging field of viscous electronics.


title: Comparing Perturbation Models for Evaluating Stability of Neuroimaging Pipelines

abstract: A lack of software reproducibility has become increasingly apparent in the last several years, calling into question the validity of scientific findings affected by published tools. Reproducibility issues may have numerous sources of error, including the underlying numerical stability of algorithms and implementations employed. Various forms of instability have been observed in neuroimaging, including across operating system versions, minor noise injections, and implementation of theoretically equivalent algorithms. In this paper we explore the effect of various perturbation methods on a typical neuroimaging pipeline through the use of i) targeted noise injections, ii) Monte Carlo Arithmetic, and iii) varying operating systems to identify the quality and severity of their impact. The work presented here demonstrates that even low order computational models such as the connectome estimation pipeline that we used are susceptible to noise. This suggests that stability is a relevant axis upon which tools should be compared, developed, or improved, alongside more commonly considered axes such as accuracy/biological feasibility or performance. The heterogeneity observed across participants clearly illustrates that stability is a property of not just the data or tools independently, but their interaction. Characterization of stability should therefore be evaluated for specific analyses and performed on a representative set of subjects for consideration in subsequent statistical testing. Additionally, identifying how this relationship scales to higher-order models is an exciting next step which will be explored. Finally, the joint application of perturbation methods with post-processing approaches such as bagging or signal normalization may lead to the development of more numerically stable analyses while maintaining sensitivity to meaningful variation.


title: Symmetry group factorization reveals the structure-function relation in the neural connectome of Caenorhabditis elegans

abstract: The neural connectome of the nematode Caenorhabditis elegans has been completely mapped, yet in spite of being one of the smallest connectomes (302 neurons), the design principles that explain how the connectome structure determines its function remain unknown. Here, we find symmetries in the locomotion neural circuit of C. elegans, each characterized by its own symmetry group which can be factorized into the direct product of normal subgroups. The action of these normal subgroups partitions the connectome into sectors of neurons that match broad functional categories. Furthermore, symmetry principles predict the existence of novel finer structures inside these normal subgroups forming feedforward and recurrent networks made of blocks of imprimitivity. These blocks constitute structures made of circulant matrices nested in a hierarchy of block-circulant matrices, whose functionality is understood in terms of neural processing filters responsible for fast processing of information.


title: Solving Math Word Problems with Double-Decoder Transformer

abstract: This paper proposes a Transformer-based model to generate equations for math word problems. It achieves much better results than RNN models when copy and align mechanisms are not used, and can outperform complex copy and align RNN models. We also show that training a Transformer jointly in a generation task with two decoders, left-to-right and right-to-left, is beneficial. Such a Transformer performs better than the one with just one decoder not only because of the ensemble effect, but also because it improves the encoder training procedure. We also experiment with adding reinforcement learning to our model, showing improved performance compared to MLE training.


title: Multimodal Neuroimaging Data Integration and Pathway Analysis

abstract: With fast advancements in technologies, the collection of multiple types of measurements on a common set of subjects is becoming routine in science. Some notable examples include multimodal neuroimaging studies for the simultaneous investigation of brain structure and function, and multi-omics studies for combining genetic and genomic information. Integrative analysis of multimodal data allows scientists to interrogate new mechanistic questions. However, the data collection and generation of integrative hypotheses is outpacing available methodology for joint analysis of multimodal measurements. In this article, we study high-dimensional multimodal data integration in the context of mediation analysis. We aim to understand the roles different data modalities play as possible mediators in the pathway between an exposure variable and an outcome. We propose a mediation model framework with two data types serving as separate sets of mediators, and develop a penalized optimization approach for parameter estimation. We study both the theoretical properties of the estimator through an asymptotic analysis, and its finite-sample performance through simulations. We illustrate our method with a multimodal brain pathway analysis having both structural and functional connectivities as mediators in the association between sex and language processing.


title: Performance Analysis of Zippers

abstract: A zipper is a powerful technique of representing a purely functional data structure in a way that allows fast access to a specific element. It is often used in cases where the imperative data structures would use a mutable pointer. However, the efficiency of zippers as a replacement for mutable pointers is not sufficiently explored. We attempt to address this issue by comparing the performance of zippers and mutable pointers in two common scenarios and three different languages: C++, C#, and Haskell.


title: Impact of anharmonicity on sound wave velocities at extreme conditions

abstract: Theoretical calculations of sound-wave velocities of materials at extreme conditions are of great importance to various fields, in particular geophysics. For example, the seismic data on sound-wave propagation through the solid iron-rich Earth's inner core have been the main source for elucidating its properties and building models. As the laboratory experiments at very high temperatures and pressures are non-trivial, ab initio predictions are invaluable. The latter, however, tend to disagree with experiment. We notice that many attempts to calculate sound-wave velocities of matter at extreme conditions in the framework of quantum-mechanics based methods have not been taking into account the effect of anharmonic atomic vibrations. We show how anharmonic effects can be incorporated into ab initio calculations and demonstrate that in particular they might be non-negligible for iron in Earth's core. Therefore, we open an avenue to reconcile experiment and ab initio theory.


title: Low Cost Sensor Networks; How Do We Know the Data are Reliable?

abstract: Plausibility of data from networks of low-cost measurement devices is a growing and important contentious issue. Informal networks of low-cost devices have particularly come to prominence for air quality monitoring. The contentious point is the believability of data without regular on-site calibration since that is a specialist task and the costs very quickly become very much larger than the cost of installation in the first place. This article suggests that approaches to the problem that involve appropriate use of independent information have the potential to resolve the contention. Ideas are illustrated particularly with reference to low-cost sensor networks for air quality measurement.


title: Physics-Informed Machine Learning Models for Predicting the Progress of Reactive-Mixing

abstract: This paper presents a physics-informed machine learning (ML) framework to construct reduced-order models (ROMs) for reactive-transport quantities of interest (QoIs) based on high-fidelity numerical simulations. QoIs include species decay, product yield, and degree of mixing. The ROMs for QoIs are applied to quantify and understand how the chemical species evolve over time. First, high-resolution datasets for constructing ROMs are generated by solving anisotropic reaction-diffusion equations using a non-negative finite element formulation for different input parameters. Non-negative finite element formulation ensures that the species concentration is non-negative (which is needed for computing QoIs) on coarse computational grids even under high anisotropy. The reactive-mixing model input parameters are a time-scale associated with flipping of velocity, a spatial-scale controlling small/large vortex structures of velocity, a perturbation parameter of the vortex-based velocity, anisotropic dispersion strength/contrast, and molecular diffusion. Second, random forests, F-test, and mutual information criterion are used to evaluate the importance of model inputs/features with respect to QoIs. Third, Support Vector Machines (SVM) and Support Vector Regression (SVR) are used to construct ROMs based on the model inputs. Then, SVR-ROMs are used to predict scaling of QoIs. Qualitatively, SVR-ROMs are able to describe the trends observed in the scaling law associated with QoIs. Fourth, the scaling law's exponent dependence on model inputs/features are evaluated using $k$-means clustering. Finally, in terms of the computational cost, the proposed SVM-ROMs and SVR-ROMs are $\mathcal{O}(10^7)$ times faster than running a high-fidelity numerical simulation for evaluating QoIs.