title: PDE-constrained optimal control problems with uncertain parameters using SAGA

abstract: We consider an optimal control problem (OCP) for a partial differential equation (PDE) with random coefficients. The optimal control function is a deterministic, distributed forcing term that minimizes an expected quadratic regularized loss functional. For the numerical approximation of this PDE-constrained OCP, we replace the expectation in the objective functional by a suitable quadrature formula and, eventually, discretize the PDE by a Galerkin method. To practically solve such approximate OCP, we propose an importance sampling version the SAGA algorithm, a type of Stochastic Gradient algorithm with a fixed-length memory term, which computes at each iteration the gradient of the loss functional in only one quadrature point, randomly chosen from a possibly non-uniform distribution. We provide a full error and complexity analysis of the proposed numerical scheme. In particular we compare the complexity of the generalized SAGA algorithm with importance sampling, with that of the Stochastic Gradient (SG) and the Conjugate Gradient (CG) algorithms, applied to the same discretized OCP.We show that SAGA converges exponentially in the number of iterations as for a CG algorithm and has a similar asymptotic computational complexity, in terms of computational cost versus accuracy (proportional with the time required if no parallel computing is used). Moreover, it features good pre-asymptotic properties, as shown by our numerical experiments, which makes it appealing in a limited budget context.


title: Revisiting the scaling of citations for research assessment

abstract: Over the past decade, national research evaluation exercises, traditionally conducted using the peer review method, have begun opening to bibliometric indicators. The citations received by a publication are assumed as proxy for its quality, but they require standardization prior to use in comparative evaluation of organizations or individual scientists: the citation data must be standardized, due to the varying citation behavior across research fields. The objective of this paper is to compare the effectiveness of the different methods of normalizing citations, in order to provide useful indications to research assessment practitioners. Simulating a typical national research assessment exercise, he analysis is conducted for all subject categories in the hard sciences and is based on the Thomson Reuters Science Citation Index-Expanded. Comparisons show that the citations average is the most effective scaling parameter, when the average is based only on the publications actually cited.


title: A consistent view of the samarium hexaboride terminations to resolve the nature of its surface states

abstract: The research effort prompted by the prediction that SmB$_6$ could be the first topological Kondo insulator has produced a wealth of new results, though not all of these seem compatible. A major discrepancy exists between scanning tunneling microscopy / spectroscopy (STM/S) and angle-resolved photoemission spectroscopy (ARPES), because the two experimental methods suggest a very different number of terminations of the (100) surface with different properties. Here we tackle this issue in a combined STM/S and ARPES study. We find that two of the well-ordered topographies reported in earlier STM studies can be associated with the crystal terminations identified using photoemission. We further observe a reversal of the STM contrast with bias voltage for one of the topographies. We ascribe this result to a different energy dependence of Sm and B-derived states, and show that it can be used to obtain element specific images of SmB$_6$ and identify which topography belongs to which termination. We finally find STS results to support a modification of the low-energy electronic structure at the surface that has been proposed as the trivial origin of surface metallicity in this material.


title: Maintaining Grasps within Slipping Bound by Monitoring Incipient Slip

abstract: In this paper, we propose an approach to detect incipient slip, i.e. predict slip, by using a high-resolution vision-based tactile sensor, GelSlim. The sensor dynamically captures the tactile imprints of the contact object and their changes with a soft gel pad. The method assumes the object is mostly rigid and treats the motion of object's imprint on sensor surface as a 2D rigid-body motion. We use the deviation of the true motion field from that of a 2D planar rigid transformation as a measure of slip. The output is a dense slip field which we use to detect when small areas of the contact patch start to slip (incipient slip). The method can detect both translational and rotational incipient slip without any prior knowledge of the object at 24 Hz. We test the method on 10 objects 240 times and achieve 86.25% detection accuracy. We further show how the slip feedback can be used to monitor the gripping force to avoid slip with a closed-loop bottle-cap screwing and unscrewing experiment with incipient slip detection feedback. The method was demonstrated to be useful for the robot to apply proper gripping force and stop screwing at the right point before breaking objects. The method can be applied to many manipulation tasks in both structured and unstructured environments.


title: An individual-level assessment of the relationship between spin-off activities and research performance in universities

abstract: One of the most problematic aspects in the creation of spin-offs by university personnel concerns the relationship between entrepreneurial activity and research activity by researcher-entrepreneurs. The literature has expressed varying and opposing views as to the nature of the relationship but very little has been produced to empirically legitimate one position or another. The present work proposes to address this shortcoming by exploring the relationship existing between academic spin-off generation and the research performance of enterprise founders. The study investigates whether, and to what extent, scientific performance by academic entrepreneurs is different than that of their colleagues, and if the involvement in entrepreneurial activity has an influence on the individual's research activity. The research questions are answered by considering all spin-offs generated by Italian universities over the period 2001-2008 and evaluating, through a bibliometric approach, the scientific performance of founders relative to that of their colleagues who carry out research in the same field. The data show better scientific performance by the researcher-entrepreneurs than that of their colleagues and in addition, although there are some variations across fields, the creation of a spin-off does not seem, on average, to have negative effects on the scientific performance of the founders.


title: Bottomonium spectroscopy motivated by general features of pNRQCD

abstract: The bottomonium mass spectra is computed in the framework of potential non-relativistic quantum chromodynamics. The potential consists of a static term incorporating Coulombic plus confinement part along with a correction term added non-perturbatively from pNRQCD, which is classified in powers of the inverse of heavy quark mass \textit{O}$(1/m)$. The masses of excited bottomonia are calculated by perturbatively adding spin-hyperfine, spin-orbit and tensor components of one gluon exchange interactions in powers of \textit{O}$(1/m^2)$. Calculated masses are found to be consistent with other theoretical studies and experimental data. The Regge trajectories of the calculated mass spectra are also constructed. The values of the wave-functions are extracted and employed to calculate the electromagnetic transition widths and $\gamma\gamma$, $e^+e^-$, light hadron and $\gamma\gamma\gamma$ decay widths of several states at various leading orders, within the non-relativistic QCD formalism. Some of the experimentally reported states of bottomonium family like $\Upsilon$(10860), $\Upsilon$(11020) and X(10610) are identified as mixed $S-D$ wave and $P$-wave states.


title: The role of temperature dependent string-inspired CPT violating backgrounds in leptogenesis and the chiral magnetic effect

abstract: In a temperature dependent CPT-Violating (CPTV) axial time-like background (induced by the Kalb-Ramond tensor field of string theory) we discuss leptogenesis by solving the Boltzmann equation.The current work non-trivially modifies the framework of a previous phenomenological approach by the authors where the CPTV axial background was considered to be a constant (with no microscopic justification). The constant background approximation though is shown to capture the main phenomenological features of leptogenesis. On comparing our analysis to the related chiral magnetic effect for axial current condensates, we conclude that the Kalb-Ramond field does not play the role of the chiral chemical potential needed for that effect.


title: The Statistical Model of Nuclear Reactions: Open Problems

abstract: Several experiments [1-3] show significant deviations from predictions of the statistical model of nuclear reactions. We summarize unsuccessful recent theoretical efforts to account for such disagreement in terms of a violation of orthogonal invariance caused by the Thomas-Ehrman shift. We report on numerical simulations involving a large number of gamma decay channels that also give rise to violation of orthogonal invariance but likewise do not account for the discrepancies. We discuss the statistical model in the light of these results.


title: Twin beam quantum-enhanced correlated interferometry for testing fundamental physics

abstract: Quantum metrology deals with improving the resolution of instruments that are otherwise limited by shot noise and it is therefore a promising avenue for enabling scientific breakthroughs. The advantage can be even more striking when quantum enhancement is combined with correlation techniques among several devices. Here, we present and realize a correlation interferometry scheme exploiting bipartite quantum correlated states injected in two independent interferometers. The scheme outperforms classical analogues in detecting a faint signal that may be correlated/uncorrelated between the two devices. We also compare its sensitivity with that obtained for a pair of two independent squeezed modes, each addressed to one interferometer, for detecting a correlated stochastic signal in the MHz frequency band. Being the simpler solution, it may eventually find application to fundamental physics tests, e.g., searching for the effects predicted by some Planck scale theories.


title: A Stochastic Penalty Model for Convex and Nonconvex Optimization with Big Constraints

abstract: The last decade witnessed a rise in the importance of supervised learning applications involving {\em big data} and {\em big models}. Big data refers to situations where the amounts of training data available and needed causes difficulties in the training phase of the pipeline. Big model refers to situations where large dimensional and over-parameterized models are needed for the application at hand. Both of these phenomena lead to a dramatic increase in research activity aimed at taming the issues via the design of new sophisticated optimization algorithms. In this paper we turn attention to the {\em big constraints} scenario and argue that elaborate machine learning systems of the future will necessarily need to account for a large number of real-world constraints, which will need to be incorporated in the training process. This line of work is largely unexplored, and provides ample opportunities for future work and applications. To handle the {\em big constraints} regime, we propose a {\em stochastic penalty} formulation which {\em reduces the problem to the well understood big data regime}. Our formulation has many interesting properties which relate it to the original problem in various ways, with mathematical guarantees. We give a number of results specialized to nonconvex loss functions, smooth convex functions, strongly convex functions and convex constraints. We show through experiments that our approach can beat competing approaches by several orders of magnitude when a medium accuracy solution is required.


title: NU Ori: a hierarchical triple system with a strongly magnetic B-type star

