title: Using Multiple Imputation to Classify Potential Outcomes Subgroups

abstract: With medical tests becoming increasingly available, concerns about over-testing and over-treatment dramatically increase. Hence, it is important to understand the influence of testing on treatment selection in general practice. Most statistical methods focus on average effects of testing on treatment decisions. However, this may be ill-advised, particularly for patient subgroups that tend not to benefit from such tests. Furthermore, missing data are common, representing large and often unaddressed threats to the validity of statistical methods. Finally, it is desirable to conduct analyses that can be interpreted causally. We propose to classify patients into four potential outcomes subgroups, defined by whether or not a patient's treatment selection is changed by the test result and by the direction of how the test result changes treatment selection. This subgroup classification naturally captures the differential influence of medical testing on treatment selections for different patients, which can suggest targets to improve the utilization of medical tests. We can then examine patient characteristics associated with patient potential outcomes subgroup memberships. We used multiple imputation methods to simultaneously impute the missing potential outcomes as well as regular missing values. This approach can also provide estimates of many traditional causal quantities. We find that explicitly incorporating causal inference assumptions into the multiple imputation process can improve the precision for some causal estimates of interest. We also find that bias can occur when the potential outcomes conditional independence assumption is violated; sensitivity analyses are proposed to assess the impact of this violation. We applied the proposed methods to examine the influence of 21-gene assay, the most commonly used genomic test, on chemotherapy selection among breast cancer patients.


title: Unusual Galactic HII Regions at the Intersection of the Central Molecular Zone and the Far Dust Lane

abstract: Sgr E is a massive star formation complex found toward the Galactic center that consists of numerous discrete, compact HII regions. It is located at the intersection between the Central Molecular Zone (CMZ) and the far dust lane of the Galactic bar, similar to "hot spots" seen in external galaxies. Compared with other Galactic star formation complexes, the Sgr E complex is unusual because its HII regions all have similar radio luminosities and angular extents, and they are deficient in ~10micron emission from their photodissociation regions (PDRs). Our Green Bank Telescope (GBT) radio recombination line observations increase the known membership of Sgr E to 19 HII regions. There are 43 additional HII region candidates in the direction of Sgr E, 26 of which are detected for the first time here using MeerKAT 1.28GHz data. Therefore, the true HII region population of Sgr E may number >60. Using APEX SEDIGISM 13CO 2-1 data we discover a 3.0x10^5 Solar mass molecular cloud associated with Sgr E, but find few molecular or far-infrared concentrations at the locations of the Sgr E HII regions. Comparison with simulations and an analysis of its radio continuum properties indicate that Sgr E formed upstream in the far dust lane of the Galactic bar a few Myr ago and will overshoot the CMZ, crashing into the near dust lane. We propose that the unusual infrared properties of the Sgr E HII regions are caused by their orbit about the Galactic center, which have possibly stripped their PDRs.


title: A Perceptually-Motivated Approach for Low-Complexity, Real-Time Enhancement of Fullband Speech

abstract: Over the past few years, speech enhancement methods based on deep learning have greatly surpassed traditional methods based on spectral subtraction and spectral estimation. Many of these new techniques operate directly in the the short-time Fourier transform (STFT) domain, resulting in a high computational complexity. In this work, we propose PercepNet, an efficient approach that relies on human perception of speech by focusing on the spectral envelope and on the periodicity of the speech. We demonstrate high-quality, real-time enhancement of fullband (48 kHz) speech with less than 5% of a CPU core.


title: Physical Modeling and Performance of Spatial-Spectral Filters for CT Material Decomposition

abstract: Material decomposition for imaging multiple contrast agents in a single acquisition has been made possible by spectral CT: a modality which incorporates multiple photon energy spectral sensitivities into a single data collection. This work presents an investigation of a new approach to spectral CT which does not rely on energy-discriminating detectors or multiple x-ray sources. Instead, a tiled pattern of K-edge filters are placed in front of the x-ray to create spatially encoded spectra data. For improved sampling, the spatial-spectral filter is moved continuously with respect to the source. A model-based material decomposition algorithm is adopted to directly reconstruct multiple material densities from projection data that is sparse in each spectral channel. Physical effects associated with the x-ray focal spot size and motion blur for the moving filter are expected to impact overall performance. In this work, those physical effects are modeled and a performance analysis is conducted. Specifically, experiments are presented with simulated focal spot widths between 0.2 mm and 4.0 mm. Additionally, filter motion blur is simulated for a linear translation speeds between 50 mm/s and 450 mm/s. The performance differential between a 0.2 mm and a 1.0 mm focal spot is less than 15% suggesting feasibility of the approach with realistic x-ray tubes. Moreover, for reasonable filter actuation speeds, higher speeds are shown to decrease error (due to improved sampling) despite motion-based spectral blur.


title: Rapid, $B_1$-insensitive, dual-band quasi-adiabatic saturation transfer with optimal control for complete quantification of myocardial ATP flux

abstract: Purpose: Phosphorus saturation-transfer experiments can quantify metabolic fluxes non-invasively. Typically, the forward flux through the creatine-kinase reaction is investigated by observing the decrease in phosphocreatine (PCr) after saturation of $\gamma$-ATP. The quantification of total ATP utilisation is currently under-explored, as it requires simultaneous saturation of inorganic phosphate (Pi) and PCr. This is challenging, as currently available saturation pulses reduce the already-low $\gamma$-ATP signal present. Methods: Using a hybrid optimal-control and Shinnar-Le-Roux method, a quasi-adiabatic RF pulse was designed for the dual-saturation of PCr and Pi to enable determination of total ATP utilisation. The pulses were evaluated in Bloch equation simulations, compared with a conventional hard-cosine DANTE saturation sequence, before application to perfused rat hearts at 11.7 Tesla. Results: The quasi-adiabatic pulse was insensitive to a $>2.5$-fold variation in $B_1$, producing equivalent saturation with a 53% reduction in delivered pulse power and a 33-fold reduction in spillover at the minimum effective $B_1$. This enabled the complete quantification of the synthesis and degradation fluxes for ATP in 30-45 minutes in the perfused rat heart. While the net synthesis flux ($4.24\pm0.8$ mM/s, SEM) was not significantly different from degradation flux ($6.88\pm2$ mM/s, $p=0.06$) and both measures are consistent with prior work, nonlinear error analysis highlights uncertainties in the Pi-to-ATP measurement that may explain a trend suggesting a possible imbalance. Conclusion: This work demonstrates a novel quasi-adiabatic dual-saturation RF pulse with significantly improved performance that can be used to measure ATP turnover in the heart in vivo.


title: Progress on a perimeter surveillance problem

abstract: We consider a perimeter surveillance problem introduced by Kingston, Beard, and Holt in 2008 and studied by Davis, Humphrey, and Kingston in 2019. In this problem, $n$ drones surveil a finite interval, moving at uniform speed and exchanging information only when they meet another drone. Kingston et al. described a particular online algorithm for coordinating their behavior and asked for an upper bound on how long it can take before the drones are fully synchronized. They divided the algorithm's behavior into two phases, and presented upper bounds on the length of each phase based on conjectured worst-case configurations. Davis et al. presented counterexamples to the conjecture for phase 1. We present sharp upper bounds on phase 2 which show that in this case the conjectured worst case is correct. We also present new lower bounds on phase 1 and the total time to synchronization, and report partial progress towards obtaining an upper bound.


title: Efficiently-coupled microring circuit for on-chip cavity QED with trapped atoms

abstract: We present a complete fabrication study of an efficiently-coupled microring optical circuit tailored for cavity quantum electrodynamics (QED) with trapped atoms. The microring structures are fabricated on a transparent membrane with high in-vacuum fiber edge-coupling efficiency in a broad frequency band. In addition, a bus waveguide pulley coupler realizes critical coupling to the microrings at both of the cesium D-line frequencies, while high coupling efficiency is achieved at the cesium 'magic' wavelengths for creating a lattice of two-color evanescent field traps above a microring. The presented platform holds promises for realizing a robust atom-nanophotonics hybrid quantum device.


title: Low-rank tensor reconstruction of concentrated densities with application to Bayesian inversion

abstract: Transport maps have become a popular mechanic to express complicated probability densities using sample propagation through an optimized push-forward. Beside their broad applicability and well-known success, transport maps suffer from several drawbacks such as numerical inaccuracies induced by the optimization process and the fact that sampling schemes have to be employed when quantities of interest, e.g. moments are to compute. This paper presents a novel method for the accurate functional approximation of probability density functions (PDF) that copes with those issues. By interpreting the pull-back result of a target PDF through an inexact transport map as a perturbed reference density, a subsequent functional representation in a more accessible format allows for efficient and more accurate computation of the desired quantities. We introduce a layer-based approximation of the perturbed reference density in an appropriate coordinate system to split the high-dimensional representation problem into a set of independent approximations for which separately chosen orthonormal basis functions are available. This effectively motivates the notion of h- and p-refinement (i.e. ``mesh size'' and polynomial degree) for the approximation of high-dimensional PDFs. To circumvent the curse of dimensionality and enable sampling-free access to certain quantities of interest, a low-rank reconstruction in the tensor train format is employed via the Variational Monte Carlo method. An a priori convergence analysis of the developed approach is derived in terms of Hellinger distance and the Kullback-Leibler divergence. Applications comprising Bayesian inverse problems and several degrees of concentrated densities illuminate the (superior) convergence in comparison to Monte Carlo and Markov-Chain Monte Carlo methods.


title: Data Efficient Voice Cloning from Noisy Samples with Domain Adversarial Training

abstract: Data efficient voice cloning aims at synthesizing target speaker's voice with only a few enrollment samples at hand. To this end, speaker adaptation and speaker encoding are two typical methods based on base model trained from multiple speakers. The former uses a small set of target speaker data to transfer the multi-speaker model to target speaker's voice through direct model update, while in the latter, only a few seconds of target speaker's audio directly goes through an extra speaker encoding model along with the multi-speaker model to synthesize target speaker's voice without model update. Nevertheless, the two methods need clean target speaker data. However, the samples provided by user may inevitably contain acoustic noise in real applications. It's still challenging to generating target voice with noisy data. In this paper, we study the data efficient voice cloning problem from noisy samples under the sequence-to-sequence based TTS paradigm. Specifically, we introduce domain adversarial training (DAT) to speaker adaptation and speaker encoding, which aims to disentangle noise from speech-noise mixture. Experiments show that for both speaker adaptation and encoding, the proposed approaches can consistently synthesize clean speech from noisy speaker samples, apparently outperforming the method adopting state-of-the-art speech enhancement module.


