title: The Lorentzian Lichnerowicz Conjecture for real-analytic, three-dimensional manifolds

abstract: Given a compact, three-dimensional, real-analytic Lorentzian manifold $(M,g)$, we prove that the identity component of the conformal group preserves a metric in the conformal class $[g]$, or $(M,g)$ is conformally flat.


title: Partially Supervised Named Entity Recognition via the Expected Entity Ratio Loss

abstract: We study learning named entity recognizers in the presence of missing entity annotations. We approach this setting as tagging with latent variables and propose a novel loss, the Expected Entity Ratio, to learn models in the presence of systematically missing tags. We show that our approach is both theoretically sound and empirically useful. Experimentally, we find that it meets or exceeds performance of strong and state-of-the-art baselines across a variety of languages, annotation scenarios, and amounts of labeled data. In particular, we find that it significantly outperforms the previous state-of-the-art methods from Mayhew et al. (2019) and Li et al. (2021) by +12.7 and +2.3 F1 score in a challenging setting with only 1,000 biased annotations, averaged across 7 datasets. We also show that, when combined with our approach, a novel sparse annotation scheme outperforms exhaustive annotation for modest annotation budgets.


title: On Kostant's weight $q$-multiplicity formula for $\mathfrak{sp}_6(\mathbb{C})$

abstract: Kostant's weight $q$-multiplicity formula is an alternating sum over a finite group known as the Weyl group, whose terms involve the $q$-analog of Kostant's partition function. The $q$-analog of the partition function is a polynomial-valued function defined by $\wp_q(\xi)=\sum_{i=0}^k c_i q^i$, where $c_i$ is the number of ways the weight $\xi$ can be written as a sum of exactly $i$ positive roots of a Lie algebra $\mathfrak{g}$. The evaluation of the $q$-multiplicity formula at $q = 1$ recovers the multiplicity of a weight in an irreducible highest weight representation of $\mathfrak{g}$. In this paper, we specialize to the Lie algebra $\mathfrak{sp}_6(\mathbb{C})$ and we provide a closed formula for the $q$-analog of Kostant's partition function, which extends recent results of Shahi, Refaghat, and Marefat. We also describe the supporting sets of the multiplicity formula (known as the Weyl alternation sets of $\mathfrak{sp}_6(\mathbb{C})$), and use these results to provide a closed formula for the $q$-multiplicity for any pair of dominant integral weights of $\mathfrak{sp}_6(\mathbb{C})$. Throughout this work, we provide code to facilitate these computations.


title: Strategic Exploration for Innovation

abstract: This paper introduces a framework to study innovation in a strategic setting, in which innovators allocate their resources between exploration and exploitation in continuous time. Exploration creates public knowledge, while exploitation delivers private benefits. Through the analysis of a class of Markov equilibria, we demonstrate that knowledge spillovers accelerate knowledge creation and expedite its availability, thereby encouraging innovators to increase exploration. The prospect of the ensuing superior long-term innovations further motivates exploration, giving rise to a positive feedback loop. This novel feedback loop can substantially mitigate the free-riding problem arising from knowledge spillovers.


title: General Formulas of the Structure Constants in the $\mathfrak{su}(N)$ Lie Algebra

abstract: We provide the analytic expressions of the totally symmetric and anti-symmetric structure constants in the $\mathfrak{su}(N)$ Lie algebra. The derivation is based on a relation linking the index of a generator to the indexes of its non-null elements. The closed formulas obtained to compute the values of the structure constants are simple expressions involving those indexes and can be analytically evaluated without any need of the expression of the generators. We hope that these expressions can be widely used for analytical and computational interest in Physics.

categories: math-ph math.MP physics.chem-ph quant-ph


title: Part 1: Spatially Dispersive Metasurfaces: Zero Thickness Surface Susceptibilities & Extended GSTCs

abstract: A simple method to describe spatially dispersive metasurfaces is proposed where the angle-dependent surface susceptibilities are explicitly used to formulate the zero thickness sheet model of practical metasurface structures. It is shown that if the surface susceptibilities of a given metasurface are expressed as a ratio of two polynomials of tangential spatial frequencies, $\boldsymbol{k_{||}}$ with complex coefficients, they can be conveniently expressed as spatial derivatives of the difference and average fields around the metasurface in the space domain, leading to extended forms of the standard Generalized Sheet Transition Conditions (GSTCs) accounting for the spatial dispersion. Using two simple examples of a short electric dipole and an all-dielectric cylindrical puck unit cells, which exhibit purely tangential surface susceptibilities and reciprocal/symmetric transmission and reflection characteristics, the proposed concept is numerically confirmed in 2D. A single Lorentzian has been found to describe the spatio-temporal frequency behavior of a short dipole unit cell, while a multi-Lorentzian description is developed to capture the complex multiple angular resonances of the dielectric puck. For both cases, the appropriate spatial boundary conditions are derived.


title: Confinement in bilayer graphene via intra- and inter-layer interactions

abstract: We consider confinement of Dirac fermions in $AB$-stacked bilayer graphene by inhomogeneous on-site interactions, (pseudo-)magnetic field or inter-layer interaction. Working within the framework of four-band approximation, we focus on the systems where the stationary equation is reducible into two stationary equations with $2\times2$ Dirac-type Hamiltonians and auxiliary interactions. We show that it is possible to find localized states by solving an effective Schr\"odinger equation with energy-dependent potential. We consider several scenarios where bilayer graphene is subject to inhomogneous (pseudo-)magnetic field, on-site interactions or inter-layer coupling. In explicit examples, we provide analytical solutions for the states localized by local fluctuations or periodicity defects of the interactions.


title: Deep learning analysis of polaritonic waves images

abstract: Deep learning (DL) is an emerging analysis tool across sciences and engineering. Encouraged by the successes of DL in revealing quantitative trends in massive imaging data, we applied this approach to nano-scale deeply sub-diffractional images of propagating polaritonic waves in complex materials. We developed a practical protocol for the rapid regression of images that quantifies the wavelength and the quality factor of polaritonic waves utilizing the convolutional neural network (CNN). Using simulated near-field images as training data, the CNN can be made to simultaneously extract polaritonic characteristics and materials parameters in a timescale that is at least three orders of magnitude faster than common fitting/processing procedures. The CNN-based analysis was validated by examining the experimental near-field images of charge-transfer plasmon polaritons at Graphene/{\alpha}-RuCl3 interfaces. Our work provides a general framework for extracting quantitative information from images generated with a variety of scanning probe methods.


title: Metall: A Persistent Memory Allocator For Data-Centric Analytics

abstract: Data analytics applications transform raw input data into analytics-specific data structures before performing analytics. Unfortunately, such data ingestion step is often more expensive than analytics. In addition, various types of NVRAM devices are already used in many HPC systems today. Such devices will be useful for storing and reusing data structures beyond a single process life cycle. We developed Metall, a persistent memory allocator built on top of the memory-mapped file mechanism. Metall enables applications to transparently allocate custom C++ data structures into various types of persistent memories. Metall incorporates a concise and high-performance memory management algorithm inspired by Supermalloc and the rich C++ interface developed by Boost.Interprocess library. On a dynamic graph construction workload, Metall achieved up to 11.7x and 48.3x performance improvements over Boost.Interprocess and memkind (PMEM kind), respectively. We also demonstrate Metall's high adaptability by integrating Metall into a graph processing framework, GraphBLAS Template Library. This study's outcomes indicate that Metall will be a strong tool for accelerating future large-scale data analytics by allowing applications to leverage persistent memory efficiently.


title: Multimodal Information Gain in Bayesian Design of Experiments

abstract: One of the well-known challenges in optimal experimental design is how to efficiently estimate the nested integrations of the expected information gain. The Gaussian approximation and associated importance sampling have been shown to be effective at reducing the numerical costs. However, they may fail due to the non-negligible biases and the numerical instabilities. A new approach is developed to compute the expected information gain, when the posterior distribution is multimodal - a situation previously ignored by the methods aiming at accelerating the nested numerical integrations. Specifically, the posterior distribution is approximated using a mixture distribution constructed by multiple runs of global search for the modes and weighted local Laplace approximations. Under any given probability of capturing all the modes, we provide an estimation of the number of runs of searches, which is dimension independent. It is shown that the novel global-local multimodal approach can be significantly more accurate and more efficient than the other existing approaches, especially when the number of modes is large. The methods can be applied to the designs of experiments with both calibrated and uncalibrated observation noises.


title: The VMC survey -- XLIII. The spatially resolved star formation history across the Large Magellanic Cloud

abstract: We derive the spatially-resolved star formation history (SFH) for a $96$ deg$^2$ area across the main body of the Large Magellanic Cloud (LMC), using the near-infrared photometry from the VISTA survey of the Magellanic Clouds (VMC). The data and analyses are characterised by a great degree of homogeneity and a low sensitivity to the interstellar extinction. 756 subregions of size $0.125$ deg$^2$ -- corresponding to projected sizes of about $296\times322\,\mathrm{pc}^{2}$ in the LMC -- are analysed. The resulting SFH maps, with typical resolution of $0.2$--$0.3$ dex in logarithm of age, reveal main features in the LMC disc at different ages: the patchy star formation at recent ages, the concentration of star formation on three spiral arms and on the Bar up to ages of $\sim\!1.6$ Gyr, and the wider and smoother distribution of older populations. The period of most intense star formation occurred roughly between 4 and 0.5 Gyr ago, at rates of $\sim\!0.3\,\mathrm{M}_{\odot}\mathrm{yr}^{-1}$. We compare young and old star formation rates with the observed numbers of RR Lyrae and Cepheids. We also derive a mean extinction and mean distance for every subregion, and the plane that best describes the spatial distribution of the mean distances. Our results cover an area about 50 per cent larger than the classical SFH maps derived from optical data by Harris & Zaritsky (2009). Main differences with respect to those maps are lower star formation rates at young ages, and a main peak of star formation being identified at ages slightly younger than $1$ Gyr.


title: Unifying Futures and Spot Market: Overbooking-Enabled Resource Trading in Mobile Edge Networks