abstract: NU Ori is a massive spectroscopic and visual binary in the Orion Nebula Cluster, with 4 components: Aa, Ab, B, and C. The B0.5 primary (Aa) is one of the most massive B-type stars reported to host a magnetic field. We report the detection of a spectroscopic contribution from the C component in high-resolution ESPaDOnS spectra, which is also detected in a Very Large Telescope Interferometer (VLTI) dataset. Radial velocity (RV) measurements of the inner binary (designated Aab) yield an orbital period of 14.3027(7) d. The orbit of the third component (designated C) was constrained using both RVs and interferometry. We find C to be on a mildly eccentric 476(1) d orbit. Thanks to spectral disentangling of mean line profiles obtained via least-squares deconvolution we show that the Zeeman Stokes $V$ signature is clearly associated with C, rather than Aa as previously assumed. The physical parameters of the stars were constrained using both orbital and evolutionary models, yielding $M_{\rm Aa} = 14.9 \pm 0.5 M_\odot$, $M_{\rm Ab} = 3.9 \pm 0.7 M_\odot$, and $M_{\rm C} = 7.8 \pm 0.7 M_\odot$. The rotational period obtained from longitudinal magnetic field $\langle B_z \rangle$ measurements is $P_{\rm rot} = 1.09468(7)$ d, consistent with previous results. Modeling of $\langle B_z \rangle$ indicates a surface dipole magnetic field strength of $\sim 8$ kG. NU Ori C has a magnetic field strength, rotational velocity, and luminosity similar to many other stars exhibiting magnetospheric H$\alpha$ emission, and we find marginal evidence of emission at the expected level ($\sim$1% of the continuum).


title: A bibliometric tool to assess the regional dimension of university-industry research collaborations

abstract: The present work proposes a bibliometric methodology for measuring the grade of correspondence between regional industry's demand for research collaboration and supply from public laboratories. The methodology also permits measurement of the intensity and direction of the regional flows of knowledge in public-private collaborations. The aim is to provide a diagnostic instrument for regional and national policy makers, which could add to existing ones to plan interventions for re-balancing sectorial public supply of knowledge with industrial absorptive capacity, and maximizing appropriability of knowledge spillovers. The methodology is applied to university-industry collaborations in the hard sciences in all Italian administrative regions.


title: Unexpected Zero Bias Conductance Peak on the Topological Semimetal Sb(111) with a Broken Bilayer

abstract: The long-sought Majorana fermion is expected to manifest in a topological-superconductor heterostructure as a zero bias conductance peak (ZBCP). As one promising platform for such heterostructures, we investigate the cleaved surface of the topological semimetal Sb(111) using scanning tunneling microscopy and spectroscopy. Remarkably, we find a robust ZBCP on some terraces of the cleaved surface, although no superconductor is present. Using quasiparticle interference imaging, Landau level spectroscopy and density functional theory, we show that the ZBCP originates from a van Hove singularity pushed up to the Fermi level by a sub-surface stacking fault. Amidst the sprint to stake claims on new Majorana fermion systems, our finding highlights the importance of using a local probe together with detailed modeling to check thoroughly for crystal imperfections that may give rise to a trivial ZBCP unrelated to Majorana physics.


title: Picking Apart Story Salads

abstract: During natural disasters and conflicts, information about what happened is often confusing, messy, and distributed across many sources. We would like to be able to automatically identify relevant information and assemble it into coherent narratives of what happened. To make this task accessible to neural models, we introduce Story Salads, mixtures of multiple documents that can be generated at scale. By exploiting the Wikipedia hierarchy, we can generate salads that exhibit challenging inference problems. Story salads give rise to a novel, challenging clustering task, where the objective is to group sentences from the same narratives. We demonstrate that simple bag-of-words similarity clustering falls short on this task and that it is necessary to take into account global context and coherence.


title: Multicritical Fermi surface topological transitions

abstract: A wide variety of complex phases in quantum materials are driven by electron-electron interactions, which are enhanced through density of states peaks. A well known example occurs at van Hove singularities where the Fermi surface undergoes a topological transition. Here we show that higher order singularities, where multiple disconnected leaves of Fermi surface touch all at once, naturally occur at points of high symmetry in the Brillouin zone. Such multicritical singularities can lead to stronger divergences in the density of states than canonical van Hove singularities, and critically boost the formation of complex quantum phases via interactions. As a concrete example of the power of these Fermi surface topological transitions, we demonstrate how they can be used in the analysis of experimental data on Sr$_3$Ru$_2$O$_7$. Understanding the related mechanisms opens up new avenues in material design of complex quantum phases.


title: Application and adaptation of the truncated Newton method for non-convex plasticity problems

abstract: Small-scale plasticity problems are often characterised by different patterning behaviours ranging from macroscopic down to the atomistic scale. In successful models of such complex behaviour, its origin lies within non-convexity of the governing free energy functional. A common approach to solve such non-convex problems numerically is by regularisation through a viscous formulation. This, however, may require the system to be overdamped and hence potentially has a strong impact on the obtained results. To avoid this side-effect, this paper addresses the treatment of the full non-convexity by an appropriate numerical solution algorithm -- the truncated Newton method. The presented method is a double iterative approach which successively generates quadratic approximations of the energy landscape and minimises these by an inner iterative scheme, based on the conjugate gradient method. The inner iterations are terminated when either a sufficient energy decrease is achieved or, to incorporate the treatment of non-convexity, a direction of negative curvature is encountered. If the latter never happens, the method reduces to Newton--Raphson iterations, solved by the conjugate gradient method, with a subsequent line search. However, in the case of a non-convex energy it avoids convergence to a saddle point and adds robustness. The stability of the truncated Newton method is demonstrated for the, highly non-convex, Peierls-Nabarro model, solved within a Finite Element framework. A potential drop in efficiency due to an occasional near singular Hessian is remedied by a trust region like extension, which is physically based on the Peierls-Nabarro disregistry profile. The result is an efficient numerical scheme with a high stability that is independent of any regularisation.


title: The wavefunction reconstruction effects in calculation of DM-induced electronic transition in semiconductor targets

abstract: The physics of the electronic excitation in semiconductors induced by sub-GeV dark matter (DM) have been extensively discussed in literature, under the framework of the standard plane wave (PW) and pseudopotential calculation scheme. In this paper, we investigate the implication of the all-electron (AE) reconstruction on estimation of the DM-induced electronic transition event rates. As a benchmark study, we first calculate the wavefunctions in silicon and germanium bulk crystals based on both the AE and pseudo (PS) schemes within the projector augmented wave (PAW) framework, and then make comparisons between the calculated excitation event rates obtained from these two approaches. It turns out that in process where large momentum transfer is kinetically allowed, the two calculated event rates can differ by a factor of a few. Such discrepancies are found to stem from the high-momentum components neglected in the PS scheme. It is thus implied that the correction from the AE wavefunction in the core region is necessary for an accurate estimate of the DM-induced transition event rate in semiconductors.


title: Accelerating SGD with momentum for over-parameterized learning

abstract: Nesterov SGD is widely used for training modern neural networks and other machine learning models. Yet, its advantages over SGD have not been theoretically clarified. Indeed, as we show in our paper, both theoretically and empirically, Nesterov SGD with any parameter selection does not in general provide acceleration over ordinary SGD. Furthermore, Nesterov SGD may diverge for step sizes that ensure convergence of ordinary SGD. This is in contrast to the classical results in the deterministic scenario, where the same step size ensures accelerated convergence of the Nesterov's method over optimal gradient descent. To address the non-acceleration issue, we introduce a compensation term to Nesterov SGD. The resulting algorithm, which we call MaSS, converges for same step sizes as SGD. We prove that MaSS obtains an accelerated convergence rates over SGD for any mini-batch size in the linear setting. For full batch, the convergence rate of MaSS matches the well-known accelerated rate of the Nesterov's method. We also analyze the practically important question of the dependence of the convergence rate and optimal hyper-parameters on the mini-batch size, demonstrating three distinct regimes: linear scaling, diminishing returns and saturation. Experimental evaluation of MaSS for several standard architectures of deep networks, including ResNet and convolutional networks, shows improved performance over SGD, Nesterov SGD and Adam.


title: Basis-invariant conditions for CP symmetry of order 4

abstract: Three-Higgs-doublet models (3HDM) allow for a novel, physically distinct form of CP invariance: CP symmetry of order 4 (CP4). Due to the large basis change freedom in 3HDM, it is imperative to recognize the presence of a possibly hidden CP4 in a basis-invariant way. In the present work, we solve this problem and establish basis-invariant necessary and sufficient conditions for a 3HDM to possess a CP4 symmetry. We also derive a basis-invariant criterion to decide whether or not a CP4 symmetric 3HDM possesses any additional CP symmetry, as well as a criterion to decide whether or not CP4 is spontaneously broken.


title: The geometric challenge of testing gravity with wide binaries

abstract: Wide binaries provide promising laboratories for testing general relativity (GR) in the low-acceleration regime. Recent observational studies have found that the difference in the proper motions and/or radial velocities of the components of nearby wide binaries appear larger than predicted by Kepler's law's, indicating a potential breakdown of GR at low accelerations. These studies have not accounted for projection effects owing to the different position of the two stars on the celestial sphere. I show that two stars in a wide binary with identical 3D space velocities often have significantly different proper motions and radial velocities purely due to projection effects. I construct a sample of simulated binaries that follow Kepler's laws and have similar phase-space distributions to the observed samples of nearby binaries. Beyond separations of $\sim$0.1 pc, direct comparison of the components' proper motions would suggest strong tensions with GR, even though the simulated binaries follow Kepler's laws by construction. The magnitude of the apparent disagreement is similar to that found observationally, suggesting that the apparent tension between observations and GR may largely be due to projection effects. I discuss prospects for constraining gravity at low accelerations with wide binaries. Robust tests of GR are possible with current data but require measurements of 3D velocities. Further work is also needed to model contamination from unbound moving groups and unrecognized hierarchical triples.


title: Stability of synchronous slowly oscillating periodic solutions for systems of delay differential equations with coupled nonlinearity

