title: Analytical modeling of micelle growth. 4. Molecular thermodynamics of wormlike micelles from ionic surfactants: theory vs. experiment

abstract: Hypotheses: The aggregation number and length of spherocylindrical (rodlike, wormlike) micelles in solutions of an ionic surfactant and salt can be predicted knowing the molecular parameters and the input concentrations of the species. This can be achieved by upgrading the quantitative molecular thermodynamic model from the previous parts of this series with an expression for the electrostatic component of micelle scission energy that is the excess free energy of the spherical endcaps with respect to the cylindrical part of the micelle. Theory: The thermodynamics of micellization is extended to the case of multicomponent system, which may contain several surfactants (both ionic and nonionic) and salts, taking into account the effect of counterion binding in the Stern layer on the micellar surface. Furthermore, the considerations are focused on a system that consists of single ionic surfactant plus salt. Findings: Excellent agreement was achieved between the theoretical model and experimental data for wormlike micelles from anionic and cationic surfactants at various concentrations of salt and temperatures. In accord with the experimental observations, at high salt concentrations, the model predicts loss of chemical equilibrium between the endcaps and cylindrical part of the wormlike micelles, which implies transition to self-assemblies of other, e.g. branched, morphology or the onset of crystallization and phase separation.


title: MPEG Media Enablers For Richer XR Experiences

abstract: With the advent of immersive media applications, the requirements for the representation and the consumption of such content has dramatically increased. The ever-increasing size of the media asset combined with the stringent motion-to-photon latency requirement makes the equation of a high quality of experience for XR streaming services difficult to solve. The MPEG-I standards aim at facilitating the wide deployment of immersive applications. This paper describes part 13, Video Decoding Interface, and part 14, Scene Description for MPEG Media of MPEG-I which address decoder management and the virtual scene composition, respectively. These new parts intend to make complex media rendering operations and hardware resources management hidden from the application, hence lowering the barrier for XR application to become mainstream and accessible to XR experience developers and designers. Both parts are expected to be published by ISO at the end of 2021.


title: Deep RL With Information Constrained Policies: Generalization in Continuous Control

abstract: Biological agents learn and act intelligently in spite of a highly limited capacity to process and store information. Many real-world problems involve continuous control, which represents a difficult task for artificial intelligence agents. In this paper we explore the potential learning advantages a natural constraint on information flow might confer onto artificial agents in continuous control tasks. We focus on the model-free reinforcement learning (RL) setting and formalize our approach in terms of an information-theoretic constraint on the complexity of learned policies. We show that our approach emerges in a principled fashion from the application of rate-distortion theory. We implement a novel Capacity-Limited Actor-Critic (CLAC) algorithm and situate it within a broader family of RL algorithms such as the Soft Actor Critic (SAC) and Mutual Information Reinforcement Learning (MIRL) algorithm. Our experiments using continuous control tasks show that compared to alternative approaches, CLAC offers improvements in generalization between training and modified test environments. This is achieved in the CLAC model while displaying the high sample efficiency of similar methods.


title: Learning Invariant Representations and Risks for Semi-supervised Domain Adaptation

abstract: The success of supervised learning hinges on the assumption that the training and test data come from the same underlying distribution, which is often not valid in practice due to potential distribution shift. In light of this, most existing methods for unsupervised domain adaptation focus on achieving domain-invariant representations and small source domain error. However, recent works have shown that this is not sufficient to guarantee good generalization on the target domain, and in fact, is provably detrimental under label distribution shift. Furthermore, in many real-world applications it is often feasible to obtain a small amount of labeled data from the target domain and use them to facilitate model training with source data. Inspired by the above observations, in this paper we propose the first method that aims to simultaneously learn invariant representations and risks under the setting of semi-supervised domain adaptation (Semi-DA). First, we provide a finite sample bound for both classification and regression problems under Semi-DA. The bound suggests a principled way to obtain target generalization, i.e. by aligning both the marginal and conditional distributions across domains in feature space. Motivated by this, we then introduce the LIRR algorithm for jointly \textbf{L}earning \textbf{I}nvariant \textbf{R}epresentations and \textbf{R}isks. Finally, extensive experiments are conducted on both classification and regression tasks, which demonstrates LIRR consistently achieves state-of-the-art performance and significant improvements compared with the methods that only learn invariant representations or invariant risks.


title: On microscopic entropy production, heat and fluctuation theorem

abstract: We demonstrate that the Gibbs-Shannon entropy is applicable to non-equilibrium systems of any size and boundary conditions. The change in microscopic entropy can be attributed to the stochastic nature of dynamic processes and to the inherent uncertainties of thermodynamic systems. The latter predicts that entropy production is nonnegative on average and varies with different trajectories according to the fluctuation theorem. By contrast, heat is affiliated with stochastic processes underlying particle motions and the ensemble average over all possible trajectories leads to the Clausius inequality. The Jarzynski/Crooks equations can be readily derived by applying the fluctuation theorem to heat variation over different trajectories linking equilibrium states.


title: A surrogate-based optimal likelihood function for the Bayesian calibration of catalytic recombination in atmospheric entry protection materials

abstract: This work deals with the inference of catalytic recombination parameters from plasma wind tunnel experiments for reusable thermal protection materials. One of the critical factors affecting the performance of such materials is the contribution to the heat flux of the exothermic recombination reactions at the vehicle surface. The main objective of this work is to develop a dedicated Bayesian framework that allows us to compare uncertain measurements with model predictions which depend on the catalytic parameter values. Our framework accounts for uncertainties involved in the model definition and incorporates all measured variables with their respective uncertainties. The physical model used for the estimation consists of a 1D boundary layer solver along the stagnation line. The chemical production term included in the surface mass balance depends on the catalytic recombination efficiency. As not all the different quantities needed to simulate a reacting boundary layer can be measured or known (such as the flow enthalpy at the inlet boundary), we propose an optimization procedure built on the construction of the likelihood function to determine their most likely values based on the available experimental data. This procedure avoids the need to introduce any a priori estimates on the nuisance quantities, namely, the boundary layer edge enthalpy, wall temperatures, static and dynamic pressures, which would entail the use of very wide priors. We substitute the optimal likelihood of the experimental data with a surrogate model to make the inference procedure both faster and more robust. We show that the resulting Bayesian formulation yields meaningful and accurate posterior distributions of the catalytic parameters with a reduction of more than 20% of the standard deviation with respect to previous works. We also study the implications of an extension of the experimental procedure.

categories: physics.data-an cs.NA math.NA physics.flu-dyn


title: LSTMs Compose (and Learn) Bottom-Up

abstract: Recent work in NLP shows that LSTM language models capture hierarchical structure in language data. In contrast to existing work, we consider the \textit{learning} process that leads to their compositional behavior. For a closer look at how an LSTM's sequential representations are composed hierarchically, we present a related measure of Decompositional Interdependence (DI) between word meanings in an LSTM, based on their gate interactions. We connect this measure to syntax with experiments on English language data, where DI is higher on pairs of words with lower syntactic distance. To explore the inductive biases that cause these compositional representations to arise during training, we conduct simple experiments on synthetic data. These synthetic experiments support a specific hypothesis about how hierarchical structures are discovered over the course of training: that LSTM constituent representations are learned bottom-up, relying on effective representations of their shorter children, rather than learning the longer-range relations independently from children.


title: Wildfire Smoke and Air Quality: How Machine Learning Can Guide Forest Management

abstract: Prescribed burns are currently the most effective method of reducing the risk of widespread wildfires, but a largely missing component in forest management is knowing which fuels one can safely burn to minimize exposure to toxic smoke. Here we show how machine learning, such as spectral clustering and manifold learning, can provide interpretable representations and powerful tools for differentiating between smoke types, hence providing forest managers with vital information on effective strategies to reduce climate-induced wildfires while minimizing production of harmful smoke.


title: Towards Social HRI for Improving Children's Healthcare Experiences

abstract: This paper describes a new research project that aims to develop a social robot designed to help children cope with painful and distressing medical procedures in a clinical setting. While robots have previously been trialled for this task, with promising initial results, the systems have tended to be teleoperated, limiting their flexibility and robustness. This project will use epistemic planning techniques as a core component for action selection in the robot system, in order to generate plans that include physical, sensory, and social actions for interacting with humans. The robot will operate in a task environment where appropriate and safe interaction with children, parents/caregivers, and healthcare professionals is required. In addition to addressing the core technical challenge of building an autonomous social robot, the project will incorporate co-design techniques involving all participant groups, and the final robot system will be evaluated in a two-site clinical trial.


title: Quantifying the multi-objective cost of uncertainty

abstract: Various real-world applications involve modeling complex systems with immense uncertainty and optimizing multiple objectives based on the uncertain model. Quantifying the impact of the model uncertainty on the given operational objectives is critical for designing optimal experiments that can most effectively reduce the uncertainty that affect the objectives pertinent to the application at hand. In this paper, we propose the concept of mean multi-objective cost of uncertainty (multi-objective MOCU) that can be used for objective-based quantification of uncertainty for complex uncertain systems considering multiple operational objectives. We provide several illustrative examples that demonstrate the concept and strengths of the proposed multi-objective MOCU. Furthermore, we present a real-world example based on the mammalian cell cycle network to demonstrate how the multi-objective MOCU can be used for quantifying the operational impact of model uncertainty when there are multiple, possibly competing, objectives.


title: Coupling electrons and vibrations in molecular quantum chemistry

abstract: We derive an electron-vibration model Hamiltonian in a quantum chemical framework, and explore the extent to which such a Hamiltonian can capture key effects of nonadiabatic dynamics. The model Hamiltonian is a simple two-body operator, and we make preliminary steps at applying standard quantum chemical methods to evaluating its properties, including mean-field theory, linear response, and a primitive correlated model. The Hamiltonian can be compared to standard vibronic Hamiltonians, but is constructed without reference to potential energy surfaces, through direct differentiation of the one- and two-electron integrals at a single reference geometry. The nature of the model Hamiltonian in the harmonic and linear-coupling regime is investigated for pyrazine, where a simple time-dependent calculation including electron-vibration correlation is demonstrated to exhibit the well-studied population transfer between the S$_2$ and S$_1$ excited states.


