title: Quantum system dynamics with a weakly nonlinear Josephson junction bath

abstract: We investigate the influence of a weakly nonlinear Josephson bath consisting of a chain of Josephson junctions on the dynamics of a small quantum system (LC oscillator). Focusing on the regime where the charging energy is the largest energy scale, we perturbatively calculate the correlation function of the Josephson bath to the leading order in the Josephson energy divided by the charging energy while keeping the cosine potential exactly. When the variation of the charging energy along the chain ensures fast decay of the bath correlation function, the dynamics of the LC oscillator that is weakly and capacitively coupled to the Josephson bath can be solved through the Markovian master equation. We establish a duality relation for the Josephson bath between the regimes of large charging and Josephson energies respectively. The results can be applied to cases where the charging energy either is nonuniformly engineered or disordered in the chain. Furthermore, we find that the Josephson bath may become non-Markovian when the temperature is increased beyond the zero-temperature limit in that the bath correlation function gets shifted by a constant and does not decay with time.


title: Dimensional analysis on superconducting plasmonics

abstract: Plasmonic response of superconductors at various dimensions are addressed in this paper. All possible parameter space is discussed and considered for theoretical demonstration towards possible future experiments. The most critical parameters in each dimensional state are found to be carrier concentration (3D), wavevector (2D), surrounding permittivity (1D) and particle radius (0D). However, it is also demonstrated that unless all parameters combined with optimal conditions, the superconducting plasmonic system does not result in an ideal lossless environment. Furthermore, under all conditions, lower temperature decreases the loss content affiliated with electron and phonon scattering.


title: Multi-Scale Merge-Split Markov Chain Monte Carlo for Redistricting

abstract: We develop a Multi-Scale Merge-Split Markov chain on redistricting plans. The chain is designed to be usable as the proposal in a Markov Chain Monte Carlo (MCMC) algorithm. Sampling the space of plans amounts to dividing a graph into a partition with a specified number of elements which each correspond to a different district. The districts satisfy a collection of hard constraints and the measure may be weighted with regard to a number of other criteria. The multi-scale algorithm is similar to our previously developed Merge-Split proposal, however, this algorithm provides improved scaling properties and may also be used to preserve nested communities of interest such as counties and precincts. Both works use a proposal which extends the ReCom algorithm which leveraged spanning trees merge and split districts. In this work we extend the state space so that each district is defined by a hierarchy of trees. In this sense, the proposal step in both algorithms can be seen as a "Forest ReCom." We also expand the state space to include edges that link specified districts, which further improves the computational efficiency of our algorithm. The collection of plans sampled by the MCMC algorithm can serve as a baseline against which a particular plan of interest is compared. If a given plan has different racial or partisan qualities than what is typical of the collection of plans, the given plan may have been gerrymandered and is labeled as an outlier.


title: Communicative Reinforcement Learning Agents for Landmark Detection in Brain Images

abstract: Accurate detection of anatomical landmarks is an essential step in several medical imaging tasks. We propose a novel communicative multi-agent reinforcement learning (C-MARL) system to automatically detect landmarks in 3D brain images. C-MARL enables the agents to learn explicit communication channels, as well as implicit communication signals by sharing certain weights of the architecture among all the agents. The proposed approach is evaluated on two brain imaging datasets from adult magnetic resonance imaging (MRI) and fetal ultrasound scans. Our experiments show that involving multiple cooperating agents by learning their communication with each other outperforms previous approaches using single agents.


title: Panchromatic calibration of Ca II triplet luminosity dependence

abstract: (ABRIDGED) Context. The line strength of the Ca II triplet (CaT) lines are a proxy to measure metallicity from individual stellar spectra of bright red giant stars. It is a mandatory step to remove the magnitude (proxy for gravity, temperature and luminosity) dependence from the equivalent width (EW) of the lines before converting them into metallicities. The working empirical procedure used for decades is to use the relative magnitude with respect to the horizontal branch level. Aims. The V filter is broadly adopted as the reference magnitude, although a few works have used different filters (I and Ks, for example). In this work we investigate the dependence of the CaT calibration using griz filters from the DECam and the GMOS, G from Gaia, BVI filters from the MCPS, YJKs filters from VIRCAM. We use as a reference FORS2 V filter used in the original analysis of the sample. Methods. Red giant stars from clusters with known metallicity and available CaT equivalent widths are used as reference. Public photometric catalogues are taken from SMASH DR2, VMC, Gaia, MCPS surveys plus VISCACHA-GMOS data, for a selection of Small Magellanic Cloud clusters. The slopes are fitted using two and three lines to be applicable to most of the metallicity scales. Results. The magnitude dependence of the CaT EWs is well described by a linear relation using any filter analysed in this work. The slope increases with wavelength of the filters. The zero point (a.k.a. reduced equivalent width), that is the metallicity indicator, remains the same. Conclusions. If the same line profile function is used with the same bandpasses and continuum regions, and the total EW comes from the same number of lines (2 or 3), then the reduced EW is the same regardless the filter used. Therefore, any filter can be used to convert the CaT equivalent widths into metallicity for a given CaT calibration.


title: Benchmarking network fabrics for data distributed training of deep neural networks

abstract: Artificial Intelligence/Machine Learning applications require the training of complex models on large amounts of labelled data. The large computational requirements for training deep models have necessitated the development of new methods for faster training. One such approach is the data parallel approach, where the training data is distributed across multiple compute nodes. This approach is simple to implement and supported by most of the commonly used machine learning frameworks. The data parallel approach leverages MPI for communicating gradients across all nodes. In this paper, we examine the effects of using different physical hardware interconnects and network-related software primitives for enabling data distributed deep learning. We compare the effect of using GPUDirect and NCCL on Ethernet and OmniPath fabrics. Our results show that using Ethernet-based networking in shared HPC systems does not have a significant effect on the training times for commonly used deep neural network architectures or traditional HPC applications such as Computational Fluid Dynamics.


title: Probing the Deuteron at Very Large Internal Momenta

