title: Supporting Expert Close Analysis of Historical Scientific Writings: A Case Study for Near-by Reading

abstract: Distant reading methodologies make use of computational processes to aid in the analysis of large text corpora which might not be pliable to traditional methods of scholarly analysis due to their volume. While these methods have been applied effectively to a variety of types of texts and contexts, they can leave unaddressed the needs of scholars in the humanities disciplines like history, who often engage in close reading of sources. Complementing the close analysis of texts with some of the tools of distant reading, such as visualization, can resolve some of the issues. We focus on a particular category of this intersection---which we refer to as near-by reading---wherein an expert engages in a computer-mediated analysis of a text with which they are familiar. We provide an example of this approach by developing a visual analysis application for the near-by reading of 19th-century scientific writings by J. W. von Goethe and A. P. de Candolle. We show that even the most formal and public texts, such as scientific treatises, can reveal unexpressed personal biases and philosophies that the authors themselves might not have recognized.


title: Polarization-independent single-photon switch based on a fiber-optical Sagnac interferometer for quantum communication networks

abstract: An essential component of future quantum networks is an optical switch capable of dynamically routing single-photons. Here we implement such a switch, based on a fiber-optical Sagnac interferometer design. The routing is implemented with a pair of fast electro-optical telecom phase modulators placed inside the Sagnac loop, such that each modulator acts on an orthogonal polarization component of the single-photons, in order to yield polarization-independent capability that is crucial for several applications. We obtain an average extinction ratio of more than 19 dB between both outputs of the switch. Our experiment is built exclusively with commercial off-the-shelf components, thus allowing direct compatibility with current optical communication systems.


title: ACDC: Weight Sharing in Atom-Coefficient Decomposed Convolution

abstract: Convolutional Neural Networks (CNNs) are known to be significantly over-parametrized, and difficult to interpret, train and adapt. In this paper, we introduce a structural regularization across convolutional kernels in a CNN. In our approach, each convolution kernel is first decomposed as 2D dictionary atoms linearly combined by coefficients. The widely observed correlation and redundancy in a CNN hint a common low-rank structure among the decomposed coefficients, which is here further supported by our empirical observations. We then explicitly regularize CNN kernels by enforcing decomposed coefficients to be shared across sub-structures, while leaving each sub-structure only its own dictionary atoms, a few hundreds of parameters typically, which leads to dramatic model reductions. We explore models with sharing across different sub-structures to cover a wide range of trade-offs between parameter reduction and expressiveness. Our proposed regularized network structures open the door to better interpreting, training and adapting deep models. We validate the flexibility and compatibility of our method by image classification experiments on multiple datasets and underlying network structures, and show that CNNs now maintain performance with dramatic reduction in parameters and computations, e.g., only 5\% parameters are used in a ResNet-18 to achieve comparable performance. Further experiments on few-shot classification show that faster and more robust task adaptation is obtained in comparison with models with standard convolutions.


title: The neutral hydrogen distribution in large-scale haloes from 21-cm intensity maps

abstract: We detect the neutral hydrogen (HI) radial brightness temperature profile in large-scale haloes by stacking 48,430 galaxies selected from the 2dFGRS catalogue onto a set of 21-cm intensity maps obtained with the Parkes radio telescope, spanning a total area of $\sim$1,300 $\text{deg}^2$ on the sky and covering the redshift range $0.06\lesssim z\lesssim 0.10$. Maps are obtained by removing both 10 and 20 foreground modes in the principal component analysis. We perform the stack at the map level and extract the profile from a circularly symmetrised version of the halo emission. We detect the HI halo emission with the significance $12.5\sigma$ for the 10-mode and $13.5\sigma$ for the 20-mode removed maps at the profile peak. We jointly fit for the observed halo mass $M_{\rm v}$ and the normalisation $c_{0,\rm HI}$ for the HI concentration parameter against the reconstructed profiles, using functional forms for the HI halo abundance proposed in the literature. We find $\log_{10}{(M_{\rm v}/\text{M}_{\odot})} = 16.1^{+0.1}_{-0.2}$, $c_{0,\rm HI}=3.5^{+0.7}_{-1.0}$ for the 10-mode and $\log_{10}{(M_{\rm v}/\text{M}_{\odot})} = 16.5^{+0.1}_{-0.2}$, $c_{0,\rm HI}=5.3^{+1.1}_{-1.7}$ for the 20-mode removed maps. These estimates show the detection of the integrated contribution from multiple galaxies located inside very massive haloes. We also consider sub-samples of 13,979 central and 34,361 satellite 2dF galaxies separately, and obtain marginal differences suggesting satellite galaxies are HI-richer. This work shows for the first time the feasibility of testing theoretical models for the HI halo content directly on profiles extracted from 21-cm maps and opens future possibilities for exploiting upcoming HI intensity-mapping data.


title: On Communication Compression for Distributed Optimization on Heterogeneous Data

abstract: Lossy gradient compression, with either unbiased or biased compressors, has become a key tool to avoid the communication bottleneck in centrally coordinated distributed training of machine learning models. We analyze the performance of two standard and general types of methods: (i) distributed quantized SGD (D-QSGD) with arbitrary unbiased quantizers and (ii) distributed SGD with error-feedback and biased compressors (D-EF-SGD) in the heterogeneous (non-iid) data setting. Our results indicate that D-EF-SGD is much less affected than D-QSGD by non-iid data, but both methods can suffer a slowdown if data-skewness is high. We further study two alternatives that are not (or much less) affected by heterogenous data distributions: first, a recently proposed method that is effective on strongly convex problems, and secondly, we point out a more general approach that is applicable to linear compressors only but effective in all considered scenarios.


title: Poset topology of $s$-weak order via SB-labelings

abstract: Ceballos and Pons generalized weak order on permutations to a partial order on certain labeled trees, thereby introducing a new class of lattices called $s$-weak order. They also generalized the Tamari lattice by defining a particular sublattice of $s$-weak order called the $s$-Tamari lattice. We prove that the homotopy type of each open interval in $s$-weak order and in the $s$-Tamari lattice is either a ball or sphere. We do this by giving $s$-weak order and the $s$-Tamari lattice a type of edge labeling known as an SB-labeling. We characterize which intervals are homotopy equivalent to spheres and which are homotopy equivalent to balls; we also determine the dimension of the spheres for the intervals yielding spheres.


title: Online Learning of Parameterized Uncertain Dynamical Environments with Finite-sample Guarantees

abstract: We present a novel online learning algorithm for a class of unknown and uncertain dynamical environments that are fully observable. First, we obtain a novel probabilistic characterization of systems whose mean behavior is known but which are subject to additive, unknown subGaussian disturbances. This characterization relies on recent concentration of measure results and is given in terms of ambiguity sets. Second, we extend the results to environments whose mean behavior is also unknown but described by a parameterized class of possible mean behaviors. Our algorithm adapts the ambiguity set dynamically by learning the parametric dependence online, and retaining similar probabilistic guarantees with respect to the additive, unknown disturbance. We illustrate the results on a differential-drive robot subject to environmental uncertainty.


title: Visualizing the Loss Landscape of Actor Critic Methods with Applications in Inventory Optimization

abstract: Continuous control is a widely applicable area of reinforcement learning. The main players of this area are actor-critic methods that utilize policy gradients of neural approximators as a common practice. The focus of our study is to show the characteristics of the actor loss function which is the essential part of the optimization. We exploit low dimensional visualizations of the loss function and provide comparisons for loss landscapes of various algorithms. Furthermore, we apply our approach to multi-store dynamic inventory control, a notoriously difficult problem in supply chain operations, and explore the shape of the loss function associated with the optimal policy. We modelled and solved the problem using reinforcement learning while having a loss landscape in favor of optimality.


title: Three Dimensional Simulations of Solar Wind Preconditioning and the 23 July 2012 Interplanetary Coronal Mass Ejection

abstract: Predicting the large-scale eruptions from the solar corona and their propagation through interplanetary space remains an outstanding challenge in solar- and helio-physics research. In this article, we describe three dimensional magnetohydrodynamic simulations of the inner heliosphere leading up to and including the extreme interplanetary coronal mass ejection (ICME) of 23 July 2012, developed using the code PLUTO. The simulations are driven using the output of coronal models for Carrington rotations 2125 and 2126 and, given the uncertainties in the initial conditions, are able to reproduce an event of comparable magnitude to the 23 July ICME, with similar velocity and density profiles at 1 au. The launch-time of this event is then varied with regards to an initial 19 July ICME and the effects of solar wind preconditioning are found to be significant for an event of this magnitude and to decrease over a time-window consistent with the ballistic refilling of the depleted heliospheric sector. These results indicate that the 23 July ICME was mostly unaffected by events prior, but would have travelled even faster had it erupted closer in time to the 19 July event where it would have experienced even lower drag forces. We discuss this systematic study of solar wind preconditioning in the context of space weather forecasting.

categories: astro-ph.SR astro-ph.EP physics.plasm-ph physics.space-ph


title: Crystal orientation and detector distance effects on resolving pseudosymmetry by electron backscatter diffraction

abstract: Accurately indexing pseudosymmetric materials has long proven challenging for electron backscatter diffraction. The recent emergence of intensity-based indexing approaches promises an enhanced ability to resolve pseudosymmetry compared to traditional Hough-based indexing approaches. However, little work has been done to understand the effects of sample position and orientation on the ability to resolve pseudosymmetry, especially for intensity-based indexing approaches. Thus, in this work we quantitatively investigate the effects of crystal orientation and detector distance in a model tetragonal ZrO2 (c/a=1.0185) material. We identify orientations that are easiest and most difficult to correctly index, characterize the effect of detector distance on indexing confidence, and analyze these trends based on the appearance of specific zone axes in the diffraction patterns. Our findings also point to the clear benefit of shorter detector distances for resolving pseudosymmetry using intensity-based indexing approaches.


title: Note on the Labelled tree graphs