title: Simultaneous Multicolor Observations of Starlink's Darksat by The Murikabushi Telescope with $MITSuME$

abstract: We present the $g'$-, $R_{\rm c}$-, and $I_{\rm c}$-band magnitudes and associated colors of the Starlink's STARLINK-1113 (one of the standard Starlink satellites) and 1130 (Darksat) with a darkening treatment to its surface. By the 105 cm Murikabushi telescope/$\it{MITSuME}$, simultaneous multicolor observations for the above satellites were conducted four times: on April 10, 2020, May 18, 2020 (for Darksat), and June 11, 2020 (for Darksat and STARLINK-1113). We found that (1) $g'$-band apparent magnitudes of Darksat ($6.95~\pm~0.11-7.65~\pm~0.11$ mag) are comparable to or brighter than that of STARLINK-1113 ($7.69~\pm~0.16$ mag), (2) the shorter the observed wavelength is, the fainter the satellite magnitudes tend to become, (3) the reflected flux by STARLINK-1113 is extremely ($> 1.0$ mag) redder than that of Darksat, (4) there is no clear correlation between the solar phase angle and orbital-altitude-scaled magnitude, (5) by flux-model fitting of the satellite trails with the blackbody radiation, it is found that the albedo of Darksat is about a half of that of STARLINK-1113. Especially, the result (1) is inconsistent with the previous studies. However, considering both solar and observer phase angles and atmospheric extinction, the brightness of STARLINK-1113 can be drastically reduced in the $g'$ and $R_{\rm c}$ band. Simultaneous multicolor-multispot observations more than three colors would give us more detailed information regarding the impact by the low Earth orbit satellite constellations.


title: Time-Variable Radio Recombination Line Emission in W49A

abstract: We present new Jansky Very Large Array (VLA) images of the central region of the W49A star-forming region at 3.6~cm and at 7~mm at resolutions of 0\farcs15 (1650 au) and 0\farcs04 (440 au), respectively. The 3.6~cm data reveal new morphological detail in the ultracompact \ion{H}{2} region population, as well as several previously unknown and unresolved sources. In particular, source A shows elongated, edge-brightened, bipolar lobes, indicative of a collimated outflow, and source E is resolved into three spherical components. We also present VLA observations of radio recombination lines at 3.6~cm and 7~mm, and IRAM Northern Extended Millimeter Array (NOEMA) observations at 1.2~mm. Three of the smallest ultracompact \ion{H}{2} regions (sources A, B2 and G2) all show broad kinematic linewidths, with $\Delta$V$_{FWHM}\gtrsim$40~km~s$^{-1}$. A multi-line analysis indicates that broad linewidths remain after correcting for pressure broadening effects, suggesting the presence of supersonic flows. Substantial changes in linewidth over the 21 year time baseline at both 3.6 cm and 7 mm are found for source G2. At 3.6 cm, the linewidth of G2 changed from 31.7$\pm$1.8 km s$^{-1}$ to 55.6$\pm$2.7 km s$^{-1}$, an increase of $+$23.9$\pm$3.4 km s$^{-1}$. The G2 source was previously reported to have shown a 3.6~cm continuum flux density decrease of 40\% between 1994 and 2015. This source sits near the center of a very young bipolar outflow whose variability may have produced these changes.


title: Epitaxy of new layered materials: 2D chalcogenides and challenges of weak van der Waals interactions

abstract: The application of new materials in nanotechnology opens new perspectives and enables ground-breaking innovations. Two-dimensional van der Waals materials and more specific, 2D chalcogenides are a promising class of new materials awaiting their usage in the semiconductor industry. However, the integration of van der Waals materials relying on industry-compatible manufacturing processes is still a major challenge. This is currently restricting the application of these new materials to the research laboratories environment only. The large-area and single-crystalline growth of van der Waals materials is one of the most important requirements to meet the challenging demands implied by the semiconductor industry. This review contributes to a more generalized understanding on the integration of van der Waals materials - and in more specific 2D chalcogenides - through the growth process of epitaxy. This, can pursue further the aspiration of large-area, single-crystalline and defect-free epitaxial integration of (quasi) van der Waals homo- and heterostructures into the great world of the semiconductor industry.


title: Case Study: Deontological Ethics in NLP

abstract: Recent work in natural language processing (NLP) has focused on ethical challenges such as understanding and mitigating bias in data and algorithms; identifying objectionable content like hate speech, stereotypes and offensive language; and building frameworks for better system design and data handling practices. However, there has been little discussion about the ethical foundations that underlie these efforts. In this work, we study one ethical theory, namely deontological ethics, from the perspective of NLP. In particular, we focus on the generalization principle and the respect for autonomy through informed consent. We provide four case studies to demonstrate how these principles can be used with NLP systems. We also recommend directions to avoid the ethical issues in these systems.


title: On nonnegative solutions for the Functionalized Cahn-Hilliard equation with degenerate mobility

abstract: The Functionalized Cahn-Hilliard equation has been proposed as a model for the interfacial energy of phase-separated mixtures of amphiphilic molecules. We study the existence of a nonnegative weak solutions of a gradient flow of the Functionalized Cahn-Hilliard equation subject to a degenerate mobility M(u) that is zero for u<=0. Assuming the initial data u0(x) is positive, we construct a weak solution as the limit of solutions corresponding to nondegenerate mobilities and verify that it satisfies an energy dissipation inequality.


title: Research, Develop, Deploy: Building a Full Spectrum Software Engineering and Research Department

abstract: At Sandia National Laboratories, the Software Engineering and Research Department seeks to provide sustainable career pathways for research software engineers (RSEs). The conceptual model for our organization follows what we call a Research, Develop, and Deploy (RDD) workflow pattern, enabling RSEs to partner with research and deployment specialists. We argue that this interdisciplinary model allows our department to act as an incubator and an accelerator for impactful ideas. We describe these tactics and our experiences as a RSE team in a scientific computing center.


title: Using Graph Neural Networks for Mass Spectrometry Prediction

abstract: Detecting and quantifying products of cellular metabolism using Mass Spectrometry (MS) has already shown great promise in many biological and biomedical applications. The biggest challenge in metabolomics is annotation, where measured spectra are assigned chemical identities. Despite advances, current methods provide limited annotation for measured spectra. Here, we explore using graph neural networks (GNNs) to predict the spectra. The input to our model is a molecular graph. The model is trained and tested on the NIST 17 LC-MS dataset. We compare our results to NEIMS, a neural network model that utilizes molecular fingerprints as inputs. Our results show that GNN-based models offer higher performance than NEIMS. Importantly, we show that ranking results heavily depend on the candidate set size and on the similarity of the candidates to the target molecule, thus highlighting the need for consistent, well-characterized evaluation protocols for this domain.


title: Deterministic computation of the characteristic polynomial in the time of matrix multiplication

abstract: This paper describes an algorithm which computes the characteristic polynomial of a matrix over a field within the same asymptotic complexity, up to constant factors, as the multiplication of two square matrices. Previously, this was only achieved by resorting to genericity assumptions or randomization techniques, while the best known complexity bound with a general deterministic algorithm was obtained by Keller-Gehrig in 1985 and involves logarithmic factors. Our algorithm computes more generally the determinant of a univariate polynomial matrix in reduced form, and relies on new subroutines for transforming shifted reduced matrices into shifted weak Popov matrices, and shifted weak Popov matrices into shifted Popov matrices.


title: Physical invariance in neural networks for subgrid-scale scalar flux modeling

abstract: In this paper we present a new strategy to model the subgrid-scale scalar flux in a three-dimensional turbulent incompressible flow using physics-informed neural networks (NNs). When trained from direct numerical simulation (DNS) data, state-of-the-art neural networks, such as convolutional neural networks, may not preserve well known physical priors, which may in turn question their application to real case-studies. To address this issue, we investigate hard and soft constraints into the model based on classical transformation invariances and symmetries derived from physical laws. From simulation-based experiments, we show that the proposed transformation-invariant NN model outperforms both purely data-driven ones as well as parametric state-of-the-art subgrid-scale models. The considered invariances are regarded as regularizers on physical metrics during the a priori evaluation and constrain the distribution tails of the predicted subgrid-scale term to be closer to the DNS. They also increase the stability and performance of the model when used as a surrogate during a large-eddy simulation. Moreover, the transformation-invariant NN is shown to generalize to regimes that have not been seen during the training phase.


title: Detection of small magnetic flux ropes from the third and fourth Parker Solar Probe encounters

abstract: We systematically search for magnetic flux rope structures in the solar wind to within the closest distance to the Sun of 0.13 AU, using data from the third and fourth orbits of the Parker Solar Probe. We extend our previous magnetic helicity based technique of identifying magnetic flux rope structures. The method is improved upon to incorporate the azimuthal flow, which becomes larger as the spacecraft approaches the Sun. A total of 21 and 34 magnetic flux ropes are identified during the third (21 days period) and fourth (17 days period) orbits of the Parker Solar Probe, respectively. We provide a statistical analysis of the identified structures, including their relation to the streamer belt and heliospheric current sheet crossing.


title: Scaling Systematic Literature Reviews with Machine Learning Pipelines

abstract: Systematic reviews, which entail the extraction of data from large numbers of scientific documents, are an ideal avenue for the application of machine learning. They are vital to many fields of science and philanthropy, but are very time-consuming and require experts. Yet the three main stages of a systematic review are easily done automatically: searching for documents can be done via APIs and scrapers, selection of relevant documents can be done via binary classification, and extraction of data can be done via sequence-labelling classification. Despite the promise of automation for this field, little research exists that examines the various ways to automate each of these tasks. We construct a pipeline that automates each of these aspects, and experiment with many human-time vs. system quality trade-offs. We test the ability of classifiers to work well on small amounts of data and to generalise to data from countries not represented in the training data. We test different types of data extraction with varying difficulty in annotation, and five different neural architectures to do the extraction. We find that we can get surprising accuracy and generalisability of the whole pipeline system with only 2 weeks of human-expert annotation, which is only 15% of the time it takes to do the whole review manually and can be repeated and extended to new data with no additional effort.


title: Asymmetric prior in wavelet shrinkage

