title: Testing effects of Lorentz invariance violation in the propagation of astroparticles with the Pierre Auger Observatory

abstract: Lorentz invariance violation (LIV) is often described by dispersion relations of the form $E_i^2=m_i^2+p_i^2+\delta_{i,n} E^{2+n}$ with delta different based on particle type $i$, with energy $E$, momentum $p$ and rest mass $m$. Kinematics and energy thresholds of interactions are modified once the LIV terms become comparable to the squared masses of the particles involved. Thus, the strongest constraints on the LIV coefficients $\delta_{i,n}$ tend to come from the highest energies. At sufficiently high energies, photons produced by cosmic ray interactions as they propagate through the Universe could be subluminal and unattenuated over cosmological distances. Cosmic ray interactions can also be modified and lead to detectable fingerprints in the energy spectrum and mass composition observed on Earth. The data collected at the Pierre Auger Observatory are therefore possibly sensitive to both the electromagnetic and hadronic sectors of LIV. In this article, we explore these two sectors by comparing the energy spectrum and the composition of cosmic rays and the upper limits on the photon flux from the Pierre Auger Observatory with simulations including LIV. Constraints on LIV parameters depend strongly on the mass composition of cosmic rays at the highest energies. For the electromagnetic sector, while no constraints can be obtained in the absence of protons beyond $10^{19}$ eV, we obtain $\delta_{\gamma,0} > -10^{-21}$, $\delta_{\gamma,1} > -10^{-40}$ eV$^{-1}$ and $\delta_{\gamma,2} > -10^{-58}$ eV$^{-2}$ in the case of a subdominant proton component up to $10^{20}$ eV. For the hadronic sector, we study the best description of the data as a function of LIV coefficients and we derive constraints in the hadronic sector such as $\delta_{\mathrm{had},0} < 10^{-19}$, $\delta_{\mathrm{had},1} < 10^{-38}$ eV$^{-1}$ and $\delta_{\mathrm{had},2}< 10^{-57}$ eV$^{-2}$ at 5$\sigma$ CL.


title: Mean-square-error-based secondary source placement in sound field synthesis with prior information on desired field

abstract: A method of optimizing secondary source placement in sound field synthesis is proposed. Such an optimization method will be useful when the allowable placement region and available number of loudspeakers are limited. We formulate a mean-square-error-based cost function, incorporating the statistical properties of possible desired sound fields, for general linear-least-squares-based sound field synthesis methods, including pressure matching and (weighted) mode matching, whereas most of the current methods are applicable only to the pressure-matching method. An efficient greedy algorithm for minimizing the proposed cost function is also derived. Numerical experiments indicated that a high reproduction accuracy can be achieved by the placement optimized by the proposed method compared with the empirically used regular placement.


title: On the Value of ML Models

abstract: We argue that, when establishing and benchmarking Machine Learning (ML) models, the research community should favour evaluation metrics that better capture the value delivered by their model in practical applications. For a specific class of use cases -- selective classification -- we show that not only can it be simple enough to do, but that it has import consequences and provides insights what to look for in a ``good'' ML model.


title: Keyphrase Generation Beyond the Boundaries of Title and Abstract

abstract: Keyphrase generation aims at generating important phrases (keyphrases) that best describe a given document. In scholarly domains, current approaches have largely used only the title and abstract of the articles to generate keyphrases. In this paper, we comprehensively explore whether the integration of additional information from the full text of a given article or from semantically similar articles can be helpful for a neural keyphrase generation model or not. We discover that adding sentences from the full text, particularly in the form of the extractive summary of the article can significantly improve the generation of both types of keyphrases that are either present or absent from the text. Experimental results with three widely used models for keyphrase generation along with one of the latest transformer models suitable for longer documents, Longformer Encoder-Decoder (LED) validate the observation. We also present a new large-scale scholarly dataset FullTextKP for keyphrase generation. Unlike prior large-scale datasets, FullTextKP includes the full text of the articles along with the title and abstract. We release the source code at https://github.com/kgarg8/FullTextKP.


title: Identifying Chern numbers of superconductors from local measurements

abstract: Fascination in topological materials originates from their remarkable response properties and exotic quasiparticles which can be utilized in quantum technologies. In particular, large-scale efforts are currently focused on realizing topological superconductors and their Majorana excitations. However, determining the topological nature of superconductors with current experimental probes is an outstanding challenge. This shortcoming has become increasingly pressing due to rapidly developing designer platforms which are theorized to display very rich topology and are better accessed by local probes rather than transport experiments. We introduce a robust machine-learning protocol for classifying the topological states of two-dimensional (2D) chiral superconductors and insulators from local density of states (LDOS) data. Since the LDOS can be measured with standard experimental techniques, our protocol contributes to overcoming the almost three decades standing problem of identifying the topological phase of 2D superconductors with broken time-reversal symmetry.


title: On the detection of a cosmic dawn signal in the radio background

abstract: The astrophysics of cosmic dawn, when star formation commenced in the first collapsed objects, is predicted to be revealed as spectral and spatial signatures in the cosmic radio background at long wavelengths. The sky-averaged redshifted 21-cm absorption line of neutral hydrogen is a probe of cosmic dawn. The line profile is determined by the evolving thermal state of the gas, radiation background, Lyman-$\alpha$ radiation from stars scattering off cold primordial gas and the relative populations of the hyperfine spin levels in neutral hydrogen atoms. We report a radiometer measurement of the spectrum of the radio sky in the 55--85~MHz band, which shows that the profile found by Bowman et al. in data taken with the Experiment to Detect the Global Epoch of Reionization Signature (EDGES) low-band instrument is not of astrophysical origin; their best-fitting profile is rejected with 95.3\% confidence. The profile was interpreted to be a signature of cosmic dawn; however, its amplitude was substantially higher than that predicted by standard cosmological models. Explanations for the amplitude of the profile included non-standard cosmology, additional mechanisms for cooling the baryons, perhaps via interactions with millicharged dark matter and an excess radio background at redshifts beyond 17. Our non-detection bears out earlier concerns and suggests that the profile found by Bowman et al. is not evidence for new astrophysics or non-standard cosmology.


title: CP-violation in oscillations of three neutrino generations: the case of degenerate masses

abstract: CP-symmetry violates in oscillations of neutrinos in vacuum in presence of three generations. In this case there is an inherent complex phase in PMNS-matrix. When any two neutrino masses degenerate the CP-violation vanishes. In the present article different mechanisms of CP-symmetry recovery are demonstrated in case of two degenerate neutrino mass-states.


title: Explanation Container in Case-Based Biomedical Question-Answering

abstract: The National Center for Advancing Translational Sciences(NCATS) Biomedical Data Translator (Translator) aims to attenuate problems faced by translational scientists. Translator is a multi-agent architecture consisting of six autonomous relay agents (ARAs) and eight knowledge providers (KPs). In this paper, we present the design of the Explanatory Agent (xARA), a case-based ARA that answers biomedical queries by accessing multiple KPs, ranking results, and explaining the ranking of results. The Explanatory Agent is designed with five knowledge containers that include the four original knowledge containers and one additional container for explanation - the Explanation Container. The Explanation Container is case-based and designed with its own knowledge containers.


title: Gromov Hyperbolicity, Geodesic Defect, and Apparent Pairs in Vietoris-Rips Filtrations

abstract: Motivated by computational aspects of persistent homology for Vietoris-Rips filtrations, we generalize a result of Eliyahu Rips on the contractibility of Vietoris-Rips complexes of geodesic spaces for a suitable parameter depending on the hyperbolicity of the space. We consider the notion of geodesic defect to extend this result to general metric spaces in a way that is also compatible with the filtration. We further show that for finite tree metrics the Vietoris-Rips complexes collapse to their corresponding subforests. We relate our result to modern computational methods by showing that these collapses are induced by the apparent pairs gradient, which is used as an algorithmic optimization in Ripser, explaining its particularly strong performance on tree-like metric data.


title: GCNDepth: Self-supervised Monocular Depth Estimation based on Graph Convolutional Network

abstract: Depth estimation is a challenging task of 3D reconstruction to enhance the accuracy sensing of environment awareness. This work brings a new solution with a set of improvements, which increase the quantitative and qualitative understanding of depth maps compared to existing methods. Recently, a convolutional neural network (CNN) has demonstrated its extraordinary ability in estimating depth maps from monocular videos. However, traditional CNN does not support topological structure and they can work only on regular image regions with determined size and weights. On the other hand, graph convolutional networks (GCN) can handle the convolution on non-Euclidean data and it can be applied to irregular image regions within a topological structure. Therefore, in this work in order to preserve object geometric appearances and distributions, we aim at exploiting GCN for a self-supervised depth estimation model. Our model consists of two parallel auto-encoder networks: the first is an auto-encoder that will depend on ResNet-50 and extract the feature from the input image and on multi-scale GCN to estimate the depth map. In turn, the second network will be used to estimate the ego-motion vector (i.e., 3D pose) between two consecutive frames based on ResNet-18. Both the estimated 3D pose and depth map will be used for constructing a target image. A combination of loss functions related to photometric, projection, and smoothness is used to cope with bad depth prediction and preserve the discontinuities of the objects. In particular, our method provided comparable and promising results with a high prediction accuracy of 89% on the publicly KITTI and Make3D datasets along with a reduction of 40% in the number of trainable parameters compared to the state of the art solutions. The source code is publicly available at https://github.com/ArminMasoumian/GCNDepth.git


title: A sea surface-based drag model for Large Eddy Simulation of wind-wave interaction

