title: Holographic modeling of nuclear matter and neutron stars

abstract: I review holographic models for (dense and cold) nuclear matter, neutron stars, and their mergers. I start by a brief general discussion on current knowledge of cold QCD matter and neutron stars, and go on discussing various approaches to model cold nuclear and quark matter by using gauge/gravity duality, pointing out their strengths and weaknesses. Then I focus on recent results for a complex bottom-up holographic framework (V-QCD), which also takes input from lattice QCD results, effective field theory, and perturbative QCD. Dense nuclear matter is modeled in V-QCD through a homogeneous non-Abelian bulk gauge field. Feasible "hybrid" equations of state for cold nuclear (and quark) matter can be constructed by using traditional methods (e.g., effective field theory) at low densities and the holographic V-QCD model at higher densities. I discuss the constraints from this approach to the properties of the nuclear to quark matter transition as well as to properties of neutron stars. Using such hybrid equations of state as an input for numerical simulations of neutron star mergers, I also derive predictions for the spectrum of produced gravitational waves.


title: Diffusion of cosmic rays in MHD turbulence

abstract: We review some recent findings on diffusion of cosmic rays (CRs) in magnetohydrodynamic (MHD) turbulence obtained by adopting the numerically-tested model of MHD turbulence, including perpendicular superdiffusion of CRs, inefficient gyroresonant scattering by Alfv\'{e}n and slow modes with scale-dependent turbulence anisotropy, resonance-broadened Transit Time Damping (TTD) interaction, and mirror diffusion. As the diffusion behavior of CRs strongly depends on the properties of MHD turbulence, theoretical modeling of CR diffusion, its numerical testing, and interpretation of CR-related observations require proper modeling of MHD turbulence.


title: Effects of Quantum Pair Creation and Annihilation on a Classical Exclusion Process: the transverse XY model with TASEP

abstract: We investigate how particle pair creation and annihilation, within the quantum transverse XY model, affects the non-equilibrium steady state (NESS) and Liouvillian gap of the stochastic Totally Asymmetric Exclusion Process (TASEP). By utilising operator quantization we formulate a perturbative description of the NESS. Furthermore, we estimate the Liouvillian gap by exploiting a Majorana canonical basis as the basis of super-operators. In this manner we show that the Liouvillian gap can remain finite in the thermodynamic limit provided the XY model anisotropy parameter remains non-zero. Additionally, we show that the character of the gap with respect to the anisotropy parameter differs depending on the phase of the XY model. The change of character corresponds to the quantum phase transition of the XY model.


title: Machine learning assisted GaAsN circular polarimeter

abstract: We demonstrate the application of a two stage machine learning algorithm that enables to correlate the electrical signals from a GaAs$_x$N$_{1-x}$ circular polarimeter with the intensity, degree of circular polarization and handedness of an incident light beam. Specifically, we employ a multimodal logistic regression to discriminate the handedness of light and a 6-layer neural network to establish the relationship between the input voltages, the intensity and degree of circular polarization. We have developed a particular neural network training strategy that substantially improves the accuracy of the device. The algorithm was trained and tested on theoretically generated photoconductivity and on photoluminescence experimental results. Even for a small training experimental dataset (70 instances), it is shown that the proposed algorithm correctly predicts linear, right and left circularly polarized light misclassifying less than $1.5\%$ of the cases and attains an accuracy larger than $97\%$ in the vast majority of the predictions ($92\%$) for intensity and degree of circular polarization. These numbers are significantly improved for the larger theoretically generated datasets (4851 instances). The algorithm is versatile enough that it can be easily adjusted to other device configurations where a map needs to be established between the input parameters and the device response. Training and testing data files as well as the algorithm are provided as supplementary material.


title: Accretion onto black holes inside neutron stars with piecewise-polytropic equations of state: analytic and numerical treatments

abstract: We consider spherically symmetric accretion onto a small, possibly primordial, black hole residing at the center of a neutron star governed by a cold nuclear equation of state (EOS). We generalize the relativistic Bondi solution for such EOSs, approximated by piecewise polytropes, and thereby obtain analytical expressions for the steady-state matter profiles and accretion rates. We compare these rates with those found by time-dependent, general relativistic hydrodynamical simulations upon relaxation and find excellent agreement. We consider several different candidate EOSs, neutron star masses and central densities and find that the accretion rates vary only little, resulting in an accretion rate that depends primarily on the black hole mass, and only weakly on the properties of the neutron star.


title: Imaging the water snowline around protostars with water and HCO$^+$ isotopologues

abstract: The water snowline location in protostellar envelopes provides crucial information about the thermal structure and the mass accretion process as it can inform about the occurrence of recent ($\lesssim$1,000 yr) accretion bursts. In addition, the ability to image water emission makes these sources excellent laboratories to test indirect snowline tracers such as H$^{13}$CO$^+$. We study the water snowline in five protostellar envelopes in Perseus using a suite of molecular line observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA) at $\sim$0.2$^{\prime\prime}-$0.7$^{\prime\prime}$ (60--210 au) resolution. B1-c provides a textbook example of compact H$_2^{18}$O ($3_{1,3}-2_{2,0}$) and HDO ($3_{1,2}-2_{2,1}$) emission surrounded by a ring of H$^{13}$CO$^+$ ($J=2-1$) and HC$^{18}$O$^+$ ($J=3-2$). Compact HDO surrounded by H$^{13}$CO$^+$ is also detected toward B1-bS. The optically thick main isotopologue HCO$^+$ is not suited to trace the snowline and HC$^{18}$O$^+$ is a better tracer than H$^{13}$CO$^+$ due to a lower contribution from the outer envelope. However, since a detailed analysis is needed to derive a snowline location from H$^{13}$CO$^+$ or HC$^{18}$O$^+$ emission, their true value as snowline tracer will lie in the application in sources where water cannot be readily detected. For protostellar envelopes, the most straightforward way to locate the water snowline is through observations of H$_2^{18}$O or HDO. Including all sub-arcsecond resolution water observations from the literature, we derive an average burst interval of $\sim$10,000 yr, but high-resolution water observations of a larger number of protostars is required to better constrain the burst frequency.


title: Predictions for Complex Distributions of Stellar Elemental Abundances in Low-Mass Galaxies

abstract: We investigate stellar elemental abundance patterns at z = 0 in 8 low-mass (M_* = 10^6 - 10^9 M_sun) galaxies in the Feedback in Realistic Environments (FIRE-2) cosmological simulations. Using magnesium (Mg) as a representative alpha-element, we explore stellar abundance patterns in [Mg/Fe] versus [Fe/H], which follow an overall monotonic trend that evolved slowly over time. Beyond this, we explore 3 notable secondary features in enrichment (found in three different case-study galaxies) that arise from a galaxy merger or bursty star formation. First, we observe a secondary track with a lower [Mg/Fe] than the main trend. At z = 0, stars from this track are predominantly found within 2-6 kpc of the center; they were accreted in a 1:3 total-mass-ratio merger ~ 0.4 Gyr ago. Second, we find a distinct elemental bi-modality that forms following a strong burst in star formation in a galaxy at t_lookback ~ 10 Gyr. This burst quenched star formation for ~ 0.66 Gyr, allowing Ia supernovae to enrich the system with iron before star formation resumed. Third, we examine stripes in enrichment that run roughly orthogonal to the dominant [Mg/Fe] versus [Fe/H] trend; these stripes correspond to short bursts of star formation during which core-collapse supernovae enrich the surrounding medium with Mg (and Fe) on short timescales. If observed, these features would substantiate the utility of elemental abundances in revealing the assembly and star formation histories of dwarf galaxies. We explore the observability of these features for upcoming spectroscopic studies. Our results show that precise measurements of elemental abundance patterns can reveal critical events in the formation histories of low-mass galaxies.


title: Explaining deep learning of galaxy morphology with saliency mapping

abstract: We successfully demonstrate the use of explainable artificial intelligence (XAI) techniques on astronomical datasets in the context of measuring galactic bar lengths. The method consists of training convolutional neural networks on human classified data from Galaxy Zoo in order to predict general galaxy morphologies, and then using SmoothGrad (a saliency mapping technique) to extract the bar for measurement by a bespoke algorithm. We contrast this to another method of using a convolutional neural network to directly predict galaxy bar lengths. These methods achieved correlation coefficients of 0.76 and 0.59, and root mean squared errors of 1.69 and 2.10 respective to human measurements. We conclude that XAI methods outperform conventional deep learning in this case, which could be reasonably explained by the larger datasets available when training the models. We suggest that our XAI method can be used to extract other galactic features (such as the bulge-to-disk ratio) without needing to collect new datasets or train new models. We also suggest that these techniques can be used to refine deep learning models as well as identify and eliminate bias within training datasets.


title: Growth modes of partially fluorinated organic molecules on amorphous silicon dioxide

abstract: We study the influence of fluorination on nucleation and growth of the organic para-sexiphenyl molecule (p-6P) on amorphous silicon dioxide ($\alpha$-SiO$_2$) by means of atomistically resolved classical molecular dynamics computer simulations. We use a simulation model that mimics the experimental deposition from the vapor and subsequent self-assembly onto the underlying surface. Our model reproduces the experimentally observed orientational changes from lying to upright standing configurations of the grown layers. We demonstrate that the increase in the number of fluorinated groups inside the p-6P leads to a smoother, layer-by-layer growth on the $\alpha$-SiO$_2$ surface: We observe that in the first layers, due to strong molecule-substrate interactions the molecules first grow in chiral (fan-like) structures, where each consecutive molecule has a higher angle, supported by molecules lying underneath. Subsequently deposited molecules bind to the already standing molecules of the chiral structures until all molecules are standing. The growth of chiral islands is the main mechanism for growth of the fluorinated p-6P derivative, while the p-6P, due to the lower interaction with the underlying substrate, forms less chiral structures. This leads to a lower energy barrier for step-edge crossing for the fluorinated molecules. We find that partial fluorination of the p-6P molecule can in this way significantly alter its growth behaviour by modifying the rough, 3D growth into a smooth, layer-by-layer growth. This has implications for the rational design of molecules and their functionalized forms which could be tailored for a desired growth behavior and structure formation.