abstract: $^{2}\mathrm{H}(e,e'p)n$ cross sections have been measured at 4-momentum transfers of $Q^{2} = 4.5 \pm 0.5$ (GeV/c)$^{2}$ over a range of neutron recoil momenta, $p_{\mathrm{r}}$, reaching up to $\sim1.0$ GeV/c. The data were obtained at fixed neutron recoil angles $\theta_{nq} = 35^\circ$, $45^\circ$ and $75^{\circ}$ with respect to the 3-momentum transfer $\vec q$. The new data agree well with previous data which reached $p_{\mathrm{r}}\sim500$ MeV/c. At $\theta_{nq} = 35^\circ$ and $45^\circ$, final state interactions (FSI), meson exchange currents (MEC) and isobar currents (IC) are suppressed and the plane wave impulse approximation (PWIA) provides the dominant cross section contribution. The new data are compared to recent theoretical calculations, where we observe a significant discrepancy for missing momenta $p_{\mathrm{r}}>700$ MeV/c.


title: When Hardness of Approximation Meets Hardness of Learning

abstract: A supervised learning algorithm has access to a distribution of labeled examples, and needs to return a function (hypothesis) that correctly labels the examples. The hypothesis of the learner is taken from some fixed class of functions (e.g., linear classifiers, neural networks etc.). A failure of the learning algorithm can occur due to two possible reasons: wrong choice of hypothesis class (hardness of approximation), or failure to find the best function within the hypothesis class (hardness of learning). Although both approximation and learnability are important for the success of the algorithm, they are typically studied separately. In this work, we show a single hardness property that implies both hardness of approximation using linear classes and shallow networks, and hardness of learning using correlation queries and gradient-descent. This allows us to obtain new results on hardness of approximation and learnability of parity functions, DNF formulas and $AC^0$ circuits.


title: Personalized Deep Learning for Ventricular Arrhythmias Detection on Medical IoT Systems

abstract: Life-threatening ventricular arrhythmias (VA) are the leading cause of sudden cardiac death (SCD), which is the most significant cause of natural death in the US. The implantable cardioverter defibrillator (ICD) is a small device implanted to patients under high risk of SCD as a preventive treatment. The ICD continuously monitors the intracardiac rhythm and delivers shock when detecting the life-threatening VA. Traditional methods detect VA by setting criteria on the detected rhythm. However, those methods suffer from a high inappropriate shock rate and require a regular follow-up to optimize criteria parameters for each ICD recipient. To ameliorate the challenges, we propose the personalized computing framework for deep learning based VA detection on medical IoT systems. The system consists of intracardiac and surface rhythm monitors, and the cloud platform for data uploading, diagnosis, and CNN model personalization. We equip the system with real-time inference on both intracardiac and surface rhythm monitors. To improve the detection accuracy, we enable the monitors to detect VA collaboratively by proposing the cooperative inference. We also introduce the CNN personalization for each patient based on the computing framework to tackle the unlabeled and limited rhythm data problem. When compared with the traditional detection algorithm, the proposed method achieves comparable accuracy on VA rhythm detection and 6.6% reduction in inappropriate shock rate, while the average inference latency is kept at 71ms.


title: The Near-Infrared Spectrum of the Nuclear Star Cluster: Looking Below the Tip of the Iceberg and Comparisons with Extragalactic Nuclei

abstract: Long-slit near-infrared (NIR) spectra of the Galactic nuclear star cluster (NSC) are discussed. The spectra sample the major axis of the NSC out to its half light radius. The equivalent widths of NaI2.21um and CaI2.26um outside of the center of the NSC, where light from nuclear bulge stars contributes a large fraction to the total flux, are consistent with solar chemical mixtures. In contrast, the equivalent widths of NaI2.21um and CaI2.26um near the center of the NSC are larger than expected from models with solar chemical mixtures, even after light from the brightest evolved stars is removed. That NaI2.21um and CaI2.26um are deeper than predicted for solar chemical mixtures over a range of evolutionary states is consistent with previous studies that find that the majority of stars near the center of the NSC formed from material that had non-solar chemical mixtures. A spectrum of nebular emission and the youngest stars near the GC is also extracted. The equivalent widths of emission features in the extracted NIR spectrum are similar to those in the nuclear spectrum of NGC 253, and it is argued that this agreement is best achieved if the current episode of star formation near the center of the NSC has been in progress for at least a few Myr.


title: Compute, Time and Energy Characterization of Encoder-Decoder Networks with Automatic Mixed Precision Training

abstract: Deep neural networks have shown great success in many diverse fields. The training of these networks can take significant amounts of time, compute and energy. As datasets get larger and models become more complex, the exploration of model architectures becomes prohibitive. In this paper we examine the compute, energy and time costs of training a UNet based deep neural network for the problem of predicting short term weather forecasts (called precipitation Nowcasting). By leveraging a combination of data distributed and mixed-precision training, we explore the design space for this problem. We also show that larger models with better performance come at a potentially incremental cost if appropriate optimizations are used. We show that it is possible to achieve a significant improvement in training time by leveraging mixed-precision training without sacrificing model performance. Additionally, we find that a 1549% increase in the number of trainable parameters for a network comes at a relatively smaller 63.22% increase in energy usage for a UNet with 4 encoding layers.


title: AB3DMOT: A Baseline for 3D Multi-Object Tracking and New Evaluation Metrics

abstract: 3D multi-object tracking (MOT) is essential to applications such as autonomous driving. Recent work focuses on developing accurate systems giving less attention to computational cost and system complexity. In contrast, this work proposes a simple real-time 3D MOT system with strong performance. Our system first obtains 3D detections from a LiDAR point cloud. Then, a straightforward combination of a 3D Kalman filter and the Hungarian algorithm is used for state estimation and data association. Additionally, 3D MOT datasets such as KITTI evaluate MOT methods in 2D space and standardized 3D MOT evaluation tools are missing for a fair comparison of 3D MOT methods. We propose a new 3D MOT evaluation tool along with three new metrics to comprehensively evaluate 3D MOT methods. We show that, our proposed method achieves strong 3D MOT performance on KITTI and runs at a rate of $207.4$ FPS on the KITTI dataset, achieving the fastest speed among modern 3D MOT systems. Our code is publicly available at http://www.xinshuoweng.com/projects/AB3DMOT.


title: Embedded Fracture Model for Coupled Flow and Geomechanics

abstract: Fluid injection and production cause changes in reservoir pressure, which result in deformations in the subsurface. This phenomenon is particularly important in reservoirs with abundant fractures and faults because the induced slip and opening of the fractures may significantly alter their hydraulic properties. Modeling strongly coupled poro-mechanical processes in naturally fractured reservoirs is a challenging problem. The Discrete Fracture Model (DFM) is a state-of-art method for modeling coupled flow and mechanics in fractured reservoirs. This method requires constructing computational grids that comform to fractures, which is very challenging in complex 3D settings. The objective of this study is to develop a numerical method that does not require gridding near fractures and can efficiently model hydromechanical interactions in fractured reservoirs. We utilize formulations based on the Strong Discontinuity Approach (SDA) for mechanics and Embedded Discrete Fracture Model (EDFM) for flow. We first present a mathematical formulation and emphasize the kinematic aspects of fracture slip and opening. We then introduce a series of mechanical tests that investigate the spatial convergence of the model and compare its accuracy with the Discrete Fracture Model (DFM). We finally consider a synthetic coupled case of a reservoir with several fractures and compare the performance of the SDA and DFM methods. Our results indicate super-linear spatial convergence of the proposed SDA algorithm. Numerical simulations confirm the applicability of the proposed method to modeling the coupling effects in subsurface applications.


title: Pseudo-differential operators associated to general type I locally compact groups

abstract: In a recent paper by M. Mantoiu and M. Ruzhansky, a global pseudo-differential calculus has been developed for unimodular groups of type I. In the present article we generalize the main results to arbitrary locally compact groups of type I. Our methods involve the use of Plancherel's theorem for non-unimodular groups. We also make connections with a $C^*$-algebraic formalism, involving dynamical systems, and give explicit constructions for the group of affine transformations of the real line.


title: Generalized Spontaneous Symmetry Breaking

abstract: We propose a new type of symmetry breaking mechanism that takes boundaries into account, and show how it can detect surface modes by interpreting them as the order parameter associated with a generalized symmetry breaking. We argue that this mechanism is analogous to, and yet distinct from, the standard Ginzburg-Landau scenarios, and show how the two communicate through the specification of boundary conditions.


title: Examples of equivariant Lagrangian mean curvature flow

abstract: In this expository note we describe important examples of Lagrangian mean curvature flow in $\mathbb{C}^2$ which are invariant under a circle action. Through these examples, we see compact and non-compact situations, long-time existence, singularities forming via explicit models, and significant objects in Riemannian and symplectic geometry, including the Clifford torus, Chekanov torus, Whitney sphere and Lawlor necks.


title: Energy-Optimal Control of a Submarine-Launched Cruise Missile

abstract: A typical mission profile of submarine-launched cruise missiles begins with the launch phase which covers the motion of the missile from the launch to the water-exit and continues with the boost phase which lasts from the water-exit to the beginning of the cruise phase. In order to achieve the desired range of the launch and boost phases, efficient utilization of available energy which carries the missile to the beginning of the cruise phase is necessary. For this purpose, this study presents a new approach for energy-optimal control of the underwater and air motion of a submarine-launched cruise missile. In this approach, the aforementioned problem is modeled and solved as a minimum-effort optimal control problem. Then, the effects of initial and final conditions on energy need are investigated, and the optimal conditions that result with the minimum energy need are determined. Prior to the guidance and control design steps, six degrees of freedom (6 DOF) motion equations are derived and the hydrodynamic and aerodynamic parameters are retrieved. The nonlinear 6 DOF motion model is simplified and linearized before minimum-effort optimal control design part. Results of the designed guidance and control strategies are presented through the nonlinear 6 DOF simulations. Finally, some comments are made and future studies are mentioned based on theoretical and simulation studies.


title: The science enabled by a dedicated solar system space telescope

abstract: The National Academy Committee on Astrobiology and Planetary Science (CAPS) made a recommendation to study a large/medium-class dedicated space telescope for planetary science, going beyond the Discovery-class dedicated planetary space telescope endorsed in Visions and Voyages. Such a telescope would observe targets across the entire solar system, engaging a broad spectrum of the science community. It would ensure that the high-resolution, high-sensitivity observations of the solar system in visible and UV wavelengths revolutionized by the Hubble Space Telescope (HST) could be extended. A dedicated telescope for solar system science would: (a) transform our understanding of time-dependent phenomena in our solar system that cannot be studied currently under programs to observe and visit new targets and (b) enable a comprehensive survey and spectral characterization of minor bodies across the solar system, which requires a large time allocation not supported by existing facilities. The time-domain phenomena to be explored are critically reliant on high spatial resolution UV-visible observations. This paper presents science themes and key questions that require a long-lasting space telescope dedicated to planetary science that can capture high-quality, consistent data at the required cadences that are free from effects of the terrestrial atmosphere and differences across observing facilities. Such a telescope would have excellent synergy with astrophysical facilities by placing planetary discoveries made by astrophysics assets in temporal context, as well as triggering detailed follow-up observations using larger telescopes. The telescope would support future missions to the Ice Giants, Ocean Worlds, and minor bodies across the solar system by placing the results of such targeted missions in the context of longer records of temporal activities and larger sample populations.


title: Cohen-Macaulayness of absolute integral closures

abstract: We prove that, modulo any power of a prime $p$, the absolute integral closure of an excellent noetherian domain is Cohen-Macaulay. A graded analog is also established, yielding variants of Kodaira vanishing "up to finite covers" in mixed characteristic. Our main tools are (log) prismatic cohomology (which yields a Frobenius action in mixed characteristic) and the $p$-adic Riemann-Hilbert functor for constructible \'etale $\mathbf{F}_p$-sheaves on varieties over a $p$-adic field (which almost controls perfectified prismatic cohomology).


title: Robust Mean Estimation on Highly Incomplete Data with Arbitrary Outliers

abstract: We study the problem of robustly estimating the mean of a $d$-dimensional distribution given $N$ examples, where most coordinates of every example may be missing and $\varepsilon N$ examples may be arbitrarily corrupted. Assuming each coordinate appears in a constant factor more than $\varepsilon N$ examples, we show algorithms that estimate the mean of the distribution with information-theoretically optimal dimension-independent error guarantees in nearly-linear time $\widetilde O(Nd)$. Our results extend recent work on computationally-efficient robust estimation to a more widely applicable incomplete-data setting.


title: AssembleNet++: Assembling Modality Representations via Attention Connections

abstract: We create a family of powerful video models which are able to: (i) learn interactions between semantic object information and raw appearance and motion features, and (ii) deploy attention in order to better learn the importance of features at each convolutional block of the network. A new network component named peer-attention is introduced, which dynamically learns the attention weights using another block or input modality. Even without pre-training, our models outperform the previous work on standard public activity recognition datasets with continuous videos, establishing new state-of-the-art. We also confirm that our findings of having neural connections from the object modality and the use of peer-attention is generally applicable for different existing architectures, improving their performances. We name our model explicitly as AssembleNet++. The code will be available at: https://sites.google.com/corp/view/assemblenet/


title: On the Evolution of Subjective Experience

abstract: Subjective Experience (SE) is part of the ancient mind-body problem, which continues to be one of deepest mysteries of science. Despite major advances in many fields, there is still no plausible causal link between SE and its realization in the body. The core issue is the incompatibility of objective (3rd person) public science with subjective (1st person) private experience. Any scientific approach to SE assumes that it arose from extended evolutionary processes and that examining evolutionary history should help us understand it. While the core mystery remains, converging evidence from theoretical, experimental, and computational studies yields strong constraints on SE and some suggestions for further research. All animals confront many of the same fitness challenges. They all need some kind of internal model to relate their life goals and actionable sensed information to action. We understand the evolution of the bodily aspects of human perception and emotion, but not the SE. The first evolutionary evidence for SE appears in vertebrates and much of its neural substrate and simulation mechanism is preserved in mammals and humans. People exhibit the same phenomena, but there are remaining mysteries of everyday experience that are demonstrably incompatible with current neuroscience. In spite of this limitation, there is considerable progress on understanding the role of SE in the success of prostheses.


title: Identification of the K-most Vulnerable Entities in a Smart Grid System

abstract: A smart grid system can be considered as a multi-layered network with power network in one layer and communication network in the other. The entities in both the layers exhibit complex intra-and-interdependencies between them. A reliable decision making by the smart grid operator is contingent upon correct analysis of such dependencies between its entities and also on accurate identification of the most critical entities in the system. The Modified Implicative Interdependency Model (MIIM) [1] successfully captures such dependencies using multi-valued Boolean Logic based equations called Interdependency Relations (IDRs) after most of the existing models made failed attempts in doing that. In this paper, for any given integer K, this model is used to identify the K-most vulnerable entities in a smart grid, failure of which can maximize the network damage. Owing to the problem being NP complete, an Integer Linear Programming (ILP) based solution is given here. Validation of the model [1] and the results of the ILP based solution is done by simulating a smart grid system of IEEE 14-Bus using MATPOWER and Java Network Simulator (JNS). Simulation results prove that not only the model MIIM [1] is correct but also it can predict the network damage for failure of K-most vulnerable entities more accurately than its predecessor Implicative Interdependency Model (IIM) [2].


title: Correspondence between dissipative phase transitions of light and time crystals

abstract: We predict the emergence of a time crystal generated by an incoherently driven and dissipative nonlinear optical oscillator, where the nonlinearity also comes from dissipation. We show that a second-order dissipative phase transition of light occurs in the frame rotating at the cavity frequency, while a boundary (dissipative) time crystal emerges in the laboratory frame. We relate these two phenomena by using the Liouvillian superoperator associated with the Lindblad master equation and its symmetries. These results connect the emergence of a second-order dissipative phase transition and a dissipative time crystal in the thermodynamic limit, allowing to interpret them as the same phenomenon in terms of the Liouvillian spectrum, but just in different frames.


title: Deploying Lifelong Open-Domain Dialogue Learning

abstract: Much of NLP research has focused on crowdsourced static datasets and the supervised learning paradigm of training once and then evaluating test performance. As argued in de Vries et al. (2020), crowdsourced data has the issues of lack of naturalness and relevance to real-world use cases, while the static dataset paradigm does not allow for a model to learn from its experiences of using language (Silver et al., 2013). In contrast, one might hope for machine learning systems that become more useful as they interact with people. In this work, we build and deploy a role-playing game, whereby human players converse with learning agents situated in an open-domain fantasy world. We show that by training models on the conversations they have with humans in the game the models progressively improve, as measured by automatic metrics and online engagement scores. This learning is shown to be more efficient than crowdsourced data when applied to conversations with real users, as well as being far cheaper to collect.


title: Constraints on Stupendously Large Black Holes

abstract: We consider the observational constraints on stupendously large black holes (SLABs) in the mass range $M \gtrsim 10^{11}\,M_{\odot}$. These have attracted little attention hitherto and we are aware of no published constraints on a SLAB population in the range $(10^{12}$ - $10^{18})\,M_{\odot}$. However, there is already evidence for black holes of up to nearly $10^{11}\,M_{\odot}$ in galactic nuclei, so it is conceivable that SLABs exist and they may even have been seeded by primordial black holes. We focus on limits associated with (i) dynamical and lensing effects, (ii) the generation of background radiation through the accretion of gas during the pregalactic epoch, and (iii) the gamma-ray emission from the annihilation of the halo of weakly interacting massive particles (WIMPs) expected to form around each SLAB if these provide the dark matter. Finally, we comment on the constraints on the mass of ultra-light bosons from future measurements of the mass and spin of SLABs.


title: Closed-Loop Design of Proton Donors for Lithium-Mediated Ammonia Synthesis with Interpretable Models and Molecular Machine Learning

abstract: In this work, we experimentally determined the efficacy of several classes of proton donors for lithium-mediated electrochemical nitrogen reduction in a tetrahydrofuran-based electrolyte, an attractive alternative method for producing ammonia. We then built an interpretable data-driven classification model which identified solvatochromic Kamlet-Taft parameters as important for distinguishing between active and inactive proton donors. After curating a dataset for the Kamlet-Taft parameters, we trained a deep learning model to predict the Kamlet-Taft parameters. The combination of classification model and deep learning model provides a predictive mapping from a given proton donor to the ability to produce ammonia. We demonstrate that this combination of classification model with deep learning is superior to a purely mechanistic or data-driven approach in accuracy and experimental data efficiency.


title: Harmonic analysis of little $q$-Legendre polynomials

abstract: Many classes of orthogonal polynomials satisfy a specific linearization property giving rise to a polynomial hypergroup structure, which offers an elegant and fruitful link to Fourier analysis, harmonic analysis and functional analysis. From the opposite point of view, this allows regarding certain Banach algebras as $L^1$-algebras, associated with underlying orthogonal polynomials. The individual behavior strongly depends on these underlying polynomials. We study the little $q$-Legendre polynomials, which are orthogonal with respect to a discrete measure. We will show that their $L^1$-algebras have the property that every element can be approximated by linear combinations of idempotents. This particularly implies that these $L^1$-algebras are weakly amenable (i. e., every bounded derivation into the dual module is an inner derivation), which is known to be shared by any $L^1$-algebra of a locally compact group; in the polynomial hypergroup context, weak amenability is rarely satisfied and of particular interest because it corresponds to a certain property of the derivatives of the underlying polynomials and their (Fourier) expansions w.r.t. the polynomial basis. To our knowledge, the little $q$-Legendre polynomials yield the first example of a polynomial hypergroup whose $L^1$-algebra is weakly amenable and right character amenable but fails to be amenable. As a crucial tool, we establish certain uniform boundedness properties of the characters. Our strategy relies on the Fourier transformation on hypergroups, the Plancherel isomorphism, continued fractions, character estimations and asymptotic behavior.


title: mlr3proba: An R Package for Machine Learning in Survival Analysis

abstract: As machine learning has become increasingly popular over the last few decades, so too has the number of machine learning interfaces for implementing these models. Whilst many R libraries exist for machine learning, very few offer extended support for survival analysis. This is problematic considering its importance in fields like medicine, bioinformatics, economics, engineering, and more. mlr3proba provides a comprehensive machine learning interface for survival analysis and connects with mlr3's general model tuning and benchmarking facilities to provide a systematic infrastructure for survival modeling and evaluation.


title: Preparation of Isotopically enriched $^{112,116,120,124}$Sn targets at VECC

abstract: Resistive heating and mechanical rolling methods have been employed to prepare isotopically enriched thin target foils of 116Sn (~380 \mu g/cm2), 124Sn(~400 \mu g/cm2) and thicker foils of 112Sn (1.7 mg/cm2),120Sn (1.6 mg/cm2),respectively. Preparation of enriched targets with small amount of material, selection of releasing agent for thin targets and separation of deposited material insolvent were among the several challenges while fabrication of the thin targets. Uniformity of the targets has been measured using 241Am {\alpha}-source. NaCl has been used as releasing agent in preparation of the thin targets. These targets have been successfully used in nuclear physics experiments at VECC.


title: Fractional powers of Bessel operator and its numerical calculation

abstract: The article discusses the fractional powers of the Bessel operator and their numerical implementation. An extensive literature is devoted to the study of fractional powers of the Laplace operator and their applications. Such degrees are used in the construction of functional spaces, in the natural generalization of the Schr\"odinger equation in quantum theory, in the construction of the models of acoustic wave propagation in complex media (for example, biological tissues) and space-time models of anomalous (very slow or very fast) diffusion, in spectral theory etc. If we assume the radiality of the function on which the Laplace operator acts, then we receive the problem of constructing the fractional power of the Bessel operator. We propose to use a compositional method for constructing the operators mentioned earlier, which leads to constructions similar in their properties to the Riesz derivatives. The Hankel transform is considered as a basic integral transformation. On its basis, the compositional method proposed by V.V. Katrakhov and S.M. Sitnik, negative powers of the Bessel operator are constructed. The resulting operator contains the Gaussian hypergeometric function in the kernel. For further study, the generalized translation operator is considered in the article, and its properties are proved. For constructing a positive fractional power of the Bessel operator known methods of regularization of the integral are considered. Then, a scheme for the numerical calculation of fractional powers of the Bessel operator is proposed. This scheme is based on the Taylor--Delsarte formula obtained by B.M. Levitan.


title: Returning radiation in strong gravity around black holes: Reverberation from the accretion disc

abstract: We study reflected X-ray emission that returns to the accretion disc in the strong gravitational fields around black holes using General Relativistic ray tracing and radiative transfer calculations. Reflected X-rays that are produced when the inner regions of the disc are illuminated by the corona are subject to strong gravitational light bending, causing up to 47 per cent of the reflected emission to be returned to the disc around a rapidly spinning black hole, depending upon the scale height of the corona. The iron K line is enhanced relative to the continuum by 25 per cent, and the Compton hump by up to a factor of three. Additional light travel time between primary and secondary reflections increases the reverberation time lag measured in the iron K band by 49 per cent, while the soft X-ray lag is increased by 25 per cent and the Compton hump response time is increased by 60 per cent. Measured samples of X-ray reverberation lags are shown to be consistent with X-rays returning to the accretion disc in strong gravity. Understanding the effects of returning radiation is important in interpreting reverberation observations to probe black holes. Reflected X-rays returning to the disc can be uniquely identified by blueshifted returning iron K line photons that are Compton scattered from the inner disc, producing excess, delayed emission in the 3.5-4.5keV energy range that will be detectable with forthcoming X-ray observatories, representing a unique test of General Relativity in the strong field limit.


title: Determining the Hubble Constant without the Sound Horizon: Measurements from Galaxy Surveys

abstract: Two sources of geometric information are encoded in the galaxy power spectrum: the sound horizon at recombination and the horizon at matter-radiation equality. Analyzing the BOSS DR12 galaxy power spectra using perturbation theory with $\Omega_m$ priors from Pantheon supernovae but no priors on $\Omega_b$, we obtain constraints on $H_0$ from the second scale, finding $H_0 = 65.1^{+3.0}_{-5.4}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$; this differs from the best-fit of SH0ES at 95\% confidence. Similar results are obtained if $\Omega_m$ is constrained from uncalibrated BAO: $H_0 = 65.6^{+3.4}_{-5.5}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$. Adding the analogous lensing results from Baxter \& Sherwin 2020, the posterior shifts to $70.6^{+3.7}_{-5.0}\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$. Using mock data, Fisher analyses, and scale-cuts, we demonstrate that our constraints do not receive significant information from the sound horizon scale. Since many models resolve the $H_0$ controversy by adding new physics to alter the sound horizon, our measurements are a consistency test for standard cosmology before recombination. A simple forecast indicates that such constraints could reach $\sigma_{H_0} \simeq 1.6\,\mathrm{km}\,\mathrm{s}^{-1}\mathrm{Mpc}^{-1}$ in the era of Euclid.

categories: astro-ph.CO astro-ph.GA astro-ph.IM gr-qc hep-ph


title: Long-Lived Eccentric Modes in Circumbinary Disks

abstract: Hydrodynamical simulations show that circumbinary disks become eccentric, even when the binary is circular. Here we demonstrate that, in steady state, the disk's eccentricity behaves as a long-lived free mode trapped by turning points that naturally arise from a continuously truncated density profile. Consequently, both the disk's precession rate and eccentricity profile may be calculated via the simple linear theory for perturbed pressure-supported disks. By formulating and solving the linear theory we find that (i) surprisingly, the precession rate is roughly determined by the binary's quadrupole, even when the quadrupole is very weak relative to pressure; (ii) the eccentricity profile is largest near the inner edge of the disk, and falls exponentially outwards; and (iii) the results from linear theory indeed agree with what is found in simulations. Understanding the development of eccentric modes in circumbinary disks is a crucial first step for understanding the long term (secular) exchange of eccentricity, angular momentum and mass between the binary and the gas. Potential applications include the search for a characteristic kinematic signature in disks around candidate binaries and precession-induced modulation of accretion over long timescales.


title: Optimal detuning for quantum filter cavities

abstract: Vacuum quantum fluctuations impose a fundamental limit on the sensitivity of gravitational-wave interferometers, which rank among the most sensitive precision measurement devices ever built. The injection of conventional squeezed vacuum reduces quantum noise in one quadrature at the expense of increasing noise in the other. While this approach improved the sensitivity of the Advanced LIGO and Advanced Virgo interferometers during their third observing run (O3), future improvements in arm power and squeezing levels will bring radiation pressure noise to the forefront. Installation of a filter cavity for frequency-dependent squeezing provides broadband reduction of quantum noise through the mitigation of this radiation pressure noise, and it is the baseline approach planned for all of the future gravitational-wave detectors currently conceived. The design and operation of a filter cavity requires careful consideration of interferometer optomechanics as well as squeezing degradation processes. In this paper, we perform an in-depth analysis to determine the optimal operating point of a filter cavity. We use our model alongside numerical tools to study the implications for filter cavities to be installed in the upcoming "A+" upgrade of the Advanced LIGO detectors.


title: An Intensity Mapping Detection of Aggregate CO Line Emission at 3 mm

abstract: We present a detection of molecular gas emission at $z\sim1-5$ using the technique of line intensity mapping. We make use of a pair of 3 mm interferometric data sets, the first from the ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS), and the second from a series of Atacama Compact Array (ACA) observations conducted between 2016 and 2018, targeting the COSMOS field. At 100 GHz, we measure non-zero power at 97.8% and 99.9% confidence in the ACA and ALMA data sets, respectively. In the joint result, we reject the zero-power hypothesis at 99.99% confidence, finding $\tilde{I}^{2}_{s}(\nu)=770\pm210\ \mu\textrm{K}^2\ \textrm{Hz}\ \textrm{sr}$. After accounting for sample variance effects, the estimated spectral shot power is $\tilde{I}^{2}_{s}(\nu)=1010_{-390}^{+550}\ \mu\textrm{K}^2\ \textrm{Hz}\ \textrm{sr}$. We derive a model for the various line species our measurement is expected to be sensitive to, and estimate the shot power to be $120_{-40}^{+80}\ \mu\textrm{K}^2\ h^{-3}\,\textrm{Mpc}^{3}$, $200^{+120}_{-70}\ \mu\textrm{K}^2\ h^{-3}\,\textrm{Mpc}^{3}$, and $90^{+70}_{-40}\ \mu\textrm{K}^2\ h^{-3}\,\textrm{Mpc}^{3}$ for CO(2-1) at $z=1.3$, CO(3-2) at $z=2.5$, and CO(4-3) at $z=3.6$, respectively. Using line ratios appropriate for high-redshift galaxies, we find these results to be in good agreement with those from the CO Power Spectrum Survey (COPSS). Adopting $\alpha_{\rm CO}=3.6\ M_{\odot}\ (\textrm{K}\ \textrm{km}\ \textrm{s}^{-1}\ \textrm{pc}^{2})^{-1}$, we estimate a cosmic molecular gas density of $\rho_{\textrm{H}_2}(z)\sim 10^{8}\ M_{\odot}\ \textrm{Mpc}^{-3}$ between $z=1-3$.