abstract: In bayesian wavelet shrinkage, the already proposed priors to wavelet coefficients are assumed to be symmetric around zero. Although this assumption is reasonable in many applications, it is not general. The present paper proposes the use of an asymmetric shrinkage rule based on the discrete mixture of a point mass function at zero and an asymmetric beta distribution as prior to the wavelet coefficients in a non-parametric regression model. Statistical properties such as bias, variance, classical and bayesian risks of the associated asymmetric rule are provided and performances of the proposed rule are obtained in simulation studies involving artificial asymmetric distributed coefficients and the Donoho-Johnstone test functions. Application in a seismic real dataset is also analyzed.


title: Impact of LHC jet data on Parton Distribution Functions

abstract: In this talk we discuss recent results on the impact of LHC jet data on global fits of parton distribution functions using theoretical predictions at NNLO in QCD supplemented by electroweak (EW) corrections.


title: Low field nano-NMR via three-level system control

abstract: Conventional control strategies for NV centers in quantum sensing are based on a two-level model of their triplet ground state. However, this approach fails in regimes of weak bias magnetic fields or strong microwave pulses, as we demonstrate. To overcome this limitation, we propose a novel control sequence that exploits all three levels by addressing a hidden Raman configuration with microwave pulses tuned to the zero-field transition. We report excellent performance in typical dynamical decoupling sequences, opening up the possibility for nano-NMR operation in low-field environments.


title: An exploratory view of water quality and sediment characteristics in Musa Estuary, the Persian Gulf

abstract: The results of mentioned physicochemical and biological analysis show that in more polluted sampling points, the number of resistance benthos (Bivalves) are more than other points of the study area. Heavy metal concentration in sediment is between 10 to 1000 times more than water samples especially for Arsenic, Mercury, and Vanadium. Dissolved Oxygen is acceptable all over the study area (more than 4 mg/L) based on environmental quality standards for dissolved oxygen which present the range of 2-9 mg/l for sensitive aquatic lives to less sensitive ones. However the presence of chemicals causes increase in the value of Chemical Oxygen Demand (around 500 to 1000 mg/L) and consequent decrease in water purification ability Nutrient uptake capacity which shown as the ratio of Nitrogen to phosphorus is appropriate (around 6) in parts of Musa estuary with more than 3 meters depth . The study area is a semi-closed tide-dominant water body with shallow depth, in comparison to the oceans. A 2-D modelling is applied to simulate current speed and direction. The results of this modelling shows except in the creeks, we have relatively suitable circulation of water in the estuary and current speed may exceed 0.8 m/s. It means Musa estuary is not a dead zone and in some parts which have enough water circulation and nutrient uptake for aquatic life, but huge amount of heavy metal (especially Mercury) in the sediment and water samples, indicates that all aquatic bodies overwhelmed by dangerous heavy metals.


title: Octagonal continued fraction and diagonal changes

abstract: In this short note we show that the octagon Farey map introduced by Smillie and Ulcigrai is an acceleration of the diagonal changes algorithm introduced by Delecroix and Ulcigrai.


title: Nifty Web Apps: Build a Web App for Any Text-Based Programming Assignment

abstract: While many students now interact with web apps across a variety of smart devices, the vast majority of our Nifty Assignments still present traditional user interfaces such as console input/output and desktop GUI. In this tutorial session, participants will learn to build simple web apps for programming assignments that execute student-written code to dynamically respond to user interactions resulting in a more modern app experience. Our approach requires up to 75% less code than similar desktop GUI apps while requiring few (if any) modifications to existing assignments. Instructors and students alike can run and modify these web apps on their own computers or deploy their apps online for access from any smart device at no cost. The tutorial presents examples from CS1 and CS2 courses in Python and Java, but the ideas apply generally.


title: BCS-BEC Crossover Effects and Pseudogap in Neutron Matter

abstract: Due to the large neutron-neutron scattering length, dilute neutron matter resembles the unitary Fermi gas, which lies half-way in the crossover from the BCS phase of weakly coupled Cooper pairs to the Bose-Einstein condensate of dimers. We discuss crossover effects in analogy with the T-matrix theory used in the physics of ultracold atoms, which we generalize to the case of a non-separable finite-range interaction. A problem of the standard Nozieres-Schmitt-Rink approach and different ways to solve it are discussed. It is shown that in the strong-coupling regime, the spectral function exhibits a pseudo-gap at temperatures above the critical temperature T_c. The effect of the correlated density on the density dependence of T_c is found to be rather weak, but a possibly important effect due to the reduced quasiparticle weight is identified.


title: The $U$ boson, interpolating between a generalized dark photon or dark $Z$, an axial boson and an axionlike particle