abstract: Securing necessary resources for edge computing processes via effective resource trading becomes a critical technique in supporting computation-intensive mobile applications. Conventional onsite spot trading could facilitate this paradigm with proper incentives, which, however, incurs excessive decision-making latency/energy consumption, and further leads to underutilization of dynamic resources. Motivated by this, a hybrid market unifying futures and spot is proposed to facilitate resource trading among an edge server (seller) and multiple smart devices (buyers) by encouraging some buyers to sign a forward contract with seller in advance, while leaving the remaining buyers to compete for available resources with spot trading. Specifically, overbooking is adopted to achieve substantial utilization and profit advantages owing to dynamic resource demands. By integrating overbooking into futures market, mutually beneficial and risk-tolerable forward contracts with appropriate overbooking rate can be achieved relying on analyzing historical statistics associated with future resource demand and communication quality, which are determined by an alternative optimization-based negotiation scheme. Besides, spot trading problem is studied via considering uniform/differential pricing rules, for which two bilateral negotiation schemes are proposed by addressing both non-convex optimization and knapsack problems. Experimental results demonstrate that the proposed mechanism achieves mutually beneficial players' utilities, while outperforming baseline methods on critical indicators, e.g., decision-making latency, resource usage, etc.


title: Revisiting Empirical Bayes Methods and Applications to Special Types of Data

abstract: Empirical Bayes methods have been around for a long time and have a wide range of applications. These methods provide a way in which historical data can be aggregated to provide estimates of the posterior mean. This thesis revisits some of the empirical Bayesian methods and develops new applications. We first look at a linear empirical Bayes estimator and apply it on ranking and symbolic data. Next, we consider Tweedie's formula and show how it can be applied to analyze a microarray dataset. The application of the formula is simplified with the Pearson system of distributions. Saddlepoint approximations enable us to generalize several results in this direction. The results show that the proposed methods perform well in applications to real data sets.


title: Thermoelectric phenomenon in hollow blocks

abstract: The work presented in this article describes thermoelectric effect in hollow blocks for heat waste harvesting purposes. The study consists of developing a numerical model formed by a heat transfer equation coupled to thermoelectric effects equations to study thermoelectric generators(TEG) incorporated inside Lebanese hollow blocks through two simulations using finite difference scheme and using finite element scheme. Results showed a voltage of 5.85mV produced from a single 8.6 x 0.4 x 0.4 cm3 thermoelectric leg made of Bismuth Antimony Telluride for {\Delta}T=30K. A design with 3 TEGs incorporated inside a hollow block was tested and validated numerically using both methods, the main results obtained for {\Delta}T=30K, showed a voltage {\Delta}V=0.72V, a current I=0.06 A and a figure of merit ZT=0.55. The design was then optimized for economic purposes.


title: Patch Attack Invariance: How Sensitive are Patch Attacks to 3D Pose?

abstract: Perturbation-based attacks, while not physically realizable, have been the main emphasis of adversarial machine learning (ML) research. Patch-based attacks by contrast are physically realizable, yet most work has focused on 2D domain with recent forays into 3D. Characterizing the robustness properties of patch attacks and their invariance to 3D pose is important, yet not fully elucidated, and is the focus of this paper. To this end, several contributions are made here: A) we develop a new metric called mean Attack Success over Transformations (mAST) to evaluate patch attack robustness and invariance; and B), we systematically assess robustness of patch attacks to 3D position and orientation for various conditions; in particular, we conduct a sensitivity analysis which provides important qualitative insights into attack effectiveness as a function of the 3D pose of a patch relative to the camera (rotation, translation) and sets forth some properties for patch attack 3D invariance; and C), we draw novel qualitative conclusions including: 1) we demonstrate that for some 3D transformations, namely rotation and loom, increasing the training distribution support yields an increase in patch success over the full range at test time. 2) We provide new insights into the existence of a fundamental cutoff limit in patch attack effectiveness that depends on the extent of out-of-plane rotation angles. These findings should collectively guide future design of 3D patch attacks and defenses.


title: Simple Formulas for Output Interception Power Estimation of Uni-Traveling Carrier Photodiodes

abstract: Simple analytical expressions for estimation of second order output intercept point (OIP2) and third order output intercept point (OIP3) of surface normal uni-traveling carrier (UTC) and modified uni-traveling carrier (MUTC) photodiode (PD) are derived. These equations are valuable for estimation of OIP for high power (M)UTC-PDs during the design phase.


title: Brian\c{c}on-Skoda exponents and the maximal root of reduced Bernstein-Sato polynomials

abstract: For a holomorphic function $f$ on a complex manifold $X$, the Brian\c{c}on-Skoda exponent $e^{\rm BS}(f)$ is the smallest integer $k$ with $f^k\in(\partial f)$ (replacing $X$ with a neighborhood of $f^{-1}(0)$), where $(\partial f)$ denotes the Jacobian ideal of $f$. It is shown that $e^{\rm BS}(f)\le d_X$ $(:=\dim X)$ by Brian\c con-Skoda. We prove that $e^{\rm BS}(f)\le[d_X-2\widetilde{\alpha}_f]+1$ with $-\widetilde{\alpha}_f$ the maximal root of the reduced Bernstein-Sato polynomial $b_f(s)/(s+1)$, assuming the latter exists (shrinking $X$ if necessary). This implies for instance that $e^{\rm BS}(f)\le d_X-2$ in the case $f^{-1}(0)$ has only rational singularities, that is, if $\widetilde{\alpha}_f>1$.


title: Better GPU Hash Tables

abstract: We revisit the problem of building static hash tables on the GPU and design and build three bucketed hash tables that use different probing schemes. Our implementations are lock-free and offer efficient memory access patterns; thus, only the probing scheme is the factor affecting the performance of the hash table's different operations. Our results show that a bucketed cuckoo hash table that uses three hash functions (BCHT) outperforms alternative methods that use power-of-two choices, iceberg hashing, and a cuckoo hash table that uses a bucket size one. At high load factors as high as 0.99, BCHT enjoys an average probe count of 1.43 during insertion. Using three hash functions only, positive and negative queries require at most 1.39 and 2.8 average probes per key, respectively.


title: Fast predictions of lattice energies by continuous isometry invariants of crystal structures

abstract: Crystal Structure Prediction (CSP) aims to discover solid crystalline materials by optimizing periodic arrangements of atoms, ions or molecules. CSP takes weeks of supercomputer time because of slow energy minimizations for millions of simulated crystals. The lattice energy is a key physical property, which determines thermodynamic stability of a crystal but has no simple analytic expression. Past machine learning approaches to predict the lattice energy used slow crystal descriptors depending on manually chosen parameters. The new area of Periodic Geometry offers much faster isometry invariants that are also continuous under perturbations of atoms. Our experiments on simulated crystals confirm that a small distance between the new invariants guarantees a small difference of energies. We compare several kernel methods for invariant-based predictions of energy and achieve the mean absolute error of less than 5kJ/mole or 0.05eV/atom on a dataset of 5679 crystals.


title: Flexibility and regularity of the hydration structure of a sodium ion. Nonempirical insight

abstract: The stationary nonempirical simulations of Na+(H2O)n clusters with n in a range of 28 to 51 carried out at the density functional level with a hybrid B3LYP functional and the Born-Oppenheimer molecular dynamics modeling of the size selected clusters reveal the interrelated structural and energetic peculiarities of sodium hydration structures. Surface, bulk, and transient configurations of the clusters are distinguished with the different location of the sodium nucleus (close to either the spatial center of the structure or one of its side faces) and its consistently changing coordination number (which typically equals five or six). The <rNaO> mean Na-O distances for the first-shell water molecules are found to depend both on the spatial character of the structure and the local coordination of sodium. The <rNaO> values are compared to different experimental estimates, and the virtual discrepancy of the latter is explained based on the results of the cluster simulations. Different coordination neighborhoods of sodium are predicted depending on its local fraction in the actual specimens.

categories: physics.chem-ph cond-mat.mtrl-sci physics.atm-clus physics.comp-ph


title: Exactly solving the Kitaev chain and generating Majorana-zero-modes out of noisy qubits

abstract: Majorana-zero-modes (MZMs) were predicted to exist as edge states of a physical system called the Kitaev chain. MZMs should host particles that are their own antiparticles and could be used as a basis for a qubit which is robust-to-noise. However, all attempts to prove their existence gave inconclusive results. Here, the Kitaev chain is exactly solved with a quantum computing methodology and properties of MZMs are probed by generating eigenstates of the Kitev Hamiltonian on 3 noisy qubits of a publicly available quantum computer. After an ontological elaboration I show that two eigenstates of the Kitaev Hamiltonian exhibit eight signatures attributed to MZMs. The results presented here are a most comprehensive set of validations of MZMs ever conducted in an actual physical system. Furthermore, the findings of this manuscript are easily reproducible for any user of publicly available quantum computers, solving another important problem of research with MZMs-the result reproducibility crisis.


title: Multihop RIS-Assisted FSO-RF System Over Double Generalized Gamma Fading

abstract: Reconfigurable intelligent surface (RIS) is a promising technology to avoid signal blockage by creating virtual line-of-sight (LOS) connectivity for free-space optical (FSO) and radio frequency (RF) wireless systems. This paper considers a mixed FSO-RF system by employing multiple RISs in both the links for multihop transmissions to extend the communication range. We develop probability density function (PDF) and cumulative density function (CDF) of the signal-to-noise ratio (SNR) for the cascaded channels by considering double generalized gamma (dGG) turbulence with pointing errors for the FSO link and the dGG distribution to model the signal fading for the RF. We derive exact closed-form expressions of the outage probability, average bit-error-rate (BER), and ergodic capacity using the decode-and-forward (DF) relaying for the mixed system. We also present asymptotic analysis on the performance in the high SNR regime depicting the impact of channel parameters on the diversity order of the system. We use computer simulations to demonstrate the effect of system and channel parameters on the RIS-aided multihop transmissions.


title: A Novel Near-field Photonic Thermal Diode with hBN and InSb