abstract: In the CHY-frame for the tree-level amplitudes, the bi-adjoint scalar theory has played a fundamental role because it gives the on-shell Feynman diagrams for all other theories. Recently, an interesting generalization of the bi-adjoint scalar theory has been given in arXiv:1708.08701 by the "Labelled tree graphs", which carries a lot of similarity comparing to the bi-adjoint scalar theory. In this note, we have investigated the Labelled tree graphs from two different angels. In the first part of the note, we have shown that we can organize all cubic Feynman diagrams produces by the Labelled tree graphs to the "effective Feynman diagrams". In the new picture, the pole structure of the whole theory is more manifest. In the second part, we have generalized the action of "picking pole" in the bi-adjoint scalar theory to general CHY-integrands which produce only simple poles.


title: Role of Nanoscale Interfacial Proximity in Contact Freezing in Water

abstract: Contact freezing is a mode of atmospheric ice nucleation in which a collision between a dry ice nucleating particle (INP) and a water droplet results in considerably faster heterogeneous nucleation. The molecular mechanism of such enhancement is, however, still a mystery. While earlier studies had attributed it to collision-induced transient perturbations, recent experiments point to the pivotal role of nanoscale proximity of the INP and the free interface. By simulating heterogeneous nucleation of ice within INP-supported nanofilms of two model water-like tetrahedral liquids, we demonstrate that such nanoscale proximity is sufficient for inducing rate increases commensurate with those observed in contact freezing experiments, but only if the free interface has a tendency to enhance homogeneous nucleation. Water is suspected of possessing this latter property, known as surface freezing propensity. Our findings therefore establish a connection between surface freezing propensity and kinetic enhancement during contact nucleation. We also observe that faster nucleation proceeds through a mechanism markedly distinct from classical heterogeneous nucleation, involving the formation of hourglass-shaped crystalline nuclei that conceive at either interface, and that have a lower free energy of formation due to the nanoscale proximity of the interfaces and the modulation of the free interfacial structure by the INP. In addition to providing valuable insights into the physics of contact nucleation, our findings can assist in improving the accuracy of heterogeneous nucleation rate measurements in experiments, and in advancing our understanding of ice nucleation on nonuniform surfaces such as organic, polymeric and biological materials.

categories: cond-mat.soft cond-mat.stat-mech physics.ao-ph physics.chem-ph physics.comp-ph


title: Class Interference Regularization

abstract: Contrastive losses yield state-of-the-art performance for person re-identification, face verification and few shot learning. They have recently outperformed the cross-entropy loss on classification at the ImageNet scale and outperformed all self-supervision prior results by a large margin (SimCLR). Simple and effective regularization techniques such as label smoothing and self-distillation do not apply anymore, because they act on multinomial label distributions, adopted in cross-entropy losses, and not on tuple comparative terms, which characterize the contrastive losses. Here we propose a novel, simple and effective regularization technique, the Class Interference Regularization (CIR), which applies to cross-entropy losses but is especially effective on contrastive losses. CIR perturbs the output features by randomly moving them towards the average embeddings of the negative classes. To the best of our knowledge, CIR is the first regularization technique to act on the output features. In experimental evaluation, the combination of CIR and a plain Siamese-net with triplet loss yields best few-shot learning performance on the challenging tieredImageNet. CIR also improves the state-of-the-art technique in person re-identification on the Market-1501 dataset, based on triplet loss, and the state-of-the-art technique in person search on the CUHK-SYSU dataset, based on a cross-entropy loss. Finally, on the task of classification CIR performs on par with the popular label smoothing, as demonstrated for CIFAR-10 and -100.


title: A Deep Learning Approach to Tongue Detection for Pediatric Population

abstract: Children with severe disabilities and complex communication needs face limitations in the usage of access technology (AT) devices. Conventional ATs (e.g., mechanical switches) can be insufficient for nonverbal children and those with limited voluntary motion control. Automatic techniques for the detection of tongue gestures represent a promising pathway. Previous studies have shown the robustness of tongue detection algorithms on adult participants, but further research is needed to use these methods with children. In this study, a network architecture for tongue-out gesture recognition was implemented and evaluated on videos recorded in a naturalistic setting when children were playing a video-game. A cascade object detector algorithm was used to detect the participants' faces, and an automated classification scheme for tongue gesture detection was developed using a convolutional neural network (CNN). In evaluation experiments conducted, the network was trained using adults and children's images. The network classification accuracy was evaluated using leave-one-subject-out cross-validation. Preliminary classification results obtained from the analysis of videos of five typically developing children showed an accuracy of up to 99% in predicting tongue-out gestures. Moreover, we demonstrated that using only children data for training the classifier yielded better performance than adult's one supporting the need for pediatric tongue gesture datasets.


title: Phoenix and Peregrine Ion Traps

abstract: The Phoenix and Peregrine ion traps are micro-fabricated surface-electrode ion traps based on silicon technology. Both are linear traps using a symmetric 6-rail design with segmented inner and outer control electrodes. The traps are fabricated on Sandia's High Optical Access (HOA) platform to provide good optical access skimming the trap surface. They are packaged in custom ceramic pin or land grid array packages using a 2.54 mm pitch. The Peregrine trap is a surface trap with all electrodes in one plane. The Phoenix trap has the same layout, but with a central through-substrate slot and its inner control electrodes are at a lower metal level. Both traps provide means to measure the substrate temperature and to heat the device by means of integrated aluminum and tungsten wires.

categories: physics.app-ph physics.atom-ph physics.ins-det quant-ph


title: Loss rates for high-$n$, $49\lesssim n \lesssim150$, 5sns($^{3}$S$_{1}$) Rydberg atoms excited in an $^{84}$Sr Bose-Einstein condensate

abstract: Measurements of the loss rates for strontium n$^{3}$S$_{1}$ Rydberg atoms excited in a dense BEC are presented for values of principal quantum number $n$ in the range $49\lesssim n \lesssim 150$ and local atom densities of $\sim1$ to $3\times10^{14}$cm$^{-3}$. Two main processes contribute to loss, associative ionization and state-changing. The relative importance of these two loss channels is investigated and their $n$- and density-dependences are discussed using a model in which Rydberg atom loss is presumed to involve a close collision between the Rydberg core ion and ground-state atoms. The present measurements are compared to earlier results obtained using rubidium Rydberg atoms. For both species the observed loss rates are sizable, $\sim10^{5}-10^{6}$s$^{-1}$, and limit the time scales over which measurements involving Rydberg atoms immersed in quantum degenerate gases can be conducted.


title: The Area Under the ROC Curve as a Measure of Clustering Quality

abstract: The Area Under the the Receiver Operating Characteristics (ROC) Curve, referred to as AUC, is a well-known performance measure in the supervised learning domain. Due to its compelling features, it has been employed in a number of studies to evaluate and compare the performance of different classifiers. In this work, we explore AUC as a performance measure in the unsupervised learning domain, more specifically, in the context of cluster analysis. In particular, we elaborate on the use of AUC as an internal/relative measure of clustering quality, which we refer to as Area Under the Curve for Clustering (AUCC). We show that the AUCC of a given candidate clustering solution has an expected value under a null model of random clustering solutions, regardless of the size of the dataset and, more importantly, regardless of the number or the (im)balance of clusters under evaluation. In addition, we elaborate on the fact that, in the context of internal/relative clustering validation as we consider, AUCC is actually a linear transformation of the Gamma criterion from Baker and Hubert (1975), for which we also formally derive a theoretical expected value for chance clusterings. We also discuss the computational complexity of these criteria and show that, while an ordinary implementation of Gamma can be computationally prohibitive and impractical for most real applications of cluster analysis, its equivalence with AUCC actually unveils a much more efficient algorithmic procedure. Our theoretical findings are supported by experimental results. These results show that, in addition to an effective and robust quantitative evaluation provided by AUCC, visual inspection of the ROC curves themselves can be useful to further assess a candidate clustering solution from a broader, qualitative perspective as well.


title: Are multi-quasiparticle interactions important in molecular ionization?

abstract: Photo-emission spectroscopy directly probes individual electronic states, ranging from single excitations to high-energy satellites, which simultaneously represent multiple quasiparticles (QPs) and encode information about electronic correlation. First-principles description of the spectra requires an efficient and accurate treatment of all many-body effects. This is especially challenging for inner valence excitations where the single QP picture breaks down. Here, we provide the full valence spectra of small closed-shell molecules, exploring the independent and interacting quasiparticle regimes, computed with the fully-correlated adaptive sampling configuration interaction (ASCI) method. We critically compare these results to calculations with the many-body perturbation theory, based on the $GW$ and vertex corrected $GW\Gamma$ approaches. The latter explicitly accounts for two-QP quantum interactions, which have been often neglected. We demonstrate that for molecular systems, the vertex correction universally improves the theoretical spectra, and it is crucial for accurate prediction of QPs as well as capturing the rich satellite structures of high-energy excitations. $GW\Gamma$ offers a unified description across all relevant energy scales. Our results suggest that the multi-QP regime corresponds to dynamical correlations, which can be described via perturbation theory.


title: Qualitative properties for solutions to subcritical fourth order systems

abstract: We prove some qualitative properties for singular solutions to a class of strongly coupled system involving a Gross--Pitaevskii-type nonlinearity. Our main theorems are vectorial fourth order counterparts of the classical results of [J. Serrin, Acta Math. (1964)], [P.-L. Lions, J. Differential Equations (1980)], [P. Aviles, Comm. Math. Phys. (1987)], and [B. Gidas and J. Spruck, Comm. Pure Appl. Math. (1981)]. On the technical level, we use the moving sphere method to classify the limit blow-up solutions to our system. Besides, applying asymptotic analysis techniques, we show that these solutions are indeed the local models near the isolated singularity. We also introduce a new fourth order nonautonomous Pohozaev functional, whose monotonicity properties yield improvement for the asymptotics results due to [R. Soranzo, Potential Anal. (1997)].


title: Exact Statistical Thermodynamics of the Pseudospin-1 System on the Diced Lattice

abstract: In this work we analyze the thermodynamic properties of the pseudospin-1 Hamiltonian on the two-dimensional {\cal{T}} -3 or Diced Lattice. Starting from the Partition function, we obtain the Grand ensemble thermodynamic potential, entropy and specific heat exactly and in the degenerate and non-degenerate regimes.


title: Tri-modality Cavitation Mapping in Shock Wave Lithotripsy