abstract: Monin-Obukhov similarity theory (MOST) is a well-tested approach for specifying the fluxes when the roughness surfaces are homogeneous. For flow over waves (inhomogeneous surfaces), phase-averaged roughness length scales are often prescribed through models based on the wave characteristics and the wind speed. However, such approaches lack generalizability over different wave ages and steepnesses due to the reliance on model coefficients tuned to specific datasets. In this paper, a sea surface-based hydrodynamic drag model applicable to moving surfaces is developed to model the pressure-based surface drag felt by the wind due to the waves. The model is based on the surface gradient approach of \citet{Anderson2010} applicable to stationary obstacles and extended here to the wind-wave problem. The wave drag model proposed specifies the hydrodynamic force based on the incoming momentum flux, wave phase speed, and the surface frontal area. The drag coefficient associated with the wind-wave momentum exchange is determined based on the wave steepness. The wave drag model is used to simulate turbulent airflow above a monochromatic wave train with different wave ages and wave steepnesses. The mean velocity profiles and model form stresses are validated with available lab-scale experimental data and show good agreement across a wide range of wave steepnesses and wave ages. The drag force is correlated with the wave surface gradient and out-of-phase with the wave height distribution by a factor $\pi/2$ for the sinusoidal wave-train considered. These results demonstrate that the current approach is sufficiently general over a wide parameter space compared to wave phase-averaged models with a minimal increase in computational cost.


title: Incompatible measurements in quantum information science

abstract: Some measurements in quantum mechanics disturb each other. This has puzzled physicists since the formulation of the theory, but only in recent decades has the incompatibility of measurements been analyzed in depth and detail, using the notion of joint measurability of generalized measurements. In this Colloquium joint measurability and incompatibility are reviewed from the perspective of quantum information science. The Colloquium starts by discussing the basic definitions and concepts. An overview on applications of incompatibility, such as in measurement uncertainty relations, the characterization of quantum correlations, or information processing tasks like quantum state discrimination, is then presented. Finally, emerging directions of research, such as a resource theory of incompatibility as well as other concepts to grasp the nature of measurements in quantum mechanics, are discussed.


title: Mass estimates of the SU(2) $0^{++}$ glueball from spectral methods

abstract: The estimation of the K\"all\'en-Lehmann spectral density from gauge invariant lattice QCD two point correlation functions is proposed, and explored via an inversion strategy based on Tikhonov regularisation. We test the method on a mesonic toy model, showing that our methodology is competitive with the traditional Maximum Entropy Method. As proof of concept the SU(2) glueball spectrum for the quantum numbers $J^{PC}=0^{++}$ is investigated, for various values of the lattice spacing, using the published data of arXiv:1910.07756. Our estimates for the ground state mass are in good agreement with the traditional approach, which is based on the large time exponential behaviour of the correlation functions. Furthermore, the spectral density also contains hints of excites states in the spectrum. Spectroscopic analysis of glueball two-point functions therefore provides a straightforward and insightful alternative to the traditional method based on the large time exponential behaviour of the correlation functions.


title: On using the complex step method for the approximation of Fr\'echet derivatives of matrix functions in automorphism groups

abstract: We show, that the Complex Step approximation to the Fr\'echet derivative of real matrix functions is applicable to the matrix sign, square root and polar mapping using iterative schemes. While this property was already discovered for the matrix sign using Newton's method, we extend the research to the family of Pad\'e iterations, that allows us to introduce iterative schemes for finding function and derivative values while approximately preserving automorphism group structure.


title: On Stability, Ancillary Services, Operation, and Security of Smart Inverters

abstract: This paper presents some recent trends in the research of grid-interactive inverters. Particularly, this paper focuses on stability, ancillary services, operation, and security of single and multiple inverters in the modern power grid. A grid-interactive inverter performs as a controllable interface between the distributed energy sources and the power grid. High-penetration of inverters in a power distribution system can create some technical challenges on the power quality, as well as voltage and frequency control of the system. In particular, a weak grid can lead to voltage oscillation and consequently instability. Moreover, the power grid is moving toward becoming a cyber-physical system in which smart inverters can exchange information for power marketing and economic dispatching. This puts the inverters at the risk of insecure operations. Hence, security enhancement has become another primary concern. Finally, the grid-interactive inverters are operated proportional power-sharing while operating together with many inverters. Recent research on coordinated operations is also discussed in this paper.


title: Locality properties of standard homogenization commutator

abstract: In the present work we study how the standard homogenization commutator, a random field that plays a central role in the theory of fluctuations, quantitatively decorrelates on large scales.


title: Matching next-to-leading-order and high-energy-resummed calculations of heavy-quarkonium-hadroproduction cross sections

abstract: The energy dependence of the total hadroproduction cross section of pseudo-scalar quarkonia is computed via matching Next-to-Leading Order (NLO) Collinear-Factorisation (CF) results with resummed higher-order corrections, proportional to $\alpha_s^{n}\ln^{n-1}(1/z)$, to the CF hard-scattering coefficient, where $z=M^2/\hat{s}$ with $M$ and $\hat{s}$ being the quarkonium mass and the partonic center-of-mass energy squared. The resummation is performed using High-Energy Factorisation (HEF) in the Doubly-Logarithmic (DL) approximation, which is a subset of the leading logarithmic $\ln (1/z)$ approximation. Doing so, one remains strictly consistent with the NLO and NNLO DGLAP evolution of the PDFs. By improving the treatment of the small-$z$ asymptotics of the CF coefficient function, the resummation cures the unphysical results of the NLO CF calculation. The matching is directly performed in the $z$-space and, for the first time, by using the Inverse-Error Weighting (InEW) matching procedure. As a by-product of the calculation, the NNLO term of the CF hard-scattering coefficient proportional to $\alpha_s^2\ln(1/z)$ is predicted from HEF.


title: Empirical atmosphere model in a mini flare during magnetic reconnection

abstract: A spatio-temporal analysis of IRIS spectra of MgII, CII, and SiIV ions allows us to study the dynamics and the stratification of the flare atmosphere along the line of sight during the magnetic reconnection phase at the jet base. Strong asymmetric MgII and CII line profiles with extended blue wings observed at the reconnection site are interpreted by the presence of two chromospheric temperature clouds: one explosive cloud with blueshifts at 290 km/s and one cloud with smaller Doppler shift (around 36 km/s). Simultaneously at the same location a mini flare was observed with strong emission in multi temperatures (AIA), in several spectral IRIS lines (e.g. Oiv and Siiv, Mgii), absorption of identified chromospheric lines in Siiv line profile, enhancement of the Balmer continuum and X-ray emission by FERMI/GBM. With the standard thick-target flare model we calculate the energy of non thermal electrons observed by FERMI and compare it to the energy radiated by the Balmer continuum emission. We show that the low energy input by non thermal electrons above 20 keV was still sufficient to produce the excess of Balmer continuum.


title: XY model with competing higher-order interactions

abstract: We study effects of competing pairwise higher-order interactions (HOI) with alternating signs and exponentially decreasing intensity on critical behavior of the XY model. It is found that critical properties of such a generalized model can be very different from the standard XY model and can strongly depend on whether the number of HOI terms is odd or even. Inclusion of any odd number of HOI terms results in two consecutive phase transitions to distinct ferromagnetic quasi-long-range order phases. Even number of HOI terms leads to two phase transitions only if the decay of the HOI intensities is relatively slow. Then the high-temperature transition to the ferromagnetic phase is followed by another transition to a peculiar competition-induced canted ferromagnetic phase. In the limit of an infinite number of HOI terms only one phase transition is confirmed, and under the conditions of fierce competition between the even and odd terms the transition temperature can be suppressed practically to zero.


title: Nonlinear pile-up separation with LSTM neural networks for cryogenic particle detectors

abstract: In high-background or calibration measurements with cryogenic particle detectors, a significant share of the exposure is lost due to pile-up of recoil events. We propose a method for the separation of pile-up events with an LSTM neural network and evaluate its performance on an exemplary data set. Despite a non-linear detector response function, we can reconstruct the ground truth of a severely distorted energy spectrum reasonably well.


title: Modeling the debonding process of osseointegrated implants due to coupled adhesion and friction

abstract: Cementless implants have become widely used for total hip replacement surgery. The long-term stability of these implants is achieved by bone growing around and into the porous surface of the implant, a process called osseointegration. However, debonding of the bone-implant interface can still occur due to aseptic implant loosening and insufficient osseointegration, which may have dramatic consequences. The aim of this work is to describe a new 3D finite element frictional contact formulation for the debonding of partially osseointegrated implants. The contact model is based on a modified Coulomb's friction law (Immel et al. 2020, Biomech. Model. Mechanobiol.) that takes into account the tangential debonding of the bone-implant interface. This model is extended in the direction normal to the bone-implant interface by considering a cohesive zone model, to account for adhesion phenomena in the normal direction and for adhesive friction of partially bonded interfaces. The model is applied to simulate the debonding of an acetabular cup implant. The influence of partial osseointegration and adhesive effects on the long-term stability of the implant is assessed. The influence of different patient- and implant-specific parameters such as the friction coefficient, the trabecular Young's modulus and the interference fit is also analyzed, in order to determine the optimal stability for different configurations. Furthermore, this work provides guidelines for future experimental and computational studies, that are necessary for further parameter calibration.


title: Small-noise approximation for Bayesian optimal experimental design with nuisance uncertainty

abstract: Calculating the expected information gain in optimal Bayesian experimental design typically relies on nested Monte Carlo sampling. When the model also contains nuisance parameters, which are parameters that contribute to the overall uncertainty of the system but are of no interest in the Bayesian design framework, this introduces a second inner loop. We propose and derive a small-noise approximation for this additional inner loop. The computational cost of our method can be further reduced by applying a Laplace approximation to the remaining inner loop. Thus, we present two methods, the small-noise Double-loop Monte Carlo and small-noise Monte Carlo Laplace methods. Moreover, we demonstrate that the total complexity of these two approaches remains comparable to the case without nuisance uncertainty. To assess the efficiency of these methods, we present three examples, and the last example includes the partial differential equation for the electrical impedance tomography experiment for composite laminate materials.