abstract: Similar to the diode in electronics, a thermal diode is a two-terminal device that allows heat to transfer easier in one direction (forward bias) than in the opposite direction (reverse bias). Unlike conductive and convective thermal diodes, a photonic thermal diode operates in a contactless mode and may afford a large operating temperature range. In this work, a novel near-field photonic thermal diode with hexagonal boron nitride (hBN) and indium antimonide (InSb) is proposed and theoretically demonstrated. The temperature dependence of the interband absorption of InSb is used to couple (or decouple) with the hyperbolic phonon polaritons in hBN. The numerical analysis predicts a rectification ratio greater than 17 for a 10 nm vacuum gap when operating at an average temperature of 300 K and a temperature difference of 200 K. The calculated rectification ratio exceeds 35 with higher average temperatures and larger temperature differences. The mechanism proposed here for achieving photonic thermal rectification provides a new way of controlling radiative heat transfer.


title: Active fault tolerant control for twin wind turbine subject to asymmetric fault

abstract: This paper addresses the problem of control of a twin wind turbine which is subject to an electrical fault affecting only one stator phase of one turbine. An active fault tolerant control is proposed. The performance and robustness of the proposed control, comparing to a passive fault tolerant one developed in the literature, are shown through numerical simulations.


title: Rotational symmetry breaking potential for two-dimensional magnets

abstract: Here we present a new perspective to the breakdown of ferromagnetic order in two-dimensional spin-lattice models employing the rotation of the underlying lattice. Using an Ising spin system on a square lattice as a prototype, we demonstrate that an additional low-symmetry interaction may lead to the absence of the truly long-range order and forms aperiodic structure, such as magnetic stripes. Employing annealing and entropic Monte Carlo simulations, we show that our model allows tuning between different phases, magnetically ordered as well as more exotic nonmagnetic phases such as Ising-nematic by changing only one control parameter, which is responsible for the arising of magnetic frustration. In addition, our methodology of considering the coupling between the magnetic structure and the host material can be extended to the study of any type of spin-exchange model in two dimensions and has many potential interesting ramifications and applications.


title: An analytical approach to the Mean--Return-Time~Phase of isotropic stochastic oscillators

abstract: One notion of phase for stochastic oscillators is based on the mean return-time (MRT): a set of points represents a certain phase if the mean time to return from any point in this set to this set after one rotation is equal to the mean rotation period of the oscillator (irrespective of the starting point). For this so far only algorithmically defined phase, we derive here analytical expressions for the important class of isotropic stochastic oscillators. This allows us to evaluate cases from the literature explicitly and to study the behavior of the MRT phase in the limits of strong noise. We also use the same formalism to show that lines of constant return time variance (instead of constant mean return time) can be defined, and that they in general differ from the MRT-isochrons.


title: Large genus bounds for the distribution of triangulated surfaces in moduli space

abstract: Triangulated surfaces are compact Riemann surfaces equipped with a conformal triangulation by equilateral triangles. In 2004, Brooks and Makover asked how triangulated surfaces are distributed in the moduli space of Riemann surfaces as the genus tends to infinity. Mirzakhani raised this question in her 2010 ICM address. We show that in the large genus case, triangulated surfaces are well distributed in moduli space in a fairly strong sense. We do this by proving upper and lower bounds for the number of triangulated surfaces lying in a Teichm\"uller ball in moduli space. In particular, we show that the number of triangulated surfaces lying in a Teichm\"uller unit ball is at most exponential in the number of triangles, independent of the genus.


title: Doubly heavy tetraquarks in an extended chromomagnetic model

abstract: Using an extended chromomagnetic model, we perform a systematic study of the masses of the doubly heavy tetraquarks. We find that the ground states of the doubly heavy tetraquarks are dominated by color-triplet $\ket{(qq)^{\bar{3}_{c}}(\bar{Q}\bar{Q})^{3_{c}}}$ configuration, which is opposite to that of the fully heavy tetraquarks. The combined results suggest that the color-triplet configuration becomes more important when the mass difference between the quarks and antiquarks increases. We find three stable states which lie below the thresholds of two pseudoscalar mesons. They are the $IJ^{P}=01^{+}$ $nn\bar{b}\bar{b}$ tetraquark, the $IJ^{P}=00^{+}$ $nn\bar{c}\bar{b}$ tetraquark and the $J^{P}=1^{+}$ $ns\bar{b}\bar{b}$ tetraquark.


title: A Physics Informed Neural Network Approach to Solution and Identification of Biharmonic Equations of Elasticity

abstract: We explore an application of the Physics Informed Neural Networks (PINNs) in conjunction with Airy stress functions and Fourier series to find optimal solutions to a few reference biharmonic problems of elasticity and elastic plate theory. Biharmonic relations are fourth-order partial differential equations (PDEs) that are challenging to solve using classical numerical methods, and have not been addressed using PINNs. Our work highlights a novel application of classical analytical methods to guide the construction of efficient neural networks with the minimal number of parameters that are very accurate and fast to evaluate. In particular, we find that enriching feature space using Airy stress functions can significantly improve the accuracy of PINN solutions for biharmonic PDEs.


title: Physics-Informed Machine Learning of the Lagrangian Dynamics of Velocity Gradient Tensor

abstract: Reduced models describing the Lagrangian dynamics of the Velocity Gradient Tensor (VGT) in Homogeneous Isotropic Turbulence (HIT) are developed under the Physics-Informed Machine Learning (PIML) framework. We consider VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and sub-filter contributions, which is accomplished by constructing them using the integrity bases and invariants of VGT. The developed models can be expressed using the extended Tensor Basis Neural Network (TBNN). Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number Direct Numerical Simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the "ground-truth" DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into TBNN.


title: Tensor distributions with covariance tensor or correlation tensor

abstract: In this article, we define the matricization of a tensor and we present some properties of the matricization. After that, we define the determinant of a tensor and we present some properties of the determinant. We define the covariance tensor and we present some properties of the covariance tensor. In a similar way, we define the correlation tensor. We define the tensor normal distribution. In a similar way, we define the tensor elliptical distributions. We prove the equivalence of the tensor elliptical distribution representations.


title: Drifting Electrons: Nonreciprocal Plasmonics and Thermal Photonics

abstract: Light propagates symmetrically in opposite directions in most materials and structures. This fact -- a consequence of the Lorentz reciprocity principle -- has tremendous implications for science and technology across the electromagnetic spectrum. Here, we investigate an emerging approach to break reciprocity that does not rely on magneto-optical effects or spacetime modulations, but is instead based on biasing a plasmonic material with a direct electric current. Using a 3D Green function formalism and microscopic considerations, we elucidate the propagation properties of surface plasmon-polaritons (SPPs) supported by a generic nonreciprocal platform of this type, revealing some previously overlooked, anomalous, wave-propagation effects. We show that SPPs can propagate in the form of steerable, slow-light, unidirectional beams associated with inflexion points in the modal dispersion. We also clarify the impact of dissipation (due to collisions and Landau damping) on nonreciprocal effects and shed light on the connections between inflexion points, exceptional points at band-edges, and complex modal transitions in leaky-wave structures. We then apply these concepts to the important area of thermal photonics, and provide the first theoretical demonstration of drift-induced nonreciprocal near-field radiative heat transfer between two planar bodies. Our findings may open new opportunities toward the development of nonreciprocal magnet-free devices that combine the benefits of plasmonics and nonreciprocal photonics for wave-guiding and energy applications.


title: Robust Hierarchical Clustering for Directed Networks: An Axiomatic Approach

abstract: We provide a complete taxonomic characterization of robust hierarchical clustering methods for directed networks following an axiomatic approach. We begin by introducing three practical properties associated with the notion of robustness in hierarchical clustering: linear scale preservation, stability, and excisiveness. Linear scale preservation enforces imperviousness to change in units of measure whereas stability ensures that a bounded perturbation in the input network entails a bounded perturbation in the clustering output. Excisiveness refers to the local consistency of the clustering outcome. Algorithmically, excisiveness implies that we can reduce computational complexity by only clustering a subset of our data while theoretically guaranteeing that the same hierarchical outcome would be observed when clustering the whole dataset. In parallel to these three properties, we introduce the concept of representability, a generative model for describing clustering methods through the specification of their action on a collection of networks. Our main result is to leverage this generative model to give a precise characterization of all robust -- i.e., excisive, linear scale preserving, and stable -- hierarchical clustering methods for directed networks. We also address the implementation of our methods and describe an application to real data.


title: Environment-assisted strong coupling regime

abstract: Strong coupling regime takes place in open hybrid systems consisting of two or more physical subsystems when the coupling strength between subsystems exceeds the relaxation rate. The relaxation arises due to the interaction of the system with environment. For this reason, it is usually believed that the enhancement of the interaction with environment inevitably leads to a transition of the system from the strong to weak coupling regime. In this paper, we refute this common opinion. We demonstrate the interaction of the coupled system with environment induces an additional coupling between the subsystems that contribute to retention the system in the strong coupling regime. We show that the environmental-induced coupling strength is proportional to the product of the Rabi coupling strength by the gradient of the density of states of the reservoir. There is a critical Rabi coupling strength above which the environmental-induced coupling ensures that the system remains in the strong coupling regime at any relaxation rate. In this case, the strong coupling regime takes place even when the relaxation rate is significantly above the Rabi coupling strength between the subsystems. The critical coupling depends on the gradient of the reservoir density of states. We demonstrate that managing this gradient can serve as an additional tool to control the properties of the coupled systems.


title: BloomNet: A Robust Transformer based model for Bloom's Learning Outcome Classification

abstract: Bloom taxonomy is a common paradigm for categorizing educational learning objectives into three learning levels: cognitive, affective, and psychomotor. For the optimization of educational programs, it is crucial to design course learning outcomes (CLOs) according to the different cognitive levels of Bloom Taxonomy. Usually, administrators of the institutions manually complete the tedious work of mapping CLOs and examination questions to Bloom taxonomy levels. To address this issue, we propose a transformer-based model named BloomNet that captures linguistic as well semantic information to classify the course learning outcomes (CLOs). We compare BloomNet with a diverse set of basic as well as strong baselines and we observe that our model performs better than all the experimented baselines. Further, we also test the generalization capability of BloomNet by evaluating it on different distributions which our model does not encounter during training and we observe that our model is less susceptible to distribution shift compared to the other considered models. We support our findings by performing extensive result analysis. In ablation study we observe that on explicitly encapsulating the linguistic information along with semantic information improves the model on IID (independent and identically distributed) performance as well as OOD (out-of-distribution) generalization capability.