abstract: A light boson $U$ from an extra $U(1)$ interpolates between a generalized dark photon coupled to $Q,\ B$ and $L_i$ (or $B-L$), plus possibly dark matter, a dark $Z$ coupled to the $Z$ current, and one axially coupled to quarks and leptons. We identify the corresponding $U(1)_F$ symmetries, with $F= \gamma_Y Y+\gamma_B B+\gamma_{L_i}L_i+\gamma_A F_A+\gamma_{F'}F'\!+\gamma_d F_d$, $F_d$ acting in a dark sector and $F'$ on possible semi-inert BEH doublets uncoupled to quarks and leptons. The $U$ current is obtained from the $U(1)_F $ and $Z$ currents, with a mixing determined by the spin-0 BEH fields. The charge $Q_U$ of chiral quarks and leptons is a combination of $Q,\ B,\ L_i$ and $T_{3L}$ with the axial $F_A$. It involves in general isovector and isoscalar axial terms, in the presence of two BEH doublets. A longitudinal $U$ with axial couplings has enhanced interactions, and behaves much as an axionlike particle. Its axial couplings $g_A$, usually restricted to $< 2\times 10^{-7} m_U$(MeV), lead to effective pseudoscalar ones $g_P= g_A\times 2m_{q,l}\,/m_U=2^{1/4}\, G_F^{1/2}\, m_{q,l} \ A_\pm $. $\,A_\pm$ is proportional to an invisibility parameter $r=\cos\theta_A$ induced by a singlet v.e.v., possibly large and allowing the $U$ to be very weakly interacting. This allows for a very small gauge coupling, expressed with two doublets and a singlet as $g"\!/4\simeq 2\times 10^{-6} m_U$(MeV) $ r/\sin 2\beta$. We discuss phenomenological implications for meson decays, neutrino interactions, atomic-physics parity violation, naturally suppressed $\pi^0\to\gamma U$ decays, etc.. The $U$ boson fits within the grand-unification framework, in symbiosis with a $SU(4)_{\rm es }$ electrostrong symmetry broken at the GUT scale, with $Q_U$ depending on $Q,\ B-L,\ F_A$ and $T_{3A}$ through three parameters $\gamma_Y,\ \gamma_A$ and $\eta$.


title: Learning Context-Free Languages with Nondeterministic Stack RNNs

abstract: We present a differentiable stack data structure that simultaneously and tractably encodes an exponential number of stack configurations, based on Lang's algorithm for simulating nondeterministic pushdown automata. We call the combination of this data structure with a recurrent neural network (RNN) controller a Nondeterministic Stack RNN. We compare our model against existing stack RNNs on various formal languages, demonstrating that our model converges more reliably to algorithmic behavior on deterministic tasks, and achieves lower cross-entropy on inherently nondeterministic tasks.


title: DFT insights into new B-containing 212 MAX phases: Hf2AB2 (A = In, Sn)

abstract: 212 MAX phase borides are new additions to the MAX phase family with enhanced thermo-mechanical properties. In this article, first-principles calculations have been carried out to explore the mechanical properties, Vickers hardness, elastic anisotropy, thermal properties, and optical properties of predicted thermodynamically stable MAX compounds Hf2AB2 (A = In, Sn). The structural properties are compared with the available data to assess the validity of calculations. The mechanical stability of Hf2AB2 (A = In, Sn) compounds is established with the help of the computed stiffness constants (Cij). The possible reason for enhanced mechanical properties and Vickers hardness of Hf2SnB2 is explained based on the analysis of bonding strength, followed by the electronic density of states. Higher mechanical strength and Vickers hardness of Hf2AB2 (A = In, Sn) compared to Hf2AC (A = In, Sn) are also indicated in the light of charge density mapping. The values of Pugh ratio, Poissons ratio and Cauchy pressure predict brittle character of the studied compounds. Besides, the anisotropic nature of the titled borides is investigated by 2D and 3D plots of elastic moduli along with some well established anisotropy indices. Thermal properties were investigated by calculating the Debye temperature, minimum thermal conductivity, Gr\"uneisen parameter, and melting temperature. The thermal properties of Hf2AB2 (A = In, Sn) are also superior to Hf2AC (A = In, Sn). The optical constants such as real and imaginary parts of the dielectric function, refractive index, extinction coefficient, absorption coefficient, photoconductivity, reflectivity, and loss function are investigated.


title: A Clustering-Based Method for Automatic Educational Video Recommendation Using Deep Face-Features of Lecturers

abstract: Discovering and accessing specific content within educational video bases is a challenging task, mainly because of the abundance of video content and its diversity. Recommender systems are often used to enhance the ability to find and select content. But, recommendation mechanisms, especially those based on textual information, exhibit some limitations, such as being error-prone to manually created keywords or due to imprecise speech recognition. This paper presents a method for generating educational video recommendation using deep face-features of lecturers without identifying them. More precisely, we use an unsupervised face clustering mechanism to create relations among the videos based on the lecturer's presence. Then, for a selected educational video taken as a reference, we recommend the ones where the presence of the same lecturers is detected. Moreover, we rank these recommended videos based on the amount of time the referenced lecturers were present. For this task, we achieved a mAP value of 99.165%.


title: Improved Complexity Bounds in Wasserstein Barycenter Problem

abstract: In this paper, we focus on computational aspects of the Wasserstein barycenter problem. We propose two algorithms to compute Wasserstein barycenters of $m$ discrete measures of size $n$ with accuracy $\e$. The first algorithm, based on mirror prox with a specific norm, meets the complexity of celebrated accelerated iterative Bregman projections (IBP), namely $\widetilde O(mn^2\sqrt n/\e)$, however, with no limitations in contrast to the (accelerated) IBP, which is numerically unstable under small regularization parameter. The second algorithm, based on area-convexity and dual extrapolation, improves the previously best-known convergence rates for the Wasserstein barycenter problem enjoying $\widetilde O(mn^2/\e)$ complexity.


title: Concurrent Alternating Least Squares for multiple simultaneous Canonical Polyadic Decompositions

abstract: Tensor decompositions, such as CANDECOMP/PARAFAC (CP), are widely used in a variety of applications, such as chemometrics, signal processing, and machine learning. A broadly used method for computing such decompositions relies on the Alternating Least Squares (ALS) algorithm. When the number of components is small, regardless of its implementation, ALS exhibits low arithmetic intensity, which severely hinders its performance and makes GPU offloading ineffective. We observe that, in practice, experts often have to compute multiple decompositions of the same tensor, each with a small number of components (typically fewer than 20), to ultimately find the best ones to use for the application at hand. In this paper, we illustrate how multiple decompositions of the same tensor can be fused together at the algorithmic level to increase the arithmetic intensity. Therefore, it becomes possible to make efficient use of GPUs for further speedups; at the same time the technique is compatible with many enhancements typically used in ALS, such as line search, extrapolation, and non-negativity constraints. We introduce the Concurrent ALS algorithm and library, which offers an interface to Matlab, and a mechanism to effectively deal with the issue that decompositions complete at different times. Experimental results on artificial and real datasets demonstrate a shorter time to completion due to increased arithmetic intensity.


title: A graph-theoretic approach to a conjecture of Dixon and Pressman

abstract: Given $n \times n$ matrices, $A_1, \dots, A_k$, consider the linear operator $L(A_1,\dots,A_k) \, \colon \; \operatorname{M}_n \to \operatorname{M}_n$ given by \[ L(A_1,\dots,A_k)(A_{k+1})= \sum_{\sigma\in S_{k+1}} \operatorname{sign}(\sigma) A_{\sigma(1)}A_{\sigma(2)} \cdots A_{\sigma(k+1)}. \] The Amitsur-Levitzki theorem asserts that $L(A_1, \ldots, A_k)$ is identically $0$ for every $k \geq 2n-1$. Dixon and Pressman conjectured that if $k$ is an even number between $2$ and $2n - 2$, then the kernel of $L(A_1, \ldots, A_k)$ is of dimension $k$ for $A_1, \ldots, A_k\in \operatorname{M}_n(\mathbb{R})$ in general position. We prove this conjecture using graph-theoretic techniques.


title: False Discovery Rate Computation: Illustrations and Modifications

abstract: False discovery rates (FDR) are an essential component of statistical inference, representing the propensity for an observed result to be mistaken. FDR estimates should accompany observed results to help the user contextualize the relevance and potential impact of findings. This paper introduces a new user-friendly R package for computing FDRs and adjusting p-values for FDR control. These tools respect the critical difference between the adjusted p-value and the estimated FDR for a particular finding, which are sometimes numerically identical but are often confused in practice. Newly augmented methods for estimating the null proportion of findings - an important part of the FDR estimation procedure - are proposed and evaluated. The package is broad, encompassing a variety of methods for FDR estimation and FDR control, and includes plotting functions for easy display of results. Through extensive illustrations, we strongly encourage wider reporting of false discovery rates for observed findings.


title: Triphoton production at hadron colliders in NNLO QCD

abstract: We present next-to-next-to-leading-order (NNLO) QCD corrections to the production of three isolated photons in hadronic collisions at the fully differential level. We employ qT subtraction within MATRIX and an efficient implementation of analytic two-loop amplitudes in the leading-colour approximation to achieve the first on-the-fly calculation for this process at NNLO accuracy. Numerical results are presented for proton-proton collisions at energies ranging from 7 TeV to 100 TeV. We find full agreement with the 8 TeV results of arXiv:1911.00479 and confirm that NNLO corrections are indispensable to describe ATLAS 8 TeV data. In addition, we demonstrate the significance of NNLO corrections for future precision studies of triphoton production at higher collision energies.


title: Reformulation of a transmission and reflection problems in terms of a quantum wave impedance function

abstract: On the base of a 1D Shr\"{o}dinger equation the non-linear first-order differential equation (Ricatti type) for a quantum wave impedance function was derived. The advantages of this approach were discussed and demonstrated for a case of a single rectangular barrier. Both the scattering and the bound states problem were reformulated in terms of a quantum wave impedance and its application for solving both these problems was considered. The expressions for a reflection and a transmission coefficient were found on the base of a quantum wave impedance approach.


title: Smooth Variational Graph Embeddings for Efficient Neural Architecture Search

abstract: Neural architecture search (NAS) has recently been addressed from various directions, including discrete, sampling-based methods and efficient differentiable approaches. While the former are notoriously expensive, the latter suffer from imposing strong constraints on the search space. Architecture optimization from a learned embedding space for example through graph neural network based variational autoencoders builds a middle ground and leverages advantages from both sides. Such approaches have recently shown good performance on several benchmarks. Yet, their stability and predictive power heavily depends on their capacity to reconstruct networks from the embedding space. In this paper, we propose a two-sided variational graph autoencoder, which allows to smoothly encode and accurately reconstruct neural architectures from various search spaces. We evaluate the proposed approach on neural architectures defined by the ENAS approach, the NAS-Bench-101 and the NAS-Bench-201 search space and show that our smooth embedding space allows to directly extrapolate the performance prediction to architectures outside the seen domain (e.g. with more operations). Thus, it facilitates to predict good network architectures even without expensive Bayesian optimization or reinforcement learning.


title: $L_1$-norm regularized $L_1$-norm best-fit line problem

abstract: This work develops a sparse and outlier-insensitive method to fit a one-dimensional subspace that can be used as a replacement for eigenvector methods such as principal component analysis (PCA). The method is insensitive to outlier observations by formulating procedures as optimization problems that seek the best-fit line according to the $\ell^1$ norm. It is also capable of producing sparse principal components by leveraging an additional penalty term induce sparsity. The algorithm has a worst-case time complexity of $O{(m^2n \log n)}$ and, under certain conditions, produces a globally optimal solution. An implementation of this algorithm in the parallel and heterogeneous environment NVIDIA CUDA is tested on synthetic and real world datasets and compared to existing methods. The results demonstrate the scalability and efficiency of the proposed approach.


title: The AGN fuelling/feedback cycle in nearby radio galaxies III. 3D relative orientations of radio jets and CO discs and their interaction

abstract: This is the third paper of a series exploring the multi-frequency properties of a sample of eleven nearby low-excitation radio galaxies (LERGs) in the southern sky. We are conducting an extensive study of different galaxy components (stars, dust, warm and cold gas, radio jets) with the aim of better understanding the AGN fuelling/feedback cycle in LERGs. Here we present new, deep, sub-kpc resolution Karl G. Jansky Very Large Array (JVLA) data for five sample sources at 10 GHz. Coupling these data with previously-acquired Atacama Large Millimeter/submillimeter Array (ALMA) CO(2-1) observations and measurements of comparable quality from the literature, we carry out for the first time a full 3D analysis of the relative orientations of jet and disc rotation axes in six FR I LERGs. This analysis shows (albeit with significant uncertainties) that the relative orientation angles span a wide range ($\approx$30$^{\circ}-60^{\circ}$). There is no case where both axes are accurately aligned and there is a marginally significant tendency for jets to avoid the disc plane. Our study also provides further evidence for the presence of a jet-CO disc interaction (already inferred from other observational indicators) in at least one source, NGC 3100. In this case, the limited extent of the radio jets, along with distortions in both the molecular gas and the jet components, suggest that the jets are young, interacting with the surrounding matter and rapidly decelerating.


title: Converging to a Desired Orientation in a Flock of Agents

abstract: This work concentrates on different aspects of the \textit{consensus problem}, when applying it to a swarm of flocking agents. We examine the possible influence an external agent, referred to as {\em influencing agent} has on the flock. We prove that even a single influencing agent with a \textit{Face Desired Orientation behaviour} that is injected into the flock is sufficient for guaranteeing desired consensus of the flock of agents which follow a Vicsek-inspired Model. We further show that in some cases this can be guaranteed also in dynamic environments.


title: A Series of Unfortunate Counterfactual Events: the Role of Time in Counterfactual Explanations

abstract: Counterfactual explanations are a prominent example of post-hoc interpretability methods in the explainable Artificial Intelligence research domain. They provide individuals with alternative scenarios and a set of recommendations to achieve a sought-after machine learning model outcome. Recently, the literature has identified desiderata of counterfactual explanations, such as feasibility, actionability and sparsity that should support their applicability in real-world contexts. However, we show that the literature has neglected the problem of the time dependency of counterfactual explanations. We argue that, due to their time dependency and because of the provision of recommendations, even feasible, actionable and sparse counterfactual explanations may not be appropriate in real-world applications. This is due to the possible emergence of what we call "unfortunate counterfactual events." These events may occur due to the retraining of machine learning models whose outcomes have to be explained via counterfactual explanation. Series of unfortunate counterfactual events frustrate the efforts of those individuals who successfully implemented the recommendations of counterfactual explanations. This negatively affects people's trust in the ability of institutions to provide machine learning-supported decisions consistently. We introduce an approach to address the problem of the emergence of unfortunate counterfactual events that makes use of histories of counterfactual explanations. In the final part of the paper we propose an ethical analysis of two distinct strategies to cope with the challenge of unfortunate counterfactual events. We show that they respond to an ethically responsible imperative to preserve the trustworthiness of credit lending organizations, the decision models they employ, and the social-economic function of credit lending.


title: The noncommutative fractional Fourier law in bounded and unbounded domains

abstract: Using the spectral theory on the $S$-spectrum it is possible to define the fractional powers of a large class of vector operators. This possibility leads to new fractional diffusion and evolution problems that are of particular interest for nonhomogeneous materials where the Fourier law is not simply the negative gradient operator but it is a nonconstant coefficients differential operator of the form $$ T=\sum_{\ell=1}^3e_\ell a_\ell(x)\partial_{x_\ell}, \ \ \ x=(x_1,x_2,x_3)\in \bar{\Omega}, $$ where, $\Omega$ can be either a bounded or an unbounded domain in $\mathbb{R}^3$ whose boundary $\partial\Omega$ is considered suitably regular, $\bar{\Omega}$ is the closure of $\Omega$ and $e_\ell$, for $\ell=1,2,3$ are the imaginary units of the quaternions $\mathbb{H}$. The operators $T_\ell:=a_\ell(x)\partial_{x_\ell}$, for $\ell=1,2,3$, are called the components of $T$ and $a_1$, $a_2$, $a_3: \bar{\Omega} \subset\mathbb{R}^3\to \mathbb{R}$ are the coefficients of $T$. In this paper we study the generation of the fractional powers of $T$, denoted by $P_{\alpha}(T)$ for $\alpha\in(0,1)$, when the operators $T_\ell$, for $\ell=1,2,3$ do not commute among themselves. To define the fractional powers $P_{\alpha}(T)$ of $T$ we have to consider the weak formulation of a suitable boundary value problem associated with the pseudo $S$-resolvent operator of $T$. In this paper we consider two different boundary conditions. If $\Omega$ is unbounded we consider Dirichlet boundary conditions. If $\Omega$ is bounded we consider the natural Robin-type boundary conditions associated with the generation of the fractional powers of $T$.


title: LaND: Learning to Navigate from Disengagements

abstract: Consistently testing autonomous mobile robots in real world scenarios is a necessary aspect of developing autonomous navigation systems. Each time the human safety monitor disengages the robot's autonomy system due to the robot performing an undesirable maneuver, the autonomy developers gain insight into how to improve the autonomy system. However, we believe that these disengagements not only show where the system fails, which is useful for troubleshooting, but also provide a direct learning signal by which the robot can learn to navigate. We present a reinforcement learning approach for learning to navigate from disengagements, or LaND. LaND learns a neural network model that predicts which actions lead to disengagements given the current sensory observation, and then at test time plans and executes actions that avoid disengagements. Our results demonstrate LaND can successfully learn to navigate in diverse, real world sidewalk environments, outperforming both imitation learning and reinforcement learning approaches. Videos, code, and other material are available on our website https://sites.google.com/view/sidewalk-learning


title: Robust Isometric Non-Rigid Structure-from-Motion

abstract: Non-Rigid Structure-from-Motion (NRSfM) reconstructs a deformable 3D object from the correspondences established between monocular 2D images. Current NRSfM methods lack statistical robustness, which is the ability to cope with correspondence errors.This prevents one to use automatically established correspondences, which are prone to errors, thereby strongly limiting the scope of NRSfM. We propose a three-step automatic pipeline to solve NRSfM robustly by exploiting isometry. Step 1 computes the optical flow from correspondences, step 2 reconstructs each 3D point's normal vector using multiple reference images and integrates them to form surfaces with the best reference and step 3 rejects the 3D points that break isometry in their local neighborhood. Importantly, each step is designed to discard or flag erroneous correspondences. Our contributions include the robustification of optical flow by warp estimation, new fast analytic solutions to local normal reconstruction and their robustification, and a new scale-independent measure of 3D local isometric coherence. Experimental results show that our robust NRSfM method consistently outperforms existing methods on both synthetic and real datasets.


title: A graph theoretical framework for the strong Gram classification of non-negative unit forms of Dynkin type A

abstract: In the context of signed line graphs, this article introduces a modified inflation technique to study strong Gram congruence of non-negative (integral quadratic) unit forms, and uses it to show that weak and strong Gram congruence coincide among positive unit forms of Dynkin type A. The concept of inverse of a quiver is also introduced, and is used to obtain and analyze the Coxeter matrix of non-negative unit forms of Dynkin type A. Connected principal unit forms of such type are also classified.


title: Recent results on hard and rare probes from ALICE

abstract: In high-energy hadronic collisions, hard parton scatterings with large momentum transfers are prerequisites for the formation of hard and rare probes. In heavy-ion collisions, these probes---final state particles related to the early hard-parton scatterings---serve as a powerful tool to explore the whole evolution of the medium including the quark--gluon plasma (QGP) stage. Hard probes in proton-proton (pp) collisions test perturbative quantum chromodynamics (pQCD) processes and hadronization, and provide a reference for the nuclear collision systems. High-multiplicity pp collisions provide a bridge to heavy-ion collisions, due to their large event activity. Heavy-flavor hadrons containing at least one charm or beauty quark belong to the hard and rare probes. Recently, ALICE has measured a broad palette of heavy-flavour baryons allowing to shed more light on charm fragmentation. Measurements of charmonium production as a function of the event multiplicity can provide insight into the interplay between charmonium-production processes and soft processes driving the multiplicity. ALICE has investigated the multiplicity dependence of J/$\psi$ production at midrapidity and the ratio of $\psi \rm (2S)$ to J/$\psi$ yields at forward rapidity. The charm and beauty cross sections can be constrained by fits of Monte-Carlo generators to the measured dielectron continuum.


title: Reshaping consumption habits by exploiting energy-related micro-moment recommendations: A case study

abstract: The environmental change and its effects, caused by human influences and natural ecological processes over the last decade, prove that it is now more prudent than ever to transition to more sustainable models of energy consumption behaviors. User energy consumption is inductively derived from the time-to-time standards of living that shape the user's everyday consumption habits. This work builds on the detection of repeated usage consumption patterns from consumption logs. It presents the structure and operation of an energy consumption reduction system, which employs a set of sensors, smart-meters and actuators in an office environment and targets specific user habits. Using our previous research findings on the value of energy-related micro-moment recommendations, the implemented system is an integrated solution that avoids unnecessary energy consumption. With the use of a messaging API, the system recommends to the user the proper energy saving action at the right moment and gradually shapes user's habits. The solution has been implemented on the Home Assistant open source platform, which allows the definition of automations for controlling the office equipment. Experimental evaluation with several scenarios shows that the system manages first to reduce energy consumption, and second, to trigger users' actions that could potentially urge them to more sustainable energy consumption habits.


title: Wavefunctions, electronic localization and bonding properties for correlated materials beyond the Kohn-Sham formalism

abstract: Many-body theories such as dynamical mean field theory (DMFT) have enabled the description of the electron exchange-correlation interactions that are missing in current density functional theory (DFT) calculations. However, there has been relatively little focus on the wavefunctions from these theories. We present the methodology of the newly developed Elk-TRIQS interface and how to calculate the DFT with DMFT (DFT+DMFT) wavefunctions, which can be used to calculate DFT+DMFT wavefunction dependent quantities. We illustrate this by calculating the electron localized function (ELF) in monolayer SrVO$_3$ and CaFe$_2$As$_2$, which provides a means of visualizing their chemical bonds. Monolayer SrVO$_3$ ELFs are sensitive to the charge redistribution between the DFT, one-shot DFT+DMFT and fully charge self-consistent DFT+DMFT calculations. In both tetragonal and collapsed tetragonal CaFe$_2$As$_2$ phases, the ELF changes weakly with correlation induced charge redistribution of the hybridized As-p and Fe-d states. Nonetheless, the interlayer As-As bond in the collapsed tetragonal structure is robust to the changes at and around the Fermi level.


title: Prevalence of oxygen defects in an in-plane anisotropic transition metal dichalcogenide

abstract: Atomic scale defects in semiconductors enable their technological applications and realization of novel quantum states. Using scanning tunneling microscopy and spectroscopy complemented by ab-initio calculations we determine the nature of defects in the anisotropic van der Waals layered semiconductor ReS$_2$. We demonstrate the in-plane anisotropy of the lattice by directly visualizing chains of rhenium atoms forming diamond-shaped clusters. Using scanning tunneling spectroscopy we measure the semiconducting gap in the density of states. We reveal the presence of lattice defects and by comparison of their topographic and spectroscopic signatures with ab initio calculations we determine their origin as oxygen atoms absorbed at lattice point defect sites. These results provide an atomic-scale view into the semiconducting transition metal dichalcogenides, paving the way toward understanding and engineering their properties.


title: Quantitative rapid and finite time stabilization of the heat equation

abstract: The null controllability of the heat equation is known for decades [19,23,30]. The finite time stabilizability of the one dimensional heat equation was proved by Coron--Nguy\^en [13], while the same question for high dimensional spaces remained widely open. Inspired by Coron--Tr\'elat [14] we find explicit stationary feedback laws that quantitatively exponentially stabilize the heat equation with decay rate $\lambda$ and $Ce^{C\sqrt{\lambda}}$ estimates, where Lebeau--Robbiano's spectral inequality [30] is naturally used. Then a piecewise controlling argument leads to null controllability with optimal cost $Ce^{C/T}$, as well as finite time stabilization.


title: Moir\'e patterns in graphene -- rhenium disulfide vertical heterostructures

abstract: Vertical stacking of atomically thin materials offers a large platform for realizing novel properties enabled by proximity effects and moir\'e patterns. Here we focus on mechanically assembled heterostructures of graphene and ReS$_2$, a van der Waals layered semiconductor. Using scanning tunneling microscopy and spectroscopy (STM/STS) we image the sharp edge between the two materials as well as areas of overlap. Locally resolved topographic images revealed the presence of a striped superpattern originating in the interlayer interactions between graphene's hexagonal structure and the triclinic, low in-plane symmetry of ReS$_2$. We compare the results with a theoretical model that estimates the shape and angle dependence of the moir\'e pattern between graphene and ReS$_2$. These results shed light on the complex interface phenomena between van der Waals materials with different lattice symmetries.


title: Cross-Correlation of Planck CMB Lensing with DESI-Like LRGs

abstract: Cross-correlations between the lensing of the cosmic microwave background (CMB) and other tracers of large-scale structure provide a unique way to reconstruct the growth of dark matter, break degeneracies between cosmology and galaxy physics, and test theories of modified gravity. We detect a cross-correlation between DESI-like luminous red galaxies (LRGs) selected from DECaLS imaging and CMB lensing maps reconstructed with the Planck satellite at a significance of $S/N = 27.2$ over scales $\ell_{\rm min} = 30$, $\ell_{\rm max} = 1000$. To correct for magnification bias, we determine the slope of the LRG cumulative magnitude function at the faint limit as $s = 0.999 \pm 0.015$, and find corresponding corrections on the order of a few percent for $C^{\kappa g}_{\ell}, C^{gg}_{\ell}$ across the scales of interest. We fit the large-scale galaxy bias at the effective redshift of the cross-correlation $z_{\rm eff} \approx 0.68$ using two different bias evolution agnostic models: a HaloFit times linear bias model where the bias evolution is folded into the clustering-based estimation of the redshift kernel, and a Lagrangian perturbation theory model of the clustering evaluated at $z_{\rm eff}$. We also determine the error on the bias from uncertainty in the redshift distribution; within this error, the two methods show excellent agreement with each other and with DESI survey expectations.


title: Adaptive Robust Quadratic Programs using Control Lyapunov and Barrier Functions

abstract: This paper presents adaptive robust quadratic program (QP) based control using control Lyapunov and barrier functions for nonlinear systems subject to time-varying and state-dependent uncertainties. An adaptive estimation law is proposed to estimate the pointwise value of the uncertainties with pre-computable estimation error bounds. The estimated uncertainty and the error bounds are then used to formulate a robust QP, which ensures that the actual uncertain system will not violate the safety constraints defined by the control barrier function. Additionally, the accuracy of the uncertainty estimation can be systematically improved by reducing the estimation sampling time, leading subsequently to reduced conservatism of the formulated robust QP. The proposed approach is validated in simulations on an adaptive cruise control problem and through comparisons with existing approaches.


title: Distributed Computing in a Pandemic: A Review of Technologies Available for Tackling COVID-19

abstract: The current COVID-19 global pandemic caused by the SARS-CoV-2 betacoronavirus has resulted in over a million deaths and is having a grave socio-economic impact, hence there is an urgency to find solutions to key research challenges. Much of this COVID-19 research depends on distributed computing. In this article, I review distributed architectures -- various types of clusters, grids and clouds -- that can be leveraged to perform these tasks at scale, at high-throughput, with a high degree of parallelism, and which can also be used to work collaboratively. High-performance computing (HPC) clusters will be used to carry out much of this work. Several bigdata processing tasks used in reducing the spread of SARS-CoV-2 require high-throughput approaches, and a variety of tools, which Hadoop and Spark offer, even using commodity hardware. Extremely large-scale COVID-19 research has also utilised some of the world's fastest supercomputers, such as IBM's SUMMIT -- for ensemble docking high-throughput screening against SARS-CoV-2 targets for drug-repurposing, and high-throughput gene analysis -- and Sentinel, an XPE-Cray based system used to explore natural products. Grid computing has facilitated the formation of the world's first Exascale grid computer. This has accelerated COVID-19 research in molecular dynamics simulations of SARS-CoV-2 spike protein interactions through massively-parallel computation and was performed with over 1 million volunteer computing devices using the Folding@home platform. Grids and clouds both can also be used for international collaboration by enabling access to important datasets and providing services that allow researchers to focus on research rather than on time-consuming data-management tasks.


title: Analytic charged BHs in $f(\mathcal{R})$ gravity

abstract: In this article, we seek exact charged spherically symmetric black holes (BHs) with considering $f(\mathcal{R})$ gravitational theory. These BHs are characterized by convolution and error functions. Those two functions depend on a constant of integration which is responsible to make such a solution deviate from the Einstein general relativity (GR). The error function which constitutes the charge potential of the Maxwell field depends on the constant of integration and when this constant is vanishing we can not reproduce the Reissner-Nordstr\"om BH in the lower order of $f(\mathcal{R})$. This means that we can not reproduce Reissner-Nordstr\"om BH in lower-order-curvature theory, i.e., in GR limit $f(\mathcal{R})=\mathcal{R}$, we can not get the well known charged BH. We study the physical properties of these BHs and show that it is asymptotically approached as a flat spacetime or approach AdS/dS spacetime. Also, we calculate the invariants of the BHS and show that the singularities are milder than those of BH's of GR. Additionally, we derive the stability condition through the use of geodesic deviation. Moreover, we study the thermodynamics of our BH and investigate the impact of the higher-order-curvature theory. Finally, we show that all the BHs are stable and have radial speed equal to one through the use of odd-type mode.


title: Mechanism Design of a Bio-inspired Armwing Mechanism for Mimicking Bat Flapping Gait

abstract: The objective of this work is to design and develop a bio-inspired soft and articulated armwing structure which will be an integral component of a morphing aerial co-bot, Aerobat. In our design, we draw inspiration from bats. Bat membranous wings possess unique functions that make them a good example to take inspiration from and transform current aerial drones. In contrast with other flying vertebrates, bats have an extremely articulated musculoskeletal system, key to their body impact-survivability and deliver an impressively adaptive and multimodal locomotion behavior. Bats exclusively use this capability with structural flexibility to generate the controlled force distribution on each wing membrane. The wing flexibility, complex wing kinematics, and fast muscle actuation allow these creatures to change the body configuration within a few tens of milliseconds. These characteristics are crucial to the unrivaled agility of bats and copying them can potentially transform the state-of-the-art aerial drone design.


title: Sparse network asymptotics for logistic regression

abstract: Consider a bipartite network where $N$ consumers choose to buy or not to buy $M$ different products. This paper considers the properties of the logistic regression of the $N\times M$ array of i-buys-j purchase decisions, $\left[Y_{ij}\right]_{1\leq i\leq N,1\leq j\leq M}$, onto known functions of consumer and product attributes under asymptotic sequences where (i) both $N$ and $M$ grow large and (ii) the average number of products purchased per consumer is finite in the limit. This latter assumption implies that the network of purchases is sparse: only a (very) small fraction of all possible purchases are actually made (concordant with many real-world settings). Under sparse network asymptotics, the first and last terms in an extended Hoeffding-type variance decomposition of the score of the logit composite log-likelihood are of equal order. In contrast, under dense network asymptotics, the last term is asymptotically negligible. Asymptotic normality of the logistic regression coefficients is shown using a martingale central limit theorem (CLT) for triangular arrays. Unlike in the dense case, the normality result derived here also holds under degeneracy of the network graphon. Relatedly, when there happens to be no dyadic dependence in the dataset in hand, it specializes to recently derived results on the behavior of logistic regression with rare events and iid data. Sparse network asymptotics may lead to better inference in practice since they suggest variance estimators which (i) incorporate additional sources of sampling variation and (ii) are valid under varying degrees of dyadic dependence.


title: Recursive Top-Down Production for Sentence Generation with Latent Trees

abstract: We model the recursive production property of context-free grammars for natural and synthetic languages. To this end, we present a dynamic programming algorithm that marginalises over latent binary tree structures with $N$ leaves, allowing us to compute the likelihood of a sequence of $N$ tokens under a latent tree model, which we maximise to train a recursive neural function. We demonstrate performance on two synthetic tasks: SCAN (Lake and Baroni, 2017), where it outperforms previous models on the LENGTH split, and English question formation (McCoy et al., 2020), where it performs comparably to decoders with the ground-truth tree structure. We also present experimental results on German-English translation on the Multi30k dataset (Elliott et al., 2016), and qualitatively analyse the induced tree structures our model learns for the SCAN tasks and the German-English translation task.


title: Algorithmic Frameworks for the Detection of High Density Anomalies

abstract: This study explores the concept of high-density anomalies. As opposed to the traditional concept of anomalies as isolated occurrences, high-density anomalies are deviant cases positioned in the most normal regions of the data space. Such anomalies are relevant for various practical use cases, such as misbehavior detection and data quality analysis. Effective methods for identifying them are particularly important when analyzing very large or noisy sets, for which traditional anomaly detection algorithms will return many false positives. In order to be able to identify high-density anomalies, this study introduces several non-parametric algorithmic frameworks for unsupervised detection. These frameworks are able to leverage existing underlying anomaly detection algorithms and offer different solutions for the balancing problem inherent in this detection task. The frameworks are evaluated with both synthetic and real-world datasets, and are compared with existing baseline algorithms for detecting traditional anomalies. The Iterative Partial Push (IPP) framework proves to yield the best detection results.


title: Uncertainty over Uncertainty: Investigating the Assumptions, Annotations, and Text Measurements of Economic Policy Uncertainty

abstract: Methods and applications are inextricably linked in science, and in particular in the domain of text-as-data. In this paper, we examine one such text-as-data application, an established economic index that measures economic policy uncertainty from keyword occurrences in news. This index, which is shown to correlate with firm investment, employment, and excess market returns, has had substantive impact in both the private sector and academia. Yet, as we revisit and extend the original authors' annotations and text measurements we find interesting text-as-data methodological research questions: (1) Are annotator disagreements a reflection of ambiguity in language? (2) Do alternative text measurements correlate with one another and with measures of external predictive validity? We find for this application (1) some annotator disagreements of economic policy uncertainty can be attributed to ambiguity in language, and (2) switching measurements from keyword-matching to supervised machine learning classifiers results in low correlation, a concerning implication for the validity of the index.


title: Logical peering for interdomain networking on testbeds

abstract: Research testbed fabrics have potential to support long-lived, evolving, interdomain experiments, including opt-in application traffic across multiple campuses and edge sites. We propose abstractions and security infrastructure to facilitate multi-domain networking, and a reusable controller toolkit (ExoPlex) for network service providers (NSPs) running in testbed-hosted virtual network slices. We demonstrate the idea on the ExoGENI testbed, which allows slices to interconnect and exchange traffic over peering links by mutual consent. Each ExoPlex NSP runs a peering controller that manages its interactions with its linked peers and controls the NSP's dataplane network via SDN. Our approach expresses policies for secure peering and routing in a declarative language---logical peering. The prototype uses logic rules to verify IP prefix ownership, filter and validate route advertisements, and implement user-specified policies for connectivity and path control in networks with multiple transit NSPs.


title: Density Property and Composition Operators on $H(b)$ Spaces of Finitely Connected Planar Domains

abstract: In this work, the density in $H(b)$ spaces of finitely connected planar domains and the boundedness of composition operators on these function spaces are studied. Density of the algebra $\mathcal{A}(D)$ is considered for both in the cases where the defining function $b$ is an extreme and non-extreme point of the unit ball of $H^\infty(D)$. In the last part boundedness of composition operators on $H(b)$ spaces is considered and as well as a generalization of the unit disk case is given, the boundedness of composition operators with generalized Blaschke symbols over finitely connected domains is characterized.


title: Modelling overlap functions for one-nucleon removal: role of the effective three-nucleon force

abstract: One-nucleon overlap functions, needed for nucleon-removal reaction calculations, are solutions of an inhomogeneous equation with the source term defined by the wave functions of the initial and final nuclear states and interaction between the removed nucleon with the rest. The source term approach (STA) allows the overlaps with correct asymptotic decrease to be modelled while using nuclear many-body functions calculated in minimal model spaces. By properly choosing the removed nucleon interaction the minimum-model-space STA can reproduce reduced values of spectroscopic factors extracted from nucleon-removal reactions and predicts isospin asymmetry in the spectroscopic factor reduction. It is well-known that model space truncation leads to the appearance of higher-order induced forces, with three-nucleon force being the most important. In this paper the role of such a force on the source term calculation is studied. Applications to one-nucleon removal from double-magic nuclei show that three-nucleon force improves the description of available phenomenological overlap functions and reduces isospin asymmetry in spectroscopic factors.


title: Tentative detection of HC5NH+ in TMC-1

abstract: Using the Yebes 40m radio telescope, we report the detection of a series of seven lines harmonically related with a rotational constant B0=1295.81581 +/- 0.00026 MHz and a distortion constant D0=27.3 +/- 0.5 Hz towards the cold dense cloud TMC-1. Ab initio calculations indicate that the best possible candidates are the cations HC5NH+ and NC4NH+. From a comparison between calculated and observed rotational constants and other arguments based on proton affinities and dipole moments, we conclude that the best candidate for a carrier of the observed lines is the protonated cyanodiacetylene cation, HC5NH+. The HC5N/HC5NH+ ratio derived in TMC-1 is 240, which is very similar to the HC3N/HC3NH+ ratio. Results are discussed in the framework of a chemical model for protonated molecules in cold dense clouds.


title: Hybrid calculation of the MSSM Higgs boson masses using the complex THDM as EFT

abstract: Recently, the Higgs boson masses in the Minimal Supersymmetric Standard Model (MSSM) and their mixing have been calculated using the complex Two-Higgs-DoubletModel (cTHDM) as an effective field theory (EFT) of the MSSM. Here, we discuss the implementation of this calculation, which we improve in several aspects, into the hybrid framework of FeynHiggs by combing the cTHDM-EFT calculation with the existing fixed-order calculation. This combination allows accurate predictions also in the intermediate regime where some SUSY particles are relatively light, some relatively heavy and some in between. Moreover, the implementation provides precise predictions for the Higgs decay rates and production cross-sections.


title: Gaussian Process (GP)-based Learning Control of Selective Laser Melting Process

abstract: Selective laser melting (SLM) is one of emerging processes for effective metal additive manufacturing. Due to complex heat exchange and material phase changes, it is challenging to accurately model the SLM dynamics and design robust control of SLM process. In this paper, we first present a data-driven Gaussian process based dynamic model for SLM process and then design a model predictive control to regulate the melt pool size. Physical and process constraints are considered in the controller design. The learning model and control design are tested and validated with high-fidelity finite element simulation. The comparison results with other control design demonstrate the efficacy of the control design.


title: PathoNet: Deep learning assisted evaluation of Ki-67 and tumor infiltrating lymphocytes (TILs) as prognostic factors in breast cancer; A large dataset and baseline

abstract: The nuclear protein Ki-67 and Tumor infiltrating lymphocytes (TILs) have been introduced as prognostic factors in predicting tumor progression and its treatment response. The value of the Ki-67 index and TILs in approach to heterogeneous tumors such as Breast cancer (BC), known as the most common cancer in women worldwide, has been highlighted in the literature. Due to the indeterminable and subjective nature of Ki-67 as well as TILs scoring, automated methods using machine learning, specifically approaches based on deep learning, have attracted attention. Yet, deep learning methods need considerable annotated data. In the absence of publicly available benchmarks for BC Ki-67 stained cell detection and further annotated classification of cells, we propose SHIDC-BC-Ki-67 as a dataset for the aforementioned purpose. We also introduce a novel pipeline and a backend, namely PathoNet for Ki-67 immunostained cell detection and classification and simultaneous determination of intratumoral TILs score. Further, we show that despite facing challenges, our proposed backend, PathoNet, outperforms the state of the art methods proposed to date in the harmonic mean measure.


title: From hard spheres to hard-core spins

abstract: A system of hard spheres exhibits physics that is controlled only by their density. This comes about because the interaction energy is either infinite or zero, so all allowed configurations have exactly the same energy. The low density phase is liquid, while the high density phase is crystalline, an example of "order by disorder" as it is driven purely by entropic considerations. Here we study a family of hard spin models, which we call hardcore spin models, where we replace the translational degrees of freedom of hard spheres with the orientational degrees of freedom of lattice spins. Their hardcore interaction serves analogously to divide configurations of the many spin system into allowed and disallowed sectors. We present detailed results on the square lattice in $d=2$ for a set of models with $\mathbb{Z}_n$ symmetry, which generalize Potts models, and their $U(1)$ limits, for ferromagnetic and antiferromagnetic senses of the interaction, which we refer to as exclusion and inclusion models. As the exclusion/inclusion angles are varied, we find a Kosterlitz-Thouless phase transition between a disordered phase and an ordered phase with quasi-long-ranged order, which is the form order by disorder takes in these systems. These results follow from a set of height representations, an ergodic cluster algorithm, and transfer matrix calculations.


title: Short-Term Solar Irradiance Forecasting Using Calibrated Probabilistic Models

abstract: Advancing probabilistic solar forecasting methods is essential to supporting the integration of solar energy into the electricity grid. In this work, we develop a variety of state-of-the-art probabilistic models for forecasting solar irradiance. We investigate the use of post-hoc calibration techniques for ensuring well-calibrated probabilistic predictions. We train and evaluate the models using public data from seven stations in the SURFRAD network, and demonstrate that the best model, NGBoost, achieves higher performance at an intra-hourly resolution than the best benchmark solar irradiance forecasting model across all stations. Further, we show that NGBoost with CRUDE post-hoc calibration achieves comparable performance to a numerical weather prediction model on hourly-resolution forecasting.


title: Spatially resolved CIII]$\lambda$1909 emission in Haro 11