title: Strong-field Gravity Tests with the Double Pulsar

abstract: Continued observations of the Double Pulsar, PSR J0737-3039A/B, consisting of two radio pulsars (A and B) that orbit each other with a period of 2.45hr in a mildly eccentric (e=0.088) binary system, have led to large improvements in the measurement of relativistic effects in this system. With a 16-yr data span, the results enable precision tests of theories of gravity for strongly self-gravitating bodies and also reveal new relativistic effects that have been expected but are now observed for the first time. These include effects of light propagation in strong gravitational fields which are currently not testable by any other method. We observe retardation and aberrational light-bending that allow determination of the pulsar's spin direction. In total, we have detected seven post-Keplerian (PK) parameters, more than for any other binary pulsar. For some of these effects, the measurement precision is so high that for the first time we have to take higher-order contributions into account. These include contributions of A's effective mass loss (due to spin-down) to the observed orbital period decay, a relativistic deformation of the orbit, and effects of the equation of state of super-dense matter on the observed PK parameters via relativistic spin-orbit coupling. We discuss the implications of our findings, including those for the moment of inertia of neutron stars. We present the currently most precise test of general relativity's (GR's) quadrupolar description of gravitational waves, validating GR's prediction at a level of $1.3 \times 10^{-4}$ (95% conf.). We demonstrate the utility of the Double Pulsar for tests of alternative theories by focusing on two specific examples and discuss some implications for studies of the interstellar medium and models for the formation of the Double Pulsar. Finally, we provide context to other types of related experiments and prospects for the future.


title: Depth Uncertainty Networks for Active Learning

abstract: In active learning, the size and complexity of the training dataset changes over time. Simple models that are well specified by the amount of data available at the start of active learning might suffer from bias as more points are actively sampled. Flexible models that might be well suited to the full dataset can suffer from overfitting towards the start of active learning. We tackle this problem using Depth Uncertainty Networks (DUNs), a BNN variant in which the depth of the network, and thus its complexity, is inferred. We find that DUNs outperform other BNN variants on several active learning tasks. Importantly, we show that on the tasks in which DUNs perform best they present notably less overfitting than baselines.


title: Symplectomorphisms of some Weinstein 4-manifolds

abstract: Let M be a Weinstein four-manifold mirror to Y\D for (Y,D) a log Calabi--Yau surface; intuitively, this is typically the Milnor fibre of a smoothing of a cusp singularity. We introduce two families of symplectomorphisms of M: Lagrangian translations, which we prove are mirror to tensors with line bundles; and nodal slide recombinations, which we prove are mirror to automorphisms of (Y,D). The proof uses a detailed compatibility between the homological and SYZ view-points on mirror symmetry. Together with spherical twists, these symplectomorphisms are expected to generate all autoequivalences of the wrapped Fukaya category of M which are compactly supported in a categorical sense. A range of applications is given.


title: Excitonic fine structure of epitaxial Cd(Se,Te) on ZnTe type-II quantum dots

abstract: The structure of the ground state exciton of Cd(Se,Te) quantum dots embedded in ZnTe matrix is studied experimentally using photoluminescence spectroscopy and theoretically using ${\bf k}\cdot{\bf p}$ and configuration interaction methods. The experiments reveal a considerable reduction of fine-structure splitting energy of the exciton with increase of Se content in the dots. That effect is interpreted by theoretical calculations to originate due to the transition from spatially direct (type-I) to indirect (type-II) transition between electrons and holes in the dot induced by increase of Se. The trends predicted by the theory match those of the experimental results very well.The theory identifies that the main mechanism causing elevated fine-structure energy in particular in type-I dots is due to the multipole expansion of the exchange interaction. Moreover, the theory reveals that for Se contents in the dot $>0.3$, there exist also a {\bf peculiar type of confinement showing signatures of both type~I and type~II} and which exhibits extraordinary properties, such as almost purely light hole character of exciton and toroidal shape of hole states.


title: Universal gate operations on nuclear spin qubits in an optical tweezer array of $^{171}$Yb atoms

abstract: Neutral atom arrays are a rapidly developing platform for quantum science. In this work, we demonstrate a universal set of quantum gate operations on a new type of neutral atom qubit: a nuclear spin in an alkaline earth-like atom (AEA), $^{171}$Yb. We present techniques for cooling, trapping and imaging using a newly discovered magic trapping wavelength at $\lambda = 486.78$ nm. We implement qubit initialization, readout, and single-qubit gates, and observe long qubit lifetimes, $T_1 \approx 20$ s and $T^*_2 = 1.24(5)$ s, and a single-qubit operation fidelity $\mathcal{F}_{1Q} = 0.99959(6)$. We also demonstrate two-qubit entangling gates using the Rydberg blockade, as well as coherent control of these gate operations using light shifts on the Yb$^+$ ion core transition at 369 nm. These results are a significant step towards highly coherent quantum gates in AEA tweezer arrays.


title: Revisiting PS191 limits on sterile neutrinos

abstract: We perform Monte Carlo simulations of the sterile neutrino signal at the fixed target experiment PS191 operated on a proton beam of 19.2\,GeV at CERN in the eighties. We find that the strongest bounds the PS191 could obtain are significantly lower than what they published, and now are obsolete being surpassed by recent T2K, NA62, E949, TRIUMF and PIENU experiments.


title: Adding species to chemical reaction networks: preserving rank preserves nondegenerate behaviours

abstract: We show that adding new chemical species into the reactions of a chemical reaction network (CRN) in such a way that the rank of the network remains unchanged preserves its capacity for multiple nondegenerate equilibria and/or periodic orbits. One consequence is that any bounded nondegenerate behaviours which can occur in a CRN can occur in a CRN with bounded stoichiometric classes. The main result adds to a family of theorems which tell us which enlargements of a CRN preserve its capacity for nontrivial dynamical behaviours. It generalises some earlier claims, and complements similar claims involving the addition of reactions into CRNs. The result gives us information on how ignoring some chemical species, as is common in biochemical modelling, might affect the allowed dynamics in differential equation models of CRNs. We demonstrate the scope and limitations of the main theorem via several examples. These illustrate how we can use the main theorem to predict multistationarity and oscillation in CRNs enlarged with additional species; but also how the enlargements can introduce new behaviours such as additional periodic orbits and new bifurcations.


title: Accounting for survey design in Bayesian disaggregation of survey-based areal estimates of proportions: an application to the American Community Survey

abstract: Understanding the effects of social determinants of health on health outcomes requires data on characteristics of the neighborhoods in which subjects live. However, estimates of these characteristics are often aggregated over space and time in a fashion that diminishes their utility. Take, for example, estimates from the American Community Survey (ACS), a multi-year nationwide survey administered by the U.S. Census Bureau: estimates for small municipal areas are aggregated over 5-year periods, whereas 1-year estimates are only available for municipal areas with populations $>$65,000. Researchers may wish to use ACS estimates in studies of population health to characterize neighborhood-level exposures. However, 5-year estimates may not properly characterize temporal changes or align temporally with other data in the study, while the coarse spatial resolution of the 1-year estimates diminishes their utility in characterizing neighborhood exposure. To circumvent this issue, in this paper we propose a modeling framework to disaggregate estimates of proportions derived from sampling surveys which explicitly accounts for the survey design effect. We illustrate the utility of our model by applying it to the ACS data, generating estimates of poverty for the state of Michigan at fine spatio-temporal resolution.


title: Filament coalescence and hub structure in MonR2: Implications to massive star and cluster formation