title: On The Development of Multidimensional Progenitor Models For Core-collapse Supernovae

abstract: Multidimensional hydrodynamic simulations of shell convection in massive stars suggest the development of aspherical perturbations that may be amplified during iron core-collapse. These perturbations have a crucial and qualitative impact on the delayed neutrino-driven core-collapse supernova explosion mechanism by increasing the total stress behind the stalled shock. In this paper, we investigate the properties of a 15 \msun model evolved in 1-,2-, and 3-dimensions (3D) for the final $\sim$424 seconds before gravitational instability and iron core-collapse using MESA and the FLASH simulation framework. We find that just before collapse, our initially perturbed fully 3D model reaches angle-averaged convective velocity magnitudes of $\approx$ 240-260 km s$^{-1}$ in the Si- and O-shell regions with a Mach number $\approx$ 0.06. We find the bulk of the power in the O-shell resides at large scales, characterized by spherical harmonic orders ($\ell$) of 2-4, while the Si-shell shows broad spectra on smaller scales of $\ell\approx30-40$. Both convective regions show an increase in power at $\ell=5$ near collapse. We show that the 1D \texttt{MESA} model agrees with the convective velocity profile and speeds of the Si-shell when compared to our highest resolution 3D model. However, in the O-shell region, we find that \texttt{MESA} predicts speeds approximately \emph{four} times slower than all of our 3D models suggest. All eight of the multi-dimensional stellar models considered in this work are publicly available.


title: Robust Validation: Confident Predictions Even When Distributions Shift

abstract: While the traditional viewpoint in machine learning and statistics assumes training and testing samples come from the same population, practice belies this fiction. One strategy -- coming from robust statistics and optimization -- is thus to build a model robust to distributional perturbations. In this paper, we take a different approach to describe procedures for robust predictive inference, where a model provides uncertainty estimates on its predictions rather than point predictions. We present a method that produces prediction sets (almost exactly) giving the right coverage level for any test distribution in an $f$-divergence ball around the training population. The method, based on conformal inference, achieves (nearly) valid coverage in finite samples, under only the condition that the training data be exchangeable. An essential component of our methodology is to estimate the amount of expected future data shift and build robustness to it; we develop estimators and prove their consistency for protection and validity of uncertainty estimates under shifts. By experimenting on several large-scale benchmark datasets, including Recht et al.'s CIFAR-v4 and ImageNet-V2 datasets, we provide complementary empirical results that highlight the importance of robust predictive validity.


title: Voronin universality on the abscissa of absolute convergence

abstract: We prove that the Voronin universality theorem for the Riemann zeta-function extends to the line Re(s)=1 if in addition to vertical shifts we also allow scaling and adding a sufficiently large constant.


title: Nonparametric prediction with spatial data

abstract: We describe a (nonparametric) prediction algorithm for spatial data, based on a canonical factorization of the spectral density function. We provide theoretical results showing that the predictor has desirable asymptotic properties. Finite sample performance is assessed in a Monte Carlo study that also compares our algorithm to a rival nonparametric method based on the infinite AR representation of the dynamics of the data. Finally, we apply our methodology to predict house prices in Los Angeles.


title: Sketching semidefinite programs for faster clustering

abstract: Many clustering problems enjoy solutions by semidefinite programming. Theoretical results in this vein frequently consider data with a planted clustering and a notion of signal strength such that the semidefinite program exactly recovers the planted clustering when the signal strength is sufficiently large. In practice, semidefinite programs are notoriously slow, and so speedups are welcome. In this paper, we show how to sketch a popular semidefinite relaxation of a graph clustering problem known as minimum bisection, and our analysis supports a meta-claim that the clustering task is less computationally burdensome when there is more signal.

categories: cs.IT cs.DS cs.LG math.IT math.ST stat.TH


title: Feasibility of 2D antiscatter grid concept for flat panel detectors: preliminary investigation of primary transmission properties

abstract: Suppressing the effects of scattered radiation in flat panel detector, FPD, based CBCT still remains to be a challenge. To address the scatter problem, we have been investigating the feasibility of a two dimensional antiscatter grid (2D ASG) concept for FPDs. Although a 2D ASG can potentially provide high scatter rejection capability, primary transmission characteristics of a 2D ASG and its implications in image quality plays a more critical role in implementation of the 2D ASG concept. Thus, in this work, a computational model was developed to investigate the primary transmission properties of the 2D ASG for various grid and FPD pixel geometries, and the improvement in signal to noise ratio,SNR, was calculated analytically to demonstrate the impact of 2D ASG's transmission characteristics on image quality. Computational model showed that average primary transmission fraction (Tp) strongly depends on the septal thickness of 2D ASG, and 2D ASG can provide higher Tp than existing radiographic ASGs at a septal thickness of 0.1 mm. Due to the higher Tp, 2D ASG was also predicted to provide SNR improvements in projections in low to moderate scatter environments typically observed in CBCT imaging. On the other hand, the model also indicated that the shadow or footprint of the 2D ASG leads to spatially nonuniform variations in primary signal in FPD pixels. Reduction of septal thickness and optimization of 2D ASG's pitch may play an essential role in reducing such variations in primary image signal, and avoiding potential image artifacts associated with 2D ASG's footprint.


title: Spin-driven evolution of asteroids' top-shapes at fast and slow spins seen from (101955) Bennu and (162173) Ryugu

abstract: Proximity observations by OSIRIS-REx and Hayabusa2 provided clues on the shape evolution processes of the target asteroids, (101955) Bennu and (162173) Ryugu. Their oblate shapes with equatorial ridges, or the so-called top shapes, may have evolved due to their rotational conditions at present and in the past. Different shape evolution scenarios were previously proposed; Bennu's top shape may have been driven by surface processing, while Ryugu's may have been developed due to large deformation. These two scenarios seem to be inconsistent. Here, we revisit the structural analyses in earlier works and fill a gap to connect these explanations. We also apply a semi-analytical technique for computing the cohesive strength distribution in a uniformly rotating triaxial ellipsoid to characterize the global failure of top-shaped bodies. Assuming that the structure is uniform, our semi-analytical approach describes the spatial variations in failed regions at different spin periods; surface regions are the most sensitive at longer spin periods, while interiors fail structurally at shorter spin periods. This finding suggests that the shape evolution of a top shape may vary due to rotation and internal structure, which can explain the different evolution scenarios of Bennu's and Ryugu's top shapes. We interpret our results as the indications of top shapes' various evolution processes.


title: Interfaces in spectral asymptotics and nodal sets

abstract: This is largely a survey of results obtained jointly with Boris Hanin and Peng Zhou on interfaces in spectral asymptotics, both for Schr\"odinger operators on $L^2({\mathbb R}^d)$ and for Toeplitz Hamiltonians acting on holomorphic sections of ample line bundles $L \to M$ over K\"ahler manifolds $(M, \omega)$. By an interface is meant a hypersurface, either in physical space ${\mathbb R}^d$ or in phase space, separating an allowed region where spectral asymptotics are standard and a forbidden region where they are non-standard. The main question is to give the detailed transition between the two types of asymptotics across the hypersurface (i.e. interface). In the real Schr\"odinger setting, the asymptotics are of Airy type; in the K\"ahler setting they are of Erf (Gaussian error function) type. In addition, we introduce the Bargmann-Fock space of a positive Hermitian line bundle and study interface asymptotics in that setting.


title: Two-dimensional critical systems with mixed boundary conditions: Exact Ising results from conformal invariance and boundary-operator expansions

abstract: With conformal-invariance methods, Burkhardt, Guim, and Xue studied the critical Ising model, defined on the upper half plane $y>0$ with different boundary conditions $a$ and $b$ on the negative and positive $x$ axes. For $ab=-+$ and $f+$, they determined the one and two-point averages of the spin $\sigma$ and energy $\epsilon$. Here $+$, $-$, and $f$ stand for spin-up, spin-down, and free-spin boundaries, respectively. The case $+-+-+\dots$, where the boundary conditions switch between $+$ and $-$ at arbitrary points, $\zeta_1$, $\zeta_2$, $\dots$ on the $x$ axis was also analyzed. In this paper the alternating boundary conditions $+f+f+\dots$ and the case $-f+$ of three different boundary conditions are considered. Exact results for the one and two-point averages of $\sigma$, $\epsilon$, and the stress tensor $T$ are derived. Using the results for $\langle T\rangle$, the critical Casimir interaction with the boundary of a wedge-shaped inclusion is analyzed for mixed boundary conditions. The paper also includes a comprehensive discussion of boundary-operator expansions in two-dimensional critical systems with mixed boundary conditions. Two types of expansions - away from switching points of the boundary condition and at switching points - are considered. The asymptotic behavior of two-point averages is expressed in terms of one-point averages with the help of the expansions. We also consider the strip geometry with mixed boundary conditions and derive the distant-wall corrections to one-point averages near one edge due to the other edge using the boundary-operator expansions. The predictions of the boundary-operator expansions are consistent with exact results for Ising systems.


title: Tailoring Sub-micrometer Periodic Surface Structures via Ultrashort Pulsed Direct Laser Interference Patterning

abstract: Direct laser Interference Patterning (DLIP) with ultrashort laser pulses (ULP) represents a precise and fast technique to produce tailored periodic sub-micrometer structures on various materials. In this work, an experimental and theoretical approach is presented to investigate the previously unexplored fundamental mechanisms for the formation of unprecedented laser-induced topographies on stainless steel following proper combinations of DLIP with ULP. DLIP is aimed to determine the initial conditions of the laser-matter interaction by defining an ablated region while double ULP are used to control the reorganisation of the self-assembled laser induced sub-micrometer sized structures by exploiting the interplay of different absorption and excitation levels coupled with the melt hydrodynamics induced by the first of the double pulses. A multiscale physical model is presented to correlate the interference period, polarization orientation and number of incident pulses with the induced morphologies. Special emphasis is given to electron excitation, relaxation processes and hydrodynamical effects that are crucial to the production of complex morphologies. Results are expected to derive new knowledge of laser-matter interaction in combined DLIP and ULP conditions and enable enhanced fabrication capabilities of complex hierarchical sub-micrometer sized structures for a variety of applications.