abstract: Shock wave lithotripsy (SWL) has been widely used for non-invasive treatment of kidney stones. Cavitation plays an important role in stone fragmentation, yet may also contribute to renal injury during SWL. It is therefore crucial to determine the spatiotemporal distributions of cavitation activities to maximize stone fragmentation while minimizing tissue injury. Traditional cavitation detection methods include high-speed optical imaging, active cavitation mapping (ACM), and passive cavitation mapping (PCM). While each of the three methods provides unique information about the dynamics of the bubbles, PCM has most practical applications in biological tissues. To image the dynamics of cavitation bubble collapse, we previously developed a sliding-window PCM (SW-PCM) method to identify each bubble collapse with high temporal and spatial resolution. To further validate and optimize the SW-PCM method, in this work, we have developed tri-modality cavitation imaging that includes 3D high-speed optical imaging, ACM, and PCM seamlessly integrated in a single system. Using the tri-modality system, we imaged and analyzed laser-induced single cavitation bubbles in both free and constricted space and shockwave-induced cavitation clusters. Collectively, our results have demonstrated the high reliability and spatial-temporal accuracy of the SW-PCM approach, which paves the way for future in vivo applications on large animals and humans in SWL.


title: Main Sequence Scatter is Real: The Joint Dependence of Galaxy Clustering on Star Formation and Stellar Mass

abstract: We present new measurements of the clustering of stellar mass-complete samples of $\sim40,000$ SDSS galaxies at $z\sim0.03$ as a joint function of stellar mass and specific star formation rate (sSFR). Our results confirm what Coil et al. (2017) find at $z\sim0.7$: galaxy clustering is a stronger function of sSFR at fixed stellar mass than of stellar mass at fixed sSFR. We also find that galaxies above the star-forming main sequence (SFMS) with higher sSFR are less clustered than galaxies below the SFMS with lower sSFR, at a given stellar mass. A similar trend is present for quiescent galaxies. This confirms that main sequence scatter, and scatter within the quiescent sequence, is physically connected to the large-scale cosmic density field. We compare the resulting galaxy bias versus sSFR, and relative bias versus sSFR ratio, for different galaxy samples across ${0<z<1.2}$ to mock galaxy catalogs based on the empirical galaxy evolution model of Behroozi et al. (2019). This model fits PRIMUS and DEEP2 clustering data well at intermediate redshift, but agreement with SDSS is not as strong. We show that increasing the correlation between galaxy SFR and halo accretion rate at $z\sim0$ in the model substantially improves agreement with SDSS data. Mock catalogs suggest that central galaxies contribute substantially to the dependence of clustering on sSFR at a given stellar mass and that the signal is not simply an effect of satellite galaxy fraction differences with sSFR. Our results are highly constraining for galaxy evolution models and show that the stellar-to-halo mass relation (SHMR) depends on sSFR.


title: Video Moment Retrieval via Natural Language Queries

abstract: In this paper, we propose a novel method for video moment retrieval (VMR) that achieves state of the arts (SOTA) performance on R@1 metrics and surpassing the SOTA on the high IoU metric (R@1, IoU=0.7). First, we propose to use a multi-head self-attention mechanism, and further a cross-attention scheme to capture video/query interaction and long-range query dependencies from video context. The attention-based methods can develop frame-to-query interaction and query-to-frame interaction at arbitrary positions and the multi-head setting ensures the sufficient understanding of complicated dependencies. Our model has a simple architecture, which enables faster training and inference while maintaining . Second, We also propose to use multiple task training objective consists of moment segmentation task, start/end distribution prediction and start/end location regression task. We have verified that start/end prediction are noisy due to annotator disagreement and joint training with moment segmentation task can provide richer information since frames inside the target clip are also utilized as positive training examples. Third, we propose to use an early fusion approach, which achieves better performance at the cost of inference time. However, the inference time will not be a problem for our model since our model has a simple architecture which enables efficient training and inference.


title: A note on the geometry of certain classes of linear operators

abstract: In this note we introduce a new technique to answer an issue posed in [7] concerning geometric properties of the set of non-surjective linear operators. We also extend and improve a related result from the same paper.


title: Quarkonium Semiclassical Transport in Quark-Gluon Plasma: Factorization and Quantum Correction

abstract: We study quarkonium transport in the quark-gluon plasma by using the potential nonrelativistic QCD (pNRQCD) effective field theory and the framework of open quantum systems. We argue that the coupling between quarkonium and the thermal bath is weak using separation of scales, so the initial density matrix of the total system factorizes and the time evolution of the subsystem is Markovian. We derive the semiclassical Boltzmann equation for quarkonium by applying a Wigner transform to the Lindblad equation and carrying out a semiclassical expansion. We resum relevant interactions to all orders in the coupling constant at leading power of the nonrelativistic and multipole expansions. The derivation is valid for both weakly coupled and strongly coupled quark-gluon plasmas. We find reaction rates in the transport equation factorize into a quarkonium dipole transition function and a chromoelectric gluon distribution function. For the differential reaction rate, the definition of the momentum dependent chromoelectric gluon distribution function involves staple-shaped Wilson lines. For the inclusive reaction rate, the Wilson lines collapse into a straight line along the real time axis and the distribution becomes momentum independent. The relation between the two Wilson lines is analogous to the relation between the Wilson lines appearing in the gluon parton distribution function (PDF) and the gluon transverse momentum dependent parton distribution function (TMDPDF). The centrality dependence of the quarkonium nuclear modification factor measured by experiments probes the momentum independent distribution while the transverse momentum dependence and measurements of the azimuthal angular anisotropy may be able to probe the momentum dependent one. We discuss one way to indirectly constrain the quarkonium in-medium real potential by using the factorization formula and lattice calculations. The leading quantum correction to the semiclassical transport equation of quarkonium is also worked out. The study can be easily generalized to quarkonium transport in cold nuclear matter, which is relevant for quarkonium production in eA collisions in the future Electron-Ion Collider.


title: Spectroscopic Detection of the Pre-White Dwarf Companion of Regulus

abstract: Mass transfer in an interacting binary will often strip the mass donor of its entire envelope and spin up the mass gainer to near critical rotation. The nearby B-type star Regulus represents a binary in the post-mass transfer stage: it is a rapid rotator with a very faint companion in a 40 d orbit. Here we present the results of a search for the spectral features of the stripped-down star in an extensive set of high S/N and high resolution spectra obtained with the CFHT/ESPaDOnS and TBL/NARVAL spectrographs. We first determine revised orbital elements in order to set accurate estimates of the orbital Doppler shifts at the times of observation. We then calculate cross-correlation functions of the observed and model spectra, and we search for evidence of the companion signal in the residuals after removal of the strong primary component. We detect a weak peak in the co-added residuals that has the properties expected for a faint pre-white dwarf. We use the dependence of the peak height and width on assumed secondary velocity semiamplitude to derive the semiamplitude, which yields masses of $M_1/M_\odot = 3.7 \pm 1.4$ and $M_2/M_\odot = 0.31 \pm 0.10$ (assuming orbital inclination equals the spin inclination of Regulus). We estimate the pre-white dwarf temperature $T_{\rm eff} = (20 \pm 4)$~kK through tests with differing temperature model spectra, and we find the radius $R_2/R_\odot = 0.061 \pm 0.011$ from the component temperatures and the flux ratio associated with the amplitude of the signal in the cross-correlation residuals.


title: Constraints on jet quenching from a multi-stage energy-loss approach

abstract: We present a multi-stage model for jet evolution through a quark-gluon plasma within the JETSCAPE framework. The multi-stage approach in JETSCAPE provides a unified description of distinct phases in jet shower contingent on the virtuality. We demonstrate a simultaneous description of leading hadron and integrated jet observables as well as jet $v_n$ using tuned parameters. Medium response to the jet quenching is implemented based on a weakly-coupled recoil prescription. We also explore the cone-size dependence of jet energy loss inside the plasma.


title: Prediction of aerothermal characteristics of a generic hypersonic inlet flow

abstract: The accurate prediction of aerothermal surface loading is of paramount importance for the design of high speed flight vehicles. In this work, we consider the numerical solution of hypersonic flow over a double-finned geometry, representative of the inlet of an air-breathing flight vehicle, characterized by three-dimensional intersecting shock-wave/turbulent boundary-layer interaction at Mach 8.3. High Reynolds numbers ($Re_L \approx 11.6 \times 10^6$ based on free-stream conditions) and the presence of cold walls ($T_w/T_o \approx 0.27$) leading to large near-wall temperature gradients necessitate the use of wall-modeled large-eddy simulation (WMLES) in order to make calculations computationally tractable. The comparison of the WMLES results with experimental measurements shows good agreement in the time-averaged surface heat flux and wall pressure distributions, and the WMLES predictions show reduced errors with respect to the experimental measurements than prior RANS calculations. The favorable comparisons are obtained using an LES wall model based on equilibrium boundary layer approximations despite the presence of numerous non-equilibrium conditions including three dimensionality, shock-boundary layer interactions, and flow separation. Lastly, it is also demonstrated that the use of semi-local eddy viscosity scaling (in lieu of the commonly used van Driest scaling) in the LES wall model is necessary to accurately predict the surface pressure loading and heat fluxes.


title: 2.5D Root of Trust: Secure System-Level Integration of Untrusted Chiplets

abstract: Dedicated, after acceptance and publication, in memory of the late Vassos Soteriou. For the first time, we leverage the 2.5D interposer technology to establish system-level security in the face of hardware- and software-centric adversaries. More specifically, we integrate chiplets (i.e., third-party hard intellectual property of complex functionality, like microprocessors) using a security-enforcing interposer. Such hardware organization provides a robust 2.5D root of trust for trustworthy, yet powerful and flexible, computation systems. The security paradigms for our scheme, employed firmly by design and construction, are: 1) stringent physical separation of trusted from untrusted components, and 2) runtime monitoring. The system-level activities of all untrusted commodity chiplets are checked continuously against security policies via physically separated security features. Aside from the security promises, the good economics of outsourced supply chains are still maintained; the system vendor is free to procure chiplets from the open market, while only producing the interposer and assembling the 2.5D system oneself. We showcase our scheme using the Cortex-M0 core and the AHB-Lite bus by ARM, building a secure 64-core system with shared memories. We evaluate our scheme through hardware simulation, considering different threat scenarios. Finally, we devise a physical-design flow for 2.5D systems, based on commercial-grade design tools, to demonstrate and evaluate our 2.5D root of trust.