title: Boundedness from below of SU(n) potentials

abstract: Vacuum stability requires that the scalar potential is bounded from below. Whether or not this is true depends on the scalar quartic interactions alone, but even so the analysis is arduous and has only been carried out for a limited set of models. Complementing the existing literature, this work contains the necessary and sufficient conditions for two SU(n) invariant potentials to be bounded from below. In particular, expressions are given for models with the fundamental and the 2-index (anti)symmetric representations of this group. A sufficient condition for vacuum stability is also provided for models with the fundamental and the adjoint representations. Finally, some considerations are made concerning the model with the gauge group SU(2) and the scalar representations $\boldsymbol{1}$, $\boldsymbol{2}$ and $\boldsymbol{3}$; such a setup is particularly important for neutrino mass generation and lepton number violation.


title: Prospects of fast flavor neutrino conversion in rotating core-collapse supernovae

abstract: There is mounting evidence that neutrinos undergo fast flavor conversion (FFC) in core-collapse supernova (CCSN). In this letter, we investigate the roles of stellar rotation on the occurrence of FFC by carrying out axisymmetric CCSN simulations with full Boltzmann neutrino transport. Our result suggests that electron neutrino lepton number (ELN) angular crossings, which are the necessary and sufficient condition to trigger FFC, preferably occur in the equatorial region for rotating CCSNe. By scrutinizing the neutrino--matter interaction and neutrino radiation field, we find some pieces of evidence that the stellar rotation facilitates the occurrence of FFC. The low-electron-fraction region in the post-shock layer expands by centrifugal force, enhancing the disparity of neutrino absorption between electron-type neutrinos ($\nu_{\rm e}$) and their anti-particles ($\bar{\nu}_{\rm e}$). This has a significant impact on the angular distribution of neutrinos in momentum space, in which $\nu_{\rm e}$ tends to be more isotropic than $\bar{\nu}_{\rm e}$; consequently, ELN crossings emerge. The ELN crossing found in this study is clearly associated with rotation, which motivates further investigation on how the subsequent FFC influences explosion dynamics, nucleosynthesis, and neutrino signals in rotating CCSNe.


title: Hidden quasiconservation laws in fracton hydrodynamics

abstract: We show that the simplest universality classes of fracton hydrodynamics in more than one spatial dimension, including isotropic theories of charge and dipole conservation, can exhibit hidden "quasiconservation laws", in which certain higher multipole moments can only decay due to dangerously irrelevant corrections to hydrodynamics. We present two simple examples of this phenomenon. Firstly, an isotropic dipole-conserving fluid in the infinite plane conserves an infinite number of "harmonic multipole charges" within linear response; we calculate the decay or growth of these charges due to dangerously irrelevant nonlinearities. Secondly, we consider a model with $xy$ and $x^2-y^2$ quadrupole conservation, in addition to dipole conservation, which is described by isotropic fourth-order subdiffusion, yet has dangerously irrelevant sixth-order corrections necessary to relax the harmonic multipole charges. We confirm our predictions for the anomalously slow decay of the harmonic conserved charges in each setting by using numerical simulations, both of the nonlinear hydrodynamic differential equations, and in quantum automaton circuits on a square lattice.


title: Mutually unbiased frames

abstract: In this work, the concept of mutually unbiased frames is introduced as the most general notion of unbiasedness for sets composed by linearly independent and normalized vectors. It encompasses the already existing notions of unbiasedness for orthonormal bases, regular simplices, equiangular tight frames, positive operator valued measure, and also includes symmetric informationally complete quantum measurements. After introducing the tool, its power is shown by finding the following results about the last mentioned class of constellations: (i) real fiducial states do not exist in any even dimension, and (ii) unknown $d$-dimensional fiducial states are parameterized, a priori, with roughly $3d/2$ real variables only, without loss of generality. Furthermore, multi-parametric families of pure quantum states having minimum uncertainty with regard to several choices of $d+1$ orthonormal bases are shown, in every dimension $d$. These last families contain all existing fiducial states in every finite dimension, and the bases include maximal sets of $d+1$ mutually unbiased bases, when $d$ is a prime number.


title: Coherence boosting: When your pretrained language model is not paying enough attention

abstract: Long-range semantic coherence remains a challenge in automatic language generation and understanding. We demonstrate that large language models have insufficiently learned the effect of distant words on next-token prediction. We present coherence boosting, an inference procedure that increases a LM's focus on a long context. We show the benefits of coherence boosting with pretrained models by distributional analyses of generated ordinary text and dialog responses. It is also found that coherence boosting with state-of-the-art models for various zero-shot NLP tasks yields performance gains with no additional training.


title: Nonlinear proper orthogonal decomposition for convection-dominated flows

abstract: Autoencoder techniques find increasingly common use in reduced order modeling as a means to create a latent space. This reduced order representation offers a modular data-driven modeling approach for nonlinear dynamical systems when integrated with a time series predictive model. In this letter, we put forth a nonlinear proper orthogonal decomposition (POD) framework, which is an end-to-end Galerkin-free model combining autoencoders with long short-term memory networks for dynamics. By eliminating the projection error due to the truncation of Galerkin models, a key enabler of the proposed nonintrusive approach is the kinematic construction of a nonlinear mapping between the full-rank expansion of the POD coefficients and the latent space where the dynamics evolve. We test our framework for model reduction of a convection-dominated system, which is generally challenging for reduced order models. Our approach not only improves the accuracy, but also significantly reduces the computational cost of training and testing.


title: ASPECTNEWS: Aspect-Oriented Summarization of News Documents

abstract: Generic summaries try to cover an entire document and query-based summaries try to answer document-specific questions. But real users' needs often fall in between these extremes and correspond to aspects, high-level topics discussed among similar types of documents. In this paper, we collect a dataset of realistic aspect-oriented summaries, AspectNews, which covers different subtopics about articles in news sub-domains. We annotate data across two domains of articles, earthquakes and fraud investigations, where each article is annotated with two distinct summaries focusing on different aspects for each domain. A system producing a single generic summary cannot concisely satisfy both aspects. Our focus in evaluation is how well existing techniques can generalize to these domains without seeing in-domain training data, so we turn to techniques to construct synthetic training data that have been used in query-focused summarization work. We compare several training schemes that differ in how strongly keywords are used and how oracle summaries are extracted. Our evaluation shows that our final approach yields (a) focused summaries, better than those from a generic summarization system or from keyword matching; (b) a system sensitive to the choice of keywords.


title: Strong $L^p$-error analysis of nonlinear Monte Carlo approximations for high-dimensional semilinear partial differential equations

abstract: Full-history recursive multilevel Picard (MLP) approximation schemes have been shown to overcome the curse of dimensionality in the numerical approximation of high-dimensional semilinear partial differential equations (PDEs) with general time horizons and Lipschitz continuous nonlinearities. However, each of the error analyses for MLP approximation schemes in the existing literature studies the $L^2$-root-mean-square distance between the exact solution of the PDE under consideration and the considered MLP approximation and none of the error analyses in the existing literature provides an upper bound for the more general $L^p$-distance between the exact solution of the PDE under consideration and the considered MLP approximation. It is the key contribution of this article to extend the $L^2$-error analysis for MLP approximation schemes in the literature to a more general $L^p$-error analysis with $p\in (0,\infty)$. In particular, the main result of this article proves that the proposed MLP approximation scheme indeed overcomes the curse of dimensionality in the numerical approximation of high-dimensional semilinear PDEs with the approximation error measured in the $L^p$-sense with $p \in (0,\infty)$.


title: Non-Euclidean Contraction Analysis of Continuous-Time Neural Networks

abstract: Critical questions in dynamical neuroscience and machine learning are related to the study of continuous-time neural networks and their stability, robustness, and computational efficiency. These properties can be simultaneously established via a contraction analysis. This paper develops a comprehensive non-Euclidean contraction theory for continuous-time neural networks. First, for non-Euclidean $\ell_{1}/\ell_{\infty}$ logarithmic norms, we establish quasiconvexity with respect to positive diagonal weights and closed-form worst-case expressions over certain matrix polytopes. Second, for locally Lipschitz maps (e.g., arising as activation functions), we show that their one-sided Lipschitz constant equals the essential supremum of the logarithmic norm of their Jacobian. Third and final, we apply these general results to classes of continuous-time neural networks, including Hopfield, firing rate, Persidskii, Lur'e and other models. For each model, we compute the optimal contraction rate and corresponding weighted non-Euclidean norm via a linear program or, in some special cases, via a Hurwitz condition on the Metzler majorant of the synaptic matrix. Our non-Euclidean analysis establishes also absolute, connective, and total contraction properties.


title: Characterization of plastic scintillator bars using fast neutrons from D-D and D-T reactions

abstract: We report results of fast neutron response in plastic scintillator (PS) bars from deuterium-deuterium (D-D) and deuterium-tritium (D-T) reactions using Purnima Neutron Generator Facility, BARC, Mumbai. These measurements are useful in context of Indian Scintillator Matrix for Reactor Anti-Neutrino (ISMRAN) detection, an array of 10x10 PS bars, used to measure reactor anti-neutrinos through inverse beta decay (IBD) signal. ISMRAN detector, an above-ground experiment close to the reactor core (~13m), deals with an active fast neutron background inside the reactor hall. A good understanding of fast neutron response in PS bars is an essential pre-requisite for suppression and discrimination of fast neutron background from IBD events. A monoenergetic neutron beam from the fusion reaction of D-D at 2.45 MeV and D-T at 14.1 MeV are used to characterize the energy response in these bars. The neutron energy response function has been simulated using the GEANT4 package and are compared with the measured data. A reasonable agreement of deposited energies by fast neutrons in PS bars between data and simulation are obtained for these reactions. The ratio of energy deposition in adjacent bars is used to discriminate between prompt IBD, fast neutron and neutron capture cascade gamma events.


