title: Nonexistence of k-bounded sobrification

abstract: In this paper, we will focus on k-bounded sober spaces and show the existence of a T_0 space X not admitting any k-bounded sobrification. This strengthens a result of Zhao, Lu and Wang, who proved that the canonical k-bounded sobrification does not exist. Our work provides a complete solution to a question of Zhao and Ho, and shows that unlike Sob and BSob,the category KBSob of all k-bounded sober spaces is not a reflective subcategory of the category Top_0 of all T_0 spaces. Furthermore, we introduce the notion of qk-bounded sober spaces and prove that the category KBSob is a full reflective subcategory of the category QKBSob of all qk-bounded sober spaces and continuous mappings preserving existing irreducible suprema.


title: Cellular chains on freehedra and operadic pairs

abstract: The paper is devoted to explaining the operadic meaning of freehedra, a family of polytopes originally defined to study free loop spaces. We introduce the notion of operadic pairs and algebras over them. Cellular chains on Stasheff associahedra and Stasheff multiplihedra assemble into a two-colored operadic pair that governs $A_\infty$-algebras and $A_\infty$-modules over them, with maps that are $A_\infty$ in both colors. An important quotient of this operadic pair is given by cellular chains on cubes and freehedra.


title: Extremely Sparse Co-Prime Sensing: Low Latency Estimation is not a Dream but a Reality

abstract: Co-prime sensing is a sub-Nyquist technique for signal acquisition. Several modifications to the prototype co-prime array have been proposed in the literature. Researchers have also demonstrated low latency estimation. This paper describes the functioning of a Family of Adjustable Pivot Co-Prime Arrays. However, the main focus is to introduce the reader to a concept called Extremely Sparse Co-Prime Arrays and Samplers. Adjustable pivot co-prime arrays are a special case of the extremely sparse arrays. The closed-form expressions for the weight function and the correlogram bias window are derived. Low latency estimation is demonstrated using the extremely sparse co-prime scheme. Furthermore, a multidimensional and a hybrid extremely sparse co-prime array is proposed as a straightforward extension of the 1D-theory. Finally, a generalized extremely sparse structure is developed with several design parameters. Most existing structures may be viewed as a special case of the generalized scheme.


title: Reachable Polyhedral Marching (RPM): A Safety Verification Algorithm for Robotic Systems with Deep Neural Network Components

abstract: We present a method for computing exact reachable sets for deep neural networks with rectified linear unit (ReLU) activation. Our method is well-suited for use in rigorous safety analysis of robotic perception and control systems with deep neural network components. Our algorithm can compute both forward and backward reachable sets for a ReLU network iterated over multiple time steps, as would be found in a perception-action loop in a robotic system. Our algorithm is unique in that it builds the reachable sets by incrementally enumerating polyhedral cells in the input space, rather than iterating layer-by-layer through the network as in other methods. If an unsafe cell is found, our algorithm can return this result without completing the full reachability computation, thus giving an anytime property that accelerates safety verification. In addition, our method requires less memory during execution compared to existing methods where memory can be a limiting factor. We demonstrate our algorithm on safety verification of the ACAS Xu aircraft advisory system. We find unsafe actions many times faster than the fastest existing method and certify no unsafe actions exist in about twice the time of the existing method. We also compute forward and backward reachable sets for a learned model of pendulum dynamics over a 50 time step horizon in 87s on a laptop computer. Algorithm source code: https://github.com/StanfordMSL/Neural-Network-Reach.


title: Transfer Learning for Oral Cancer Detection using Microscopic Images

abstract: Oral cancer has more than 83% survival rate if detected in its early stages, however, only 29% of cases are currently detected early. Deep learning techniques can detect patterns of oral cancer cells and can aid in its early detection. In this work, we present the first results of neural networks for oral cancer detection using microscopic images. We compare numerous state-of-the-art models via transfer learning approach and collect and release an augmented dataset of high-quality microscopic images of oral cancer. We present a comprehensive study of different models and report their performance on this type of data. Overall, we obtain a 10-15% absolute improvement with transfer learning methods compared to a simple Convolutional Neural Network baseline. Ablation studies show the added benefit of data augmentation techniques with finetuning for this task.


title: FERN: Fair Team Formation for Mutually Beneficial Collaborative Learning

abstract: Automated Team Formation is becoming increasingly important for a plethora of applications in open source community projects, remote working platforms, as well as online educational systems. The latter case, in particular, poses significant challenges that are specific to the educational domain. Indeed, teaming students aims to accomplish far more than the successful completion of a specific task. It needs to ensure that all members in the team benefit from the collaborative work, while also ensuring that the participants are not discriminated with respect to their protected attributes, such as race and gender. Towards achieving these goals, this work introduces FERN, a fair team formation approach that promotes mutually beneficial peer learning, dictated by protected group fairness as equality of opportunity in collaborative learning. We formulate the problem as a multi-objective discrete optimization problem. We show this problem to be NP-hard and propose a heuristic hill-climbing algorithm. Extensive experiments on both synthetic and real-world datasets against well-known team formation techniques show the effectiveness of the proposed method.


title: Multifold behavior of the information transmission by the quantum 3-switch

abstract: We uncover new behaviors of the transmission of information by three quantum channels in superposition of causal orders subject to some level of noise. We find that the transmission can exhibit three different behaviors as the level of noise is varied. This multifold behavior can be explained by the different equivalence classes of quantum switch matrices related to specific combinations of causal orders. We classify these matrices using their characteristic polynomials and matrix invariants, and we calculate analytical expressions for the Holevo information in three representative cases. Our results are a step forward to understand and harness quantum control of causal orders with different levels of noise. We also study the Holevo information as function of a continuous order parameter and analyse transitions at integer values.


title: Strong detection of the CMB lensingxgalaxy weak lensingcross-correlation from ACT-DR4,PlanckLegacy and KiDS-1000

abstract: We measure the cross-correlation between galaxy weak lensing data from the Kilo Degree Survey (KiDS-1000, DR4) and cosmic microwave background (CMB) lensing data from the Atacama Cosmology Telescope (ACT, DR4) and the Planck Legacy survey. We use two samples of source galaxies, selected with photometric redshifts, $(0.1<z_{\rm B}<1.2)$ and $(1.2<z_{\rm B}<2)$, which produce a combined detection significance of the CMB lensing/weak galaxy lensing cross-spectrum of $7.7\sigma$. With the lower redshift galaxy sample, for which the cross-correlation is detected at a significance of $5.3\sigma$, we present joint cosmological constraints on the matter density parameter, $\Omega_{\rm m}$, and the matter fluctuation amplitude parameter, $\sigma_8$, marginalising over three nuisance parameters that model our uncertainty in the redshift and shear calibration, and the intrinsic alignment of galaxies. We find our measurement to be consistent with the best-fitting flat $\Lambda$CDM cosmological models from both Planck and KiDS-1000. We demonstrate the capacity of CMB-weak lensing cross-correlations to set constraints on either the redshift or shear calibration, by analysing a previously unused high-redshift KiDS galaxy sample $(1.2<z_{\rm B}<2)$, with the cross-correlation detected at a significance of $7\sigma$. This analysis provides an independent assessment for the accuracy of redshift measurements in a regime that is challenging to calibrate directly owing to known incompleteness in spectroscopic surveys.


title: A framework for constrained static state estimation in unbalanced distribution networks

abstract: State estimation plays a key role in the transition from the passive to the active operation of distribution systems, as it allows to monitor these networks and, successively, to perform control actions. However, designing state estimators for distribution systems carries a significant amount of challenges. This is due to the physical complexity of the networks, e.g., phase unbalance, and limited measurements. Furthermore, the features of the distribution system present significant local variations, e.g., voltage level and number and type of customers, which makes it hard to design a "one-size-fits-all" state estimator. The present paper introduces a unifying framework that allows to easily implement and compare diverse unbalanced static state estimation models. This is achieved by formulating state estimation as a general constrained optimization problem. The advantages of this approach are described and supported by numerical illustration on a large set of distribution feeders. The framework is also implemented and made available open-source.


title: Femtoscopy of the Origin of the Nucleon Mass

abstract: We study the prospects of using femtoscopic low-momentum correlation measurements at the Large Hadron Collider to access properties of the J/psi-nucleon interaction. The QCD multipole expansion in terms of the J/psi chromopolarizability relates the forward scattering amplitude to a key matrix element to the origin of the nucleon mass problem, the average chromoelectric gluon distribution in the nucleon. We use information on the J/psi-nucleon interaction provided by lattice QCD simulations and phenomenological models to compute J/psi-nucleon correlation functions. The computed correlation functions show clear sensitivity to the interaction, in particular to the J/psi chromopolarizability.


title: The Field Substellar Mass Function Based on the Full-sky 20-pc Census of 525 L, T, and Y Dwarfs

abstract: We present final Spitzer trigonometric parallaxes for 361 L, T, and Y dwarfs. We combine these with prior studies to build a list of 525 known L, T, and Y dwarfs within 20 pc of the Sun, 38 of which are presented here for the first time. Using published photometry and spectroscopy as well as our own follow-up, we present an array of color-magnitude and color-color diagrams to further characterize census members, and we provide polynomial fits to the bulk trends. Using these characterizations, we assign each object a $T_{\rm eff}$ value and judge sample completeness over bins of $T_{\rm eff}$ and spectral type. Except for types $\ge$ T8 and $T_{\rm eff} <$ 600K, our census is statistically complete to the 20-pc limit. We compare our measured space densities to simulated density distributions and find that the best fit is a power law ($dN/dM \propto M^{-\alpha}$) with $\alpha = 0.6{\pm}0.1$. We find that the evolutionary models of Saumon & Marley correctly predict the observed magnitude of the space density spike seen at 1200K $< T_{\rm eff} <$ 1350K, believed to be caused by an increase in the cooling timescale across the L/T transition. Defining the low-mass terminus using this sample requires a more statistically robust and complete sample of dwarfs $\ge$Y0.5 and with $T_{\rm eff} <$ 400K. We conclude that such frigid objects must exist in substantial numbers, despite the fact that few have so far been identified, and we discuss possible reasons why they have largely eluded detection.


title: Die Hard Holographic Phenomenology of Cuprates