title: A Bootstrapped Model to Detect Abuse and Intent in White Supremacist Corpora

abstract: Intelligence analysts face a difficult problem: distinguishing extremist rhetoric from potential extremist violence. Many are content to express abuse against some target group, but only a few indicate a willingness to engage in violence. We address this problem by building a predictive model for intent, bootstrapping from a seed set of intent words, and language templates expressing intent. We design both an n-gram and attention-based deep learner for intent and use them as colearners to improve both the basis for prediction and the predictions themselves. They converge to stable predictions in a few rounds. We merge predictions of intent with predictions of abusive language to detect posts that indicate a desire for violent action. We validate the predictions by comparing them to crowd-sourced labelling. The methodology can be applied to other linguistic properties for which a plausible starting point can be defined.


title: SemEval-2020 Task 9: Overview of Sentiment Analysis of Code-Mixed Tweets

abstract: In this paper, we present the results of the SemEval-2020 Task 9 on Sentiment Analysis of Code-Mixed Tweets (SentiMix 2020). We also release and describe our Hinglish (Hindi-English) and Spanglish (Spanish-English) corpora annotated with word-level language identification and sentence-level sentiment labels. These corpora are comprised of 20K and 19K examples, respectively. The sentiment labels are - Positive, Negative, and Neutral. SentiMix attracted 89 submissions in total including 61 teams that participated in the Hinglish contest and 28 submitted systems to the Spanglish competition. The best performance achieved was 75.0% F1 score for Hinglish and 80.6% F1 for Spanglish. We observe that BERT-like models and ensemble methods are the most common and successful approaches among the participants.


title: Lie PCA: Density estimation for symmetric manifolds

abstract: We introduce an extension to local principal component analysis for learning symmetric manifolds. In particular, we use a spectral method to approximate the Lie algebra corresponding to the symmetry group of the underlying manifold. We derive the sample complexity of our method for a variety of manifolds before applying it to various data sets for improved density estimation.


title: Order-preserving dynamics in one dimension -- single-file diffusion and caging from the perspective of dynamical density functional theory

abstract: Dynamical density functional theory (DDFT) is a powerful variational framework to study the nonequilibrium properties of colloids by only considering a time-dependent one-body number density. Despite the large number of recent successes, properly modeling the long-time dynamics in interacting systems within DDFT remains a notoriously difficult problem, since structural information, accounting for temporary or permanent particle cages, gets lost. Here we address such a caging scenario by reducing it to a clean one-dimensional problem, where the particles are naturally ordered (arranged on a line) by perfect cages created by their two next neighbors. In particular, we construct a DDFT approximation based on an equilibrium system with an asymmetric pair potential, such that the corresponding one-body densities still carry the footprint of particle order. Applied to a system of confined hard rods, this order-preserving dynamics (OPD) yields exact results at the system boundaries, in addition to the imprinted correct long-time behavior of density profiles representing individual particles. In an open system, our approach correctly reproduces the reduced long-time diffusion coefficient and subdiffusion, characteristic for a single-file setup. These observations cannot be made using current forms of DDFT without particle order.


title: Numerical Investigation of Critical Electrochemical Factors for Localized Corrosion using a Multi-species Reactive Transport Model

abstract: A multi-species reactive transport model based on the sequential non-iterative approach is employed to investigate two major stages involved in artificial pit experiments of stainless steel: I. Stable pitting under a salt film and II. Film-free dissolution that transitions to repassivation. Data of current density and electrical potential obtained from rapid polarization scans of pits with different depths are utilized to calibrate the elapsed time and the electrode kinetics of each stage. The local chemistry near the base of pits at different temperature and bulk concentrations is simulated to determine several critical electrochemical factors at saturation and repassivation.


title: Maximally Localized Dynamical Quantum Embedding for Solving Many-Body Correlated Systems

abstract: We present a quantum embedding methodology to resolve the Anderson impurity model in the context of dynamical mean-field theory, based on an extended exact diagonalization method. Our method provides a maximally localized quantum impurity model, where the non-local components of the correlation potential remain minimal. This method comes at a large benefit, as the environment used in the quantum embedding approach is described by propagating correlated electrons and hence offers a polynomial increase $O(N^4)$ of the number of degrees of freedom for the embedding mapping without adding bath sites. We report that quantum impurity models with as few as 3 bath sites can reproduce both the Mott transition and the Kondo physics, thus opening a more accessible route to the description of time-dependent phenomena. Finally, we obtain excellent agreement for dynamical magnetic susceptibilities, poising this approach as a candidate to describe 2-particle excitations such as excitons in correlated systems. We expect that our approach will be highly beneficial for the implementation of embedding algorithms on quantum computers, as it allows for a fine description of the correlation in materials with a reduced number of required qubits.


title: The Threshold Strong Dimension of a Graph

abstract: Let $G$ be a connected graph and $u,v$ and $w$ vertices of $G$. Then $w$ is said to {\em strongly resolve} $u$ and $v$, if there is either a shortest $u$-$w$ path that contains $v$ or a shortest $v$-$w$ path that contains $u$. A set $W$ of vertices of $G$ is a {\em strong resolving set} if every pair of vertices of $G$ is strongly resolved by some vertex of $W$. A smallest strong resolving set of a graph is called a {\em strong basis} and its cardinality, denoted $\beta_s(G)$, the {\em strong dimension} of $G$. The {\em threshold strong dimension} of a graph $G$, denoted $\tau_s(G)$, is the smallest strong dimension among all graphs having $G$ as spanning subgraph. A graph whose strong dimension equals its threshold strong dimension is called $\beta_s$-{\em irreducible}. In this paper we establish a geometric characterization for the threshold strong dimension of a graph $G$ that is expressed in terms of the smallest number of paths (each of sufficiently large order) whose strong product admits a certain type of embedding of $G$. We demonstrate that the threshold strong dimension of a graph is not equal to the previously studied threshold dimension of a graph. Graphs with strong dimension $1$ and $2$ are necessarily $\beta_s$-irreducible. It is well-known that the only graphs with strong dimension $1$ are the paths. We completely describe graphs with strong dimension $2$ in terms of the strong resolving graphs introduced by Oellermann and Peters-Fransen. We obtain sharp upper bounds for the threshold strong dimension of general graphs and determine exact values for this invariant for certain subclasses of trees.


title: Transverse Single-Spin Asymmetries of Midrapidity Direct Photons, Neutral Pions, and Eta Mesons in 200 GeV Polarized Proton-Proton Collisions at PHENIX

abstract: Experimental observations of strikingly large transverse single-spin asymmetries (TSSAs) opened a window into quark and gluon dynamics present in hadronic collisions, revealing large spin-momentum correlations within nucleons and in the process of forming hadrons. Though originally measured in lower energy fixed target experiments, they have been found to persist in collisions with momentum transfer well into the perturbative regime of quantum chromodynamics (QCD) and yet their origin remains poorly understood. TSSA measurements have allowed for the development of both transverse momentum dependent and collinear twist-3 descriptions of nonperturbative spin-momentum correlations for both initial- and final-state effects. Results are presented for the TSSAs of direct photons, neutral pions, and eta mesons in the pseudorapidity range $|\eta|<0.35$ from $p^\uparrow+p$ collisions with $\sqrt{s} = 200$ GeV at PHENIX. As hadrons, $\pi^0$ and $\eta$ mesons are sensitive to both initial- and final-state effects. At midrapidity, $\pi^0$ and $\eta$ measurements are sensitive to the dynamics of gluons along with a mix of quark flavors. These results are a factor of three increase in statistical precision and extend to higher transverse momentum when compared with previous PHENIX measurements in this kinematic region. Because direct photon production does not include hadronization, the direct photon TSSA is only sensitive to spin-momentum correlations in the proton. The kinematics of this result in particular make the direct photon TSSA a clean probe of gluon dynamics in the transversely polarized proton. All three of these asymmetries will help constrain the twist-3 trigluon collinear correlation function as well as the gluon Sivers function, improving our knowledge of spin-dependent gluon dynamics in QCD.


title: Time Fused Coefficient SIR Model with Application to COVID-19 Epidemic in the United States

abstract: In this paper, we propose a Susceptible-Infected-Removal (SIR) model with time fused coefficients. In particular, our proposed model discovers the underlying time homogeneity pattern for the SIR model's transmission rate and removal rate via Bayesian shrinkage priors. MCMC sampling for the proposed method is facilitated by the nimble package in R. Extensive simulation studies are carried out to examine the empirical performance of the proposed methods. We further apply the proposed methodology to analyze different levels of COVID-19 data in the United States.


title: COVID19 Tracking: An Interactive Tracking, Visualizing and Analyzing Platform

abstract: The Coronavirus Disease 2019 (COVID-19) has now become a pandemic, inflicting millions of people and causing tens of thousands of deaths. To better understand the dynamics of COVID-19, we present a comprehensive COVID-19 tracking and visualization platform that pinpoints the dynamics of the COVID-19 worldwide. Four essential components are implemented: 1) presenting the visualization map of COVID-19 confirmed cases and total counts all over the world; 2) showing the worldwide trends of COVID-19 at multi-grained levels; 3) provide multi-view comparisons, including confirmed cases per million people, mortality rate and accumulative cure rate; 4) integrating a multi-grained view of the disease spreading dynamics in China and showing how the epidemic is taken under control in China.


title: Parameter estimation in the SIR model from early infections

abstract: A standard model for epidemics is the SIR model on a graph. We introduce a simple algorithm that uses the early infection times from a sample path of the SIR model to estimate the parameters this model, and we provide a performance guarantee in the setting of locally tree-like graphs.


title: J-PAS: Measuring emission lines with artificial neural networks