title: The Dangers of Underclaiming: Reasons for Caution When Reporting How NLP Systems Fail

abstract: Researchers in NLP often frame and discuss research results in ways that serve to deemphasize the field's successes, often in response to the field's widespread hype. Though well-meaning, this has yielded many misleading or false claims about the limits of our best technology. This is a problem, and it may be more serious than it looks: It harms our credibility in ways that can make it harder to mitigate present-day harms, like those involving biased systems for content moderation or resume screening. It also limits our ability to prepare for the potentially enormous impacts of more distant future advances. This paper urges researchers to be careful about these claims and suggests some research directions and communication strategies that will make it easier to avoid or rebut them.


title: Positive energy and non-SUSY flows in ISO(7) gauged supergravity

abstract: We consider maximal gauged supergravity in 4D with ISO(7) gauge group, which arises from a consistent truncation of massive IIA supergravity on a six-sphere. Within its G$_2$ invariant sector, the theory is known to possess a supersymmetric AdS extremum as well as two non-supersymmetric ones. In this context we provide a first order formulation of the theory by making use of the Hamilton-Jacobi (HJ) formalism. This allows us to derive a positive energy theorem for both non-supersymmetric extrema. Subsequently, we also find novel non-supersymmetric domain walls (DW) interpolating between the supersymmetric extremum and each of the other two. Finally, we discuss a perturbative HJ technique that may be used in order to solve for curved DW geometries.


title: Dynamics in Systems with Modulated Symmetries

abstract: We extend the notions of multipole and subsystem symmetries to more general {\it spatially modulated} symmetries. We uncover two instances with exponential and (quasi)-periodic modulations, and provide simple microscopic models in one, two and three dimensions. Seeking to understand their effect in the long-time dynamics, we numerically study a stochastic cellular automaton evolution that obeys such symmetries. We prove that in one dimension, the periodically modulated symmetries lead to a diffusive scaling of correlations modulated by a finite microscopic momentum. In higher dimensions, these symmetries take the form of lines and surfaces of conserved momenta. These give rise to exotic forms of sub-diffusive behavior with a rich spatial structure influenced by lattice-scale features. Exponential modulation, on the other hand, can lead to correlations that are infinitely long-lived at the boundary, while decaying exponentially in the bulk.


title: Minimum Viable Device Drivers for ARM TrustZone

abstract: While TrustZone can isolate IO hardware, it lacks drivers for modern IO devices. Rather than porting drivers, we propose a novel approach to deriving minimum viable drivers: developers exercise a full driver and record the driver/device interactions; the processed recordings, dubbed driverlets, are replayed in the TEE at run time to access IO devices. Driverlets address two key challenges: correctness and expressiveness, for which they build on a key construct called interaction template. The interaction template ensures faithful reproduction of recorded IO jobs (albeit on new IO data); it accepts dynamic input values; it tolerates nondeterministic device behaviors. We demonstrate driverlets on a series of sophisticated devices, making them accessible to TrustZone for the first time to our knowledge. Our experiments show that driverlets are secure, easy to build, and incur acceptable overhead (1.4x -2.7x compared to native drivers). Driverlets fill a critical gap in the TrustZone TEE, realizing its long-promised vision of secure IO.


title: Nonlinear Landauer formula: Nonlinear response theory of disordered and topological materials

abstract: The Landauer formula provides a general scattering formulation of electrical conduction. Despite its utility, it has been mainly applied to the linear-response regime, and a scattering theory of nonlinear response has yet to be fully developed. Here, we extend the Landauer formula to the nonlinear-response regime. We show that while the linear conductance is directly related to the transmission probability, the nonlinear conductance is given by its derivatives with respect to energy. This sensitivity to the energy derivatives is shown to produce unique nonlinear transport phenomena of mesoscopic systems including disordered and topological materials. By way of illustration, we investigate nonlinear conductance of disordered chains and identify their universal behavior according to symmetry. In particular, we find large singular nonlinear conductance for zero modes, including Majorana zero modes in topological superconductors. We also show the critical behavior of nonlinear response around the mobility edges due to the Anderson transitions. Moreover, we study nonlinear response of graphene as a prime example of topological materials featuring quantum anomaly. Furthermore, considering the geometry of electronic wave functions, we develop a scattering theory of the nonlinear Hall effect. We establish a new connection between the nonlinear Hall response and the nonequilibrium quantum fluctuations. We also discuss the influence of disorder and Anderson localization on the nonlinear Hall effect. Our work opens a new avenue in quantum physics beyond the linear-response regime.


title: Universal finite-size amplitude and anomalous entanglement entropy of $z=2$ quantum Lifshitz criticalities in topological chains

abstract: We consider Lifshitz criticalities with dynamical exponent $z=2$ that emerge in a class of topological chains. There, such a criticality plays a fundamental role in describing transitions between symmetry-enriched conformal field theories (CFTs). We report that, at such critical points in one spatial dimension, the finite-size correction to the energy scales with system size, $L$, as $\sim L^{-2}$, with universal and anomalously large coefficient. The behavior originates from the specific dispersion around the Fermi surface, $\epsilon \propto \pm k^2$. We also show that the entanglement entropy exhibits at the criticality a non-logarithmic dependence on $l/L$, where $l$ is the length of the sub-system. In the limit of $l\ll L$, the maximally-entangled ground state has the entropy, $S(l/L)=S_0+2n(l/L)\log(l/L)$. Here $S_0$ is some non-universal entropy originating from short-range correlations and $n$ is half-integer or integer depending on the degrees of freedom in the model. We show that the novel entanglement originates from the long-range correlation mediated by a zero mode in the low energy sector. The work paves the way to study finite-size effects and entanglement entropy around Lifshitz criticalities and offers an insight into transitions between symmetry-enriched criticalities.


title: Memory-augmented Adversarial Autoencoders for Multivariate Time-series Anomaly Detection with Deep Reconstruction and Prediction

abstract: Detecting anomalies for multivariate time-series without manual supervision continues a challenging problem due to the increased scale of dimensions and complexity of today's IT monitoring systems. Recent progress of unsupervised time-series anomaly detection mainly use deep autoencoders to solve this problem, i.e. training on normal samples and producing significant reconstruction error on abnormal inputs. However, in practice, autoencoders can reconstruct anomalies so well, due to powerful capabilites of neural networks. Besides, these approaches can be ineffective for identifying non-point anomalies, e.g. contextual anomalies and collective anomalies, since they solely utilze a point-wise reconstruction objective. To tackle the above issues, we propose MemAAE (\textit{Memory-augmented Adversarial Autoencoders with Deep Reconstruction and Prediction}), a novel unsupervised anomaly detection method for time-series. By jointly training two complementary proxy tasks, reconstruction and prediction, with a shared network architecture, we show that detecting anomalies via multiple tasks obtains superior performance rather than single-task training. Additionally, a compressive memory module is introduced to preserve normal patterns, avoiding unexpected generalization on abnormal inputs. Through extensive experiments, MemAAE achieves an overall F1 score of 0.90 on four public datasets, significantly outperforming the best baseline by 0.02.


title: GrowSpace: Learning How to Shape Plants

abstract: Plants are dynamic systems that are integral to our existence and survival. Plants face environment changes and adapt over time to their surrounding conditions. We argue that plant responses to an environmental stimulus are a good example of a real-world problem that can be approached within a reinforcement learning (RL)framework. With the objective of controlling a plant by moving the light source, we propose GrowSpace, as a new RL benchmark. The back-end of the simulator is implemented using the Space Colonisation Algorithm, a plant growing model based on competition for space. Compared to video game RL environments, this simulator addresses a real-world problem and serves as a test bed to visualize plant growth and movement in a faster way than physical experiments. GrowSpace is composed of a suite of challenges that tackle several problems such as control, multi-stage learning,fairness and multi-objective learning. We provide agent baselines alongside case studies to demonstrate the difficulty of the proposed benchmark.


title: Adaptive and Fair Transformation for Recoverable Mutual Exclusion

abstract: Mutual exclusion is one of the most commonly used techniques to handle contention in concurrent systems. Traditionally, mutual exclusion algorithms have been designed under the assumption that a process does not fail while acquiring/releasing a lock or while executing its critical section. However, failures do occur in real life, potentially leaving the lock in an inconsistent state. This gives rise to the problem of recoverable mutual exclusion (RME) that involves designing a mutual exclusion (ME) algorithm that can tolerate failures, while maintaining safety and liveness properties. In this work, we present a framework that transforms any algorithm that solves the RME problem into an algorithm that can also simultaneously adapt to (1) the number of processes competing for the lock, as well as (2) the number of failures that have occurred in the recent past, while maintaining the correctness and performance properties of the underlying RME algorithm. Additionally, the algorithm constructed as a result of this transformation adds certain desirable properties like fairness (a variation of FCFS) and bounded recovery. Assume that the worst-case RMR complexity of a critical section request in the underlying RME algorithm is $R(n)$. Then, our framework yields an RME algorithm for which the worst-case RMR complexity of a critical section request is given by $\mathcal{O}(\min \{\ddot{c}, \sqrt{F+1}, R(n)\})$, where $\ddot{c}$ denotes the point contention of the request and $F$ denotes the number of failures in the recent past of the request. We further extend our framework by presenting a novel memory reclamation algorithm to bound the worst-case space complexity of the RME algorithm. The memory reclamation techniques maintain the fairness, performance and correctness properties of our transformation and is general enough to be employed to bound the space of other RME algorithms.


title: On endomorphisms of automatic groups

abstract: We propose two geometric versions of the bounded reduction property and find conditions for them to coincide. In particular, for the natural automatic structure on a hyperbolic group, the two notions are equivalent. We study endomorphisms with $L$-quasiconvex image and prove that those with finite kernel satisfy a synchronous version of the bounded reduction property. Finally, we use these techniques to prove $L$-quasiconvexity of the equalizer of two endomorphisms under certain (strict) conditions.