abstract: Here we study the MonR2 star forming region, which has a rich network of filaments joining in a star cluster forming hub, aiming at understanding the hub structure and to examine the mass fraction residing in the hub and in the filaments, which is a key factor that influences massive star formation. We conducted a multi-scale, multi-component analysis of the Herschel column density maps (resolution of 18.2" or $\sim$0.07 pc at 830 pc) of the region using a newly developed algorithm "getsf" to identify the structural components, namely, extended cloud, filaments, and sources. We find that cascades of lower column density filaments coalesce to form higher density filaments eventually merging inside the hub (0.8 pc radius). As opposed to the previous view of the hub as a massive clump with $\sim$1 pc radius, we find it to be a network of short high-density filaments. The total mass reservoir in the MonR2 HFS (5 pc $\times$ 5 pc) is split between filaments (54%), extended cloud (37%) and sources (9%). The M/L of filaments increase from $\sim$ 10 Msun/pc at 1.5pc from the hub to $\sim$ 100 Msun/pc at its centre, while the number of filaments per annulus of 0.2pc width decreases from 20 to 2 in the same range. The observed radial column density structure of the HFS (filament component only) displays a power-law dependence of $N_{\mathrm{H}_2} \propto r^{-2.17}$ up to a radius of $\sim$2.5 pc from the central hub, resembling a global collapse of the HFS. We present a scenario where the HFS can be supported by magnetic fields which interact, merge and reorganize themselves as the filaments coalesce. In the new view of the hub as a network of high-density filaments, we suggest that only the stars located in the network can benefit from the longitudinal flows of gas to become massive, which may explain the reason for the formation of many low-mass stars in cluster centres.


title: Open or Sneaky? Fast or Slow? Light or Heavy?: Investigating Security Releases of Open Source Packages

abstract: Vulnerabilities in open source packages can be a security risk for the client projects that use these packages as dependencies. When a new vulnerability is discovered in a package, the package should quickly release a fix in a new version, referred to as security release in this study. The security release should be well-documented and require minimal migration effort to facilitate fast adoption by the client projects. However, to what extent the open source packages follow these recommendations is not known. The goal of this study is to aid software practitioners and researchers in understanding the current practice of releasing security fixes by open source packages and identifying areas for improvement through an empirical study of security releases. Specifically, in this paper, we study (1) the time lag between fix and release; (2) how security fixes are documented in the release notes; (3) code change characteristics (size and semantic versioning) of the release; and (4) the time lag between the release and an advisory publication on Snyk or NVD (two popular vulnerability databases) for security releases over a dataset of 4,377 security advisories across seven package ecosystems. We find that the median security release is available in under 4 days of the corresponding fix and contains 134 lines of code (LOC) change. Further, we find that 61.5% of the security releases come with a release note that documents the corresponding security fix. However, Snyk and NVD may take a median of 25 days (from the release) to publish an advisory for these security releases, possibly resulting in delayed notification to the client projects. Based on our findings, we make four recommendations for the package maintainers and the ecosystem administrators, such as using private fork for security fixes and standardizing the practice for announcing security releases.


title: Spontaneous anomalous Hall effect arising from an unconventional compensated magnetic phase in a semiconductor

abstract: The anomalous Hall effect, commonly observed in metallic magnets, has been established to originate from the time-reversal symmetry breaking by an internal macroscopic magnetization in ferromagnets or by a non-collinear magnetic order. Here we observe a spontaneous anomalous Hall signal in the absence of an external magnetic field in an epitaxial film of MnTe, which is a semiconductor with a collinear antiparallel magnetic ordering of Mn moments and a vanishing net magnetization. The anomalous Hall effect arises from an unconventional phase with strong time-reversal symmetry breaking and alternating spin polarization in real-space crystal structure and momentum-space electronic structure. The anisotropic crystal environment of magnetic Mn atoms due to the non-magnetic Te atoms is essential for establishing the unconventional phase and generating the anomalous Hall effect.


title: Automatic Multi-Class Cardiovascular Magnetic Resonance Image Quality Assessment using Unsupervised Domain Adaptation in Spatial and Frequency Domains

abstract: Population imaging studies rely upon good quality medical imagery before downstream image quantification. This study provides an automated approach to assess image quality from cardiovascular magnetic resonance (CMR) imaging at scale. We identify four common CMR imaging artefacts, including respiratory motion, cardiac motion, Gibbs ringing, and aliasing. The model can deal with images acquired in different views, including two, three, and four-chamber long-axis and short-axis cine CMR images. Two deep learning-based models in spatial and frequency domains are proposed. Besides recognising these artefacts, the proposed models are suitable to the common challenges of not having access to data labels. An unsupervised domain adaptation method and a Fourier-based convolutional neural network are proposed to overcome these challenges. We show that the proposed models reliably allow for CMR image quality assessment. The accuracies obtained for the spatial model in supervised and weakly supervised learning are 99.41+0.24 and 96.37+0.66 for the UK Biobank dataset, respectively. Using unsupervised domain adaptation can somewhat overcome the challenge of not having access to the data labels. The maximum achieved domain gap coverage in unsupervised domain adaptation is 16.86%. Domain adaptation can significantly improve a 5-class classification task and deal with considerable domain shift without data labels. Increasing the speed of training and testing can be achieved with the proposed model in the frequency domain. The frequency-domain model can achieve the same accuracy yet 1.548 times faster than the spatial model. This model can also be used directly on k-space data, and there is no need for image reconstruction.


title: Hedging Cryptocurrency Options

abstract: The cryptocurrency market is volatile, non-stationary and non-continuous. Together with liquid derivatives markets, this poses a unique opportunity to study risk management, especially the hedging of options, in a turbulent market. We study the hedge behaviour and effectiveness for the class of affine jump diffusion models and infinite activity Levy processes. First, market data is calibrated to stochastic volatility inspired (SVI)-implied volatility surfaces to price options. To cover a wide range of market dynamics, we generate Monte Carlo price paths using an SVCJ model (stochastic volatility with correlated jumps), a close-to-actual-market GARCH-filtered kernel density estimation as well as a historical backtest. In all three settings, options are dynamically hedged with Delta, Delta-Gamma, Delta-Vega and Minimum Variance strategies. Including a wide range of market models allows to understand the trade-off in the hedge performance between complete, but overly parsimonious models, and more complex, but incomplete models. The calibration results reveal a strong indication for stochastic volatility, low jump frequency and evidence of infinite activity. Short-dated options are less sensitive to volatility or Gamma hedges. For longer-dated options, tail risk is consistently reduced by multiple-instrument hedges, in particular by employing complete market models with stochastic volatility.


title: Qualitative properties of space-dependent SIR models with constant delay and their numerical solutions

abstract: In this article a space-dependent epidemic model equipped with a constant latency period is examined. We construct a delay partial integro-differential equation and show that its solution possesses some biologically reasonable features. We propose some numerical schemes and show that by choosing the time step to be sufficiently small the schemes preserve the qualitative properties of the original continuous model. Finally, some numerical experiments are presented that confirm the aforementioned theoretical results.


title: Persistent Animal Identification Leveraging Non-Visual Markers

abstract: Our objective is to locate and provide a unique identifier for each mouse in a cluttered home-cage environment through time, as a precursor to automated behaviour recognition for biological research. This is a very challenging problem due to (i) the lack of distinguishing visual features for each mouse, and (ii) the close confines of the scene with constant occlusion, making standard visual tracking approaches unusable. However, a coarse estimate of each mouse's location is available from a unique RFID implant, so there is the potential to optimally combine information from (weak) tracking with coarse information on identity. To achieve our objective, we make the following key contributions: (a) the formulation of the object identification problem as an assignment problem (solved using Integer Linear Programming), and (b) a novel probabilistic model of the affinity between tracklets and RFID data. The latter is a crucial part of the model, as it provides a principled probabilistic treatment of object detections given coarse localisation. Our approach achieves 77% accuracy on this animal identification problem, and is able to reject spurious detections when the animals are hidden.


title: Bento and the Art of Repeated Research

abstract: Bento provides a new approach to developing file systems, with safety and high-velocity development in mind. This is achieved by using Rust, a modern and memory-safe systems programming language, and by providing a framework to run a single file system implementation in kernel space with the VFS or in user space with FUSE. In this paper, the benchmarking experiments from the Bento paper are repeated. We fail to exactly reproduce the results of the Bento paper, but more or less find the same patterns albeit with more outlying results. Additionally we unsuccessfully run a standardized test suite, and expand the set of experiments with latency benchmarks and throughput benchmarks using a RAM block device. The latency benchmarks show that ext4 with journaling consistently outperforms Bento-fs and the RAM throughput benchmarks show no additional consistent performance pattern. During this experimentation, a set of 12 bugs was encountered and analyzed. We find that the ratio of memory related bugs is lower than other systems programming projects that use C as opposed to Rust, thus supporting the claims of the Bento framework.


title: Ex-post moral hazard and manipulation-proof contracts

abstract: We examine the trade-off between the provision of incentives to exert costly effort (ex-ante moral hazard) and the incentives needed to prevent the agent from manipulating the profit observed by the principal (ex-post moral hazard). Formally, we build a model of two-stage hidden actions where the agent can both influence the expected revenue of a business and manipulate its observed profit. We show that manipulation-proofness is sensitive to the interaction between the manipulation technology and the probability distribution of the stochastic output. The optimal contract is manipulation-proof whenever the manipulation technology is linear. However, a convex manipulation technology sometimes leads to contracts with manipulations in equilibrium. Whenever the distribution satisfies the monotone likelihood ratio property, we can always find a manipulation technology for which the optimal contract is not manipulation-proof.


title: Jet quenching in glasma

abstract: We discuss the transverse momentum broadening of hard probes traversing an evolving glasma, which is the earliest phase of the matter produced in relativistic heavy-ion collisions. The coefficient $\hat q$ is calculated using the Fokker-Planck equation, and an expansion in the proper time $\tau$ which is applied to describe the temporal evolution of the glasma. The correlators of the chromodynamic fields that determine the Fokker-Planck collision term, which in turn provides $\hat q$, are computed to fifth order in $\tau$. The momentum broadening is shown to rapidly grow in time and reach a magnitude of several ${\rm GeV^2/fm}$. We show that the transient pre-equilibrium phase provides a contribution to the energy loss of hard probes which is comparable to that of the long lasting, hydrodynamically evolving, equilibrium phase.


title: Induction of Large Magnetic Anisotropy Energy and Formation of Multiple Dirac States in SrIrO$_3$ Films: Role of correlation and spin-orbit coupling

abstract: The 5$d$ transition metal oxides, in particular iridates, host novel electronic and magnetic phases due to the interplay between onsite Coulomb repulsion ($U$) and spin-orbit coupling (SOC). The reduced dimensionality brings another degree of freedom to increase the functionality of these systems. By taking the example of ultrathin films of SrIrO$_3$,theoretically, we demonstrate that confinement led localization can introduce large magnetic anisotropy energy (MAE) in the range 2-7 meV/Ir which is one to two orders higher than that of the traditional MAE compounds formed out of transition metals and their multilayers. Furthermore, in the weak correlation limit, tailored terminations can yield multiple Dirac states across a large energy window of 2 eV around the Fermi energy which is a rare phenomena in correlated oxides and upon experimental realization it will give rise to unique transport properties with excitation and doping.


title: Post-Quantum Era in V2X Security: Convergence of Orchestration and Parallel Computation

abstract: Along with the potential emergence of quantum computing, safety and security of new and complex communication services such as automated driving need to be redefined in the post-quantum era. To ensure reliable, continuous and secure operation of these scenarios, quantum resistant security algorithms (QRSAs) that enable secure connectivity must be integrated into the network management and orchestration systems of mobile networks. This paper explores a roadmap study of post-quantum era convergence with cellular connectivity using the Service & Computation Orchestrator (SCO) framework for enhanced data security in radio access and backhaul transmission with a particular focus on Vehicle-to-Everything (V2X) services. Using NTRU as a QSRA, we have shown that the parallelization performance of Toom-Cook and Karatsuba computation methods can vary based on different CPU load conditions through extensive simulations and that the SCO framework can facilitate the selection of the most efficient computation for a given QRSA. Finally, we discuss the evaluation results, identify the current standardization efforts, and possible directions for the coexistence of post-quantum and mobile network connectivity through a SCO framework that leverages parallel computing.


title: Envelope theorem and discontinuous optimisation: the case of positioning choice problems

abstract: This article examines differentiability properties of the value function of positioning choice problems, a class of optimisation problems in finite-dimensional Euclidean spaces. We show that positioning choice problems' value function is always almost everywhere differentiable even when the objective function is discontinuous. To obtain this result we first show that the Dini superdifferential is always well-defined for the maxima of positioning choice problems. This last property allows to state first-order necessary conditions in terms of Dini supergradients. We then prove our main result, which is an ad-hoc envelope theorem for positioning choice problems. Lastly, after discussing necessity of some key assumptions, we conjecture that similar theorems might hold in other spaces as well.


title: A Fully Fiber-Integrated Ion Trap for Portable Optical Atomic Clocks

abstract: We present a novel, single-ion trap with integrated optical fibers directly embedded within the trap structure to deliver laser light as well as collect the ion's fluorescence. This eliminates the need for optical windows. We characterise the system's performance and measure signal-to-background ratios in the ion's fluorescence on the order of 50, which allows us to perform state readout with a fidelity over 99% in 600 $\mu$s. We test the system's resilience to thermal variations in the range between 22{\deg}C and 53{\deg}C, and the system's vibration resilience at 34 Hz and 300 Hz and find no effect on its performance. The combination of compactness and robustness of our fiber-coupled trap makes it well suited for applications in, as well as outside, research laboratory environments and in particular for highly compact portable optical atomic clocks. While our system is designed for trapping $^{40}$Ca$^{+}$ ions the fundamental design principles can be applied to other ion species.


title: Insurance design and arson-type risks

abstract: We design the insurance contract when the insurer faces arson-type risks. The optimal contract must be manipulation-proof. It is therefore continuous, it has a bounded slope, and it satisfies the no-sabotage condition when arson-type actions are free. Any contract that mixes a deductible, coinsurance and an upper limit is manipulation-proof. We also show that the ability to perform arson-type actions reduces the insured's welfare as less coverage is offered in equilibrium.


title: Composable constraints

abstract: We introduce a notion of compatibility between constraint encoding and compositional structure. Phrased in the language of category theory, it is given by a "composable constraint encoding". We show that every composable constraint encoding can be used to construct an equivalent notion of a constrained category in which morphisms are supplemented with the constraints they satisfy. We further describe how to express the compatibility of constraints with additional categorical structures of their targets, such as parallel composition, compactness, and time-symmetry. We present a variety of concrete examples. Some are familiar in the study of quantum protocols and quantum foundations, such as signalling and sectorial constraints; others arise by construction from basic categorical notions. We use the language developed to discuss the notion of intersectability of constraints and the simplifications it allows for when present, and to show that any time-symmetric theory of relational constraints admits a faithful notion of intersection.


title: Contracting on average iterated function systems by metric change

abstract: We study contraction conditions for an iterated function system of continuous maps on a metric space which are chosen randomly, identically and independently. We investigate metric changes, preserving the topological structure of the space, which turn the IFS into one which is contracting on average. For the particular case of a system of $C^1$-diffeomorphisms of the circle which is proximal and does not have a probability measure simultaneously invariant by every map, we derive a strongly equivalent metric which contracts on average.


title: Dynamical photon-photon interaction mediated by a quantum emitter

abstract: Single photons constitute a main platform in quantum science and technology: they carry quantum information over extended distances in the future quantum internet and can be manipulated in advanced photonic circuits enabling scalable photonic quantum computing. The main challenge in quantum photonics is how to generate advanced entangled resource states and efficient light-matter interfaces. Here we utilize the efficient and coherent coupling of a single quantum emitter to a nanophotonic waveguide for realizing quantum nonlinear interaction between single-photon wavepackets. This inherently multimode quantum system constitutes a new research frontier in quantum optics. We demonstrate control of a photon with another photon and experimentally unravel the dynamical response of two-photon interactions mediated by a quantum emitter, and show that the induced quantum correlations are controlled by the pulse duration. The work will open new avenues for tailoring complex photonic quantum resource states.


title: Power-law holographic dark energy and cosmology

abstract: We formulate power-law holographic dark energy, which is a modified holographic dark energy model based on the extended entropy relation arising from the consideration of state mixing between the ground and the excited ones in the calculation of the entanglement entropy. We construct two cases of the scenario, imposing the usual future event horizon choice, as well as the Hubble one. Thus, the former model is a one-parameter extension of standard holographic dark energy, recovering it in the limit where power-law extended entropy recovers Bekenstein-Hawking one, while the latter belongs to the class of running vacuum models, a feature that may reveal the connection between holography and the renormalization group running. For both models we extract the differential equation that determines the evolution of the dark-energy density parameter and we provide the expression for the corresponding equation-of-state parameter. We find that the scenario can describe the sequence of epochs in the Universe evolution, namely the domination of matter followed by the domination of dark energy. Moreover, the dark-energy equation of state presents a rich behavior, lying in the quintessence regime or passing into the phantom one too, depending on the values of the two model parameters, a behavior that is richer than the one of standard holographic dark energy.


title: Quantile Regression under Limited Dependent Variable

abstract: A new Stata command, ldvqreg, is developed to estimate quantile regression models for the cases of censored (with lower and/or upper censoring) and binary dependent variables. The estimators are implemented using a smoothed version of the quantile regression objective function. Simulation exercises show that it correctly estimates the parameters and it should be implemented instead of the available quantile regression methods when censoring is present. An empirical application to women's labor supply in Uruguay is considered.


title: Multi-Asset Spot and Option Market Simulation

abstract: We construct realistic spot and equity option market simulators for a single underlying on the basis of normalizing flows. We address the high-dimensionality of market observed call prices through an arbitrage-free autoencoder that approximates efficient low-dimensional representations of the prices while maintaining no static arbitrage in the reconstructed surface. Given a multi-asset universe, we leverage the conditional invertibility property of normalizing flows and introduce a scalable method to calibrate the joint distribution of a set of independent simulators while preserving the dynamics of each simulator. Empirical results highlight the goodness of the calibrated simulators and their fidelity.

categories: q-fin.CP cs.LG q-fin.MF q-fin.ST stat.ML


title: A short review of the Casimir effect with emphasis on dynamical boundary conditions

abstract: We give a short review on the static and dynamical Casimir effects, recalling their historical prediction, as well as their more recent experimental verification. We emphasise on the central role played by so-called {\it dynamical boundary conditions} (for which the boundary condition depends on a second time derivative of the field) in the experimental verification of the dynamical Casimir effect by Wilson et al. We then go on to review our previous work on the static Casimir effect with dynamical boundary conditions, providing an overview on how to compute the so-called local Casimir energy, the total Casimir energy and the Casimir force. We give as a future perspective the direction in which this work should be generalised to put the theoretical predictions of the dynamical Casimir effect experiments on a rigorous footing.


title: VL-Adapter: Parameter-Efficient Transfer Learning for Vision-and-Language Tasks

abstract: Recently, fine-tuning language models pre-trained on large text corpora have provided huge improvements on vision-and-language (V&L) tasks as well as on pure language tasks. However, fine-tuning the entire parameter set of pre-trained models becomes impractical since the model size is growing rapidly. Hence, in this paper, we introduce adapter-based parameter-efficient transfer learning techniques to V&L models such as VL-BART and VLT5. We evaluate our methods in a unified multi-task setup on both image-text and video-text benchmarks. For the image-text tasks, we use four diverse V&L datasets: VQAv2, GQA, NLVR2 , and MSCOCO image captioning. For video-text tasks, we use TVQA, How2QA, TVC, and YC2C. With careful training and thorough experiments, we benchmark three popular adapter-based methods (Adapter, Hyperformer, Compacter) against the standard full fine-tuning and the recently proposed prompt-tuning approach. We also enhance the efficiency and performance of adapters by sharing their weights to attain knowledge across tasks. Our results demonstrate that training the adapter with the weight-sharing technique (4.18% of total parameters for image-text tasks and 3.39% for video-text tasks) can match the performance of fine-tuning the entire model. Lastly, we present a comprehensive analysis including the combination of adapter and task-specific prompts and the impact of V&L pre-training on adapters. Our code is available at: https://github.com/ylsung/VL_adapter.


title: BScNets: Block Simplicial Complex Neural Networks

abstract: Simplicial neural networks (SNN) have recently emerged as the newest direction in graph learning which expands the idea of convolutional architectures from node space to simplicial complexes on graphs. Instead of pre-dominantly assessing pairwise relations among nodes as in the current practice, simplicial complexes allow us to describe higher-order interactions and multi-node graph structures. By building upon connection between the convolution operation and the new block Hodge-Laplacian, we propose the first SNN for link prediction. Our new Block Simplicial Complex Neural Networks (BScNets) model generalizes the existing graph convolutional network (GCN) frameworks by systematically incorporating salient interactions among multiple higher-order graph structures of different dimensions. We discuss theoretical foundations behind BScNets and illustrate its utility for link prediction on eight real-world and synthetic datasets. Our experiments indicate that BScNets outperforms the state-of-the-art models by a significant margin while maintaining low computation costs. Finally, we show utility of BScNets as the new promising alternative for tracking spread of infectious diseases such as COVID-19 and measuring the effectiveness of the healthcare risk mitigation strategies.


title: Making Fanaroff-Riley I radio sources III. The effects of the magnetic field on relativistic jets' propagation and source morphologies

abstract: Extragalactic radio sources appear under different morphologies, the most frequent ones are classified as Fanaroff-Riley type I (FR I), typically with lower luminosities, and Fanaroff-Riley type II, (FR II), typically more luminous. This simple classification, however, has many exceptions that we intend to investigate. Following previous analyses in the three-dimensional Hydrodynamic and Magneto-Hydrodynamic limits, we extend the numerical investigation to the Relativistic Magneto-Hydrodynamic regime, to include sources whose jet kinetic power sets in the range that separates FR Is from FR IIs. We consider weakly and mildly relativistic, underdense, supersonic jets that propagate in a stratified medium. In the model, the ambient temperature increases with distance from the jet origin maintaining constant pressure. We present three cases with low, high and intermediate kinetic luminosity that evolve into different morphologies. We find that the resulting morphology can be highly time dependent and that, apart from the jet power, the jet-to-ambient density ratio and the magnetization parameter play a crucial role in the jet evolution as well.


title: Statures and Sobrification Ranks of Noetherian Spaces

abstract: There is a rich theory of maximal order types of well-partial-orders (wpos), pioneered by de Jongh and Parikh (1977) and Schmidt (1981). Every wpo is Noetherian in its Alexandroff topology, and there are more; this prompts us to investigate an analogue of that theory in the wider context of Noetherian spaces. The notion of maximal order type does not seem to have a direct analogue in Noetherian spaces per se, but the equivalent notion of stature, investigated by Blass and Gurevich (2008) does: we define the stature $||X||$ of a Noetherian space $X$ as the ordinal rank of its poset of proper closed subsets. We obtain formulas for statures of sums, of products, of the space of words on a space $X$, of the space of finite multisets on $X$, in particular. They confirm previously known formulas on wpos, and extend them to Noetherian spaces. The proofs are, by necessity, rather different from their wpo counterparts, and rely on explicit characterizations of the sobrifications of the corresponding spaces, as obtained by Finkel and the first author (2020). We also give formulas for the statures of some natural Noetherian spaces that do not arise from wpos: spaces with the cofinite topology, Hoare powerspaces, powersets, and spaces of words on $X$ with the so-called prefix topology. Finally, because our proofs require it, and also because of its independent interest, we give formulas for the ordinal ranks of the sobrifications of each of those spaces, which we call their sobrification ranks.


title: Stabilized finite element methods for a fully-implicit logarithmic reformulation of the Oldroyd-B constitutive law

abstract: Logarithmic conformation reformulations for viscoelastic constitutive laws have alleviated the high Weissenberg number problem, and the exploration of highly elastic flows became possible. However, stabilized formulations for logarithmic conformation reformulations in the context of finite element methods have not yet been widely studied. We present stabilized formulations for the logarithmic reformulation of the Oldroyd-B model by Saramito (2014) based on the Variational Multiscale framework and Galerkin/Least-Squares method. The reformulation allows the use of Newton's method due to its fully-implicit nature for solving the steady-state problem directly. The proposed stabilization methods cure instabilities of convection and compatibility while preserving a three-field problem. The spatial accuracy of the formulations is assessed with the four-roll periodic box, revealing comparable accuracy between the methods. The formulations are validated with benchmark flows past a cylinder and through a 4:1 contraction. We found an excellent agreement to benchmark results in the literature. The algebraic sub-grid scale method is highly robust, indicated by the high limiting Weissenberg numbers in the benchmark flows.


title: DSPS: Differentiable Stellar Population Synthesis

abstract: Models of stellar population synthesis (SPS) are the fundamental tool that relates the physical properties of a galaxy to its spectral energy distribution (SED). In this paper, we present DSPS: a python package for stellar population synthesis. All of the functionality in DSPS is implemented natively in the JAX library for automatic differentiation, and so our predictions for galaxy photometry are fully differentiable, and directly inherit the performance benefits of JAX, including portability onto GPUs. DSPS also implements several novel features, such as i) a flexible empirical model for stellar metallicity that incorporates correlations with stellar age, and ii) support for the diffstar model that provides a physically-motivated connection between the star formation history of a galaxy (SFH) and the mass assembly of its underlying dark matter halo. We detail a set of theoretical techniques for using autodiff to calculate gradients of predictions for galaxy SEDs with respect to SPS parameters that control a range of physical effects, including SFH, stellar metallicity, nebular emission, and dust attenuation. When forward modeling the colors of a synthetic galaxy population, we find that DSPS can provide a factor of 5 speedup over standard SPS codes on a CPU, and a factor of 300-400 on a modern GPU. When coupled with gradient-based techniques for optimization and inference, DSPS makes it practical to conduct expansive likelihood analyses of simulation-based models of the galaxy--halo connection that fully forward model galaxy spectra and photometry.


