title: Infrared Variability due to Magnetic Pressure Driven Jets, Dust Ejection and Quasi-Puffed-Up Inner Rims

abstract: The interaction between a YSO stellar magnetic field and its protostellar disc can result in stellar accretional flows and outflows from the inner disc rim. Gas flows with a velocity component perpendicular to disc midplane subject particles to centrifugal acceleration away from the protostar, resulting in particles being catapulted across the face of the disc. The ejected material can produce a "dust fan", which may be dense enough to mimic the appearance of a "puffed-up" inner disc rim. We derive analytic equations for the time dependent disc toroidal field, the disc magnetic twist, the size of the stable toroidal disc region, the jet speed and the disc region of maximal jet flow speed. We show how the observed infrared variability of the pre-transition disc system LRLL~31 can be modelled by a dust ejecta fan from the inner-most regions of the disc whose height is partially dependent on the jet flow speed. The greater the jet flow speed, the higher is the potential dust fan scale height. An increase in mass accretion onto the star tends to increase the height and optical depth of the dust ejection fan, increasing the amount of 1--8~$\mu$m radiation. The subsequent shadow reduces the amount of light falling on the outer disc and decreases the 8-- 40~$\mu$m radiation. A decrease in the accretion rate reverses this scenario, thereby producing the observed "see-saw" infrared variability.


title: Modeling microbial cross-feeding at intermediate scale portrays community dynamics and species coexistence

abstract: Social interaction between microbes can be described at many levels of details, ranging from the biochemistry of cell-cell interactions to the ecological dynamics of populations. Choosing the best level to model microbial communities without losing generality remains a challenge. Here we propose to model cross-feeding interactions at an intermediate level between genome-scale metabolic models of individual species and consumer-resource models of ecosystems, which is suitable to empirical data. We applied our method to three published examples of multi-strain Escherichia coli communities with increasing complexity consisting of uni-, bi-, and multi-directional cross-feeding of either substitutable metabolic byproducts or essential nutrients. The intermediate-scale model accurately described empirical data and could quantify exchange rates elusive by other means, such as the byproduct secretions, even for a complex community of 14 amino acid auxotrophs. We used the three models to study each community's limits of robustness to perturbations such as variations in resource supply, antibiotic treatments and invasion by other "cheaters" species. Our analysis provides a foundation to quantify cross-feeding interactions from experimental data, and highlights the importance of metabolic exchanges in the dynamics and stability of microbial communities.


title: Interactive Natural Language-based Person Search

abstract: In this work, we consider the problem of searching people in an unconstrained environment, with natural language descriptions. Specifically, we study how to systematically design an algorithm to effectively acquire descriptions from humans. An algorithm is proposed by adapting models, used for visual and language understanding, to search a person of interest (POI) in a principled way, achieving promising results without the need to re-design another complicated model. We then investigate an iterative question-answering (QA) strategy that enable robots to request additional information about the POI's appearance from the user. To this end, we introduce a greedy algorithm to rank questions in terms of their significance, and equip the algorithm with the capability to dynamically adjust the length of human-robot interaction according to model's uncertainty. Our approach is validated not only on benchmark datasets but on a mobile robot, moving in a dynamic and crowded environment.


title: The G-stable rank for tensors

abstract: We introduce the $G$-stable rank of a higher order tensors over perfect fields. The $G$-stable rank is related to the Hilbert-Mumford criterion for stability in Geometric Invariant Theory. We will relate the $G$-stable rank to the tensor rank and slice rank. For numerical applications, we express the $G$-stable rank as a solution to an optimization problem. Over the field ${\mathbb F}_3$ we discuss an application to the Cap Set Problem.


title: Residual Bootstrap Exploration for Bandit Algorithms

abstract: In this paper, we propose a novel perturbation-based exploration method in bandit algorithms with bounded or unbounded rewards, called residual bootstrap exploration (\texttt{ReBoot}). The \texttt{ReBoot} enforces exploration by injecting data-driven randomness through a residual-based perturbation mechanism. This novel mechanism captures the underlying distributional properties of fitting errors, and more importantly boosts exploration to escape from suboptimal solutions (for small sample sizes) by inflating variance level in an \textit{unconventional} way. In theory, with appropriate variance inflation level, \texttt{ReBoot} provably secures instance-dependent logarithmic regret in Gaussian multi-armed bandits. We evaluate the \texttt{ReBoot} in different synthetic multi-armed bandits problems and observe that the \texttt{ReBoot} performs better for unbounded rewards and more robustly than \texttt{Giro} \cite{kveton2018garbage} and \texttt{PHE} \cite{kveton2019perturbed}, with comparable computational efficiency to the Thompson sampling method.


title: CopyCat: Controlled Instruction-Level Attacks on Enclaves

abstract: The adversarial model presented by trusted execution environments (TEEs) has prompted researchers to investigate unusual attack vectors. One particularly powerful class of controlled-channel attacks abuses page-table modifications to reliably track enclave memory accesses at a page-level granularity. In contrast to noisy microarchitectural timing leakage, this line of deterministic controlled-channel attacks abuses indispensable architectural interfaces and hence cannot be mitigated by tweaking microarchitectural resources. We propose an innovative controlled-channel attack, named CopyCat, that deterministically counts the number of instructions executed within a single enclave code page. We show that combining the instruction counts harvested by CopyCat with traditional, coarse-grained page-level leakage allows the accurate reconstruction of enclave control flow at a maximal instruction-level granularity. CopyCat can identify intra-page and intra-cache line branch decisions that ultimately may only differ in a single instruction, underscoring that even extremely subtle control flow deviations can be deterministically leaked from secure enclaves. We demonstrate the improved resolution and practicality of CopyCat on Intel SGX in an extensive study of single-trace and deterministic attacks against cryptographic implementations, and give novel algorithmic attacks to perform single-trace key extraction that exploit subtle vulnerabilities in the latest versions of widely-used cryptographic libraries. Our findings highlight the importance of stricter verification of cryptographic implementations, especially in the context of TEEs.


title: Fine tuning U-Net for ultrasound image segmentation: which layers?

abstract: Fine-tuning a network which has been trained on a large dataset is an alternative to full training in order to overcome the problem of scarce and expensive data in medical applications. While the shallow layers of the network are usually kept unchanged, deeper layers are modified according to the new dataset. This approach may not work for ultrasound images due to their drastically different appearance. In this study, we investigated the effect of fine-tuning different layers of a U-Net which was trained on segmentation of natural images in breast ultrasound image segmentation. Tuning the contracting part and fixing the expanding part resulted in substantially better results compared to fixing the contracting part and tuning the expanding part. Furthermore, we showed that starting to fine-tune the U-Net from the shallow layers and gradually including more layers will lead to a better performance compared to fine-tuning the network from the deep layers moving back to shallow layers. We did not observe the same results on segmentation of X-ray images, which have different salient features compared to ultrasound, it may therefore be more appropriate to fine-tune the shallow layers rather than deep layers. Shallow layers learn lower level features (including speckle pattern, and probably the noise and artifact properties) which are critical in automatic segmentation in this modality.


title: AdvMS: A Multi-source Multi-cost Defense Against Adversarial Attacks

abstract: Designing effective defense against adversarial attacks is a crucial topic as deep neural networks have been proliferated rapidly in many security-critical domains such as malware detection and self-driving cars. Conventional defense methods, although shown to be promising, are largely limited by their single-source single-cost nature: The robustness promotion tends to plateau when the defenses are made increasingly stronger while the cost tends to amplify. In this paper, we study principles of designing multi-source and multi-cost schemes where defense performance is boosted from multiple defending components. Based on this motivation, we propose a multi-source and multi-cost defense scheme, Adversarially Trained Model Switching (AdvMS), that inherits advantages from two leading schemes: adversarial training and random model switching. We show that the multi-source nature of AdvMS mitigates the performance plateauing issue and the multi-cost nature enables improving robustness at a flexible and adjustable combination of costs over different factors which can better suit specific restrictions and needs in practice.


title: Cooperative LIDAR Object Detection via Feature Sharing in Deep Networks

abstract: The recent advancements in communication and computational systems has led to significant improvement of situational awareness in connected and autonomous vehicles. Computationally efficient neural networks and high speed wireless vehicular networks have been some of the main contributors to this improvement. However, scalability and reliability issues caused by inherent limitations of sensory and communication systems are still challenging problems. In this paper, we aim to mitigate the effects of these limitations by introducing the concept of feature sharing for cooperative object detection (FS-COD). In our proposed approach, a better understanding of the environment is achieved by sharing partially processed data between cooperative vehicles while maintaining a balance between computation and communication load. This approach is different from current methods of map sharing, or sharing of raw data which are not scalable. The performance of the proposed approach is verified through experiments on Volony dataset. It is shown that the proposed approach has significant performance superiority over the conventional single-vehicle object detection approaches.


title: Diagnostics of plasma ionisation using torsional Alf\'en waves

abstract: Using the recently observed torsional Alfv\'en waves in solar prominences, we determine the ionisation state of the plasma by taking into account that Alfv\'en waves propagate in a partially ionised prominence plasma. We derive the evolutionary equation of waves and compare the analytical solutions to observations to determine the number density of neutrals. Using a single fluid plasma approximation, where the wave damping is provided by the Cowling resistivity, we study the temporal evolution of waves. By comparing the solution of equations with observational data (period, amplitude, propagation speed), we determined the value of the Cowling resistivity that led us to draw a conclusion on the amount of neutrals in the partially ionised plasma, a quantity that cannot be measured directly or indirectly. Our results show that damped torsional Alfv\'en waves are an ideal diagnostic tool for the ionisation state of the plasma. Using a simple model, we find that at the observational temperature of torsional Alfv\'en waves, the number of neutrals is of the order of $5\times 10^{10}$ cm$^{-3}$.


title: Detection of the Red Supergiant Wind from the Progenitor of Cassiopeia A