title: Bias of Root Numbers for Hilbert Newforms of Cubic Level

abstract: We give a general formula of the bias of root numbers for Hilbert modular newforms of cubic level. Explicit calculation is given when the base field is $\mathbb{Q}, \mathbb{Q}(\sqrt{2}), \mathbb{Q}(\sqrt{5})$ and the level is the cube of certain rational integers. This complements a previous result of the second author and extends the bias phenomenon to the number fields. Our method is based on Jacquet-Zagier's trace formula, and the explicit calculation works generally for all real quadratic fields of narrow class number one and for rational cubic levels.


title: Non-existing and ill-behaved coequalizers of locally ordered spaces

abstract: Categories of locally ordered spaces are especially well-adapted to the realization of most precubical sets, though their colimits are not so easy to determine (in comparison with colimits in the category of d-spaces for example). We use the plural here, as the notion of a locally ordered space vary from an author to another, only differing according to seemingly anodyne technical details. As we explain in this article, these differences have dramatic consequences on colimits. In particular, we show that most categories of locally ordered spaces are not cocomplete, thus answering a question that was neglected so far. The strategy is the following: given a directed loop {\gamma} on a locally ordered space X, we try to identify the image of {\gamma} with a single point. If it were taken in the category of d-spaces, such an identification would be likely to create a vortex, while locally ordered spaces have no vortices. Concretely, the antisymmetry of local orders gets more points to be identified than in a mere topological quotient. However, the effect of this phenomenon is in some sense limited to the neighbourhood of (the image of) {\gamma}. So the existence and the nature of the corresponding coequalizer strongly depends on the topology around the image of {\gamma}. As an extreme example, if the latter forms a connected component, the coequalizer exists and its underlying space matches with the topological coequalizer.


title: Hanner's Inequality For Positive Semidefinite Matrices

abstract: We prove an analogous Hanner's Inequality of $L^p$ spaces for positive semidefinite matrices. Let $||X||_p=\text{Tr}[(X^\ast X)^{p/2}]^{1/p}$ denote the $p$-Schatten norm of a matrix $X\in M_{n\times n}(\mathbb{C})$. We show that the inequality $||X+Y||_p^p+||X-Y||_p^p\geq (||X||_p+||Y||_p)^p+(|||X||_p-||Y||_p|)^p$ holds for $1\leq p\leq 2$ and reverses for $p\geq 2$ when $X,Y\in M_{n\times n}(\mathbb{C})^+$. This was previously known in the $1<p\leq 4/3$, $p=2$, and $p\geq 4$ cases, or with additional special assumptions. We outline these previous methods, and comment on their failure to extend to the general case. We further show that there is equality if and only if $Y=cX$, which is analogous to the equality case in $L^p$. With the general inequality, it is confirmed that the unit ball in $C^{p}_+$ has the same moduli of smoothness and convexity as the unit ball in $L^p$.


title: Reduced Order Dynamical Models For Complex Dynamics in Manufacturing and Natural Systems Using Machine Learning

abstract: Dynamical analysis of manufacturing and natural systems provides critical information about production of manufactured and natural resources respectively, thus playing an important role in assessing sustainability of these systems. However, current dynamic models for these systems exist as mechanistic models, simulation of which is computationally intensive and does not provide a simplified understanding of the mechanisms driving the overall dynamics. For such systems, lower-order models can prove useful to enable sustainability analysis through coupled dynamical analysis. There have been few attempts at finding low-order models of manufacturing and natural systems, with existing work focused on model development of individual mechanism level. This work seeks to fill this current gap in the literature of developing simplified dynamical models for these systems by developing reduced-order models using a machine learning (ML) approach. The approach is demonstrated on an entire soybean-oil to soybean-diesel process plant and a lake system. We use a grey-box ML method with a standard nonlinear optimization approach to identify relevant models of governing dynamics as ODEs using the data simulated from mechanistic models. Results show that the method identifies a high accuracy linear ODE models for the process plant, reflective of underlying linear stoichiometric mechanisms and mass balance driving the dynamics. For the natural systems, we modify the ML approach to include the effect of past dynamics, which gives non-linear ODE. While the modified approach provides a better match to dynamics of stream flow, it falls short of completely recreating the dynamics. We conclude that the proposed ML approach work well for systems where dynamics is smooth, such as in manufacturing plant whereas does not work perfectly well in case of chaotic dynamics such as water stream flow.


title: 3D Human Pose Estimation for Free-form Activity Using WiFi Signals

abstract: WiFi human sensing has become increasingly attractive in enabling emerging human-computer interaction applications. The corresponding technique has gradually evolved from the classification of multiple activity types to more fine-grained tracking of 3D human poses. However, existing WiFi-based 3D human pose tracking is limited to a set of predefined activities. In this work, we present Winect, a 3D human pose tracking system for free-form activity using commodity WiFi devices. Our system tracks free-form activity by estimating a 3D skeleton pose that consists of a set of joints of the human body. In particular, we combine signal separation and joint movement modeling to achieve free-form activity tracking. Our system first identifies the moving limbs by leveraging the two-dimensional angle of arrival of the signals reflected off the human body and separates the entangled signals for each limb. Then, it tracks each limb and constructs a 3D skeleton of the body by modeling the inherent relationship between the movements of the limb and the corresponding joints. Our evaluation results show that Winect is environment-independent and achieves centimeter-level accuracy for free-form activity tracking under various challenging environments including the none-line-of-sight (NLoS) scenarios.


title: Categorical blow-up formula for Hilbert schemes of points

abstract: Let $S$ be a smooth projective surface, and $\hat{S}$ be its blow-up at a point. In this paper, we study the derived category of the Hilbert scheme of points on the blow-up $\hat{S}$. We obtain a semi-orthogonal decomposition consisting of the derived categories of the Hilbert schemes on the original surface $S$, which recovers the blow-up formula for the Euler characteristics obtained by G\"ottsche and Nakajima-Yoshioka. The proof uses the Quot formula, which was conjectured by Jiang and recently proved by Toda.


title: Stability of an Ionization Front in Bondi Accretion

abstract: Spherical Bondi accretion is used in astrophysics as an approximation to investigate many types of accretion processes. Two-phase accretion flows that transition from neutral to ionized have observational support in high-mass star formation, and have application to accretion flows around any ionizing source, but the hydrodynamic stability of two-phase Bondi accretion is not understood. With both semi-analytic and fully numerical methods we find that these flows may be stable, conditionally stable or unstable depending on the initial conditions. The transition from an R-type to a D-type ionization front plays a key role in conditionally stable and unstable flows.


title: Coverage Probability of Double-IRS Assisted Communication Systems

abstract: In this paper, we focus on the coverage probability of a double-intelligent reflecting surface (IRS) assisted wireless network and study the impact of multiplicative beamforming gain and correlated Rayleigh fading. In particular, we obtain a novel closed-form expression of the coverage probability of a single-input single-output (SISO) system assisted by two large IRSs while being dependent on the corresponding arbitrary reflecting beamforming matrices (RBMs) and large-scale statistics in terms of correlation matrices. Taking advantage of the large-scale statistics, i.e., statistical channel state information (CSI), we perform optimization of the RBMs of both IRSs once per several coherence intervals rather than at each interval. Hence, we achieve a reduction of the computational complexity, otherwise increased in multi-IRS-assisted networks during their RBM optimization. Numerical results validate the analytical expressions even for small IRSs, confirm enhanced performance over the conventional single-IRS counterpart, and reveal insightful properties.


title: Dynamic probabilistic logic models for effective abstractions in RL

abstract: State abstraction enables sample-efficient learning and better task transfer in complex reinforcement learning environments. Recently, we proposed RePReL (Kokel et al. 2021), a hierarchical framework that leverages a relational planner to provide useful state abstractions for learning. We present a brief overview of this framework and the use of a dynamic probabilistic logic model to design these state abstractions. Our experiments show that RePReL not only achieves better performance and efficient learning on the task at hand but also demonstrates better generalization to unseen tasks.


title: A Proof of the Optimal Leapfrogging Conjecture

abstract: Suppose we place checkers in the lower left corner of a Go board and wish to move them to the upper right corner in as few moves as possible, where the pieces move as in the game of Chinese checkers. Auslander, Benjamin, and Wilkerson in 1993 generalized this game for integer lattices and defined a measure of speed for a starting configuration of pieces. They proved that the maximum speed of any configuration is 1, and only three configurations, called "speed-of-light" configurations, attain this speed. We prove their conjecture that the maximum speed of a non-speed-of-light configuration is 2/3 in the 2-dimensional case, and present a framework that should extend to higher dimensions.


title: Semimartingale and continuous-time Markov chain approximation for rough stochastic local volatility models

abstract: Rough volatility models have recently been empirically shown to provide a good fit to historical volatility time series and implied volatility smiles of SPX options. They are continuous-time stochastic volatility models, whose volatility process is driven by a fractional Brownian motion with Hurst parameter less than half. Due to the challenge that it is neither a semimartingale nor a Markov process, there is no unified method that not only applies to all rough volatility models, but also is computationally efficient. This paper proposes a semimartingale and continuous-time Markov chain (CTMC) approximation approach for the general class of rough stochastic local volatility (RSLV) models. In particular, we introduce the perturbed stochastic local volatility (PSLV) model as the semimartingale approximation for the RSLV model and establish its existence , uniqueness and Markovian representation. We propose a fast CTMC algorithm and prove its weak convergence. Numerical experiments demonstrate the accuracy and high efficiency of the method in pricing European, barrier and American options. Comparing with existing literature, a significant reduction in the CPU time to arrive at the same level of accuracy is observed.


title: Efficient privacy-preserving inference for convolutional neural networks