abstract: We study stability of so-called synchronous slowly oscillating periodic solutions (SOPSs) for a system of identical delay differential equations (DDEs) with linear decay and nonlinear delayed negative feedback that are coupled through their nonlinear term. Under a row sum condition on the coupling matrix, existence of a unique SOPS for the corresponding scalar DDE implies existence of a unique synchronous SOPS for the coupled DDEs. However, stability of the SOPS for the scalar DDE does not generally imply stability of the synchronous SOPS for the coupled DDEs. We obtain an explicit formula, depending only on the spectrum of the coupling matrix, the strength of the linear decay and the values of the nonlinear negative feedback function near plus/minus infinity, that determines the stability of the synchronous SOPS in the asymptotic regime where the nonlinear term is heavily weighted. We also treat the special cases of so-called weakly coupled systems, near uniformly coupled systems, and doubly nonnegative coupled systems, in the aforementioned asymptotic regime. Our approach is to estimate the characteristic (Floquet) multipliers for the synchronous SOPS. We first reduce the analysis of the multidimensional variational equation to the analysis of a family of scalar variational-type equations, and then establish limits for an associated family of monodromy-type operators. We illustrate our results with examples of systems of DDEs with mean-field coupling and systems of DDEs arranged in a ring.


title: Spin molecular-orbit coupling and magnetic properties of the decorated honeycomb layers of Mo3S7(dmit)3 crystals

abstract: We explore the magnetic properties of isolated a-b planes of trinuclear organometallic crystals, Mo3S7(dmit)3, in which an interplay of strong electronic correlations and spin molecular-orbital coupling (SMOC) occurs. The magnetic properties can be described by a XXZ+120, S=1 Heisenberg model on a honeycomb lattice with single-spin anisotropy, D, which depends strongly on SMOC. Based on ab initio estimates of SMOC in Mo3S7(dmit)3 crystals, we predict that the honeycomb layers of Mo3S7(dmit)3 are Neel ordered. However, in materials with a greater degree of magnetic frustration, Neel order can give way to the large-D phase.


title: Differentiable MPC for End-to-end Planning and Control

abstract: We present foundations for using Model Predictive Control (MPC) as a differentiable policy class for reinforcement learning in continuous state and action spaces. This provides one way of leveraging and combining the advantages of model-free and model-based approaches. Specifically, we differentiate through MPC by using the KKT conditions of the convex approximation at a fixed point of the controller. Using this strategy, we are able to learn the cost and dynamics of a controller via end-to-end learning. Our experiments focus on imitation learning in the pendulum and cartpole domains, where we learn the cost and dynamics terms of an MPC policy class. We show that our MPC policies are significantly more data-efficient than a generic neural network and that our method is superior to traditional system identification in a setting where the expert is unrealizable.


title: FAQBism

abstract: We answer several questions that have been Frequently Asked about QBism. These remarks (many of them lighthearted) should be considered supplements to more systematic treatments by the authors and others.


title: Bounding bias due to selection

abstract: When epidemiologic studies are conducted in a subset of the population, selection bias can threaten the validity of causal inference. This bias can occur whether or not that selected population is the target population, and can occur even in the absence of exposure-outcome confounding. However, it is often difficult to quantify the extent of selection bias, and sensitivity analysis can be challenging to undertake and to understand. In this article we demonstrate that the magnitude of the bias due to selection can be bounded by simple expressions defined by parameters characterizing the relationships between unmeasured factor(s) responsible for the bias and the measured variables. No functional form assumptions are necessary about those unmeasured factors. Using knowledge about the selection mechanism, researchers can account for the possible extent of selection bias by specifying the size of the parameters in the bounds. We also show that the bounds, which differ depending on the target population, result in summary measures that can be used to calculate the minimum magnitude of the parameters required to shift a risk ratio to the null. The summary measure can be used to determine the overall strength of selection that would be necessary to explain away a result. We then show that the bounds and summary measures can be simplified in certain contexts or with certain assumptions. Using examples with varying selection mechanisms, we also demonstrate how researchers can implement these simple sensitivity analyses.


title: Dirac quantisation condition: a comprehensive review

abstract: In most introductory courses on electrodynamics, one is taught the electric charge is quantised but no theoretical explanation related to this law of nature is offered. Such an explanation is postponed to graduate courses on electrodynamics, quantum mechanics and quantum field theory, where the famous Dirac quantisation condition is introduced, which states that a single magnetic monopole in the Universe would explain the electric charge quantisation. Even when this condition assumes the existence of a not-yet-detected magnetic monopole, it provides the most accepted explanation for the observed quantisation of the electric charge. However, the usual derivation of the Dirac quantisation condition involves the subtle concept of an "unobservable" semi-infinite magnetised line, the so-called "Dirac string," which may be difficult to grasp in a first view of the subject. The purpose of this review is to survey the concepts underlying the Dirac quantisation condition, in a way that may be accessible to advanced undergraduate and graduate students. Some of the discussed concepts are gauge invariance, singular potentials, single-valuedness of the wave function, undetectability of the Dirac string and quantisation of the electromagnetic angular momentum. Five quantum-mechanical and three semi-classical derivations of the Dirac quantisation condition are reviewed. In addition, a simple derivation of this condition involving heuristic and formal arguments is presented.


title: Unsupervised Identification of Disease Marker Candidates in Retinal OCT Imaging Data

abstract: The identification and quantification of markers in medical images is critical for diagnosis, prognosis, and disease management. Supervised machine learning enables the detection and exploitation of findings that are known a priori after annotation of training examples by experts. However, supervision does not scale well, due to the amount of necessary training examples, and the limitation of the marker vocabulary to known entities. In this proof-of-concept study, we propose unsupervised identification of anomalies as candidates for markers in retinal Optical Coherence Tomography (OCT) imaging data without a constraint to a priori definitions. We identify and categorize marker candidates occurring frequently in the data, and demonstrate that these markers show predictive value in the task of detecting disease. A careful qualitative analysis of the identified data driven markers reveals how their quantifiable occurrence aligns with our current understanding of disease course, in early- and late age-related macular degeneration (AMD) patients. A multi-scale deep denoising autoencoder is trained on healthy images, and a one-class support vector machine identifies anomalies in new data. Clustering in the anomalies identifies stable categories. Using these markers to classify healthy-, early AMD- and late AMD cases yields an accuracy of 81.40%. In a second binary classification experiment on a publicly available data set (healthy vs. intermediate AMD) the model achieves an area under the ROC curve of 0.944.


title: Demonstration of resonant inelastic X-ray scattering as a probe of exciton-phonon coupling

abstract: Resonant inelastic X-ray scattering (RIXS) is a promising technique for obtaining electron-phonon coupling constants. However, the ability to extract these coupling constants throughout the Brillouin zone for crystalline materials remains limited. To address this need, we developed a Green's function formalism to capture electron-phonon contributions to core-level spectroscopies without explicitly solving the full vibronic problem. Our approach is based on the cumulant expansion of the Green's function combined with many-body theory calculated vibrational coupling constants. The methodology is applied to X-ray photoemission spectroscopy, X-ray absorption spectroscopy (XAS), and RIXS. In the case of the XAS and RIXS, we use a 2-particle exciton Green's function, which accounts implicitly for particle-hole interference effects that have previously proved difficult. To demonstrate the methodology and gain a deeper understanding of the experimental technique, we apply our formalism to small molecules, for which unambiguous experimental data exist. This comparison reveals that the vibronic coupling constant probed by RIXS is in fact related to exciton-phonon coupling rather than electron-phonon coupling, challenging the conventional interpretation of the experiment.


title: A link invariant related to Khovanov homology and knot Floer homology

abstract: In this paper we introduce a chain complex $C_{1 \pm 1}(D)$ where D is a plat braid diagram for a knot K. This complex is inspired by knot Floer homology, but it the construction is purely algebraic. It is constructed as an oriented cube of resolutions with differential d=d_0+d_1. We show that the E_2 page of the associated spectral sequence is isomorphic to the Khovanov homology of K, and that the total homology is a link invariant which we conjecture is isomorphic to \delta-graded knot Floer homology. The complex can be refined to a tangle invariant for braids on 2n strands, where the associated invariant is a bimodule over an algebra A_n. We show that A_n is isomorphic to B'(2n+1, n), the algebra used for the DA-bimodule constructed by Ozsvath and Szabo in their algebraic construction of knot Floer homology.


title: On The Inductive Bias of Words in Acoustics-to-Word Models

abstract: Acoustics-to-word models are end-to-end speech recognizers that use words as targets without relying on pronunciation dictionaries or graphemes. These models are notoriously difficult to train due to the lack of linguistic knowledge. It is also unclear how the amount of training data impacts the optimization and generalization of such models. In this work, we study the optimization and generalization of acoustics-to-word models under different amounts of training data. In addition, we study three types of inductive bias, leveraging a pronunciation dictionary, word boundary annotations, and constraints on word durations. We find that constraining word durations leads to the most improvement. Finally, we analyze the word embedding space learned by the model, and find that the space has a structure dominated by the pronunciation of words. This suggests that the contexts of words, instead of their phonetic structure, should be the future focus of inductive bias in acoustics-to-word models.


title: A Symmetric Two Higgs Doublet Model

abstract: After the discovery of a SM-like scalar state with mass of 125 GeV in 2012, the electroweak sector demands to be studied with all possible theoretical and experimental efforts. In this work we present a symmetric two Higgs doublet model with a discrete interchange symmetry between the two Higgs doublets ($\Phi_1$ $\leftrightarrow$ $\Phi_2$). Apart from the SM-like scalar state (h) with $m_h = 125 $ GeV, the model has several distinguishing features including the pseudoscalar (A), the charged scalars($H^\pm$) and the neutral scalar H, which do not have any direct coupling to fermions. The neutral scalar $H$ can have mass lighter than the 125 GeV SM-like Higgs state $h$. Due to the presence of a residual $Z_2$ symmetry after the Spontaneous Symmetry Breaking (SSB), the neutral scalar $H$ emerges as the Dark Matter candidate in this scenario. As an effect of this possibility, the SM-like scalar $h$ will have an extra invisible decay mode of $h \rightarrow H H$ in this framework. We propose the model and discuss some of the interesting features with a guideline of possible phenomenological searches at the LHC present in this scenario.


title: You May Not Need Attention

abstract: In NMT, how far can we get without attention and without separate encoding and decoding? To answer that question, we introduce a recurrent neural translation model that does not use attention and does not have a separate encoder and decoder. Our eager translation model is low-latency, writing target tokens as soon as it reads the first source token, and uses constant memory during decoding. It performs on par with the standard attention-based model of Bahdanau et al. (2014), and better on long sentences.