title: Sensitivity estimation for calculated phase equilibria

abstract: The development of a consistent framework for Calphad model sensitivity is necessary for the rational reduction of uncertainty via new models and experiments. In the present work, a sensitivity theory for Calphad was developed, and a closed-form expression for the log-likelihood gradient and Hessian of a multi-phase equilibrium measurement was presented. The inherent locality of the defined sensitivity metric was mitigated through the use of Monte Carlo averaging. A case study of the Cr-Ni system was used to demonstrate visualizations and analyses enabled by the developed theory. Criteria based on the classical Cram\'er-Rao bound were shown to be a useful diagnostic in assessing the accuracy of parameter covariance estimates from Markov Chain Monte Carlo. The developed sensitivity framework was applied to estimate the statistical value of phase equilibria measurements in comparison with thermochemical measurements, with implications for Calphad model uncertainty reduction.


title: On the restrictiveness of the hazard rate order

abstract: Every element $\theta=(\theta_1,\ldots,\theta_n)$ of the probability $n$-simplex induces a probability distribution $P_\theta$ of a random variable $X$ that can assume only a finite number of real values $x_1 < \cdots < x_n$ by defining $P_\theta(X=x_i) = \theta_i, 1\leq i \leq n$. We show that if $\Theta$ and $\Theta'$ are two random vectors uniformly distributed on $\Delta^n$, then $P(P_\Theta\leq_{\rm hr} P_{\Theta'})=\frac{1}{2^{n-1}}$ where $\leq_{\rm hr}$ denotes the hazard rate order.


title: Discrete surfaces with length and area and minimal fillings of the circle

abstract: We propose to imagine that every Riemannian metric on a surface is discrete at the small scale, made of curves called walls. The length of a curve is its number of wall crossings, and the area of the surface is the number of crossings of the walls themselves. We show how to approximate a Riemannian (or self-reverse Finsler) metric by a wallsystem. This work is motivated by Gromov's filling area conjecture (FAC) that the hemisphere minimizes area among orientable Riemannian surfaces that fill a circle isometrically. We introduce a discrete FAC: every square-celled surface that fills isometrically a $2n$-cycle graph has at least $n(n-1)/2$ squares. We prove that our discrete FAC is equivalent to the FAC for surfaces with self-reverse metric. If the surface is a disk, the discrete FAC follows from Steinitz's algorithm for transforming curves into pseudolines. This gives a new proof of the FAC for disks with self-reverse metric. We also imitate Ivanov's proof of the same fact, using discrete differential forms. And we prove that the FAC holds for M\"obius bands with self-reverse metric. For this we use a combinatorial curve shortening flow developed by de Graaf--Schrijver and Hass--Scott. With the same method we prove the systolic inequality for Klein bottles with self-reverse metric, conjectured by Sabourau--Yassine. Self-reverse metrics can be discretized using walls because every normed plane satisfies Crofton's formula: the length of every segment equals the symplectic measure of the set of lines that it crosses. Directed 2-dimensional metrics have no Crofton formula, but can be discretized as well. Their discretization is a triangulation where the length of each edge is 1 in one way and 0 in the other, and the area of the surface is the number of triangles. This structure is a simplicial set, dual to a plabic graph. The role of the walls is played by Postnikov's strands.


title: Relative Tur\'{a}n Problems for Uniform Hypergraphs

abstract: For two graphs $F$ and $H$, the relative Tur\'{a}n number $\mathrm{ex}(H,F)$ is the maximum number of edges in an $F$-free subgraph of $H$. Foucaud, Krivelevich, and Perarnau \cite{FKP} and Perarnau and Reed \cite{PR} studied these quantities as a function of the maximum degree of $H$. In this paper, we study a generalization for uniform hypergraphs. If $F$ is a complete $r$-partite $r$-uniform hypergraph with parts of sizes $s_1,s_2,\dots,s_r$ with each $s_{i + 1}$ sufficiently large relative to $s_i$, then with $1/\beta = \sum_{i = 2}^r \prod_{j = 1}^{i - 1} s_j$ we prove that for any $r$-uniform hypergraph $H$ with maximum degree $\Delta$, \[\mathrm{ex}(H,F)\ge \Delta^{-\beta - o(1)} \cdot e(H).\] This is tight as $\Delta \rightarrow \infty$ up to the $o(1)$ term in the exponent, since we show there exists a $\Delta$-regular $r$-graph $H$ such that $\mathrm{ex}(H,F)=O(\Delta^{-\beta}) \cdot e(H)$. Similar tight results are obtained when $H$ is the random $n$-vertex $r$-graph $H_{n,p}^r$ with edge-probability $p$, extending results of Balogh and Samotij \cite{BS} and Morris and Saxton \cite{MS}.


title: Interplay Between Resilience and Accuracy in Resilient Vector Consensus in Multi-Agent Networks

abstract: In this paper, we study the relationship between resilience and accuracy in the resilient distributed multi-dimensional consensus problem. We consider a network of agents, each of which has a state in $\mathbb{R}^d$. Some agents in the network are adversarial and can change their states arbitrarily. The normal (non-adversarial) agents interact locally and update their states to achieve consensus at some point in the convex hull $\calC$ of their initial states. This objective is achievable if the number of adversaries in the neighborhood of normal agents is less than a specific value, which is a function of the local connectivity and the state dimension $d$. However, to be resilient against adversaries, especially in the case of large $d$, the desired local connectivity is large. We discuss that resilience against adversarial agents can be improved if normal agents are allowed to converge in a bounded region $\calB\supseteq\calC$, which means normal agents converge at some point close to but not necessarily inside $\calC$ in the worst case. The accuracy of resilient consensus can be measured by the Hausdorff distance between $\calB$ and $\calC$. As a result, resilience can be improved at the cost of accuracy. We propose a resilient bounded consensus algorithm that exploits the trade-off between resilience and accuracy by projecting $d$-dimensional states into lower dimensions and then solving instances of resilient consensus in lower dimensions. We analyze the algorithm, present various resilience and accuracy bounds, and also numerically evaluate our results.


title: Explanation of Unintended Radiated Emission Classification via LIME

abstract: Unintended radiated emissions arise during the use of electronic devices. Identifying and mitigating the effects of these emissions is a key element of modern power engineering and associated control systems. Signal processing of the electrical system can identify the sources of these emissions. A dataset known as Flaming Moes includes captured unintended radiated emissions from consumer electronics. This dataset was analyzed to construct next-generation methods for device identification. To this end, a neural network based on applying the ResNet-18 image classification architecture to the short time Fourier transforms of short segments of voltage signatures was constructed. Using this classifier, the 18 device classes and background class were identified with close to 100 percent accuracy. By applying LIME to this classifier and aggregating the results over many classifications for the same device, it was possible to determine the frequency bands used by the classifier to make decisions. Using ensembles of classifiers trained on very similar datasets from the same parent data distribution, it was possible to recover robust sets of features of device output useful for identification. The additional understanding provided by the application of LIME enhances the trainability, trustability, and transferability of URE analysis networks.


title: Assess the impacts of human mobility change on COVID-19 dynamics in Arizona, U.S.: a modeling study incorporating Google Community Mobility Reports

abstract: In June 2020, Arizona, U.S., emerged as one of the world's worst coronavirus disease 2019(COVID-19) spots after the stay-at-home order was lifted in the middle of May. However, with the decisions to reimpose restrictions, the number of COVID-19 cases has been declining, and Arizona is considered to be a good model in slowing the epidemic. In this paper, we aimed to examine the COVID-19 situation in Arizona and assess the impact of human mobility change. We constructed the mobility integrated metapopulation susceptible-infectious-removed model and fitted to publicly available datasets on COVID-19 cases and mobility changes in Arizona. Our simulations showed that by reducing human mobility, the peak time was delayed, and the final size of the epidemic was decreased in all three regions. Our analysis suggests that rapid and effective decision making is crucial to control human mobility and, therefore, COVID-19 epidemics. Until a vaccine is available, reimplementations of mobility restrictions in response to the increase of new COVID-19 cases might need to be considered in Arizona and beyond.


title: New Families of Solutions for the Space Time Fractional Burgers Equation

abstract: In this paper, the hyperbolic tangent function method is applied for constructing exact solutions for space-time conformal fractional Burgers equation. Furthermore, the space-time conformal fractional Burgers equation is tested for the Painleve property and consequently new numerous exact solutions are generated via Backlund transform.


title: Evidence for top quark production in nucleus-nucleus collisions

abstract: Using $1.7 \pm 0.1\,\mathrm{nb^{-1}}$ of lead-lead ($A = 208$) collision data recorded by the CMS experiment at a nucleon-nucleon center-of-mass energy of 5.02 TeV, we report evidence of top quark pair ($\mathrm{t\bar{t}}$) production. The $\mathrm{t\bar{t}}$ cross section ($\sigma_\mathrm{t\bar{t}}$) is extracted from likelihood fits to a multivariate discriminator using lepton kinematic variables in dilepton final states and two methods. One method relies on the leptonic information alone, and the second one exploits, in addition, the presence of bottom quarks. The measured $\sigma_\mathrm{t\bar{t}}$ is $2.54^{+0.84}_{-0.74}$ and $2.03^{+0.71}_{-0.64}\, \mu\mathrm{b}$ in the two cases, respectively, consistent with predictions from perturbative quantum chromodynamics. We demonstrate, for the first time, that top quark decay products (leptonically decaying W bosons and bottom quarks) can be identified, irrespective of any possible final-state interactions with the quark-gluon plasma.


title: Unicritical polynomial maps with rational multipliers

abstract: In this article, we prove that every unicritical polynomial map that has only rational multipliers is either a power map or a Chebyshev map. This provides some evidence in support of a conjecture by Milnor concerning rational maps whose multipliers are all integers.


title: A General Framework for Fairness in Multistakeholder Recommendations