abstract: Throughout this paper we present a new method to detect and measure emission lines in J-PAS up to $z = 0.35$. J-PAS will observe $8000$~deg$^2$ of the northern sky in the upcoming years with 56 photometric bands. The release of such amount of data brings us the opportunity to employ machine learning methods in order to overcome the difficulties associated with photometric data. We used Artificial Neural Networks (ANNs) trained and tested with synthetic J-PAS photometry from CALIFA, MaNGA, and SDSS spectra. We carry out two tasks: firstly, we cluster galaxies in two groups according to the values of the equivalent width (EW) of $H\alpha$, $H\beta$, $[NII]{\lambda 6584}$, and $ [OIII]{\lambda 5007}$ lines measured in the spectra. Then, we train an ANN to assign to each galaxy a group. We are able to classify them with the uncertainties typical of the photometric redshift measurable in J-PAS. Secondly, we utilize another ANN to determine the values of those EWs. Subsequently, we obtain the $[NII]/H\alpha$, $[OIII]/H\beta$, and \ion{O}{3}\ion{N}{2} ratios recovering the BPT diagram . We study the performance of the ANN in two training samples: one is only composed of synthetic J-PAS photo-spectra (J-spectra) from MaNGA and CALIFA (CALMa set) and the other one is composed of SDSS galaxies. We can reproduce properly the main sequence of star forming galaxies from the determination of the EWs. With the CALMa training set we reach a precision of 0.093 and 0.081 dex for the $[NII]/H\alpha$ and $[OIII]/H\beta$ ratios in the SDSS testing sample. Nevertheless, we find an underestimation of those ratios at high values in galaxies hosting an AGN. We also show the importance of the dataset used for both training and testing the model. ANNs are extremely useful to overcome the limitations previously expected concerning the detection and measurements of the emission lines in surveys like J-PAS.


title: An early warning system for electromagnetic follow-up of gravitational-wave events

abstract: Binary neutron stars (BNSs) will spend $\simeq 10$ -- 15 minutes in the band of Advanced LIGO and Virgo detectors at design sensitivity. Matched-filtering of gravitational-wave (GW) data could in principle accumulate enough signal-to-noise ratio (SNR) to identify a forthcoming event tens of seconds before the companions collide and merge. Here we report on the design and testing of an early warning gravitational-wave detection pipeline. Early warning alerts can be produced for sources that are at low enough redshift so that a large enough SNR accumulates $\sim 10 - 60\,\rm s$ before merger. We find that about 7% (respectively, 49%) of the total detectable BNS mergers will be detected $60\, \rm s$ ($10\, \rm s$) before the merger. About 2% of the total detectable BNS mergers will be detected before merger and localized to within $100\, \rm \text{deg}^2$ (90% credible interval). Coordinated observing by several wide-field telescopes could capture the event seconds before or after the merger. LIGO-Virgo detectors at design sensitivity could facilitate observing at least one event at the onset of merger.


title: $f(R)$ gravity and Tsallis holographic dark energy

abstract: The $f(R)$ gravity theory is considered. It is a gravitational theory that generalizes the Einstein-Hilbert action. In this context, a holographic dark energy model is studied. Tsallis non-extensive entropy is used to introduce the dark energy density based on the holographic principle. Then the Friedmann equation, deceleration parameter, and equation of state for this model are investigated.


title: Geometry of the Gaussian multiplicative chaos in the Wiener space

abstract: We develop an approach for investigating geometric properties of Gaussian multiplicative chaos (GMC) in an infinite dimensional set up. The base space is chosen to be the space of continuous functions endowed with Wiener measure, and the random field is a space-time white noise integrated against Brownian paths. In this set up, we show that in any dimension $d\geq 1$ and for any inverse temperature, the GMC-volume of a ball, uniformly around all paths, decays exponentially with an explicit decay rate. The exponential rate reflects the balance between two competing terms, namely the principal eigenvalue of the Dirichlet Laplacian and an energy functional defined over a certain compactification developed earlier in [MV14]. For $d\geq 3$ and high temperature, the underlying Gaussian field is also shown to attain very high values under the GMC -- that is, all paths are "GMC-thick" in this regime. Both statements are natural infinite dimensional extensions of similar behavior captured by $2d$ Liouville quantum gravity and reflect a certain "atypical behavior" of the GMC: while the GMC volume decays exponentially fast uniformly over all paths, the field itself attains atypically large values on all paths when sampled according to the GMC. It is also shown that, despite the exponential decay of volume for any temperature, for small enough temperature, the normalized overlap of two independent paths tends to follow one of only a finite number of independent paths for most of its allowed time horizon, allowing the GMC probability to accumulate most of its mass along such trajectories.


title: Learning Spectral Templates for Photometric Redshift Estimation from Broadband Photometry

abstract: Estimating redshifts from broadband photometry is often limited by how accurately we can map the colors of galaxies to an underlying spectral template. Current techniques utilize spectrophotometric samples of galaxies or spectra derived from spectral synthesis models. Both of these approaches have their limitations, either the sample sizes are small and often not representative of the diversity of galaxy colors or the model colors can be biased (often as a function of wavelength) which introduces systematics in the derived redshifts. In this paper we learn the underlying spectral energy distributions from an ensemble of $\sim$100K galaxies with measured redshifts and colors. We show that we are able to reconstruct emission and absorption lines at a significantly higher resolution than the broadband filters used to measure the photometry for a sample of 20 spectral templates. We find that our training algorithm reduces the fraction of outliers in the derived photometric redshifts by up to 28%, bias up to 91%, and scatter up to 25%, when compared to estimates using a standard set of spectral templates. We discuss the current limitations of this approach and its applicability for recovering the underlying properties of galaxies. Our derived templates and the code used to produce these results are publicly available in a dedicated Github repository: https://github.com/dirac-institute/photoz_template_learning.


title: Quadric surfaces in the Pfaffian hypersurface in $\mathbb{P}^{14}$

abstract: We study smooth quadric surfaces in the Pfaffian hypersurface in $\mathbb{P}^{14}$ parameterising $6 \times 6$ skew-symmetric matrices of rank at most 4, not intersecting the Grassmannian $\mathbb{G}(1,5)$. Such surfaces correspond to quadratic systems of skew-symmetric matrices of size 6 and constant rank 4, and give rise to a globally generated vector bundle $E$ on the quadric. We analyse these bundles and their geometry, relating them to linear congruences of lines in $\mathbb{P}^5$.


title: Two-stage building energy consumption clustering based on temporal and peak demand patterns

abstract: Analyzing smart meter data to understand energy consumption patterns helps utilities and energy providers perform customized demand response operations. Existing energy consumption segmentation techniques use assumptions that could result in reduced quality of clusters in representing their members. We address this limitation by introducing a two-stage clustering method that more accurately captures load shape temporal patterns and peak demands. In the first stage, load shapes are clustered by allowing a large number of clusters to accurately capture variations in energy use patterns and cluster centroids are extracted by accounting for shape misalignments. In the second stage, clusters of similar centroid and power magnitude range are merged by using Dynamic Time Warping. We used three datasets consisting of ~250 households (~15000 profiles) to demonstrate the performance improvement, compared to baseline methods, and discuss the impact on energy management.


title: Skein theories for virtual tangles

abstract: In this paper, we use skein-theoretic techniques to classify all virtual knot polynomials and trivalent graph invariants with certain smallness conditions. The first half of the paper classifies all virtual knot polynomials giving non-trivial invariants strictly smaller than the one given by the Higman-Sims spin model. In particular, we exhibit a family of skein theories coming from $\text{Rep}(O(2))$ with an interesting braiding. In addition, all skein theories of oriented virtual tangles with some smallness conditions are classified. In the second half of the paper, we classify all non-trivial invariants of (perhaps non-planar) trivalent graphs coming from symmetric trivalent categories. For each of these categories, we also classify when the sub-category generated by only the trivalent vertex is braided. An interesting example of this arise from the tensor category Deligne's $S_t$.


title: Segmentation analysis and the recovery of queuing parameters via the Wasserstein distance: a study of administrative data for patients with chronic obstructive pulmonary disease

abstract: This work uses a data-driven approach to analyse how the resource requirements of patients with chronic obstructive pulmonary disease (COPD) may change, quantifying how those changes impact the hospital system with which the patients interact. This approach is composed of a novel combination of often distinct modes of analysis: segmentation, operational queuing theory, and the recovery of parameters from incomplete data. By combining these methods as presented here, this work demonstrates that potential limitations around the availability of fine-grained data can be overcome. Thus, finding useful operational results despite using only administrative data. The paper begins by finding a useful clustering of the population from this granular data that feeds into a multi-class M/M/c model, whose parameters are recovered from the data via parameterisation and the Wasserstein distance. This model is then used to conduct an informative analysis of the underlying queuing system and the needs of the population under study through several what-if scenarios. The analyses used to form and study this model consider, in effect, all types of patient arrivals and how those types impact the system. With that, this study finds that there are no quick solutions to reduce the impact of COPD patients on the system, including adding capacity to the system. In this analysis, the only effective intervention to reduce the strain caused by those presenting with COPD is to enact external policies which directly improve the overall health of the COPD population before they arrive at the hospital.


title: Concurrent Fixed-Size Allocation and Free in Constant Time

abstract: Our goal is to efficiently solve the dynamic memory allocation problem in a concurrent setting where processes run asynchronously. On $p$ processes, we can support allocation and free for fixed-sized blocks with $O(1)$ worst-case time per operation, $\Theta(p^2)$ additive space overhead, and using only single-word read, write, and CAS. While many algorithms rely on having constant-time fixed-size allocate and free, we present the first implementation of these two operations that is constant time with reasonable space overhead.


title: A residual a posteriori error estimate for the time-domain boundary element method

abstract: This article investigates residual a posteriori error estimates and adaptive mesh refinements for time-dependent boundary element methods for the wave equation. We obtain reliable estimates for Dirichlet and acoustic boundary conditions which hold for a large class of discretizations. Efficiency of the error estimate is shown for a natural discretization of low order. Numerical examples confirm the theoretical results. The resulting adaptive mesh refinement procedures in 3d recover the adaptive convergence rates known for elliptic problems.


title: Probing New Physics in Dimension-8 Neutral Gauge Couplings at $e^+e^-$ Colliders

abstract: Neutral triple gauge couplings (nTGCs) are absent in the standard model effective theory up to dimension-6 operators, but could arise from dimension-8 effective operators. In this work, we study the pure gauge operators of dimension-8 that contribute to nTGCs and are independent of the dimension-8 operator involving the Higgs doublet. We show that the pure gauge operators generate both $Z\gamma Z^*$ and $Z\gamma\gamma^*$ vertices with rapid energy dependence $\propto E^5$, which can be probed sensitively via the reaction $e^+e^- \to Z\gamma$. We demonstrate that measuring the nTGCs via the reaction $e^+e^- \to Z\gamma$ followed by $Z \to q\bar{q}$ decays can probe the new physics scales of dimension-8 pure gauge operators up to the range $(1-5)$TeV at the CEPC, FCC-ee and ILC colliders with $\sqrt{s}=(0.25-1)$TeV, and up to the range $(10-16)$TeV at CLIC with $\sqrt{s}=(3-5)$TeV, assuming in each case an integrated luminosity of 5/ab. We compare these sensitivities with the corresponding probes of the dimension-8 nTGC operators involving Higgs doublets and the dimension-8 fermionic contact operators that contribute to the $e^+e^-Z\gamma$ vertex.