title: Mechanism for Switchability in Electron-Doped Ferroelectric Interfaces

abstract: With the recent experimental verification that ferroelectric lattice distortions survive in the metallic phase of some materials, there is a desire to create devices that are both switchable and take advantage of the novel functionalities afforded by polar interfaces. In this work, we explore a simple model for such an interface and demonstrate a mechanism by which a metallic ferroelectric substrate may be switched by a bias voltage. This finding is in contrast to the reasonable expectation that hysteresis is prevented by screening of external fields in ferroelectric metals. Instead, the electron gas binds to polarization gradients to form a compensated state. Uncompensated electrons, which may screen external fields, are generated either when the electron density exceeds the ferroelectric polarization or when the bias voltage exceeds a "spillover" threshold. We propose that switchable thin films may be optimized by choosing an electron density that is slightly less than the lattice polarization. In addition to the high-polarization states, we find that thin metallic films also have a low-polarization state with average polarization near zero. Unlike in insulating films, where the polarization is small everywhere in this state, the low-polarization state in the metallic films consists of two head-to-head domains of opposite polarization. This domain formation is enabled by the screening of depolarizing fields by the electron gas.


title: Molecular Beam Epitaxy of a 2D material nearly lattice matched to a 3D substrate: $NiTe_{2}$ on $GaAs$

abstract: The lattice mismatch between interesting 2D materials and commonly available 3D substrates is one of the obstacles in the epitaxial growth of monolithic 2D/3D heterostructures, but a number of 2D materials have not yet been considered for epitaxy. Here we present the first molecular beam epitaxy growth of NiTe$_{2}$ 2D transition metal dichalcogenide. Importantly, the growth is realized on a nearly lattice matched GaAs(111)B substrate. Structural properties of the grown layers are investigated by electron diffraction, X-ray diffraction, and scanning tunnelling microscopy. Surface coverage and atomic scale order is evidenced by images obtained with atomic force, scanning electron, and transmission electron microscopy. Basic transport properties were measured confirming that NiTe$_{2}$ layers are metallic, with the Hall concentration of $10^{20}$cm$^{-3}$ to $10^{23}$cm$^{-3}$, depending on the growth conditions.


title: Circuit complexity in proca theory

abstract: In this paper, we study circuit complexity in Proca theory with Nielsen's approach and Fubini-Study (FS) metric approach. We place the fields on a lattice to gain a regularized theory, and obtain the ground state by adopting proper coordinates. We calculate complexities of the ground and thermofield double (TFD) states with Nielsen's approach, complexity of the TFD state is found to grows like a logarithmic function. We quantize the Proca fields and give the approximate ground state and TFD state by acting unitary circuit operators on the associated reference states. The circuit lengths are calculated with FS metric, the minimal lengths are given according to the associated geometric spaces. The complexity of TFD state is found to grows linearly with time.


title: Who's Waldo? Linking People Across Text and Images

abstract: We present a task and benchmark dataset for person-centric visual grounding, the problem of linking between people named in a caption and people pictured in an image. In contrast to prior work in visual grounding, which is predominantly object-based, our new task masks out the names of people in captions in order to encourage methods trained on such image-caption pairs to focus on contextual cues (such as rich interactions between multiple people), rather than learning associations between names and appearances. To facilitate this task, we introduce a new dataset, Who's Waldo, mined automatically from image-caption data on Wikimedia Commons. We propose a Transformer-based method that outperforms several strong baselines on this task, and are releasing our data to the research community to spur work on contextual models that consider both vision and language.


title: Coexistence of extended flat band and Kekul\'e order in Li-intercalated graphene

abstract: Doping graphene near the 1/4 filling to shift the extended flat band and van Hove singularity below E$_F$ has been highly desirable. Here we report the experimental observation of an extended flat band below E$_F$ in Li-intercalated graphene. Strong electron-phonon interaction is clearly identified by notable kinks in the band dispersion. Moreover, the evolution of the band structure upon Li intercalation shows that the extended flat band and the Kekul\'e order emerge simultaneously. Our work provides opportunities for investigating flat band related instabilities and its interplay with the Kekul\'e order


title: X-ray spectra, light-curves and SEDs of blazars frequently observed by Swift

abstract: Blazars research is one of the hot topics of contemporary extra-galactic astrophysics. That is because these sources are the most abundant type of extra-galactic gamma-ray sources and are suspected to play a central role in multi-messenger astrophysics. We have used swift_xrtproc, a tool to carry out an accurate spectral and photometric analysis of the Swift-XRT data of all blazars observed by Swift at least 50 times between December 2004 and the end of 2020. We present a database of X-ray spectra, best-fit parameter values, count-rates and flux estimations in several energy bands of over 31,000 X-ray observations and single snapshots of 65 blazars. The results of the X-ray analysis have been combined with other multi-frequency archival data to assemble the broad-band Spectral Energy Distributions (SEDs) and the long-term light-curves of all sources in the sample. Our study shows that large X-ray luminosity variability on different timescales is present in all objects. Spectral changes are also frequently observed with a "harder-when-brighter" or "softer-when-brighter" behavior depending on the SED type of the blazars. The peak energy of the synchrotron component nu_peak in the SED of HBL blazars, estimated from the log-parabolic shape of their X-ray spectra, also exhibits very large changes in the same source, spanning a range of over two orders of magnitude in Mrk421 and Mrk501, the objects with the best data sets in our sample.


title: NeuraCrypt is not private

abstract: NeuraCrypt (Yara et al. arXiv 2021) is an algorithm that converts a sensitive dataset to an encoded dataset so that (1) it is still possible to train machine learning models on the encoded data, but (2) an adversary who has access only to the encoded dataset can not learn much about the original sensitive dataset. We break NeuraCrypt privacy claims, by perfectly solving the authors' public challenge, and by showing that NeuraCrypt does not satisfy the formal privacy definitions posed in the original paper. Our attack consists of a series of boosting steps that, coupled with various design flaws, turns a 1% attack advantage into a 100% complete break of the scheme.


title: Resummation of the Sivers asymmetry in heavy flavor dijet production at the Electron-Ion Collider

abstract: We review our transverse momentum dependent factorization and resummation formalism for heavy flavor dijet production at the EIC. In this formalism, we have calculated the heavy flavor mass corrections in the collinear-soft and jet functions, and in the resummed expression for the cross section. By establishing this formalism, we then study the effects of the mass corrections by providing predictions at the EIC for the massive case and for the case where the mass is neglected. We find that the heavy flavor mass effects can give sizable corrections to the predicted asymmetry.


title: On the Opportunities and Risks of Foundation Models

abstract: AI is undergoing a paradigm shift with the rise of models (e.g., BERT, DALL-E, GPT-3) that are trained on broad data at scale and are adaptable to a wide range of downstream tasks. We call these models foundation models to underscore their critically central yet incomplete character. This report provides a thorough account of the opportunities and risks of foundation models, ranging from their capabilities (e.g., language, vision, robotics, reasoning, human interaction) and technical principles(e.g., model architectures, training procedures, data, systems, security, evaluation, theory) to their applications (e.g., law, healthcare, education) and societal impact (e.g., inequity, misuse, economic and environmental impact, legal and ethical considerations). Though foundation models are based on standard deep learning and transfer learning, their scale results in new emergent capabilities,and their effectiveness across so many tasks incentivizes homogenization. Homogenization provides powerful leverage but demands caution, as the defects of the foundation model are inherited by all the adapted models downstream. Despite the impending widespread deployment of foundation models, we currently lack a clear understanding of how they work, when they fail, and what they are even capable of due to their emergent properties. To tackle these questions, we believe much of the critical research on foundation models will require deep interdisciplinary collaboration commensurate with their fundamentally sociotechnical nature.


title: APReL: A Library for Active Preference-based Reward Learning Algorithms

abstract: Reward learning is a fundamental problem in human-robot interaction to have robots that operate in alignment with what their human user wants. Many preference-based learning algorithms and active querying techniques have been proposed as a solution to this problem. In this paper, we present APReL, a library for active preference-based reward learning algorithms, which enable researchers and practitioners to experiment with the existing techniques and easily develop their own algorithms for various modules of the problem. APReL is available at https://github.com/Stanford-ILIAD/APReL.


title: Reassessing the Limitations of CNN Methods for Camera Pose Regression

abstract: In this paper, we address the problem of camera pose estimation in outdoor and indoor scenarios. In comparison to the currently top-performing methods that rely on 2D to 3D matching, we propose a model that can directly regress the camera pose from images with significantly higher accuracy than existing methods of the same class. We first analyse why regression methods are still behind the state-of-the-art, and we bridge the performance gap with our new approach. Specifically, we propose a way to overcome the biased training data by a novel training technique, which generates poses guided by a probability distribution from the training set for synthesising new training views. Lastly, we evaluate our approach on two widely used benchmarks and show that it achieves significantly improved performance compared to prior regression-based methods, retrieval techniques as well as 3D pipelines with local feature matching.


title: $K$-core analysis of shear-thickening suspensions