abstract: Contemporary recommender systems act as intermediaries on multi-sided platforms serving high utility recommendations from sellers to buyers. Such systems attempt to balance the objectives of multiple stakeholders including sellers, buyers, and the platform itself. The difficulty in providing recommendations that maximize the utility for a buyer, while simultaneously representing all the sellers on the platform has lead to many interesting research problems.Traditionally, they have been formulated as integer linear programs which compute recommendations for all the buyers together in an \emph{offline} fashion, by incorporating coverage constraints so that the individual sellers are proportionally represented across all the recommended items. Such approaches can lead to unforeseen biases wherein certain buyers consistently receive low utility recommendations in order to meet the global seller coverage constraints. To remedy this situation, we propose a general formulation that incorporates seller coverage objectives alongside individual buyer objectives in a real-time personalized recommender system. In addition, we leverage highly scalable submodular optimization algorithms to provide recommendations to each buyer with provable theoretical quality bounds. Furthermore, we empirically evaluate the efficacy of our approach using data from an online real-estate marketplace.


title: The Effects of Primordial Black Holes on Dark Matter Models

abstract: We investigate the effects of producing dark matter by Hawking evaporation of primordial black holes (PBHs) in scenarios that may have a second well-motivated dark matter production mechanism, such as freeze-out, freeze-in, or gravitational production. We show that the interplay between PBHs and the alternative sources of dark matter can give rise to model-independent modifications to the required dark matter abundance from each production mechanism, which in turn affect the prospects for dark matter detection. In particular, we demonstrate that for the freeze-out mechanism, accounting for evaporation of PBHs after freeze-out demands a larger annihilation cross section of dark matter particles than its canonical value for a thermal dark matter. For mechanisms lacking thermalization due to a feeble coupling to the thermal bath, we show that the PBH contribution to the dark matter abundance leads to the requirement of an even feebler coupling. Moreover, we show that when a large initial abundance of PBHs causes an early matter-dominated epoch, PBH evaporation alone cannot explain the whole abundance of dark matter today. In this case, an additional production mechanism is required, in contrast to the case when PBHs are formed and evaporate during a radiation-dominated epoch.


title: Chiral topological superconductivity in Josephson junction

abstract: We consider a heterostructure of semiconductor layers sandwiched between two superconductors, forming a two-dimensional Josephson junction. Applying a Zeeman field perpendicular to the junction can render a topological superconducting phase with chiral Majorana edge mode. We show that the phase difference between two superconductors can efficiently reduce the magnetic field required to achieve the chiral topological superconductivity, providing an experimentally feasible setup to realize chiral Majorana edge modes. We also construct a lattice Hamiltonian of the setup to demonstrate the chiral Majorana edge mode and the Majorana bound state localized in vortices.


title: How fast is a quantum jump?

abstract: A proposal is put forward for an estimate of the duration of a transition between atomic states. The proposal rests on the consideration that a resonance of the atomic electron with modes of the zero-point radiation field of Compton's frequency is at the core of the phenomenon. The theoretical result, given essentially by the expression (\alpha\omega_{C})^{-1}, where \alpha is the fine structure constant and \omega_{C} the Compton angular frequency for the electron, lies well within the range of the recently experimentally estimated values of the order of attoseconds (10^{-18} s).


title: A Polariton-Stabilized Spin Clock

abstract: Atom-like quantum systems in solids have been proposed as a compact alternative for atomic clocks, but realizing the potential of solid-state technology will requires an architecture design which overcomes traditional limitations such as magnetic and temperature-induced systematics. Here, we propose a solution to this problem: a `solid-state spin clock' that hybridizes a microwave resonator with a magnetic-field-insensitive spin transition within the ground state of the diamond nitrogen-vacancy center. Detailed numerical and analytical modeling of this `polariton-stabilized' spin clock (PSSC) indicates a potential fractional frequency instability below $10^{-13}$ at 1 second measurement time, assuming present-day experimental parameters. This stability would represent a significant improvement over the state-of-the-art in miniaturized atomic vapor clocks.


title: An equivalence between enriched $\infty$-categories and $\infty$-categories with weak action

abstract: We show that an $\infty$-category $\mathcal{M}$ with a closed left action of a monoidal $\infty$-category $\mathcal{V}$ is completely determined by the $\mathcal{V}$-valued graph of morphism objects equipped with the structure of a $\mathcal{V}$-enrichment in the sense of Gepner-Haugseng. We prove a similar result when $\mathcal{M}$ is a $\mathcal{V}$-enriched $\infty$-category in the sense of Lurie, an operadic generalization of the notion of $\infty$-category with closed left action. Precisely, we prove that sending a $\mathcal{V}$-enriched $\infty$-category in the sense of Lurie to the $\mathcal{V}$-valued graph of morphism objects refines to an equivalence $\chi$ between the $\infty$-category of $\mathcal{V}$-enriched $\infty$-categories in the sense of Lurie and of Gepner-Haugseng. Moreover if $\mathcal{V}$ is a presentably $\mathbb{E}_{\mathrm{k+1}}$-monoidal $\infty$-category for $1 \leq k \leq \infty$, we prove that $\chi$ restricts to a lax $\mathbb{E}_{\mathrm{k}}$-monoidal functor between the $\infty$-category of left $\mathcal{V}$-modules in $\mathrm{Pr}^L$, the symmetric monoidal $\infty$-category of presentable $\infty$-categories, endowed with the relative tensor product, and the tensor product of $\mathcal{V}$-enriched $\infty$-categories of Gepner-Haugseng. As an application of our theory we construct a lax symmetric monoidal embedding of the $\infty$-category of small stable $\infty$-categories into the $\infty$-category of small spectral $\infty$-categories. As a second application we produce a Yoneda-embedding for Lurie's notion of enriched $\infty$-categories.


title: Player Identification in Hockey Broadcast Videos

abstract: We present a deep recurrent convolutional neural network (CNN) approach to solve the problem of hockey player identification in NHL broadcast videos. Player identification is a difficult computer vision problem mainly because of the players' similar appearance, occlusion, and blurry facial and physical features. However, we can observe players' jersey numbers over time by processing variable length image sequences of players (aka 'tracklets'). We propose an end-to-end trainable ResNet+LSTM network, with a residual network (ResNet) base and a long short-term memory (LSTM) layer, to discover spatio-temporal features of jersey numbers over time and learn long-term dependencies. For this work, we created a new hockey player tracklet dataset that contains sequences of hockey player bounding boxes. Additionally, we employ a secondary 1-dimensional convolutional neural network classifier as a late score-level fusion method to classify the output of the ResNet+LSTM network. This achieves an overall player identification accuracy score over 87% on the test split of our new dataset.


title: Towards the Development of Entropy-Based Anomaly Detection in an Astrophysics Simulation

abstract: The use of AI and ML for scientific applications is currently a very exciting and dynamic field. Much of this excitement for HPC has focused on ML applications whose analysis and classification generate very large numbers of flops. Others seek to replace scientific simulations with data-driven surrogate models. But another important use case lies in the combination application of ML to improve simulation accuracy. To that end, we present an anomaly problem which arises from a core-collapse supernovae simulation. We discuss strategies and early successes in applying anomaly detection techniques from machine learning to this scientific simulation, as well as current challenges and future possibilities.


title: Accenture at CheckThat! 2020: If you say so: Post-hoc fact-checking of claims using transformer-based models

abstract: We introduce the strategies used by the Accenture Team for the CLEF2020 CheckThat! Lab, Task 1, on English and Arabic. This shared task evaluated whether a claim in social media text should be professionally fact checked. To a journalist, a statement presented as fact, which would be of interest to a large audience, requires professional fact-checking before dissemination. We utilized BERT and RoBERTa models to identify claims in social media text a professional fact-checker should review, and rank these in priority order for the fact-checker. For the English challenge, we fine-tuned a RoBERTa model and added an extra mean pooling layer and a dropout layer to enhance generalizability to unseen text. For the Arabic task, we fine-tuned Arabic-language BERT models and demonstrate the use of back-translation to amplify the minority class and balance the dataset. The work presented here was scored 1st place in the English track, and 1st, 2nd, 3rd, and 4th place in the Arabic track.


title: BP-RRT: Barrier Pair Synthesis for Temporal Logic Motion Planning

abstract: For a nonlinear system (e.g. a robot) with its continuous state space trajectories constrained by a linear temporal logic specification, the synthesis of a low-level controller for mission execution often results in a non-convex optimization problem. We devise a new algorithm to solve this type of non-convex problems by formulating a rapidly-exploring random tree of barrier pairs, with each barrier pair composed of a quadratic barrier function and a full state feedback controller. The proposed method employs a rapid-exploring random tree to deal with the non-convex constraints and uses barrier pairs to fulfill the local convex constraints. As such, the method solves control problems fulfilling the required transitions of an automaton in order to satisfy given linear temporal logic constraints. At the same time it synthesizes locally optimal controllers in order to transition between the regions corresponding to the alphabet of the automaton. We demonstrate this new algorithm on a simulation of a two linkage manipulator robot.


title: Quantitative weighted estimates for rough singular integrals on homogeneous groups

abstract: In this paper, we study weighted $L^{p}(w)$ boundedness ($1<p<\infty$ and $w$ a Muckenhoupt $A_{p}$ weight) of singular integrals with homogeneous convolution kernel $K(x)$ on an arbitrary homogeneous group $\mathbb H$ of dimension $\mathbb{Q}$, {under the assumption that $K_0$, the restriction of $K$ to the unit annulus, is mean zero and $L^{q}$ integrable for some $q_{0}<q\leq \infty$,} where $q_{0}$ is a fixed constant depending on $w$. We obtain a quantitative weighted bound, which is consistent with the one obtained by Hyt\"onen--Roncal--Tapiola in the Euclidean setting, for this operator on $L^{p}(w)$. Comparing to the previous results in the Euclidean setting, our assumptions on the kernel and on the underlying space are weaker. Moreover, we investigate the quantitative weighted bound for the bi-parameter rough singular integrals on product homogeneous Lie groups.


title: Galactic seismology: the evolving "phase spiral" after the Sagittarius dwarf impact