title: A time-dependent FEM-BEM coupling method for fluid-structure interaction in 3d

abstract: We consider the well-posedness and a priori error estimates of a 3d FEM-BEM coupling method for fluid-structure interaction in the time domain. For an elastic body immersed in a fluid, the exterior linear wave equation for the fluid is reduced to an integral equation on the boundary involving the Poincar\'{e}-Steklov operator. The resulting problem is solved using a Galerkin boundary element method in the time domain, coupled to a finite element method for the Lam\'{e} equation inside the elastic body. Based on ideas from the time-independent coupling formulation, we obtain an a priori error estimate and discuss the implementation of the proposed method. Numerical experiments illustrate the performance of our scheme for model problems.


title: On the Equivalence of Three Complete Cyclic Systems of Integers

abstract: The system of coaches by Hilton and Pedersen, the system of cyclic sequences of Schick, and Braendli-Bayne, related to diagonals in regular (2 n)-gons, and the system of modified modular doubling sequences elaborated in this paper are proved to be equivalent. The latter system employs the modified modular equivalence used by Braendli-Bayne. A sequence of Euler tours related on Schick's cycles of diagonals is also presented.


title: Enhanced interlayer interactions in Ni-doped MoS$_2$, and structural and electronic signatures of doping site

abstract: The crystal structure of MoS$_2$ with strong covalent bonds in plane and weak Van der Waals interactions out of plane gives rise to interesting properties for applications such as solid lubrication, optoelectronics, and catalysis, which can be enhanced by transition-metal doping. However, the mechanisms for improvement and even the structure of the doped material can be unclear, which we address with theoretical calculations. Building on our previous work on Ni-doping of the bulk 2H phase, now we compare to polytypes (1H monolayer and 3R bulk), to determine favorable sites for Ni and the doping effect on structure, electronic properties, and the layer dissociation energy. The most favorable intercalation/adatom sites are tetrahedral intercalation for 3R (like 2H) and Mo-atop for 1H. The relative energies indicate a possibility of phase change from 2H to 3R with substitution of Mo or S. We find structural and electronic properties that can be used to identify the doping sites, including metallic behavior in Mo-substituted 3R and 2H, and in-gap states for Mo- and S-substituted 1H, which could have interesting optoelectronic applications. We observe a large enhancement in the interlayer interactions of Ni-doped MoS$_2$, opposite to the effect of other transition metals. For lubrication applications, this increased layer dissociation energy could be the mechanism of low wear. Our systematic study shows the effect of doping concentration and we extrapolate to the low-doping limit. This work gives insight into the previously unclear structure of Ni-doped MoS$_2$ and how it can be detected experimentally, the relation of energy and structures of doped monolayers and bulk systems, the electronic properties under doping, and the effect of doping on interlayer interactions.


title: Gate-free state preparation for fast variational quantum eigensolver simulations: ctrl-VQE

abstract: The variational quantum eigensolver (VQE) is currently the flagship algorithm for solving electronic structure problems on near-term quantum computers. This hybrid quantum/classical algorithm involves implementing a sequence of parameterized gates on quantum hardware to generate a target quantum state, and then measuring the expectation value of the molecular Hamiltonian. Due to finite coherence times and frequent gate errors, the number of gates that can be implemented remains limited on current quantum devices, preventing accurate applications to systems with significant entanglement, such as strongly correlated molecules. In this work, we propose an alternative algorithm (which we refer to as ctrl-VQE) where the quantum circuit used for state preparation is removed entirely and replaced by a quantum control routine which variationally shapes a pulse to drive the initial Hartree-Fock state to the full CI target state. As with VQE, the objective function optimized is the expectation value of the qubit-mapped molecular Hamiltonian. However, by removing the quantum circuit, the coherence times required for state preparation can be drastically reduced by directly optimizing the pulses. We demonstrate the potential of this method numerically by directly optimizing pulse shapes which accurately model the dissociation curves of the hydrogen molecule (covalent bond) and helium hydride ion (ionic bond), and we compute the single point energy for LiH with four transmons.


title: Coloring Fast Without Learning Your Neighbors' Colors

abstract: We give an improved randomized CONGEST algorithm for distance-$2$ coloring that uses $\Delta^2+1$ colors and runs in $O(\log n)$ rounds, improving the recent $O(\log \Delta \cdot \log n)$-round algorithm in [Halld\'orsson, Kuhn, Maus; PODC '20]. We then improve the time complexity to $O(\log \Delta) + 2^{O(\sqrt{\log\log n})}$.


title: A Generalization of the Kramers-Kronig Relations for Linear Time-Varying Media

abstract: We explore the mathematical theory to rigorously describe the response of media with linear time-varying, generally dispersive, electromagnetic constitutive parameters. We show that even when the temporal inhomogeneity takes place on a time scale comparable or shorter than the driving fields' time period, one can still define a physically meaningful time-varying dispersion. Accordingly, a generalized set of Kramers-Kronig relations is investigated to link the real and imaginary parts of the time-varying frequency-dispersive spectra characterizing the medium's constitutive response. Among others, we study the case of a Lorentzian dielectric response with time-varying volumetric density of polarizable atoms and present the varying circuital equivalents of the governing differential equation, which in turn allow us to use the notion of generalized time-varying impedances/admittances of a time-dependent resistor, inductor and capacitor.


title: Experimental Demonstration of a Reconfigurable Coupled Oscillator Platform to Solve the Max-Cut Problem

abstract: In this work, we experimentally demonstrate an integrated circuit (IC) of 30 relaxation oscillators with reconfigurable capacitive coupling to solve the NP-Hard Maximum Cut (Max-Cut) problem. We show that under the influence of an external second-harmonic injection signal, the oscillator phases exhibit a bi-partition which can be used to calculate a high quality approximate Max-Cut solution. Leveraging the all-to-all reconfigurable coupling architecture, we experimentally evaluate the computational properties of the oscillators using randomly generated graph instances of varying size and edge density . Further, comparing the Max-Cut solutions with the optimal values, we show that the oscillators (after simple post-processing) produce a Max-Cut that is within 99% of the optimal value in 28 of the 36 measured graphs; importantly, the oscillators are particularly effective in dense graphs with the Max-Cut being optimal in seven out of nine measured graphs with edge density 0.8. Our work marks a step towards creating an efficient, room-temperature-compatible non-Boolean hardware-based solver for hard combinatorial optimization problems.


title: Instantaneous frequencies in the Kuramoto model

abstract: Using the main results of the Kuramoto theory of globally coupled phase oscillators combined with methods from probability and generalized function theory in a geometric analysis, we extend Kuramoto's results and obtain a mathematical description of the instantaneous frequency (phase-velocity) distribution. Our result is validated against numerical simulations, and we illustrate it in cases where the natural frequencies have normal and Beta distributions. In both cases, we vary the coupling strength and compare systematically the distribution of time-averaged frequencies (a known result of Kuramoto theory) to that of instantaneous frequencies, focussing on their qualitative differences near the synchronized frequency and in their tails. For a class of natural frequency distributions with power-law tails, which includes the Cauchy-Lorentz distribution, we analyze rare events by means of an asymptotic formula obtained from a power series expansion of the instantaneous frequency distribution.


title: Spectrum of twists of Cayley and Cayley sum graphs

abstract: Let $G$ be a finite group with $|G|\geq 4$ and $S$ be a subset of $G$. Given an automorphism $\sigma$ of $G$, the twisted Cayley graph $C(G, S)^\sigma$ (resp. the twisted Cayley sum graph $C_\Sigma(G, S)^\sigma$) is defined as the graph having $G$ as its set of vertices and the adjacent vertices of a vertex $g\in G$ are of the form $\sigma(gs)$ (resp. $\sigma(g^{-1} s)$) for some $s\in S$. If the twisted Cayley graph $C(G, S)^\sigma$ is undirected and connected, then we prove that the nontrivial spectrum of its normalised adjacency operator is bounded away from $-1$ and this bound depends only on its degree, the order of $\sigma$ and the vertex Cheeger constant of $C(G, S)^\sigma$. Moreover, if the twisted Cayley sum graph $C_\Sigma(G, S)^\sigma$ is undirected and connected, then we prove that the nontrivial spectrum of its normalised adjacency operator is bounded away from $-1$ and this bound depends only on its degree and the vertex Cheeger constant of $C_\Sigma(G, S)^\sigma$. We also study these twisted graphs with respect to anti-automorphisms, and obtain similar results. Further, we prove an analogous result for the Schreier graphs satisfying certain conditions.


title: CG-SENSE revisited: Results from the first ISMRM reproducibility challenge

abstract: Purpose: The aim of this work is to shed light on the issue of reproducibility in MR image reconstruction in the context of a challenge. Participants had to recreate the results of "Advances in sensitivity encoding with arbitrary k-space trajectories" by Pruessmann et al. Methods: The task of the challenge was to reconstruct radially acquired multi-coil k-space data (brain/heart) following the method in the original paper, reproducing its key figures. Results were compared to consolidated reference implementations created after the challenge, accounting for the two most common programming languages used in the submissions (Matlab/Python). Results: Visually, differences between submissions were small. Pixel-wise differences originated from image orientation, assumed field-of-view or resolution. The reference implementations were in good agreement, both visually and in terms of image similarity metrics. Discussion and Conclusion: While the description level of the published algorithm enabled participants to reproduce CG-SENSE in general, details of the implementation varied, e.g., density compensation or Tikhonov regularization. Implicit assumptions about the data lead to further differences, emphasizing the importance of sufficient meta-data accompanying open data sets. Defining reproducibility quantitatively turned out to be non-trivial for this image reconstruction challenge, in the absence of ground-truth results. Typical similarity measures like NMSE of SSIM were misled by image intensity scaling and outlier pixels. Thus, to facilitate reproducibility, researchers are encouraged to publish code and data alongside the original paper. Future methodological papers on MR image reconstruction might benefit from the consolidated reference implementations of CG-SENSE presented here, as a benchmark for methods comparison.