abstract: The CIII]1909 (hereafter, C III]) line is the strongest ultraviolet emission line after Ly$\alpha$ and is therefore of interest to high redshift studies of star-forming (SF) galaxies near the epoch of reionization. It is thought that C III] emission is strongest in galaxies with subsolar metallicity and low mass, however, spectral observations of numerous such galaxies at high and low redshift produce inconclusive or even contradictory results. We present the first ever C III ] imaging, obtained with HST/STIS for the low-redshift SF galaxy Haro 11. Cluster parameters like stellar mass, dust fraction and dust attenuation, and ionization parameter, obtained through spectral energy distribution fitting, show no correlation with the CIII] equivalent width (EW), which may be due to a combination of the limitation of the models and the age-homogeneity of the cluster population. Comparing the ratio of C III] emission line flux from individual clusters to that of H$\alpha$, [O III], and [O II] we find that the clusters with the highest EW(C III]) can be reconciled only with Cloudy models with an extremely high C/O ratio of $\ge$ 1.4(C/O)$_\odot$ for an ionizing population of single stars, binary stars, or a mixture of binary stars and active galactic nuclei. Given the point-like nature of strong C III], the integrated total strength of EW(C III]) becomes dependent on the morphology of the galaxy, which would explain the large scatter in EW(C III]) strengths, observed in galaxies with otherwise similar SF properties, and of similarly low metallicity and stellar mass.