title: Incorporating Measurement Error in Astronomical Object Classification

abstract: Most general-purpose classification methods, such as support-vector machine (SVM) and random forest (RF), fail to account for an unusual characteristic of astronomical data: known measurement error uncertainties. In astronomical data, this information is often given in the data but discarded because popular machine learning classifiers cannot incorporate it. We propose a simulation-based approach that incorporates heteroscedastic measurement error into existing classification method to better quantify uncertainty in classification. The proposed method first simulates perturbed realizations of the data from a Bayesian posterior predictive distribution of a Gaussian measurement error model. Then, a chosen classifier is fit to each simulation. The variation across the simulations naturally reflects the uncertainty propagated from the measurement errors in both labeled and unlabeled data sets. We demonstrate the use of this approach via two numerical studies. The first is a thorough simulation study applying the proposed procedure to SVM and RF, which are well-known hard and soft classifiers, respectively. The second study is a realistic classification problem of identifying high-$z$ $(2.9 \leq z \leq 5.1)$ quasar candidates from photometric data. The data are from merged catalogs of the Sloan Digital Sky Survey, the $Spitzer$ IRAC Equatorial Survey, and the $Spitzer$-HETDEX Exploratory Large-Area Survey. The proposed approach reveals that out of 11,847 high-$z$ quasar candidates identified by a random forest without incorporating measurement error, 3,146 are potential misclassifications with measurement error. Additionally, out of $1.85$ million objects not identified as high-$z$ quasars without measurement error, 936 can be considered new candidates with measurement error.