title: Nuclear Burning in Collapsar Accretion Disks

abstract: The core collapse of massive, rapidly-rotating stars are thought to be the progenitors of long-duration gamma-ray bursts (GRB) and their associated hyper-energetic supernovae (SNe). At early times after the collapse, relatively low angular momentum material from the infalling stellar envelope will circularize into an accretion disk located just outside the black hole horizon, resulting in high accretion rates necessary to power a GRB jet. Temperatures in the disk midplane at these small radii are sufficiently high to dissociate nuclei, while outflows from the disk can be neutron-rich and may synthesize r-process nuclei. However, at later times, and for high progenitor angular momentum, the outer layers of the stellar envelope can circularize at larger radii $\gtrsim 10^{7}$ cm, where nuclear reaction can take place in the disk midplane ((e.g.~$^{4}$He + $^{16}$O $\rightarrow$ $^{20}$Ne + $\gamma$).. Here we explore the effects of nuclear burning on collapsar accretion disks and their outflows by means of hydrodynamical $\alpha$-viscosity torus simulations coupled to a 19-isotope nuclear reaction network, which are designed to mimic the late infall epochs in collapsar evolution when the viscous time of the torus has become comparable to the envelope fall-back time. Our results address several key questions, such as the conditions for quiescent burning and accretion versus detonation and the generation of $^{56}$Ni in disk outflows, which we show could contribute significantly to powering GRB supernovae. Being located in the slowest, innermost layers of the ejecta, the latter could provide the radioactive heating source necessary to make the spectral signatures of r-process elements visible in late-time GRB-SNe spectra.


title: Neutrino masses, vacuum stability and quantum gravity prediction for the mass of the top quark

abstract: A general prediction from asymptotically safe quantum gravity is the approximate vanishing of all quartic scalar couplings at the UV fixed point beyond the Planck scale. A vanishing Higgs doublet quartic coupling near the Planck scale translates into a prediction for the ratio between the mass of the Higgs boson $M_H$ and the top quark $M_t$. If only the standard model particles contribute to the running of couplings below the Planck mass, the observed $M_H\sim125\,{\rm GeV}$ results in the prediction for the top quark mass $M_t\sim 171\,{\rm GeV}$, in agreement with recent measurements. In this work, we study how the asymptotic safety prediction for the top quark mass is affected by possible physics at an intermediate scale. We investigate the effect of a $SU(2)$ triplet scalar and right-handed neutrinos, needed to explain the tiny mass of left-handed neutrinos. For pure seesaw II, with no or very heavy right handed neutrinos, the top mass can increase to $M_t\sim 172.5\,{\rm GeV}$ for a triplet mass of $M_\Delta\sim 10^8{\rm GeV}$. Right handed neutrino masses at an intermediate scale increase the uncertainty of the predictions of $M_t$ due to unknown Yukawa couplings of the right-handed neutrinos and a cubic interaction in the scalar potential. For an appropriate range of Yukawa couplings there is no longer an issue of vacuum stability.


title: Cannibal domination and the matter power spectrum

abstract: Decoupled hidden sectors can easily and generically result in a period of cannibal domination, during which the dominant component of the Universe has an equation of state intermediate between radiation and matter due to self-heating by number-changing interactions. We present for the first time the consequences of a cannibal-dominated era prior to big bang nucleosynthesis for structure formation on small scales. We find that an early cannibal-dominated era imprints a characteristic peak on the dark matter power spectrum, with scale and amplitude directly determined by the mass, lifetime, and number-changing interaction strength of the cannibal field. This enhancement to the small-scale matter power spectrum will generate early-forming dark matter microhalos, and we provide a detailed and transparent map between the properties of the cannibal species and the characteristic mass and formation time of these structures. These relations demonstrate how the internal workings of a hidden sector leave a potentially observable imprint on the matter power spectrum even if dark matter has no direct couplings to the Standard Model.


title: HCN production in Titan's Atmosphere: Coupling quantum chemistry and disequilibrium atmospheric modeling

abstract: Hydrogen cyanide (HCN) is a critical reactive source of nitrogen for building key biomolecules relevant for the origin of life. Still, many HCN reactions remain uncharacterized by experiments and theory, and the complete picture of HCN production in planetary atmospheres is not fully understood. To improve this situation, we develop a novel technique making use of computational quantum chemistry, experimental data, and atmospheric numerical simulations. First, we use quantum chemistry simulations to explore the entire field of possible reactions for a list of primary species in N2-, CH4-, and H2-dominated atmospheres. In this process, we discover 33 new reactions with no previously known rate coefficients. From here, we develop a consistent reduced atmospheric hybrid chemical network (CRAHCN) containing experimental values when available, and our calculated rate coefficients otherwise. Next, we couple CRAHCN to a 1D chemical kinetic model (ChemKM) to compute the HCN abundance as a function of atmospheric depth on Titan. Our simulated atmospheric HCN profile agrees very well with the Cassini observations. CRAHCN contains 104 reactions however nearly all of the simulated atmospheric HCN profile can be obtained using a scaled down network of only 19 dominant reactions. From here, we form a complete picture of HCN chemistry in Titan's atmosphere, from the dissociation of the main atmospheric species, down to the direct production of HCN along 4 major channels. One of these channels was first discovered and characterized in Pearce et al. (2019) and this work.


title: PSF Deconvolution of the IFU Data and Restoration of Galaxy Stellar Kinematics

abstract: We present a performance test of the Point Spread Function deconvolution algorithm applied to astronomical Integral Field Unit (IFU) Spectroscopy data for restoration of galaxy kinematics. We deconvolve the IFU data by applying the Lucy-Richardson algorithm to the 2D image slice at each wavelength. We demonstrate that the algorithm can effectively recover the true stellar kinematics of the galaxy, by using mock IFU data with diverse combination of surface brightness profile, S/N, line-of-sight geometry and Line-Of-Sight Velocity Distribution (LOSVD). In addition, we show that the proxy of the spin parameter $\lambda_{R_{e}}$ can be accurately measured from the deconvolved IFU data. We apply the deconvolution algorithm to the actual SDSS-IV MaNGA IFU survey data. The 2D LOSVD, geometry and $\lambda_{R_{e}}$ measured from the deconvolved MaNGA IFU data exhibit noticeable difference compared to the ones measured from the original IFU data. The method can be applied to any other regular-grid IFU data to extract the PSF-deconvolved spatial information.


title: Discovery of a New WZ Sagittae Type Cataclysmic Variable in the Kepler/K2 Data

abstract: We identify a new, bright transient in the Kepler/K2 Campaign 11 field. Its light curve rises over seven magnitudes in a day and then declines three magnitudes over a month before quickly fading another two magnitudes. The transient was still detectable at the end of the campaign. The light curve is consistent with a WZ~Sge type dwarf nova outburst. Early superhumps with a period of 82 minutes are seen in the first 10 days and suggest that this is the orbital period of the binary which is typical for the WZ~Sge class. Strong superhump oscillations develop ten days after peak brightness with periods ranging between 83 and 84 minutes. At 25 days after the peak brightness a bump in the light curve appears to signal a subtle rebrightening phase implying that this was an unusual type-A outburst. This is the only WZ~Sge type system observed by Kepler/K2 during an outburst. The early rise of this outburst is well-fit with a broken power law. In first 10 hours the system brightened linearly and then transitioned to a steep rise with a power law index of 4.8. Looking at archival Kepler/K2 data and new TESS observations, a linear rise in the first several hours at the initiation of a superoutburst appears to be common in SU~UMa stars.


title: High-redshift radio galaxies: a potential new source of 21-cm fluctuations

abstract: Radio sources are expected to have formed at high redshifts, producing an excess radiation background above the cosmic microwave background (CMB) at low frequencies. Their effect on the redshifted 21-cm signal of neutral hydrogen is usually neglected, as it is assumed that the associated background is small. Recently, an excess radio background above the level of the CMB has been proposed as one of the possible explanations for the unusually strong 21-cm signal from redshift $z\sim 17$ reported by the EDGES collaboration. As a result, the implications of a smooth and extremely strong excess radio background on both the sky-averaged (global) 21-cm signal and its fluctuations have been considered. Here we take into account the inhomogeneity of the radio background created by a population of high-redshift galaxies, and show that it adds a new type of 21-cm fluctuations to the well-known contributions of density, velocity, Ly-${\alpha}$ coupling, heating and reionization. We find that a population of high-redshift galaxies even with a moderately-enhanced radio efficiency (unrelated to the EDGES result) can have a significant effect on the 21-cm power spectrum and global signal in models with weak X-ray heating. For models that can explain the EDGES data, we conduct a large parameter survey to explore their signatures. We show that in such models the 21-cm power spectrum at $z\sim 17$ is enhanced by up to two orders of magnitude compared to the CMB-only standard case, and the shape and time evolution of the power spectrum is significantly modified by the radio fluctuations. These fluctuations are within reach of upcoming radio interferometers. We also find that these models can be significantly constrained by current and future observations of radio sources.


title: 'Deep' Dive into $b \to c$ Anomalies: Standardized and Future-proof Model Selection Using Self-normalizing Neural Networks

abstract: Noting the erroneous proclivity of information-theoretic approaches, like the Akaike information criterion (AIC), to select simpler models while performing model selection with a small sample size, we address the problem of new physics model selection in $b\to c \tau \nu_{\tau}$ decays in this paper by employing a specific machine learning algorithm (self-normalizing neural networks, a.k.a. SNN) for supervised classification and regression, in a model-independent framework. While the outcomes of the classification with real data-set are compared with AIC, with the SNNs outperforming AIC$_c$ in all aspects of model selection, the regression-outcomes are compared with the results from Bayesian analyses; the obtained parameter spaces differ considerably while keeping maximum posterior (MAP) estimates similar. A few of the two-operator scenarios with a tensor-type interaction are found to be the most probable solution for the data. We also test the effectiveness of our trained networks with the expected, more precise data in Belle-II. The trained networks and associated functionalities are supplied for the use of the community.


title: Evidence for a jet and outflow from Sgr A*: a continuum and spectral line study

abstract: We study the environment of Sgr A* using spectral and continuum observations with the ALMA and VLA. Our analysis of sub-arcsecond H30alpha, H39alpha, H52alpha and H56alpha line emission towards Sgr A* confirm the recently published broad peak ~500 km/s~spectrum toward Sgr~A*. We also detect emission at more extreme radial velocities peaking near -2500 and 4000 km/s, within 0.2''. We then present broad band radio continuum images at multiple frequencies on scales from arcseconds to arcminutes. A number of elongated continuum structures lie parallel to the Galactic plane, extending from ~0.4'' to 10'. We note a nonthermal elongated structure on an arcminute scale emanating from Sgr A* at low frequencies between 1 and 1.4 GHz where thermal emission from the mini-spiral is depressed by optical depth effects. The position angle of this elongated structure and the sense of motion of ionized features with respect to Sgr A* suggest a symmetric, collimated jet emerging from Sgr A* with an opening angle of ~30deg and a position angle of ~60deg punching through the medium before accelerating a significant fraction of the orbiting ionized gas to high velocities. The jet with estimated mass flow rate ~1.4x10^{-5} solar mass/yr emerges perpendicular to the equatorial plane of the accretion flow near the event horizon of Sgr A* and runs along the Galactic plane. To explain a number of east-west features near Sgr A*, we also consider the possibility of an outflow component with a wider-angle launched from the accretion flow at larger radii.


title: Measuring the polarization of boosted, hadronic $W$ bosons with jet substructure observables

abstract: In this work, we present a new technique to measure the longitudinal and transverse polarization fractions of hadronic decays of boosted $W$ bosons. We introduce a new jet substructure observable denoted as $p_\theta$, which is a proxy for the parton level decay polar angle of the $W$ boson in its rest-frame. We show that the distribution of this observable is sensitive to the polarization of $W$ bosons and can therefore be used to reconstruct the $W$ polarization in a model-independent way. As a test case, we study the efficacy of our technique on vector boson scattering processes at the high luminosity Large Hadron Collider and we find that our technique can determine the longitudinal polarization fraction to within $\pm 0.15$. We also show that our technique can be used to identify the parity of beyond Standard Model scalar or pseudo-scalar resonances decaying to $W$ bosons with just 20 events.


title: Sun-as-a-star Spectral Irradiance Observations of Transiting Active Regions

abstract: Major solar flares are prone to occur in active region atmospheres associated with large, complex, dynamically-evolving sunspots. This points to the importance of monitoring the evolution of starspots, not only in visible but also in ultra violet (UV) and X-rays, in understanding the origin and occurrence of stellar flares. To this end, we perform spectral irradiance analysis on different types of transiting solar active regions by using a variety of full-disk synoptic observations. The target events are an isolated sunspot, spotless plage, and emerging flux in prolonged quiet-Sun conditions selected from the past decade. We find that the visible continuum and total solar irradiance become darkened when the spot is at the central meridian, whereas it is bright near the solar limb; UV bands sensitive to the chromosphere correlate well with the variation of total unsigned magnetic flux in the photosphere; amplitudes of EUV and soft X-ray increase with the characteristic temperature, whose light curves are flat-topped due to their sensitivity to the optically thin corona; the transiting spotless plage does not show the darkening in the visible irradiance, while the emerging flux produces an asymmetry in all light curves about the central meridian. The multi-wavelength sun-as-a-star study described here indicates that such time lags between the coronal and photospheric light curves have the potential to probe the extent of coronal magnetic fields above the starspots. In addition, EUV wavelengths that are sensitive to the transition-region temperature sometimes show anti-phased variations, which may be used for diagnosing plasmas around starspots.


title: Structure and rotation of young massive star clusters in a simulated dwarf starburst

abstract: We analyze the three-dimensional shapes and kinematics of the young star cluster population forming in a high-resolution GRIFFIN project simulation of a metal-poor dwarf galaxy starburst. The star clusters, which follow a power-law mass distribution, form from the cold ISM phase with an IMF sampled with individual stars down to 4 solar masses at sub-parsec spatial resolution. Massive stars and their important feedback mechanisms are modelled in detail. The simulated clusters follow a surprisingly tight relation between the specific angular momentum and mass with indications of two sub-populations. Massive clusters ($M_\mathrm{cl}\gtrsim 3\times 10^4 M_{\odot})$ have the highest specific angular momenta at low ellipticities ($\epsilon\sim 0.2$) and show alignment between their shapes and rotation. Lower mass clusters have lower specific angular momenta with larger scatter, show a broader range of elongations, and are typically misaligned indicating that they are not shaped by rotation. The most massive clusters $(M \gtrsim 10^5\,M_{\odot})$ accrete gas and proto-clusters from a $ \lesssim 100\,\rm pc$ scale local galactic environment on a $t \lesssim 10\,\rm Myr$ timescale, inheriting the ambient angular momentum properties. Their two-dimensional kinematic maps show ordered rotation at formation, up to $v \sim 8.5\,\rm km s^{-1}$, consistent with observed young massive clusters and old globular clusters, which they might evolve into. The massive clusters have angular momentum parameters $\lambda_R\lesssim 0.5$ and show Gauss-Hermite coefficients $h_3$ that are anti-correlated with the velocity, indicating asymmetric line-of-sight velocity distributions as a signature of a dissipative formation process.


title: The Exotic Type Ic Broad-Lined Supernova SN 2018gep: Blurring the Line Between Supernovae and Fast Optical Transients

abstract: In the last decade a number of rapidly evolving transients have been discovered that are not easily explained by traditional supernovae models. We present optical and UV data on onee such object, SN 2018gep, that displayed a fast rise with a mostly featureless blue continuum around maximum light, and evolved to develop broad features more typical of a SN Ic-bl while retaining significant amounts of blue flux throughout its observations. The blue excess is most evident in its near-UV flux that is over 4 magnitudes brighter than other stripped envelope supernovae, but also visible in optical g$-$r colors at early times. Its fast rise time of $t_{\rm rise,V} \lesssim 6.2 \pm 0.8$ days puts it squarely in the emerging class of Fast Evolving Luminous Transients, or Fast Blue Optical Transients. With a peak absolute magnitude of M$_r=-19.49 \pm 0.23 $ mag it is on the extreme end of both the rise time and peak magnitude distribution for SNe Ic-bl. Only one other SN Ic-bl has similar properties, iPTF16asu, for which less of the important early time and UV data have been obtained. We show that the objects SNe 2018gep and iPTF16asu have similar photometric and spectroscopic properties and that they overall share many similarities with both SNe Ic-bl and Fast Evolving Transients. We obtain IFU observations of the SN 2018gep host galaxy and derive a number of properties for it. We show that the derived host galaxy properties for both SN 2018gep and iPTF16asu are overall consistent with the SNe Ic-bl and GRB/SNe sample while being on the extreme edge of the observed Fast Evolving Transient sample. These photometric observations are consistent with a simple SN Ic-bl model that has an additional form of energy injection at early times that drives the observed rapid, blue rise, and we speculate that this additional power source may extrapolate to the broader Fast Evolving Transient sample.


title: Prospects for distinguishing galaxy evolution models with surveys at redshifts $z \gtrsim 4$

abstract: Many semi-empirical galaxy formation models have recently emerged to interpret high-$z$ galaxy luminosity functions and make predictions for future galaxy surveys. A common approach assumes a "universal" star formation efficiency, $f_{\ast}$, independent of cosmic time but strongly dependent on the masses of dark matter halos. Though this class of models has been very successful in matching observations over much of cosmic history, simple stellar feedback models do predict redshift evolution in $f_{\ast}$, and are commonly used in semi-analytic models. In this work, we calibrate a set of universal $f_{\ast}$ and feedback-regulated models to the same set of rest-ultraviolet $z \gtrsim 4$ observations, and find that a rapid, $\sim (1+z)^{-3/2}$ decline in both the efficiency of dust production and duty cycle of star formation are needed to reconcile feedback-regulated models with current observations. By construction, these models remain nearly identical to universal $f_{\ast}$ models in rest-ultraviolet luminosity functions and colours. As a result, the only way to distinguish these competing scenarios is either via (i) improved constraints on the clustering of galaxies -- universal and feedback-regulated models differ in predictions for the galaxy bias by $0.1 \lesssim \Delta \langle b \rangle \lesssim 0.3$ over $4 \lesssim z \lesssim 10$ -- or (ii) independent constraints on the dust contents and/or duty cycle of star formation. This suggests that improved constraints on the `dustiness' and `burstiness' of high-$z$ galaxies will not merely add clarity to a given model of star formation in high-$z$ galaxies, but rather fundamentally determine our ability to identify the correct model in the first place.


title: IceCube-Gen2: The Window to the Extreme Universe

abstract: The observation of electromagnetic radiation from radio to $\gamma$-ray wavelengths has provided a wealth of information about the universe. However, at PeV (10$^{15}$ eV) energies and above, most of the universe is impenetrable to photons. New messengers, namely cosmic neutrinos, are needed to explore the most extreme environments of the universe where black holes, neutron stars, and stellar explosions transform gravitational energy into non-thermal cosmic rays. The discovery of cosmic neutrinos with IceCube has opened this new window on the universe. In this white paper, we present an overview of a next-generation instrument, IceCube-Gen2, which will sharpen our understanding of the processes and environments that govern the universe at the highest energies. IceCube-Gen2 is designed to: 1) Resolve the high-energy neutrino sky from TeV to EeV energies; 2) Investigate cosmic particle acceleration through multi-messenger observations; 3) Reveal the sources and propagation of the highest energy particles in the universe; 4) Probe fundamental physics with high-energy neutrinos. IceCube-Gen2 will increase the annual rate of observed cosmic neutrinos by a factor of ten compared to IceCube, and will be able to detect sources five times fainter than its predecessor. Furthermore, through the addition of a radio array, IceCube-Gen2 will extend the energy range by several orders of magnitude compared to IceCube. Construction will take 8 years and cost about \$350M. The goal is to have IceCube-Gen2 fully operational by 2033. IceCube-Gen2 will play an essential role in shaping the new era of multi-messenger astronomy, fundamentally advancing our knowledge of the high-energy universe. This challenging mission can be fully addressed only in concert with the new survey instruments across the electromagnetic spectrum and gravitational wave detectors which will be available in the coming years.