title: $U(1)$ Vacuum, Chern-Simons Diffusion and Real-Time Simulations

abstract: Despite its importance for subjects ranging from cosmology to plasma physics, first principle simulations of the dynamics associated with a U(1) chiral anomaly have been started only recently. In this work, we report on the current status of these investigations. We discuss a possible set-up and highlight some results. We present a determination of the Chern-Simons diffusion rate, which shows some discrepancy with the usual effective description which is anomalous magnetohydrodynamic. We also present some exploratory results on the behaviour of an initial chiral chemical potential.


title: A quantum alternating operator ansatz with hard and soft constraints for lattice protein folding

abstract: Gate-based universal quantum computers form a rapidly evolving field of quantum computing hardware technology. In previous work, we presented a quantum algorithm for lattice protein folding on a cubic lattice, tailored for quantum annealers. In this paper, we introduce a novel approach for solving the lattice protein folding problem on universal gate-based quantum computing architectures. Lattice protein models are coarse-grained representations of proteins that have been used extensively over the past thirty years to examine the principles of protein folding and design.These models can be used to explore a vast number of possible protein conformations and to infer structural properties of more complex atomistic protein structures. We formulate the problem as a quantum alternating operator ansatz, a member of the wider class of variational quantum/classical hybrid algorithms. To increase the probability of sampling the ground state, we propose splitting the optimization problem into hard and soft constraints. This enables us to use a previously under-utilised component of the variational algorithm to constrain the search to the subspace of solutions that satisfy the hard constraints.


title: Pulsational pair-instability supernovae in very close binaries

abstract: Pair-instability and pulsational pair-instability supernovae (PPISN) have not been unambiguously observed so far. They are, however, promising candidates for the progenitors of the heaviest binary black hole (BBH) mergers detected. If these BBHs are the product of binary evolution, then PPISNe could occur in very close binaries. Motivated by this, we discuss the implications of a PPISN happening with a close binary companion, and what impact these events have on the formation of merging BBHs through binary evolution. For this, we have computed a set of models of metal-poor ($Z_\odot/10$) single helium stars using the \texttt{MESA} software instrument. For PPISN progenitors with pre-pulse masses $>50M_\odot$ we find that, after a pulse, heat deposited throughout the layers of the star that remain bound cause it to expand to more than $100R_\odot$ for periods of $10^2-10^4\;$~yrs depending on the mass of the progenitor. This results in long-lived phases of Roche-lobe overflow or even common-envelope events if there is a close binary companion, leading to additional electromagnetic transients associated to PPISN eruptions. If we ignore the effect of these interactions, we find that mass loss from PPISNe reduces the final black hole spin by $\sim 30\%$, induces eccentricities below the threshold of detectability of the LISA observatory, and can produce a double-peaked distribution of measured chirp masses in BBH mergers observed by ground-based detectors.


title: Mass Splitting in a Linear Sigma Model

abstract: We derived the tree level spectrum to an extension to the linear sigma model describing an EFT for an $SU(3)_c$ gauge theory with $N_f$ flavors of fermions and $N_1$ fermions have a mass $m_l$ and $N_2$ fermions have a mass $m_h$. We examined the effects of a small mass splitting on single mass data for 8 and 12 flavors of fermions corresponding to the unperturbed case. Our intention is to encourage more simulations of split mass theories with 8, 10 and 12 flavors of fermions.


title: Extracting Linguistic Resources from the Web for Concept-to-Text Generation

abstract: Many concept-to-text generation systems require domain-specific linguistic resources to produce high quality texts, but manually constructing these resources can be tedious and costly. Focusing on NaturalOWL, a publicly available state of the art natural language generator for OWL ontologies, we propose methods to extract from the Web sentence plans and natural language names, two of the most important types of domain-specific linguistic resources used by the generator. Experiments show that texts generated using linguistic resources extracted by our methods in a semi-automatic manner, with minimal human involvement, are perceived as being almost as good as texts generated using manually authored linguistic resources, and much better than texts produced by using linguistic resources extracted from the relation and entity identifiers of the ontology.


title: Calculation of the Uncertainties in the Planetary Precessions with the Recent EPM2017 Ephemerides and their Use in Fundamental Physics and Beyond

abstract: I tentatively compile the formal uncertainties in the secular rates of change of the orbital elements $a,~e,~I,~\Omega$ and $\varpi$ of the planets of the solar system from the recently released formal errors in $a$ and the nonsingular elements $h,~k,~p$ and $q$ estimated for the same bodies with the EPM2017 ephemerides by E.\,V. Pitjeva and N.\,P. Pitjev. The highest accuracies occur for the inner planets and Saturn in view of the extensive use of radiotechnical data collected over the last decades. For the inclination $I$, node $\Omega$ and perihelion $\varpi$ of Mercury and Mars, I obtain accuracies $\sigma_{\dot I},\,\sigma_{\dot\Omega},\,\sigma_{\dot \varpi}\simeq 1-10\,\mu\textrm{as~cty}^{-1}$, while for Saturn they are $\sigma_{\dot I},\,\sigma_{\dot\Omega},\,\sigma_{\dot \varpi}\simeq 10\,\mu\textrm{as~cty}^{-1}-1\,\textrm{mas~cty}^{-1}$. As far as the semimajor axis $a$ is concerned, its rates for the inner planets are accurate to the $\sigma_{\dot a}\simeq 1-100\,\textrm{mm~cty}^{-1}$ level, while for Saturn I obtain $\sigma_{\dot a}\simeq 17\,\textrm{m~cty}^{-1}$. In terms of the parameterized post-Newtonian (PPN) parameters $\beta$ and $\gamma$, a formal error as little as $8\,\mu\textrm{as~cty}^{-1}$ for the Hermean apsidal rate corresponds to a $\simeq 2\times 10^{-7}$ bias in the combination $\left(2 + 2\gamma- \beta\right)/3$ parameterizing the Schwarzschild-type periehlion precession of Mercury. The realistic uncertainties of the planetary precessions may be up to one order of magnitude larger. I discuss their potential multiple uses in fundamental physics, astronomy, and planetology.


title: Simulating 3-D Stellar Hydrodynamics using PPM and PPB Multifluid Gas Dynamics on CPU and CPU+GPU Nodes

abstract: The special computational challenges of simulating 3-D hydrodynamics in deep stellar interiors are discussed, and numerical algorithmic responses described. Results of recent simulations carried out at scale on the NSF's Blue Waters machine at the University of Illinois are presented, with a special focus on the computational challenges they address. Prospects for future work using GPU-accelerated nodes such as those on the DoE's new Summit machine at Oak Ridge National Laboratory are described, with a focus on numerical algorithmic accommodations that we believe will be necessary.


title: Geometric Flows of G_2 Structures

abstract: Geometric flows have proved to be a powerful geometric analysis tool, perhaps most notably in the study of 3-manifold topology, the differentiable sphere theorem, Hermitian-Yang-Mills connections and canonical Kaehler metrics. In the context of G_2 geometry, there are several geometric flows which arise. Each flow provides a potential means to study the geometry and topology associated with a given class of G_2 structures. We will introduce these flows, and describe some of the key known results and open problems in the field.


title: Sharp error estimates for spline approximation: explicit constants, $n$-widths, and eigenfunction convergence

abstract: In this paper we provide a priori error estimates in standard Sobolev (semi-)norms for approximation in spline spaces of maximal smoothness on arbitrary grids. The error estimates are expressed in terms of a power of the maximal grid spacing, an appropriate derivative of the function to be approximated, and an explicit constant which is, in many cases, sharp. Some of these error estimates also hold in proper spline subspaces, which additionally enjoy inverse inequalities. Furthermore, we address spline approximation of eigenfunctions of a large class of differential operators, with a particular focus on the special case of periodic splines. The results of this paper can be used to theoretically explain the benefits of spline approximation under $k$-refinement by isogeometric discretization methods. They also form a theoretical foundation for the outperformance of smooth spline discretizations of eigenvalue problems that has been numerically observed in the literature, and for optimality of geometric multigrid solvers in the isogeometric analysis context.


title: Imaging emergent heavy Dirac fermions of a topological Kondo insulator

abstract: Kondo insulators are primary candidates in the search for strongly correlated topological quantum phases, which may host topological order, fractionalization, and non-Abelian statistics. Within some Kondo insulators, the hybridization gap is predicted to protect a nontrivial topological invariant and to harbor emergent heavy Dirac fermion surface modes. We use high-energy-resolution spectroscopic imaging in real and momentum space on the Kondo insulator, SmB$_6$. On cooling through $T^*_{\Delta}\approx$ 35 K we observe the opening of an insulating gap that expands to $\Delta\approx$ 10 meV at 2 K. Within the gap, we image the formation of linearly dispersing surface states with effective masses reaching $m^* = (410\pm20)m_e$. We thus demonstrate existence of a strongly correlated topological Kondo insulator phase hosting the heaviest known Dirac fermions.


title: Development of a 3-D Position Sensitive Neutron Detector Based on Organic Scintillators with Double Side SiPM Readout

abstract: A 3-D position sensitive neutron detector is being developed based on a plastic scintillator array. A double side SiPM readout is used to determine the depth of interaction (DOI) in each scintillator unit. In the preliminary test, the DOI in a 254 x 6 x 6 mm3 SP101 plastic scintillator is measured at different positions using a collimated Co-60 source. The SiPM (KETEK PM6660) signals are recorded by a 2.5 GS/s digital oscilloscope. The DOI results are calculated using both the amplitude and the temporal information. Position resolutions (FWHM) of 2.5 cm and 6.6 cm are realized, respectively. A detector based on a 2-D array is capable of recording the 3-D position information of the incident neutron. The 3-D detector is to be used together with a neutron time projection chamber as a directional fast neutron detector. According to the simulation results, the angular resolution (8 degree FWHM) is much better than that of a typical neutron scatter camera.


title: Multiple Measurement Vectors Problem: A Decoupling Property and its Applications