title: Seeking the Building Blocks of Visual Imagery and Creativity in a Cognitively Inspired Neural Network

abstract: How do we imagine visual objects and combine them to create new forms? To answer this question, we need to explore the cognitive, computational and neural mechanisms underlying imagery and creativity. The body of research on deep learning models with creative behaviors is growing. However, in this paper we suggest that the complexity of such models and their training sets is an impediment to using them as tools to understand human aspects of creativity. We propose using simpler models, inspired by neural and cognitive mechanisms, that are trained with smaller data sets. We show that a standard deep learning architecture can demonstrate imagery by generating shape/color combinations using only symbolic codes as input. However, generating a new combination that was not experienced by the model was not possible. We discuss the limitations of such models, and explain how creativity could be embedded by incorporating memory mechanisms to combine the output of the network into new combinations and use that as new training data.


title: Beyond Ads: Sequential Decision-Making Algorithms in Law and Public Policy

abstract: We explore the promises and challenges of employing sequential decision-making algorithms -- such as bandits, reinforcement learning, and active learning -- in law and public policy. While such algorithms have well-characterized performance in the private sector (e.g., online advertising), the tendency to naively apply algorithms motivated by one domain, often online advertisements, can be called the "advertisement fallacy." Our main thesis is that law and public policy pose distinct methodological challenges that the machine learning community has not yet addressed. Machine learning will need to address these methodological problems to move "beyond ads." Public law, for instance, can pose multiple objectives, necessitate batched and delayed feedback, and require systems to learn rational, causal decision-making policies, each of which presents novel questions at the research frontier. We discuss a wide range of potential applications of sequential decision-making algorithms in regulation and governance, including public health, environmental protection, tax administration, occupational safety, and benefits adjudication. We use these examples to highlight research needed to render sequential decision making policy-compliant, adaptable, and effective in the public sector. We also note the potential risks of such deployments and describe how sequential decision systems can also facilitate the discovery of harms. We hope our work inspires more investigation of sequential decision making in law and public policy, which provide unique challenges for machine learning researchers with potential for significant social benefit.


title: Mapping and fixed point property theorems for inverse limits with set-valued bonding functions

abstract: Among other results, the paper gives new mapping theorems and new fixed point property theorems for inverse limits of inverse sequences of compact metric spaces with upper semicontinuous set-valued bonding functions. We also revisit the results from two papers, Mioduszewski's ''Mappings of inverse limits'' and Feuerbacher's ''Mappings of inverse limits revisited''.


title: Rotons and their damping in elongated dipolar Bose-Einstein condensates

abstract: We discuss finite temperature damping of rotons in elongated Bose-condensed dipolar gases, which are in the Thomas-Fermi regime in the tightly confined directions. The presence of many branches of excitations which can participate in the damping process, is crucial for the Landau damping and results in significant increase of the damping rate. It is found, however, that even rotons with energies close to the roton gap may remain fairly stable in systems with the roton gap as small as 1nK.


title: Loop induced $h \rightarrow \gamma \gamma, Z\gamma, gg $ decays in 3-Higgs doublet models

abstract: In this paper, we analyze the loop induced decays of neutral Higgs boson $ h $ into pairs of photons, gluons, and Z$ \gamma $ in 3-Higgs Doublet Models (3HDMs), where all the doublets have non-zero vacuum expectation values. In the 3HDM, there are contributions to these decays from loops of charged scalars $ H^{\pm}_{1,2} $ that are not present in the Standard Model. We find an enhancement in the channel $ h\rightarrow Z\gamma $ with respect to the Standard Model prediction in some parameter region while respecting constraints from the measurements of $ h\rightarrow \gamma\gamma $ at the LHC. The charged Higgs bosons $ H^{\pm}_{1,2} $ lead to the enhancement of $ Z\gamma $ decay mode in some part of the parameter space of the 3HDM. Then, we discuss the correlations between the decay modes $ h\rightarrow Z\gamma $ and $ h\rightarrow \gamma\gamma $. Moreover, we study Type-I and Type-Z of 3HDMs and find that these types of 3HDMs can be distinguished via only the gluonic decays, since the contributions of the $ H^{\pm}_{1,2} $ loops are dominant to those of fermion loops. We also investigate the effects of the model parameters on these loop induced decays. As an interesting phenomenological consequence, the 3HDM is dominantly governed by the parameters: two ratios of vacuum expectation values (VEVs); $ \tan\gamma,~\tan\beta $ and the mixing angles of the CP-even Higgs bosons, $ \alpha_{1,2}~ \text{and}~\theta $.