title: Anisotropy of Long-period Comets Explained by Their Formation Process

abstract: Long-period comets coming from the Oort cloud are thought to be planetesimals formed in the planetary region on the ecliptic plane. We have investigated the orbital evolution of these bodies due to the Galactic tide. We extended Higuchi et al. (2007) and derived the analytical solutions to the Galactic longitude and latitude of the direction of aphelion, $L$ and $B$. Using the analytical solutions, we show that the ratio of the periods of the evolution of $L$ and $B$ is very close to either 2 or $\infty$ for initial eccentricities $e_i\simeq 1$, as is true for the Oort cloud comets. From the relation between $L$ and $B$, we predict that Oort cloud comets returning to the planetary region are concentrated on the ecliptic plane and a second plane, which we call the "empty ecliptic." This consists in a rotation of the ecliptic around the Galactic pole by 180$^\circ$. Our numerical integrations confirm that the radial component of the Galactic tide, which is neglected in the derivation of the analytical solutions, is not strong enough to break the relation between $L$ and $B$ derived analytically. Brief examination of observational data shows that there are concentrations near both the ecliptic and the empty ecliptic. We also show that the anomalies of the distribution of $B$ of long-period comets mentioned by several authors are explained by the concentrations on the two planes more consistently than by previous explanations.