abstract: The processing of sensitive user data using deep learning models is an area that has gained recent traction. Existing work has leveraged homomorphic encryption (HE) schemes to enable computation on encrypted data. An early work was CryptoNets, which takes 250 seconds for one MNIST inference. The main limitation of such approaches is that of the expensive FFT-like operations required to perform operations on HE-encrypted ciphertext. Others have proposed the use of model pruning and efficient data representations to reduce the number of HE operations required. We focus on improving upon existing work by proposing changes to the representations of intermediate tensors during CNN inference. We construct and evaluate private CNNs on the MNIST and CIFAR-10 datasets, and achieve over a two-fold reduction in the number of operations used for inferences of the CryptoNets architecture.


title: Robustness of different loss functions and their impact on networks learning capability

abstract: Recent developments in AI have made it ubiquitous, every industry is trying to adopt some form of intelligent processing of their data. Despite so many advances in the field, AIs full capability is yet to be exploited by the industry. Industries that involve some risk factors still remain cautious about the usage of AI due to the lack of trust in such autonomous systems. Present-day AI might be very good in a lot of things but it is very bad in reasoning and this behavior of AI can lead to catastrophic results. Autonomous cars crashing into a person or a drone getting stuck in a tree are a few examples where AI decisions lead to catastrophic results. To develop insight and generate an explanation about the learning capability of AI, we will try to analyze the working of loss functions. For our case, we will use two sets of loss functions, generalized loss functions like Binary cross-entropy or BCE and specialized loss functions like Dice loss or focal loss. Through a series of experiments, we will establish whether combining different loss functions is better than using a single loss function and if yes, then what is the reason behind it. In order to establish the difference between generalized loss and specialized losses, we will train several models using the above-mentioned losses and then compare their robustness on adversarial examples. In particular, we will look at how fast the accuracy of different models decreases when we change the pixels corresponding to the most salient gradients.


title: On Learning the Transformer Kernel

abstract: In this work we introduce KERNELIZED TRANSFORMER, a generic, scalable, data driven framework for learning the kernel function in Transformers. Our framework approximates the Transformer kernel as a dot product between spectral feature maps and learns the kernel by learning the spectral distribution. This not only helps in learning a generic kernel end-to-end, but also reduces the time and space complexity of Transformers from quadratic to linear. We show that KERNELIZED TRANSFORMERS achieve performance comparable to existing efficient Transformer architectures, both in terms of accuracy as well as computational efficiency. Our study also demonstrates that the choice of the kernel has a substantial impact on performance, and kernel learning variants are competitive alternatives to fixed kernel Transformers, both in long as well as short sequence tasks.


title: Mitigating Membership Inference Attacks by Self-Distillation Through a Novel Ensemble Architecture

abstract: Membership inference attacks are a key measure to evaluate privacy leakage in machine learning (ML) models. These attacks aim to distinguish training members from non-members by exploiting differential behavior of the models on member and non-member inputs. The goal of this work is to train ML models that have high membership privacy while largely preserving their utility; we therefore aim for an empirical membership privacy guarantee as opposed to the provable privacy guarantees provided by techniques like differential privacy, as such techniques are shown to deteriorate model utility. Specifically, we propose a new framework to train privacy-preserving models that induces similar behavior on member and non-member inputs to mitigate membership inference attacks. Our framework, called SELENA, has two major components. The first component and the core of our defense is a novel ensemble architecture for training. This architecture, which we call Split-AI, splits the training data into random subsets, and trains a model on each subset of the data. We use an adaptive inference strategy at test time: our ensemble architecture aggregates the outputs of only those models that did not contain the input sample in their training data. We prove that our Split-AI architecture defends against a large family of membership inference attacks, however, it is susceptible to new adaptive attacks. Therefore, we use a second component in our framework called Self-Distillation to protect against such stronger attacks. The Self-Distillation component (self-)distills the training dataset through our Split-AI ensemble, without using any external public datasets. Through extensive experiments on major benchmark datasets we show that SELENA presents a superior trade-off between membership privacy and utility compared to the state of the art.


title: Minor Embedding in Broken Chimera and Pegasus Graphs is NP-complete

abstract: The embedding is an essential step when calculating on the D-Wave machine. In this work we show the hardness of the embedding problem for both types of existing hardware, represented by the Chimera and the Pegasus graphs, containing unavailable qubits. We construct certain broken Chimera graphs, where it is hard to find a Hamiltonian cycle. As the Hamiltonian cycle problem is a special case of the embedding problem, this proves the general complexity result for the Chimera graphs. By exploiting the subgraph relation between the Chimera and the Pegasus graphs, the proof is then further extended to the Pegasus graphs.


title: Comparing One-step and Two-step Scatter Correction and Density Reconstruction in X-ray CT

abstract: In this work, we compare one-step and two-step approaches for X-ray computed tomography (CT) scatter correction and density reconstruction. X-ray CT is an important imaging technique in medical and industrial applications. In many cases, the presence of scattered X-rays leads to loss of contrast and undesirable artifacts in reconstructed images. Many approaches to computationally removing scatter treat scatter correction as a preprocessing step that is followed by a reconstruction step. Treating scatter correction and reconstruction jointly as a single, more complicated optimization problem is less studied. It is not clear from the existing literature how these two approaches compare in terms of reconstruction accuracy. In this paper, we compare idealized versions of these two approaches with synthetic experiments. Our results show that the one-step approach can offer improved reconstructions over the two-step approach, although the gap between them is highly object-dependent.


title: Solving Image PDEs with a Shallow Network

abstract: Partial differential equations (PDEs) are typically used as models of physical processes but are also of great interest in PDE-based image processing. However, when it comes to their use in imaging, conventional numerical methods for solving PDEs tend to require very fine grid resolution for stability, and as a result have impractically high computational cost. This work applies BLADE (Best Linear Adaptive Enhancement), a shallow learnable filtering framework, to PDE solving, and shows that the resulting approach is efficient and accurate, operating more reliably at coarse grid resolutions than classical methods. As such, the model can be flexibly used for a wide variety of problems in imaging.


title: On the stability of homogeneous black strings in AdS

abstract: In this paper we extend the analysis of the stability of an homogeneous black string in the presence of a negative cosmological constant with minimally coupled scalar fields. We recall the linear stability of this solutions under generic perturbations on the metric and of the scalar fields. Then, we extend the study of the stability by presenting the existence of a non-generic perturbation which may lead to an unstable behavior. The later mode is fine-tuned since it requires the scalar field degree of freedom to be absent through the whole evolution of the system.


title: Control Prefixes for Parameter-Efficient Text Generation

abstract: Prefix-tuning is a powerful lightweight technique for adapting a large pre-trained language model to a downstream application. However, it uses the same dataset-level tuned prompt for all examples in the dataset. We extend this idea and propose a dynamic method, Control Prefixes, which allows for the inclusion of conditional input-dependent information, combining the benefits of prompt tuning and controlled generation. The method incorporates attribute-level learnable representations into different layers of a pre-trained transformer, allowing for the generated text to be guided in a particular direction. We provide a systematic evaluation of the technique and apply it to five datasets from the GEM benchmark for natural language generation (NLG). Although the aim is to develop a parameter-efficient model, we show Control Prefixes can even outperform full fine-tuning methods. We present state-of-the-art results on several data-to-text datasets, including WebNLG.


title: Nothing Wasted: Full Contribution Enforcement in Federated Edge Learning

abstract: The explosive amount of data generated at the network edge makes mobile edge computing an essential technology to support real-time applications, calling for powerful data processing and analysis provided by machine learning (ML) techniques. In particular, federated edge learning (FEL) becomes prominent in securing the privacy of data owners by keeping the data locally used to train ML models. Existing studies on FEL either utilize in-process optimization or remove unqualified participants in advance. In this paper, we enhance the collaboration from all edge devices in FEL to guarantee that the ML model is trained using all available local data to accelerate the learning process. To that aim, we propose a collective extortion (CE) strategy under the imperfect-information multi-player FEL game, which is proved to be effective in helping the server efficiently elicit the full contribution of all devices without worrying about suffering from any economic loss. Technically, our proposed CE strategy extends the classical extortion strategy in controlling the proportionate share of expected utilities for a single opponent to the swiftly homogeneous control over a group of players, which further presents an attractive trait of being impartial for all participants. Moreover, the CE strategy enriches the game theory hierarchy, facilitating a wider application scope of the extortion strategy. Both theoretical analysis and experimental evaluations validate the effectiveness and fairness of our proposed scheme.


title: A New Approach for Interpretability and Reliability in Clinical Risk Prediction: Acute Coronary Syndrome Scenario

abstract: We intend to create a new risk assessment methodology that combines the best characteristics of both risk score and machine learning models. More specifically, we aim to develop a method that, besides having a good performance, offers a personalized model and outcome for each patient, presents high interpretability, and incorporates an estimation of the prediction reliability which is not usually available. By combining these features in the same approach we expect that it can boost the confidence of physicians to use such a tool in their daily activity. In order to achieve the mentioned goals, a three-step methodology was developed: several rules were created by dichotomizing risk factors; such rules were trained with a machine learning classifier to predict the acceptance degree of each rule (the probability that the rule is correct) for each patient; that information was combined and used to compute the risk of mortality and the reliability of such prediction. The methodology was applied to a dataset of patients admitted with any type of acute coronary syndromes (ACS), to assess the 30-days all-cause mortality risk. The performance was compared with state-of-the-art approaches: logistic regression (LR), artificial neural network (ANN), and clinical risk score model (Global Registry of Acute Coronary Events - GRACE). The proposed approach achieved testing results identical to the standard LR, but offers superior interpretability and personalization; it also significantly outperforms the GRACE risk model and the standard ANN model. The calibration curve also suggests a very good generalization ability of the obtained model as it approaches the ideal curve. Finally, the reliability estimation of individual predictions presented a great correlation with the misclassifications rate. Those properties may have a beneficial application in other clinical scenarios as well. [abridged]


title: Topological Quantum Critical Points in Strong Coupling limits: Global Symmetries and Strongly Interacting Majorana Fermions