abstract: In 2018, the ESA \Gaia\ satellite discovered a remarkable spiral pattern ("phase spiral") in the $z-V_z$ phase plane throughout the solar neighbourhood, where $z$ and $V_z$ are the displacement and velocity of a star perpendicular to the Galactic disc. In response to Binney \& Sch\"onrich's analytic model of a disc-crossing satellite to explain the \Gaia\ data, we carry out a high-resolution, N-body simulation (N$\:\approx 10^8$ particles) of an impulsive mass ($2\times 10^{10}$ \Msun) that interacts with a cold stellar disc at a single transit point. The disc response is complex since the impulse triggers a superposition of two distinct bisymmetric ($m=2$) modes $-$ a density wave and a corrugated bending wave $-$ that wrap up at different rates. Stars in the {\it faster} density wave wrap up with time $T$ according to $\phi_D(R,T)=(\Omega_D(R) + \Omega_{\rm o})\:T$ where $\phi_D$ describes the spiral pattern and $\Omega_D =\Omega(R) -\kappa(R)/2$, where $\kappa$ is the epicyclic frequency. While the pattern speed $\Omega_{\rm o}$ is small, it is non-zero. The {\it slower} bending wave wraps up according to $\Omega_B\approx\Omega_D/2$ producing a corrugated wave. The bunching effect of the density wave triggers the phase spiral as it rolls up and down on the bending wave ("rollercoaster" model). The phase spiral emerges slowly about $\Delta T \approx 400$ Myr after impact. It appears to be a long-lived, disc-wide phenomenon that continues to evolve over most of the 2~Gyr simulation. Thus, given Sagittarius' (Sgr) low total mass today ($M_{\rm tot}\sim 3\times 10^8$ \Msun\ within 10 kpc diameter), we believe the phase spiral was excited by the disc-crossing dwarf some $1-2$ Gyr {\it before} the recent transit. For this to be true, Sgr must be losing mass at 0.5-1 dex per orbit loop.


title: Determination of the scaling characteristics of time-dependent optical properties of microalgae using electromagnetic scattering

abstract: The time-dependent optical properties of microalgae are crucial for light transfer in photobioreactor (PBR) designs. Here, the time-dependent optical properties were derived using electromagnetic scattering theory informed by the experimentally measured optical properties of Chlorella protothecoides. The temporal scaling functions (TSFs) of nearly wavelength-independent absorption and scattering cross-sections were demonstrated using the electromagnetic scattering theory, leading to the first concrete expression of the TSF. The TSF establishes the relationship between the time-dependent absorption/scattering cross-sections at the stationery and growth phases of microalgal development. The concrete expression of the TSF provides a new way to calculate the time-dependent optical properties of microalgae using electromagnetic scattering theory. The areas where the TSF of microalgae has exciting potential are in remote sensing and photobioreactor applications.


title: Communication-efficient distributed eigenspace estimation

abstract: Distributed computing is a standard way to scale up machine learning and data science algorithms to process large amounts of data. In such settings, avoiding communication amongst machines is paramount for achieving high performance. Rather than distribute the computation of existing algorithms, a common practice for avoiding communication is to compute local solutions or parameter estimates on each machine and then combine the results; in many convex optimization problems, even simple averaging of local solutions can work well. However, these schemes do not work when the local solutions are not unique. Spectral methods are a collection of such problems, where solutions are orthonormal bases of the leading invariant subspace of an associated data matrix, which are only unique up to rotation and reflections. Here, we develop a communication-efficient distributed algorithm for computing the leading invariant subspace of a data matrix. Our algorithm uses a novel alignment scheme that minimizes the Procrustean distance between local solutions and a reference solution, and only requires a single round of communication. For the important case of principal component analysis (PCA), we show that our algorithm achieves a similar error rate to that of a centralized estimator. We present numerical experiments demonstrating the efficacy of our proposed algorithm for distributed PCA, as well as other problems where solutions exhibit rotational symmetry, such as node embeddings for graph data and spectral initialization for quadratic sensing.


title: GazeMAE: General Representations of Eye Movements using a Micro-Macro Autoencoder

abstract: Eye movements are intricate and dynamic events that contain a wealth of information about the subject and the stimuli. We propose an abstract representation of eye movements that preserve the important nuances in gaze behavior while being stimuli-agnostic. We consider eye movements as raw position and velocity signals and train separate deep temporal convolutional autoencoders. The autoencoders learn micro-scale and macro-scale representations that correspond to the fast and slow features of eye movements. We evaluate the joint representations with a linear classifier fitted on various classification tasks. Our work accurately discriminates between gender and age groups, and outperforms previous works on biometrics and stimuli clasification. Further experiments highlight the validity and generalizability of this method, bringing eye tracking research closer to real-world applications.


title: A Large Scale Empirical Study of the Impact of Spaghetti Code and Blob Anti-patterns on Program Comprehension

abstract: Context: Several studies investigated the impact of anti-patterns (i.e., "poor" solutions to recurring design problems) during maintenance activities and reported that anti-patterns significantly affect the developers' effort required to edit files. However, before developers edit files, they must understand the source code of the systems. This source code must be easy to understand by developers. Objective: In this work, we provide a complete assessment of the impact of two instances of two anti-patterns, Blob or Spaghetti Code, on program comprehension. Method: We analyze the impact of these two anti-patterns through three empirical studies conducted at Polytechnique Montr\'eal (Canada) with 24 participants; at Carlton University (Canada) with 30 participants; and at University Basilicata (Italy) with 79 participants. Results: We collect data from 372 tasks obtained thanks to 133 different participants from the three universities. We use three metrics to assess the developers' comprehension of the source code: (1) the duration to complete each task; (2) their percentage of correct answers; and, (3) the NASA task load index for their effort. Conclusions: We report that, although single occurrences of Blob or Spaghetti code anti-patterns have little effect on code comprehension, two occurrences of either Blob or Spaghetti Code significantly increases the developers' time spent in their tasks, reduce their percentage of correct answers, and increase their effort. Hence, we recommend that developers act on both anti-patterns, which should be refactored out of the source code whenever possible. We also recommend further studies on combinations of anti-patterns rather than on single anti-patterns one at a time.


title: Optimizing Mode Connectivity via Neuron Alignment

abstract: The loss landscapes of deep neural networks are not well understood due to their high nonconvexity. Empirically, the local minima of these loss functions can be connected by a learned curve in model space, along which the loss remains nearly constant; a feature known as mode connectivity. Yet, current curve finding algorithms do not consider the influence of symmetry in the loss surface created by model weight permutations. We propose a more general framework to investigate the effect of symmetry on landscape connectivity by accounting for the weight permutations of the networks being connected. To approximate the optimal permutation, we introduce an inexpensive heuristic referred to as neuron alignment. Neuron alignment promotes similarity between the distribution of intermediate activations of models along the curve. We provide theoretical analysis establishing the benefit of alignment to mode connectivity based on this simple heuristic. We empirically verify that the permutation given by alignment is locally optimal via a proximal alternating minimization scheme. Empirically, optimizing the weight permutation is critical for efficiently learning a simple, planar, low-loss curve between networks that successfully generalizes. Our alignment method can significantly alleviate the recently identified robust loss barrier on the path connecting two adversarial robust models and find more robust and accurate models on the path.


title: Game Industry Problems: an Extensive Analysis of the Gray Literature

abstract: Context: Given its competitiveness, the video-game industry has a closed-source culture. Hence, little is known of the problems faced by game developers. However, game developers do share information about their games projects through postmortems, which describe informally what happened during the projects. Objective: The software-engineering research community and game developers would benefit from a state of the problems of the video game industry, in particular the problems faced by game developers, their evolution in time, and their root causes. This state of the practice would allow researchers and practitioners to work towards solving these problems. Method: We analyzed 200 postmortems from 1997 to 2019, resulting in 927 problems divided into 20 types. Through our analysis, we described the overall landscape of game industry problems in the past 23 years and how these problems evolved over the years. We also give details on the most common problems, their root causes, and possible solutions. We finally discuss suggestions for future projects. Results: We observe that (1) the game industry suffers from management and production problems in the same proportion; (2) management problems decreased over the years giving space to business problems, while production problems remained constant; (3a) technical and game design problems are decreasing over the years, the latter only after the last decade; (3b) problems related to the team increase over the last decade;(3c) marketing problems are the ones that had the biggest increase over the 23 years compared to other problem types; (4) finally, the majority of the main root causes are related to people, not technologies. Conclusions: In this paper we provide a state of the practice for researchers to understand and study video-game development problems. We also offer suggestions to help practitioners to avoid the most common problems.


title: Polyphorm: Structural Analysis of Cosmological Datasets via Interactive Physarum Polycephalum Visualization

abstract: This paper introduces Polyphorm, an interactive visualization and model fitting tool that provides a novel approach for investigating cosmological datasets. Through a fast computational simulation method inspired by the behavior of Physarum polycephalum, an unicellular slime mold organism that efficiently forages for nutrients, astrophysicists are able to extrapolate from sparse datasets, such as galaxy maps archived in the Sloan Digital Sky Survey, and then use these extrapolations to inform analyses of a wide range of other data, such as spectroscopic observations captured by the Hubble Space Telescope. Researchers can interactively update the simulation by adjusting model parameters, and then investigate the resulting visual output to form hypotheses about the data. We describe details of Polyphorm's simulation model and its interaction and visualization modalities, and we evaluate Polyphorm through three scientific use cases that demonstrate the effectiveness of our approach.


title: Monogenic pure cubics

abstract: Let $k\geq 2$ be a square-free integer. We prove that the number of square-free integers $m\in [1,N]$ such that $(k,m)=1$ and $\mathbb{Q}(\sqrt[3]{k^2m})$ is monogenic is $\gg N^{1/3}$ and $\ll N/(\log N)^{1/3-\epsilon}$ for any $\epsilon>0$. Assuming ABC, the upper bound can be improved to $O(N^{(1/3)+\epsilon})$. Let $F$ be the finite field of order $q$ with $(q,3)=1$ and let $g(t)\in F[t]$ be non-constant square-free. We prove unconditionally the analogous result that the number of square-free $h(t)\in F[t]$ such that $\deg(h)\leq N$, $(g,h)=1$ and $F(t,\sqrt[3]{g^2h})$ is monogenic is $\gg q^{N/3}$ and $\ll N^2q^{N/3}$.