title: Unsupervised 3D Brain Anomaly Detection

abstract: Anomaly detection (AD) is the identification of data samples that do not fit a learned data distribution. As such, AD systems can help physicians to determine the presence, severity, and extension of a pathology. Deep generative models, such as Generative Adversarial Networks (GANs), can be exploited to capture anatomical variability. Consequently, any outlier (i.e., sample falling outside of the learned distribution) can be detected as an abnormality in an unsupervised fashion. By using this method, we can not only detect expected or known lesions, but we can even unveil previously unrecognized biomarkers. To the best of our knowledge, this study exemplifies the first AD approach that can efficiently handle volumetric data and detect 3D brain anomalies in one single model. Our proposal is a volumetric and high-detail extension of the 2D f-AnoGAN model obtained by combining a state-of-the-art 3D GAN with refinement training steps. In experiments using non-contrast computed tomography images from traumatic brain injury (TBI) patients, the model detects and localizes TBI abnormalities with an area under the ROC curve of ~75%. Moreover, we test the potential of the method for detecting other anomalies such as low quality images, preprocessing inaccuracies, artifacts, and even the presence of post-operative signs (such as a craniectomy or a brain shunt). The method has potential for rapidly labeling abnormalities in massive imaging datasets, as well as identifying new biomarkers.