abstract: In this article, we discuss strong coupling limits of topological quantum critical points (TQCPs) where quantum phase transitions between two topological distinct superconducting states take place. We illustrate that while superconducting phases on both sides of TQCPs spontaneously break same symmetries, universality classes of critical states can be identified only when global symmetries in topological states are further specified. In dimensions $d=2,3$, we find that continuous $(d+1)$th order transitions at weakly interacting TQCPs that were pointed out previously in the presence of emergent Lorentz symmetry can be terminated by strongly interacting fixed points of majorana fields. For $2d$ time reversal symmetry breaking TQCPs, termination points are supersymmetric with ${\mathcal N}=4N_f={1}$ (where $N_f$ is the number of four-component Dirac fermions and ${\mathcal N}$ is the number of two-component real fermions) beyond which transitions are discontinuous first order ones. For $2d$ time reversal symmetric TQCPs without other global symmetries, termination points of $(d+1)$th order continuous transition lines are generically conformal invariant without supersymmetry. Beyond these strong coupling fixed points, there are first-order discontinuous transitions as far as the protecting symmetry is not spontaneously broken but no direct transitions if the protecting symmetry is spontaneously broken in the presence of strong interactions. In $3d$, strong coupling termination points can be further effectively represented by new emergent gapless real bosons weakly coupled with free gapless majorana fermions. However, in $1d$, time reversal symmetric $(d+1)$th continuous transition lines of TQCPs are terminated by simple free majorana fermion fixed points.


title: Probing the Earth's Core using Atmospheric Neutrinos at INO

abstract: The proposed 50 kt Iron Calorimeter (ICAL) detector at the India-based Neutrino Observatory (INO) aims to detect atmospheric muon neutrinos and antineutrinos separately in the multi-GeV range of energies and over a wide range of path lengths. While passing through the Earth, the upward-going neutrinos experience a density-dependent matter effect, which can be utilized to extract information about the internal structure of Earth. Since the Earth's matter effect modifies the neutrino oscillation patterns differently for neutrinos and antineutrinos, the capability of ICAL to distinguish $\mu^-$ and $\mu^+$ events plays an important role in observing this matter effect. Taking advantage of good angular resolution, ICAL would be able to observe about 331 $\mu^-$ and 146 $\mu^+$ events corresponding to the core-passing neutrinos and antineutrinos, respectively, in 10 years. We demonstrate for the first time that ICAL would be able to validate the presence of Earth's core by ruling out a two-layered profile consisting of only mantle and crust in fit with respect to the PREM profile in data with a median $\Delta \chi^2$ of 7.45 for normal mass ordering (NO) and 4.83 for inverted mass ordering (IO) using 500 kt$\cdot$yr exposure. If we do not use the charge identification capability of ICAL, these sensitivities deteriorate to a $\Delta\chi^2$ of 3.76 for NO and 1.59 for IO.


title: Gravitational Instabilities of Uniform Black Strings in AdS

abstract: Locally AdS$_{d-1}\times\mathbb{R}$ uniform black strings (UBS) in the presence of a massless scalar field are believed to avoid the onset of the Gregory-Laflamme (GL) instability in $d\ge4$ as no tachyonic modes exist in the spectrum of the Laplace-Beltrami operator. We present analytic and numerical evidence of GL modes in the Lichnerowicz spectrum indicating that AdS$_{d-1}$ UBSs are classically and thermodynamically unstable at the linear level in $d>4$. In $d=4$, we confirm that uniform BTZ$_3$ strings are indeed stable as previously suggested. We propose that linear instabilities of black strings are triggered if and only if a tachyonic mode exists in the Lichnerowicz spectrum. At the end state of the instability, AdS$_{d-1}$ UBSs of finite length may tunnel to a SAdS$_d$ black hole or converge onto a novel non-uniform AdS$_d$ black string. We conjecture that weak cosmic censorship is violated if the non-uniform solution is an exact AdS$_d$ black funnel and compute entropy estimates in $d>4$ as evidence.


title: Trigger Hunting with a Topological Prior for Trojan Detection

abstract: Despite their success and popularity, deep neural networks (DNNs) are vulnerable when facing backdoor attacks. This impedes their wider adoption, especially in mission critical applications. This paper tackles the problem of Trojan detection, namely, identifying Trojaned models -- models trained with poisoned data. One popular approach is reverse engineering, i.e., recovering the triggers on a clean image by manipulating the model's prediction. One major challenge of reverse engineering approach is the enormous search space of triggers. To this end, we propose innovative priors such as diversity and topological simplicity to not only increase the chances of finding the appropriate triggers but also improve the quality of the found triggers. Moreover, by encouraging a diverse set of trigger candidates, our method can perform effectively in cases with unknown target labels. We demonstrate that these priors can significantly improve the quality of the recovered triggers, resulting in substantially improved Trojan detection accuracy as validated on both synthetic and publicly available TrojAI benchmarks.


title: Rapid-cycle Thouless pumping in a one-dimensional optical lattice

abstract: An adiabatic cycle around a degeneracy point in the parameter space of a one-dimensional band insulator is known to result in an integer valued noiseless particle transport in the thermodynamic limit. Recently, it was shown that in the case of an infinite bipartite lattice the adiabatic Thouless protocol can be continuously deformed into a fine tuned finite-frequency cycle preserving the properties of noiseless quantized transport. In this paper, we numerically investigate the implementation of such an ideal rapid-cycle Thouless pumping protocol in a one-dimensional optical lattice. It is shown that the rapidity will cause first order corrections due to next-to-nearest-neighbour hopping and second order corrections due to the addition of a harmonic potential. Lastly, the quantization of the change in center of mass of the particle distribution is investigated, and shown to have corrections in the first order of the potential curvature.


title: On the Integrability of Pfaffian Forms on ${\mathbb R}^{n}$

abstract: This paper details the lesser known conditions on ${\mathbb {R}}^{n}$ for the integrability of pfaffian forms, or 1-forms. Emphasis is given to locality of these conditions, and proofs in some additional detail are provided for theorems due to Clairaut and Carath\'eodory. Considering the importance of the integrability of pfaffian forms, in particular in mathematical-physics, this paper shows that: there is a hidden content in Carath\'eodory's theorem in the direction of a global integrability.


title: Exploratory Lagrangian-Based Particle Tracing Using Deep Learning

abstract: Time-varying vector fields produced by computational fluid dynamics simulations are often prohibitively large and pose challenges for accurate interactive analysis and exploration. To address these challenges, reduced Lagrangian representations have been increasingly researched as a means to improve scientific time-varying vector field exploration capabilities. This paper presents a novel deep neural network-based particle tracing method to explore time-varying vector fields represented by Lagrangian flow maps. In our workflow, in situ processing is first utilized to extract Lagrangian flow maps, and deep neural networks then use the extracted data to learn flow field behavior. Using a trained model to predict new particle trajectories offers a fixed small memory footprint and fast inference. To demonstrate and evaluate the proposed method, we perform an in-depth study of performance using a well-known analytical data set, the Double Gyre. Our study considers two flow map extraction strategies as well as the impact of the number of training samples and integration durations on efficacy, evaluates multiple sampling options for training and testing and informs hyperparameter settings. Overall, we find our method requires a fixed memory footprint of 10.5 MB to encode a Lagrangian representation of a time-varying vector field while maintaining accuracy. For post hoc analysis, loading the trained model costs only two seconds, significantly reducing the burden of I/O when reading data for visualization. Moreover, our parallel implementation can infer one hundred locations for each of two thousand new pathlines across the entire temporal resolution in 1.3 seconds using one NVIDIA Titan RTX GPU.


title: A Static Analysis Framework for Data Science Notebooks

abstract: Notebooks provide an interactive environment for programmers to develop code, analyse data and inject interleaved visualizations in a single environment. Despite their flexibility, a major pitfall that data scientists encounter is unexpected behaviour caused by the unique out-of-order execution model of notebooks. As a result, data scientists face various challenges ranging from notebook correctness, reproducibility and cleaning. In this paper, we propose a framework that performs static analysis on notebooks, incorporating their unique execution semantics. Our framework is general in the sense that it accommodate for a wide range of analyses, useful for various notebook use cases. We have instantiated our framework on a diverse set of analyses and have evaluated them on 2211 real world notebooks. Our evaluation demonstrates that the vast majority (98.7%) of notebooks can be analysed in less than a second, well within the time frame required by interactive notebook clients


title: Return migration of German-affiliated researchers: Analyzing departure and return by gender, cohort, and discipline using Scopus bibliometric data 1996-2020

abstract: The international migration of researchers is an important dimension of scientific mobility, and has been the subject of considerable policy debate. However, tracking the migration life courses of researchers is challenging due to data limitations. In this study, we use Scopus bibliometric data on eight million publications from 1.1 million researchers who have published at least once with an affiliation address from Germany in 1996-2020. We construct the partial life histories of published researchers in this period and explore both their out-migration and the subsequent return of a subset of this group: the returnees. Our analyses shed light on the career stages and gender disparities between researchers who remain in Germany, those who emigrate, and those who eventually return. We find that the return migration streams are even more gender imbalanced, which points to the need for additional efforts to encourage female researchers to come back to Germany. We document a slightly declining trend in return migration among more recent cohorts of researchers who left Germany, which, for most disciplines, was associated with a decrease in the German collaborative ties of these researchers. Moreover, we find that the gender disparities for the most gender imbalanced disciplines are unlikely to be mitigated by return migration given the gender compositions of the cohorts of researchers who have left Germany and of those who have returned. This analysis uncovers new dimensions of migration among scholars by investigating the return migration of published researchers, which is critical for the development of science policy.


title: Regional Stability Analysis of Transitional Fluid Flows

abstract: A method to bound the maximum energy perturbation for which regional stability of transitional fluid flow models can be guaranteed is introduced. The proposed method exploits the fact that the fluid model's nonlinearities are both lossless and locally bounded and uses the axes lengths of the ellipsoids for the trajectory set containment as variables in the stability conditions. Compared to existing approaches, the proposed method leads to an average increase in the maximum allowable energy perturbation of 29% for the Waleffe-KimHamilton (WKH) shear flow model and of 38% for the 9-state reduced model of Couette flow.