title: The entropy production of an active particle in a box

abstract: A run-and-tumble particle in a one dimensional box (infinite potential well) is studied. The steady state is analytically solved and analyzed, revealing the emergent length scale of the boundary layer where particles accumulate near the walls. The mesoscopic steady state entropy production rate of the system is derived from coupled Fokker-Planck equations with a linear reaction term, resulting in an exact analytic expression. The entropy production density is shown to peak at the walls. Additionally, the derivative of the entropy production rate peaks at a system size proportional to the length scale of the accumulation boundary layer, suggesting that the behavior of the entropy production rate and its derivatives as a function of the control parameter may signify a qualitative behavior change in the physics of active systems, such as phase transitions.


title: Monte-Carlo neutrino transport in neutron star merger simulations

abstract: Gravitational waves and electromagnetic signals from merging neutron star binaries provide valuable information about the the properties of dense matter, the formation of heavy elements, and high-energy astrophysics. To fully leverage observations of these systems, we need numerical simulations that provide reliable predictions for the properties of the matter unbound in these mergers. An important limitation of current simulations is the use of approximate methods for neutrino transport that do not converge to a solution of the transport equations as numerical resolution increases, and thus have errors that are impossible to quantify. Here, we report on a first simulation of a binary neutron star merger that uses Monte-Carlo techniques to directly solve the transport equations in low-density regions. In high-density regions, we use approximations inspired by implicit Monte-Carlo to greatly reduce the cost of simulations, while only introducing errors quantifiable through more expensive convergence studies. We simulate an unequal mass neutron star binary merger up to $5\,{\rm ms}$ past merger, and report on the properties of the matter and neutrino outflows. Finally, we compare our results to the output of our best approximate `M1' transport scheme, demonstrating that an M1 scheme that carefully approximates the neutrino energy spectrum only leads to $\sim 10\%$ uncertainty in the composition and velocity of the ejecta, and $\sim20\%$ uncertainty in the $\nu_e$ and $\bar\nu_e$ luminosities and energies. The most significant disagreement found between M1 and Monte-Carlo results is a factor of $\sim 2$ difference in the luminosity of heavy-lepton neutrinos.


title: The Formation of Very Massive Stars in Early Galaxies and Implications for Intermediate Mass Black Holes

abstract: We investigate the ab-initio formation of super-massive stars in a pristine atomic cooling halo. The halo is extracted from a larger self-consistent parent simulation. The halo remains metal-free and star formation is suppressed due to a combination of dynamical heating from mergers and a mild ($J_{\rm LW} \sim 2 - 10 \ J_{21}$(z)) Lyman-Werner (LW) background. We find that more than 20 very massive stars form with stellar masses greater than 1000 M$_{\odot}$. The most massive star has a stellar mass of over 6000 M$_{\odot}$. However, accretion onto all stars declines significantly after the first $\sim$ 100 kyr of evolution as the surrounding material is accreted and the turbulent nature of the gas causes the stars to move to lower density regions. We post-process the impact of ionising radiation from the stars and find that ionising radiation is not a limiting factor when considering SMS formation and growth. Rather the birth environments are highly turbulent and a steady accretion flow is not maintained within the timescale (2 Myr) of our simulations. As the massive stars end their lives as direct collapse black holes this will seed these embryonic haloes with a population of black holes with masses between approximately 300 M$_{\odot}$ and 10,000 M$_{\odot}$. Afterwards they may sink to the centre of the haloes, eventually coalescing to form larger intermediate mass black holes whose in-situ mergers will be detectable by LISA.


title: NGC 2770: high supernova rate due to interaction

abstract: Galaxies which hosted many core-collapse supernovae (SN) explosions can be used to study the conditions necessary for the formation of massive stars. NGC 2770 was dubbed a SN factory, because it hosted four core-collapse SNe in 20 years (three type Ib and one type IIn). Its star formation rate (SFR) was reported not to be enhanced and therefore not compatible with such a high SN rate. We aim at explaining the high SN rate of NGC 2770. We used archival HI line data for NGC 2770 and reinterpret the Halpha and optical continuum data. Even though the continuum-based SFR indicators do not yield high values, the dust-corrected Halpha luminosity implies a high SFR, consistent with the high SN rate. Such disparity between the SFR estimators is an indication of recently enhanced star formation activity, because the continuum indicators trace long timescale of the order of 100 Myr, unlike the line indicators, which trace timescales of the order of 10 Myr. Hence, the unique feature of NGC 2770 compared to other galaxies is the fact that it is observed very recently after the enhancement of the SFR. It also has high dust extinction, E(B-V) above 1 mag. We provide support for the hypothesis that the increased SFR in NGC 2770 is due to the interaction with its companion galaxies. We report an HI bridge between NGC 2770 and its closest companion and the existence of a total of four companions within 100 kpc (one identified for the first time). There are no clear HI concentrations close to the positions of SNe in NGC 2770 such as those detected for hosts of gamma-ray bursts (GRBs) and broad-lined SNe type Ic (IcBL). This suggests that the progenitors of type Ib SNe are not born out of recently accreted atomic gas, as was suggested for GRB and IcBL SN progenitors.


title: CGM$^{2}$ I: The Extent of the Circumgalactic Medium Traced by Neutral Hydrogen

abstract: We present initial results from the \textit{COS and Gemini Mapping the Circumgalactic Medium} (\mbox{CGMCGM} $\equiv$ CGM$^{2}$) survey. The CGM$^{2}$ survey consists of 1689 galaxies, all with high-quality Gemini GMOS spectra, within 1 Mpc of twenty-two $z \lesssim 1$ quasars, all with S/N$\sim$10 {\emph{HST/COS}} G130M$+$G160M spectra. For 572 of these galaxies having stellar masses $10^{7} M_{\odot} < M_{\star} < 10^{11} M_{\odot}$ and $z \lesssim 0.5$, we show that the \ion{H}{1} covering fraction above a threshold of \NHI$>10^{14} $cm$^{-2}$ is $\gtrsim 0.5$ within 1.5 virial radii ($R_{\rm vir} \sim R_{200m}$). We examine the \ion{H}{1} kinematics and find that the majority of absorption lies within $\pm$ 250 km s$^{-1}$ of the galaxy systemic velocity. We examine \ion{H}{1} covering fractions over a range of impact parameters to infer a characteristic size of the CGM, $R^{14}_{\rm CGM}$, as a function of galaxy mass. $R^{14}_{\rm CGM}$ is the impact parameter at which the probability of observing an absorber with \NHI $>$ 10$^{14}$ cm$^{-2}$ is $>$ 50\%. In this framework, the radial extent of the CGM of $M_{\star} > 10^{9.9} M_{\odot}$ galaxies is $R^{14}_{\rm CGM} = 346^{+57}_{-53}$ kpc or $R^{14}_{\rm CGM} \simeq 1.2R_{\rm vir}$. Intermediate-mass galaxies with $10^{9.2} < M_{\star}/M_{\odot} < 10^{9.9}$ have an extent of $R^{14}_{\rm CGM} = 353^{+64}_{-50}$ kpc or $R^{14}_{\rm CGM} \simeq 2.4R_{\rm vir}$. Low-mass galaxies, $M_{\star} < 10^{9.2} M_{\odot}$, show a smaller physical scale $R^{14}_{\rm CGM} = 177_{-65}^{+70}$ kpc and extend to $R^{14}_{\rm CGM} \simeq 1.6R_{\rm vir}$. Our analysis suggests that using $R_{\rm vir}$ as a proxy for the characteristic radius of the CGM likely underestimates its extent.


title: Machine Learning Approach to Integral Field Unit Spectroscopy Observations: I. HII Region Kinematics

abstract: SITELLE is a novel integral field unit spectroscopy instrument that has an impressive spatial (11 by 11 arcmin), spectral coverage, and spectral resolution (R=1-20000). SIGNALS is anticipated to obtain deep observations (down to 3.6x10-17ergs s-1cm-2) of 40 galaxies, each needing complex and substantial time to extract spectral information. We present a method that uses Convolution Neural Networks (CNN) for estimating emission line parameters in optical spectra obtained with SITELLE as part of the SIGNALS large program. Our algorithm is trained and tested on synthetic data representing typical emission spectra for HII regions based on Mexican Million Models database(3MdB) BOND simulations. The network's activation map demonstrates its ability to extract the dynamical (broadening and velocity) parameters from a set of 5 emission lines (e.g. H{\alpha}, N[II] doublet, and S[II] doublet) in the SN3 (651-685 nm) filter of SITELLE. Once trained, the algorithm was tested on real SITELLE observations in the SIGNALS program of one of the South West fields of M33. The CNN recovers the dynamical parameters with an accuracy better than 5 km s-1 in regions with a signal-to-noise ratio greater than 15 over the H{\alpha}line. More importantly, our CNN method reduces calculation time by over an order of magnitude on the spectral cube with native spatial resolution when compared with standard fitting procedures. These results clearly illustrate the power of machine learning algorithms for the use in future IFU-based missions. Subsequent work will explore the applicability of the methodology to other spectral parameters such as the flux of key emission lines.


title: Demonstration of an amplitude filter cavity at gravitational-wave frequencies

abstract: Quantum vacuum fluctuations fundamentally limit the precision of optical measurements, such as those in gravitational-wave detectors. Injection of conventional squeezed vacuum can be used to reduce quantum noise in the readout quadrature, but this reduction is at the cost of increasing noise in the orthogonal quadrature. For detectors near the limits imposed by quantum radiation pressure noise (QRPN), both quadratures impact the measurement, and the benefits of conventional squeezing are limited. In this paper, we demonstrate the use of a critically-coupled 16m optical cavity to diminish anti-squeezing at frequencies below 90Hz where it exacerbates QRPN, while preserving beneficial squeezing at higher frequencies. This is called an amplitude filter cavity, and it is useful for avoiding degradation of detector sensitivity at low frequencies. The attenuation from the cavity also provides technical advantages such as mitigating backscatter.


title: The singularities and Kodaira dimension of unitary Shimura varieties

abstract: The Kodaira dimension of Shimura varieties has been studied by many people. Kondo and Gritsenko-Hulek-Sankaran studied the singularities of orthogonal Shimura varieties related to the moduli spaces of polarized K3 surfaces. They proved that they have canonical singularities and are of general type if the polarization degree is sufficiently large. In this paper, we work on similar problems for unitary Shimura varieties. We show that they have canonical singularities if $n > 4$. As an application, we show that certain unitary Shimura varieties associated with Hermitian forms over the rings of integers of $\mathbb{Q}(\sqrt{-1})$, $\mathbb{Q}(\sqrt{-3})$ are of general type. We use modular forms of low weight vanishing on ramification divisors, which are the restrictions of the quasi-pullbacks of the Borcherds form $\Phi_{12}$.


title: A Carbon-enhanced Lyman Limit System: Signature of the First Generation of Stars?

abstract: We present the study of a Lyman limit system (LLS) at $z_{\rm abs}$ = 1.5441 towards quasar J134122.50+185213.9 observed with VLT X-shooter. This is a very peculiar system with strong C I absorption seen associated with a neutral hydrogen column density of log $N$(H I) (cm$^{-2}$) = 18.10, too small to shield the gas from any external UV flux. The low ionization absorption lines exhibit a simple kinematic structure consistent with a single component. Using CLOUDY models to correct for ionization, we find that the ionization parameter of the gas is in the range $-$ 4.5 $<$ log $U$ $<$ $-$4.2 and the gas density $-$1.5 $<$ log $n$(H) (cm$^{-3}$) $<$ $-$1.2. The models suggest that carbon is overabundant relative to iron, [C/Fe] $>$ +2.2 at [Fe/H] $\sim$ $-$1.6. Such a metal abundance pattern is reminiscent of carbon-enhanced metal-poor stars detected in the Galaxy halo. Metal enrichment by the first generation of supernovae provides a plausible explanation for the inferred abundance pattern in this system.


title: Disc-binary interactions in depleted post-AGB binaries

abstract: Binary post-asymptotic giant branch (post-AGB) stars have orbital periods in the range of 100--2500 days in eccentric orbits. They are surrounded by circumbinary dusty discs. They are the immediate result of unconstrained binary interaction processes. Their observed orbital properties do not correspond to model predictions: Neither the periods nor the high eccentricities are expected. Our goal is to investigate if interactions between a binary and its circumbinary disc during the post-AGB phase can result in their eccentric orbits, while simultaneously explaining the chemical anomaly known as depletion. For this paper, we selected three binaries (EP Lyr, RU Cen, HD 46703) with well-constrained orbits, luminosities, and chemical abundances. We used the MESA code to evolve post-AGB models, while including the accretion of metal-poor gas. This allows us to constrain the evolution of the stars and study the impact of circumbinary discs on the orbital properties of the models. We investigate the effect of torques produced by gas inside the binary cavity and the effect of Lindblad resonances on the orbit, while also including the tidal interaction following the equilibrium tide model. We find that none of our models are able to explain the high orbital eccentricities of the binaries in our sample. The accretion torque does not significantly impact the binary orbit, while Lindblad resonances can pump the eccentricity up to only $e \approx 0.2$. At higher eccentricities, the tidal interaction becomes too strong, so the high observed eccentricities cannot be reproduced. However, even if we assume tides to be ineffective, the eccentricities in our models do not exceed $\approx 0.25$. We conclude that either our knowledge of disc-binary interactions is still incomplete, or the binaries must have left their phase of strong interaction in an eccentric orbit.


title: Two-Particle Self-Consistent method for the multi-orbital Hubbard model

abstract: One of the most challenging problems in solid state systems is the microscopic analysis of electronic correlations. A paramount minimal model that encodes correlation effects is the Hubbard Hamiltonian, which -- albeit its simplicity -- is exactly solvable only in a few limiting cases and approximate many-body methods are required for its solution. In this review we present an overview on the non-perturbative Two-Particle Self-Consistent method (TPSC) which was originally introduced to describe the electronic properties of the single-band Hubbard model. We introduce here a detailed derivation of the multi-orbital generalization of TPSC and discuss particular features of the method on exemplary interacting models in comparison to dynamical mean-field theory results.


title: The LOFAR view of FR0 radio galaxies

abstract: We explore the low-frequency radio properties of the sources in the Fanaroff-Riley class 0 catalog (FR0CAT) as seen by the LOw Frequency ARray (LOFAR) observations at 150 MHz. This sample includes 104 compact radio active galactic nuclei (AGN) associated with nearby (z<0.05) massive early-type galaxies. Sixty-six FR0CAT sources are in the sky regions observed by LOFAR and all of them are detected, usually showing point-like structures with sizes smaller than 3-6 kpc. However, 12 FR0s present resolved emission of low surface brightness which contributes between 5% and 40% of the total radio power at 150 MHz, usually with a jetted morphology extending between 15 and 50 kpc. No extended emission is detected around the other FR0s, with a typical luminosity limit of 5 x 10$^{22}$ W/Hz over an area of 100 kpc x 100 kpc. The spectral slopes of FR0s between 150 MHz and 1.4 GHz span a broad range (-0.7 < $\alpha$ < 0.8) with a median value of $\overline\alpha \sim 0.1$; 20% of them have a steep spectrum ($\alpha$ > 0.5), an indication of the presence of substantial extended emission confined within the spatial resolution limit. The fraction of FR0s showing evidence for the presence of jets, by including both spectral and morphological information, is at least ~40%. This study confirms that FR0s and FRIs can be interpreted as two extremes of a continuous population of jetted sources, with the FR0s representing the low end in size and radio power.


title: Reconciling hints on axion-like-particles from high-energy gamma rays with stellar bounds

abstract: It has been recently claimed by two different groups that the spectral modulation observed in gamma rays from Galactic pulsars and supernova remnants can be due to conversion of photons into ultra-light axion-like-particles (ALPs) in large-scale Galactic magnetic fields. While we show the required best-fit photon-ALP coupling, $g_{a\gamma} \sim 2 \times 10^{-10}$ GeV${}^{-1}$, to be consistent with constraints from observations of photon-ALPs mixing in vacuum, this is in conflict with other bounds, specifically from the CAST solar axion limit, from the helium-burning lifetime in globular clusters, and from the non-observations of gamma rays in coincidence with SN 1987A. In order to reconcile these different results, we propose that environmental effects in matter would suppress the ALP production in dense astrophysical plasma, allowing to relax previous bounds and make them compatible with photon-ALP conversions in the low-density Galactic medium. If this explanation is correct, the claimed ALP signal would be on the reach of next-generations laboratory experiments such as ALPS II.


title: ALMA 0.02"-resolution observations reveal HCN-abundance-enhanced counter-rotating and outflowing dense molecular gas at the NGC 1068 nucleus

abstract: We present ALMA ~0.02"-resolution observations of the nucleus of the nearby (~14 Mpc) type-2 AGN NGC 1068 at HCN/HCO+/HNC J=3-2 lines, as well as at their 13C isotopologue and vibrationally excited lines, to scrutinize the morphological/dynamical/chemical/physical properties of dense molecular gas in the putative dusty molecular torus around a mass-accreting supermassive black hole. We confirm almost east-west-oriented dense molecular gas emission both morphologically and dynamically, which we regard as coming from the torus. Bright emission is compact (<3 pc), and low-surface-brightness emission extends out to 5-7 pc. These dense molecular gas properties are not symmetric between the eastern and western torus. The HCN J=3-2 emission is stronger than the HCO+ J=3-2 emission within the ~7 pc torus region, with an estimated dense molecular mass of (0.4-1.0)x10^6Msun. We interpret that HCN abundance is enhanced in the torus. We detect signatures of outflowing dense molecular gas and a vibrationally excited HCN J=3-2 line. Finally, we find that in the innermost (<1 pc) part of the torus, the dense molecular line rotation velocity, relative to the systemic velocity, is the opposite of that in the outer (>2 pc) part, in both the eastern and western torus. We prefer a scenario of counter-rotating dense molecular gas with innermost almost-Keplerian-rotation and outer slowly rotating (far below Keplerian) components. Our high-spatial-resolution dense molecular line data reveal that torus properties of NGC 1068 are much more complicated than the simple axi-symmetrically rotating torus picture in the classical AGN unification paradigm.


title: Effects of defects in the XY chain with frustrated boundary conditions

abstract: It has been recently proven that new types of bulk, local order can ensue due to frustrated boundary condition, that is, periodic boundary conditions with an odd number of lattice sites and anti-ferromagnetic interactions. For the quantum XY chain in zero external fields, the usual antiferromagnetic order has been found to be replaced either by a mesoscopic ferromagnet or by an incommensurate AFM order. In this work we examine the resilience of these new types of orders against a defect that breaks the translational symmetry of the model. We find that, while a ferromagnetic defect restores the traditional, staggered order, an AFM one stabilizes the incommensurate order. The robustness of the frustrated order to certain kinds of defects paves the way for its experimental observability.

categories: cond-mat.str-el cond-mat.stat-mech hep-th quant-ph


title: On the Numerical Solution of Nonlinear Eigenvalue Problems for the Monge-Amp\`{e}re Operator

abstract: In this article, we report the results we obtained when investigating the numerical solution of some nonlinear eigenvalue problems for the Monge-Amp\`{e}re operator $v\rightarrow \det \mathbf{D}^2 v$. The methodology we employ relies on the following ingredients: (i) A divergence formulation of the eigenvalue problems under consideration. (ii) The time discretization by operator-splitting of an initial value problem (a kind of gradient flow) associated with each eigenvalue problem. (iii) A finite element approximation relying on spaces of continuous piecewise affine functions. To validate the above methodology, we applied it to the solution of problems with known exact solutions: The results we obtained suggest convergence to the exact solution when the space discretization step $h\rightarrow 0$. We considered also test problems with no known exact solutions.


title: The role of the chiral anomaly in polarized deeply inelastic scattering I: Finding the triangle graph inside the box diagram in Bjorken and Regge asymptotics

abstract: We revisit the role of the chiral "triangle" anomaly in deeply inelastic scattering (DIS) of electrons off polarized protons employing a powerful worldline formalism. We demonstrate how the triangle anomaly appears at high energies in the DIS box diagram for the polarized proton structure function $g_1(x_B,Q^2)$ in both the Bjorken limit of large $Q^2$ and in the Regge limit of small $x_B$. We show that the operator product expansion is not required to extract the anomaly in either asymptotics though it is sufficient in the Bjorken limit. Likewise, the infrared pole in the anomaly arises in both limits. The leading contribution to $g_1$, in both Bjorken and Regge asymptotics, is therefore given by the expectation value of the topological charge density, generalizing a result previously argued by Jaffe and Manohar to hold for the first moment of $g_1$. In follow-up work, we will show how our results motivate the derivation of a helicity-dependent effective action incorporating the physics of the anomaly at small $x_B$ and shall discuss the QCD evolution of $g_1(x_B,Q^2)$ in this framework.


title: On the steady motion of Navier-Stokes flows past a fixed obstacle in a three-dimensional channel under mixed boundary conditions

abstract: We analyze the steady motion of a viscous incompressible fluid in a three-dimensional channel containing an obstacle through the Navier-Stokes equations with mixed boundary conditions: the inflow is given by a fairly general datum and the flow is assumed to satisfy a constant traction boundary condition on the outlet, together with the standard no-slip assumption on the obstacle and on the remaining walls of the domain. Explicit bounds on the inflow velocity guaranteeing existence and uniqueness of such steady motion are provided after estimating some Sobolev embedding constants and constructing a suitable solenoidal extension of the inlet velocity through the Bogovskii formula. A quantitative analysis of the forces exterted by the fluid over the obstacle constitutes the main application of our results: by deriving a volume integral formula for the drag and lift, explicit upper bounds on these forces are given in terms of the geometrical constraints of the domain.


title: Temperature structures of embedded disks: young disks in Taurus are warm

abstract: The chemical composition of gas and ice in disks around young stars set the bulk composition of planets. In contrast to protoplanetary disks (Class II), young disks that are still embedded in their natal envelope (Class 0 and I) are predicted to be too warm for CO to freeze out, as has been confirmed observationally for L1527 IRS. To establish whether young disks are generally warmer than their more evolved counterparts, we observed five young (Class 0/I and Class I) disks in Taurus with the Atacama Large Millimeter/submillimeter Array (ALMA), targeting C$^{17}$O $2-1$, H$_2$CO $3_{1,2}-2_{1,1}$, HDO $3_{1,2}-2_{2,1}$ and CH$_3$OH $5_K-4_K$ transitions at $0.48^{\prime\prime} \times 0.31^{\prime\prime}$ resolution. The different freeze-out temperatures of these species allow us to derive a global temperature structure. C$^{17}$O and H$_2$CO are detected in all disks, with no signs of CO freeze-out in the inner $\sim$100 au, and a CO abundance close to $\sim$10$^{-4}$. H$_2$CO emission originates in the surface layers of the two edge-on disks, as witnessed by the especially beautiful V-shaped emission pattern in IRAS~04302+2247. HDO and CH$_3$OH are not detected, with column density upper limits more than 100 times lower than for hot cores. Young disks are thus found to be warmer than more evolved protoplanetary disks around solar analogues, with no CO freeze-out (or only in the outermost part of $\gtrsim$100 au disks) or CO processing. However, they are not as warm as hot cores or disks around outbursting sources, and therefore do not have a large gas-phase reservoir of complex molecules.


title: Impact of correlations between $a_\mu$ and $\alpha_{QED}$ on the EW fit

abstract: We study the potential impact on the electroweak (EW) fits due to the tensions between the current determinations of the hadronic vacuum polarisation (HVP) contributions to the anomalous magnetic moment of the muon ($a_\mu$), based on either phenomenological dispersion integrals using measured hadronic spectra or on Lattice QCD calculations. The impact of the current tension between the experimental measurement of $a_\mu$ and the total theoretical prediction based on the phenomenological calculations of the HVP are also studied. The correlations between the uncertainties of the theoretical predictions of $a_\mu$ and of the running of $\alpha_{QED}$ are taken into account in the studies. We conclude that the impact on the EW fit can be large in improbable scenarios involving global shifts of the full HVP contribution, while it is much smaller if the shift is restricted to a lower mass range and/or if the shift in $\alpha_{QED}$ is obtained from that in $a_\mu$ through appropriate use of the correlations. Indeed, the latter scenarios only imply at most a 2.6/16 increase in the $\chi^2/n.d.f.$ of the EW fits and relatively small changes for the resulting fit parameter values.


title: Dark Scalars and Heavy Neutral Leptons at DarkQuest

abstract: The proposed DarkQuest beam dump experiment, a modest upgrade to the existing SeaQuest/SpinQuest experiment, has great potential for uncovering new physics within a dark sector. We explore both the near-term and long-term prospects for observing two distinct, highly-motivated hidden sector benchmark models: heavy neutral leptons and Higgs-mixed scalars. We comprehensively examine the particle production and detector acceptance at DarkQuest, including an updated treatment of meson production, and light scalar production through both bremsstrahlung and gluon-gluon fusion. In both benchmark models, DarkQuest will provide an opportunity to probe previously inaccessible interesting regions of parameter space on a fairly short timescale when compared to other proposed experiments.


title: Local-density dependent Markov processes on graphons with epidemiological applications

abstract: We investigate local-density dependent Markov processes on a class of large graphs sampled from a graphon, where the transition rates of the vertices are influenced by the states of their neighbors. We show that as the average degree converges to infinity, the evolution of the process in the transient regime converges to the solution of a set of non-local integro-partial differential equations. We also provide rigorous derivation for the epidemic threshold in the case of the Susceptible-Infected-Susceptible (SIS) process on such graphons.


title: Numerical Semigroups of small and large type

abstract: A numerical semigroup is a sub-semigroup of the natural numbers that has a finite complement. Some of the key properties of a numerical semigroup are its Frobenius number F, genus g and type t. It is known that for any numerical semigroup $\frac{g}{F+1-g}\leq t\leq 2g-F$. Numerical semigroups with $t=2g-F$ are called almost symmetric, we introduce a new property that characterises them. We give an explicit characterisation of numerical semigroups with $t=\frac{g}{F+1-g}$. We show that for a fixed $\alpha$ the number of numerical semigroups with Frobenius number $F$ and type $F-\alpha$ is eventually constant for large $F$. Also the number of numerical semigroups with genus $g$ and type $g-\alpha$ is also eventually constant for large $g$.


title: Splitting methods for solution decomposition in nonstationary problems

abstract: In approximating solutions of nonstationary problems, various approaches are used to compute the solution at a new time level from a number of simpler (sub-)problems. Among these approaches are splitting methods. Standard splitting schemes are based on one or another additive splitting of the operator into "simpler" operators that are more convenient/easier for the computer implementation and use inhomogeneous (explicitly-implicit) time approximations. In this paper, a new class of splitting schemes is proposed that is characterized by an additive representation of the solution instead of the operator corresponding to the problem (called problem operator). A specific feature of the proposed splitting is that the resulting coupled equations for individual solution components consist of the time derivatives of the solution components. The proposed approaches are motivated by various applications, including multiscale methods, domain decomposition, and so on, where spatially local problems are solved and used to compute the solution. Unconditionally stable splitting schemes are constructed for a first-order evolution equation, which is considered in a finite-dimensional Hilbert space. In our splitting algorithms, we consider the decomposition of both the main operator of the system and the operator at the time derivative. Our goal is to provide a general framework that combines temporal splitting algorithms and spatial decomposition and its analysis. Applications of the framework will be studied separately.


title: A Theory Building Study of Enterprise Architecture Practices and Benefits

abstract: Academics and practitioners have made various claims regarding the benefits that Enterprise Architecture (EA) delivers for both individual projects and the organization as a whole. At the same time, there is a lack of explanatory theory regarding how EA delivers these benefits. Moreover, EA practices and benefits have not been extensively investigated by empirical research, with especially quantitative studies on the topic being few and far between. This paper therefore presents the statistical findings of a theory-building survey study (n=293). The resulting PLS model is a synthesis of current implicit and fragmented theory, and shows how EA practices and intermediate benefits jointly work to help the organization reap benefits for both the organization and its projects. The model shows that EA and EA practices do not deliver benefits directly, but operate through intermediate results, most notably compliance with EA and architectural insight. Furthermore, the research identifies the EA practices that have a major impact on these results, the most important being compliance assessments, management propagation of EA, and different types of knowledge exchange. The results also demonstrate that projects play an important role in obtaining benefits from EA, but that they generally benefit less than the organization as a whole.


title: Complementary Language Model and Parallel Bi-LRNN for False Trigger Mitigation

abstract: False triggers in voice assistants are unintended invocations of the assistant, which not only degrade the user experience but may also compromise privacy. False trigger mitigation (FTM) is a process to detect the false trigger events and respond appropriately to the user. In this paper, we propose a novel solution to the FTM problem by introducing a parallel ASR decoding process with a special language model trained from "out-of-domain" data sources. Such language model is complementary to the existing language model optimized for the assistant task. A bidirectional lattice RNN (Bi-LRNN) classifier trained from the lattices generated by the complementary language model shows a $38.34\%$ relative reduction of the false trigger (FT) rate at the fixed rate of $0.4\%$ false suppression (FS) of correct invocations, compared to the current Bi-LRNN model. In addition, we propose to train a parallel Bi-LRNN model based on the decoding lattices from both language models, and examine various ways of implementation. The resulting model leads to further reduction in the false trigger rate by $10.8\%$.


title: Commonsense Knowledge in Wikidata

abstract: Wikidata and Wikipedia have been proven useful for reason-ing in natural language applications, like question answering or entitylinking. Yet, no existing work has studied the potential of Wikidata for commonsense reasoning. This paper investigates whether Wikidata con-tains commonsense knowledge which is complementary to existing commonsense sources. Starting from a definition of common sense, we devise three guiding principles, and apply them to generate a commonsense subgraph of Wikidata (Wikidata-CS). Within our approach, we map the relations of Wikidata to ConceptNet, which we also leverage to integrate Wikidata-CS into an existing consolidated commonsense graph. Our experiments reveal that: 1) albeit Wikidata-CS represents a small portion of Wikidata, it is an indicator that Wikidata contains relevant commonsense knowledge, which can be mapped to 15 ConceptNet relations; 2) the overlap between Wikidata-CS and other commonsense sources is low, motivating the value of knowledge integration; 3) Wikidata-CS has been evolving over time at a slightly slower rate compared to the overall Wikidata, indicating a possible lack of focus on commonsense knowledge. Based on these findings, we propose three recommended actions to improve the coverage and quality of Wikidata-CS further.


title: TIDE: A General Toolbox for Identifying Object Detection Errors

abstract: We introduce TIDE, a framework and associated toolbox for analyzing the sources of error in object detection and instance segmentation algorithms. Importantly, our framework is applicable across datasets and can be applied directly to output prediction files without required knowledge of the underlying prediction system. Thus, our framework can be used as a drop-in replacement for the standard mAP computation while providing a comprehensive analysis of each model's strengths and weaknesses. We segment errors into six types and, crucially, are the first to introduce a technique for measuring the contribution of each error in a way that isolates its effect on overall performance. We show that such a representation is critical for drawing accurate, comprehensive conclusions through in-depth analysis across 4 datasets and 7 recognition models. Available at https://dbolya.github.io/tide/


title: Phase Transitions in Asymptotically Singular Anderson Hamiltonian and Parabolic Model

abstract: Let $\xi$ be a Gaussian white noise on $\mathbb R^d$ ($d=1,2,3$). Let $(\xi_\varepsilon)_{\varepsilon>0}$ be continuous Gaussian processes such that $\xi_\varepsilon\to\xi$ as $\varepsilon\to0$, defined by convolving $\xi$ against a mollifier. We consider the asymptotics of the parabolic Anderson model (PAM) with noise $\xi_{\varepsilon(t)}$ for large time $t\gg1$, and the Dirichlet eigenvalues of the Anderson Hamiltonian (AH) with potential $\xi_{\varepsilon(t)}$ on large boxes $(-t,t)^d$, where the parameter $\varepsilon(t)$ vanishes as $t\to\infty$. We prove that the asymptotics in question exhibit a phase transition in the rate at which $\varepsilon(t)$ vanishes, which distinguishes between the behavior observed in the AH/PAM with continuous Gaussian noise and white noise. By comparing our main theorems with previous results on the AH/PAM with white noise, our results show that some asymptotics of the latter can be accessed with solely elementary methods, and we obtain quantitative estimates on the difference between the AH/PAM with white noise and its continuous-noise approximations as $t\to\infty$.


title: Bounds on Distributional Treatment Effect Parameters using Panel Data with an Application on Job Displacement

abstract: This paper develops new techniques to bound distributional treatment effect parameters that depend on the joint distribution of potential outcomes -- an object not identified by standard identifying assumptions such as selection on observables or even when treatment is randomly assigned. I show that panel data and an additional assumption on the dependence between untreated potential outcomes for the treated group over time (i) provide more identifying power for distributional treatment effect parameters than existing bounds and (ii) provide a more plausible set of conditions than existing methods that obtain point identification. I apply these bounds to study heterogeneity in the effect of job displacement during the Great Recession. Using standard techniques, I find that workers who were displaced during the Great Recession lost on average 34\% of their earnings relative to their counterfactual earnings had they not been displaced. Using the methods developed in the current paper, I also show that the average effect masks substantial heterogeneity across workers.


title: Review of Machine Learning Algorithms for Brain Stroke Diagnosis and Prognosis by EEG Analysis

abstract: Currently, strokes are the leading cause of adult disability in the United States. Traditional treatment and rehabilitation options such as physical therapy and tissue plasminogen activator are limited in their effectiveness and ability to restore mobility and function to the patient. As a result, there exists an opportunity to greatly improve the treatment for strokes. Machine learning, specifically techniques that utilize Brain-Computer Interfaces (BCIs) to help the patient either restore neurologic pathways or effectively communicate with an electronic prosthetic, show promising results when applied to both stroke diagnosis and rehabilitation. In this review, sources that design and implement BCIs for treatment of stroke patients are evaluated and categorized based on their successful applications for stroke diagnosis or stroke rehabilitation. The various machine learning techniques and algorithms that are addressed and combined with BCI technology show that the use of BCIs for stroke treatment is a promising and rapidly expanding field.


title: Beyond the MSW effect: Neutrinos in a dense medium

abstract: We present a theory of neutrino oscillations in a dense medium which goes beyond the effective matter potential used in the description of the MSW effect. We show how the purity of the neutrino state is degraded by neutrino interactions with the environment and how neutrino--matter interactions can be a source of decoherence. We present new oscillation formulae for neutrinos interacting with leptons and carry out a numerical analysis which exhibits deviations from the MSW formulae for propagation through the Earth of ultra-high energy neutrinos. In particular, we show that at high density and/or high neutrino energy, the vanishing transition probabilities derived for MSW effect, are non zero when the scattering is taken into account.


title: Smooth loops and loop bundles

abstract: A loop is a rather general algebraic structure that has an identity element and division, but is not necessarily associative. Smooth loops are a direct generalization of Lie groups. A key example of a non-Lie smooth loop is the loop of unit octonions. In this paper, we study properties of smooth loops and their associated tangent algebras, including a loop analog of the Mauer-Cartan equation. Then, given a manifold, we introduce a loop bundle as an associated bundle to a particular principal bundle. Given a connection on the principal bundle, we define the torsion of a loop bundle structure and show how it relates to the curvature, and also consider the critical points of some related functionals. Throughout, we see how some of the known properties of $G_{2}$-structures can be seen from this more general setting.


title: Near-Zero-Field Spin-Dependent Recombination Current and Electrically Detected Magnetic Resonance from the Si/SiO$_2$ interface

abstract: Dielectric interfaces critical for metal-oxide-semiconductor (MOS) electronic devices, such as the Si/SiO$_2$ MOS field effect transistor (MOSFET), possess trap states that can be visualized with electrically-detected spin resonance techniques, however the interpretation of such measurements has been hampered by the lack of a general theory of the phenomena. This article presents such a theory for two electrical spin-resonance techniques, electrically detected magnetic resonance (EDMR) and the recently observed near-zero field magnetoresistance (NZFMR), by generalizing Shockley Read Hall trap-assisted recombination current calculations via stochastic Liouville equations. Spin mixing at this dielectric interface occurs via the hyperfine interaction, which we show can be treated either quantum mechanically or semiclassically, yielding distinctive differences in the current across the interface. By analyzing the bias dependence of NZFMR and EDMR, we find that the recombination in a Si/SiO$_2$ MOSFET is well understood within a semiclassical approach.


title: $L_\infty$-structure on Barzdell's complex for monomial algebras

abstract: Let $A$ be a monomial associative finite dimensional algebra over a field $\Bbbk$ of characteristic zero. It is well known that the Hochschild cohomology of $A$ can be computed using Bardzell's complex $B(A)$. The aim of this article is to describe an explict $L_\infty$-structure on $B(A)$ that induces a weak equivalence of $L_\infty$-algebras between $B(A)$ and the Hochschild complex $C(A)$ of $A$. This allows us to describe the Maurer-Cartan equation in terms of elements of degree $2$ in $B(A)$. Finally, we make concrete computations when $A$ is a truncated algebra, and we prove that Bardzell's complex for radical square zero algebras is in fact a dg-Lie algebra.


title: On the Boundedness of Globally $F$-split varieties

abstract: This paper proposes the use of $F$-split and globally $F$-regular conditions in the pursuit of BAB type results in positive characteristic. The main technical work comes in the form of a detailed study of threefold Mori fibre spaces over positive dimensional bases. As a consequence we prove the main theorem, which reduces birational boundedness for a large class of varieties to the study of prime Fano varieties.


title: WloopPHI: A tool for ab initio characterization of Weyl semimetals

abstract: WloopPHI is a Python code that expands the features of WIEN2k, a full-potential all-electron density functional theory package, by the characterization of Weyl semimetals. It enables the calculation of the chirality (or "monopole charge") associated with Weyl nodes and nodal lines. The theoretical methodology for the calculation of the chirality is based on an extended Wilson loop method and a Berry phase approach. We validate the code using TaAs, which is a well-characterized Weyl semimetal, both theoretically and experimentally. Afterwards, we applied the method to the characterization of YRh$_6$Ge$_4$ and found two sets of Weyl points (ca. 0.2 eV below the Fermi energy) together with a topological nodal line (protected by mirror symmetry) crossing the Fermi energy and mapped their chiralities.


title: Abelian Closures of Infinite Binary Words

abstract: Two finite words $u$ and $v$ are called Abelian equivalent if each letter occurs equally many times in both $u$ and $v$. The abelian closure $\mathcal{A}(\mathbf{x})$ of (the shift orbit closure of) an infinite word $\mathbf{x}$ is the set of infinite words $\mathbf{y}$ such that, for each factor $u$ of $\mathbf{y}$, there exists a factor $v$ of $\mathbf{x}$ which is abelian equivalent to $u$. The notion of an abelian closure gives a characterization of Sturmian words: among binary uniformly recurrent words, Sturmian words are exactly those words for which $\mathcal{A}(\mathbf{x})$ equals the shift orbit closure $\Omega(\mathbf{x})$. In this paper we show that, contrary to larger alphabets, the abelian closure of a uniformly recurrent aperiodic binary word which is not Sturmian contains infinitely many minimal subshifts.


title: $Z$-knotted and $Z$-homogeneous triangulations of surfaces

abstract: A triangulation is called $z$-knotted if it has a single zigzag (up to reversing). A $z$-orientation on a triangulation is a minimal collection of zigzags which double covers the set of edges. An edge is of type I if zigzags from the $z$-orientation pass through it in different directions, otherwise this edge is of type II. If all zigzags from the $z$-orientation contain precisely two edges of type I after any edge of type II, then the $z$-oriented triangulation is said to be $z$-homogeneous. We describe an algorithm transferring each $z$-homogeneous trianguation to other $z$-homogeneous triangulation which is also $z$-knotted.


title: Research on rolling friction's dependence on ball bearings' radius

abstract: There are two alternative historical laws of rolling resistance formulated by French scientist Coulomb and Dupuit. It has been decided to verify experimentally again, which of these laws describes freely rolling ball bearings on a hard surface better. An inducement to carrying out the measurements was the idea of the constant thickness of roadbed, which is consistent with Dupuit's theory. Measurements have been done using the damped oscillations in the pendulum bearings. Results have shown better consistency with the Coulomb's theory with small, but measurable deviations. These deviations were successfully explained by the so called "Cobblestones model". Parameters designated by this model have been successfully verified by the surface roughness's profile measurement. An additional theoretical aspect of this work is distinguishing two types of rolling friction force: dynamical and kinematical in an analogy to two types of specific heat capacity in the thermodynamics of gases.


title: Explicit Constructions of Halphen Pencils

abstract: We construct rational elliptic surfaces of index two by explicitly constructing their associated Halphen pencils in the projective plane $\mathbb{P}^2$. For each of the types of singular fibers that occur we construct at least one example having that type of fiber and in fact, for some, we construct all possible examples. We establish a precise dictionary between the fibers in a rational elliptic surface and the corresponding plane curves and, in particular, we study the singularities of the curves appearing in a Halphen pencil.


title: Graphene Bilayers with a Twist

abstract: Near a magic twist angle, bilayer graphene transforms from a weakly correlated Fermi liquid to a strongly correlated two-dimensional electron system with properties that are extraordinarily sensitive to carrier density and to controllable environmental factors such as the proximity of nearby gates and twist-angle variation. Among other phenomena magic-angle twisted bilayer graphene hosts superconductivity, interaction induced insulating states, magnetism, electronic nematicity, linear-in-T low-temperature resistivity, and quantized anomalous Hall states. We highlight some key research results in this field, point to important questions that remain open, and comment on the place of magic angle twisted bilayer graphene in the strongly correlated quantum matter world.


title: Single bubble and drop techniques for characterizing foams and emulsions

abstract: The physics of foams and emulsions has traditionally been studied using bulk foam/emulsion tests and single film platforms such as the Scheludko cell. Recently there has been a renewed interest in a third class of techniques that we term as single bubble/drop tests, which employ isolated whole bubbles and drops to probe the characteristics of foams and emulsions. Single bubble and drop techniques provide a convenient framework for investigating a number of important characteristics of foams and emulsions, including the rheology, stabilization mechanisms, and rupture dynamics. In this review we provide a comprehensive discussion of the various single bubble/drop platforms and the associated experimental measurement protocols including the construction of coalescence time distributions, visualization of the thin film profiles and characterization of the interfacial rheological properties. Subsequently, we summarize the recent developments in foam and emulsion science with a focus on the results obtained through single bubble/drop techniques. We conclude the review by presenting important venues for future research.


title: Planning and Operations of Mixed Fleets in Mobility-on-Demand Systems

abstract: Automated vehicles (AVs) are expected to be beneficial for Mobility-on-Demand (MoD), thanks to their ability of being globally coordinated. To facilitate the steady transition towards full autonomy, we consider the transition period of AV deployment, whereby an MoD system operates a mixed fleet of automated vehicles (AVs) and human-driven vehicles (HVs). In such systems, AVs are centrally coordinated by the operator, and the HVs might strategically respond to the coordination of AVs. We devise computationally tractable strategies to coordinate mixed fleets in MoD systems. Specifically, we model an MoD system with a mixed fleet using a Stackelberg framework where the MoD operator serves as the leader and human-driven vehicles serve as the followers. We develop two models: 1) a steady-state model to analyze the properties of the problem and determine the planning variables (e.g., compensations, prices, and the fleet size of AVs), and 2) a time-varying model to design a real-time coordination algorithm for AVs. The proposed models are validated using a case study inspired by real operational data of a MoD service in Singapore. Results show that the proposed algorithms can significantly improve system performance.


title: Subharmonic Solutions In Reversible Non-Autonomous Differential Equations

abstract: We study the existence of subharmonic solutions in the system $\ddot {u}(t)=f(t,u(t))$, where $u(t)\in\mathbb{R}^{k}$ and $f$ is an even and $p$-periodic function in time. Under some additional symmetry conditions on the function $f$, the problem of finding $mp$-periodic solutions can be reformulated in a functional space as a $\Gamma\times\mathbb{Z}_{2}\times D_{m}% $-equivariant equation, where the group $\Gamma\times\mathbb{Z}_{2}$ acts on the space $\mathbb{R}^{k}$ and $D_{m}$ acts on $u(t)$ by time-shifts and reflection. We apply Brouwer equivariant degree to prove the existence of an infinite number of subharmonic solutions for the function $f$ that satisfies additional hypothesis on linear behavior near zero and the Nagumo condition at infinity. We also discuss the bifurcation of subharmonic solutions when the system depends on an extra parameter.


title: The determinants of COVID-19 case fatality rate (CFR) in the Italian regions and provinces: an analysis of environmental, demographic, and healthcare factors

abstract: The Italian government has been one of the most responsive to COVID-19 emergency, through the adoption of quick and increasingly stringent measures to contain the outbreak. Despite this, Italy has suffered a huge human and social cost, especially in Lombardy. The aim of this paper is dual: i) first, to investigate the reasons of the case fatality rate (CFR) differences across Italian 20 regions and 107 provinces, using a multivariate OLS regression approach; and ii) second, to build a taxonomy of provinces with similar mortality risk of COVID-19, by using the Ward hierarchical agglomerative clustering method. I considered health system metrics, environmental pollution, climatic conditions, demographic variables, and three ad hoc indexes that represent the health system saturation. The results showed that overall health care efficiency, physician density, and average temperature helped to reduce the CFR. By the contrary, population aged 70 and above, car and firm density, level of air pollutants (NO2, O3, PM10, and PM2.5), relative average humidity, COVID-19 prevalence, and all three indexes of health system saturation were positively associated with the CFR. Population density, social vertical integration, and altitude were not statistically significant. In particular, the risk of dying increases with age, as 90 years old and above had a three-fold greater risk than the 80 to 89 years old and four-fold greater risk than 70 to 79 years old. Moreover, the cluster analysis showed that the highest mortality risk was concentrated in the north of the country, while the lowest risk was associated with southern provinces. Finally, since prevalence and health system saturation indexes played the most important role in explaining the CFR variability, a significant part of the latter may have been caused by the massive stress of the Italian health system.


title: Differentially Private Sketches for Jaccard Similarity Estimation

abstract: This paper describes two locally-differential private algorithms for releasing user vectors such that the Jaccard similarity between these vectors can be efficiently estimated. The basic building block is the well known MinHash method. To achieve a privacy-utility trade-off, MinHash is extended in two ways using variants of Generalized Randomized Response and the Laplace Mechanism. A theoretical analysis provides bounds on the absolute error and experiments show the utility-privacy trade-off on synthetic and real-world data. The paper ends with a critical discussion of related work.