title: Sparse Interventions in Language Models with Differentiable Masking

abstract: There has been a lot of interest in understanding what information is captured by hidden representations of language models (LMs). Typically, interpretation methods i) do not guarantee that the model actually uses the encoded information, and ii) do not discover small subsets of neurons responsible for a considered phenomenon. Inspired by causal mediation analysis, we propose a method that discovers within a neural LM a small subset of neurons responsible for a particular linguistic phenomenon, i.e., subsets causing a change in the corresponding token emission probabilities. We use a differentiable relaxation to approximately search through the combinatorial space. An $L_0$ regularization term ensures that the search converges to discrete and sparse solutions. We apply our method to analyze subject-verb number agreement and gender bias detection in LSTMs. We observe that it is fast and finds better solutions than the alternative (REINFORCE). Our experiments confirm that each of these phenomenons is mediated through a small subset of neurons that do not play any other discernible role.


title: Transformation of an energy spectrum and wave functions in the crossover from two- to three-dimensional topological insulator in HgTe quantum wells: long and thorny way

abstract: A magnetotransport and quantum capacitance of the two-dimensional electron gas in HgTe/Cd$_x$Hg$_{1-x}$Te quantum wells of a width ($20.2-46.0$)~nm are experimentally investigated. It is shown that the first energy subband of spatial quantization is split due to the spin-orbit interaction and the split branches are single-spin, therewith the splitting strength increases with the increase of the quantum well width. The electron effective masses in the branches are close to each other within the actual density range. Magneto-intersubband oscillations (MISO) observed in the structures under study exhibit the growing amplitude with the increasing electron density that contradicts to the expected decrease of wave function overlap for the rectangular quantum well. To interpret the data obtained, we have used a self-consistent approach to calculate the electron energy spectrum and the wave function within framework of the \emph{kP}-model. It has been in particular shown that the MISO amplitude increase results from the increasing overlap of the wave functions due to their shift from the gate electrode with the gate voltage increase known as phenomenon of the negative electron polarizability. The results obtained from the transport experiments are supported by quantum capacitance measurements.


title: Diffractive longitudinal structure function at the Electron Ion Collider

abstract: Possibilities for the measurement of the longitudinal structure function in diffraction $F_\mathrm{L}^\mathrm{D}$ at the future US Electron Ion Collider are investigated. The sensitivity to $F_\mathrm{L}^\mathrm{D}$ arises from the variation of the reduced diffractive cross section with centre-of-mass energy. Simulations are performed with various sets of beam energy combinations and for different assumptions on the precision of the diffractive cross section measurements. Scenarios compatible with current EIC performance expectations lead to an unprecedented precision on $F_\mathrm{L}^\mathrm{D}$ at the 5-10 % level in the best measured regions. While scenarios with data at a larger number of centre-of-mass energies allow the extraction of $F_\mathrm{L}^\mathrm{D}$ in the widest kinematic domain and with the smallest uncertainties, even the more conservative assumptions lead to precise measurements. The ratio $R^\mathrm{D}$ of photoabsorption cross sections for longitudinally to transversely polarised photons can also be obtained with high precision using a separate extraction method.


title: The $\hbar\rightarrow 0$ Limit of the Entanglement Entropy

abstract: Entangled quantum states share properties that do not have classical analogs, in particular, they show correlations that can violate Bell inequalities. It is therefore an interesting question to see what happens to entanglement measures -- such as the entanglement entropy for a pure state -- taking the semi-classical limit, where the naive expectation is that they may become singular or zero. This conclusion is however incorrect. In this paper, we determine the $\hbar\rightarrow 0$ limit of the bipartite entanglement entropy for a one-dimensional system of $N$ quantum particles in an external potential and we explicitly show that this limit is finite. Moreover, if the particles are fermionic, we show that the $\hbar\rightarrow 0$ limit of the bipartite entanglement entropy coincides with the Shannon entropy of $N$ bits.


title: Probably approximately correct quantum source coding

abstract: Information-theoretic lower bounds are often encountered in several branches of computer science, including learning theory and cryptography. In the quantum setting, Holevo's and Nayak's bounds give an estimate of the amount of classical information that can be stored in a quantum state. Previous works have shown how to combine information-theoretic tools with a counting argument to lower bound the sample complexity of distribution-free quantum probably approximately correct (PAC) learning. In our work, we establish the notion of Probably Approximately Correct Source Coding and we show two novel applications in quantum learning theory and delegated quantum computation with a purely classical client. In particular, we provide a lower bound of the sample complexity of a quantum learner for arbitrary functions under the Zipf distribution, and we improve the security guarantees of a classically-driven delegation protocol for measurement-based quantum computation (MBQC).


title: Trapping of Planar Brownian Motion: Full First Passage Time Distributions by Kinetic Monte-Carlo, Asymptotic and Boundary Integral Methods

abstract: We consider the problem of determining the arrival statistics of unbiased planar random walkers to complex target configurations. In contrast to problems posed in finite domains, simple moments of the distribution, such as the mean (MFPT) and variance, are not defined and it is necessary to obtain the full arrival statistics. We describe several methods to obtain these distributions and other associated quantities such as splitting probabilities. One approach combines a Laplace transform of the underlying parabolic equation with matched asymptotic analysis followed by numerical transform inversion. The second approach is similar, but uses a boundary integral equation method to solve for the Laplace transformed variable. To validate the results of this theory, and to obtain the arrival time statistics in very general configurations of absorbers, we introduce an efficient Kinetic Monte Carlo (KMC) method that describes trajectories as a combination of large but exactly solvable projection steps. The effectiveness of these methodologies is demonstrated on a variety of challenging examples highlighting the applicability of these methods to a variety of practical scenarios, such as source inference. A particularly useful finding arising from these results is that homogenization theories, in which complex configurations are replaced by equivalent simple ones, are remarkably effective at describing arrival time statistics.


title: Toric Promotion

abstract: We introduce toric promotion as a cyclic analogue of Sch\"utzenberger's promotion operator. Toric promotion acts on the set of labelings of a graph $G$. We discuss connections between toric promotion and previously-studied notions such as toric posets and friends-and-strangers graphs. Our main theorem provides a surprisingly simple description of the orbit structure of toric promotion when $G$ is a forest.


title: Optimal friction matrix for underdamped Langevin sampling

abstract: A systematic procedure for optimising the friction coefficient in underdamped Langevin dynamics as a sampling tool is given by taking the gradient of the associated asymptotic variance with respect to friction. We give an expression for this gradient in terms of the solution to an appropriate Poisson equation and show that it can be approximated by short simulations of the associated first variation/tangent process under concavity assumptions on the log density. Our algorithm is applied to the estimation of posterior means in Bayesian inference problems and reduced variance is demonstrated when compared to the original underdamped and overdamped Langevin dynamics in both full and stochastic gradient cases.


title: Artificial selection of communities drives the emergence of structured interactions

abstract: Species-rich communities, such as the microbiota or microbial ecosystems, provide key functions for human health and climatic resilience. Increasing effort is being dedicated to design experimental protocols for selecting community-level functions of interest. These experiments typically involve selection acting on populations of communities, each of which is composed of multiple species. If numerical simulations started to explore the evolutionary dynamics of this complex, multi-scale system, a comprehensive theoretical understanding of the process of artificial selection of communities is still lacking. Here, we propose a general model for the evolutionary dynamics of communities composed of a large number of interacting species, described by disordered generalised Lotka-Volterra equations. Our analytical and numerical results reveal that selection for scalar community functions leads to the emergence, along an evolutionary trajectory, of a low-dimensional structure in an initially featureless interaction matrix. Such structure reflects the combination of the properties of the ancestral community and of the selective pressure. Our analysis determines how the speed of adaptation scales with the system parameters and the abundance distribution of the evolved communities. Artificial selection for larger total abundance is thus shown to drive increased levels of mutualism and interaction diversity. Inference of the interaction matrix is proposed as a method to assess the emergence of structured interactions from experimentally accessible measures.


title: Extremal points of total generalized variation balls in 1D: characterization and applications

abstract: The total generalized variation (TGV) is a popular regularizer in inverse problems and imaging combining discontinuous solutions and higher order smoothing. In particular, empirical observations suggest that its order two version strongly favors piecewise affine functions. In the present manuscript, we formalize this statement for the one-dimensional TGV-functional by characterizing the extremal points of its sublevel sets with respect to a suitable quotient space topology. These results imply that 1D TGV-regularized linear inverse problems with finite dimensional observations admit piecewise affine minimizers. As further applications of this characterization we include precise first-order necessary optimality conditions without requiring convexity of the fidelity term, and a simple solution algorithm for TGV-regularized minimization problems.


title: Singlet-Doublet Self-interacting Dark Matter and Radiative Neutrino Mass