title: Effect of transcription reinitiation in stochastic gene expression

abstract: Gene expression (GE) is an inherently random or stochastic or noisy process. The randomness in different steps of GE, e.g., transcription, translation, degradation, etc., leading to cell-to-cell variations in mRNA and protein levels. This variation appears in organisms ranging from microbes to metazoans. Stochastic gene expression has important consequences for cellular function. The random fluctuations in protein levels produce variability in cellular behavior. It is beneficial in some contexts and harmful to others. These situations include stress response, metabolism, development, cell cycle, circadian rhythms, and aging. Different model studies e.g., constitutive, two-state, etc., reveal that the fluctuations in mRNA and protein levels arise from different steps of gene expression among which the steps in transcription have the maximum effect. The pulsatile mRNA production through RNAP-II based reinitiation of transcription is an important part of gene transcription. Though, the effect of that process on mRNA and protein levels is very little known. The addition of any biochemical step in the constitutive or two-state process generally decreases the mean and increases the Fano factor. In this study, we have shown that the RNAP-II based reinitiation process in gene transcription can have different effects on both mean and Fano factor at mRNA levels in different model systems. It decreases the mean and Fano factor both at the mRNA levels in the constitutive network whereas in other networks it can simultaneously increase or decrease both quantities or it can have mixed-effect at mRNA levels. We propose that a constitutive network with reinitiation behaves like a product independent negative feedback circuit whereas other networks behave as either product independent positive or negative or mixed feedback circuit.


title: Cross-domain Adaptation with Discrepancy Minimization for Text-independent Forensic Speaker Verification

abstract: Forensic audio analysis for speaker verification offers unique challenges due to location/scenario uncertainty and diversity mismatch between reference and naturalistic field recordings. The lack of real naturalistic forensic audio corpora with ground-truth speaker identity represents a major challenge in this field. It is also difficult to directly employ small-scale domain-specific data to train complex neural network architectures due to domain mismatch and loss in performance. Alternatively, cross-domain speaker verification for multiple acoustic environments is a challenging task which could advance research in audio forensics. In this study, we introduce a CRSS-Forensics audio dataset collected in multiple acoustic environments. We pre-train a CNN-based network using the VoxCeleb data, followed by an approach which fine-tunes part of the high-level network layers with clean speech from CRSS-Forensics. Based on this fine-tuned model, we align domain-specific distributions in the embedding space with the discrepancy loss and maximum mean discrepancy (MMD). This maintains effective performance on the clean set, while simultaneously generalizes the model to other acoustic domains. From the results, we demonstrate that diverse acoustic environments affect the speaker verification performance, and that our proposed approach of cross-domain adaptation can significantly improve the results in this scenario.


title: Learning from the past: A process recommendation system for video game projects using postmortems experiences

abstract: Context: The video game industry is a billion dollar industry that faces problems in the way games are developed. One method to address these problems is using developer aid tools, such as Recommendation Systems. These tools assist developers by generating recommendations to help them perform their tasks. Objective: This article describes a systematic approach to recommend development processes for video game projects, using postmortem knowledge extraction and a model of the context of the new project, in which "postmortems" are articles written by video game developers at the end of projects, summarizing the experience of their game development team. This approach aims to provide reflections about development processes used in the game industry as well as guidance to developers to choose the most adequate process according to the contexts they're in. Method: Our approach is divided in three separate phases: in the the first phase, we manually extracted the processes from the postmortems analysis; in the second one, we created a video game context and algorithm rules for recommendation; and finally in the third phase, we evaluated the recommended processes by using quantitative and qualitative metrics, game developers feedback, and a case study by interviewing a video game development team. Contributions: This article brings three main contributions. The first describes a database of developers' experiences extracted from postmortems in the form of development processes. The second defines the main attributes that a video game project contain, which it uses to define the contexts of the project. The third describes and evaluates a recommendation system for video game projects, which uses the contexts of the projects to identify similar projects and suggest a set of activities in the form of a process.


title: Transcendental Series of Reciprocals of Fibonacci and Lucas Numbers

abstract: Let $F_1=1,F_2=1,\ldots$ be the Fibonacci sequence. Motivated by the identity $\displaystyle\sum_{k=0}^{\infty}\frac{1}{F_{2^k}}=\frac{7-\sqrt{5}}{2}$, Erd\"os and Graham asked whether $\displaystyle\sum_{k=1}^{\infty}\frac{1}{F_{n_k}}$ is irrational for any sequence of positive integers $n_1,n_2,\ldots$ with $\frac{n_{k+1}}{n_k}\geq c>1$. We resolve the transcendence counterpart of their question: as a special case of our main theorem, we have that $\displaystyle\sum_{k=1}^{\infty}\frac{1}{F_{n_k}}$ is transcendental when $\frac{n_{k+1}}{n_k}\geq c>2$. The bound $c>2$ is best possible thanks to the identity at the beginning. This paper provides a new way to apply the Subspace Theorem to obtain transcendence results and extends previous non-trivial results obtainable by only Mahler's method for special sequences of the form $n_k=d^k+r$.


title: Oceanic non-Kolmogorov optical turbulence and spherical wave propagation

abstract: Light propagation in turbulent media is conventionally studied with the help of the spatio-temporal power spectra of the refractive index fluctuations. In particular, for natural water turbulence several models for the spatial power spectra have been developed based on the classic, Kolmogorov postulates. However, as currently widely accepted, non-Kolmogorov turbulent regime is also common in the stratified flow fields, as suggested by recent developments in atmospheric optics. Until now all the models developed for the non-Kolmogorov optical turbulence were pertinent to atmospheric research and, hence, involved only one advected scalar, e.g., temperature. We generalize the oceanic spatial power spectrum, based on two advected scalars, temperature and salinity concentration, to the non-Kolmogorov turbulence regime, with the help of the so-called "Upper-Bound Limitation" and by adopting the concept of spectral correlation of two advected scalars. The proposed power spectrum can handle general non-Kolmogorov, anisotropic turbulence but reduces to Kolmogorov, isotropic case if the power law exponents of temperature and salinity are set to 11/3 and anisotropy coefficient is set to unity. To show the application of the new spectrum, we derive the expression for the second-order mutual coherence function of a spherical wave and examine its coherence radius (in both scalar and vector forms) to characterize the turbulent disturbance. Our numerical calculations show that the statistics of the spherical wave vary substantially with temperature and salinity non-Kolmogorov power law exponents and temperature-salinity spectral correlation coefficient. The introduced spectrum is envisioned to become of significance for theoretical analysis and experimental measurements of non-classic natural water double-diffusion turbulent regimes.


title: Are the Old Days Gone? A Survey on Actual Software Engineering Processes in Video Game Industry

abstract: In the past 10 years, several researches studied video game development process who proposed approaches to improve the way how games are developed. These approaches usually adopt agile methodologies because of claims that traditional practices and the waterfall process are gone. However, are the "old days" really gone in the game industry? In this paper, we present a survey of software engineering processes in video game industry from postmortem project analyses. We analyzed 20 postmortems from Gamasutra Portal. We extracted their processes and modelled them through using the Business Process Model and Notation (BPMN). This work presents three main contributions. First, a postmortem analysis methodology to identify and extract project processes. Second, the study main result: \textbf{the "old days" are gone, but not completely}. \textbf{Iterative practices} are increasing and are applied to at least \textbf{65\% of projects} in which \textbf{45\% of this projects} explicitly adopted Agile practices. However, \textbf{waterfall} process is still applied at least \textbf{30\% of projects}. Finally, we discuss some implications, directions and opportunities for video game development community.


title: Hierarchical Roofline Analysis: How to Collect Data using Performance Tools on Intel CPUs and NVIDIA GPUs

abstract: This paper surveys a range of methods to collect necessary performance data on Intel CPUs and NVIDIA GPUs for hierarchical Roofline analysis. As of mid-2020, two vendor performance tools, Intel Advisor and NVIDIA Nsight Compute, have integrated Roofline analysis into their supported feature set. This paper fills the gap for when these tools are not available, or when users would like a more customized workflow for certain analysis. Specifically, we will discuss how to use Intel Advisor, RRZE LIKWID, Intel SDE and Intel Amplifier on Intel architectures, and nvprof, Nsight Compute metrics, and Nsight Compute section files on NVIDIA architectures. These tools will be used to collect information for as many memory/cache levels in the memory hierarchy as possible in order to provide insights into application's data reuse and cache locality characteristics.


title: Chiral Rings, Futaki Invariants, Plethystics, and Groebner Bases

abstract: We study chiral rings of 4d $\mathcal{N}=1$ supersymmetric gauge theories via the notion of K-stability. We show that when using Hilbert series to perform the computations of Futaki invariants, it is not enough to only include the test symmetry information in the former's denominator. We discuss a way to modify the numerator so that K-stability can be correctly determined, and a rescaling method is also applied to simplify the calculations involving test configurations. All of these are illustrated with a host of examples, by considering vacuum moduli spaces of various theories. Using Gr\"obner basis and plethystic techniques, many non-complete intersections can also be addressed, thus expanding the list of known theories in the literature.


title: Corrosion Resistance of Sulfur-Selenium Alloy Coatings

abstract: Despite decades of research, metallic corrosion remains a long-standing challenge in many engineering applications. Specifically, designing a material that can resist corrosion both in abiotic as well as biotic environments remains elusive. Here we design a lightweight sulfur-selenium (S-Se) alloy with high stiffness and ductility that can serve as a universal corrosion-resistant coating with protection efficiency of ~99.9% for steel in a wide range of diverse environments. S-Se coated mild steel shows a corrosion rate that is 6-7 orders of magnitude lower than bare metal in abiotic (simulated seawater and sodium sulfate solution) and biotic (sulfate-reducing bacterial medium) environments. The coating is strongly adhesive and mechanically robust. We attribute the high corrosion resistance of the alloy in diverse environments to its semi-crystalline, non-porous, anti-microbial, and viscoelastic nature with superior mechanical performance, enabling it to successfully block a variety of diffusing species.


title: A Lower Bound on Determinantal Complexity