title: Translation of Tschebotarow's "The Resolvent Problem"

abstract: This is an English translation of "The Resolvent Problem" by Tschebotarow/Chebotarev. In this article, Chebotarev summarizes the history of the resolvent problem from compass and ruler constructions to Klein and Hilbert' formlutions of the problems. He also describes his work on the problem up until that time. The ideas of resolvent degree and essential dimension arose out of the resolvent problems Chebotarev describes here.


title: The relationship between driving volatility in time to collision and crash injury severity in a naturalistic driving environment

abstract: As a key indicator of unsafe driving, driving volatility characterizes the variations in microscopic driving decisions. This study characterizes volatility in longitudinal and lateral driving decisions and examines the links between driving volatility in time to collision and crash injury severity. By using a unique real-world naturalistic driving database from the 2nd Strategic Highway Research Program (SHRP), a test set of 671 crash events featuring around 0.2 million temporal samples of real world driving are analyzed. Based on different driving performance measures, 16 different volatility indices are created. To explore the relationships between crash-injury severity outcomes and driving volatility, the volatility indices are then linked with individual crash events including information on crash severity, drivers' pre crash maneuvers and behaviors, secondary tasks and durations, and other factors. As driving volatility prior to crash involvement can have different components, an indepth analysis is conducted using the aggregate as well as segmented (based on time to collision) real world driving data. To account for the issues of observed and unobserved heterogeneity, fixed and random parameter logit models with heterogeneity in parameter means and variances are estimated. The empirical results offer important insights regarding how driving volatility in time to collision relates to crash severity outcomes. Overall, statistically significant positive correlations are found between the aggregate (as well as segmented) volatility measures and crash severity outcomes. The findings suggest that greater driving volatility (both in longitudinal and lateral direction) in time to collision increases the likelihood of police reportable or most severe crash events... ...