abstract: We study a Compressed Sensing (CS) problem known as Multiple Measurement Vectors (MMV) problem, which arises in joint estimation of multiple signal realizations when the signal samples have a common (joint) sparse support over a fixed known dictionary. Although there is a vast literature on the analysis of MMV, it is not yet fully known how the number of signal samples and their statistical correlations affects the performance of the joint estimation in MMV. Moreover, in many instances of MMV the underlying sparsifying dictionary may not be precisely known, and it is still an open problem to quantify how the dictionary mismatch may affect the estimation performance. In this paper, we focus on $\ell_{2,1}$-norm regularized least squares ($\ell_{2,1}$-LS) as a well-known and widely-used MMV algorithm in the literature. We prove an interesting decoupling property for $\ell_{2,1}$-LS, where we show that it can be decomposed into two phases: i) use all the signal samples to estimate the signal covariance matrix (coupled phase), ii) plug in the resulting covariance estimate as the true covariance matrix into the Minimum Mean Squared Error (MMSE) estimator to reconstruct each signal sample individually (decoupled phase). As a consequence of this decomposition, we are able to provide further insights on the performance of $\ell_{2,1}$-LS for MMV. In particular, we address how the signal correlations and dictionary mismatch affects its performance. Moreover, we show that by using the decoupling property one can obtain a variety of MMV algorithms with performances even better than that of $\ell_{2,1}$-LS. We also provide numerical simulations to validate our theoretical results.


title: A fast, sensitive, room-temperature graphene nanomechanical bolometer

abstract: Bolometers are a powerful and vital means of detecting light in the IR to THz frequencies, and they have been adopted for a range of uses from astronomical observation to thermal imaging. As uses diversify, there is an increasing demand for faster, more sensitive room-temperature bolometers. To this end, graphene has generated interest because of its miniscule heat capacity and its intrinsic ultra-broadband absorption, properties that would allow it to quickly detect low levels of light of nearly any wavelength. Yet, graphene has disappointed its expectations in traditional electrical bolometry at room temperature, because of its weakly temperature-dependent resistivity and exceptionally high thermal conductivity. Here, we overcome these challenges with a new approach that detects light by tracking the resonance frequency of a graphene nanomechanical resonator. The absorbed light heats up and thermally tensions the resonator, thereby changing its frequency. Using this approach, we achieve a room-temperature noise-equivalent power of 7 pW/Hz^1/2, a value 100 times more sensitive than electrical graphene bolometers, and speeds (1.3 MHz) that greatly surpass state-of-the-art microbolometers.


title: KSZ tomography and the bispectrum

abstract: Several statistics have been proposed for measuring the kSZ effect by combining the small-scale CMB with galaxy surveys. We review five such statistics, and show that they are all mathematically equivalent to the optimal bispectrum estimator of type $\langle ggT \rangle$. Reinterpreting these kSZ statistics as special cases of bispectrum estimation makes many aspects transparent, for example optimally weighting the estimator, or incorporating photometric redshift errors. We analyze the information content of the bispectrum and show that there are two observables: the small-scale galaxy-electron power spectrum $P_{ge}(k_S)$, and the large-scale galaxy-velocity power spectrum $P_{gv}(k)$. The cosmological constraining power of the kSZ arises from its sensitivity to fluctuations on large length scales, where its effective noise level can be much better than galaxy surveys.


title: Constraining local non-Gaussianities with kSZ tomography

abstract: Kinetic Sunyaev Zel'dovich (kSZ) tomography provides a powerful probe of the radial velocity field of matter in the Universe. By cross-correlating a high resolution CMB experiment like CMB S4 and a galaxy survey like DESI or LSST, one can measure the radial velocity field with very high signal to noise over a large volume of the universe. In this paper we show how this measurement can be used to improve constraints on primordial non-Gaussianities of the local type. The velocity field provides a measurement of the unbiased matter perturbations on large scales, which can be cross-correlated with the biased large-scale galaxy density field. This results in sample variance cancellation for a measurement of scale-dependent bias due to a non-zero $f_{NL}$. Using this method we forecast that CMB S4 and LSST combined reach a sensitivity $\sigma_{f_{NL}} \sim 0.5$, which is a factor of three improvement over the sensitivity using LSST alone (without internal sample variance cancellation). We take into account critical systematics like photometric redshifts, the kSZ optical depth degeneracy, and systematics affecting the shape of the galaxy auto-power spectrum and find that these have negligible impact, thus making kSZ tomography a robust probe for primordial non-Gaussianities. We also forecast the impact of mass binning on our constraints. The techniques proposed in this paper could be an important component of achieving the theoretically important threshold of $\sigma_{f_{NL}} \lesssim 1$ with future surveys.


title: Understanding Deep Neural Networks through Input Uncertainties

abstract: Techniques for understanding the functioning of complex machine learning models are becoming increasingly popular, not only to improve the validation process, but also to extract new insights about the data via exploratory analysis. Though a large class of such tools currently exists, most assume that predictions are point estimates and use a sensitivity analysis of these estimates to interpret the model. Using lightweight probabilistic networks we show how including prediction uncertainties in the sensitivity analysis leads to: (i) more robust and generalizable models; and (ii) a new approach for model interpretation through uncertainty decomposition. In particular, we introduce a new regularization that takes both the mean and variance of a prediction into account and demonstrate that the resulting networks provide improved generalization to unseen data. Furthermore, we propose a new technique to explain prediction uncertainties through uncertainties in the input domain, thus providing new ways to validate and interpret deep learning models.


title: Optimal constants in nontrapping resolvent estimates and applications in numerical analysis

abstract: We study the resolvent for nontrapping obstacles on manifolds with Euclidean ends. It is well known that for such manifolds, the outgoing resolvent satisfies $\|\chi R(k) \chi\|_{L^2\to L^2}\leq C{k}^{-1}$ for ${k}>1$, but the constant $C$ has been little studied. We show that, for high frequencies, the constant is bounded above by $2/\pi$ times the length of the longest generalized bicharacteristic of $|\xi|_g^2-1$ remaining in the support of $\chi.$ We show that this estimate is optimal in the case of manifolds without boundary. We then explore the implications of this result for the numerical analysis of the Helmholtz equation.


title: Unsupervised Dimension Selection using a Blue Noise Spectrum

abstract: Unsupervised dimension selection is an important problem that seeks to reduce dimensionality of data, while preserving the most useful characteristics. While dimensionality reduction is commonly utilized to construct low-dimensional embeddings, they produce feature spaces that are hard to interpret. Further, in applications such as sensor design, one needs to perform reduction directly in the input domain, instead of constructing transformed spaces. Consequently, dimension selection (DS) aims to solve the combinatorial problem of identifying the top-$k$ dimensions, which is required for effective experiment design, reducing data while keeping it interpretable, and designing better sensing mechanisms. In this paper, we develop a novel approach for DS based on graph signal analysis to measure feature influence. By analyzing synthetic graph signals with a blue noise spectrum, we show that we can measure the importance of each dimension. Using experiments in supervised learning and image masking, we demonstrate the superiority of the proposed approach over existing techniques in capturing crucial characteristics of high dimensional spaces, using only a small subset of the original features.


title: Topological Quantum Field Theory for Abelian Topological Phases and Loop Braiding Statistics in $(3+1)$-Dimensions

abstract: Topological qauntum field theory(TQFT) is a very powerful theoretical tool to study topological phases and phase transitions. In $2+1$D, it is well known that the Chern-Simons theory captures all the universal topological data of topological phases, e.g., quasi-particle braiding statistics, chiral central charge and even provides us a deep insight for the nature of topological phase transitions. Recently, topological phases of quantum matter are also intensively studied in $3+1$D and it has been shown that loop like excitation obeys the so-called three-loop-braiding statistics. In this paper, we will try to establish a TQFT framework to understand the quantum statistics of particle and loop like excitation in $3+1$D. We will focus on Abelian topological phases for simplicity, however, the general framework developed here is not limited to Abelian topological phases.


title: Angular momentum of fission fragments

abstract: The strong deformation present immediately after scission has consequences for the angular momentum population of the fragments as well as the angular distribution of their decay radiation. We find that the usual spin-cutoff parameterization describes very well the angular momentum distribution associated with the deformation of the fragments at the scission point. Depending on the deformation, its contribution can comparable to the thermal contribution to the angular momentum of the newly formed fragments. The $M$-distribution of the angular momentum is highly polarized and gives rise to large anisotropies in the subsequent gamma cascade. We treat in detail a typical gamma cascade in a daughter nucleus, following usual model assumptions except for the anisotropy of the initial state. In principle, the observed anisotropy can provide information on the relative amounts of deformation and thermal energy present at the scission point.


title: Exceptionally Monadic Error Handling

abstract: We notice that the type of catch :: c a -> (e -> c a) -> c a operator is a special case of monadic bind operator (>>=) :: m a -> (a -> m b) -> m b, the semantics (surprisingly) matches, and this observation has many interesting consequences. For instance, the reader is probably aware that the monadic essence of the (>>=) operator of the error monad $\lambda A.E \lor A$ is to behave like identity monad for "normal" values and to stop on "errors". The unappreciated fact is that handling of said "errors" with a catch operator of the "flipped" "conjoined" error monad $\lambda E.E \lor A$ is, too, a monadic computation that treats still unhandled "errors" as "normal" values and stops when an "error" is finally handled. We show that for an appropriately indexed type of computations such a "conjoined" structure naturally follows from the conventional operational semantics of throw and catch operators. Consequently, we show that this structure uniformly generalizes all conventional monadic error handling mechanisms we are aware of. We also demonstrate several more interesting instances of this structure of which at least bi-indexed monadic parser combinators and conventional exceptions implemented via continuations have immediate practical applications. Finally, we notice that these observations provide surprising perspectives on error handling in general and point to a largely unexplored trail in programming language design space.


title: Targeted stochastic gradient Markov chain Monte Carlo for hidden Markov models with rare latent states