title: Simple Periodic Boundary Conditions for Molecular Simulation of Uniaxial Flow

abstract: We present rotating periodic boundary conditions (PBCs) for the simulation of nonequilibrium molecular dynamics (NEMD) under uniaxial stretching flow (USF) or biaxial stretching flow (BSF). Such nonequilibrium flows need specialized PBCs since the simulation box deforms with the background flow. The technique builds on previous models using one or lattice remappings, and is simpler than the PBCs developed for the general three dimensional flow. For general three dimensional flows, Dobson \cite{Dobson} and Hunt \cite{Hunt} proposed schemes which are not time-periodic since they use more than one automorphism remapping. This paper presents a single automorphism remapping PBCs for USF and BSF which is time periodic up to a rotation matrix and has a better minimum lattice spacing properties.


title: Hyperseed: Unsupervised Learning with Vector Symbolic Architectures

abstract: Motivated by recent innovations in biologically-inspired neuromorphic hardware, this article presents a novel unsupervised machine learning algorithm named Hyperseed that draws on the principles of Vector Symbolic Architectures (VSA) for fast learning of a topology preserving feature map of unlabelled data. It relies on two major operations of VSA, binding and bundling. The algorithmic part of Hyperseed is expressed within Fourier Holographic Reduced Representations model, which is specifically suited for implementation on spiking neuromorphic hardware. The two primary contributions of the Hyperseed algorithm are, few-shot learning and a learning rule based on single vector operation. These properties are empirically evaluated on synthetic datasets as well as on illustrative benchmark use-cases, IRIS classification, and a language identification task using n-gram statistics. The results of these experiments confirm the capabilities of Hyperseed and its applications in neuromorphic hardware.


title: Growth and Dissipation of Be Star Discs in Misaligned Binary Systems

abstract: We use a three-dimensional smoothed particle hydrodynamics code to simulate growth and dissipation of Be star discs in systems where the binary orbit is misaligned with respect to the spin axis of the primary star. We investigate six different scenarios of varying orbital period and misalignment angle, feeding the disc at a constant rate for 100 orbital periods, and then letting the disc dissipate for 100 orbital periods. During the disc growth phase, we find that the binary companion tilts the disc away from its initial plane at the equator of the primary star before settling to a constant orientation after 40 to 50 orbital periods. While the mass-injection into the disc is ongoing, the tilting of the disc can cause material to reaccrete onto the primary star prematurely. Once disc dissipation begins, usually the disc precesses about the binary companion's orbital axis with precession periods ranging from 20 to 50 orbital periods. In special cases we detect phenomena of disc tearing, as well as Kozai-Lidov oscillations of the disc. These oscillations reach a maximum eccentricity of about 0.6, and a minimum inclination of about \ang{20} with respect to the binary's orbit. We also find the disc material to have highly eccentric orbits beyond the transition radius, where the disc changes from being dominated by viscous forces, to heavily controlled by the companion star, in contrast to its nearly circular motion inward of the transition radius. Finally, we offer predictions to how these changes will affect Be star observables.


title: Towards Transparent Interactive Semantic Parsing via Step-by-Step Correction

abstract: Existing studies on semantic parsing focus primarily on mapping a natural-language utterance to a corresponding logical form in one turn. However, because natural language can contain a great deal of ambiguity and variability, this is a difficult challenge. In this work, we investigate an interactive semantic parsing framework that explains the predicted logical form step by step in natural language and enables the user to make corrections through natural-language feedback for individual steps. We focus on question answering over knowledge bases (KBQA) as an instantiation of our framework, aiming to increase the transparency of the parsing process and help the user appropriately trust the final answer. To do so, we construct INSPIRED, a crowdsourced dialogue dataset derived from the ComplexWebQuestions dataset. Our experiments show that the interactive framework with human feedback has the potential to greatly improve overall parse accuracy. Furthermore, we develop a pipeline for dialogue simulation to evaluate our framework w.r.t. a variety of state-of-the-art KBQA models without involving further crowdsourcing effort. The results demonstrate that our interactive semantic parsing framework promises to be effective across such models.


title: Implementation of a Multiple Target Tracking Filter on an Adiabatic Quantum Computer

abstract: Recent work at Fraunhofer FKIE shows that Morefield's method for multiple target data association can in theory be solved on an adiabatic quantum computer. The present paper validates the theory and examines the significant limitations of currently available adiabatic quantum computers for solving the data association problem. The limitations of such architectures are both theoretical and practical in nature, and both are discussed. The data association problem is formulated as a quadratic unconstrained binary optimization (QUBO) problem; consequently, much of the discussion is relevant to other applications which are, or can be, posed as QUBO problems.


title: Shapes, Softness and Non-Yrast Collectivity in 186W

abstract: Non-yrast, excited states in neutron-rich $^{186}$W were populated via inelastic-scattering reactions using beams of $^{136}$Xe nuclei accelerated to 725 and 800 MeV. Levels populated in the reactions were investigated via particle-$ \gamma $ coincidence techniques using the Gammasphere array of High-Purity Germanium detectors and the compact heavy-ion counter, CHICO2. The $ K^{\pi} = 2 ^{+} $ ($\gamma$), $ K^{\pi} = 0^{+}$ and $ K^{\pi} = 2^{-} $ (octupole) rotational side bands were extended to spins $ 14\hbar $, $ 12\hbar $, and $ 13\hbar $, respectively. A staggering pattern observed in the energies of levels in the $ K^{\pi} = 2^{+} $ band was found to be consistent with a potential that gets softer to vibration in the $ \gamma $ degree of freedom with increasing spin. The odd-even staggering of states in the $ K^{\pi} = 2^{-}$ band was found to exhibit a phase opposite to that seen in the $ \gamma $ band; an effect most probably associated with Coriolis coupling to other, unobserved octupole vibrational bands in $^{186}$W.


title: A central limit theorem concerning uncertainty in estimates of individual admixture

abstract: The concept of individual admixture (IA) assumes that the genome of individuals is composed of alleles inherited from $K$ ancestral populations. Each copy of each allele has the same chance $q_k$ to originate from population $k$, and together with the allele frequencies $p$ in all populations at all $M$ markers, comprises the admixture model. Here, we assume a supervised scheme, i.e.\ allele frequencies $p$ are given through a reference database of size $N$, and $q$ is estimated via maximum likelihood for a single sample. We study laws of large numbers and central limit theorems describing effects of finiteness of both, $M$ and $N$, on the estimate of $q$. We recall results for the effect of finite $M$, and provide a central limit theorem for the effect of finite $N$, introduce a new way to express the uncertainty in estimates in standard barplots, give simulation results, and discuss applications in forensic genetics.


title: QNMs of scalar fields on small Reissner-Nordstr\"{o}m-AdS$\mathbf{_5}$ black holes

abstract: We study the quasinormal modes (QNMs) of a charged scalar field on a Reissner-Nordstr\"{o}m-anti-de Sitter (RN-AdS$_{5}$) black hole in the small radius limit by using the isomonodromic method. We also derive the low-temperature expansion of the fundamental QNM frequency. Finally, we provide numerical evidence that instabilities appear in the small radius limit for large values of the charge of the scalar field.


title: Differentiable Network Pruning for Microcontrollers

abstract: Embedded and personal IoT devices are powered by microcontroller units (MCUs), whose extreme resource scarcity is a major obstacle for applications relying on on-device deep learning inference. Orders of magnitude less storage, memory and computational capacity, compared to what is typically required to execute neural networks, impose strict structural constraints on the network architecture and call for specialist model compression methodology. In this work, we present a differentiable structured network pruning method for convolutional neural networks, which integrates a model's MCU-specific resource usage and parameter importance feedback to obtain highly compressed yet accurate classification models. Our methodology (a) improves key resource usage of models up to 80x; (b) prunes iteratively while a model is trained, resulting in little to no overhead or even improved training time; (c) produces compressed models with matching or improved resource usage up to 1.4x in less time compared to prior MCU-specific methods. Compressed models are available for download.


title: Robustness of Stark many-body localization in the $J_1$-$J_2$ Heisenberg model

abstract: Stark many-body localization (SMBL) is a phenomenon observed in interacting systems with a nearly uniform spatial gradient applied field. Contrasting to the traditional many-body localization phenomenon, SMBL does not require disorder. Here we investigate SMBL in a spin-$1/2$ described by a Heisenberg model including a next-nearest-neighbor exchange coupling. By employing an exact diagonalization approach and time evolution calculation we analyze both level spacing ratio (LSR) statistics of the Hamiltonian model as well as the dynamics of the system from a given initial state. Our results reveals that for zero field in our finite system, LSR statistics suggest localization while the dynamics shows thermalization, which has been attributed to a finite-size effect. Slightly nonuniform field gradient, LSR statistic predictions agree very well with the dynamics of the physical quantities indicating delocalization and localization for small and large field gradient, respectively. More interestingly, we find that localization is robust in the presence of next-nearest-neighbor coupling in the Hamiltonian. Moreover, this coupling can be tuned to enhance SMBL in the system, meaning that localized regimes can be obtained for smaller field gradient as compared to the traditional nearest-neighbor isotropic Heisenberg model.

categories: cond-mat.dis-nn cond-mat.soft cond-mat.stat-mech cond-mat.str-el


title: Omni-sparsity DNN: Fast Sparsity Optimization for On-Device Streaming E2E ASR via Supernet

abstract: From wearables to powerful smart devices, modern automatic speech recognition (ASR) models run on a variety of edge devices with different computational budgets. To navigate the Pareto front of model accuracy vs model size, researchers are trapped in a dilemma of optimizing model accuracy by training and fine-tuning models for each individual edge device while keeping the training GPU-hours tractable. In this paper, we propose Omni-sparsity DNN, where a single neural network can be pruned to generate optimized model for a large range of model sizes. We develop training strategies for Omni-sparsity DNN that allows it to find models along the Pareto front of word-error-rate (WER) vs model size while keeping the training GPU-hours to no more than that of training one singular model. We demonstrate the Omni-sparsity DNN with streaming E2E ASR models. Our results show great saving on training time and resources with similar or better accuracy on LibriSpeech compared to individually pruned sparse models: 2%-6.6% better WER on Test-other.