title: The Mysterious Bursts observed by Telescope Array and Axion Quark Nuggets

abstract: Telescope Array (TA) experiment has recorded \cite{Abbasi:2017rvx,Okuda_2019} several short time bursts of air shower like events. These bursts are very distinct from conventional single showers, and are found to be strongly correlated with lightnings. We propose that these bursts represent the direct manifestation of the dark matter annihilation events within the so-called axion quark nugget (AQN) model, which was originally invented for completely different purpose to explain the observed similarity between the dark and the visible components in the Universe, i.e. $\Omega_{\rm DM}\sim \Omega_{\rm visible}$ without any fitting parameters. We support this proposal by demonstrating that the observations \cite{Abbasi:2017rvx,Okuda_2019}, including the frequency of appearance, temporal and spatial distributions, intensity, and other related observables are nicely match the emission features of the AQNs propagating in the atmosphere under thunderstorm. We propose to test these ideas by reanalyzing the existing data by increasing the cutoff time scale $\Delta t=1$ ms for the bursts. We also suggest to test this proposal by analyzing the correlations with proper infrasound and seismic instruments. We also suggest to search for the radio signal with frequency $\nu\in (3-300)$ MHz which must be synchronized with bursts.

categories: hep-ph astro-ph.EP astro-ph.HE physics.ao-ph


title: Coulomb interactions and renormalization of semi-Dirac fermions near a topological Lifshitz transition

abstract: We aim to understand how the spectrum of semi-Dirac fermions is renormalized due to long-range Coulomb electron-electron interactions at a topological Lifshitz transition, where two Dirac cones merge. At the transition, the electronic spectrum is characterized by massive quadratic dispersion in one direction, while it remains linear in the other. We have found that, to lowest order, the unconventional log squared (double logarithmic) correction to the quasiparticle mass in bare perturbation theory leads to resummation into strong mass renormalization in the exact full solution of the perturbative renormalization group equations. This behavior effectively wipes out the curvature of the dispersion and leads to Dirac cone restoration at low energy: the system flows towards Dirac dispersion which is anisotropic but linear in momentum, with interaction-depended logarithmic modulation. The Berry phase associated with the restored critical Dirac spectrum is zero - a property guaranteed by time-reversal symmetry and unchanged by renormalization. Our results are in contrast with the behavior that has been found within the large-$N$ approach.


title: The MOSDEF Survey: Differences in SFR and Metallicity for Morphologically-Selected Mergers at z~2

abstract: We study the properties of 55 morphologically-identified merging galaxy systems at z~2. These systems are flagged as mergers based on features such as tidal tails, double nuclei, and asymmetry. Our sample is drawn from the MOSFIRE Deep Evolution Field (MOSDEF) survey, along with a control sample of isolated galaxies at the same redshift. We consider the relationships between stellar mass, star formation rate (SFR), and gas-phase metallicity for both merging and non-merging systems. In the local universe, merging systems are characterized by an elevated SFR and depressed metallicity compared to isolated systems at a given mass. Our results indicate SFR enhancement and metallicity deficit for merging systems relative to non-merging systems for a fixed stellar mass at z~2, though larger samples are required to establish these preliminary results with higher statistical significance. In future work, it will be important to establish if the enhanced SFR and depressed metallicity in high-redshift mergers deviate from the "fundamental metallicity relation," as is observed in mergers in the local universe, and therefore shed light on gas flows during galaxy interactions.


title: Quantum Entangled-Probe Scattering Theory

abstract: We develop an entangled-probe scattering theory, including quantum detection, that extends the scope of standard scattering approaches. We argue that these probes may be revolutionary in studying entangled matter such as unconventional phases of strongly correlated systems. Our presentation focuses on a neutron beam probe that is mode-entangled in spin and path as is experimentally realized in [1], although similar ideas also apply to photon probes. We generalize the traditional van Hove theory [2] whereby the response is written as a properly-crafted combination of two-point correlation functions. Tuning the probe's entanglement length allows us to interrogate spatial scales of interest by analyzing interference patterns in the differential cross-section. Remarkably, for a spin dimer target we find that the typical Young-like interference pattern observed if the target state is un-entangled gets quantum erased when that state becomes maximally entangled.

categories: quant-ph cond-mat.mes-hall cond-mat.str-el nucl-th


title: The MOSDEF Survey: Calibrating the relationship between H$\alpha$ star-formation rate and radio continuum luminosity at $1.4 < z < 2.6$

abstract: The observed empirical relation between the star-formation rates (SFR) of low-redshift galaxies and their radio continuum luminosity offers a potential means of measuring SFR in high redshift galaxies that is unaffected by dust obscuration. In this study, we make the first test for redshift evolution in the SFR-radio continuum relation at high redshift using dust-corrected H$\alpha$ SFR. Our sample consists of 178 galaxies from the MOSFIRE Deep Evolution Field (MOSDEF) Survey at $1.4 < z < 2.6$ with rest-frame optical spectroscopy and deep 1.5 GHz radio continuum observations from the Karl G. Jansky Very Large Array (VLA) GOODS North field. Using a stacking analysis we compare the observed radio continuum luminosities with those predicted from the dust-corrected H$\alpha$ SFR assuming a range of $z\sim0$ relations. We find no evidence for a systematic evolution with redshift, when stacking the radio continuum as a function of dust-corrected H$\alpha$ SFR and when stacking both optical spectroscopy and radio continuum as a function of stellar mass. We conclude that locally calibrated relations between SFR and radio continuum luminosity remain valid out to $z\sim 2$.


title: MHV Graviton Scattering Amplitudes and Current Algebra on the Celestial Sphere

abstract: The Cachazo-Strominger subleading soft graviton theorem for a positive helicity soft graviton is equivalent to the Ward identities for $\overline{SL(2,\mathbb C)}$ currents. This naturally gives rise to a $\overline{SL(2,\mathbb C)}$ current algebra living on the celestial sphere. The generators of the $\overline{SL(2,\mathbb C)}$ current algebra and the supertranslations, coming from a positive helicity leading soft graviton, form a closed algebra. We find that the OPE of two graviton primaries in the Celestial CFT, extracted from MHV amplitudes, is completely determined in terms of this algebra. To be more precise, 1) The subleading terms in the OPE are determined in terms of the leading OPE coefficient if we demand that both sides of the OPE transform in the same way under this local symmetry algebra. 2) Positive helicity gravitons have null states under this local algebra whose decoupling leads to differential equations for MHV amplitudes. An $n$ point MHV amplitude satisfies two systems of $(n-2)$ linear first order PDEs corresponding to $(n-2)$ positive helicity gravitons. We have checked, using Hodges' formula, that one system of differential equations is satisfied by any MHV amplitude, whereas the other system has been checked up to six graviton MHV amplitude. 3) One can determine the leading OPE coefficients from these differential equations. This points to the existence of an autonomous sector of the Celestial CFT which holographically computes the MHV graviton scattering amplitudes and is completely defined by this local symmetry algebra. The MHV-sector of the Celestial CFT is like a minimal model of $2$-D CFT.


title: Nonperturbative Topological Current in Weyl and Dirac Semimetals in Laser Fields

abstract: We study non-perturbatively the anomalous Hall current and its high harmonics generated in Weyl and Dirac semimetals by strong elliptically polarized laser fields, in the context of kinetic theory. We find a novel crossover between perturbative and non-perturbative regimes characterized by the electric field strength $\mathcal{E}^{*}= \frac{\mu \omega}{ 2 e v_F}$ ($\omega$: laser frequency, $\mu$: Fermi energy, $v_F$: Fermi velocity). In the perturbative regime, the anomalous Hall current quadratically depends on the field strength ($\mathcal{E}$), whereas the higher order corrections, as well as high harmonics, vanish at zero temperature. In the non-perturbative regime, the anomalous Hall current saturates and decays as $(\log{\mathcal{E}})/\mathcal{E}$, while even-order high harmonics are generated when inplane rotational symmetry is broken. Based on the analytical solution of the Boltzmann equation, we reveal the topological origin of the sharp crossover: the Weyl monopole stays inside or moves outside of the Fermi sphere, respectively, during its fictitious motion in the pertubative or non-pertubative regimes. Our findings establish a new non-linear response intrinsically connected to topology, characteristic to Weyl and Dirac semimetals.


title: The Gilmore-NASG model to predict single-bubble cavitation in compressible liquids

abstract: The Gilmore model is combined with the Noble-Abel-stiffened-gas (NASG) equation of state to yield a simple model to predict the expansion and collapse of spherical bubbles based on real gas thermodynamics. The NASG equation of state resolves the temperature inaccuracy associated with the commonly employed Tait equation of state for liquids and, thus, can provide a consistent description of compressible and thermal effects of the bubble content and the surrounding liquid during cavitation. After a detailed derivation of the proposed Gilmore-NASG model, the differences between the classical Gilmore-Tait model and the proposed model are highlighted with results of single-bubble cavitation related to bubble collapse and driven by an acoustic excitation in frequency and amplitude regimes relevant to sonoluminescence, high-intensity focused ultrasound and shock wave lithotripsy. Especially for rapidly and violently collapsing bubbles, substantial differences in the bubble behaviour can be observed between the proposed Gilmore-NASG model and the classical Gilmore-Tait model. The ability of the Gilmore-NASG model to simultaneously predict reliable pressure and temperature values in gas, vapour and liquid, makes the proposed model particularly attractive for sonochemistry and biomedical applications.