abstract: Markov chain Monte Carlo (MCMC) algorithms for hidden Markov models often rely on the forward-backward sampler. This makes them computationally slow as the length of the time series increases, motivating the recent development of sub-sampling-based approaches. These approximate the full posterior by using small random subsequences of the data at each MCMC iteration within stochastic gradient MCMC. In the presence of imbalanced data resulting from rare latent states, subsequences often exclude rare latent state data, leading to inaccurate inference and prediction/detection of rare events. We propose a targeted sub-sampling (TASS) approach that over-samples observations corresponding to rare latent states when calculating the stochastic gradient of parameters associated with them. TASS uses an initial clustering of the data to construct subsequence weights that reduce the variance in gradient estimation. This leads to improved sampling efficiency, in particular in settings where the rare latent states correspond to extreme observations. We demonstrate substantial gains in predictive and inferential accuracy on real and synthetic examples.


title: Making root cause analysis feasible for large code bases: a solution approach for a climate model

abstract: For large-scale simulation codes with huge and complex code bases, where bit-for-bit comparisons are too restrictive, finding the source of statistically significant discrepancies (e.g., from a previous version, alternative hardware or supporting software stack) in output is non-trivial at best. Although there are many tools for program comprehension through debugging or slicing, few (if any) scale to a model as large as the Community Earth System Model (CESM; trademarked), which consists of more than 1.5 million lines of Fortran code. Currently for the CESM, we can easily determine whether a discrepancy exists in the output using a by now well-established statistical consistency testing tool. However, this tool provides no information as to the possible cause of the detected discrepancy, leaving developers in a seemingly impossible (and frustrating) situation. Therefore, our aim in this work is to provide the tools to enable developers to trace a problem detected through the CESM output to its source. To this end, our strategy is to reduce the search space for the root cause(s) to a tractable size via a series of techniques that include creating a directed graph of internal CESM variables, extracting a subgraph (using a form of hybrid program slicing), partitioning into communities, and ranking nodes by centrality. Runtime variable sampling then becomes feasible in this reduced search space. We demonstrate the utility of this process on multiple examples of CESM simulation output by illustrating how sampling can be performed as part of an efficient parallel iterative refinement procedure to locate error sources, including sensitivity to CPU instructions. By providing CESM developers with tools to identify and understand the reason for statistically distinct output, we have positively impacted the CESM software development cycle and, in particular, its focus on quality assurance.


title: Two types of axisymmetric helical magnetorotational instability in rotating flows with positive shear

abstract: We reveal and investigate a new type of linear axisymmetric helical magnetorotational instability which is capable of destabilizing viscous and resistive rotational flows with radially increasing angular velocity, or positive shear. This instability is double-diffusive by nature and is different from the more familiar helical magnetorotational instability, operating at positive shear above the Liu limit, in that it works instead for a wide range of the positive shear when ${\rm (i)}$ a combination of axial/poloidal and azimuthal/toroidal magnetic fields is applied and ${\rm (ii)}$ the magnetic Prandtl number is not too close to unity. We study this instability first with radially local WKB analysis and then confirm its existence using a global stability analysis of the magnetized flow between two rotating cylinders with conducting or insulating boundaries. From an experimental point of view, we also demonstrate the presence of the new instability in a magnetized viscous and resistive Taylor-Couette flow with positive shear for such values of the flow parameters, which can be realized in upcoming experiments at the DRESDYN facility. Finally, this instability might have implications for the dynamics of the equatorial parts of the solar tachocline and dynamo action there, since the above two necessary conditions for the instability to take place are satisfied in this region. Our global stability calculations for the tachocline-like configuration, representing a thin rotating cylindrical layer with the appropriate boundary conditions -- conducting inner and insulating outer cylinders -- and the values of the flow parameters, indicate that it can indeed arise in this case with a characteristic growth time comparable to the solar cycle period.

categories: physics.flu-dyn astro-ph.EP astro-ph.SR physics.plasm-ph


title: Scientific Domain Knowledge Improves Exoplanet Transit Classification with Deep Learning

abstract: Space-based missions such as Kepler, and soon TESS, provide large datasets that must be analyzed efficiently and systematically. Recent work by Shallue & Vanderburg (2018) successfully used state-of-the-art deep learning models to automatically classify Kepler transit signals as either exoplanets or false positives; our application of their model yielded 95.8% accuracy and 95.5% average precision. Here we expand upon that work by including additional scientific domain knowledge into the network architecture and input representations to significantly increase overall model performance to 97.5% accuracy and 98.0% average precision. Notably, we achieve 15-20% gains in recall for the lowest signal-to-noise transits that can correspond to rocky planets in the habitable zone. We input into the network centroid time-series information derived from Kepler data plus key stellar parameters taken from the Kepler DR25 catalogue. We also implement data augmentation techniques to alleviate model over-fitting. These improvements allow us to drastically reduce the size of the model, while still maintaining improved performance; smaller models are better for generalization, for example from Kepler to TESS data. This work illustrates the importance of including expert domain knowledge in even state-of-the-art deep learning models when applying them to scientific research problems that seek to identify weak signals in noisy data. This classification tool will be especially useful for upcoming space-based photometry missions focused on finding small planets, such as TESS and PLATO.


title: Anisotropies in the astrophysical gravitational-wave background: The impact of black hole distributions

abstract: We use population inference to explore the impact that uncertainties in the distribution of binary black holes (BBH) have on the astrophysical gravitational-wave background (AGWB). Our results show that the AGWB monopole is sensitive to the nature of the BBH population (particularly the local merger rate), while the anisotropic $C_\ell$ spectrum is only modified to within a few percent, at a level which is insignificant compared to other sources of uncertainty (such as cosmic variance). This is very promising news for future observational studies of the AGWB, as it shows that (i) the monopole can be used as a new probe of the population of compact objects throughout cosmic history, complementary to direct observations by LIGO and Virgo and (ii) we are able to make surprisingly robust predictions for the $C_\ell$ spectrum, even with only very approximate knowledge of the black hole population. As a result, the AGWB anisotropies have enormous potential as a new probe of the large-scale structure of the Universe, and of late-Universe cosmology in general.


title: Gaia-DR2 extended kinematical maps. Part I: Method and application

abstract: CONTEXT. The Gaia Collaboration has used Gaia-DR2 sources with six-dimensional (6D) phase space information to derive kinematical maps within 5 kpc of the Sun, which is a reachable range for stars with relative error in distance lower than 20%. AIMS. Here we aim to extend the range of distances by a factor of two to three, thus adding the range of Galactocentric distances between 13 kpc and 20 kpc to the previous maps, with their corresponding error and root mean square values. METHODS. We make use of the whole sample of stars of Gaia-DR2 including radial velocity measurements, which consists in more than seven million sources, and we apply a statistical deconvolution of the parallax errors based on the Lucy's inversion method of the Fredholm integral equations of the first kind, without assuming any prior. RESULTS. The new extended maps provide lots of new and corroborated information about the disk kinematics: significant departures of circularity in the mean orbits with radial Galactocentric velocities between -20 and +20 km/s and vertical velocities between -10 and +10 km/s; variations of the azimuthal velocity with position; asymmetries between the northern and the southern Galactic hemispheres, especially towards the anticenter that includes a larger azimuthal velocity in the south; and others. CONCLUSIONS. These extended kinematical maps can be used to investigate the different dynamical models of our Galaxy, and we will present our own analyses in the forthcoming second part of this paper. At present, it is evident that the Milky Way is far from a simple stationary configuration in rotational equilibrium, but is characterized by streaming motions in all velocity components with conspicuous velocity gradients.


title: Production of a degenerate Fermi-Fermi mixture of dysprosium and potassium atoms

abstract: We report on the realization of a mixture of fermionic $^{161}$Dy and fermionic $^{40}$K where both species are deep in the quantum-degenerate regime. Both components are spin-polarized in their absolute ground states, and the low temperatures are achieved by means of evaporative cooling of the dipolar dysprosium atoms together with sympathetic cooling of the potassium atoms. We describe the trapping and cooling methods, in particular the final evaporation stage, which leads to Fermi degeneracy of both species. Analyzing cross-species thermalization we obtain an estimate of the magnitude of the inter-species $s$-wave scattering length at low magnetic field. We demonstrate magnetic levitation of the mixture as a tool to ensure spatial overlap of the two components. The properties of the Dy-K mixture make it a very promising candidate to explore the physics of strongly interacting mass-imbalanced Fermi-Fermi mixtures.


title: Particle dynamics simulation of wet granulation in a rotating drum

abstract: We simulate the granulation process of solid spherical particles in the presence of a viscous liquid in a horizontal rotating drum by using molecular dynamics simulations in three dimensions. The numerical approach accounts for the cohesive and viscous effects of the binding liquid, which is assumed to be transported by wet particles and re-distributed homogeneously between wet particles in contact. We investigate the growth of a single granule introduced into the granular bed and the cumulative numbers of accreted and eroded particles as a function of time for a range of values of material parameters such as mean particle size, size polydispersity, friction coefficient, and liquid viscosity. We find that the granule growth is an exponential function of time, reflecting the decrease of the number of free wet particles. The influence of material parameters on the accretion and erosion rates reveals the nontrivial dynamics of the granulation process. It opens the way to a granulation model based on the realistic determination of particle-scale mechanisms of granulation.


title: RELICS: Strong Lensing Analysis of MACS J0417.5-1154 and Predictions for Observing the Magnified High-Redshift Universe with JWST

abstract: Strong gravitational lensing by clusters of galaxies probes the mass distribution at the core of each cluster and magnifies the universe behind it. MACS J0417.5-1154 at z=0.443 is one of the most massive clusters known based on weak lensing, X-ray, and Sunyaev-Zel'dovich analyses. Here we compute a strong lens model of MACS J0417 based on Hubble Space Telescope imaging observations collected, in part, by the Reionization Lensing Cluster Survey (RELICS), and recently reported spectroscopic redshifts from the MUSE instrument on the Very Large Telescope (VLT). We measure an Einstein radius of ThetaE=36'' at z = 9 and a mass projected within 200 kpc of M(200 kpc) = 1.78+0.01-0.03x10**14Msol. Using this model, we measure a ratio between the mass attributed to cluster-member galaxy halos and the main cluster halo of order 1:100. We assess the probability to detect magnified high-redshift galaxies in the field of this cluster, both for comparison with RELICS HST results and as a prediction for the James Webb Space Telescope (JWST) Guaranteed Time Observations upcoming for this cluster. Our lensing analysis indicates that this cluster has similar lensing strength to other clusters in the RELICS program. Our lensing analysis predicts a detection of at least a few z~6-8 galaxies behind this cluster, at odds with a recent analysis that yielded no such candidates in this field. Reliable strong lensing models are crucial for accurately predicting the intrinsic properties of lensed galaxies. As part of the RELICS program, our strong lensing model produced with the Lenstool parametric method is publicly available through the Mikulski Archive for Space Telescopes (MAST).