abstract: Cassiopeia A (Cas A) is one of the best studied young Galactic supernova remnants. While providing a rare opportunity to study in detail the remnant of a Type IIb supernova, questions remain regarding the nature of its progenitor, its mass-loss history, and its pre-SN evolution. Here we present an optical investigation of the circumstellar environment around Cas A and find clumpy and filamentary Halpha emission nebulosities concentrated 10-15 pc (10-15 arcminutes) to the north and east. First reported by Minkowski as a faint H II region, these nebulosities exhibit distinct morphological and spectroscopic properties relative to the surrounding diffuse emissions. Compared to neighboring H II regions, these nebulae show stronger [N II] 6548, 6583 A and [S II] 6716, 6731 A emissions relative to Halpha. We show that Cas A's highest-velocity ejecta knots are interacting with some of the closest projected emission nebulae, thus providing strong evidence that these nebulae lie at the same distance as the remnant. We interpret these surrounding nebulosities to be the remains of the progenitor's red supergiant wind which accumulated against the southern edge of a large extended H II region located north of Cas A. Our findings are consistent with the view that Cas A's progenitor underwent considerable mass-loss, first from a fast main-sequence wind, then from a slower, clumpy red supergiant wind, and finally from a brief high-velocity wind, like that from a yellow supergiant.


title: Simultaneous Inference for Massive Data: Distributed Bootstrap

abstract: In this paper, we propose a bootstrap method applied to massive data processed distributedly in a large number of machines. This new method is computationally efficient in that we bootstrap on the master machine without over-resampling, typically required by existing methods \cite{kleiner2014scalable,sengupta2016subsampled}, while provably achieving optimal statistical efficiency with minimal communication. Our method does not require repeatedly re-fitting the model but only applies multiplier bootstrap in the master machine on the gradients received from the worker machines. Simulations validate our theory.


title: On the universal unfolding of vector fields in one variable: A proof of Kostov's theorem

abstract: In this note we present variants of Kostov's theorem on a versal deformation of a parabolic point of a complex analytic $1$-dimensional vector field. First we provide a self-contained proof of Kostov's theorem, together with a proof that this versal deformation is indeed universal. We then generalize to the real analytic and formal cases, where we show universality, and to the $C^\infty$ case, where we show that only versality is possible.


title: Potential Coexistence of Exciton and Fermion Pair Condensations

abstract: Extensive theoretical and experimental investigation has been conducted on fermion pair condensation and exciton condensation as distinct classes of Bose-Einstein-like condensation. In this work, the existence of a fermion-exciton condensate---a single quantum state in which character of both fermion pair and exciton condensates coexist---is established computationally in the low-particle-number ($N$) limit and theoretically in the large-$N$ thermodynamic limit. The trade-off between the fermion pair and excitonic character of the fermion-exciton condensate is shown to be elliptic in nature. The possibility that the properties of fermion-exciton condensates could be a hybrid of the properties of fermion pair condensates and exciton condensates is discussed, and future experimental and computational exploration of this new class of condensate, which may potentially be realizable in a bilayer of superconductors, is anticipated.


title: Counterexamples to $L^p$ collapsing estimates

abstract: We show that certain $L^2$ space-time estimates for generalized density matrices which have been used by several authors in recent years to study equations of BBGKY or Hartree-Fock type, do not have non-trivial $L^pL^q$ generalizations.


title: Cataclysmic Variables in the First Year of the Zwicky Transient Facility

abstract: Using selection criteria based on amplitude, time and color, we have identified 329 objects as known or candidate cataclysmic variable (CVs) during the first year of testing and operation of the Zwicky Transient Facility (ZTF). Of these, 90 are previously confirmed CVs, 218 are strong candidates based on the shape and color of their light curves obtained during 3-562 days of observations, and the remaining 21 are possible CVs but with too few data points to be listed as good candidates. Almost half the strong candidates are within 10 deg of the galactic plane, in contrast to most other large surveys which have avoided crowded fields. The available Gaia parallaxes are consistent with sampling the low mass transfer CVs, as predicted by population models. Our followup spectra have confirmed Balmer/helium emission lines in 27 objects, with four showing high excitation HeII emission, including candidates for an AM CVn, a polar and an intermediate polar. Our results demonstrate that a complete survey of the galactic plane is needed to accomplish an accurate determination of the number of CVs existing in the Milky Way.


title: SD-GAN: Structural and Denoising GAN reveals facial parts under occlusion

abstract: Certain facial parts are salient (unique) in appearance, which substantially contribute to the holistic recognition of a subject. Occlusion of these salient parts deteriorates the performance of face recognition algorithms. In this paper, we propose a generative model to reconstruct the missing parts of the face which are under occlusion. The proposed generative model (SD-GAN) reconstructs a face preserving the illumination variation and identity of the face. A novel adversarial training algorithm has been designed for a bimodal mutually exclusive Generative Adversarial Network (GAN) model, for faster convergence. A novel adversarial "structural" loss function is also proposed, comprising of two components: a holistic and a local loss, characterized by SSIM and patch-wise MSE. Ablation studies on real and synthetically occluded face datasets reveal that our proposed technique outperforms the competing methods by a considerable margin, even for boosting the performance of Face Recognition.


title: Spin 2 Quasinormal Modes in Generalized Nariai Spacetimes

abstract: In this work we analytically obtain the quasinormal spectrum for the gravitational perturbation on a higher-dimensional generalization of the Nariai spacetime that is comprised of the direct product of the two-dimensional de Sitter space with several two-spheres. A key step in order to attain this result is to use a suitable basis for the angular functions depending on the rank of the tensorial degree of freedom that one needs to describe. Here we define such a basis, which is a generalization of the tensor spherical harmonics that is suited for spaces that are the product of several spaces of constant curvature.


title: Unitary Quantum Lattice Simulations for Maxwell Equations in Vacuum and in Dielectric Media

abstract: Utilizing the similarity between the spinor representation of the Dirac equation and the Maxwell equations that has been recognized since the early days of relativistic quantum mechanics, a quantum lattice (QLA) representation of unitary collision stream operators of Maxwell equations is derived for both homogeneous and inhomogeneous media. A second order accurate 4 spinor scheme is developed and tested successfully for two dimensional (2D) propagation of a Gaussian pulse in a uniform medium while for normal (1D) incidence of an electromagnetic Gaussian pulse onto a dielectric interface requires 8 component spinors. In particular, the well-known phase change, field amplitudes and profile widths are recovered by the QLA asymptotic profiles. The QLA simulations yield the time dependent electromagnetic fields as the pulse enters and straddles the dielectric boundary. QLA involves unitary interleaved noncommuting collision and streaming operators that can be coded onto a quantum computer. The noncommutation being the only reason why one perturbatively recovers the Maxwell equations.


title: Sr$_2$RuO$_4$, like doped cuprates and barium bismuthate, is a negative charge-transfer gap even parity superconductor with $\frac{3}{4}$-filled oxygen band

abstract: A comprehensive theory of superconductivity in Sr$_2$RuO$_4$ must explain experiments that suggest even parity superconducting order and others that suggest broken time reversal symmetry. Completeness further requires that the theory applies to Ca$_2$RuO$_4$, a Mott-Hubbard semiconductor that exhibits an unprecedented insulator-to-metal transition driven by very small electric field, and also by doping with very small concentration of electrons, leading to a metallic state proximate to ferromagnetism. A valence transition model, previously proposed for superconducting cuprates [Phys. Rev. B {\bf 98}, 205153] is extended to Sr$_2$RuO$_4$ and Ca$_2$RuO$_4$. The insulator to metal transition is distinct from that expected from the simple melting of the Mott-Hubbard semiconductor. Rather, the Ru ions occur as low spin Ru$^{4+}$ in the semiconductor, and as high spin Ru$^{3+}$ in the metal, the driving force behind the valence transition being the strong spin-charge coupling and consequent large ionizaton energy in the low charge state. Metallic and superconducting ruthenates are two-component systems in which the half-filled high spin Ru$^{3+}$ ions determine the magnetic behavior but not transport, while the charge carriers are entirely on on the layer oxygen ions, which have an average charge -1.5. Spin singlet superconductivity evolves from the correlated lattice frustrated 3/4 filled band of layer oxygen ions alone, in agreement with quantum many body calculations that have demonstrated enhancement by electron-electron interactions of superconducting pair-pair correlations tions uniquely at or very close to this filling [Phys. Rev. B {\bf 93}, 165110 and {\bf 93}, 205111]. Several model specific experimental predictions are made, including that spin susceptibility due to Ru ions will remain unchanged as Sr$_2$RuO$_4$ is taken through superconducting Tc.


title: Finiteness results concerning algebraic power series

abstract: We construct an explicit filtration of the ring of algebraic power series by finite dimensional constructible sets, measuring the complexity of these series. As an application, we give a bound on the dimension of the set of algebraic power series of bounded complexity lying on an algebraic variety defined over the field of power series.


title: q-Measurement Scheme set in Quantum Photonic Frameworks: Linking Abstract Space and Lab Spaces

abstract: Beginning in abstract space and dislodging the representational form paves a way to formulate a version of a quantum physical measurement scheme. With materiality playing sustainment roles with respect to q-states, these latter control system response towards external probes. Sustainment concerns all possible q-states. The particle character usually assigned to eigenstates and eigenvalues dies out, while energy gaps relate to possible energy quanta exchanges; the occupancy concept no longer required. Abstract and laboratory spaces are incommensurate, linking them leads to a possible physical q-measurement theory. Connectors (links) are referred to as gangplanks and/or gangways, they perform in frontiers arbitrating quantum to classical regions. Inertial frames facilitate linking q-state domains via projection to wave function space. A physical quantum measurement theory obtains grounded on these elements including quantum information concepts.


title: A Recurrent Neural Network Based Patch Recommender for Linux Kernel Bugs