abstract: Shear thickening of suspensions is studied by discrete-particle simulation, accounting for hydrodynamic, repulsive, and contact forces. The contact forces, including friction, are activated when the imposed shear stress $\sigma$ is able to overcome the repulsive force. The simulation method captures strong continuous and discontinuous shear thickening (CST and DST) in the range of solid volume fraction $0.54 \le \phi\le 0.56$ studied here. This work presents characteristics of the contact force network developed in the suspension under shear. The number of frictional contacts per particle $Z$ is shown to have a one-to-one relationship with the suspension stress, and the conditions for simple percolation of frictional contacts are found to deviate strongly from those of a random network model. The stress is shown to have important correlations with topological invariant metrics of the contact network known as $k$-cores; the $k$-cores are maximal subgraphs (`clusters') in which all member particles have $k$ or more frictional contacts to other members of the same subgraph. Only $k\le 3$ is found in this work at solid volume fractions $\phi \le 0.56$. Distinct relationships between the suspension rheology and the $k$-cores are found. One is that the stress susceptibility, defined as $\partial \sigma/\partial \dot{\gamma}$ where $\dot{\gamma}$ is the shear rate, is found to peak at the condition of onset of the $3$-core, regardless of whether the system exhibits CST or DST. A second is that the stress per particle within cores of different $k$ increases sharply with increase of $k$ at the onset of DST; in CST, the difference is mild.


title: Attractor mechanisms of moduli spaces of Calabi-Yau 3-folds

abstract: We investigate the complex and K\"ahler attractor mechanisms of moduli spaces of Calabi-Yau 3-folds. The complex attractor mechanism was previously studied by Ferrara-Kallosh-Strominger, Moore and others in string theory. It is concerned with the minimizing problems of the normalized central charges of 3-cycles and defines a new interesting class of Calabi-Yau 3-folds called, the complex attractor varieties. In light of mirror symmetry, we introduce the K\"ahler attractor mechanism and define the K\"ahler attractor varieties. The complex and K\"ahler attractor varieties are expected to possess very rich structures, in particular certain complex and K\"ahler rigidities.


title: Variational Inference at Glacier Scale

abstract: We characterize the complete joint posterior distribution over spatially-varying basal traction and and ice softness parameters of an ice sheet model from observations of surface speed by using stochastic variational inference combined with natural gradient descent to find an approximating variational distribution. By placing a Gaussian process prior over the parameters and casting the problem in terms of eigenfunctions of a kernel, we gain substantial control over prior assumptions on parameter smoothness and length scale, while also rendering the inference tractable. In a synthetic example, we find that this method recovers known parameters and accounts for mutual indeterminacy, both of which can influence observed surface speed. In an application to Helheim Glacier in Southeast Greenland, we show that our method scales to glacier-sized problems. We find that posterior uncertainty in regions of slow flow is high regardless of the choice of observational noise model.


title: A model problem for multiplicative chaos in number theory

abstract: Resolving a conjecture of Helson, Harper recently established that partial sums of random multiplicative functions typically exhibit more than square-root cancellation. Harper's work gives an example of a problem in number theory that is closely linked to ideas in probability theory connected with multiplicative chaos; another such closely related problem is the Fyodorov-Hiary-Keating conjecture on the maximum size of the Riemann zeta function in intervals of bounded length on the critical line. In this paper we consider a problem that might be thought of as a simplified function field version of Helson's conjecture. We develop and simplify the ideas of Harper in this context, with the hope that the simplified proof would be of use to readers seeking a gentle entry-point to this fascinating area.


title: The Integrated Probabilistic Data Association Filter Adapted to Lie Groups

abstract: The Integrated Probabilistic Data Association Filter (IPDAF) is a target tracking algorithm based on the Probabilistic Data Association Filter that calculates a statistical measure that indicates if an estimated representation of the target properly represents the target or is generated from non-target-originated measurements. The main contribution of this paper is to adapt the IPDAF to constant velocity target models that evolve on connected, unimodular Lie groups, and where the measurements are also defined on a Lie group. We present an example where the methods developed in the paper are applied to the problem of tracking a ground vehicle on the special Euclidean group SE(2).


title: The tree world

abstract: Tree amplitudes of any gauge theory and gravity can be factorized into primitive three-particle amplitudes by the BCFW recursion relations. We show that the amplitudes at any perturbation order are given by tree amplitudes with additional pairs of particles in the forward limit. As an example we compute the four-gluon all-plus helicity amplitude at next-to-leading order.


title: Relativistic Mechanics Theory for Electrons that Exhibits Spin, Zitterbewegung, Superposition and Produces Dirac's Wave Equation

abstract: A neo-classical relativistic mechanics theory is presented where the spin of an electron is a natural part of its space-time path as a point particle. The fourth-order equation of motion corresponds to the same Lagrangian function in proper time as in special relativity except for an additional spin energy term. The total motion can be decomposed into a sum of a local spin motion about a point and a global motion of this point, called the spin center. The global motion is sub-luminal and obeys Newton's second law in proper time, the time for a clock fixed at the spin center, while the total motion occurs at the speed of light c, consistent with the eigenvalues of Dirac's velocity operators having magnitude c. The local spin motion corresponds to Schr\"odinger's zitterbewegung and is a perpetual motion, which for a free electron has a circular path in the spin-center frame. In an electro-magnetic field, this spin motion generates magnetic and electric dipole energies through the Lorentz force on the electron's point charge. The corresponding electric dipole energy is consistent with the spin-orbit coupling term in the corrected Pauli non-relativistic Hamiltonian but the magnetic dipole energy is one half of that in Dirac's theory. By defining a spin tensor as the angular momentum of the electron's total motion about its spin center, the fundamental equations of motion can be re-written in an identical form to those of the Barut-Zanghi electron theory. These equations of motion can then be expressed using operators applied to a state function of proper time satisfying a Dirac-Schr\"odinger spinor equation. The operators produce dynamic variables without any probability implications. For the free electron, the state function satisfies Dirac's relativistic wave equation when the Lorentz transformation is applied to express proper time in terms of an observer's space-time coordinates.


title: The angular correlation between the fission fragment intrinsic spins

abstract: The generation of fission fragments (FF) spins and of their relative orbital angular momentum has been debated for more than six decades and no consensus has been yet achieved so far. (The published abstract, being too long, can be found in the manuscript.)


title: The Influence of Approximation in Generalized Uncertainty Principle on Black Hole Evaporation

abstract: The generalized uncertainty principle is often used to modify various thermodynamics systems by regarding the greater-than-equal relation as an approximate relation. We give a method to improve this approximation and compare the differences between the original and improved methods during the evaporation of black hole from two aspects of positive and negative parameters. Finally, we prove the rationality of the improved method and give some guiding opinions.


title: Information Exchange in Causally Nonseparable Processes

abstract: For a system of two parties, the process matrix framework predicts the existence of causally nonseparable structures. We characterize the information exchanged, showing that the total entropy of the two parties acts as a measure for the nonseparability.


title: Disorder induced 2D metal-insulator transition in moir\'e transition metal dichalcogenide multilayers

abstract: We develop a minimal theory for the recently observed metal-insulator transition (MIT) in two-dimensional (2D) moir\'e multilayer transition metal dichalcogenides (mTMD) using Coulomb disorder in the environment as the underlying mechanism. In particular, carrier scattering by random charged impurities leads to an effective 2D MIT approximately controlled by the Ioffe-Regel criterion, which is qualitatively consistent with the experiments. We find the necessary disorder to be around $5$-$10\times10^{10}$cm$^{-2}$ random charged impurities in order to quantitatively explain much, but not all, of the observed MIT phenomenology as reported by two different experimental groups. Our estimate is consistent with the known disorder content in TMDs.


title: Classical approaches to prethermal discrete time crystals in one, two, and three dimensions

abstract: We provide a comprehensive account of prethermal discrete time crystals within classical Hamiltonian dynamics, complementing and extending our recent work [Phys. Rev. Lett. 127, 140602 (2021)]. Considering power-law interacting spins on one-, two-, and three-dimensional hypercubic lattices, we investigate the interplay between dimensionality and interaction range in the stabilization of these non-equilibrium phases of matter that break the discrete time-translational symmetry of a periodic drive.

categories: quant-ph cond-mat.dis-nn cond-mat.stat-mech cond-mat.str-el nlin.CD


title: Driven-dissipative criticality within the discrete truncated Wigner approximation

abstract: We present an approach to the numerical simulation of open quantum many-body systems based on the semiclassical framework of the discrete truncated Wigner approximation. We establish a quantum jump formalism to integrate the quantum master equation describing the dynamics of the system, which we find to be exact in both the noninteracting limit and the limit where the system is described by classical rate equations. We apply our method to simulation of the paradigmatic dissipative Ising model, where we are able to capture the critical fluctuations of the system beyond the level of mean-field theory.


title: The time traveler's guide to the quantization of zero modes

abstract: We study the relationship between the quantization of a massless scalar field on the two-dimensional Einstein cylinder and in a spacetime with a time machine. We find that the latter picks out a unique prescription for the state of the zero mode in the Einstein cylinder. We show how this choice arises from the computation of the vacuum Wightman function and the vacuum renormalized stress-energy tensor in the time-machine geometry. Finally, we relate the previously proposed regularization of the zero mode state as a squeezed state with the time-machine warp parameter, thus demonstrating that the quantization in the latter regularizes the quantization in an Einstein cylinder.


title: Superfluid Effective Field Theory for Dark Matter Direct Detection

abstract: We develop an effective field theory (EFT) framework for superfluid ${}^4$He to model the interactions among quasiparticles, helium atoms and probe particles. Our effective field theory approach brings together symmetry arguments and power-counting and matches to classical fluid dynamics. We then present the decay and scattering rates for the relevant processes involving quasiparticles and helium atoms. The presented EFT framework and results can be used to understand the dynamics of thermalization in the superfluid, and can be further applied to sub-GeV dark matter direct detection with superfluid ${}^4$He.


title: On the single-event-based identification of primordial black hole mergers at cosmological distances

abstract: The existence of primordial black holes (PBHs), which may form from the collapse of matter overdensities shortly after the Big Bang, is still under debate. Among the potential signatures of PBHs are gravitational waves (GWs) emitted from binary black hole (BBH) mergers at redshifts $z\gtrsim 30$, where the formation of astrophysical black holes is unlikely. Future ground-based GW detectors, Cosmic Explorer and Einstein Telescope, will be able to observe equal-mass BBH mergers with total mass of $\mathcal{O}(10-100)~M_{\odot}$ at such distances. In this work, we investigate whether the redshift measurement of a single BBH source can be precise enough to establish its primordial origin. We simulate BBHs of different masses, mass ratios and orbital orientations. We show that for BBHs with total masses between $20~M_{\odot}$ and $40~M_{\odot}$ merging at $z \geq 40$ one can infer $z>30$ at up to 97\% credibility, with a network of one Einstein Telescope, one 40-km Cosmic Explorer in the US and one 20-km Cosmic Explorer in Australia. A smaller network made of one Einstein Telescope and one 40-km Cosmic Explorer in the US measures $z>30$ at larger than 90\% credibility for roughly half of the sources than the larger network. We then assess the dependence of this result on the Bayesian redshift priors used for the analysis, specifically on the relative abundance of the BBH mergers originated from the first stars, and the primordial BBH mergers.


title: Data-driven expectations for electromagnetic counterpart searches based on LIGO/Virgo public alerts

abstract: Searches for electromagnetic counterparts of gravitational-wave signals have redoubled since the first detection in 2017 of a binary neutron star merger with a gamma-ray burst, optical/infrared kilonova, and panchromatic afterglow. Yet, one LIGO/Virgo observing run later, there has not yet been a second, secure identification of an electromagnetic counterpart. This is not surprising given that the localization uncertainties of events in LIGO and Virgo's third observing run, O3, were much larger than predicted. We explain this by showing that improvements in data analysis that now allow LIGO/Virgo to detect weaker and hence more poorly localized events have increased the overall number of detections, of which well-localized, gold-plated events make up a smaller proportion overall. We present simulations of the next two LIGO/Virgo/KAGRA observing runs, O4 and O5, that are grounded in the statistics of O3 public alerts. To illustrate the significant impact that the updated predictions can have, we study the follow-up strategy for the Zwicky Transient Facility. Realistic and timely forecasting of gravitational-wave localization accuracy is paramount given the large commitments of telescope time and the need to prioritize which events are followed up. We include a data release of our simulated localizations as a public proposal planning resource for astronomers.


title: SN 2021csp -- the explosion of a stripped envelope star within a H and He-poor circumstellar medium

abstract: We present observations of SN 2021csp, a unique supernova (SN) which displays evidence for interaction with H- and He- poor circumstellar material (CSM) at early times. Using high-cadence spectroscopy taken over the first week after explosion, we show that the spectra of SN 2021csp are dominated by C III lines with a velocity of 1800 km s$^{-1}$. We associate this emission with CSM lost by the progenitor prior to explosion. Subsequently, the SN displays narrow He lines before metamorphosing into a broad-lined Type Ic SN. We model the bolometric light curve of SN 2021csp, and show that it is consistent with the energetic ($4\times10^{51}$ erg) explosion of a stripped star, producing 0.4 M$_\odot$ of 56Ni within a $\sim$1 M$_\odot$ shell of CSM extending out to 400 R$_\odot$.


title: Improving constraints on the reionization parameters using 21-cm bispectrum

abstract: Radio interferometric experiments aim to constrain the reionization model parameters by measuring the 21-cm signal statistics, primarily the power spectrum. However the Epoch of Reionization (EoR) 21-cm signal is highly non-Gaussian, and this non-Gaussianity encodes important information about this era. The bispectrum is the lowest order statistic able to capture this inherent non-Gaussianity. Here we are the first to demonstrate that bispectra for large and intermediate length scales and for all unique $k$-triangle shapes provide tighter constraints on the EoR parameters compared to the power spectrum or the bispectra for a limited number of shapes of $k$-triangles. We use the Bayesian inference technique to constrain EoR parameters. We have also developed an Artificial Neural Network (ANN) based emulator for the EoR 21-cm power spectrum and bispectrum which we use to remarkably speed up our parameter inference pipeline. Here we have considered the sample variance and the system noise uncertainties corresponding to $1000$ hrs of SKA-Low observations for estimating errors in the signal statistics. We find that using all unique $k$-triangle bispectra improves the constraints on parameters by a factor of $2-4$ (depending on the stage of reionization) over the constraints that are obtained using power spectrum alone.


title: Quantum Extremal Surfaces and the Holographic Entropy Cone

abstract: Quantum states with geometric duals are known to satisfy a stricter set of entropy inequalities than those obeyed by general quantum systems. The set of allowed entropies derived using the Ryu-Takayanagi (RT) formula defines the Holographic Entropy Cone (HEC). These inequalities are no longer satisfied once general quantum corrections are included by employing the Quantum Extremal Surface (QES) prescription. Nevertheless, the structure of the QES formula allows for a controlled study of how quantum contributions from bulk entropies interplay with HEC inequalities. In this paper, we initiate an exploration of this problem by relating bulk entropy constraints to boundary entropy inequalities. In particular, we show that requiring the bulk entropies to satisfy the HEC implies that the boundary entropies also satisfy the HEC. Further, we also show that requiring the bulk entropies to obey monogamy of mutual information (MMI) implies the boundary entropies also obey MMI.


title: Where, when and why: occurrence of fast-pairwise collective neutrino oscillation in three-dimensional core-collapse supernova models

abstract: Fast-pairwise collective neutrino oscillation represents a key uncertainty in the theory of core-collapse supernova (CCSN). Despite the potentially significant impact on CCSN dynamics, it is usually neglected in numerical models of CCSN because of the formidable technical difficulties of self-consistently incorporating this physics. In this paper, we investigate the prospects for the occurrence of fast flavor conversion by diagnosing electron neutrino lepton number (ELN) crossing in more than a dozen state-of-the-art three-dimensional CCSN models. ELN crossings is a necessary condition for triggering flavor conversion. Although only zeroth and first angular moments are available from the simulations, our new method enables us to look into the angular distributions of neutrinos in momentum space and provide accurate insight into ELN crossings. Our analysis suggests that fast flavor conversion generally occurs in the post-shock region of CCSNe, and that explosive models provide more favorable conditions for the flavor conversion than failed CCSNe. We also find that there are both common and progenitor-dependent characteristics. Classifying ELN crossings into two types, we analyze the generation mechanism of each case by scrutinizing the neutrino radiation field and matter interactions. We find key ingredients of CCSN dynamics driving the ELN crossings: proto-neutron star (PNS) convection, asymmetric neutrino emission, neutrino absorptions and scatterings. This study suggests that we need to accommodate fast flavor conversions in realistic CCSN models.


title: HERA Phase I Limits on the Cosmic 21-cm Signal: Constraints on Astrophysics and Cosmology During the Epoch of Reionization

abstract: Recently, the Hydrogen Epoch of Reionization Array (HERA) collaboration has produced the experiment's first upper limits on the power spectrum of 21-cm fluctuations at z~8 and 10. Here, we use several independent theoretical models to infer constraints on the intergalactic medium (IGM) and galaxies during the epoch of reionization (EoR) from these limits. We find that the IGM must have been heated above the adiabatic cooling threshold by z~8, independent of uncertainties about the IGM ionization state and the nature of the radio background. Combining HERA limits with galaxy and EoR observations constrains the spin temperature of the z~8 neutral IGM to 27 K < T_S < 630 K (2.3 K < T_S < 640 K) at 68% (95%) confidence. They therefore also place a lower bound on X-ray heating, a previously unconstrained aspects of early galaxies. For example, if the CMB dominates the z~8 radio background, the new HERA limits imply that the first galaxies produced X-rays more efficiently than local ones (with soft band X-ray luminosities per star formation rate constrained to L_X/SFR = { 10^40.2, 10^41.9 } erg/s/(M_sun/yr) at 68% confidence), consistent with expectations of X-ray binaries in low-metallicity environments. The z~10 limits require even earlier heating if dark-matter interactions (e.g., through millicharges) cool down the hydrogen gas. Using a model in which an extra radio background is produced by galaxies, we rule out (at 95% confidence) the combination of high radio and low X-ray luminosities of L_{r,\nu}/SFR > 3.9 x 10^24 W/Hz/(M_sun/yr) and L_X/SFR<10^40 erg/s/(M_sun/yr). The new HERA upper limits neither support nor disfavor a cosmological interpretation of the recent EDGES detection. The analysis framework described here provides a foundation for the interpretation of future HERA results.


title: Sub-arcsecond imaging with the International LOFAR Telescope I. Foundational calibration strategy and pipeline

abstract: [abridged] The International LOFAR Telescope is an interferometer with stations spread across Europe. With baselines of up to ~2,000 km, LOFAR has the unique capability of achieving sub-arcsecond resolution at frequencies below 200 MHz, although this is technically and logistically challenging. Here we present a calibration strategy that builds on previous high-resolution work with LOFAR. We give an overview of the calibration strategy and discuss the special challenges inherent to enacting high-resolution imaging with LOFAR, and describe the pipeline, which is publicly available, in detail. We demonstrate the calibration strategy by using the pipeline on P205+55, a typical LOFAR Two-metre Sky Survey (LoTSS) pointing. We perform in-field delay calibration, solution referencing to other calibrators, self-calibration, and imaging of example directions of interest in the field. For this specific field and these ionospheric conditions, dispersive delay solutions can be transferred between calibrators up to ~1.5 degrees away, while phase solution transferral works well over 1 degree. We demonstrate a check of the astrometry and flux density scale. Imaging in 17 directions, the restoring beam is typically 0.3" x 0.2" although this varies slightly over the entire 5 square degree field of view. We achieve ~80 to 300 $\mu$Jy/bm image rms noise, which is dependent on the distance from the phase centre; typical values are ~90 $\mu$Jy/bm for the 8 hour observation with 48 MHz of bandwidth. Seventy percent of processed sources are detected, and from this we estimate that we should be able to image ~900 sources per LoTSS pointing. This equates to ~3 million sources in the northern sky, which LoTSS will entirely cover in the next several years. Future optimisation of the calibration strategy for efficient post-processing of LoTSS at high resolution (LoTSS-HR) makes this estimate a lower limit.


title: Sub-arcsecond imaging with the International LOFAR Telescope: II. Completion of the LOFAR Long-Baseline Calibrator Survey

abstract: The Low-Frequency Array (LOFAR) Long-Baseline Calibrator Survey (LBCS) was conducted between 2014 and 2019 in order to obtain a set of suitable calibrators for the LOFAR array. In this paper we present the complete survey, building on the preliminary analysis published in 2016 which covered approximately half the survey area. The final catalogue consists of 30006 observations of 24713 sources in the northern sky, selected for a combination of high low-frequency radio flux density and flat spectral index using existing surveys (WENSS, NVSS, VLSS, and MSSS). Approximately one calibrator per square degree, suitable for calibration of $\geq$ 200 km baselines is identified by the detection of compact flux density, for declinations north of 30 degrees and away from the Galactic plane, with a considerably lower density south of this point due to relative difficulty in selecting flat-spectrum candidate sources in this area of the sky. Use of the VLBA calibrator list, together with statistical arguments by comparison with flux densities from lower-resolution catalogues, allow us to establish a rough flux density scale for the LBCS observations, so that LBCS statistics can be used to estimate compact flux densities on scales between 300 mas and 2 arcsec, for sources observed in the survey. The LBCS can be used to assess the structures of point sources in lower-resolution surveys, with significant reductions in the degree of coherence in these sources on scales between 2 arcsec and 300 mas. The LBCS survey sources show a greater incidence of compact flux density in quasars than in radio galaxies, consistent with unified schemes of radio sources. Comparison with samples of sources from interplanetary scintillation (IPS) studies with the Murchison Widefield Array (MWA) shows consistent patterns of detection of compact structure in sources observed both interferometrically with LOFAR and using IPS.


title: NLO electroweak potentials for minimal dark matter and beyond

abstract: We calculate the one-loop correction to the static potential induced by $\gamma$, $W$ and $Z$-exchange at tree-level for arbitrary heavy standard model multiplets. We find that the result obeys a "Casimir-like" scaling, making the NLO correction to the potential a "low-energy" property of the electroweak gauge bosons. Furthermore, we discuss the phenomenology of the NLO potentials, the analytically known asymptotic limits and provide fitting functions in position space for easy use of the results.


title: Pushing subarcsecond resolution imaging down to 30 MHz with the trans-European International LOFAR Telescope

abstract: Relatively little information is available about the Universe at ultra-low radio frequencies, i.e. below 50 MHz (ULF), although the ULF spectral window contains a wealth of unique diagnostics for studying galactic and extragalactic phenomena. Sub-arcsecond resolution imaging at these frequencies is extremely difficult, due to the long baselines (>1000 km) required and large ionospheric perturbations. We have conducted a pilot project to investigate the ULF performance and potential of the International LOFAR Telescope (ILT), a trans-European interferometric array with baselines up to ~2000 km and observing frequencies down to 10 MHz. We have successfully produced images with sub-arcsecond resolution for 6 radio sources at frequencies down to 30 MHz. This is more than an order of magnitude better resolution than pre-ILT observations at similar frequencies. The six targets that we have imaged (3C 196, 3C 225, 3C 273, 3C 295, 3C 298 and 3C 380) are bright radio sources with compact structures. By comparing our data of 3C 196 and 3C 273 with observations at higher frequencies, we investigate their spatially resolved radio spectral properties. Our success shows that at frequencies down to 30 MHz, sub-arcsecond imaging with the ILT is possible. Further analysis is needed to determine the feasibility of observations of fainter sources or sources with less compact emission.


title: Origin of the ring structures in Hercules A -- Sub-arcsecond 144 MHz to 7 GHz observations

abstract: The prominent radio source Hercules A features complex structures in its radio lobes. Although it is one of the most comprehensively studied sources in the radio sky, the origin of the ring structures in the Hercules A radio lobes remains an open question. We present the first sub-arcsecond angular resolution images at low frequencies (<300 MHz) of Hercules A, made with the International LOFAR Telescope. With the addition of data from the Karl G. Jansky Very Large Array, we mapped the structure of the lobes from 144 MHz to 7 GHz. We explore the origin of the rings within the lobes of Hercules A, and test whether their properties are best described by a shock model, where shock waves are produced by the jet propagating in the radio lobe, or by an inner-lobe model, where the rings are formed by decelerated jetted plasma. From spectral index mapping our large frequency coverage reveals that the curvature of the different ring spectra increases with distance away from the central active galactic nucleus. We demonstrate that the spectral shape of the rings is consistent with synchrotron aging, which speaks in favor of an inner-lobe model where the rings are formed from the deposition of material from past periods of intermittent core activity.


title: The resolved jet of 3C 273 at 150 MHz

abstract: Since its discovery in 1963, 3C273 has become one of the most widely studied quasars with investigations spanning the electromagnetic spectrum. While much has been discovered about this historically notable source, its low-frequency emission is far less well understood. Observations in the MHz regime have traditionally lacked the resolution required to explore small-scale structures that are key to understanding the processes that result in the observed emission. In this paper we use the first sub-arcsecond images of 3C273 at MHz frequencies to investigate the morphology of the compact jet structures and the processes that result in the observed spectrum. Using the full complement of LOFAR's international stations, we produce $0.31 \times 0.21$ arcsec images of 3C273 at 150 MHz to determine the jet's kinetic power, place constraints on the bulk speed and inclination angle of the jets, and look for evidence of the elusive counter-jet at 150 MHz. Using ancillary data at GHz frequencies, we fit free-free absorption (FFA) and synchrotron self-absorption (SSA) models to determine their validity in explaining the observed spectra. The images presented display for the first time that robust, high-fidelity imaging of low-declination complex sources is now possible with the LOFAR international baselines. We show that the main small-scale structures of 3C273 match those seen at higher frequencies and that absorption is present in the observed emission. We determine the kinetic power of the jet to be in the range of $3.5 \times 10^{43}$ - $1.5 \times 10^{44}$ erg s$^{-1}$ which agrees with estimates made using higher frequency observations. We derive lower limits for the bulk speed and Lorentz factor of $\beta \gtrsim 0.55$ and $\Gamma \geq 1.2$ respectively. The counter-jet remains undetected at $150$ MHz, placing a limit on the peak brightness of $S_\mathrm{cj\_150} < 40$ mJy beam$^{-1}$.


title: Unmasking the history of 3C 293 with LOFAR sub-arcsecond imaging

abstract: Active galactic nuclei (AGNs) show episodic activity, evident in galaxies that exhibit restarted radio jets. These restarted jets can interact with their environment, leaving signatures on the radio spectral energy distribution. Tracing these signatures requires resolved spectral index measurements over a broad frequency range including low frequencies. We present such a study for the radio galaxy 3C 293. Using the International LOFAR telescope (ILT) we probed spatial scales as fine as ~0.2" at 144 MHz, and to constrain the spectrum we combined these data with Multi-Element Radio Linked Interferometer Network (MERLIN) and Very Large Array (VLA) archival data. In the inner lobes (~2 kpc), we detect the presence of a spectral turnover that peaks at ~225 MHz and is most likely caused by free-free absorption from the rich surrounding medium. We confirm that these inner lobes are part of a jet-dominated young radio source (spectral age $\lesssim$0.17 Myr), which is strongly interacting with the rich interstellar medium (ISM) of the host galaxy. The outer lobes (~100 kpc) have a spectral index of $\alpha$~0.6-0.8 from 144-4850 MHz with a remarkably uniform spatial distribution and only mild spectral curvature ($\Delta\alpha\lesssim$ 0.2). We propose that intermittent fuelling and jet flow disruptions are powering the mechanisms that keep the spectral index in the outer lobes from steepening and maintain the spatial uniformity of the spectral index. Overall, it appears that 3C 293 has gone through multiple (two to three) epochs of activity. This study adds 3C 293 to the new sub-group of restarted galaxies with short interruption time periods. This is the first time a spatially resolved study simultaneously studies a young source as well as the older outer lobes at such low frequencies. This illustrates the potential of the ILT to expand such studies to a larger sample of radio galaxies.


title: High-resolution international LOFAR observations of 4C~43.15 -- Spectral ages and injection indices in a high-z radio galaxy

abstract: Radio sources with steep spectra are preferentially associated with the most distant galaxies, the $\alpha-z$ relation, but the reason for this relation is an open question. The spatial distribution of spectra in high-z radio sources can be used to study this relation, and low-frequency observations are particularly important in understanding the particle acceleration and injection mechanisms. However, the small angular sizes of high-z sources together with the inherently low resolution of low-frequency radio telescopes until now has prevented high angular resolution low-frequency observations of distant objects. Here we present subarcsecond observations of a $z = 2.4$ radio galaxy at frequencies between $121$ MHz and $166$ MHz. We measure the spatial distribution of spectra, and discuss the implications for models of the $\alpha-z$ relation. We targeted 4C 43.15 with the High Band Antennas (HBAs) of the \textit{International LOFAR Telescope} (ILT) with a range of baselines up to $1300\ \mathrm{km}$. At the central frequency of $143$ MHz we achieve an angular resolution of $\sim 0.3''$. By complementing our data with archival \textit{Very Large Array} (VLA) data we study the spectral index distribution across 4C 43.15 between $55\ \mathrm{MHz}$ and $8.4\ \mathrm{GHz}$ at resolutions of $0.4''$ and $0.9''$. With a magnetic field strength of $B = 5.2$ nT and fitted injection indices of $\alpha^\mathrm{north}_\mathrm{inj} = -0.8$ and $\alpha^\mathrm{south}_\mathrm{inj} = -0.6$, fitting a Tribble spectral ageing model results in a spectral age of $\tau_\mathrm{spec} = 1.1 \pm 0.1$ Myr. We conclude that our data on 4C 43.15 indicates that inverse Compton losses could become comparable to or exceed synchrotron losses at higher redshifts and that inverse Compton losses could be a viable explanation for the $\alpha-z$ relation.


title: Subarcsecond LOFAR imaging of Arp299 at 150 MHz. Tracing the nuclear and diffuse extended emission of a bright LIRG

abstract: We study for the first time the low-frequency ($\sim$150 MHz) radio brightness distribution of Arp~299 at subarcsecond resolution, tracing in both compact and extended emission regions the local spectral energy distribution (SED) in order to characterize the dominant emission and absorption processes. We analysed the spatially resolved emission of Arp 299 revealed by 150 MHz international baseline Low-Frequency Array (LOFAR) and 1.4, 5.0, and 8.4 GHz Very Large Array (VLA) observations. We present the first subarcsecond (0.4"$\sim$100~pc) image of the whole Arp~299 system at 150~MHz. The high surface brightness sensitivity of our LOFAR observations ($\sim$100 $\mu$Jy/beam) allowed us to detect all of the nuclear components detected at higher frequencies, as well as the extended steep-spectrum emission surrounding the nuclei. We obtained spatially resolved, two-point spectral index maps for the whole galaxy: the compact nuclei show relatively flat spectra, while the extended, diffuse component shows a steep spectrum. We fitted the radio SED of the nuclear regions using two different models: a continuous free-free medium model and a clumpy model. The continuous model can explain the SED of the nuclei assuming a population of relativistic electrons subjected to synchrotron, bremsstrahlung, and ionization losses. The clumpy model fits assuming relativistic electrons with negligible energy losses, and thermal fractions that are more typical of star-forming galaxies than those required for the continuous model. Our results confirm the usefulness of combining spatially resolved radio imaging at both MHz and GHz frequencies to characterize in detail the radio emission properties of LIRGs from the central 100 pc out to the kiloparsec galaxy-wide scales.


title: Supermeasured: Violating Bell-Statistical Independence without violating physical statistical independence

abstract: Bell's theorem is often said to imply that quantum mechanics violates local causality, and that local causality cannot be restored with a hidden-variables theory. This however is only correct if the hidden-variables theory fulfils an assumption called Statistical Independence. Violations of Statistical Independence are commonly interpreted as correlations between the measurement settings and the hidden variables (which determine the measurement outcomes). Such correlations have been discarded as ``fine-tuning'' or a ``conspiracy''. We here point out that the common interpretation is at best physically ambiguous and at worst incorrect. The problem with the common interpretation is that Statistical Independence might be violated because of a non-trivial measure in state space, a possibility we propose to call ``supermeasured''. We use Invariant Set Theory as an example of a supermeasured theory that violates the Statistical Independence assumption in Bell's theorem without requiring correlations between hidden variables and measurement settings (physical statistical independence).


title: High-resolution imaging with the International LOFAR Telescope: Observations of the gravitational lenses MG 0751+2716 and CLASS B1600+434

abstract: We present Low-Frequency Array (LOFAR) telescope observations of the radio-loud gravitational lens systems MG 0751+2716 and CLASS B1600+434. These observations produce images at 300 milliarcseconds (mas) resolution at 150 MHz. In the case of MG 0751+2716, lens modelling is used to derive a size estimate of around 2 kpc for the low-frequency source, which is consistent with a previous 27.4 GHz study in the radio continuum with Karl G. Jansky Very Large Array (VLA). This consistency implies that the low-frequency radio source is cospatial with the core-jet structure that forms the radio structure at higher frequencies, and no significant lobe emission or further components associated with star formation are detected within the magnified region of the lens. CLASS B1600+434 is a two-image lens where one of the images passes through the edge-on spiral lensing galaxy, and the low radio frequency allows us to derive limits on propagation effects, namely scattering, in the lensing galaxy. The observed flux density ratio of the two lensed images is 1.19 +/- 0.04 at an observed frequency of 150 MHz. The widths of the two images give an upper limit of 0.035 kpc m^-20/3 on the integrated scattering column through the galaxy at a distance approximately 1 kpc above its plane, under the assumption that image A is not affected by scattering. This is relatively small compared to limits derived through very long baseline interferometry (VLBI) studies of differential scattering in lens systems. These observations demonstrate that LOFAR is an excellent instrument for studying gravitational lenses. We also report on the inability to calibrate three further lens observations: two from early observations that have less well determined station calibration, and a third observation impacted by phase transfer problems.


title: Spectral analysis of spatially-resolved 3C295 (sub-arcsecond resolution) with the International LOFAR Telescope

abstract: 3C295 is a bright, compact steep spectrum source with a well-studied integrated radio spectral energy distribution (SED) from 132 MHz to 15 GHz. However, spatially resolved spectral studies have been limited due to a lack of high resolution images at low radio frequencies. These frequencies are crucial for measuring absorption processes, and anchoring the overall spectral modelling of the radio SED. In this paper, we use International LOFAR (LOw-Frequency ARray) Telescope (ILT) observations of 3C295 to study its spatially resolved spectral properties with sub-arcsecond resolution at 132 MHz. Combining our new 132 MHz observation with archival data at 1.6 GHz, 4.8 GHz, and 15 GHz, we are able to carry out a resolved radio spectral analysis. The spectral properties of the hotspots provides evidence for low frequency flattening. In contrast, the spectral shape across the lobes is consistent with a JP spectral ageing model. Using the integrated spectral information for each component, we then fit low-frequency absorption models to the hotspots, finding that both free-free absorption and synchrotron self-absorption models provide a better fit to the data than a standard power law. Although we can say there is low-frequency absorption present in the hotspots of 3C295, future observations with the Low Band Antenna of the ILT at 55 MHz may allow us to distinguish the type of absorption.


title: The old conformal bootstrap revisited

abstract: The "old" conformal bootstrap was originally formulated by Migdal and Polyakov (MP) as a method for calculating conformal dimensions self-consistently. In this work we revisit the MP bootstrap and apply efficient multi-loop Feynman integral techniques in order to solve the corresponding equations. We obtain solutions at first order in the skeleton expansion for finite coupling and for integer values of the spacetime dimension. We focus in particular on the $\phi^3$ theory in seven dimensions and the $O(N)$ vector models in three dimensions.


title: A critical appraisal of the singularity theorems

abstract: The 2020 Nobel prize in Physics has revived the interest in the singularity theorems and, in particular, in the Penrose theorem published in 1965. In this short paper I briefly review the main ideas behind the theorems and then proceed to an evaluation of their hypotheses and implications. I will try to dispel some common misconceptions about the theorems and their conclusions, as well as to convey some of their rarely mentioned consequences. In particular, a discussion of spacetime extensions in relation to the theorems is provided. The nature of the singularity inside black holes is also analyzed.


title: K_4-free graphs have sparse halves

abstract: Every $K_4$-free graph on $n$ vertices has a set of $\lfloor n/2\rfloor$ vertices spanning at most $n^2/18$ edges.


title: Cosmological functional renormalisation group, extended Galilean invariance and the sweeping effect

abstract: The functional renormalisation group is employed to study the non-linear regime of late-time cosmic structure formation. This framework naturally allows for non-perturbative approximation schemes, usually guided by underlying symmetries or a truncation of the theory space. An extended symmetry that is related to Galilean invariance is studied and corresponding Ward identities are derived. These are used to obtain (formally) closed renormalisation group flow equations for two-point correlation functions in the limit of large wave numbers (small scales). The flow equations are analytically solved in an approximation that is connected to the 'sweeping effect' previously described in the context of fluid turbulence.


title: Dirac type relativistic quantum mechanics for massive photons

abstract: Constructing a relativistic quantum mechanics (RQM) for a massive photon, without appealing to the quantum field theoretical approach such as the massive Proca model, we find a new theoretical particle solution which allows the massive photon having either positive or negative energy as solutions. In particular, we predict the existence of the so-called anti-photon corresponding to the negative energy solution, similar to the positron in the Dirac RQM for an electron. Note that the anti-photons could be an intense radiation flare of the gamma ray burst. In this RQM for the massive photon, we construct a positive definite probability density and a nontrivial diagonal Hamiltonian, and also discuss a massless photon. Moreover we confirm the covariances of the relativistic equation of motion for the massive photon and the corresponding probability continuity equation under the Lorentz transformation.


title: A Numerical Method for a Nonlocal Diffusion Equation with Additive Noise

abstract: We consider a nonlocal evolution equation representing the continuum limit of a large ensemble of interacting particles on graphs forced by noise. The two principle ingredients of the continuum model are a nonlocal term and Q-Wiener process describing the interactions among the particles in the network and stochastic forcing respectively. The network connectivity is given by a square integrable function called a graphon. We prove that the initial value problem for the continuum model is well-posed. Further, we construct a semidiscrete (discrete in space and continuous in time) and a fully discrete schemes for the nonlocal model. The former is obtained by a discontinuous Galerkin method and the latter is based on further discretizing time using the Euler-Maruyama method. We prove convergence and estimate the rate of convergence in each case. For the semidiscrete scheme, the rate of convergence estimate is expressed in terms of the regularity of the graphon, Q-Wiener process, and the initial data. We work in generalized Lipschitz spaces, which allows to treat models with data of lower regularity. This is important for applications as many interesting types of connectivity including small-world and power-law are expressed by graphons that are not smooth. The error analysis of the fully discrete scheme, on the other hand, reveals that for some models common in applied science, one has a higher speed of convergence than that predicted by the standard estimates for the Euler-Maruyama method. The rate of convergence analysis is supplemented with detailed numerical experiments, which are consistent with our analytical results. As a by-product, this work presents a rigorous justification for taking continuum limit for a large class of interacting dynamical systems on graphs subject to noise.


title: Understanding the factors driving the opioid epidemic using machine learning

abstract: In recent years, the US has experienced an opioid epidemic with an unprecedented number of drugs overdose deaths. Research finds such overdose deaths are linked to neighborhood-level traits, thus providing opportunity to identify effective interventions. Typically, techniques such as Ordinary Least Squares (OLS) or Maximum Likelihood Estimation (MLE) are used to document neighborhood-level factors significant in explaining such adverse outcomes. These techniques are, however, less equipped to ascertain non-linear relationships between confounding factors. Hence, in this study we apply machine learning based techniques to identify opioid risks of neighborhoods in Delaware and explore the correlation of these factors using Shapley Additive explanations (SHAP). We discovered that the factors related to neighborhoods environment, followed by education and then crime, were highly correlated with higher opioid risk. We also explored the change in these correlations over the years to understand the changing dynamics of the epidemic. Furthermore, we discovered that, as the epidemic has shifted from legal (i.e., prescription opioids) to illegal (e.g.,heroin and fentanyl) drugs in recent years, the correlation of environment, crime and health related variables with the opioid risk has increased significantly while the correlation of economic and socio-demographic variables has decreased. The correlation of education related factors has been higher from the start and has increased slightly in recent years suggesting a need for increased awareness about the opioid epidemic.