title: Emergent classical spacetime from microstates of an incipient black hole

abstract: Black holes have an enormous underlying space of microstates, but universal macroscopic physics characterized by mass, charge and angular momentum as well as a causally disconnected interior. This leads two related puzzles: (1) How does the effective factorization of interior and exterior degrees of freedom emerge in gravity?, and (2) How does the underlying degeneracy of states wind up having a geometric realization in the horizon area and in properties of the singularity? We explore these puzzles in the context of an incipient black hole in the AdS/CFT correspondence, the microstates of which are dual to half-BPS states of the $\mathcal{N}=4$ super-Yang-Mills theory. First, we construct a code subspace for this black hole and show how to organize it as a tensor product of a universal macroscopic piece (describing the exterior), and a factor corresponding to the microscopic degrees of freedom (describing the interior). We then study the classical phase space and symplectic form for low-energy excitations around the black hole. On the AdS side, we find that the symplectic form has a new physical degree of freedom at the stretched horizon of the black hole, reminiscent of soft hair, which is absent in the microstates. We explicitly show how such a soft mode emerges from the microscopic phase space in the dual CFT via a canonical transformation and how it encodes partial information about the microscopic degrees of freedom of the black hole.


title: Improving Machine Reading Comprehension with General Reading Strategies

abstract: Reading strategies have been shown to improve comprehension levels, especially for readers lacking adequate prior knowledge. Just as the process of knowledge accumulation is time-consuming for human readers, it is resource-demanding to impart rich general domain knowledge into a deep language model via pre-training. Inspired by reading strategies identified in cognitive science, and given limited computational resources -- just a pre-trained model and a fixed number of training instances -- we propose three general strategies aimed to improve non-extractive machine reading comprehension (MRC): (i) BACK AND FORTH READING that considers both the original and reverse order of an input sequence, (ii) HIGHLIGHTING, which adds a trainable embedding to the text embedding of tokens that are relevant to the question and candidate answers, and (iii) SELF-ASSESSMENT that generates practice questions and candidate answers directly from the text in an unsupervised manner. By fine-tuning a pre-trained language model (Radford et al., 2018) with our proposed strategies on the largest general domain multiple-choice MRC dataset RACE, we obtain a 5.8% absolute increase in accuracy over the previous best result achieved by the same pre-trained model fine-tuned on RACE without the use of strategies. We further fine-tune the resulting model on a target MRC task, leading to an absolute improvement of 6.2% in average accuracy over previous state-of-the-art approaches on six representative non-extractive MRC datasets from different domains (i.e., ARC, OpenBookQA, MCTest, SemEval-2018 Task 11, ROCStories, and MultiRC). These results demonstrate the effectiveness of our proposed strategies and the versatility and general applicability of our fine-tuned models that incorporate these strategies. Core code is available at https://github.com/nlpdata/strategy/.


title: Consistency between ARPES and STM measurements on SmB$_6$

abstract: Strongly correlated topological surface states are promising platforms for next-generation quantum applications, but they remain elusive in real materials. The correlated Kondo insulator SmB$_6$ is one of the most promising candidates, with theoretically predicted heavy Dirac surface states supported by transport and scanning tunneling microscopy (STM) experiments. However, a puzzling discrepancy appears between STM and angle-resolved photoemission (ARPES) experiments on SmB$_6$. Although ARPES detects spin-textured surface states, their velocity is an order of magnitude higher than expected, while the Dirac point -- the hallmark of any topological system -- can only be inferred deep within the bulk valence band. A significant challenge is that SmB$_6$ lacks a natural cleavage plane, resulting in ordered surface domains limited to 10s of nanometers. Here we use STM to show that surface band bending can shift energy features by 10s of meV between domains. Starting from our STM spectra, we simulate the full spectral function as an average over multiple domains with different surface potentials. Our simulation shows excellent agreement with ARPES data, and thus resolves the apparent discrepancy between large-area measurements that average over multiple band-shifted domains and atomically-resolved measurements within a single domain.


title: The inverse Born problem in contextual probability theories: quantum spin and continuous random variables

abstract: We revise contextual probability theory and the associated inverse Born problem, solved by the Quantum-Like Representation Algorithm (QLRA) in case of two discrete random variables. After pointing out the special feature and limitations of QLRA for two binary random variables, we generalize the QLRA procedure to solve the inverse Born problem to the case of three binary random variables. We analyze the quantum spin model in the light of our results and showed that can be easily obtained from a contextual probability model. We furthermore study the inverse Born problem in the case of two continuous random variables, exploiting the general idea underlying QLRA.


title: Model parameter estimation using coherent structure coloring

abstract: Lagrangian data assimilation is a complex problem in oceanic and atmospheric modeling. Tracking drifters in large-scale geophysical flows can involve uncertainty in drifter location, complex inertial effects, and other factors which make comparing them to simulated Lagrangian trajectories from numerical models extremely challenging. Temporal and spatial discretization, factors necessary in modeling large scale flows, also contribute to separation between real and simulated drifter trajectories. The chaotic advection inherent in these turbulent flows tends to separate even closely spaced tracer particles, making error metrics based solely on drifter displacements unsuitable for estimating model parameters. We propose to instead use error in the coherent structure coloring (CSC) field to assess model skill. The CSC field provides a spatial representation of the underlying coherent patterns in the flow, and we show that it is a more robust metric for assessing model accuracy. Through the use of two test cases, one considering spatial uncertainty in particle initialization, and one examining the influence of stochastic error along a trajectory and temporal discretization, we show that error in the coherent structure coloring field can be used to accurately determine single or multiple simultaneously unknown model parameters, whereas a conventional error metric based on error in drifter displacement fails. Because the CSC field enhances the difference in error between correct and incorrect model parameters, error minima in model parameter sweeps become more distinct. The effectiveness and robustness of this method for single and multi-parameter estimation in analytical flows suggests that Lagrangian data assimilation for real oceanic and atmospheric models would benefit from a similar approach.

categories: physics.flu-dyn cs.LG physics.ao-ph physics.data-an stat.ML


title: Numerical calculation of the full two-loop electroweak corrections to muon (g-2)

abstract: Numerical calculation of two-loop electroweak corrections to the muon anomalous magnetic moment ($g$-2) is done based on, on shell renormalization scheme (OS) and free quark model (FQM). The GRACE-FORM system is used to generate Feynman diagrams and corresponding amplitudes. Total 1780 two-loop diagrams and 70 one-loop diagrams composed of counter terms are calculated to get the renormalized quantity. As for the numerical calculation, we adopt trapezoidal rule with Double Exponential method (DE). Linear extrapolation method (LE) is introduced to regularize UV- and IR-divergences and to get finite values. The reliability of our result is guaranteed by several conditions. The sum of one and two loop electroweak corrections in this renormalization scheme becomes $a_\mu^{EW:OS}[1{\rm+}2{\rm -loop}]= 151.2 (\pm 1.0)\times 10^{-11}$, where the error is due to the numerical integration and the uncertainty of input mass parameters and of the hadronic corrections to electroweak loops. By taking the hadronic corrections into account, we get $a_\mu^{EW}[1{\rm+}2 {\rm -loop}]= 152.9 (\pm 1.0)\times 10^{-11}$. It is in agreement with the previous works given in PDG within errors.


title: Structured Parallel Programming

abstract: Based on our previous work on algebraic laws for true concurrency, we design a structured parallel programming language for true concurrency called PPL. Different to most programming languages, PPL has an explicit parallel operator as an essential operator, including its operational, denotational and axiomatic semantics. PPL can structure a truly concurrent graph to a normal form which means that it is possible to implement a compiler for PPL.


title: High-speed standard magneto-rotational instability

abstract: The large Reynolds number asymptotic approximation of the neutral curve of Taylor-Couette flow subject to axial uniform magnetic field is analysed. The flow has been extensively studied since early 90's as the magneto-rotational instability (MRI) occurring in the flow could possibly explain the origin of the instability observed in certain astrophysical objects. Elsewhere the ideal approximation has been used to study high-speed flows, whilst it sometimes produces paradoxical results. For example, ideal flows must be completely stabilised for strong enough applied magnetic field, but on the other hand the vanishing magnetic Prandtl number limit of the stability should be purely hydrodynamic so instability must occur when Rayleigh's stability condition is violated. The first our discovery is that this apparent contradiction can be resolved by showing the abrupt appearance of the hydrodynamic instability at certain critical value of magnetic Prandtl number, which can be found by asymptotically large Reynolds number limit but with long enough wavelength to retain some diffusive effects. The second our finding concerns so-called Velikhov-Chandrasekhar paradox, namely the mismatch of the zero external magnetic field limit of the Velikhov-Chandrasekhar stability criterion and Rayleigh's stability criterion. We show for fully wide gap cases that the high Reynolds number asymptotic analysis of the MRI naturally yields the simple stability condition that describes smooth transition from Rayleigh to Velikhov-Chandrasekhar stability criteria with increasing Lundquist number.


title: WaveGlow: A Flow-based Generative Network for Speech Synthesis

abstract: In this paper we propose WaveGlow: a flow-based network capable of generating high quality speech from mel-spectrograms. WaveGlow combines insights from Glow and WaveNet in order to provide fast, efficient and high-quality audio synthesis, without the need for auto-regression. WaveGlow is implemented using only a single network, trained using only a single cost function: maximizing the likelihood of the training data, which makes the training procedure simple and stable. Our PyTorch implementation produces audio samples at a rate of more than 500 kHz on an NVIDIA V100 GPU. Mean Opinion Scores show that it delivers audio quality as good as the best publicly available WaveNet implementation. All code will be made publicly available online.