abstract: Software bugs in a production environment have an undesirable impact on quality of service, unplanned system downtime, and disruption in good customer experience, resulting in loss of revenue and reputation. Existing approaches to automated software bug repair focuses on known bug templates detected using static code analysis tools and test suites, and in automatic generation of patch code for these bugs. We describe the typical bug fixing process employed in the Linux kernel, and motivate the need for a new automated tool flow to fix bugs. We present an initial design of such an automated tool that uses Recurrent Neural Network (RNN) based Natural Language Processing to generate patch recommendations from user generated bug reports. At the 50th percentile of the test bugs, the correct patch occurs within the top 11.5 patch recommendations output by the model. Further, we present a Linux kernel developer's assessment of the quality of patches recommended for new unresolved kernel bugs.


title: EyeTAP: A Novel Technique using Voice Inputs to Address the Midas Touch Problem for Gaze-based Interactions

abstract: One of the main challenges of gaze-based interactions is the ability to distinguish normal eye function from a deliberate interaction with the computer system, commonly referred to as 'Midas touch'. In this paper we propose, EyeTAP (Eye tracking point-and-select by Targeted Acoustic Pulse) a hands-free interaction method for point-and-select tasks. We evaluated the prototype in two separate user studies, each containing two experiments with 33 participants and found that EyeTAP is robust even in presence of ambient noise in the audio input signal with tolerance of up to 70 dB, results in a faster movement time, and faster task completion time, and has a lower cognitive workload than voice recognition. In addition, EyeTAP has a lower error rate than the dwell-time method in a ribbon-shaped experiment. These characteristics make it applicable for users for whom physical movements are restricted or not possible due to a disability. Furthermore, EyeTAP has no specific requirements in terms of user interface design and therefore it can be easily integrated into existing systems with minimal modifications. EyeTAP can be regarded as an acceptable alternative to address the Midas touch.


title: From Poincar\'e Recurrence to Convergence in Imperfect Information Games: Finding Equilibrium via Regularization

abstract: In this paper we investigate the Follow the Regularized Leader dynamics in sequential imperfect information games (IIG). We generalize existing results of Poincar\'e recurrence from normal-form games to zero-sum two-player imperfect information games and other sequential game settings. We then investigate how adapting the reward (by adding a regularization term) of the game can give strong convergence guarantees in monotone games. We continue by showing how this reward adaptation technique can be leveraged to build algorithms that converge exactly to the Nash equilibrium. Finally, we show how these insights can be directly used to build state-of-the-art model-free algorithms for zero-sum two-player Imperfect Information Games (IIG).


title: ivmodel: An R Package for Inference and Sensitivity Analysis of Instrumental Variables Models with One Endogenous Variable

abstract: We present a comprehensive R software ivmodel for analyzing instrumental variables with one endogenous variable. The package implements a general class of estimators called k- class estimators and two confidence intervals that are fully robust to weak instruments. The package also provides power formulas for various test statistics in instrumental variables. Finally, the package contains methods for sensitivity analysis to examine the sensitivity of the inference to instrumental variables assumptions. We demonstrate the software on the data set from Card (1995), looking at the causal effect of levels of education on log earnings where the instrument is proximity to a four-year college.


title: Caveats in Eliciting Mobile App Requirements

abstract: Factors such as app stores or platform choices heavily affect functional and non-functional mobile app requirements. We surveyed 45 companies and interviewed ten experts to explore how factors that impact mobile app requirements are understood by requirements engineers in the mobile app industry. We observed a lack of knowledge in several areas. For instance, we observed that all practitioners were aware of data privacy concerns, however, they did not know that certain third-party libraries, usage aggregators, or advertising libraries also occasionally leak sensitive user data. Similarly, certain functional requirements may not be implementable in the absence of a third-party library that is either banned from an app store for policy violations or lacks features, for instance, missing desired features in ARKit library for iOS made practitioners turn to Android. We conclude that requirements engineers should have adequate technical experience with mobile app development as well as sufficient knowledge in areas such as privacy, security and law, in order to make informed decisions during requirements elicitation.


title: Regularity of Lyapunov Exponents for Diffeomorphisms with Dominated Splitting

abstract: We consider families of diffeomorphisms with dominated splittings and preserving a Borel probability measure, and we study the regularity of the Lyapunov exponents associated to the invariant bundles with respect to the parameter. We obtain that the regularity is at least the sum of the regularities of the two invariant bundles (for regularities in $[0,1]$), and under suitable conditions we obtain formulas for the derivatives. Similar results are obtained for families of flows, and for the case when the invariant measure depends on the map. We also obtain several applications. Near the time one map of a geodesic flow of a surface of negative curvature the metric entropy of the volume is Lipschitz with respect to the parameter. At the time one map of a geodesic flow on a manifold of constant negative curvature the topological entropy is differentiable with respect to the parameter, and we give a formula for the derivative. Under some regularity conditions, the critical points of the Lyapunov exponent function are non-flat (the second derivative is nonzero for some families). Also, again under some regularity conditions, the criticality of the Lyapunov exponent function implies some rigidity of the map, in the sense that the volume decomposes as a product along the two complimentary foliations. In particular for area preserving Anosov diffeomorphisms, the only critical points are the maps smoothly conjugated to the linear map, corresponding to the global extrema.


title: Absence of $\mu$SR evidence for magnetic order in the pseudogap phase of Bi$_{2+x}$Sr$_{2-x}$CaCu$_2$O$_{8+\delta}$

abstract: We present an extended zero-field muon spin relaxation (ZF-$\mu$SR) study of overdoped Bi$_{2+x}$Sr$_{2-x}$CaCu$_2$O$_{8+\delta}$ (Bi2212) single crystals, intended to elucidate the origin of weak quasistatic magnetism previously detected by $\mu$SR in the superconducting and normal states of optimally-doped and overdoped samples. New results on heavily-overdoped single crystals show a similar monotonically decreasing ZF-$\mu$SR relaxation rate with increasing temperature that persists above the pseudogap (PG) temperature $T^*$ and does not evolve with hole doping ($p$). Additional measurements using an ultra-low background apparatus confirm that this behavior is an intrinsic property of Bi2212, which cannot be due to magnetic order associated with the PG phase. Instead we show that the temperature-dependent relaxation rate is most likely caused by structural changes that modify the contribution of the nuclear dipole fields to the ZF-$\mu$SR signal. Our results for Bi2212 emphasize the importance of not assuming the nuclear-dipole field contribution is independent of temperature in ZF-$\mu$SR studies of high-temperature (high-$T_c$) cuprate superconductors, and do not support a recent $\mu$SR study of YBa$_2$Cu$_3$O$_{6+x}$ that claims to detect magnetic order in the PG phase.


title: An Alternating Direction Explicit Method for Time Evolution Equations with Applications to Fractional Differential Equations

abstract: We derive and analyze the alternating direction explicit (ADE) method for time evolution equations with the time-dependent Dirichlet boundary condition and with the zero Neumann boundary condition. The original ADE method is an additive operator splitting (AOS) method, which has been developed for treating a wide range of linear and nonlinear time evolution equations with the zero Dirichlet boundary condition. For linear equations, it has been shown to achieve the second order accuracy in time yet is unconditionally stable for an arbitrary time step size. For the boundary conditions considered in this work, we carefully construct the updating formula at grid points near the boundary of the computational domain and show that these formulas maintain the desired accuracy and the property of unconditional stability. We also construct numerical methods based on the ADE scheme for two classes of fractional differential equations. We will give numerical examples to demonstrate the simplicity and the computational efficiency of the method.


title: Inter-numerology Interference Management with Adaptive Guards: A Cross-layer Approach

abstract: The next-generation communication technologies are evolving towards increased flexibility in various aspects. Although orthogonal frequency division multiplexing (OFDM) remains as the waveform of the upcoming fifth-generation (5G) standard, the new radio provides flexibility in waveform parametrization (a.k.a. numerology) to address diverse requirements. However, managing the peaceful coexistence of mixed numerologies is challenging due to inter-numerology interference (INI). This paper proposes the utilization of adaptive guards in both time and frequency domains as a solution along with a multi-window operation in the physical (PHY) layer. The adaptive windowing operation needs a guard duration to reduce the unwanted emissions, and a guard band is required to handle the INI level on the adjacent band. The guards in both domains are jointly optimized with respect to the subcarrier spacing, use case (i.e., service requirement), and power offset between the numerologies. Also, the multi-window approach provides managing each side of the spectrum independently in case of an asymmetric interference scenario. Since the allowed interference level depends on the numerologies operating in the adjacent bands, the potential of adaptive guards is further increased and exploited with a medium access control (MAC) layer scheduling technique. The proposed INI-based scheduling algorithm decreases the need for guards by allocating the numerologies to the available bands, considering their subcarrier spacing, power level, and SIR requirements. Therefore, INI management is performed with a cross-layer (PHY and MAC) approach in this study. The results show that the precise design that accommodates such flexibility reduces the guards significantly and improves the spectral efficiency of mixed numerology systems.


title: Tricking Johnny into Granting Web Permissions

abstract: We studied the web permission API dialog box in popular mobile and desktop browsers, and found that it typically lacks measures to protect users from unwittingly granting web permission when clicking too fast. We developed a game that exploits this issue, and tricks users into granting webcam permission. We conducted three experiments, each with 40 different participants, on both desktop and mobile browsers. The results indicate that in the absence of a prevention mechanism, we achieve a considerably high success rate in tricking 95% and 72% of participants on mobile and desktop browsers, respectively. Interestingly, we also tricked 47% of participants on a desktop browser where a prevention mechanism exists.


title: A fitted finite volume method for stochastic optimal control Problems

abstract: In this article, we provide a numerical method based on fitted finite volume method to approximate the Hamilton-Jacobi-Bellman (HJB) equation coming from stochastic optimal control problems. The computational challenge is due to the nature of the HJB equation, which may be a second-order degenerated partial differential equation coupled with optimization. In the work, we discretize the HJB equation using the fitted finite volume method and show that matrix resulting from spatial discretization is an M-matrix. The optimization problem is solved at every time step using iterative method. Numerical results are presented to show the robustness of the fitted finite volume numerical method comparing to the standard finite difference method.