title: Compensated Convex Based Transforms for Image Processing and Shape Interrogation

abstract: This paper reviews some recent applications of the theory of the compensated convex transforms or of the proximity hull as developed by the authors to image processing and shape interrogation with special attention given to the Hausdorff stability and multiscale properties. The paper contains also numerical experiments that demonstrate the performance of our methods compared to the state-of-art ones.


title: Thermal boundaries in kinetic and hydrodynamic limits

abstract: We investigate how a thermal boundary, modelled by a Langevin dynamics, affect the macroscopic evolution of the energy at different space-time scales.


title: Dynamic Green's functions in discrete flexural systems

abstract: The paper presents an analysis of the dynamic behaviour of discrete flexural systems composed of Euler--Bernoulli beams. The canonical object of study is the discrete Green's function, from which information regarding the dynamic response of the lattice under point loading by forces and moments can be obtained. Special attention is devoted to the interaction between flexural and torsional waves in a square lattice of Euler--Bernoulli beams, which is shown to yield a range of novel effects, including extreme dynamic anisotropy, non-reciprocity, wave-guiding, filtering, and the ability to create localised defect modes, all without the need for additional resonant elements or interfaces. The analytical study is complimented by numerical computations and finite element simulations, both of which are used to illustrate the effects predicted. A general algorithm is provided for constructing Green's functions as well as defect modes. This algorithm allows the tuning of the lattice to produce pass bands, band gaps, resonant modes, wave-guides, and defect modes, over any desired frequency range.


title: Signs, Spin, SMEFT: Sum Rules at Dimension Six

abstract: We place theoretical constraints on the leading deviations to four-fermion standard model interactions. Invoking S-matrix analyticity and partial wave unitarity, we develop new dispersion relations that yield either spin-dependent sum rules on dimension-six fermionic operators or information about the amplitude's behavior at large momentum. The pattern of SMEFT inequalities we find for theories obeying certain large-momentum constraints enables a diagnosis of properties of emerging new physics. These relations form a bridge between new physics searches: discovery of flavor-violating $\tau$ decays at Belle II would motivate the search for flavor-conserving new physics below 25 TeV, as the results would provide definitive information about the UV.


title: Kozai Migration Naturally Explains the White Dwarf Planet WD1856b

abstract: The Jovian-sized object WD~1856~b transits a white dwarf (WD) in a compact $1.4$-day orbit. Unlikely to have endured stellar evolution in its current orbit, WD~1856~b is thought to have migrated from much wider separations. Because the WD is old, and a member of a well-characterized hierarchical multiple, the well-known Kozai mechanism provides an effective migration channel for WD~1856~b. Moreover, the lack of tides in the star allows us to directly connect the current semi-major axis to the pre-migration one, from which we can infer the initial conditions of the system. By further demanding that successful migrators survive all previous phases of stellar evolution, we are able to constrain the mass of WD~1856~b to be $\simeq0.7-3M_{\rm J}$ and its main sequence semi-major axis to be $\simeq 2-2.5$ au. These properties imply that WD~1856~b was born a typical gas giant. We further estimate the occurrence rate of Kozai-migrated planets around WDs to be ${\cal O}(10^{-3}{-}10^{-4})$, suggesting that WD~1856~b is the only one in the {\it TESS} sample, but implying ${\cal O}(10^2)$ future detections by LSST. In a sense, WD~1856~b was an ordinary Jovian planet that underwent an extraordinary dynamical history.


title: The Number Densities and Stellar Populations of Massive Galaxies at 3 < z < 6: A Diverse, Rapidly Forming Population in the Early Universe

abstract: We present the census of massive (log(M$_{*}$/M$_{\odot}$)$\geq 11$) galaxies at $3<z<6$ identified over the COSMOS/UltraVISTA Ultra-Deep field stripes: consisting of $\approx100$ and $\approx20$ high-confidence candidates at $3<z<4$ and $4<z<6$, respectively. The $3<z<4$ population is comprised of post-starburst, UV star-forming and dusty-star forming galaxies in roughly equal fractions, while UV-star-forming galaxies dominate at $4<z<6$ . We account for various sources of biases in SED modelling, finding that the treatment of emission line contamination is essential for understanding the number densities and mass growth histories of massive galaxies at $z>3$. The significant increase in observed number densities at $z\sim4$ ($>\times$ 5 in $\lesssim600$ Myrs) implies that this is the epoch at which log(M$_{*}$/M$_{\odot}$)$\geq 11$ galaxies emerge in significant numbers, with stellar ages ($\approx500-900$ Myrs) indicating rapid formation epochs as early as $z\sim7$. Leveraging ancillary multi-wavelength datasets, we perform panchromatic SED modelling to constrain the total star-formation activity of the sample. The star-formation activity of the sample is generally consistent with being on the star-formation main sequence at the considered redshifts, with $\approx15-25\%$ of the population showing evidence of suppressed star-formation rates, indicating that quenching mechanisms are already at play by $z\sim4$. We stack available HST imaging, confirming their compact nature ($r_{e}\lesssim2.2$ kpc), consistent with expected sizes of high-$z$ star-forming galaxies. Finally, we discuss how our results are in-line with the early formation epochs and short formation timescales inferred from the fossil records of the most massive galaxies in the Universe.


title: Transverse spin in the light-ray OPE

abstract: We study a product of null-integrated local operators $\mathcal{O}_1$ and $\mathcal{O}_2$ on the same null plane in a CFT. Such null-integrated operators transform like primaries in a fictitious $d-2$ dimensional CFT in the directions transverse to the null integrals. We give a complete description of the OPE in these transverse directions. The terms with low transverse spin are light-ray operators with spin $J_1+J_2-1$. The terms with higher transverse spin are primary descendants of light-ray operators with higher spins $J_1+J_2-1+n$, constructed using special conformally-invariant differential operators that appear precisely in the kinematics of the light-ray OPE. As an example, the OPE between average null energy operators contains light-ray operators with spin $3$ (as described by Hofman and Maldacena), but also novel terms with spin $5,7,9,$ etc.. These new terms are important for describing energy two-point correlators in non-rotationally-symmetric states, and for computing multi-point energy correlators. We check our formulas in a non-rotationally-symmetric energy correlator in $\mathcal{N}=4$ SYM, finding perfect agreement.


title: The Scalar Chemical Potential in Cosmological Collider Physics

abstract: Non-analyticity in co-moving momenta within the non-Gaussian bispectrum is a distinctive sign of on-shell particle production during inflation, presenting a unique opportunity for the "direct detection" of particles with masses as large as the inflationary Hubble scale ($H$). However, the strength of such non-analyticity ordinarily drops exponentially by a Boltzmann-like factor as masses exceed $H$. In this paper, we study an exception provided by a dimension-5 derivative coupling of the inflaton to heavy-particle currents, applying it specifically to the case of two real scalars. The operator has a "chemical potential" form, which harnesses the large kinetic energy scale of the inflaton, $\dot{\phi}_{0}^{1/2} \approx 60H$, to act as an efficient source of scalar particle production. Derivative couplings of inflaton ensure radiative stability of the slow-roll potential, which in turn maintains (approximate) scale-invariance of the inflationary correlations. We show that a signal not suffering Boltzmann suppression can be obtained in the bispectrum with strength $f_{\mathrm{NL}} \sim \mathcal{O}(0.01-10)$ for an extended range of scalar masses, $M \lesssim \dot{\phi}_{0}^{1/2}$, potentially as high as $10^{15}$ GeV, within the sensitivity of upcoming LSS and more futuristic 21-cm experiments. The mechanism does not invoke any particular fine-tuning of parameters or breakdown of perturbation-theoretic control. The leading contribution appears at tree-level, which makes the calculation analytically tractable and removes the loop-suppression as compared to earlier chemical potential studies of non-zero spins. The steady particle production allows us to infer the effective mass of the heavy particles and the chemical potential from the variation in bispectrum oscillations as a function of co-moving momenta. Our analysis sets the stage for generalization to heavy bosons with non-zero spin.


title: Stellar Collapse Diversity and the Diffuse Supernova Neutrino Background

abstract: The diffuse cosmic supernova neutrino background (DSNB) is observational target of the gadolinium-loaded Super-Kamiokande (SK) detector and the forthcoming JUNO and Hyper-Kamiokande detectors. Current predictions are hampered by our still incomplete understanding of the supernova (SN) explosion mechanism and of the neutron star (NS) equation of state and maximum mass. In our comprehensive study we revisit this problem on grounds of the landscapes of successful and failed SN explosions obtained by Sukhbold et al. and Ertl et al. with parametrized one-dimensional neutrino engines for large sets of single-star and helium-star progenitors, with the latter serving as proxy of binary evolution effects. Besides considering engines of different strengths, leading to different fractions of failed SNe with black-hole (BH) formation, we also vary the NS mass limit, the spectral shape of the neutrino emission, and include contributions from poorly understood alternative NS-formation channels such as accretion-induced or merger-induced collapse events. Since the neutrino signals of our large model sets are approximate, we calibrate the associated degrees of freedom by using state-of-the-art simulations of proto-neutron star cooling. Our predictions are higher than other recent ones because of a large fraction of failed SNe with long delay to BH formation. Our best-guess model predicts a DSNB electron-antineutrino-flux of 28.8^{+24.6}_{-10.9} cm^{-2}s^{-1} with 6.0^{+5.1}_{-2.1} cm^{-2}s^{-1} in the favorable measurement interval of [10,30] MeV, and 1.3^{+1.1}_{-0.4} cm^{-2}s^{-1} with electron-antineutrino energies > 17.3 MeV, which is roughly a factor of two below the current SK limit. The uncertainty range is dominated by the still insufficiently constrained cosmic rate of stellar core-collapse events.