title: Deep Net Features for Complex Emotion Recognition

abstract: This paper investigates the influence of different acoustic features, audio-events based features and automatic speech translation based lexical features in complex emotion recognition such as curiosity. Pretrained networks, namely, AudioSet Net, VoxCeleb Net and Deep Speech Net trained extensively for different speech based applications are studied for this objective. Information from deep layers of these networks are considered as descriptors and encoded into feature vectors. Experimental results on the EmoReact dataset consisting of 8 complex emotions show the effectiveness, yielding highest F1 score of 0.85 as against the baseline of 0.69 in the literature.


title: Modelling N2O dynamics of activated sludge biomass under nitrifying and denitrifying conditions: pathway contributions and uncertainty analysis

abstract: Nitrous oxide (N2O) is a potent greenhouse gas emitted during biological wastewater treatment. A pseudo-mechanistic model describing three biological pathways for nitric oxide (NO) and N2O production was calibrated for mixed culture biomass from an activated sludge process using laboratory-scale experiments. The model (NDHA) comprehensively describes N2O producing pathways by both autotrophic ammonium oxidizing bacteria and heterotrophic bacteria. Extant respirometric assays and anaerobic batch experiments were designed to calibrate endogenous and exogenous processes (heterotrophic denitrification and autotrophic ammonium/nitrite oxidation) together with the associated net N2O production. Ten parameters describing heterotrophic processes and seven for autotrophic processes were accurately estimated (variance/mean < 25%). The model predicted NO and N2O dynamics at varying dissolved oxygen, ammonium and nitrite levels and was validated against an independent set of experiments with the same biomass. Aerobic ammonium oxidation experiments at two oxygen levels used for model evaluation (2 and 0.5 mg/L) indicated that both the nitrifier denitrification (42, 64%) and heterotrophic denitrification (7, 17%) pathways increased and dominated N2O production at high nitrite and low oxygen concentrations; while the nitrifier nitrification pathway showed the largest contribution at high dissolved oxygen levels (51, 19%). The uncertainty of the biological parameter estimates was propagated to N2O model outputs via Monte Carlo simulations as 95% confidence intervals. The accuracy of the estimated parameters resulted in a low uncertainty of the N2O emission factors (4.6 +- 0.6% and 1.2 +- 0.1%).


title: Properties of turbulence in the reconnection exhaust: numerical simulations compared with observations

abstract: The properties of the turbulence which develops in the outflows of magnetic reconnection have been investigated using self-consistent plasma simulations, in three dimensions. As commonly observed in space plasmas, magnetic reconnection is characterized by the presence of turbulence. Here we provide a direct comparison of our simulations with reported observations of reconnection events in the magnetotail investigating the properties of the electromagnetic field and the energy conversion mechanisms. In particular, simulations show the development of a turbulent cascade consistent with spacecraft observations, statistics of the the dissipation mechanisms in the turbulent outflows similar to the one observed in reconnection jets in the magnetotail, and that the properties of turbulence vary as a function of the distance from the reconnecting X-line.


title: Low-Dimensional Bottleneck Features for On-Device Continuous Speech Recognition

abstract: Low power digital signal processors (DSPs) typically have a very limited amount of memory in which to cache data. In this paper we develop efficient bottleneck feature (BNF) extractors that can be run on a DSP, and retrain a baseline large-vocabulary continuous speech recognition (LVCSR) system to use these BNFs with only a minimal loss of accuracy. The small BNFs allow the DSP chip to cache more audio features while the main application processor is suspended, thereby reducing the overall battery usage. Our presented system is able to reduce the footprint of standard, fixed point DSP spectral features by a factor of 10 without any loss in word error rate (WER) and by a factor of 64 with only a 5.8% relative increase in WER.


title: Robustly Disentangled Causal Mechanisms: Validating Deep Representations for Interventional Robustness

abstract: The ability to learn disentangled representations that split underlying sources of variation in high dimensional, unstructured data is important for data efficient and robust use of neural networks. While various approaches aiming towards this goal have been proposed in recent times, a commonly accepted definition and validation procedure is missing. We provide a causal perspective on representation learning which covers disentanglement and domain shift robustness as special cases. Our causal framework allows us to introduce a new metric for the quantitative evaluation of deep latent variable models. We show how this metric can be estimated from labeled observational data and further provide an efficient estimation algorithm that scales linearly in the dataset size.


title: Magnetic field constraint in the outskirts of spiral galaxies

abstract: Based on CO(2-1) public data, we study the monoxide oxygen gas excitation conditions and the magnetic field strength of NGC 2841, NGC 3077, NGC 3184, NGC 3351 spiral galaxies. For their galaxy outskirts, we found kinetic temperatures in the range of $35 - 38 ~K$, CO column densities $10^{15} - 10^{16} ~cm^{-2}$ and H2 masses of $4 ~10^6 - 6 ~ 10^8 M_\odot$. An H2 density 1$0^3 cm^{-3}$ is suitable to explain the 2 sigma upper limits of the CO(2-1) line intensity. We constrain the magnetic field strength for our sample of spiral galaxies and their outskirts, evaluating a simplified expression of the magneto-hydrodynamic force equation. Our estimations provide values for the magnetic field strength in the order of $6-31 ~\mu G$.


title: An easily-implemented, block-based fast marching method with superior sequential and parallel performance

abstract: The fast marching method is well-known for its worst-case optimal computational complexity in solving the Eikonal equation, and has been employed in numerous scientific and engineering fields. However, it has barely benefited from fast-advancing multi-core and many-core architectures, due to the challenges presented by its apparent sequential nature. In this paper, we present a straightforward block-based approach for a highly scalable parallelization of the fast marching method on shard-memory computers. Central to our new algorithm is a simplified restarted narrow band strategy, with which the global bound for terminating the front marching is replaced by an incremental one, increasing by a given stride in each restart. It greatly reduces load imbalance among blocks through a synchronized exchanging step after the marching step. Furthermore, simple activation mechanisms are introduced to skip blocks not involved in either step. Thus the new algorithm mainly consists of two loops, each for performing one step on a group of blocks. Notably, it does not contain any data race conditions at all, ideal for a direct loop level parallelization using multiple threads. A systematic performance study is carried out for several point source problems with various speed functions on four grids with up to $1024^3$ points using two different computers. Substantial parallel speedups, e.g., up to 3--5 times the number of cores on a 16-core/32-thread computer and up to two orders of magnitude on a 64-core/256-thread computer, are demonstrated with all cases. As an extra bonus, our new algorithm also gives improved sequential performance in most scenarios. Detailed pseudo-codes are provided to illustrate the modification procedure from the single-block sequential algorithm to the multi-block one in a step-by-step manner.


title: The Sloan Digital Sky Survey Reverberation Mapping Project: Systematic Investigations of Short-Timescale CIV Broad Absorption Line Variability

abstract: We systematically investigate short-timescale ($<$10-day rest-frame) CIV broad absorption-line (BAL) variability to constrain quasar-wind properties and provide insights into BAL-variability mechanisms in quasars. We employ data taken by the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, as the rapid cadence of these observations provides a novel opportunity to probe BAL variability on shorter rest-frame timescales than have previously been explored. In a sample of 27 quasars with a median of 58 spectral epochs per quasar, we have identified 15 quasars ($55^{+18}_{-14}$%), 19 of 37 CIV BAL troughs ($51^{+15}_{-12}$%), and 54 of 1460 epoch pairs ($3.7 \pm 0.5$%) that exhibit significant CIV BAL equivalent-width variability on timescales of less than 10 days in the quasar rest frame. These frequencies indicate that such variability is common among quasars and BALs, though somewhat rare among epoch pairs. Thus, models describing BALs and their behavior must account for variability on timescales down to less than a day in the quasar rest frame. We also examine a variety of spectral characteristics and find that in some cases, BAL variability is best described by ionization-state changes, while other cases are more consistent with changes in covering fraction or column density. We adopt a simple model to constrain the density and radial distance of two outflows appearing to vary by ionization-state changes, yielding outflow density lower limits consistent with previous work.


title: H$_2$O abundances and cloud properties in ten hot giant exoplanets

abstract: Transmission spectroscopy of exoplanets has the potential to provide precise measurements of atmospheric chemical abundances, in particular of hot Jupiters whose large sizes and high temperatures make them conducive to such observations. To date, several transmission spectra of hot Jupiters have revealed low amplitude features of water vapour compared to expectations from cloud-free atmospheres of solar metallicity. The low spectral amplitudes in such atmospheres could either be due to the presence of aerosols that obscure part of the atmosphere or due to inherently low abundances of H$_2$O in the atmospheres. A recent survey of transmission spectra of ten hot Jupiters used empirical metrics to suggest atmospheres with a range of cloud/haze properties but with no evidence for H$_2$O depletion. Here, we conduct a detailed and homogeneous atmospheric retrieval analysis of the entire sample and report the H$_2$O abundances, cloud properties, terminator temperature profiles, and detection significances of the chemical species. Our present study finds that the majority of hot Jupiters have atmospheres consistent with sub-solar H$_2$O abundances at their day-night terminators. The best constrained abundances range from $\mathrm{log(H_2O)}$ of $-5.04^{+0.46}_{-0.30}$ to $-3.16^{+0.66}_{-0.69}$, which compared to expectations from solar-abundance equilibrium chemistry correspond to $0.018^{+0.035}_{-0.009}\times$ solar to $1.40^{+4.97}_{-1.11}\times$ solar. Besides H$_2$O we report statistical constraints on other chemical species and cloud/haze properties, including cloud/haze coverage fractions which range from $0.18^{+0.26}_{-0.12}$ to $0.76^{+0.13}_{-0.15}$. The retrieved H$_2$O abundances suggest sub-solar oxygen and/or super-solar C/O ratios, and can provide important constraints on the formation and migration pathways of hot giant exoplanets.