abstract: This note discusses the attempts of fitting a number of the approximate power-law dependencies observed in the cuprates into one consistent holographic or holographically inspired hydrodynamic framework. Contrary to the expectations, the goal of reproducing as many as possible of the established behaviors of the thermodynamic and transport coefficients appears to be achievable within the simplistic picture of a non-degenerate fermion fluid with quadratic dispersion. While not immediately elucidating the essential physics of the cuprates, this observation suggests a possible reason for which the previous attempts towards that goal have remained inconclusive.


title: Analysis of the Harmonic Function Used in Bin-Packing

abstract: The harmonic function was first introduced by Lee and Lee (JACM 1985) for analyzing their online bin-packing algorithm. Subsequently, it has been used to obtain approximation algorithms for many different packing problems. Here we slightly generalize the harmonic function and give alternative proofs of its important properties.


title: Neural collapse with unconstrained features

abstract: Neural collapse is an emergent phenomenon in deep learning that was recently discovered by Papyan, Han and Donoho. We propose a simple "unconstrained features model" in which neural collapse also emerges empirically. By studying this model, we provide some explanation for the emergence of neural collapse in terms of the landscape of empirical risk.


title: Many-Body Quantum Zeno Effect and Measurement-Induced Subradiance Transition

abstract: It is well known that by repeatedly measuring a quantum system it is possible to completely freeze its dynamics into a well defined state, a signature of the quantum Zeno effect. Here we show that for a many-body system evolving under competing unitary evolution and variable-strength measurements the onset of the Zeno effect takes the form of a sharp phase transition. Using the Quantum Ising chain with continuous monitoring of the transverse magnetization as paradigmatic example we show that for weak measurements the entanglement produced by the unitary dynamics remains protected, and actually enhanced by the monitoring, while only above a certain threshold the system is sharply brought into an uncorrelated Zeno state. We show that this transition is invisible to the average dynamics, but encoded in the rare fluctuations of the stochastic measurement process, which we show to be perfectly captured by a non-Hermitian Hamiltonian which takes the form of a Quantum Ising model in an imaginary valued transverse field. We provide analytical results based on the fermionization of the non-Hermitian Hamiltonian in supports of our exact numerical calculations.


title: Quantum versus thermal fluctuations in the fcc antiferromagnet: alternative routes to order by disorder

abstract: In frustrated magnetic systems with competing interactions fluctuations can lift the residual accidental degeneracy. We argue that the state selection may have different outcomes for quantum and thermal order by disorder. As an example, we consider the semiclassical Heisenberg fcc antiferromagnet with only the nearest-neighbor interactions. Zero-point oscillations select the type 3 collinear antiferromagnetic state at T=0. Thermal fluctuations favor instead the type 1 antiferromagnetic structure. The opposite tendencies result in a finite-temperature transition between the two collinear states. Competition between effects of quantum and thermal order by disorder is a general phenomenon and is also realized in the J1-J2 square-lattice antiferromagnet at the critical point J2 = 0.5 J1.


title: On the Large Charge Sector in the Critical $O(N)$ Model at Large $N$

abstract: We study operators in the rank-$j$ totally symmetric representation of $O(N)$ in the critical $O(N)$ model in arbitrary dimension $d$, in the limit of large $N$ and large charge $j$ with $j/N\equiv \hat{j}$ fixed. The scaling dimensions of the operators in this limit may be obtained by a semiclassical saddle point calculation. Using the standard Hubbard-Stratonovich description of the critical $O(N)$ model at large $N$, we solve the relevant saddle point equation and determine the scaling dimensions as a function of $d$ and $\hat{j}$, finding agreement with all existing results in various limits. In $4<d<6$, we observe that the scaling dimension of the large charge operators becomes complex above a critical value of the ratio $j/N$, signaling an instability of the theory in that range of $d$. Finally, we also derive results for the correlation functions involving two "heavy" and one or two "light" operators. In particular, we determine the form of the "heavy-heavy-light" OPE coefficients as a function of the charges and $d$.


title: Left orderability of cyclic branched covers of rational knots $C(2n+1,2m,2)$

abstract: We compute the nonabelian $\mathrm{SL_2}(\mathbb{C})$-character varieties of the rational knots $C(2n+1,2m,2)$ in the Conway notation, where $m$ and $n$ are non-zero integers. By studying real points on these varieties, we determine the left orderability of the fundamental groups of the cyclic branched covers of $C(2n+1,2m,2)$.


title: Spatiotemporal dynamics of particle collisions in quantum spin chains

abstract: Recent developments have highlighted the potential of quantum spin models to realize the phenomenology of confinement leading to the formation of bound states such as mesons. In this work we show that Ising chains also provide a platform to realize and probe particle collisions in pristine form with the key advantage that one can not only monitor the asymptotic particle production, but also the whole spatiotemporal dynamics of the collision event. We study both elastic and inelastic collisions between different kinds of mesons and also more complex bound states of mesons, which one can interpret as an analog of exotic particles such as the tetraquark in quantum chromodynamics. We argue that our results not only apply to the specific studied spin model, but can be readily extended to lattice gauge theories in a more general context. As the considered Ising chains admit a natural realization in various quantum simulator platforms, it is a key implication of this work that particle collisions therefore become amenable within current experimental scope. Concretely, we discuss a potentially feasible implementation in systems of Rydberg atoms.

categories: cond-mat.quant-gas hep-lat hep-th quant-ph


title: Non-Conformality, Subregion Complexity and Meson Binding

abstract: We study holographically the zero and finite temperature behavior of the potential energy and holographic subregion complexity corresponding to a probe meson in a non-conformal model. We observe that in zero and low temperature non-conformality has a decreasing effect on the dimensionless meson potential energy. However, non-conformal corrections increase absolute value of the dimensionless holographic subregion complexity in both zero and finite temperature which means the non-conformal state needs less information to be specified. In other words, considering the effect of non-conformality, the less bounded meson state needs less information to be specified. In low temperature limits, thermal corrections decrease meson potential energy and do not have a specific effect on holographic subregion complexity. We find that in the vicinity of the phase transition, the zero temperature meson state is more favorable than the finite temperature state, from the holographic subregion complexity point of view.


title: A theoretical framework of BL Her stars. I. Effect of metallicity and convection parameters on period-luminosity and period-radius relations