title: Revisiting Popularity and Demographic Biases in Recommender Evaluation and Effectiveness

abstract: Recommendation algorithms are susceptible to popularity bias: a tendency to recommend popular items even when they fail to meet user needs. A related issue is that the recommendation quality can vary by demographic groups. Marginalized groups or groups that are under-represented in the training data may receive less relevant recommendations from these algorithms compared to others. In a recent study, Ekstrand et al. investigate how recommender performance varies according to popularity and demographics, and find statistically significant differences in recommendation utility between binary genders on two datasets, and significant effects based on age on one dataset. Here we reproduce those results and extend them with additional analyses. We find statistically significant differences in recommender performance by both age and gender. We observe that recommendation utility steadily degrades for older users, and is lower for women than men. We also find that the utility is higher for users from countries with more representation in the dataset. In addition, we find that total usage and the popularity of consumed content are strong predictors of recommender performance and also vary significantly across demographic groups.


title: Faster Modular Composition

abstract: A new Las Vegas algorithm is presented for the composition of two polynomials modulo a third one, over an arbitrary field. When the degrees of these polynomials are bounded by $n$, the algorithm uses $O(n^{1.43})$ field operations, breaking through the $3/2$ barrier in the exponent for the first time. The previous fastest algebraic algorithms, due to Brent and Kung in 1978, require $O(n^{1.63})$ field operations in general, and ${n^{3/2+o(1)}}$ field operations in the special case of power series over a field of large enough characteristic. If cubic-time matrix multiplication is used, the new algorithm runs in ${n^{5/3+o(1)}}$ operations, while previous ones run in $O(n^2)$ operations. Our approach relies on the computation of a matrix of algebraic relations that is typically of small size. Randomization is used to reduce arbitrary input to this favorable situation.


title: Clean or Annotate: How to Spend a Limited Data Collection Budget

abstract: Crowdsourcing platforms are often used to collect datasets for training machine learning models, despite higher levels of inaccurate labeling compared to expert labeling. There are two common strategies to manage the impact of such noise. The first involves aggregating redundant annotations, but comes at the expense of labeling substantially fewer examples. Secondly, prior works have also considered using the entire annotation budget to label as many examples as possible and subsequently apply denoising algorithms to implicitly clean the dataset. We find a middle ground and propose an approach which reserves a fraction of annotations to explicitly clean up highly probable error samples to optimize the annotation process. In particular, we allocate a large portion of the labeling budget to form an initial dataset used to train a model. This model is then used to identify specific examples that appear most likely to be incorrect, which we spend the remaining budget to relabel. Experiments across three model variations and four natural language processing tasks show our approach outperforms or matches both label aggregation and advanced denoising methods designed to handle noisy labels when allocated the same finite annotation budget.


title: Operator Complexity for Quantum Scalar Fields and Cosmological Perturbations

abstract: We calculate the operator complexity for the displacement, squeeze and rotation operators of a quantum harmonic oscillator. The complexity of the time-dependent displacement operator is constant, equal to the magnitude of the coherent state parameter, while the complexity of unitary evolution by a generic quadratic Hamiltonian is proportional to the amount of squeezing and is sensitive to the time-dependent phase of the unitary operator. We apply these results to study the complexity of a free massive scalar field, finding that the complexity has a period of rapid linear growth followed by a saturation determined by the UV cutoff and the number of spatial dimensions. We also study the complexity of the unitary evolution of quantum cosmological perturbations in de Sitter space, which can be written as time-dependent squeezing and rotation operators on individual Fourier mode pairs. The complexity of a single mode pair at late times grows linearly with the number of e-folds, while the complexity at early times oscillates rapidly due to the sensitivity of operator complexity to the phase of unitary time evolution. Integrating over all modes, the total complexity of cosmological perturbations scales as the square root of the (exponentially) growing volume of de Sitter space, suggesting that inflation leads to an explosive growth in complexity of the Universe.


title: Sensitivity of the least stable modes to passive control for a flow around an elastically-mounted circular cylinder

abstract: In this paper, a methodology to calculate the sensitivity of the least stable modes of fluid-structure interaction systems with respect to local forces is presented. We make use of the adjoint equations of the flow-structure coupled system to calculate the gradients, and the algorithms were implemented using the spectral/\emph{hp} element method for the spatial discretization. The methodology was applied to two-dimensional incompressible laminar steady flows around an elastically-mounted circular cylinder, and we obtained the gradients of the real and imaginary parts of the least stable eigenvalues with respect to forces located at arbitrary points in the flow domain. Selected values of mass ratio and reduced velocity were considered in the simulations, and the results were compared to those obtained for a fixed cylinder at the same Reynolds number. We noticed that the sensitivity fields of the fluid-structure interaction system can be very different from its fixed structure counterpart, and amongst the cases, with an elastic structure, the fields vary greatly according to the reduced velocity. Finally, the sensitivity results were verified against linear and nonlinear simulations of flows with small control cylinders placed at locations selected according to the sensitivity fields. The agreement between the predictions made with the sensitivity analyses and the linear and nonlinear results of the forced flows was excellent. In some cases, it was possible to completely suppress the cylinder vibration.


title: 1321+045: a Compact Steep Spectrum radio source in a cool-core galaxy cluster

abstract: Cluster-central gigahertz peak and compact steep spectrum (CSS) sources offer an opportunity to study the earliest phases of AGN feedback, but few have yet been examined in detail. We present results from radio and X-ray observations of 1321+045, a CSS source in a 4.4 keV cluster at z=0.263. The cluster has a strongly cooling core, and disturbances from a minor cluster merger may have triggered a period of jet activity which formed the 16 kpc radio lobes 2.0 [+0.3,-0.2] Myr ago. However, new VLBA imaging shows a ~20 pc jet on a different projected axis, which is probably only a few hundred years old. We consider possible histories for the system, with either one or two periods of jet activity. While this single system is informative, a broader study of the youngest cluster-central radio sources is desirable.


title: Projected-Search Methods for Bound-Constrained Optimization

abstract: Projected-search methods for bound-constrained optimization are based on performing a search along a piecewise-linear continuous path obtained by projecting a search direction onto the feasible region. A benefit of these methods is that many changes to the active set can be made at the cost of computing a single search direction. As the objective function is not differentiable along the search path, it is not possible to use a projected-search method with a step that satisfies the Wolfe conditions, which require the directional derivative of the objective function at a point on the path. Thus, methods based on a simple backtracking procedure must be used to give a step that satisfies an "Armijo-like" sufficient decrease condition. As a consequence, conventional projected-search methods are unable to exploit sophisticated safeguarded polynomial interpolation techniques that have been shown to be effective for the unconstrained case. This paper concerns the formulation and analysis of projected-search methods based on a new quasi-Wolfe line search that is appropriate for piecewise differentiable functions. The behavior of the line search is similar to conventional Wolfe line search, except that a step is accepted under a wider range of conditions. These conditions take into consideration steps at which the restriction of the objective function on the search path is not differentiable. Two new classes of method are proposed that may be broadly categorized as active-set and interior-point methods. Computational results are given for two specific methods from these general classes: a projected-search active-set method that uses a limited-memory quasi-Newton approximation of the Hessian; and a projected-search primal-dual interior-point method. The results show that in these contexts, a quasi-Wolfe line search is substantially more efficient and reliable than an Armijo line search.


title: Type Ia Supernova Models: Asymmetric Remnants and Supernova Remnant G1.9+0.3

abstract: The youngest Galactic supernova remnant G1.9+0.3, probably the result of a Type Ia supernova, shows surprising anomalies in the distribution of its ejecta in space and velocity. In particular, high-velocity shocked iron is seen in several locations far from the remnant center, in some cases beyond prominent silicon and sulfur emission. These asymmetries strongly suggest a highly asymmetric explosion. We present high-resolution hydrodynamic simulations in two and three dimensions of the evolution from ages of 100 seconds to hundreds of years of two asymmetric Type Ia models, expanding into a uniform medium. At the age of G1.9+0.3 (about 100 years), our 2D model shows almost no iron shocked to become visible in X-rays. Only in a much higher-density environment could significant iron be shocked, at which time the model's expansion speed is completely inconsistent with the observations of G1.9+0.3. Our 3D model, evolving the most asymmetric of a suite of Type Ia SN models from Seitenzahl et al.~(2013), shows some features resembling G1.9+0.3. We characterize its evolution with images of composition in three classes: C and O, intermediate-mass elements (IMEs), and iron-group elements (IGEs). From ages of 13 to 1800 years, we follow the evolution of the highly asymmetric initial remnant as the explosion asymmetries decrease in relative strength to be replaced by asymmetries due to evolutionary hydrodynamic instabilities. At an age of about 100 years, our 3D model has comparable shocked masses of C+O, IMEs, and IGEs, with about 0.03 $M_\odot$ each. Evolutionary changes appear to be rapid enough that continued monitoring with the Chandra X-ray Observatory may show significant variations.


title: Electronic pumping of heat without charge transfer

abstract: A mechanism of electron-mediated pumping of heat in the absence of net charge transfer is proposed. It may be realized in charge-neutral electron systems, such as graphene, coupled to an external electric potential. The flow of heat in this pumping cycle is not accompanied by a buildup of voltage along the system, which offers advantages over traditional thermoelectric cooling setups. Efficiency of heat pumping and magnitude of heat flux are studied in the hydrodynamic regime for weak disorder. In a pristine system, even for an infinitesimal pumping potential the heat flux remains finite. In particular, for a potential in the form of a traveling wave moving with velocity $c$ the pumping is perfect; the entire heat content of the electron liquid is advected with velocity $c$. For a general pumping cycle the heat flux is determined by the cycle geometry and disorder strength.