title: Descriptive and Predictive Analysis of Euroleague Basketball Games and the Wisdom of Basketball Crowds

abstract: In this study we focus on the prediction of basketball games in the Euroleague competition using machine learning modelling. The prediction is a binary classification problem, predicting whether a match finishes 1 (home win) or 2 (away win). Data is collected from the Euroleague's official website for the seasons 2016-2017, 2017-2018 and 2018-2019, i.e. in the new format era. Features are extracted from matches' data and off-the-shelf supervised machine learning techniques are applied. We calibrate and validate our models. We find that simple machine learning models give accuracy not greater than 67% on the test set, worse than some sophisticated benchmark models. Additionally, the importance of this study lies in the "wisdom of the basketball crowd" and we demonstrate how the predicting power of a collective group of basketball enthusiasts can outperform machine learning models discussed in this study. We argue why the accuracy level of this group of "experts" should be set as the benchmark for future studies in the prediction of (European) basketball games using machine learning.


title: Criptocurrencies, Fiat Money, Blockchains and Databases

abstract: Two taxonomies of money that include cryptocurrencies are analyzed. A definition of the term cryptocurrency is given and a taxonomy of them is presented, based on how its price is fixed. The characteristics of the use of current fiat money and the operation of two-level banking systems are discussed. Cryptocurrencies are compared with fiat money and the aspects in which the latter cannot be overcome are indicated. The characteristics of blockchains and databases are described. The possible cases of use of both technologies are compared, and it is noted that blockchains, in addition to cryptocurrencies and certain records, have not yet shown their usefulness, while databases constitute the foundation of most of the automated systems in operation.


title: Analysis on Complete Set of Fock States with Explicit Wavefunctions for the Covariant Harmonic Oscillator Problem