abstract: Self-interacting dark matter (SIDM) with a light mediator is a promising scenario to alleviate the small-scale problems of the cold dark matter paradigm while being consistent with the latter at large scales, as suggested by astrophysical observations. This, however, leads to an under-abundant SIDM relic due to large annihilation rates into mediator particles, often requiring an extension of the simplest thermal or non-thermal relic generation mechanism. In this work, we consider a singlet-doublet fermion dark matter scenario where the singlet fermion with a light scalar mediator gives rise to the velocity-dependent dark matter self-interaction through a Yukawa type attractive potential. The doublet fermion, by virtue of its tiny mixing with the singlet, can be long-lived and can provide a non-thermal contribution to the singlet relic at late epochs, filling the deficit in the thermal relic of singlet SIDM. The light scalar mediator, due to its mixing with the standard model Higgs, paves the path for detecting such SIDM at terrestrial laboratories leading to constraints on model parameters from CRESST-III and XENON1T experiments. Enlarging the dark sector particles by two more singlet fermions and one scalar doublet, all odd under an unbroken $\mathcal{Z}_2$ symmetry can also explain non-zero neutrino mass in scotogenic fashion.


title: Peer-to-Peer Communication Trade-Offs for Smart Grid Applications

abstract: Virtual topologies in peer-to-peer networks can reduce the traffic consumed by altering the logical connectivity of peers without altering the underlying network. However, such sparsely connected virtual topologies do not focus on the needs for smart grid applications, which is information dissemination throughout the network, and in turn degrade the performance of distributed control algorithms running on peer-to-peer networks. This paper provides a flexible solution for application developers to prototype and deploy different virtual topologies that balances these trade-offs. First, it introduces a configurable virtual communication topology framework, TopLinkMgr, which enables users to specify any chosen connectivity configuration and deploy peer-to-peer applications using it. Second, it proposes a novel fault-tolerant self-adaptive virtual topology management algorithm, Bounded Path Dissemination, that can ensure the dissemination of information to all peers within a specified threshold. Experiments show that the algorithm improves on convergence speed and accuracy over state-of-the-art methods and is also robust against node failures while consuming significantly less communication bandwidth.


title: Infinite geodesics, competition interfaces and the second class particle in the scaling limit

abstract: We establish fundamental properties of infinite geodesics and competition interfaces in the directed landscape. We construct infinite geodesics in the directed landscape, establish their uniqueness and coalescence, and define Busemann functions. We then define competition interfaces in the directed landscape. We prove the second class particle in tasep converges under KPZ scaling to a competition interface. Under suitable conditions, we show the competition interface has an asymptotic direction, analogous to the speed of a second class particle, and determine its law. Moreover, we prove the competition interface has an absolutely continuous law on compact sets with respect to infinite geodesics.


title: Chiralspin symmetry and confinement

abstract: Interesting lattice QCD simulations at high temperature in QCD and particular truncated studies have shown the emergence of an unexpected group symmetry, so called \textit{chiralspin}. However this is not a symmetry of the QCD action for free quarks, which makes unclear the transition to deconfinement at high temperature in QCD. Therefore we try to redefine this group so that is a symmetry of free quark action and it is consistent with the presence of deconfinement in QCD.


title: The Taylor-Wiles method for coherent cohomology, II

abstract: We show that the Taylor-Wiles method can be applied to the cohomology of a Shimura variety $S$ of PEL type attached to a unitary similitude group $G$, with coefficients in the coherent sheaf attached to an automorphic vector bundle $\CF$ , when $S$ has a smooth model over a $p$-adic integer ring. This generalizes the main results of the article \cite{H13}, which treated the case when $S$ is compact. As in the previous article, the starting point is a theorem of Lan and Suh that proves the vanishing of torsion in the cohomology under certain conditions on the parameters of the bundle $\CF$ and the prime $p$. Most of the additional difficulty in the non-compact case is related to showing that the contributions of boundary cohomology are all of Eisenstein type. We also need to show that the coverings giving rise to the diamond operators can be extended to \'etale coverings of appropriate toroidal compactifications. The result is applied to show that, when the Taylor-Wiles method applies, the congruence ideal attached to a coherent cohomological realization of an automorphic Galois representation is independent of the signatures of the hermitian form to which $G$ is attached. We also show that the Gorenstein hypothesis used to construct $p$-adic $L$-functions in \cite{EHLS} is valid under rather general hypotheses.


title: Symmetry constraints on spin order transfer in parahydro-gen-induced polarization (PHIP)

abstract: It is well known that the association of parahydrogen (pH2) with an unsaturated molecule or a transient metalorganic complex can enhance the intensity of NMR signals; the effect is known as parahydrogen-induced polarization (PHIP). During the last decades, numerous methods were proposed for converting pH2-derived nuclear spin order to the observable magnetization of pro-tons or other nuclei of interest, usually 13C or 15N. Here, we analyze the constraints imposed by the topological symmetry of the spin systems on the amplitude of transferred polarization. In asymmetric systems, heteronuclei can be polarized to 100%. However, the amplitude drops to 75% in A2BX systems and further to 50% in A3B2X systems. The latter case is of primary importance for biological applications of PHIP using sidearm hydrogenation (PHIP-SAH). If the polarization is transferred to the same type of nuclei, i.e. 1H, symmetry constraints impose significant boundaries on the spin-order distribution. For AB, A2B, A3B, A2B2, AA'(AA') systems, the maxi-mum average polarization for each spin is 100%, 50%, 33.3%, 25%, and 0, respectively, when A and B (or A') came from pH2 We also discuss the effect of dipole-dipole induced pH2 spin-order distribution in heterogeneous catalysis or nematic liquid crystals. Practical examples from the literature illustrate our theoretical analysis.


title: The whole and the parts: the MDL principle and the a-contrario framework

abstract: This work explores the connections between the Minimum Description Length (MDL) principle as developed by Rissanen, and the a-contrario framework for structure detection proposed by Desolneux, Moisan and Morel. The MDL principle focuses on the best interpretation for the whole data while the a-contrario approach concentrates on detecting parts of the data with anomalous statistics. Although framed in different theoretical formalisms, we show that both methodologies share many common concepts and tools in their machinery and yield very similar formulations in a number of interesting scenarios ranging from simple toy examples to practical applications such as polygonal approximation of curves and line segment detection in images. We also formulate the conditions under which both approaches are formally equivalent.


title: Scheduled Relaxation Jacobi schemes for non-elliptic partial differential equations

abstract: The Scheduled Relaxation Jacobi (SRJ) method is a linear solver algorithm which greatly improves the convergence of the Jacobi iteration through the use of judiciously chosen relaxation factors (an SRJ scheme) which attenuate the solution error. Until now, the method has primarily been used to accelerate the solution of elliptic PDEs (e.g. Laplace, Poisson's equation) as the currently available schemes are restricted to solving this class of problems. The goal of this paper is to present a methodology for constructing SRJ schemes which are suitable for solving non-elliptic PDEs (or equivalent, nonsymmetric linear systems arising from the discretization of these PDEs), thereby extending the applicability of this method to a broader class of problems. These schemes are obtained by numerically solving a constrained minimization problem which guarantees the solution error will not grow as long as the linear system has eigenvalues which lie in certain regions of the complex plane. We demonstrate that these schemes are able to accelerate the convergence of standard Jacobi iteration for the nonsymmetric linear systems arising from discretization of the 1D and 2D steady advection-diffusion equations.


title: A mean-field theory approach to 3D nematic phase transitions in microtubules

abstract: Microtubules are dynamic intracellular fibers that have been observed experimentally to undergo spontaneous self-alignment. We formulate a 3D mean-field theory model to analyze the nematic phase transition of microtubules growing and interacting within a 3D space then make a comparison with computational simulations. We identify a control parameter $G_\text{eff}$ and predict a unique critical value $G_\text{eff}=1.56$ for which a phase transition can occur. Furthermore, we show both analytically and using simulations that this predicted critical value does not depend on the presence of zippering. The mean-field theory developed here provides an analytical estimate of microtubule patterning characteristics without running time-consuming simulations and is a step towards bridging scales from microtubule behavior to multicellular simulations.


title: Thermodynamics to infer the astrophysics of binary black hole mergers

abstract: We introduce the Merger Entropy Index ($\mathcal{I}_\mathrm{BBH}$), a new parameter to measure the efficiency of entropy transfer for any generic binary black hole merger in General Relativity. We find that $\mathcal{I}_\mathrm{BBH}$ is bounded between an asymptotic maximum and minimum. For the observed population of mergers detected by LIGO and Virgo, we find that $\mathcal{I}_\mathrm{BBH}$ is $\lesssim30\%$ of its theoretical maximum. By imposing the thermodynamical consistency between the pre- and post-merger states through $\mathcal{I}_\mathrm{BBH}$, we showcase BRAHMA -- a novel framework to infer the properties and astrophysical implications of gravitational-wave detections. For GW190521 -- the heaviest confirmed binary black hole merger observed so far -- our framework rules out high mass-ratio, negative effective inspiral spin, and electromagnetic counterpart claims. Furthermore, our analysis provides an independent confirmation that GW190521 belongs to a separate population.


title: Software solutions for numerical modeling of wide-field telescopes

abstract: This paper presents an integrated modeling software to analyze the PSF of wide-field telescopes affected by misalignments. Even relatively small misalignments in the optical system of a telescope can significantly deteriorate the image quality by introducing large aberrations. In particular, wide-field telescopes are critically affected by these errors, insomuch that usually a closed-loop active optics system is adopted for a continuous correction, rather than for sporadic alignment procedures. Typically, a ray-tracing software such as Zemax OpticStudio is employed to accurately analyze the system during the optical design. However, an analytical model of the optical system is preferable when the PSF of the telescope must be reconstructed quickly for algorithmic purposes. Here the analytical model is derived through a hybrid approach and developed in a custom software package, designed to be general and flexible in order to be tailored to different optical configurations. First, leveraging on the Zemax OpticStudio API, the ray-tracing software is integrated into a Matlab pipeline. This allows to perform a statistical analysis by automatically simulating the system response in a variety of misaligned working conditions. Then, the resulting dataset is employed to populate a database of parameters describing the model.