abstract: The determinantal complexity of a polynomial $P \in \mathbb{F}[x_1, \ldots, x_n]$ over a field $\mathbb{F}$ is the dimension of the smallest matrix $M$ whose entries are affine functions in $\mathbb{F}[x_1, \ldots, x_n]$ such that $P = Det(M)$. We prove that the determinantal complexity of the polynomial $\sum_{i = 1}^n x_i^n$ is at least $1.5n - 3$. For every $n$-variate polynomial of degree $d$, the determinantal complexity is trivially at least $d$, and it is a long standing open problem to prove a lower bound which is super linear in $\max\{n,d\}$. Our result is the first lower bound for any explicit polynomial which is bigger by a constant factor than $\max\{n,d\}$, and improves upon the prior best bound of $n + 1$, proved by Alper, Bogart and Velasco [ABV17] for the same polynomial.


title: Incidence bounds on Edge Partitions of $K_n$

abstract: We solve a problem conjectured by Cheriyan, giving sharp bounds for incidence of certain edge partitions of the connected graph on $n$-vertices. We briefly discuss the history of the problem and relation to node connectivity of strongly regular graphs. We show that the bound cannot be made sharper.


title: Reflected anticipated backward stochastic differential equations with nonlinear resistance

abstract: In this paper, we consider reflected anticipated backward stochastic differential equations (RABSDEs, for short) with an additional resistance in the generators. Firstly, we study the existence and uniqueness results. In Luo (2020), the condition of a small time horizon is needed. Compared with the proving method in Luo (2020), we use a different proving method to avoid requiring the Lipschitz coefficients of generators $f(t,y,z,\theta,\vartheta,m,\bar{m})$ for $y,z,\theta,\vartheta$ to be small enough. We only require the Lipschitz coefficient for resistance in generator is small enough. Moreover, a probabilistic structure for solution is specified. Secondly, we give a comparison theorem for this type of equation. At last, under the linear growth condition and some other conditions on resistance , we derive the minimal solution.


title: User-Guided Domain Adaptation for Rapid Annotation from User Interactions: A Study on Pathological Liver Segmentation

abstract: Mask-based annotation of medical images, especially for 3D data, is a bottleneck in developing reliable machine learning models. Using minimal-labor user interactions (UIs) to guide the annotation is promising, but challenges remain on best harmonizing the mask prediction with the UIs. To address this, we propose the user-guided domain adaptation (UGDA) framework, which uses prediction-based adversarial domain adaptation (PADA) to model the combined distribution of UIs and mask predictions. The UIs are then used as anchors to guide and align the mask prediction. Importantly, UGDA can both learn from unlabelled data and also model the high-level semantic meaning behind different UIs. We test UGDA on annotating pathological livers using a clinically comprehensive dataset of 927 patient studies. Using only extreme-point UIs, we achieve a mean (worst-case) performance of 96.1%(94.9%), compared to 93.0% (87.0%) for deep extreme points (DEXTR). Furthermore, we also show UGDA can retain this state-of-the-art performance even when only seeing a fraction of available UIs, demonstrating an ability for robust and reliable UI-guided segmentation with extremely minimal labor demands.


title: Sparse domination and weighted estimates for rough bilinear singular integrals

abstract: Let $r>\frac{4}{3}$ and let $\Omega \in L^{r}(\mathbb{S}^{2n-1})$ have vanishing integral. We show that the bilinear rough singular integral $$T_{\Omega}(f,g)(x)= \textrm{p.v.} \int_{\mathbb{R}^n}\int_{\mathbb{R}^n}\frac{\Omega((y,z)/|(y,z)|)}{|(y,z)|^{2n}}f(x-y)g(x-z)\,dydz,$$ satisfies a sparse bound by $(p,p,p)$-averages, where $p$ is bigger than a certain number explicitly related to $r$ and $n$. As a consequence we deduce certain quantitative weighted estimates for bilinear homogeneous singular integrals associated with rough homogeneous kernels.


title: Unleashing In-network Computing on Scientific Workloads

abstract: Many recent efforts have shown that in-network computing can benefit various datacenter applications. In this paper, we explore a relatively less-explored domain which we argue can benefit from in-network computing: scientific workloads in high-performance computing. By analyzing canonical examples of HPC applications, we observe unique opportunities and challenges for exploiting in-network computing to accelerate scientific workloads. In particular, we find that the dynamic and demanding nature of scientific workloads is the major obstacle to the adoption of in-network approaches which are mostly open-loop and lack runtime feedback. In this paper, we present NSinC (Network-accelerated ScIeNtific Computing), an architecture for fully unleashing the potential benefits of in-network computing for scientific workloads by providing closed-loop runtime feedback to in-network acceleration services. We outline key challenges in realizing this vision and a preliminary design to enable acceleration for scientific applications.


title: A Visual Analytics Approach for Exploratory Causal Analysis: Exploration, Validation, and Applications

abstract: Using causal relations to guide decision making has become an essential analytical task across various domains, from marketing and medicine to education and social science. While powerful statistical models have been developed for inferring causal relations from data, domain practitioners still lack effective visual interface for interpreting the causal relations and applying them in their decision-making process. Through interview studies with domain experts, we characterize their current decision-making workflows, challenges, and needs. Through an iterative design process, we developed a visualization tool that allows analysts to explore, validate, and apply causal relations in real-world decision-making scenarios. The tool provides an uncertainty-aware causal graph visualization for presenting a large set of causal relations inferred from high-dimensional data. On top of the causal graph, it supports a set of intuitive user controls for performing what-if analyses and making action plans. We report on two case studies in marketing and student advising to demonstrate that users can effectively explore causal relations and design action plans for reaching their goals.


title: Bio-inspired Structure Identification in Language Embeddings

abstract: Word embeddings are a popular way to improve downstream performances in contemporary language modeling. However, the underlying geometric structure of the embedding space is not well understood. We present a series of explorations using bio-inspired methodology to traverse and visualize word embeddings, demonstrating evidence of discernible structure. Moreover, our model also produces word similarity rankings that are plausible yet very different from common similarity metrics, mainly cosine similarity and Euclidean distance. We show that our bio-inspired model can be used to investigate how different word embedding techniques result in different semantic outputs, which can emphasize or obscure particular interpretations in textual data.


title: Near-inertial wave critical layers over sloping bathymetry

abstract: This study describes a specific type of critical layer for near-inertial waves (NIWs) that forms when isopycnals run parallel to sloping bathymetry. Upon entering this slantwise critical layer, the group velocity of the waves decreases to zero and the NIWs become trapped and amplified, which can enhance mixing. A realistic simulation of anticyclonic eddies on the Texas-Louisiana shelf reveals that such critical layers can form where the eddies impinge onto the sloping bottom. Velocity shear bands in the simulation indicate that wind-forced NIWs are radiated downward from the surface in the eddies, bend upward near the bottom, and enter critical layers over the continental shelf, resulting in inertially-modulated enhanced mixing. Idealized simulations designed to capture this flow reproduce the wave propagation and enhanced mixing. The link between the enhanced mixing and wave trapping in the slantwise critical layer is made using ray-tracing and an analysis of the waves' energetics in the idealized simulations. An ensemble of simulations is performed spanning the relevant parameter space that demonstrates that the strength of the mixing is correlated with the degree to which NIWs are trapped in the critical layers. While the application here is for a shallow coastal setting, the mechanisms could be active in the open ocean as well where isopycnals align with bathymetry.


title: Profiling US Restaurants from Billions of Payment Card Transactions

abstract: A payment card (such as debit or credit) is one of the most convenient payment methods for purchasing goods and services. Hundreds of millions of card transactions take place across the globe every day, generating a massive volume of transaction data. The data render a holistic view of cardholder-merchant interactions, containing insights that can benefit various applications, such as payment fraud detection and merchant recommendation. However, utilizing these insights often requires additional information about merchants missing from the data owner's (i.e., payment company's) perspective. For example, payment companies do not know the exact type of product a merchant serves. Collecting merchant attributes from external sources for commercial purposes can be expensive. Motivated by this limitation, we aim to infer latent merchant attributes from transaction data. As proof of concept, we concentrate on restaurants and infer the cuisine types of restaurants from transactions. To this end, we present a framework for inferring the cuisine types of restaurants from transaction data. Our proposed framework consists of three steps. In the first step, we generate cuisine labels for a limited number of restaurants via weak supervision. In the second step, we extract a wide variety of statistical features and neural embeddings from the restaurant transactions. In the third step, we use deep neural networks (DNNs) to infer the remaining restaurants' cuisine types. The proposed framework achieved a 76.2% accuracy in classifying the US restaurants. To the best of our knowledge, this is the first framework to infer the cuisine types of restaurants by analyzing transaction data as the only source.


title: Complementary Dynamics

abstract: We present a novel approach to enrich arbitrary rig animations with elastodynamic secondary effects. Unlike previous methods which pit rig displacements and physical forces as adversaries against each other, we advocate that physics should complement artists intentions. We propose optimizing for elastodynamic displacements in the subspace orthogonal to displacements that can be created by the rig. This ensures that the additional dynamic motions do not undo the rig animation. The complementary space is high dimensional, algebraically constructed without manual oversight, and capable of rich high-frequency dynamics. Unlike prior tracking methods, we do not require extra painted weights, segmentation into fixed and free regions or tracking clusters. Our method is agnostic to the physical model and plugs into non-linear FEM simulations, geometric as-rigid-as-possible energies, or mass-spring models. Our method does not require a particular type of rig and adds secondary effects to skeletal animations, cage-based deformations, wire deformers, motion capture data, and rigid-body simulations.


title: Unifying Clustered and Non-stationary Bandits

abstract: Non-stationary bandits and online clustering of bandits lift the restrictive assumptions in contextual bandits and provide solutions to many important real-world scenarios. Though the essence in solving these two problems overlaps considerably, they have been studied independently. In this paper, we connect these two strands of bandit research under the notion of test of homogeneity, which seamlessly addresses change detection for non-stationary bandit and cluster identification for online clustering of bandit in a unified solution framework. Rigorous regret analysis and extensive empirical evaluations demonstrate the value of our proposed solution, especially its flexibility in handling various environment assumptions.