abstract: The earlier treatments of Lorentz covariant harmonic oscillator have brought to light various difficulties, such as reconciling Lorentz symmetry with the full Fock space, and divergence issues with their functional representations. We present here a full solution avoiding those problems. The complete set of Fock states is obtained, together with the corresponding explicit wavefunction and their inner product integrals free from any divergence problem and the Lorentz symmetry fully maintained without additional constraints imposed. By a simple choice of the pseudo-unitary representation of the underlying symmetry group, motivated from the perspective of the Minkowski spacetime as a representation for the Lorentz group, we obtain the natural non-unitary Fock space picture commonly considered though not formulated and presented in the careful details given here. From a direct derivation of the appropriate basis state wavefunctions of the finite-dimensional irreducible representations of the Lorentz symmetry, the relation between the latter and the Fock state wavefunctions is also explicitly shown. Moreover, the full picture including the states with non-positive norm may give consistent physics picture as a version of Lorentz covariant quantum mechanics. Probability interpretation for the usual von Neumann measurements is not a problem as all wavefunctions restricted to a definite value for the `time' variable are just like those of the usual time independent quantum mechanics. A further understanding from a perspective of the dynamics from the symplectic geometry of the phase space is shortly discussed.


title: Effective Correlates of Motor Imagery Performance based on Default Mode Network in Resting-State

abstract: Motor imagery based brain-computer interfaces (MI-BCIs) allow the control of devices and communication by imagining different muscle movements. However, most studies have reported a problem of "BCI-illiteracy" that does not have enough performance to use MI-BCI. Therefore, understanding subjects with poor performance and finding the cause of performance variation is still an important challenge. In this study, we proposed predictors of MI performance using effective connectivity in resting-state EEG. As a result, the high and low MI performance groups had a significant difference as 23% MI performance difference. We also found that connection from right lateral parietal to left lateral parietal in resting-state EEG was correlated significantly with MI performance (r = -0.37). These findings could help to understand BCI-illiteracy and to consider alternatives that are appropriate for the subject.


title: Prediction of Individual Propofol Requirements based on Preoperative EEG Signals

abstract: The patient must be given an adequate amount of propofol for safe surgery since overcapacity and low capacity cause accidents. However, the sensitivity of propofol varies from patient to patient, making it very difficult to determine the propofol requirements for anesthesia. This paper aims to propose a neurophysiological predictor of propofol requirements based on the preoperative electroencephalogram (EEG). We exploited the canonical correlation analysis that infers the amount of information on the propofol requirements. The results showed that the preoperative EEG included the factor that could explain the propofol requirements. Specifically, the frontal and posterior regions had crucial information on the propofol requirements. Moreover, there was a significantly different power in the frontal and posterior regions between baseline and unconsciousness periods, unlike the alpha power in the central region. These findings showed the potential that preoperative EEG could predict the propofol requirements.


title: The Information in Emotion Communication

abstract: How much information is transmitted when animals use emotions to communicate? It is clear that emotions are used as communication systems in humans and other species. The quantitative theory of emotion information presented here is based on Shannon's mathematical theory of information in communication systems. The theory explains myriad aspects of emotion communication and offers dozens of new directions for research. It is superior to the "contagion" theory of emotion spreading, which is currently dominant. One important application of the information theory of emotion communication is that it permits the development of emotion security systems for social networks to guard against the widespread emotion manipulation we see online today.


title: Scaled Fixed Point Algorithm for Computing the Matrix Square Root

abstract: This paper addresses the numerical solution of the matrix square root problem. Two fixed point iterations are proposed by rearranging the nonlinear matrix equation $A - X^2 = 0$ and incorporating a positive scaling parameter. The proposals only need to compute one matrix inverse and at most two matrix multiplications per iteration. A global convergence result is established. The numerical comparisons versus some existing methods from the literature, on several test problems, demonstrate the efficiency and effectiveness of our proposals.


title: Prescriptions for Correcting Ultraviolet-Based Redshifts for Luminous Quasars at High Redshift

abstract: High-redshift quasars typically have their redshift determined from rest-frame ultraviolet (UV) emission lines. However, these lines, and more specifically the prominent C IV $\lambda 1549$ emission line, are typically blueshifted yielding highly uncertain redshift estimates compared to redshifts determined from rest-frame optical emission lines. We present near-infrared spectroscopy of 18 luminous quasars at $2.15 < z < 3.70$ that allows us to obtain reliable systemic redshifts for these sources. Together with near-infrared spectroscopy of an archival sample of 44 quasars with comparable luminosities and redshifts, we provide prescriptions for correcting UV-based redshifts. Our prescriptions reduce velocity offsets with respect to the systemic redshifts by $\sim140$ km s$^{-1}$ and reduce the uncertainty on the UV-based redshift by $\sim25\%$ with respect to the best method currently used for determining such values. We also find that the redshifts determined from the Sloan Digital Sky Survey Pipeline for our sources suffer from significant uncertainties, which cannot be easily mitigated. We discuss the potential of our prescriptions to improve UV-based redshift corrections given a much larger sample of high redshift quasars with near-infrared spectra.


title: Revisiting Training Strategies and Generalization Performance in Deep Metric Learning

abstract: Deep Metric Learning (DML) is arguably one of the most influential lines of research for learning visual similarities with many proposed approaches every year. Although the field benefits from the rapid progress, the divergence in training protocols, architectures, and parameter choices make an unbiased comparison difficult. To provide a consistent reference point, we revisit the most widely used DML objective functions and conduct a study of the crucial parameter choices as well as the commonly neglected mini-batch sampling process. Under consistent comparison, DML objectives show much higher saturation than indicated by literature. Further based on our analysis, we uncover a correlation between the embedding space density and compression to the generalization performance of DML models. Exploiting these insights, we propose a simple, yet effective, training regularization to reliably boost the performance of ranking-based DML models on various standard benchmark datasets. Code and a publicly accessible WandB-repo are available at https://github.com/Confusezius/Revisiting_Deep_Metric_Learning_PyTorch.


title: Online Policies for Efficient Volunteer Crowdsourcing

abstract: Nonprofit crowdsourcing platforms such as food recovery organizations rely on volunteers to perform time-sensitive tasks. To encourage volunteers to complete a task, platforms use nudging mechanisms to notify a subset of volunteers with the hope that at least one of them responds positively. However, since excessive notifications may reduce volunteer engagement, the platform faces a trade-off between notifying more volunteers for the current task and saving them for future ones. Motivated by these applications, we introduce the online volunteer notification problem, a generalization of online stochastic bipartite matching where tasks arrive following a known time-varying distribution over task types. Upon arrival of a task, the platform notifies a subset of volunteers with the objective of minimizing the number of missed tasks. To capture each volunteer's adverse reaction to excessive notifications, we assume that a notification triggers a random period of inactivity, during which she will ignore all notifications. However, if a volunteer is active and notified, she will perform the task with a given pair-specific match probability that captures her preference for the task. We develop an online randomized policy that achieves a constant-factor guarantee close to the upper bound we establish for the performance of any online policy. Our policy as well as hardness results are parameterized by the minimum discrete hazard rate of the inter-activity time distribution. The design of our policy relies on sparsifying an ex-ante feasible solution by solving a sequence of dynamic programs. Further, in collaboration with Food Rescue U.S., a volunteer-based food recovery platform, we demonstrate the effectiveness of our policy by testing them on the platform's data from various locations across the U.S.


title: Fast Multi-Subset Transform and Weighted Sums Over Acyclic Digraphs

abstract: The zeta and Moebius transforms over the subset lattice of $n$ elements and the so-called subset convolution are examples of unary and binary operations on set functions. While their direct computation requires $O(3^n)$ arithmetic operations, less naive algorithms only use $2^n \mathrm{poly}(n)$ operations, nearly linear in the input size. Here, we investigate a related $n$-ary operation that takes $n$ set functions as input and maps them to a new set function. This operation, we call multi-subset transform, is the core ingredient in the known inclusion--exclusion recurrence for weighted sums over acyclic digraphs, which extends Robinson's recurrence for the number of labelled acyclic digraphs. Prior to this work the best known complexity bound was the direct $O(3^n)$. By reducing the task to multiple instances of rectangular matrix multiplication, we improve the complexity to $O(2.985^n)$.


title: A non-inferiority test for R-squared with random regressors

abstract: Determining the lack of association between an outcome variable and a number of different explanatory variables is frequently necessary in order to disregard a proposed model. This paper proposes a non-inferiority test for the coefficient of determination (or squared multiple correlation coefficient), R-squared, in a linear regression analysis with random predictors. The test is derived from inverting a one-sided confidence interval based on a scaled central F distribution.


title: Optimally Guarding Perimeters and Regions with Mobile Range Sensors

abstract: We investigate the problem of using mobile robots equipped with 2D range sensors to optimally guard perimeters or regions, i.e., 1D or 2D sets. Given such a set of arbitrary shape to be guarded, and $k$ mobile sensors where the $i$-th sensor can guard a circular region with a variable radius $r_i$, we seek the optimal strategy to deploy the $k$ sensors to fully cover the set such that $\max r_i$ is minimized. On the side of computational complexity, we show that computing a $1.152$-optimal solution for guarding a perimeter or a region is NP-hard, i.e., the problem is hard to approximate. The hardness result on perimeter guarding holds when each sensor may guard at most two disjoint perimeter segments. On the side of computational methods, for the guarding perimeters, we develop a fully polynomial time approximation scheme (FPTAS) for the special setting where each sensor may only guard a single continuous perimeter segment, suggesting that the aforementioned hard-to-approximate result on the two-disjoint-segment sensing model is tight. For the general problem, we first describe a polynomial-time (2+$\epsilon)$-approximation algorithm as an upper bound, applicable to both perimeter guarding and region guarding. This is followed by a high-performance integer linear programming (ILP) based method that computes near-optimal solutions. Thorough computational benchmarks as well as evaluation on potential application scenarios demonstrate the effectiveness of these algorithmic solutions.


title: Non-Hermitian quantum impurity systems in and out of equilibrium: noninteracting case

abstract: We provide systematic analysis on a non-Hermitian PT -symmetric quantum impurity system both in and out of equilibrium, based on exact computations. In order to understand the interplay between non-Hermiticity and Kondo physics, we focus on a prototypical noninteracting impurity system, the resonant level model, with complex coupling constants. Explicitly constructing biorthogonal basis, we study its thermodynamic properties as well as the Loschmidt echo starting from the initially disconnected two free fermion chains. Remarkably, we observe the universal crossover physics in the Loschmidt echo, both in the PT broken and unbroken regimes. We also find that the ground state quantities we compute in the PT broken regime can be obtained by analytic continuation. It turns out that Kondo screening ceases to exist in the PT broken regime, which was also previously predicted in the non-hermitian Kondo model. All the analytical results are corroborated against biorthogonal free fermion numerics.


title: Periodic orbits of Linear and invariant flows on connected Lie groups

abstract: Our main is to study periodic orbits of linear and invariant flows on a real, connected Lie group. Since each linear flow $\varphi_t$ has a derivation associated $\mathcal{D}$, we show that the existence of periodic orbits of $\varphi_t$ is based on the eigenvalues of the derivation $\mathcal{D}$. From this, we study periodic orbits of a linear flow on noncompact, semisimple Lie groups, and we work with periodic orbits of a linear flow on connected, simply connected, solvable Lie groups of dimension 2 or 3.


title: On contact loci of hyperplane arrangements

abstract: We give an explicit expression for the contact loci of hyperplane arrangements and show that their cohomology rings are combinatorial invariants. We also give an expression for the restricted contact loci in terms of Milnor fibers of associated hyperplane arrangements. We prove the degeneracy of a spectral sequence related to the restricted contact loci of a hyperplane arrangement and which conjecturally computes algebraically the Floer cohomology of iterates of the Milnor monodromy. We give formulas for the Betti numbers of contact loci and restricted contact loci in generic cases.


title: A novel approach to building micromegas detectors enabled by precision manufacturing

abstract: Micromegas detectors are a relatively modern design concept for micropattern gas detectors, designed to handle high particle flux while providing high gain, high spatial resolution, and fast response times for a variety of radiation detection applications. Due to the advancement of industry, building a micromegas detector without advanced in-house capabilities is now possible. In this work, we present an innovative method to build micromegas detectors utilizing precision manufacturers to fabricate the core components. Two detectors were built using the newly described method and are experimentally validated.


title: Assessment of the beta-delayed proton decay rate of $^{11}$Be

abstract: The $^{11}$Be neutron halo nucleus appears to decay into $^{10}$Be with a rate that exceeds expectations. Neutron disappearance into dark matter, beta decay of a halo neutron, or beta-delayed proton decay have been offered as explanations. In this work we study the latter option; we carry out shell model calculations and sequential decay analysis examining the beta-delayed proton decay going through a resonance in $^{11}$B. The narrow energy window, lack of states with sufficient spectroscopic strength, overwhelming alpha decay branch, all make reconciling the observed rate with beta-delayed proton decay difficult.


title: Strength from Weakness: Fast Learning Using Weak Supervision

abstract: We study generalization properties of weakly supervised learning. That is, learning where only a few "strong" labels (the actual target of our prediction) are present but many more "weak" labels are available. In particular, we show that having access to weak labels can significantly accelerate the learning rate for the strong task to the fast rate of $\mathcal{O}(\nicefrac1n)$, where $n$ denotes the number of strongly labeled data points. This acceleration can happen even if by itself the strongly labeled data admits only the slower $\mathcal{O}(\nicefrac{1}{\sqrt{n}})$ rate. The actual acceleration depends continuously on the number of weak labels available, and on the relation between the two tasks. Our theoretical results are reflected empirically across a range of tasks and illustrate how weak labels speed up learning on the strong task.


title: Estimating Training Data Influence by Tracing Gradient Descent

abstract: We introduce a method called TracIn that computes the influence of a training example on a prediction made by the model. The idea is to trace how the loss on the test point changes during the training process whenever the training example of interest was utilized. We provide a scalable implementation of TracIn via: (a) a first-order gradient approximation to the exact computation, (b) saved checkpoints of standard training procedures, and (c) cherry-picking layers of a deep neural network. In contrast with previously proposed methods, TracIn is simple to implement; all it needs is the ability to work with gradients, checkpoints, and loss functions. The method is general. It applies to any machine learning model trained using stochastic gradient descent or a variant of it, agnostic of architecture, domain and task. We expect the method to be widely useful within processes that study and improve training data.


title: The Minimal Genus of Homology Classes in a Finite 2-Complex

abstract: We study surface representatives of homology classes of finite complexes which minimize certain complexity measures, including its genus and Euler characteristic. Our main result is that up to surgery at nullhomotopic curves minimizers are homotopic to cellwise coverings to the 2-skeleton. From this we conclude that the minimizing problem is in general algorithmically undecidable, but can be solved for 2-dimensional CAT(-1)-complexes.


title: The accuracy of post-processing nucleosynthesis

abstract: The computational requirements posed by multi-dimensional simulations of type Ia supernovae make it difficult to incorporate complex nuclear networks to follow the release of nuclear energy along with the propagation of the flame. Instead, these codes usually model the flame and use simplified nuclear kinetics, with the goal of determining a sufficiently accurate rate of nuclear energy generation and, afterwards, post-processing the thermodynamic trajectories with a large nuclear network to obtain more reliable nuclear yields. In this work, I study the performance of simplified nuclear networks with respect to reproduction of the nuclear yields obtained with a one-dimensional supernova code equipped with a large nuclear network. I start by defining a strategy to follow the properties of matter in nuclear statistical equilibrium (NSE). I propose to use published tables of NSE properties, together with a careful interpolation routine. Short networks (iso7 and 13{\alpha}) are able to give an accurate yield of 56Ni, after post-processing, but can fail by order of magnitude in predicting the ejected mass of even mildly abundant species (> 0.001 solar masses). A network of 21 species reproduces the nucleosynthesis of the Chandrasekhar and sub-Chandrasekhar explosions studied here with average errors better than 20% for the whole set of stable elements and isotopes followed in the models.


title: More on Quantum Chiral Higher Spin Gravity

abstract: Chiral Higher Spin Gravity is unique in being the smallest higher spin extension of gravity and in having a simple local action both in flat and (anti)-de Sitter spaces. It must be a closed subsector of any other higher spin theory in four dimensions, which makes it an important building block and benchmark. Using the flat space version for simplicity, we perform a thorough study of quantum corrections in Chiral Theory, which strengthen our earlier results arXiv:1805.00048. Even though the interactions are naively non-renormalizable, we show that there are no UV-divergences in two-, three- and four-point amplitudes at one loop thanks to the higher spin symmetry. We also give arguments that the AdS Chiral Theory should exhibit similar properties. It is shown that Chiral Theory admits Yang-Mills gaugings with $U(N)$, $SO(N)$ and $USp(N)$ groups, which is reminiscent of the Chan-Paton symmetry in string theory.


title: Massive MIMO with Multi-Antenna Users under Jointly Correlated Ricean Fading

abstract: We study the uplink performance of massive multiple-input multiple-output (MIMO) when users are equipped with multiple antennas. To this end, we consider a generalized channel model that accounts for line-of-sight propagation and spatially correlated multipath fading. Most importantly, we employ the Weichselberger correlation model, which has been shown to alleviate the deficiencies of the popular Kronecker model. The main contribution of this paper is a rigorous closed-form expression for the uplink spectral efficiency using maximum-ratio combining and minimum mean square error channel estimation. Our result is a non-trivial generalization of previous results on massive MIMO with spatially correlated channels, thereby enabling us to have suitable designs for future massive MIMO systems. Numerical simulations corroborate our analysis and provide useful insights on how different propagation conditions affect system performance.


title: On the Versatility of Open Logical Relations: Continuity, Automatic Differentiation, and a Containment Theorem

abstract: Logical relations are one of the most powerful techniques in the theory of programming languages, and have been used extensively for proving properties of a variety of higher-order calculi. However, there are properties that cannot be immediately proved by means of logical relations, for instance program continuity and differentiability in higher-order languages extended with real-valued functions. Informally, the problem stems from the fact that these properties are naturally expressed on terms of non-ground type (or, equivalently, on open terms of base type), and there is no apparent good definition for a base case (i.e. for closed terms of ground types). To overcome this issue, we study a generalization of the concept of a logical relation, called \emph{open logical relation}, and prove that it can be fruitfully applied in several contexts in which the property of interest is about expressions of first-order type. Our setting is a simply-typed $\lambda$-calculus enriched with real numbers and real-valued first-order functions from a given set, such as the one of continuous or differentiable functions. We first prove a containment theorem stating that for any such a collection of functions including projection functions and closed under function composition, any well-typed term of first-order type denotes a function belonging to that collection. Then, we show by way of open logical relations the correctness of the core of a recently published algorithm for forward automatic differentiation. Finally, we define a refinement-based type system for local continuity in an extension of our calculus with conditionals, and prove the soundness of the type system using open logical relations.


title: Modelling response to trypophobia trigger using intermediate layers of ImageNet networks

abstract: In this paper, we approach the problem of detecting trypophobia triggers using Convolutional neural networks. We show that standard architectures such as VGG or ResNet are capable of recognizing trypophobia patterns. We also conduct experiments to analyze the nature of this phenomenon. To do that, we dissect the network decreasing the number of its layers and parameters. We prove, that even significantly reduced networks have accuracy above 91% and focus their attention on the trypophobia patterns as presented on the visual explanations.


title: Entrywise convergence of iterative methods for eigenproblems

abstract: Several problems in machine learning, statistics, and other fields rely on computing eigenvectors. For large scale problems, the computation of these eigenvectors is typically performed via iterative schemes such as subspace iteration or Krylov methods. While there is classical and comprehensive analysis for subspace convergence guarantees with respect to the spectral norm, in many modern applications other notions of subspace distance are more appropriate. Recent theoretical work has focused on perturbations of subspaces measured in the $\ell_{2 \to \infty}$ norm, but does not consider the actual computation of eigenvectors. Here we address the convergence of subspace iteration when distances are measured in the $\ell_{2 \to \infty}$ norm and provide deterministic bounds. We complement our analysis with a practical stopping criterion and demonstrate its applicability via numerical experiments. Our results show that one can get comparable performance on downstream tasks while requiring fewer iterations, thereby saving substantial computational time.


title: Equal Risk Pricing of Derivatives with Deep Hedging

abstract: This article presents a deep reinforcement learning approach to price and hedge financial derivatives. This approach extends the work of Guo and Zhu (2017) who recently introduced the equal risk pricing framework, where the price of a contingent claim is determined by equating the optimally hedged residual risk exposure associated respectively with the long and short positions in the derivative. Modifications to the latter scheme are considered to circumvent theoretical pitfalls associated with the original approach. Derivative prices obtained through this modified approach are shown to be arbitrage-free. The current paper also presents a general and tractable implementation for the equal risk pricing framework inspired by the deep hedging algorithm of Buehler et al. (2019). An $\epsilon$-completeness measure allowing for the quantification of the residual hedging risk associated with a derivative is also proposed. The latter measure generalizes the one presented in Bertsimas et al. (2001) based on the quadratic penalty. Monte Carlo simulations are performed under a large variety of market dynamics to demonstrate the practicability of our approach, to perform benchmarking with respect to traditional methods and to conduct sensitivity analyses.


title: Stochastic Regret Minimization in Extensive-Form Games

abstract: Monte-Carlo counterfactual regret minimization (MCCFR) is the state-of-the-art algorithm for solving sequential games that are too large for full tree traversals. It works by using gradient estimates that can be computed via sampling. However, stochastic methods for sequential games have not been investigated extensively beyond MCCFR. In this paper we develop a new framework for developing stochastic regret minimization methods. This framework allows us to use any regret-minimization algorithm, coupled with any gradient estimator. The MCCFR algorithm can be analyzed as a special case of our framework, and this analysis leads to significantly-stronger theoretical on convergence, while simultaneously yielding a simplified proof. Our framework allows us to instantiate several new stochastic methods for solving sequential games. We show extensive experiments on three games, where some variants of our methods outperform MCCFR.


title: A quasi-static particle-in-cell algorithm based on an azimuthal Fourier decomposition for highly efficient simulations of plasma-based acceleration: QPAD

abstract: The 3D quasi-static particle-in-cell (PIC) algorithm is a very efficient method for modeling short-pulse laser or relativistic charged particle beam-plasma interactions. In this algorithm, the plasma response to a non-evolving laser or particle beam is calculated using Maxwell's equations based on the quasi-static approximate equations that exclude radiation. The plasma fields are then used to advance the laser or beam forward using a large time step. The algorithm is many orders of magnitude faster than a 3D fully explicit relativistic electromagnetic PIC algorithm. It has been shown to be capable to accurately model the evolution of lasers and particle beams in a variety of scenarios. At the same time, an algorithm in which the fields, currents and Maxwell equations are decomposed into azimuthal harmonics has been shown to reduce the complexity of a 3D explicit PIC algorithm to that of a 2D algorithm when the expansion is truncated while maintaining accuracy for problems with near azimuthal symmetry. This hybrid algorithm uses a PIC description in r-z and a gridless description in $\phi$. We describe a novel method that combines the quasi-static and hybrid PIC methods. This algorithm expands the fields, charge and current density into azimuthal harmonics. A set of the quasi-static field equations are derived for each harmonic. The complex amplitudes of the fields are then solved using the finite difference method. The beam and plasma particles are advanced in Cartesian coordinates using the total fields. Details on how this algorithm was implemented using a similar workflow to an existing quasi-static code, QuickPIC, are presented. The new code is called QPAD for QuickPIC with Azimuthal Decomposition. Benchmarks and comparisons between a fully 3D explicit PIC code, a full 3D quasi-static code, and the new quasi-static PIC code with azimuthal decomposition are also presented.


title: Eccentricity terrain of $\delta$-hyperbolic graphs

abstract: A graph $G=(V,E)$ is $\delta$-hyperbolic if for any four vertices $u,v,w,x$, the two larger of the three distance sums $d(u,v)+d(w,x)$, $d(u,w)+d(v,x)$, and $d(u,x)+d(v,w)$ differ by at most $2\delta \geq 0$. Recent work shows that many real-world graphs have small hyperbolicity $\delta$. This paper describes the eccentricity terrain of a $\delta$-hyperbolic graph. The eccentricity function $e_G(v)=\max\{d(v,u) : u \in V\}$ partitions the vertex set of $G$ into eccentricity layers $C_{k}(G) = \{v \in V : e(v)=rad(G)+k\}$, $k \in \mathbb{N}$, where $rad(G)=\min\{e_G(v): v\in V\}$ is the radius of $G$. The paper studies the eccentricity layers of vertices along shortest paths, identifying such terrain features as hills, plains, valleys, terraces, and plateaus. It introduces the notion of $\beta$-pseudoconvexity, which implies Gromov's $\epsilon$-quasiconvexity, and illustrates the abundance of pseudoconvex sets in $\delta$-hyperbolic graphs. In particular, it shows that all sets $C_{\leq k}(G)=\{v\in V : e_G(v) \leq rad(G) + k\}$, $k\in \mathbb{N}$, are $(2\delta-1)$-pseudoconvex. Additionally, several bounds on the eccentricity of a vertex are obtained which yield a few approaches to efficiently approximating all eccentricities. An $O(\delta |E|)$ time eccentricity approximation $\hat{e}(v)$, for all $v\in V$, is presented that uses distances to two mutually distant vertices and satisfies $e_G(v)-2\delta \leq \hat{e}(v) \leq {e_G}(v)$. It also shows existence of two eccentricity approximating spanning trees $T$, one constructible in $O(\delta |E|)$ time and the other in $O(|E|)$ time, which satisfy ${e}_G(v) \leq e_T(v) \leq {e}_G(v)+4\delta+1$ and ${e}_G(v) \leq e_T(v) \leq {e}_G(v)+6\delta$, respectively. Thus, the eccentricity terrain of a tree gives a good approximation (up-to an additive error $O(\delta))$ of the eccentricity terrain of a $\delta$-hyperbolic graph.


title: Brownian absolute continuity of the KPZ fixed point with arbitrary initial condition

abstract: We show that the law of the KPZ fixed point starting from arbitrary initial condition is absolutely continuous with respect to the law of Brownian motion $B$ on every compact interval. In particular, the Airy$_1$ process is absolutely continuous with respect to $B$ on any compact interval.


title: Quantum phase estimation for a class of generalized eigenvalue problems

abstract: Quantum phase estimation provides a path to quantum computation of solutions to Hermitian eigenvalue problems $Hv = \lambda v$, such as those occurring in quantum chemistry. It is natural to ask whether the same technique can be applied to generalized eigenvalue problems $Av = \lambda B v$, which arise in many areas of science and engineering. We answer this question affirmatively. A restricted class of generalized eigenvalue problems could be solved as efficiently as standard eigenvalue problems. A paradigmatic example is provided by Sturm--Liouville problems. Another example comes from linear ideal magnetohydrodynamics, where phase estimation could be used to determine the stability of magnetically confined plasmas in fusion reactors.


title: Space Efficient Deterministic Approximation of String Measures

abstract: We study approximation algorithms for the following three string measures that are widely used in practice: edit distance (ED), longest common subsequence (LCS), and longest increasing sequence (LIS). All three problems can be solved exactly by standard algorithms that run in polynomial time with roughly $\Theta(n)$ space, where $n$ is the input length, and our goal is to design deterministic approximation algorithms that run in polynomial time with significantly smaller space. Towards this, we design several algorithms that achieve $1+\epsilon$ or $1-\epsilon$ approximation for all three problems, where $\epsilon>0$ can be any constant and even slightly sub constant. Our algorithms are flexible and can be adjusted to achieve the following two regimes of parameters: 1) space $n^\delta$ for any constant $\delta>0$ with running time essentially the same as or slightly more than the standard algorithms; and 2) space $\mathsf{polylog}(n)$ with (a larger) polynomial running time, which puts the approximation versions of the three problems in Steve's class (\textbf{SC}). Our algorithms significantly improve previous results in terms of space complexity, where all known results need to use space at least $\Omega(\sqrt{n})$. Some of our algorithms can also be adapted to work in the asymmetric streaming model [SS13], and output the corresponding sequence. Furthermore, our results can be used to improve a recent result by Farhadi et. al. [FHRS20] about approximating ED in the asymmetric streaming model, reducing the running time from being exponential in [FHRS20] to a polynomial. Our algorithms are based on the idea of using recursion as in Savitch's theorem [Sav70], and a careful adaption of previous techniques to make the recursion work. Along the way we also give a new logspace reduction from longest common subsequence to longest increasing sequence, which may be of independent interest.


title: On the generation and propagation of solitary waves in integrable and non-integrable nonlinear lattices

abstract: We investigate the generation and propagation of solitary waves in the context of the Hertz chain and Toda lattice, with the aim to highlight the similarities, as well as differences between these systems. We begin by discussing the kinetic and potential energy of a solitary wave in these systems, and show that under certain circumstances the kinetic and potential energy profiles in these systems (i.e., their spatial distribution) look reasonably close to each other. While this and other features, such as the connection between the amplitude and the total energy of the wave bear similarities between the two models, there are also notable differences, such as the width of the wave. We then study the dynamical behavior of these systems in response to an initial velocity impulse. For the Toda lattice, we do so by employing the inverse scattering transform, and we obtain analytically the ratio between the energy of the resulting solitary wave and the energy of the impulse, as a function of the impulse velocity; we then compare the dynamics of the Toda system to that of the Hertz system, for which the corresponding quantities are obtained through numerical simulations. In the latter system, we obtain a universality in the fraction of the energy stored in the resulting solitary traveling wave irrespectively of the size of the impulse. This fraction turns out to only depend on the nonlinear exponent. Finally, we investigate the relation between the velocity of the resulting solitary wave and the velocity of the impulse. In particular, we provide an alternative proof for the numerical scaling rule of Hertz type systems.


title: Processing topical queries on images of historical newspaper pages

abstract: Historical newspapers are a source of research for the human and social sciences. However, these image collections are difficult to read by machine due to the low quality of the print, the lack of standardization of the pages in addition to the low quality photograph of some files. This paper presents the processing model of a topic navigation system in historical newspaper page images. The general procedure consists of four modules which are: segmentation of text sub-images and text extraction, preprocessing and representation, induced topic extraction and representation, and document viewing and retrieval interface. The algorithmic and technological approaches of each module are described and the initial test results about a collection covering a range of 28 years are presented.


title: Simple post-translational circadian clock models from selective sequestration

abstract: It is possible that there are post-translational circadian oscillators that continue functioning in the absence of negative feedback transcriptional repression in many cell types from diverse organisms. Apart from the KaiABC system from cyanobacteria, the molecular components and interactions required to create in-vitro ("test-tube") circadian oscillations in different cell types are currently unknown. Inspired by the KaiABC system, I provide "proof-of-principle" mathematical models that a protein with 2 (or more) modification sites which selectively sequesters an effector/cofactor molecule can function as a circadian time-keeper. The 2-site mechanism can be implemented using two relatively simple coupled non-linear ODEs in terms of site occupancy; the models do not require overly special fine-tuning of parameters for generating stable limit cycle oscillations.


title: Adinkra Foundation of Component Decomposition and the Scan for Superconformal Multiplets in 11D, N = 1 Superspace

abstract: For the first time in the physics literature, the Lorentz representations of all 2,147,483,648 bosonic degrees of freedom and 2,147,483,648 fermionic degrees of freedom in an unconstrained eleven dimensional scalar superfield are presented. Comparisons of the conceptual bases for this advance in terms of component field, superfield, and adinkra arguments, respectively, are made. These highlight the computational efficiency of the adinkra-based approach over the others. It is noted at level sixteen in the 11D, N = 1 scalar superfield, the {65} representation of SO(1,10), the conformal graviton, is present. Thus, Adinkra-based arguments suggest the surprising possibility that the 11D, N = 1 scalar superfield alone might describe a Poincare supergravity prepotential in analogy to one of the off-shell versions of 4D, N = 1 superfield supergravity.


title: The metric dimension of critical Galton-Watson trees and linear preferential attachment trees

abstract: The metric dimension of a graph $G$ is the minimal size of a subset $R$ of vertices of $G$ that, upon reporting their graph distance from a distingished (source) vertex $v^\star$, enable unique identification of the source vertex $v^\star$ among all possible vertices of $G$. In this paper we show a Law of Large Numbers (LLN) for the metric dimension of some classes of trees: critical Galton-Watson trees conditioned to have size $n$, and growing general linear preferential attachment trees. The former class includes uniform random trees, the latter class includes Yule-trees (also called random recursive trees), $m$-ary increasing trees, binary search trees, and positive linear preferential attachment trees. In all these cases, we are able to identify the limiting constant in the LLN explicitly. Our result relies on the insight that the metric dimension can be related to subtree properties, and hence we can make use of the powerful fringe-tree literature developed by Aldous and Janson et al.


title: Minimal rational curves on the moduli spaces of symplectic and orthogonal bundles

abstract: Let $C$ be an algebraic curve of genus $g$ and $L$ a line bundle over $C$. Let $\mathcal{MS}_C(n,L)$ and $\mathcal{MO}_C(n,L)$ be the moduli spaces of $L$-valued symplectic and orthogonal bundles respectively, over $C$ of rank $n$. We construct rational curves on these moduli spaces which generalize Hecke curves on the moduli space of vector bundles. As a main result, we show that these Hecke type curves have the minimal degree among the rational curves passing through a general point of the moduli spaces. As its byproducts, we show the non-abelian Torelli theorem and compute the automorphism group of moduli spaces.


title: A Bayes Factor Approach with Informative Prior for Rare Genetic Variant Analysis from Next Generation Sequencing Data

abstract: The discovery of rare genetic variants through Next Generation Sequencing is a very challenging issue in the field of human genetics. We propose a novel region-based statistical approach based on a Bayes Factor (BF) to assess evidence of association between a set of rare variants (RVs) located on the same genomic region and a disease outcome in the context of case-control design. Marginal likelihoods are computed under the null and alternative hypotheses assuming a binomial distribution for the RV count in the region and a beta or mixture of Dirac and beta prior distribution for the probability of RV. We derive the theoretical null distribution of the BF under our prior setting and show that a Bayesian control of the False Discovery Rate (BFDR) can be obtained for genome-wide inference. Informative priors are introduced using prior evidence of association from a Kolmogorov-Smirnov test statistic. We use our simulation program, sim1000G, to generate RV data similar to the 1,000 genomes sequencing project. Our simulation studies showed that the new BF statistic outperforms standard methods (SKAT, SKAT-O, Burden test) in case-control studies with moderate sample sizes and is equivalent to them under large sample size scenarios. Our real data application to a lung cancer case-control study found enrichment for RVs in known and novel cancer genes. It also suggests that using the BF with informative prior improves the overall gene discovery compared to the BF with non-informative prior.


title: Causal Inference under Networked Interference and Intervention Policy Enhancement

abstract: Estimating individual treatment effects from data of randomized experiments is a critical task in causal inference. The Stable Unit Treatment Value Assumption (SUTVA) is usually made in causal inference. However, interference can introduce bias when the assigned treatment on one unit affects the potential outcomes of the neighboring units. This interference phenomenon is known as spillover effect in economics or peer effect in social science. Usually, in randomized experiments or observational studies with interconnected units, one can only observe treatment responses under interference. Hence, how to estimate the superimposed causal effect and recover the individual treatment effect in the presence of interference becomes a challenging task in causal inference. In this work, we study causal effect estimation under general network interference using GNNs, which are powerful tools for capturing the dependency in the graph. After deriving causal effect estimators, we further study intervention policy improvement on the graph under capacity constraint. We give policy regret bounds under network interference and treatment capacity constraint. Furthermore, a heuristic graph structure-dependent error bound for GNN-based causal estimators is provided.


title: Adiabatic quantum pumping in buckled graphene nanoribbon driven by a kink

abstract: We propose a new type of quantum pump in buckled graphene nanoribbon, adiabatically driven by a kink moving along the ribbon. From a practical point of view, pumps with moving scatterers present advantages as compared to gate-driven pumps, like enhanced charge transfer per cycle per channel. The kink geometry is simplified by truncating the spatial arrangement of carbon atoms with the classical $\phi^4$ model solution, including a width renormalization following from the Su-Schrieffer-Heeger model for carbon nanostructures. We demonstrate numerically, as a proof of concept, that for moderate deformations a stationary kink at the ribbon center may block the current driven by the external voltage bias. In the absence of a bias, a moving kink leads to highly-effective pump effect, with a charge per unit cycle dependent on the gate voltage.


title: How far actually is the Galactic Center IRS 13E3 from Sagittarius A$^{\ast}$?

abstract: The Galactic Center IRS 13E cluster is a very intriguing IR object located at $\sim0.13$ pc from Sagittarius A$^\ast$ (Sgr A$^\ast$) in projection distance. There are both arguments for and against the hypothesis that a dark mass like an intermediate mass black hole (IMBH) exists in the cluster. Recently we have detected the rotating ionized gas ring around IRS 13E3, which belongs to the cluster, in the H30$\alpha$ recombination line using ALMA. The enclosed mass is derived to be $M_{\mathrm{encl.}}\simeq2\times10^4$ $M_\odot$, which agrees with an IMBH and is barely less than the astrometric upper limit mass of the IMBH around Sgr A$^\ast$. Because the limit mass depends on the true three-deminsional (3D) distance from Sgr A$^\ast$, it is very important to determine it observationally. However, the 3D distance is indefinite because it is hard to determine the line-of-sight (LOS) distance by usual methods. We would attempt to estimate the LOS distance by spectroscopic informations. The CH$_3$OH molecule is easily destroyed by cosmic ray around Sgr A$^{\ast}$. However, we detected a highly excited CH$_3$OH emission line in the ionized gas stream associated with IRS 13E3. This indicates that IRS 13E3 is located at $r\gtrsim 0.4$ pc from Sgr A$^{\ast}$.


title: Twisting Noncommutative Geometries with Applications to High Energy Physics

abstract: With the bare essentials of noncommutative geometry (defined by a spectral triple), we first describe how it naturally gives rise to gauge theories. Then, we quickly review the notion of twisting (in particular, minimally) noncommutative geometries and how it induces a Wick rotation, that is, a transition of the metric signature from euclidean to Lorentzian. We focus on comparatively more tractable examples of spectral triples; such as the ones corresponding to a closed Riemannian spin manifold, $U(1)$ gauge theory, and electrodynamics. By minimally twisting these examples and computing their associated fermionic actions, we demonstrate how to arrive at physically relevant actions (such as the Weyl and Dirac actions) in Lorentz signature, even though starting from euclidean spectral triples. In the process, not only do we extract a physical interpretation of the twist, but we also capture exactly how the Wick rotation takes place at the level of the fermionic action.


title: Expressing Objects just like Words: Recurrent Visual Embedding for Image-Text Matching

abstract: Existing image-text matching approaches typically infer the similarity of an image-text pair by capturing and aggregating the affinities between the text and each independent object of the image. However, they ignore the connections between the objects that are semantically related. These objects may collectively determine whether the image corresponds to a text or not. To address this problem, we propose a Dual Path Recurrent Neural Network (DP-RNN) which processes images and sentences symmetrically by recurrent neural networks (RNN). In particular, given an input image-text pair, our model reorders the image objects based on the positions of their most related words in the text. In the same way as extracting the hidden features from word embeddings, the model leverages RNN to extract high-level object features from the reordered object inputs. We validate that the high-level object features contain useful joint information of semantically related objects, which benefit the retrieval task. To compute the image-text similarity, we incorporate a Multi-attention Cross Matching Model into DP-RNN. It aggregates the affinity between objects and words with cross-modality guided attention and self-attention. Our model achieves the state-of-the-art performance on Flickr30K dataset and competitive performance on MS-COCO dataset. Extensive experiments demonstrate the effectiveness of our model.


title: Tuning of LANSCE 805-MHz High-Energy Linear Accelerator with Reduced Beam Losses

abstract: Suppression of beam losses is essential for successful operation of high-intensity linac. Historically, the values of the field amplitudes and phases of the side-coupled, 805-MHz LANSCE linac modules are maintained using a well-known delta-t tuning procedure. Transverse matching of the beam with accelerator is performed through adjustments of beam ellipses in 4D phase space with accelerator lattice using matching quadrupoles. Control of the beam-energy ramp along the length of a proton linear accelerator is required to improve tune of accelerator and decrease beam losses. Time-of-flight measurements of the H- beam energy are now being used to confirm and improve the overall control of the energy ramp along the linac. The time-of-flight method utilizes absolute measurements of beam energy using direct signals from beam at an oscilloscope, as well as the difference in RF phases measured as the beam passes installed delta-t pickup loops. A newly developed BPPM data acquisition system is used. Beam energy measurement along accelerator together with phase scans, klystrons output power control, and delta-t method, allow tuning of accelerator much more close to original design and reduce losses generated by linac. Details of the upgraded tuning procedure and results of tuning and operation are presented.


title: The Problem with Metrics is a Fundamental Problem for AI

abstract: Optimizing a given metric is a central aspect of most current AI approaches, yet overemphasizing metrics leads to manipulation, gaming, a myopic focus on short-term goals, and other unexpected negative consequences. This poses a fundamental contradiction for AI development. Through a series of real-world case studies, we look at various aspects of where metrics go wrong in practice and aspects of how our online environment and current business practices are exacerbating these failures. Finally, we propose a framework towards mitigating the harms caused by overemphasis of metrics within AI by: (1) using a slate of metrics to get a fuller and more nuanced picture, (2) combining metrics with qualitative accounts, and (3) involving a range of stakeholders, including those who will be most impacted.


title: A Trust-Region Method For Nonsmooth Nonconvex Optimization

abstract: We propose a trust-region type method for a class of nonsmooth nonconvex optimization problems where the objective function is a summation of a (probably nonconvex) smooth function and a (probably nonsmooth) convex function. The model function of our trust-region subproblem is always quadratic and the linear term of the model is generated using abstract descent directions. Therefore, the trust-region subproblems can be easily constructed as well as efficiently solved by cheap and standard methods. When the accuracy of the model function at the solution of the subproblem is not sufficient, we add a safeguard on the stepsizes for improving the accuracy. For a class of functions that can be "truncated", an additional truncation step is defined and a stepsize modification strategy is designed. The overall scheme converges globally and we establish fast local convergence under suitable assumptions. In particular, using a connection with a smooth Riemannian trust-region method, we prove local quadratic convergence for partly smooth functions under a strict complementary condition. Preliminary numerical results on a family of $\ell_1$-optimization problems are reported and demonstrate the efficiency of our approach.


title: Queueing Subject To Action-Dependent Server Performance: Utilization Rate Reduction

abstract: We consider a discrete-time system comprising a first-come-first-served queue, a non-preemptive server, and a stationary non-work-conserving scheduler. New tasks enter the queue according to a Bernoulli process with a pre-specified arrival rate. At each instant, the server is either busy working on a task or is available. When the server is available, the scheduler either assigns a new task to the server or allows it to remain available (to rest). In addition to the aforementioned availability state, we assume that the server has an integer-valued activity state. The activity state is non-decreasing during work periods, and is non-increasing otherwise. In a typical application of our framework, the server performance (understood as task completion probability) worsens as the activity state increases. In this article, we build on and transcend recent stabilizability results obtained for the same framework. Specifically, we establish methods to design scheduling policies that not only stabilize the queue but also reduce the utilization rate - understood as the infinite-horizon time-averaged portion of time the server is working. This article has a main theorem leading to two key results: (i) We put forth a tractable method to determine, using a finite-dimensional linear program (LP), the infimum of all utilization rates that can be achieved by scheduling policies that are stabilizing, for a given arrival rate. (ii) We propose a design method, also based on finite-dimensional LPs, to obtain stabilizing scheduling policies that can attain a utilization rate arbitrarily close to the aforementioned infimum. We also establish structural and distributional convergence properties, which are used throughout the article, and are significant in their own right.


title: Evaporation Induced Rayleigh-Taylor Instabilities in Polymer Solutions

abstract: Understanding the mechanics of detrimental convective instabilities in drying polymer solutions is crucial in many applications such as the production of film coatings. It is well known that solvent evaporation in polymer solutions can lead to Rayleigh-B\'enard or Marangoni-type instabilities. Here we reveal another mechanism, namely that evaporation can cause the interface to display Rayleigh-Taylor instabilities due to the build-up of a dense layer at the air-liquid interface. We study experimentally the onset time ($t_p$) of the instability as a function of the macroscopic properties of aqueous polymer solutions, which we tune by varying the polymer concentration ($c_0$), molecular weight and polymer type. In dilute solutions, $t_p$ shows two limiting behaviors depending on the polymer diffusivity. For high diffusivity polymers (low molecular weight), the pluming time scales as $c_0^{-2/3}$. This result agrees with previous studies on gravitational instabilities in miscible systems where diffusion stabilizes the system. On the other hand, in low diffusivity polymers the pluming time scales as $c_0^{-1}$. The stabilizing effect of an effective interfacial tension, similar to those in immiscible systems, explains this strong concentration dependence. Above a critical concentration, $\hat{c}$, viscosity delays the growth of the instability, allowing time for diffusion to act as the dominant stabilizing mechanism. This results in $t_p$ scaling as $(\nu/c_0)^{2/3}$.


title: The CGM at Cosmic Noon with KCWI: Outflows from a Star-forming Galaxy at z=2.071

abstract: We present the first results from our CGM at Cosmic Noon with KCWI program to study gas flows in the circumgalactic medium (CGM) at $z=2-3$. Combining the power of a high-resolution VLT/UVES quasar spectrum, an HST/ACS image, and integral field spectroscopy with Keck/KCWI, we detected Lya emission from a $1.7L_{\ast}$ galaxy at $z_{\rm gal}=2.0711$ associated with a Lyman limit system with weak MgII ($W_r(2796)=0.24$ Ang) in quasar field J143040$+$014939. The galaxy is star-forming (${\rm SFR}_{\rm FUV}=37.8$ M$_{\odot}$ yr$^{-1}$) and clumpy: either an edge-on disk ($i=85^{\circ}$) or, less likely, a major merger. The background quasar probes the galaxy at an impact parameter of $D=66$ kpc along the projected galaxy minor axis ($\Phi=89^{\circ}$). From photoionization modeling of the absorption system, we infer a total line-of-sight CGM metallicity of ${\rm [Si/H]}=-1.5^{+0.4}_{-0.3}$. The absorption system is roughly kinematically symmetric about $z_{\rm gal}$, with a full MgII velocity spread of $\sim210$ km s$^{-1}$. Given the galaxy-quasar orientation, CGM metallicity, and gas kinematics, we interpret this gas as an outflow that has likely swept-up additional material. By modeling the absorption as a polar outflow cone, we find the gas is decelerating with average radial velocity $V_{\rm out}=109-588$ km s$^{-1}$ for half opening angles of $\theta_0=14^{\circ}-75^{\circ}$. Assuming a constant $V_{\rm out}$, it would take on average $t_{\rm out}\sim111-597$ Myr for the gas to reach 66 kpc. The outflow is energetic, with a mass outflow rate of $\dot{M}_{\rm out}<52{\pm37}$ M$_{\odot}$ yr$^{-1}$ and mass loading factor of $\eta<1.4{\pm1.0}$. We aim to build a sample of $\sim50$ MgII absorber--galaxy pairs at this epoch to better understand gas flows when they are most actively building galaxies.