abstract: We present a new grid of convective BL Herculis models using the state-of-the-art 1D non-linear radial stellar pulsation tool MESA-RSP. We investigate the impact of metallicity and four sets of different convection parameters on multi-wavelength properties. Non-linear models were computed for periods typical for BL Her stars, i.e. $1 \leq \mathrm{P (days)} \leq 4$ covering a wide range of input parameters - metallicity ($-2.0\; \mathrm{dex} \leq \mathrm{[Fe/H]} \leq 0.0\; \mathrm{dex}$), stellar mass (0.5M$_{\odot}$-0.8M$_{\odot}$), luminosity (50L$_{\odot}$-300L$_{\odot}$) and effective temperature (full extent of the instability strip; in steps of 50K). The total number of BL Her models with full-amplitude stable pulsations used in this study is 10280 across the four sets of convection parameters. We obtain their multiband ($UBVRIJHKLL'M$) light curves and derive new theoretical period-luminosity ($PL$), period-Wesenheit ($PW$) and period-radius ($PR$) relations at mean light. We find that the models computed with radiative cooling show statistically similar slopes for $PL$, $PW$ and $PR$ relations. Most empirical relations match well with the theoretical $PL$, $PW$ and $PR$ relations from the BL Her models computed using the four sets of convection parameters. However, $PL$ slopes of the models with radiative cooling provide a better match to empirical relations for BL Her stars in the LMC in the $HK_S$ bands. For each set of convection parameters, the effect of metallicity is significant in $U$ and $B$-bands and negligible in infrared bands, which is consistent with empirical results. No significant metallicity effects are seen in the $PR$ relations.


title: Cycle-consistent Generative Adversarial Networks for Neural Style Transfer using data from Chang'E-4

abstract: Generative Adversarial Networks (GANs) have had tremendous applications in Computer Vision. Yet, in the context of space science and planetary exploration the door is open for major advances. We introduce tools to handle planetary data from the mission Chang'E-4 and present a framework for Neural Style Transfer using Cycle-consistency from rendered images. The experiments are conducted in the context of the Iris Lunar Rover, a nano-rover that will be deployed in lunar terrain in 2021 as the flagship of Carnegie Mellon, being the first unmanned rover of America to be on the Moon.


title: Canonical geometrization of orientable $3$-manifolds defined by vector-colourings of $3$-polytopes

abstract: In short geometrization conjecture of W.\,Thurston (finally proved by G.~Perelman) says that any oriented $3$-manifold can be canonically partitioned into pieces, which have a geometric structure of one of the eight types. In the seminal paper (1991) M.\,W.\,Davis and T.\,Januszkiewicz introduced a wide class of $n$-dimensional manifolds -- small covers over simple $n$-polytopes. We give a complete answer to the following problem: to build an explicit canonical decomposition for any orientable $3$-manifold defined by a vector-colouring of a simple $3$-polytope, in particular for a small cover. The proof is based on analysis of results in this direction obtained before by different authors.


title: Bulge formation through disc instability -- I

abstract: We use simulations to study the growth of a pseudobulge in an isolated thin exponential stellar disc embedded in a static spherical halo. We observe a transition from later to earlier morphological types and an increase in bar prominence for higher disc-to-halo mass ratios, for lower disc-to-halo size ratios, and for lower halo concentrations. We compute bulge-to-total stellar mass ratios $B/T$ by fitting a two-component S\'ersic-exponential surface-density distribution. The final $B/T$ is strongly related to the disc's fractional contribution $f_{\rm d}$ to the total gravitational acceleration at the optical radius. The formula $B/T=0.5\,f_{\rm }^{1.8}$ fits the simulations to an accuracy of $30\%$, is consistent with observational measurements of B/T and f_d as a function of luminosity, and reproduces the observed relation between $B/T$ and stellar mass when incorporated into the GalICS~2.0 semi-analytic model of galaxy formation.


title: Betrayed by Motion: Camouflaged Object Discovery via Motion Segmentation

abstract: The objective of this paper is to design a computational architecture that discovers camouflaged objects in videos, specifically by exploiting motion information to perform object segmentation. We make the following three contributions: (i) We propose a novel architecture that consists of two essential components for breaking camouflage, namely, a differentiable registration module to align consecutive frames based on the background, which effectively emphasises the object boundary in the difference image, and a motion segmentation module with memory that discovers the moving objects, while maintaining the object permanence even when motion is absent at some point. (ii) We collect the first large-scale Moving Camouflaged Animals (MoCA) video dataset, which consists of over 140 clips across a diverse range of animals (67 categories). (iii) We demonstrate the effectiveness of the proposed model on MoCA, and achieve competitive performance on the unsupervised segmentation protocol on DAVIS2016 by only relying on motion.


title: Explainable Multivariate Time Series Classification: A Deep Neural Network Which Learns To Attend To Important Variables As Well As Informative Time Intervals

abstract: Time series data is prevalent in a wide variety of real-world applications and it calls for trustworthy and explainable models for people to understand and fully trust decisions made by AI solutions. We consider the problem of building explainable classifiers from multi-variate time series data. A key criterion to understand such predictive models involves elucidating and quantifying the contribution of time varying input variables to the classification. Hence, we introduce a novel, modular, convolution-based feature extraction and attention mechanism that simultaneously identifies the variables as well as time intervals which determine the classifier output. We present results of extensive experiments with several benchmark data sets that show that the proposed method outperforms the state-of-the-art baseline methods on multi-variate time series classification task. The results of our case studies demonstrate that the variables and time intervals identified by the proposed method make sense relative to available domain knowledge.


title: MacLeR: Machine Learning-based Run-Time Hardware Trojan Detection in Resource-Constrained IoT Edge Devices

abstract: Traditional learning-based approaches for run-time Hardware Trojan detection require complex and expensive on-chip data acquisition frameworks and thus incur high area and power overhead. To address these challenges, we propose to leverage the power correlation between the executing instructions of a microprocessor to establish a machine learning-based run-time Hardware Trojan (HT) detection framework, called MacLeR. To reduce the overhead of data acquisition, we propose a single power-port current acquisition block using current sensors in time-division multiplexing, which increases accuracy while incurring reduced area overhead. We have implemented a practical solution by analyzing multiple HT benchmarks inserted in the RTL of a system-on-chip (SoC) consisting of four LEON3 processors integrated with other IPs like vga_lcd, RSA, AES, Ethernet, and memory controllers. Our experimental results show that compared to state-of-the-art HT detection techniques, MacLeR achieves 10\% better HT detection accuracy (i.e., 96.256%) while incurring a 7x reduction in area and power overhead (i.e., 0.025% of the area of the SoC and <0.07% of the power of the SoC). In addition, we also analyze the impact of process variation and aging on the extracted power profiles and the HT detection accuracy of MacLeR. Our analysis shows that variations in fine-grained power profiles due to the HTs are significantly higher compared to the variations in fine-grained power profiles caused by the process variations (PV) and aging effects. Moreover, our analysis demonstrates that, on average, the HT detection accuracy drop in MacLeR is less than 1% and 9% when considering only PV and PV with worst-case aging, respectively, which is ~10x less than in the case of the state-of-the-art ML-based HT detection technique.


title: Proportional integral derivative, modeling and ways of stabilization for the spark plasma sintering process

abstract: The stability ofthe proportional--integral--derivative (PID)controlof temperature in the spark plasma sintering (SPS) process is investigated.ThePID regulationsof this process are tested fordifferent SPS toolingdimensions, physical parameters conditions,andareas of temperature control. It isshown thatthe PID regulation quality strongly depends on the heating time lag between the area of heat generation and the area of the temperature control. Tooling temperature rate maps arestudied to revealpotential areas forhighlyefficientPID control.The convergence of the model and experiment indicatesthat even with non-optimal initial PIDcoefficients, it is possible to reduce the temperature regulation inaccuracy to less than 4K by positioning the temperature control location in highlyresponsiveareas revealed by the finite-element calculationsof the temperature spatial distribution.


title: Master sintering curve with dissimilar grain growth trajectories: A case study on MgAl2O4

abstract: Sintering is a key step in the processing of high performance ceramics. Both the density and the grain size play a crucial role on the ceramic sintering kinetics and the final material properties. The master sintering curve (MSC) is a well-known tool for exploring sintering models kinetics. However, the conventional MSC theory assumes a unique sintering trajectory, while our study on MgAl2O4 spinel shows dissimilar growth response. Park's MSC theory has been applied and compared with the conventional MSC approach for obtaining the activation energy with and without dissimilar grain growth trajectories.


title: An elastic framework for ensemble-based large-scale data assimilation

abstract: Prediction of chaotic systems relies on a floating fusion of sensor data (observations) with a numerical model to decide on a good system trajectory and to compensate nonlinear feedback effects. Ensemble-based data assimilation (DA) is a major method for this concern depending on propagating an ensemble of perturbed model realizations.In this paper we develop an elastic, online, fault-tolerant and modular framework called Melissa-DA for large-scale ensemble-based DA. Melissa-DA allows elastic addition or removal of compute resources for state propagation at runtime. Dynamic load balancing based on list scheduling ensuresefficient execution. Online processing of the data produced by ensemble members enables to avoid the I/O bottleneck of file-based approaches. Our implementation embeds the PDAF parallel DA engine, enabling the use of various DA methods. Melissa-DA can support extra ensemble-based DAmethods by implementing the transformation of member background states into analysis states. Experiments confirm the excellent scalability of Melissa-DA, running on up to 16,240 cores, to propagate 16,384 members for a regional hydrological critical zone assimilation relying on theParFlow model on a domain with about 4 M grid cells.


title: Blade Envelopes Part I: Concept and Methodology

abstract: Blades manufactured through flank and point milling will likely exhibit geometric variability. Gauging the aerodynamic repercussions of such variability, prior to manufacturing a component, is challenging enough, let alone trying to predict what the amplified impact of any in-service degradation will be. While rules of thumb that govern the tolerance band can be devised based on expected boundary layer characteristics at known regions and levels of degradation, it remains a challenge to translate these insights into quantitative bounds for manufacturing. In this work, we tackle this challenge by leveraging ideas from dimension reduction to construct low-dimensional representations of aerodynamic performance metrics. These low-dimensional models can identify a subspace which contains designs that are invariant in performance -- the inactive subspace. By sampling within this subspace, we design techniques for drafting manufacturing tolerances and for quantifying whether a scanned component should be used or scrapped. We introduce the blade envelope as a computational manufacturing guide for a blade that is also amenable to qualitative visualizations. In this paper, the first of two parts, we discuss its underlying concept and detail its computational methodology, assuming one is interested only in the single objective of ensuring that the loss of all manufactured blades remains constant. To demonstrate the utility of our ideas we devise a series of computational experiments with the Von Karman Institute's LS89 turbine blade.


title: Nudge Attacks on Point-Cloud DNNs

abstract: The wide adaption of 3D point-cloud data in safety-critical applications such as autonomous driving makes adversarial samples a real threat. Existing adversarial attacks on point clouds achieve high success rates but modify a large number of points, which is usually difficult to do in real-life scenarios. In this paper, we explore a family of attacks that only perturb a few points of an input point cloud, and name them nudge attacks. We demonstrate that nudge attacks can successfully flip the results of modern point-cloud DNNs. We present two variants, gradient-based and decision-based, showing their effectiveness in white-box and grey-box scenarios. Our extensive experiments show nudge attacks are effective at generating both targeted and untargeted adversarial point clouds, by changing a few points or even a single point from the entire point-cloud input. We find that with a single point we can reliably thwart predictions in 12--80% of cases, whereas 10 points allow us to further increase this to 37--95%. Finally, we discuss the possible defenses against such attacks, and explore their limitations.


title: Fungal bioremediation of diuron-contaminated waters: evaluation of its degradation and the effect of amendable factors on its removal in a trickle-bed reactor under non-sterile conditions

abstract: The occurrence of the extensively used herbicide diuron in the environment poses a severe threat to the ecosystem and human health. Four different ligninolytic fungi were studied as biodegradation candidates for the removal of diuron. Among them, T. versicolor was the most effective species, degrading rapidly not only diuron (83%) but also the major metabolite 3,4-dichloroaniline (100%), after 7-day incubation. During diuron degradation, five transformation products (TPs) were found to be formed and the structures for three of them are tentatively proposed. According to the identified TPs, a hydroxylated intermediate 3-(3,4-dichlorophenyl)-1-hydroxymethyl-1-methylurea (DCPHMU) was further metabolized into the N-dealkylated compounds 3-(3,4-dichlorophenyl)-1-methylurea (DCPMU) and 3,4-dichlorophenylurea (DCPU). The discovery of DCPHMU suggests a relevant role of hydroxylation for subsequent N-demethylation, helping to better understand the main reaction mechanisms of diuron detoxification. Experiments also evidenced that degradation reactions may occur intracellularly and be catalyzed by the cytochrome P450 system. A response surface method, established by central composite design, assisted in evaluating the effect of operational variables in a trickle-bed bioreactor immobilized with T. versicolor on diuron removal. The best performance was obtained at low recycling ratios and influent flow rates. Furthermore, results indicate that the contact time between the contaminant and immobilized fungi plays a crucial role in diuron removal. This study represents a pioneering step forward amid techniques for bioremediation of pesticides-contaminated waters using fungal reactors at a real scale.


title: Improved fully automated method for the determination of medium to highly polar pesticides in surface and groundwater and application in two distinct agriculture-impacted areas

abstract: Water is an essential resource for all living organisms. The continuous and increasing use of pesticides in agricultural and urban activities results in the pollution of water resources and represents an environmental risk. To control and reduce pesticide pollution, reliable multi-residue methods for the detection of these compounds in water are needed. In this context, the present work aimed at providing an analytical method for the simultaneous determination of trace levels of 51 target pesticides in water and applying it to the investigation of target pesticides in two agriculture-impacted areas of interest. The method developed, based on an isotopic dilution approach and on-line solid-phase extraction-liquid chromatography-tandem mass spectrometry, is fast, simple, and to a large extent automated, and allows the analysis of most of the target compounds in compliance with European regulations. Further application of the method to the analysis of selected water samples collected at the lowest stretches of the two largest river basins of Catalonia (NE Spain), Llobregat and Ter, revealed the presence of a wide suite of pesticides, and some of them at concentrations above the water quality standards (irgarol and dichlorvos) or the acceptable method detection limits (methiocarb, imidacloprid, and thiacloprid), in the Llobregat, and much cleaner waters in the Ter River basin. Risk assessment of the pesticide concentrations measured in the Llobregat indicated high risk due to the presence of irgarol, dichlorvos, methiocarb, azinphos ethyl, imidacloprid, and diflufenican (hazard quotient (HQ) values>10), and an only moderate potential risk in the Ter River associated to the occurrence of bentazone and irgarol (HQ>1).


title: Heterogeneous Multipartite Entanglement Purification for Size-Constrained Quantum Devices

abstract: The entanglement resource required for quantum information processing comes in a variety of forms, from Bell states to multipartite GHZ states or cluster states. Purifying these resources after their imperfect generation is an indispensable step towards using them in quantum architectures. While this challenge, both in the case of Bell pairs and more general multipartite entangled states, is mostly overcome in the presence of perfect local quantum hardware with unconstrained qubit register sizes, devising optimal purification strategies for finite-size realistic noisy hardware has remained elusive. Here we depart from the typical purification paradigm for multipartite states explored in the last twenty years. We present cases where the hardware limitations are taken into account, and surprisingly find that smaller `sacrificial' states, like Bell pairs, can be more useful in the purification of multipartite states than additional copies of these same states. This drastically simplifies the requirements and presents a fundamentally new pathway to leverage near term networked quantum hardware.


title: Inconsistencies arising from the coupling of galaxy formation sub-grid models to Pressure-Smoothed Particle Hydrodynamics

abstract: Smoothed Particle Hydrodynamics (SPH) is a Lagrangian method for solving the fluid equations that is commonplace in astrophysics, prized for its natural adaptivity and stability. The choice of variable to smooth in SPH has been the topic of contention, with smoothed pressure (P-SPH) being introduced to reduce errors at contact discontinuities relative to smoothed density schemes. Smoothed pressure schemes produce excellent results in isolated hydrodynamics tests; in more complex situations however, especially when coupling to the `sub-grid' physics and multiple time-stepping used in many state-of-the-art astrophysics simulations, these schemes produce large force errors that can easily evade detection as they do not manifest as energy non-conservation. Here two scenarios are evaluated: the injection of energy into the fluid (common for stellar feedback) and radiative cooling. In the former scenario, force and energy conservation errors manifest (of the same order as the injected energy), and in the latter large force errors that change rapidly over a few timesteps lead to instability in the fluid (of the same order as the energy lost to cooling). Potential ways to remedy these issues are explored with solutions generally leading to large increases in computational cost. Schemes using a Density-based formulation do not create these instabilities and as such it is recommended that they are preferred over Pressure-based solutions when combined with an energy diffusion term to reduce errors at contact discontinuities.


title: Measuring the vertical response of the Galactic disc to an infalling satellite

abstract: Using N-body simulations of the Milky Way interacting with a satellite similar to the Sagittarius dwarf galaxy, we quantitatively analyse the vertical response of the Galactic disc to the satellite's repeated impacts. We approximate the vertical distortion of the Galactic disc as the sum of the first three Fourier azimuthal terms m = 0, 1 and 2, and observe their evolution in different dynamical regimes of interaction. After the first interaction, the m=0 term manifests itself as outgoing ring-like vertical distortions. The m=1 term (S-shape warp) is prograde when the impacts of the satellite are more frequent, or in general close to an interaction, whereas it is slowly retrograde in the most quiescent phases. The m=2 term is typically prograde, and close to an interaction it couples with the m=1 term. Finally, we find that the vertical response of the disc can be recovered in an unbiased way using the instantaneous positions and velocities of stars in a limited volume of the Galactic disc, analogous to real data, and that the measured vertical pattern speeds have a constraining power in the context of a Milky Way-satellite interaction.


title: Exploring the epoch of hydrogen reionization using FRBs

abstract: We describe three different methods for exploring the hydrogen reionization epoch using fast radio bursts (FRBs) and provide arguments for the existence of FRBs at high redshift (z). The simplest way, observationally, is to determine the maximum dispersion measure (DM$_{\rm max}$) of FRBs for an ensemble that includes bursts during the reionization. The DM$_{\rm max}$ provides information regarding reionization much like the optical depth of the CMB to Thomson scattering does, and it has the potential to be more accurate than constraints from Planck, if DM$_{\rm max}$ can be measured to a precision better than 500 $\mbox{pc cm}^{-3}$. Another method is to measure redshifts of about 40 FRBs between z of 6-10 with$\sim10\%$ accuracy to obtain the average electron density in 4 different z-bins with $\sim4\%$ accuracy. These two methods don't require knowledge of the FRB luminosity function and its possible redshift evolution. Finally, we show that the reionization history is reflected in the number of FRBs per unit DM, given a fluence limited survey of FRBs that includes bursts during the reionization epoch; we show using FIRE simulations that the contributions to DM from the FRB host galaxy $\&$ CGM during the reionization era is a small fraction of the observed DM. This third method requires no redshift information but does require knowledge of the FRB luminosity function.


title: Purification and Entanglement Routing on Quantum Networks

abstract: We present an approach to purification and entanglement routing on complex quantum network architectures, that is, how a quantum network equipped with imperfect channel fidelities and limited memory storage time can distribute entanglement between users. We explore how network parameters influence the performance of path-finding algorithms necessary for optimizing routing and, in particular, we explore the interplay between the bandwidth of a quantum channels and the choice of purification protocol. Finally, we demonstrate multi-path routing on various network topologies with resource constraints, in an effort to inform future design choices for quantum network configurations. Our work optimizes both the choice of path over the quantum network and the choice of purification schemes used between nodes. We consider not only pair-production rate, but optimize over the fidelity of the delivered entangled state. We introduce effective heuristics enabling fast path-finding algorithms for maximizing entanglement shared between two nodes on a quantum network, with performance comparable to that of a computationally-expensive brute-force path search.


title: Eppur \`e piatto? The cosmic chronometer take on spatial curvature and cosmic concordance

abstract: The question of whether Cosmic Microwave Background (CMB) temperature and polarization data from Planck favor a spatially closed Universe with curvature parameter $\Omega_K<0$ has been the subject of recent intense discussions. Attempts to break the geometrical degeneracy combining Planck data with external datasets such as Baryon Acoustic Oscillation (BAO) measurements all point towards a spatially flat Universe, at the cost of significant tensions with Planck, which make the resulting dataset combination problematic. Settling this issue requires identifying a dataset which can break the geometrical degeneracy while not incurring in these tensions. We argue that cosmic chronometers (CC), measurements of the expansion rate $H(z)$ from the relative ages of massive early-type passively evolving galaxies, are the dataset we are after. Furthermore, CC come with the additional advantage of being virtually free of cosmological model assumptions. Combining Planck 2018 CMB temperature and polarization data with the latest CC measurements, we break the geometrical degeneracy and find $\Omega_K=-0.0054 \pm 0.0055$, consistent with a spatially flat Universe and competitive with the Planck+BAO constraint. Our results are stable against minimal parameter space extensions and CC systematics, and we find no substantial tension between Planck and CC data within a non-flat Universe, making the resulting combination reliable. Our results allow us to assert with confidence that the Universe is spatially flat to the ${\cal O}(10^{-2})$ level, a finding which might possibly settle the ongoing spatial curvature debate, and lends even more support to the already very successful inflationary paradigm.


title: Black hole superradiance of self-interacting scalar fields

abstract: Black hole superradiance is a powerful probe of light, weakly-coupled hidden sector particles. Many candidate particles, such as axions, generically have self-interactions that can influence the evolution of the superradiant instability. As pointed out in arXiv:1604.06422 in the context of a toy model, much of the existing literature on spin-0 superradiance does not take into account the most important self-interaction-induced processes. These processes lead to energy exchange between quasi-bound levels and particle emission to infinity; for large self-couplings, superradiant growth is saturated at a quasi-equilibrium configuration of reduced level occupation numbers. In this paper, we perform a detailed analysis of the rich dynamics of spin-0 superradiance with self-interactions, and the resulting observational signatures. We focus on quartic self-interactions, which dominate the evolution for most models of interest. We explore multiple distinct regimes of parameter space introduced by a non-zero self-interaction, including the simultaneous population of two or more bound levels; at large coupling, we confirm the basic picture of quasi-equilibrium saturation and provide evidence that the "bosenova" collapse does not occur in most of the astrophysical parameter space. Compared to gravitational superradiance, we find that gravitational wave "annihilation" signals and black hole spin-down are parametrically suppressed with increasing interactions, while new gravitational wave "transition" signals can take place for moderate interactions. The novel phenomenon of scalar wave emission is less suppressed at large couplings, and if the particle has Standard Model interactions, then coherent, monochromatic axion wave signals from black hole superradiance may be detectable in proposed axion dark matter experiments.


title: Disformal Dark Matter

abstract: We generalize dark matter production to a two-metric framework whereby the physical metric, which couples to the Standard Model (SM), is conformally and/or disformally related to the metric governing the gravitational dynamics. We show that this setup is naturally present in many Ultra Violet (UV) constructions, from K\"ahler moduli fields to tensor-portal models, and from emergent gravity to supergravity models. In this setting we study dark matter production in the early Universe resulting from both scatterings off the thermal bath and the radiative decay of the inflaton. We also take into account non-instantaneous reheating effects at the end of inflation. In this context, dark matter emerges from the production of the scalar field mediating the conformal/disformal interactions with the SM, i.e. realising a Feebly Interacting Matter Particle (FIMP) scenario where the suppression scale of the interaction between the scalar and the SM can be taken almost as high as the Planck scale in the deep UV.


title: Unraveling the origin of magnetic fields in galaxies

abstract: Despite their ubiquity, there are many open questions regarding galactic and cosmic magnetic fields. Specifically, current observational constraints cannot rule out if magnetic fields observed in galaxies were generated in the Early Universe or are of astrophysical nature. Motivated by this we use our magnetic tracers algorithm to investigate whether the signatures of primordial magnetic fields persist in galaxies throughout cosmic time. We simulate a Milky Way-like galaxy in four scenarios: magnetised solely by primordial magnetic fields, magnetised exclusively by SN-injected magnetic fields, and two combined primordial + SN magnetisation cases. We find that once primordial magnetic fields with a comoving strength $B_0 >10^{-12}$ G are considered, they remain the primary source of galaxy magnetisation. Our magnetic tracers show that, even combined with galactic sources of magnetisation, when primordial magnetic fields are strong, they source the large-scale fields in the warm metal-poor phase of the simulated galaxy. In this case, the circumgalactic and intergalactic medium can be used to probe $B_0$ without risk of pollution by magnetic fields originated in the galaxy. Furthermore, whether magnetic fields are primordial or astrophysically-sourced can be inferred by studying local gas metallicity. As a result, we predict that future state-of-the-art observational facilities of magnetic fields in galaxies will have the potential to unravel astrophysical and primordial magnetic components of our Universe.


title: Primordial non-Gaussianity from the Effects of the Standard Model Higgs during Reheating after Inflation

abstract: We propose a new way of studying the Higgs potential at extremely high energies. The Standard Model (SM) Higgs boson, as a light spectator field during inflation in the early Universe, can acquire large field values from its quantum fluctuations which vary among different causal (Hubble) patches. Such a space dependence of the Higgs after the end of inflation leads to space-dependent SM particle masses and hence variable efficiency of reheating, when the inflaton decays to Higgsed SM particles. Inhomogeneous reheating results in (observable) temperature anisotropies. Further, the resulting temperature anisotropy spectrum acquires a significant non-Gaussian component, which is constrained by $\textit{Planck}$ observations of the Cosmic Microwave Background (CMB) and potentially detectable in next-generation experiments. Constraints on this non-Gaussian signal largely exclude the possibility of the observed temperature anisotropies arising primarily from Higgs effects. Hence, in principle, observational searches for non-Gaussianity in the CMB can be used to constrain the dynamics of the Higgs boson at very high (inflationary) energies.


title: Formation of Massive Star Clusters by Fast HI Gas Collision

abstract: Young massive clusters (YMCs) are dense aggregates of young stars, which are essential to galaxy evolution, owing to their ultraviolet radiation, stellar winds, and supernovae. The typical mass and radius of YMCs are M~10^4 M_sun and R~1 pc, respectively, indicating that many stars are located in a small region. The formation of YMC precursor clouds may be difficult because a very compact massive cloud should be formed before stellar feedback blows off the cloud. Recent observational studies suggest that YMCs can be formed as a consequence of the fast HI gas collision with a velocity of ~100 km s^-1, which is the typical velocity of the galaxy-galaxy interaction. In this study, we examine whether the fast HI gas collision triggers YMC formation using three-dimensional magnetohydrodynamics simulations, which include the effects of self-gravity, radiative cooling/heating, and chemistry. We demonstrate that massive gravitationally bound gas clumps with M >10^4 M_sun and L~4 pc are formed in the shock compressed region induced by the fast HI gas collision, which massive gas clumps can evolve into YMCs. Our results show that the YMC precursors are formed by the global gravitational collapse of molecular clouds, and YMCs can be formed even in low-metal environments, such as the Magellanic Clouds. Additionally, the very massive YMC precursor cloud, with M >10^5 M_sun, can be created when we consider the fast collision of HI clouds, which may explain the origin of the very massive stellar cluster R136 system in the Large Magellanic Cloud.


title: The axiomatic and the operational approaches to resource theories of magic do not coincide

abstract: Stabiliser operations occupy a prominent role in fault-tolerant quantum computing. They are defined operationally: by the use of Clifford gates, Pauli measurements and classical control. These operations can be efficiently simulated on a classical computer, a result which is known as the Gottesman-Knill theorem. However, an additional supply of magic states is enough to promote them to a universal, fault-tolerant model for quantum computing. To quantify the needed resources in terms of magic states, a resource theory of magic has been developed. Stabiliser operations (SO) are considered free within this theory, however they are not the most general class of free operations. From an axiomatic point of view, these are the completely stabiliser-preserving (CSP) channels, defined as those that preserve the convex hull of stabiliser states. It has been an open problem to decide whether these two definitions lead to the same class of operations. In this work, we answer this question in the negative, by constructing an explicit counter-example. This indicates that recently proposed stabiliser-based simulation techniques of CSP maps are strictly more powerful than Gottesman-Knill-like methods. The result is analogous to a well-known fact in entanglement theory, namely that there is a gap between the operationally defined class of local operations and classical communication (LOCC) and the axiomatically defined class of separable channels.


title: Massless Positivity in Graviton Exchange

abstract: We formulate Positivity Bounds for scattering amplitudes including exchange of massless particles. We generalize the standard construction through dispersion relations to include the presence of a branch cut along the real axis in the complex plane for the Maldestam variable $s$. In general, validity of these bounds require the cancellation of divergences in the forward limit of the amplitude, proportional to $t^{-1}$ and $\log(t)$. We show that this is possible in the case of gravitons if one assumes a Regge behavior of the amplitude at high energies below the Planck scale, as previously suggested in the literature, and that the concrete UV behaviour of the amplitude is uniquely determined by the structure of IR divergences. We thus extend previous results by including a sub-leading logarithmic term, which we show to be universal. The bounds that we present here have the potential of constraining very general models of modified gravity and EFTs of matter coupled to gravitation.


title: Symmetries of stationary points of the $G$-invariant potential and the framework of the auxiliary group

abstract: We classify the constraints on a stationary point of the potential invariant under a finite group into intrinsic and extrinsic based on whether they are independent of the coefficients in the potential or not. We find that the symmetry group of a set of stationary points can be larger than that of the potential and the stabilizer under this group generates intrinsic constraints. By applying these findings in the framework of the auxiliary group, we show that the constraints that can only be obtained extrinsically in an elementary theory can be generated intrinsically in an effective theory.


title: Quantum speedups for convex dynamic programming

abstract: We present a quantum algorithm to solve dynamic programming problems with convex value functions. For linear discrete-time systems with a $d$-dimensional state space of size $N$, the proposed algorithm outputs a quantum-mechanical representation of the value function in time $O(T \gamma^{dT}\mathrm{polylog}(N,(T/\varepsilon)^{d}))$, where $\varepsilon$ is the accuracy of the solution, $T$ is the time horizon, and $\gamma$ is a problem-specific parameter depending on the condition numbers of the cost functions. This allows us to evaluate the value function at any fixed state in time $O(T \gamma^{dT}\sqrt{N}\,\mathrm{polylog}(N,(T/\varepsilon)^{d}))$, and the corresponding optimal action can be recovered by solving a convex program. The class of optimization problems to which our algorithm can be applied includes provably hard stochastic dynamic programs. Finally, we show that the algorithm obtains a quadratic speedup (up to polylogarithmic factors) compared to the classical Bellman approach on some dynamic programs with continuous state space that have $\gamma=1$.


title: On the relation between the $\overline{\mathrm{MS}}$ and the kinetic mass of heavy quarks

abstract: We compute the relation between the pole mass and the kinetic mass of a heavy quark to three loops. Using the known relation between the pole and the $\overline{\rm MS}$ mass we obtain precise conversion relations between the $\overline{\rm MS}$ and kinetic masses. The kinetic mass is defined via the moments of the spectral function for the scattering involving a heavy quark close to threshold. This requires the computation of the imaginary part of a forward scattering amplitude up to three-loop order. We discuss in detail the expansion procedure and the reduction to master integrals. For the latter analytic results are provided. We apply our result both to charm and bottom quark masses. In the latter case we compute and include finite charm quark mass effects. Furthermore, we determine the large-$\beta_0$ result for the conversion formula at four-loop order. For the bottom quark we estimate the uncertainty in the conversion between the $\overline{\rm MS}$ and kinetic masses to about 15 MeV which is an improvement by a factor two to three as compared to the two-loop formula. The improved precision is crucial for the extraction of the Cabibbo-Kobayashi-Maskawa matrix element $|V_{cb}|$ at Belle II.


title: AGN on the move: A search for off-nuclear AGN from recoiling SMBHs and ongoing galaxy mergers with the Zwicky Transient Facility

abstract: A supermassive black hole (SMBH) ejected from the potential well of its host galaxy via gravitational wave recoil carries important information about the mass ratio and spin alignment of the pre-merger SMBH binary. Such a recoiling SMBH may be detectable as an active galactic nucleus (AGN) broad line region offset by up to 10\,kpc from a disturbed host galaxy. We describe a novel methodology using forward modeling with \texttt{The Tractor} to search for such offset AGN in a sample of 5493 optically variable AGN detected with the Zwicky Transient Facility (ZTF). We present the discovery of 9 AGN which may be spatially offset from their host galaxies and are candidates for recoiling SMBHs. Five of these offset AGN exhibit double-peaked broad Balmer lines which may arise from unobscured accretion disk emission and four show radio emission indicative of a relativistic jet. The fraction of double-peaked emitters in our spatially offset AGN sample is significantly larger than the 16\% double-peaked emitter fraction observed for ZTF AGN overall. In our sample of variable AGN we also identified 52 merging galaxies, including a new spectroscopically confirmed dual AGN. Finally, we detected the dramatic rebrightening of SDSS1133, a previously discovered variable object and recoiling SMBH candidate, in ZTF. The flare was accompanied by the re-emergence of strong P-Cygni line features indicating that it may be an outbursting luminous blue variable star.


title: A modular characterization of supersolvable lattices

abstract: We characterize supersolvable lattices in terms of a certain modular type relation. McNamara and Thomas earlier characterized this class of lattices as those graded lattices having a maximal chain that consists of left-modular elements. Our characterization replaces the condition of gradedness with a second modularity condition on the maximal chain of left-modular elements.


title: Kaluza-Klein spectrometry for ${\rm AdS_{3}}$ vacua

abstract: We use exceptional field theory to compute Kaluza-Klein mass spectra around ${\rm AdS_{3}}$ vacua that sit in half-maximal gauged supergravity in three dimensions. The formalism applies to any vacuum that arises from a consistent truncation of higher-dimensional supergravity, no matter what symmetries are preserved. We illustrate its efficiency by computing the spectra of ${\cal N}=(2,0)$ and ${\cal N}=(1,1)$ six-dimensional supergravities on ${\rm AdS_{3}}\times S^{3}$ and of type II supergravity on ${\rm AdS_{3}}\times S^{3}\times S^{3}\times S^{1}$.


title: Heun operator of Lie type and the modified algebraic Bethe ansatz

abstract: The generic Heun operator of Lie type is identified as a certain $BC$-Gaudin magnet Hamiltonian in a magnetic field. By using the modified algebraic Bethe ansatz introduced to diagonalize such Gaudin models, we obtain the spectrum of the generic Heun operator of Lie type in terms of the Bethe roots of inhomogeneous Bethe equations. We show also that these Bethe roots are intimately associated to the roots of polynomial solutions of the differential Heun equation. We illustrate the use of this approach in two contexts: the representation theory of $O(3)$ and the computation of the entanglement entropy for free Fermions on the Krawtchouk chain.


title: Differentially Private Learning Needs Better Features (or Much More Data)

abstract: We demonstrate that differentially private machine learning has not yet reached its "AlexNet moment" on many canonical vision tasks: linear models trained on handcrafted features significantly outperform end-to-end deep neural networks for moderate privacy budgets. To exceed the performance of handcrafted features, we show that private learning requires either much more private data, or access to features learned on public data from a similar domain. Our work introduces simple yet strong baselines for differentially private learning that can inform the evaluation of future progress in this area.


title: Macroscopic Superposition States in Isolated Quantum Systems

abstract: For any choice of initial state and weak assumptions about the Hamiltonian, large isolated quantum systems undergoing Schrodinger evolution spend most of their time in macroscopic superposition states. The result follows from von Neumann's 1929 Quantum Ergodic Theorem. As a specific example, we consider a box containing a solid ball and some gas molecules. Regardless of the initial state, the system will evolve into a quantum superposition of states with the ball in macroscopically different positions. Thus, despite their seeming fragility, macroscopic superposition states are ubiquitous consequences of quantum evolution. We discuss the connection to many worlds quantum mechanics.


title: Spin-orbit-enhanced magnetic surface second-harmonic generation in Sr$_2$IrO$_4$

abstract: An anomalous optical second-harmonic generation (SHG) signal was previously reported in Sr$_2$IrO$_4$ and attributed to a hidden odd-parity bulk magnetic state. Here we investigate the origin of this SHG signal using a combination of bulk magnetic susceptibility, magnetic-field-dependent SHG rotational anisotropy, and overlapping wide-field SHG imaging and atomic force microscopy measurements. We find that the anomalous SHG signal exhibits a two-fold rotational symmetry as a function of in-plane magnetic field orientation that is associated with a crystallographic distortion. We also show a change in SHG signal across step edges that tracks the bulk antiferromagnetic stacking pattern. While we do not rule out the existence of hidden order in Sr$_2$IrO$_4$, our results altogether show that the anomalous SHG signal in parent Sr$_2$IrO$_4$ originates instead from a surface-magnetization-induced electric-dipole process that is enhanced by strong spin-orbit coupling.


title: FOREVER22: galaxy formation in protocluster regions

abstract: We present results from a new cosmological hydrodynamics simulation campaign of protocluster (PC) regions, FOREVER22: FORmation and EVolution of galaxies in Extremely-overdense Regions motivated by SSA22. The simulations cover a wide range of cosmological scales using three different zoom set-ups in a parent volume of $(714.2~\rm cMpc)^{3}$: PCR (Proto-Cluster Region; $V= (28.6~{\rm cMpc})^{3} $, SPH particle mass, $m_{\rm{SPH}} = 4.1 \times 10^{6}~\rm M_{\odot}$ and final redshift, $z_{\rm end}=2.0$), BCG (Brightest proto-Cluster Galaxy; $V \sim (10~{\rm cMpc})^{3} $, $m_{\rm SPH} = 5.0\times10^{5}~\rm M_{\odot}$ and $z_{\rm end}=4.0$ ), and First ( $V \sim (3~{\rm cMpc})^{3} $, $m_{\rm SPH} = 7.9 \times 10^{3}~\rm M_{\odot}$ and $z_{\rm end}=9.5$) runs, that allow to focus on different aspects of galaxy formation. In the PCR runs, we follow 10 PCs, each harbouring 1 - 4 SMBHs with $M_{\rm BH} \ge 10^{9}~\rm M_{\odot}$. One of the PC cores shows a spatially close arrangement of seven starburst galaxies with ${\rm SFR} \gtrsim 100~\rm M_{\odot}~yr^{-1}$ each, that are dust-obscured and would appear as submillimeter galaxies with flux $\gtrsim 1~$ mJy at $1.1~ \rm mm$ in observations. The BCG runs show that the total SFRs of haloes hosting BCGs are affected by AGN feedback, but exceed $1000~\rm M_{\odot}~yr^{-1}$ at $z \lesssim 6$. The First runs resolve mini-haloes hosting population (Pop) III stars and we show that, in PC regions, the dominant stellar population changes from Pop III to Pop II at $z \gtrsim 20$, and the first galaxies with ${\rm SFR} \gtrsim 18~\rm M_{\odot}~yr^{-1}$ form at $z \sim 10$. These can be prime targets for future observations with the James Webb Space Telescope. Our simulations successfully reproduce the global star formation activities in observed PCs and suggest that PCs can kickstart cosmic reionization.


title: Equations of linear subvarieties of strata of differentials

abstract: For a linear subvariety $M$ of a stratum of meromorphic differentials, we investigate its closure in the multi-scale compactification constructed by Bainbridge-Chen-Gendron-Grushevsky-M\"oller. We prove various restrictions on the type of defining linear equations in period coordinates for $M$ near its boundary, and prove that the closure is locally a toric variety. As applications, we give a fundamentally new proof of a generalization of the cylinder deformation theorem of Wright to the case of meromorphic strata, and construct a smooth compactification of the Hurwitz space of covers of the Riemann sphere.


title: Vanishing of Avramov Obstructions for Products of Sequentially Transverse Ideals

abstract: Two ideals $I$ and $J$ are called transverse if $I \cap J = IJ$. We show that the obstructions defined by Avramov for classes of (sequentially) transverse ideals in regular local rings are always $0$. In particular, we compute an explicit free resolution and Koszul homology for all such ideals. Moreover, we construct an explicit trivial Massey operation on the associated Koszul complex and hence (by Golod's construction) a minimal free resolution of the residue field over the quotient defined by the product of transverse ideals. We conclude with questions about the existence of associative multiplicative structures on the minimal free resolution of the quotient defined by products of transverse ideals.


title: The Environment of FRB 121102 and Possible Relation to SGR/PSR J1745-2900

abstract: Variations of the dispersion (DM) and rotation (RM) measures of FRB 121102 indicate magnetic fields $\sim$ 3--17 mG in the dispersing plasma. The electron density may be $\sim 10^4\,$cm$^{-3}$. The observed time scales $\sim 1$ year constrain the size of the plasma cloud. Increasing DM excludes simple models involving an expanding supernova remnant, and the non-zero RM excludes spherical symmetry. The varying DM and RM may be attributable to the motion of plasma into or out of the line of sight or changing electron density within slower-moving plasma. The extraordinarily large RM of FRB 121102 implies an environment, and possibly also a formation process and source, qualitatively different from those of other FRB. The comparable and comparably varying RM of SGR/PSR J1745-2900 suggests it as an FRB candidate. An appendix discusses the age of FRB 121102 in the context of a "Copernican Principle".


title: What is nonlocal in counterfactual quantum communication?

abstract: We revisit the "counterfactual quantum communication" of Salih et al. [1], who claim that an observer "Bob" can send one bit of information to a second observer "Alice" without any physical particle traveling between them. We show that a locally conserved, massless current - specifically, a current of modular angular momentum, $L_z$ mod 2$\hbar$ - carries the one bit of information. We integrate the flux of $L_z$ mod 2$\hbar$ from Bob to Alice and show that it equals one of the two eigenvalues of $L_z$ mod 2$\hbar$, either 0 or $\hbar$, thus precisely accounting for the one bit of information he sends her. We previously [2] obtained this result using weak values of $L_z$ mod 2$\hbar$; here we do not use weak values.


title: State selection in frustrated magnets

abstract: Magnets with frustration often show accidental degeneracies, characterized by a large classical ground-state space (CGSS). Quantum fluctuations may `select' one of these ground states -- a phenomenon labeled `order by (quantum) disorder' in literature. In this article, we examine the mechanism(s) by which such state selection takes place. We argue that a magnet, at low energies, maps to a particle moving on the CGSS. State selection corresponds to localization of the particle at a certain point on this space. We distinguish two mechanisms that can bring about localization. In the first, quantum fluctuations generate a potential on the CGSS. If the potential has a deep enough minimum, then the particle localizes in its vicinity. We denote this as `order by potential' (ObP). In the second scenario, the particle localizes at a self-intersection point due to bound-state formation -- a consequence of geometry and quantum interference. Following recent studies by the present authors, we denote this scenario as `order by singularity' (ObS). We place our discussion within the context of the one-dimensional spin-$S$ Kitaev model. We map out its CGSS which grows systematically with increasing system size. It resembles a network where the number of nodes increases exponentially. In addition, the number of wires that cross at each node also grows exponentially. This self-intersecting structure leads to ObS, with the low-energy physics determined by a small subset of the CGSS, consisting of `Cartesian' states. A contrasting picture emerges when an additional XY antiferromagnetic coupling is introduced. The CGSS simplifies dramatically, taking the form of a circle. Spin-wave fluctuations generate a potential on this space, giving rise to state selection by ObP under certain conditions. Apart from contrasting ObS and ObP, we discuss the possibility of ObS in macroscopic magnets.


title: On piecewise hyperdefinable groups

abstract: The aim of this paper is to generalize and improve two of the main model-theoretic results of "Stable group theory and approximate subgroups" by E. Hrushovski to the context of piecewise hyperdefinable sets. The first one is the existence of Lie models. The second one is the stabilizer theorem. In the process, a systematic study of the structure of piecewise hyperdefinable sets is developed. In particular, we show the most significant properties of their logic topologies.


title: Recognizing Proper Tree-Graphs

abstract: We investigate the parameterized complexity of the recognition problem for the proper $H$-graphs. The $H$-graphs are the intersection graphs of connected subgraphs of a subdivision of a multigraph $H$, and the properness means that the containment relationship between the representations of the vertices is forbidden. The class of $H$-graphs was introduced as a natural (parameterized) generalization of interval and circular-arc graphs by Bir\'o, Hujter, and Tuza in 1992, and the proper $H$-graphs were introduced by Chaplick et al. in WADS 2019 as a generalization of proper interval and circular-arc graphs. For these graph classes, $H$ may be seen as a structural parameter reflecting the distance of a graph to a (proper) interval graph, and as such gained attention as a structural parameter in the design of efficient algorithms. We show the following results. - For a tree $T$ with $t$ nodes, it can be decided in $ 2^{\mathcal{O}(t^2 \log t)} \cdot n^3 $ time, whether an $n$-vertex graph $ G $ is a proper $ T $-graph. For yes-instances, our algorithm outputs a proper $T$-representation. This proves that the recognition problem for proper $H$-graphs, where $H$ required to be a tree, is fixed-parameter tractable when parameterized by the size of $T$. Previously only NP-completeness was known. - Contrasting to the first result, we prove that if $H$ is not constrained to be a tree, then the recognition problem becomes much harder. Namely, we show that there is a multigraph $H$ with 4 vertices and 5 edges such that it is NP-complete to decide whether $G$ is a proper $H$-graph.


title: Streaming Multi-speaker ASR with RNN-T

abstract: Recent research shows end-to-end ASR systems can recognize overlapped speech from multiple speakers. However, all published works have assumed no latency constraints during inference, which does not hold for most voice assistant interactions. This work focuses on multi-speaker speech recognition based on a recurrent neural network transducer (RNN-T) that has been shown to provide high recognition accuracy at a low latency online recognition regime. We investigate two approaches to multi-speaker model training of the RNN-T: deterministic output-target assignment and permutation invariant training. We show that guiding separation with speaker order labels in the former case enhances the high-level speaker tracking capability of RNN-T. Apart from that, with multistyle training on single- and multi-speaker utterances, the resulting models gain robustness against ambiguous numbers of speakers during inference. Our best model achieves a WER of 10.2% on simulated 2-speaker LibriSpeech data, which is competitive with the previously reported state-of-the-art nonstreaming model (10.3%), while the proposed model could be directly applied for streaming applications.


title: Characterization and Computation of Feasible Trajectories for an Articulated Probe with a Variable-Length End Segment

abstract: An articulated probe is modeled in the plane as two line segments, $ab$ and $bc$, joined at $b$, with $ab$ being very long, and $bc$ of some small length $r$. We investigate a trajectory planning problem involving the articulated two-segment probe where the length $r$ of $bc$ can be customized. Consider a set $P$ of simple polygonal obstacles with a total of $n$ vertices, a target point $t$ located in the free space such that $t$ cannot see to infinity, and a circle $S$ centered at $t$ enclosing $P$. The probe initially resides outside $S$, with $ab$ and $bc$ being collinear, and is restricted to the following sequence of moves: a straight line insertion of $abc$ into $S$ followed by a rotation of $bc$ around $b$. The goal is to compute a feasible obstacle-avoiding trajectory for the probe so that, after the sequence of moves, $c$ coincides with $t$. We prove that, for $n$ line segment obstacles, the smallest length $r$ for which there exists a feasible probe trajectory can be found in $O(n^{2+\epsilon})$ time using $O(n^{2+\epsilon})$ space, for any constant $\epsilon > 0$. Furthermore, we prove that all values $r$ for which a feasible probe trajectory exists form $O(n^2)$ intervals, and can be computed in $O(n^{5/2})$ time using $O(n^{2+\epsilon})$ space. We also show that, for a given $r$, the feasible trajectory space of the articulated probe can be characterized by a simple arrangement of complexity $O(n^2)$, which can be constructed in $O(n^2)$ time. To obtain our solutions, we design efficient data structures for a number of interesting variants of geometric intersection and emptiness query problems.


title: Does BERT Understand Sentiment? Leveraging Comparisons Between Contextual and Non-Contextual Embeddings to Improve Aspect-Based Sentiment Models

abstract: When performing Polarity Detection for different words in a sentence, we need to look at the words around to understand the sentiment. Massively pretrained language models like BERT can encode not only just the words in a document but also the context around the words along with them. This begs the questions, "Does a pretrain language model also automatically encode sentiment information about each word?" and "Can it be used to infer polarity towards different aspects?". In this work we try to answer this question by showing that training a comparison of a contextual embedding from BERT and a generic word embedding can be used to infer sentiment. We also show that if we finetune a subset of weights the model built on comparison of BERT and generic word embedding, it can get state of the art results for Polarity Detection in Aspect Based Sentiment Classification datasets.


title: Low-Resolution Face Recognition In Resource-Constrained Environments

abstract: A non-parametric low-resolution face recognition model for resource-constrained environments with limited networking and computing is proposed in this work. Such environments often demand a small model capable of being effectively trained on a small number of labeled data samples, with low training complexity, and low-resolution input images. To address these challenges, we adopt an emerging explainable machine learning methodology called successive subspace learning (SSL).SSL offers an explainable non-parametric model that flexibly trades the model size for verification performance. Its training complexity is significantly lower since its model is trained in a one-pass feedforward manner without backpropagation. Furthermore, active learning can be conveniently incorporated to reduce the labeling cost. The effectiveness of the proposed model is demonstrated by experiments on the LFW and the CMU Multi-PIE datasets.


title: Scaling Wide Residual Networks for Panoptic Segmentation

abstract: The Wide Residual Networks (Wide-ResNets), a shallow but wide model variant of the Residual Networks (ResNets) by stacking a small number of residual blocks with large channel sizes, have demonstrated outstanding performance on multiple dense prediction tasks. However, since proposed, the Wide-ResNet architecture has barely evolved over the years. In this work, we revisit its architecture design for the recent challenging panoptic segmentation task, which aims to unify semantic segmentation and instance segmentation. A baseline model is obtained by incorporating the simple and effective Squeeze-and-Excitation and Switchable Atrous Convolution to the Wide-ResNets. Its network capacity is further scaled up or down by adjusting the width (i.e., channel size) and depth (i.e., number of layers), resulting in a family of SWideRNets (short for Scaling Wide Residual Networks). We demonstrate that such a simple scaling scheme, coupled with grid search, identifies several SWideRNets that significantly advance state-of-the-art performance on panoptic segmentation datasets in both the fast model regime and strong model regime.


title: Physics of ULIRGs with MUSE and ALMA: The PUMA project I. Properties of the survey and first MUSE data results

abstract: Ultraluminous infrared galaxies (ULIRGs) are characterised by extreme starburst (SB) and AGN activity, and are therefore ideal laboratories for studying the outflow phenomena. We have recently started a project called Physics of ULIRGs with MUSE and ALMA (PUMA), which is a survey of 25 nearby (z < 0.165) ULIRGs observed with the integral field spectrograph MUSE and the interferometer ALMA. This sample includes systems with both AGN and SB nuclear activity in the pre- and post-coalescence phases of major mergers. The main goals of the project are to study the prevalence of multi-phase outflows as a function of the galaxy properties, to constrain the driving mechanisms of the outflows (e.g. distinguish between SB and AGN winds), and to identify feedback effects on the host galaxy. In this first paper, we present details on the sample selection, MUSE observations, and derive first data products. MUSE data were analysed to study the dynamical status of each of the 21 ULIRGs observed so far, taking the stellar kinematics and the morphological properties inferred from MUSE narrow-band images into account. We also located the ULIRG nuclei, using near-IR (HST) and mm (ALMA) data, and studied their optical spectra to infer the ionisation state through BPT diagnostics, and outflows in both ionised and neutral gas. We show that the morphological and stellar kinematic classifications are consistent: post-coalescence systems are more likely associated with ordered motions, while interacting (binary) systems are dominated by non-ordered and streaming motions. We also find broad and asymmetric [OIII] and NaID profiles in almost all nuclear spectra, with line widths in the range 300-2000 km/s, possibly associated with AGN- and SB-driven winds. This result reinforces previous findings that indicated that outflows are ubiquitous during the pre- and post-coalescence phases of major mergers.


title: Toward Quantum Simulations of $\mathbb{Z}_2$ Gauge Theory Without State Preparation

abstract: Preparing strongly-coupled particle states on quantum computers requires large resources. In this work, we show how classical sampling coupled with projection operators can be used to compute Minkowski matrix elements without explicitly preparing these states on the quantum computer. We demonstrate this for the 2+1d $\mathbb{Z}_2$ lattice gauge theory on small lattices with a quantum simulator.


title: Singular solutions to parabolic equations in nondivergence form

abstract: For any $\alpha \in (0,1)$, we construct an example of a solution to a parabolic equation with measurable coefficients in two space dimensions which has an isolated singularity and is not better that $C^\alpha$. We prove that there exists no solution to a fully nonlinear uniformly parabolic equation, in any dimension, which has an isolated singularity where it is not $C^2$ while it is analytic elsewhere, and it is homogeneous in $x$ at the time of the singularity. We build an example of a non homogeneous solution to a fully nonlinear uniformly parabolic equation with an isolated singularity, which we verify with the aid of a numerical computation.


title: Ensemble- and Distance-Based Feature Ranking for Unsupervised Learning

abstract: In this work, we propose two novel (groups of) methods for unsupervised feature ranking and selection. The first group includes feature ranking scores (Genie3 score, RandomForest score) that are computed from ensembles of predictive clustering trees. The second method is URelief, the unsupervised extension of the Relief family of feature ranking algorithms. Using 26 benchmark data sets and 5 baselines, we show that both the Genie3 score (computed from the ensemble of extra trees) and the URelief method outperform the existing methods and that Genie3 performs best overall, in terms of predictive power of the top-ranked features. Additionally, we analyze the influence of the hyper-parameters of the proposed methods on their performance, and show that for the Genie3 score the highest quality is achieved by the most efficient parameter configuration. Finally, we propose a way of discovering the location of the features in the ranking, which are the most relevant in reality.


title: Bridging the Planet Radius Valley: Stellar Clustering as a Key Driver for Turning Sub-Neptunes into Super-Earths

abstract: Extrasolar planets with sizes between that of the Earth and Neptune ($R_{\rm p}=1{-}4~{\rm R}_\oplus$) have a bimodal radius distribution. This 'planet radius valley' separates compact, rocky super-Earths ($R_{\rm p}=1.0{-}1.8~{\rm R}_\oplus$) from larger sub-Neptunes ($R_{\rm p}=1.8{-}3.5~{\rm R}_\oplus$) hosting a gaseous hydrogen-helium envelope around their rocky core. Various hypotheses for this radius valley have been put forward, which all rely on physics internal to the planetary system: photoevaporation by the host star, long-term mass loss driven by the cooling planetary core, or the transition between two fundamentally different planet formation modes as gas is lost from the protoplanetary disc. Here we report the discovery that the planet radius distribution exhibits a strong dependence on ambient stellar clustering, characterised by measuring the position-velocity phase space density with \textit{Gaia}. When dividing the planet sample into 'field' and 'overdensity' sub-samples, we find that planetary systems in the field exhibit a statistically significant ($p=5.5\times10^{-3}$) dearth of planets below the radius valley compared to systems in phase space overdensities. This implies that the large-scale stellar environment of a planetary system is a key factor setting the planet radius distribution. We discuss how models for the radius valley might be revised following our findings and conclude that a multi-scale, multi-physics scenario is needed, connecting planet formation and evolution, star and stellar cluster formation, and galaxy evolution.


title: Carbon stars as standard candles: II. The median J magnitude as a distance indicator

abstract: We introduce a new distance determination method using carbon-rich asymptotic giant branch stars (CS) as standard candles and the Large and Small Magellanic Clouds (LMC and SMC) as the fundamental calibrators. We select the samples of CS from the ($(J-K_{s})_0$, $J_0$) colour-magnitude diagrams, as, in this combination of filters, CS are bright and easy to identify. We fit the CS $J$-band luminosity functions using a Lorentzian distribution modified to allow the distribution to be asymmetric. We use the parameters of the best-fit distribution to determine if the CS luminosity function of a given galaxy resembles that of the LMC or SMC. Based on this resemblance, we use either the LMC or SMC as the calibrator and estimate the distance to the given galaxy using the median $J$ magnitude ($\overline{J}$) of the CS samples. We apply this new method to the two Local Group galaxies NGC 6822 and IC 1613. We find that NGC 6822 has an "LMC-like" CS luminosity function while IC 1613 is more "SMC-like". Using the values for the median absolute $J$ magnitude for the LMC and SMC found in Paper I we find a distance modulus of $\mu_{0}=23.54\pm0.03$ (stat) for NGC 6822 and $\mu_{0}=24.34\pm0.05$ (stat) for IC 1613.


title: Efficient Construction of Nonlinear Models over Normalized Data

abstract: Machine Learning (ML) applications are proliferating in the enterprise. Relational data which are prevalent in enterprise applications are typically normalized; as a result, data has to be denormalized via primary/foreign-key joins to be provided as input to ML algorithms. In this paper, we study the implementation of popular nonlinear ML models, Gaussian Mixture Models (GMM) and Neural Networks (NN), over normalized data addressing both cases of binary and multi-way joins over normalized relations. For the case of GMM, we show how it is possible to decompose computation in a systematic way both for binary joins and for multi-way joins to construct mixture models. We demonstrate that by factoring the computation, one can conduct the training of the models much faster compared to other applicable approaches, without any loss in accuracy. For the case of NN, we propose algorithms to train the network taking normalized data as the input. Similarly, we present algorithms that can conduct the training of the network in a factorized way and offer performance advantages. The redundancy introduced by denormalization can be exploited for certain types of activation functions. However, we demonstrate that attempting to explore this redundancy is helpful up to a certain point; exploring redundancy at higher layers of the network will always result in increased costs and is not recommended. We present the results of a thorough experimental evaluation, varying several parameters of the input relations involved and demonstrate that our proposals for the training of GMM and NN yield drastic performance improvements typically starting at 100%, which become increasingly higher as parameters of the underlying data vary, without any loss in accuracy.


title: Existence and uniqueness of global weak solutions to strain-limiting viscoelasticity with Dirichlet boundary data

abstract: We consider a system of evolutionary equations that is capable of describing certain viscoelastic effects in linearized yet nonlinear models of solid mechanics. The essence of the paper is that the constitutive relation, involving the Cauchy stress, the small strain tensor and the symmetric velocity gradient, is given in an implicit form. For a large class of implicit constitutive relations we establish the existence and uniqueness of a global-in-time large-data weak solution. We then focus on the class of so-called limiting strain models, i.e., models for which the magnitude of the strain tensor is known to remain small a priori, regardless of the magnitude of the Cauchy stress tensor. For this class of models, a new technical difficulty arises, which is that the Cauchy stress is only an integrable function over its domain of definition, resulting in the underlying function spaces being nonreflexive and thus the weak compactness of bounded sequences of elements of these spaces is lost. Nevertheless, even for problems of this type we are able to provide a satisfactory existence theory, as long as the initial data have finite elastic energy and the boundary data fulfill natural compatibility conditions.


title: Time evolution of gaps in stellar streams in axisymmetric St\"ackel potentials

abstract: When a subhalo interacts with a cold stellar stream it perturbs its otherwise nearly smooth distribution of stars, and this leads to the creation of a gap. The properties of such gaps depend on the parameters of the interaction. Their characterisation could thus lead to the determination of the mass spectrum of the perturbers and possibly reveal the existence of dark subhalos orbiting the Milky Way. Our goal is to construct a fully analytical model of the formation and evolution of gaps embedded in streams orbiting in a realistic Milky Way potential. To this end, we extend the model of Helmi & Koppelman (2016) for spherical potentials, and predict the properties of gaps in streams evolving in axisymmetric St\"ackel potentials. We make use of action-angles and their simple behaviour to calculate the divergence of initially nearby orbits slightly perturbed by the interaction with a subhalo. Our model, corroborated by N-body experiments, predicts that the size of a gap grows linearly with time. We obtain analytical expressions for the dependencies of the growth rate on the orbit of the stream, the properties of the subhalo (mass, scale-radius), and the geometry of the encounter (relative velocity, impact parameter). We find that the density at the centre of the gap decreases with time as a power-law in the same way as the density of a stream. This results in the density-contrast between a pristine and a perturbed stream on the same orbit asymptotically reaching a constant value dependent only on the encounter parameters. We find that at a fixed age, smallish gaps are sensitive mostly to the mass of the subhalo, while gaps formed by subhalo flybys with a low relative velocity, or when the stream and subhalo move parallel, are degenerate to the encounter parameters.


title: Computing Feasible Trajectories for an Articulated Probe in Three Dimensions

abstract: Consider an input consisting of a set of $n$ disjoint triangular obstacles in $\mathbb{R}^3$ and a target point $t$ in the free space, all enclosed by a large sphere $S$ of radius $R$ centered at $t$. An articulated probe is modeled as two line segments $ab$ and $bc$ connected at point $b$. The length of $ab$ can be equal to or greater than $R$, while $bc$ is of a given length $r \leq R$. The probe is initially located outside $S$, assuming an unarticulated configuration, in which $ab$ and $bc$ are collinear and $b \in ac$. The goal is to find a feasible (obstacle-avoiding) probe trajectory to reach $t$, with the condition that the probe is constrained by the following sequence of moves -- a straight-line insertion of the unarticulated probe into $S$, possibly followed by a rotation of $bc$ at $b$ for at most $\pi/2$ radians, so that $c$ coincides with $t$. We prove that if there exists a feasible probe trajectory, then a set of extremal feasible trajectories must be present. Through careful case analysis, we show that these extremal trajectories can be represented by $O(n^4)$ combinatorial events. We present a solution approach that enumerates and verifies these combinatorial events for feasibility in overall $O(n^{4+\epsilon})$ time using $O(n^{4+\epsilon})$ space, for any constant $\epsilon > 0$. The enumeration algorithm is highly parallel, considering that each combinatorial event can be generated and verified for feasibility independently of the others. In the process of deriving our solution, we design the first data structure for addressing a special instance of circular sector emptiness queries among polyhedral obstacles in three dimensional space, and provide a simplified data structure for the corresponding emptiness query problem in two dimensions.


title: Deconfinement Phase Transition under Chemical Equilibrium

abstract: In this work, we investigate how the assumption of chemical equilibrium with leptons affects the deconfinement phase transition to quark matter. This is done within the framework of the Chiral Mean Field model (CMF) allowing for non-zero net strangeness, corresponding to the conditions found in astrophysical scenarios. We build 3-dimensional QCD phase diagrams with temperature, baryon chemical potential, and either charge or isospin fraction or chemical potential to show how the deconfinement region collapses to a line in the special case of chemical equilibrium, such as the one established the interior of cold catalyzed neutron stars.


title: A Comprehensive Study of the Three- and Four-Neutron Systems at Low Energies

abstract: This work presents further analysis of the three- and four-neutron systems in the low energy regime using adiabatic hyperspherical methods. In our previous Phys. Rev. Lett. article (Phys. Rev. Lett. 125, 052501 (2020)), the low-energy behavior of these neutron systems was treated in the adiabatic approximation, neglecting the off-diagonal non-adiabatic couplings. A thorough analysis of the density of states through a multi-channel treatment of the three-and four-neutron scattering near the scattering continuum threshold is performed, showing no evidence of a 4n resonance at low energy. A detailed analysis of the long-range behavior of the lowest few adiabatic hyperspherical potentials shows there is an attractive $\rho^{-3}$ universal behavior which dominates in the low-energy regime of the multi-channel scattering. This long-range behavior leads to a divergent behavior of the density of state for $E\rightarrow0$ that could account for the low-energy signal observed in the 2016 experiment by Kisamori et al. (Phys. Rev. Lett. 116, 052501 (2016)).


title: An analysis of replies to Trump's tweets

abstract: Donald Trump has tweeted thousands of times during his presidency. These public statements are an increasingly important way through which Trump communicates his political and personal views. A better understanding of the way the American public consumes and responds to these tweets is therefore critical. In the present work, we address both consumption of and response to Trump's tweets by studying replies to them on Twitter. With respect to response, we find that a small number of older, white, left-leaning, and female Americans are responsible for the vast majority of replies to Trump's tweets. These individuals also attend to a broader range of Trump's tweets than the rest of the individuals we study. With respect to consumption, we note that Trump's tweets are often viewed not in isolation, but rather in the context of a set of algorithmically-curated replies. These replies may therefore color the way Americans consume Trump's tweets. To this end, we find some evidence that Twitter accounts see replies in line with their political leanings. However, we show that this can be largely, although not entirely, attributed to the fact that Twitter is more likely to show replies by accounts a user follows. As a basis for comparison, all results for Trump are compared and contrasted with replies to Joe Biden's tweets.