title: The buildup of galaxies and their spheroids: the contributions of mergers, disc instabilities and star formation

abstract: We use the GALFORM semi-analytical model of galaxy formation and the Planck-Millennium simulation to investigate the origins of stellar mass in galaxies and their spheroids. We compare the importance of mergers and disc instabilities, as well as the starbursts that they trigger. We find that the fraction of galaxy stellar mass formed \textit{ex situ} ($f_\mathrm{ex}$) increases sharply from $M_*=10^{11}$ M$_\odot$ upwards, reaching $80\%$ at $M_*=10^{11.3}$ M$_\odot$. For low-mass galaxies we find larger \textit{\textit{ex situ}} contributions at $z=0$ than in other models ($7$-$12\%$), with a decrease towards higher redshifts. The global \textit{ex situ} fraction of all stellar mass falls sharply with redshift, from $40\%$ at $z=0$ to $3\%$ at $z=10$. Major mergers contribute roughly half of the \textit{ex situ} mass, with minor mergers and smooth accretion of satellites both accounting for $\approx25\%$, almost independent of stellar mass and redshift. Mergers dominate in building up high-mass ($M_\mathrm{*,sph}>10^{11}$ M$_\odot$) and low-mass ($M_\mathrm{*,sph}<10^{8.5}$ M$_\odot$) spheroids. Disc instabilities and their associated starbursts dominate for intermediate-mass spheroids ($10^{8.5}<M_\mathrm{*,sph}<10^{11}$ M$_\odot$) at $z=0$. The mass regime where pseudobulges dominate is in agreement with observed pseudobulge fractions, but the peak value in the pseudobulge fraction predicted by GALFORM is likely too high. The total contributions of disc instabilities and their starbursts are roughly equal at $z=0$, with the former dominating for lower-mass spheroids (peak at $M_\mathrm{*,sph}=10^{9.5}$ M$_\odot$) and the latter for higher-mass ones (peak at $M_\mathrm{*,sph}=10^{10.5}$ M$_\odot$).


title: Chandra Observations of Six Peter Pan Disks: Diversity of X-ray-driven Internal Photoevaporation Rates Doesn't Explain Their Rare Longevity

abstract: We present Chandra X-ray observations of 6 previously-identified Peter Pan objects, rare 40 Myr systems with evidence of primordial disk retention. We observe X-ray luminosities (0.8-3.0 keV) ranging from log Lx 27.7-29.1. We find that our Peter Pan sample exhibits X-ray properties similar to that of weak-lined T-Tauri stars and do not exhibit evidence of stellar accretion induced X-ray suppression. Our observed Peter Pan X-ray luminosities are consistent with that measured for field dM stars of similar spectral type and age, implying their long primordial disk lifetimes are likely not a consequence of unusually faint X-ray host stars. Our derived X-ray photoevaporative mass loss rates predict our systems have passed the point of rapid gas dispersal and call into question the impact of this internal mechanism for primordial disk dispersal around dM stars. Our qualitative assessment of the surrounding Peter Pan environments also does not predict unusually low levels of external photoevaporation relative to other respective moving group members. Overall, our results suggest Peter Pan disks may be a consequence of the low FUV flux incident on the disk in low-mass DM stars given their relatively lower levels of accretion over the course of their pre-main-sequence evolution.


title: Jet Broadening in Flowing Matter -- Resummation

abstract: In this work, we obtain the leading corrections to the jet momentum broadening distribution in a QCD medium arising from the transverse flow of the matter. We first derive the single-particle propagator of a highly energetic parton resumming its multiple interactions with the homogeneous flowing matter, explicitly keeping the leading subeikonal flow terms. Then, we use this propagator to obtain the jet broadening distribution and its leading moments. We show that this distribution becomes anisotropic in the presence of transverse flow, since its odd moments are generally non-zero and proportional to the transverse velocity of the medium. Finally, we evaluate several odd moments, which we compare to the corresponding results at first order in opacity, showing that accounting for multiple in-medium scatterings is essential to describe some observables in dense nuclear matter.


title: Massive neutrino self-interactions and inflation

abstract: Certain inflationary models like Natural inflation (NI) and Coleman-Weinberg inflation (CWI) are disfavoured by cosmological data in the standard $\Lambda\textrm{CDM}+r$ model (where $r$ is the scalar-to-tensor ratio), as these inflationary models predict the regions in the $n_s-r$ parameter space that are excluded by the cosmological data at more than 2$\sigma$ (here $n_s$ is the scalar spectral index). The same is true for single field inflationary models with an inflection point that can account for all or majority of dark matter in the form of PBHs (primordial black holes). Cosmological models incorporating strongly self-interacting neutrinos (with a heavy mediator) are, however, known to prefer lower $n_s$ values compared to the $\Lambda\rm CDM$ model. Considering such neutrino self-interactions can, thus, open up the parameter space to accommodate the above inflationary models. In this work, we implement the massive neutrino self-interactions with a heavy mediator in two different ways: flavour-universal (among all three neutrinos), and flavour-specific (involving only one neutrino species). We implement the new interaction in both scalar and tensor perturbation equations of neutrinos. Interestingly, we find that the current cosmological data can support the aforementioned inflationary models at 2$\sigma$ in the presence of such neutrino self-interactions.


title: Node Replication: Theory And Practice

abstract: We define and study a term calculus implementing higher-order node replication. It is used to specify two different (weak) evaluation strategies: call-by-name and fully lazy call-by-need, that are shown to be observationally equivalent by using type theoretical technical tools.


title: Aiolos- A multi-purpose 1-D hydrodynamics code for planetary atmospheres

abstract: We present a new 1-D multi-physics simulation code with use cases intended for, but not limited to, hydrodynamic escapeproblems of planetary atmospheres and planetary accretion models. Our formulation treats an arbitrary number of species asseparated hydrodynamic fields, couples them via friction laws, allows for a multi-band flux-limited radiation transport, and tracksionization fronts in high-energy irradiation bands. Besides coupling various known numerical solution techniques together, weimprove on the numerical stability of deep hydrostatic atmospheres by using a well-balanced scheme, hence preventing unphysicaldriving of atmospheric in- or outflow. We demonstrate the correct physical behaviour of the individual code modules and presenta few simple, new applications, such as a proof-of-concept simulations of combined core-powered mass-loss and UV-drivenatmospheric escape, along with a fully time-dependent core-collapse giant planet simulation. The multi-species nature of thecode opens up the area of exploring simulations that are agnostic towards the dominant atmospheric species and can lead toimplementations of advanced planetary evolution schemes.


title: A lab-based test of the gravitational redshift with a miniature clock network

abstract: Einstein's theory of general relativity predicts that a clock at a higher gravitational potential will tick faster than an otherwise identical clock at a lower potential, an effect known as the gravitational redshift. Here we perform a laboratory-based, blinded test of the gravitational redshift using differential clock comparisons within an evenly spaced array of 5 atomic ensembles spanning a height difference of 1 cm. We measure a fractional frequency gradient of $[-12.4\pm0.7_{\rm{(stat)}}\pm2.5_{\rm{(sys)}}]\times10^{-19}/$cm, consistent with the expected redshift gradient of $-10.9\times10^{-19}/$cm. Our results can also be viewed as relativistic gravitational potential difference measurements with sensitivity to mm scale changes in height on the surface of the Earth. These results highlight the potential of local-oscillator-independent differential clock comparisons for emerging applications of optical atomic clocks including geodesy, searches for new physics, gravitational wave detection, and explorations of the interplay between quantum mechanics and gravity.


title: Evidence for Extended Hydrogen-Poor CSM in the Three-Peaked Light Curve of Stripped Envelope Ib Supernova

abstract: We present multi-band ATLAS photometry for SN 2019tsf, a stripped-envelope Type Ib supernova (SESN). The SN shows a triple-peaked light curve and a late (re-)brightening, making it unique among stripped-envelope systems. The re-brightening observations represent the latest photometric measurements of a multi-peaked Type Ib SN to date. As late-time photometry and spectroscopy suggest no hydrogen, the potential circumstellar material (CSM) must be H-poor. Moreover, late (>150 days) spectra show no signs of narrow emission lines, further disfavouring CSM interaction. On the contrary, an extended CSM structure is seen through a follow-up radio campaign with Karl G. Jansky Very Large Array (VLA), indicating a source of bright optically thick radio emission at late times, which is highly unusual among H-poor SESNe. We attribute this phenomenology to an interaction of the supernova ejecta with spherically-asymmetric CSM, potentially disk-like, and we present several models that can potentially explain the origin of this rare Type Ib supernova. The warped disc model paints a novel picture, where the tertiary companion perturbs the progenitors CSM, that can explain the multi-peaked light curves of SNe, and here we apply it to SN 2019tsf. This SN 2019tsf is likely a member of a new sub-class of Type Ib SNe and among the recently discovered class of SNe that undergo mass transfer at the moment of explosion


title: A classification of injective FI$^m$-modules

abstract: In this paper we generalize a shift theorem, which plays a key role in studying representations of FI$^m$, the product category of the category of finite sets and injections, and classify finitely generated injective FI$^m$-modules over a field of characteristic 0.


title: Quantum entropy expansion using n-qubit permutation matrices in Galois field

abstract: Random numbers are critical for any cryptographic application. However, the data that is flowing through the internet is not secure because of entropy deprived pseudo random number generators and unencrypted IoTs. In this work, we address the issue of lesser entropy of several data formats. Specifically, we use the large information space associated with the n-qubit permutation matrices to expand the entropy of any data without increasing the size of the data. We take English text with the entropy in the range 4 - 5 bits per byte. We manipulate the data using a set of n-qubit (n $\leq$ 10) permutation matrices and observe the expansion of the entropy in the manipulated data (to more than 7.9 bits per byte). We also observe similar behaviour with other data formats like image, audio etc. (n $\leq$ 15).


title: Bug Fix Time Optimization Using Matrix Factorization and Iterative Gale-Shaply Algorithms

abstract: Bug triage is an essential task in software maintenance phase. It assigns developers (fixers) to bug reports to fix them. This process is performed manually by a triager, who analyzes developers profiles and submitted bug reports to make suitable assignments. Bug triaging process is time consuming thus automating this process is essential to improve the quality of software. Previous work addressed triaging problem either as an information retrieval or classification problem. This paper tackles this problem as a resource allocation problem, that aims at the best assignments of developers to bug reports, that reduces the total fixing time of the newly submitted bug reports, in addition to the even distribution of bug reports over developers. In this paper, a combination of matrix factorization and Gale Shapely algorithm, supported by the differential evolution is firstly introduced to optimize the total fix time and normalize developers work load. Matrix factorization is used to establish a recommendation system for Gale-Shapley to make assignment decisions. Differential evolution provides the best set of weights to build developers score profiles. The proposed approach is assessed over three repositories, Linux, Apache and Eclipse. Experimental results show that the proposed approach reduces the bug fixing time, in comparison to the manual triage, by 80.67%, 23.61% and 60.22% over Linux, Eclipse and Apache respectively. Moreover, the workload for the developers is uniform.


title: Making Linear MDPs Practical via Contrastive Representation Learning

abstract: It is common to address the curse of dimensionality in Markov decision processes (MDPs) by exploiting low-rank representations. This motivates much of the recent theoretical study on linear MDPs. However, most approaches require a given representation under unrealistic assumptions about the normalization of the decomposition or introduce unresolved computational challenges in practice. Instead, we consider an alternative definition of linear MDPs that automatically ensures normalization while allowing efficient representation learning via contrastive estimation. The framework also admits confidence-adjusted index algorithms, enabling an efficient and principled approach to incorporating optimism or pessimism in the face of uncertainty. To the best of our knowledge, this provides the first practical representation learning method for linear MDPs that achieves both strong theoretical guarantees and empirical performance. Theoretically, we prove that the proposed algorithm is sample efficient in both the online and offline settings. Empirically, we demonstrate superior performance over existing state-of-the-art model-based and model-free algorithms on several benchmarks.


title: Measuring the thermal and ionization state of the low-$z$ IGM using likelihood free inference

abstract: We present a new approach to measure the power-law temperature density relationship $T=T_0 (\rho / \bar{\rho})^{\gamma -1}$ and the UV background photoionization rate $\Gamma_{\rm HI}$ of the IGM based on the Voigt profile decomposition of the Ly$\alpha$ forest into a set of discrete absorption lines with Doppler parameter $b$ and the neutral hydrogen column density $N_{\rm HI}$. Previous work demonstrated that the shape of the $b$-$N_{\rm HI}$ distribution is sensitive to the IGM thermal parameters $T_0$ and $\gamma$, whereas our new inference algorithm also takes into account the normalization of the distribution, i.e. the line-density d$N$/d$z$, and we demonstrate that precise constraints can also be obtained on $\Gamma_{\rm HI}$. We use density-estimation likelihood-free inference (DELFI) to emulate the dependence of the $b$-$N_{\rm HI}$ distribution on IGM parameters trained on an ensemble of 624 Nyx hydrodynamical simulations at $z = 0.1$, which we combine with a Gaussian process emulator of the normalization. To demonstrate the efficacy of this approach, we generate hundreds of realizations of realistic mock HST/COS datasets, each comprising 34 quasar sightlines, and forward model the noise and resolution to match the real data. We use this large ensemble of mocks to extensively test our inference and empirically demonstrate that our posterior distributions are robust. Our analysis shows that by applying our new approach to existing Ly$\alpha$ forest spectra at $z\simeq 0.1$, one can measure the thermal and ionization state of the IGM with very high precision ($\sigma_{\log T_0} \sim 0.08$ dex, $\sigma_\gamma \sim 0.06$, and $\sigma_{\log \Gamma_{\rm HI}} \sim 0.07$ dex).


title: Simultaneous imaging of dopants and free charge carriers by STEM-EELS

abstract: Doping inhomogeneities in solids are not uncommon, but their microscopic observation and understanding are limited due to the lack of bulk-sensitive experimental techniques with high-enough spatial and spectral resolution. Here, we demonstrate nanoscale imaging of both dopants and free charge carriers in La-doped BaSnO3 (BLSO) using high-resolution electron energy-loss spectroscopy (EELS). By analyzing both high- and low-energy excitations in EELS, we reveal chemical and electronic inhomogeneities within a single BLSO nanocrystal. The inhomogeneous doping leads to distinctive localized infrared surface plasmons, including a novel plasmon mode that is highly confined between high- and low-doping regions. We further quantify the carrier density, effective mass, and dopant activation percentage from EELS data and transport measurements on the bulk single crystals of BLSO. These results represent a unique way of studying heterogeneities in solids, understanding structure-property relationships at the nanoscale, and opening the way to leveraging nanoscale doping texture in the design of nanophotonic devices.


title: A Phase-Field Study on the Effects of Nanoparticles on Solidification and Grain Growth

abstract: Nanoparticle reinforced alloys offer the potential of high strength, high temperature alloys. While promising, during rapid solidification processes, alloys suffer from nanoparticle clustering, which can discount any strength benefit. An open-source phase-field model is developed using PRISMS-PF to explore the impact of nanoparticles and clustering on alloy solidification. Heterogenous nucleation and grain boundary pinning are explicitly included, and a wide range of nanoparticle area fractions and nucleation rates are modeled. At low area fractions less than 0.05, particle clustering increases grain size between 15-45% compared to a random distribution. Our quantitative analyses inform a modified Zener grain size relationship that not only depends on nanoparticle size and area fraction, but also on the nucleation rate. Grain size first drastically decreases before plateauing at higher nucleation rates. Our simulations reveal a strong preference of nanoparticles pinning grain boundaries. Pinning fraction increases rapidly with nucleation rate before saturating between 0.85-0.90. Across the range of area fractions and nucleation rates considered, the random and clustered grain sizes each collapse to a simple analytical expression that depends only on nanoparticle radius and pinning fraction. Comparisons against experimental data reveal the expressions deduced from our analyses fit reported grain sizes better than classic Zener analysis. A simple model of strength and cost tradeoffs indicates nanoparticles can be a cost-effective way to improve alloy strength.


title: A Model for the Tropical Cyclone Wind Field Response to Idealized Landfall

abstract: The impacts of a tropical cyclone after landfall depend not only on storm intensity but also on the size and structure of the wind field. Hence, a simple predictive model for the wind field after landfall has significant potential value. This work tests existing theory for wind structure and size over the ocean against idealized axisymmetric landfall experiments in which the surface beneath a mature storm is instantaneously dried and roughened individually or simultaneously. Structure theory captures the response of the low-level wind field to different types of idealized landfalls given intensity and size response. Storm size, modeled to follow the ratio of simulated time-dependent storm intensity to the Coriolis parameter vm(t)/f, can generally predict the transient response of the storm gale wind radii r34kt to inland surface forcings, particularly for at least moderate surface roughening regardless of the level of drying. Given knowledge of the intensity evolution, the above results combine to yield a theoretical model that can predict the full tangential wind field response to idealized landfalls. An example application to a real world landfall is provided to demonstrate its potential utility for risk applications.


title: An effective criterion for finite monodromy of $\ell$-adic sheaves

abstract: We provide an effective version of Katz' criterion for finiteness of the monodromy group of a lisse, pure of weight zero, $\ell$-adic sheaf on a normal variety over a finite field, depending on the numerical complexity of the sheaf


title: Enhanced power graphs of groups are weakly perfect

abstract: A graph is weakly perfect if its clique number and chromatic number are equal. We show that the enhanced power graph of a finite group $G$ is weakly perfect: its clique number and chromatic number are equal to the maximum order of an element of $G$. The proof requires a combinatorial lemma. We give some remarks about related graphs.


title: $pp$ scattering at the LHC with the lepton pair production and one proton tagging

abstract: Analytical formulas for the cross section of the reaction $pp \to p + \ell^+ \ell^- + X$ are presented. Fiducial cross sections are compared with those measured recently by the ATLAS collaboration.


title: Streaming complexity of CSPs with randomly ordered constraints

abstract: We initiate a study of the streaming complexity of constraint satisfaction problems (CSPs) when the constraints arrive in a random order. We show that there exists a CSP, namely $\textsf{Max-DICUT}$, for which random ordering makes a provable difference. Whereas a $4/9 \approx 0.445$ approximation of $\textsf{DICUT}$ requires $\Omega(\sqrt{n})$ space with adversarial ordering, we show that with random ordering of constraints there exists a $0.48$-approximation algorithm that only needs $O(\log n)$ space. We also give new algorithms for $\textsf{Max-DICUT}$ in variants of the adversarial ordering setting. Specifically, we give a two-pass $O(\log n)$ space $0.48$-approximation algorithm for general graphs and a single-pass $\tilde{O}(\sqrt{n})$ space $0.48$-approximation algorithm for bounded degree graphs. On the negative side, we prove that CSPs where the satisfying assignments of the constraints support a one-wise independent distribution require $\Omega(\sqrt{n})$-space for any non-trivial approximation, even when the constraints are randomly ordered. This was previously known only for adversarially ordered constraints. Extending the results to randomly ordered constraints requires switching the hard instances from a union of random matchings to simple Erd\"os-Renyi random (hyper)graphs and extending tools that can perform Fourier analysis on such instances. The only CSP to have been considered previously with random ordering is $\textsf{Max-CUT}$ where the ordering is not known to change the approximability. Specifically it is known to be as hard to approximate with random ordering as with adversarial ordering, for $o(\sqrt{n})$ space algorithms. Our results show a richer variety of possibilities and motivate further study of CSPs with randomly ordered constraints.


title: Reduced Order Model for Chemical Kinetics: A case study with Primordial Chemical Network

abstract: Chemical kinetics plays an important role in governing the thermal evolution in reactive flows problems. The possible interactions between chemical species increase drastically with the number of species considered in the system. Various ways have been proposed before to simplify chemical networks with an aim to reduce the computational complexity of the chemical network. These techniques oftentimes require domain-knowledge experts to handcraftedly identify important reaction pathways and possible simplifications. Here, we propose a combination of autoencoder and neural ordinary differential equation to model the temporal evolution of chemical kinetics in a reduced subspace. We demonstrated that our model has achieved a close-to 10-fold speed-up compared to commonly used astro-chemistry solver for a 9-species primordial network, while maintaining 1 percent accuracy across a wide-range of density and temperature.


title: Case study on quantum convolutional neural network scalability

abstract: One of the crucial tasks in computer science is the processing time reduction of various data types, i.e., images, which is important for different fields -- from medicine and logistics to virtual shopping. Compared to classical computers, quantum computers are capable of parallel data processing, which reduces the data processing time. This quality of quantum computers inspired intensive research of the potential of quantum technologies applicability to real-life tasks. Some progress has already revealed on a smaller volumes of the input data. In this research effort, I aimed to increase the amount of input data (I used images from 2 x 2 to 8 x 8), while reducing the processing time, by way of skipping intermediate measurement steps. The hypothesis was that, for increased input data, the omitting of intermediate measurement steps after each quantum convolution layer will improve output metric results and accelerate data processing. To test the hypothesis, I performed experiments to chose the best activation function and its derivative in each network. The hypothesis was partly confirmed in terms of output mean squared error (MSE) -- it dropped from 0.25 in the result of classical convolutional neural network (CNN) training to 0.23 in the result of quantum convolutional neural network (QCNN) training. In terms of the training time, however, which was 1.5 minutes for CNN and 4 hours 37 minutes in the least lengthy training iteration, the hypothesis was rejected.


title: The Effects of Radiative Feedback and Supernova Induced Turbulence on Early Galaxies

abstract: The recently launched James Webb Space Telescope promises unparalleled advances in our understanding of the first stars and galaxies, but realizing this potential requires cosmological simulations that capture the key physical processes that affected these objects. Here we show that radiative transfer and subgrid turbulent mixing are two such processes. By comparing simulations with and without radiative transfer but with exactly the same physical parameters and subgrid turbulent mixing model, we show that tracking radiative transfer suppresses the Population III (Pop III) star formation density by a factor of approximately 4. In both simulations, $\gtrsim 90\%$ of Pop III stars are found in the unresolved pristine regions tracked by our subgrid model, which does a better job at modeling the regions surrounding proto-galaxy cores where metals from supernovae take tens of Myrs to mix thoroughly. At the same time, radiative transfer suppresses Pop III star formation, via the development of ionized bubbles that slows gas accretion in these regions, and it results in compact high-redshift galaxies that are surrounded by isolated low mass satellites. Thus turbulent mixing and radiative transfer are both essential processes that must be included to accurately model the morphology, composition, and growth of primordial galaxies.


title: Audio-guided Album Cover Art Generation with Genetic Algorithms

abstract: Over 60,000 songs are released on Spotify every day, and the competition for the listener's attention is immense. In that regard, the importance of captivating and inviting cover art cannot be underestimated, because it is deeply entangled with a song's character and the artist's identity, and remains one of the most important gateways to lead people to discover music. However, designing cover art is a highly creative, lengthy and sometimes expensive process that can be daunting, especially for non-professional artists. For this reason, we propose a novel deep-learning framework to generate cover art guided by audio features. Inspired by VQGAN-CLIP, our approach is highly flexible because individual components can easily be replaced without the need for any retraining. This paper outlines the architectural details of our models and discusses the optimization challenges that emerge from them. More specifically, we will exploit genetic algorithms to overcome bad local minima and adversarial examples. We find that our framework can generate suitable cover art for most genres, and that the visual features adapt themselves to audio feature changes. Given these results, we believe that our framework paves the road for extensions and more advanced applications in audio-guided visual generation tasks.


title: UV Spectropolarimetry with Polstar: Massive Star Binary Colliding Winds

abstract: The winds of massive stars are important for their direct impact on the interstellar medium, and for their influence on the final state of a star prior to it exploding as a supernova. However, the dynamics of these winds is understood primarily via their illumination from a single central source. The Doppler shift seen in resonance lines is a useful tool for inferring these dynamics, but the mapping from that Doppler shift to the radial distance from the source is ambiguous. Binary systems can reduce this ambiguity by providing a second light source at a known radius in the wind, seen from orbitally modulated directions. From the nature of the collision between the winds, a massive companion also provides unique additional information about wind momentum fluxes. Since massive stars are strong ultraviolet (UV) sources, and UV resonance line opacity in the wind is strong, UV instruments with a high resolution spectroscopic capability are essential for extracting this dynamical information. Polarimetric capability also helps to further resolve ambiguities in aspects of the wind geometry that are not axisymmetric about the line of sight, because of its unique access to scattering direction information. We review how the proposed MIDEX-scale mission Polstar can use UV spectropolarimetric observations to critically constrain the physics of colliding winds, and hence radiatively-driven winds in general. We propose a sample of 20 binary targets, capitalizing on this unique combination of illumination by companion starlight, and collision with a companion wind, to probe wind attributes over a range in wind strengths. Of particular interest is the hypothesis that the radial distribution of the wind acceleration is altered significantly, when the radiative transfer within the winds becomes optically thick to resonance scattering in multiple overlapping UV lines.


title: Pantheon+ constraints on dark energy and modified gravity: An evidence of dynamical dark energy

abstract: We use the latest Type Ia supernovae sample Pantheon+ to explore new physics on cosmic scales. Specifically, in light of this new sample, we constrain the interacting dark energy and Hu-Sawicki $f(R)$ gravity models and employ the model-independent Gaussian processes to investigate whether there is an evidence of dark energy evolution. We find that Pantheon+ alone just gives weak constraints. However, a data combination of Pantheon+ with cosmic microwave background, baryon acoustic oscillations and cosmic chronometers gives a strong constraint on the modified matter expansion rate $\epsilon=0.048\pm0.026$, which indicates that the momentum may transfer from dark energy to dark matter in the dark sector of the universe at the $1.85\,\sigma$ confidence level. In the meantime, we obtain a very tight constraint on the deviation from general relativity $\log_{10} f_{R0}< -6.32$ at the $2\,\sigma$ confidence level. Interestingly, when combining Pantheon+ with cosmic chronometers and cosmic microwave background data, we find a quintessence-like dark energy signal beyond the $2\,\sigma$ confidence level in the redshift range $z\in[0.70,1.05]$. This implies the nature of dark energy may actually be dynamical.

categories: astro-ph.CO astro-ph.HE gr-qc hep-ph hep-th


title: Estimating Emotion Contagion on Social Media via Localized Diffusion in Dynamic Graphs

abstract: We present a computational approach for estimating emotion contagion on social media networks. Built on a foundation of psychology literature, our approach estimates the degree to which the perceivers' emotional states (positive or negative) start to match those of the expressors, based on the latter's content. We use a combination of deep learning and social network analysis to model emotion contagion as a diffusion process in dynamic social network graphs, taking into consideration key aspects like causality, homophily, and interference. We evaluate our approach on user behavior data obtained from a popular social media platform for sharing short videos. We analyze the behavior of 48 users over a span of 8 weeks (over 200k audio-visual short posts analyzed) and estimate how contagious the users with whom they engage with are on social media. As per the theory of diffusion, we account for the videos a user watches during this time (inflow) and the daily engagements; liking, sharing, downloading or creating new videos (outflow) to estimate contagion. To validate our approach and analysis, we obtain human feedback on these 48 social media platform users with an online study by collecting responses of about 150 participants. We report users who interact with more number of creators on the platform are 12% less prone to contagion, and those who consume more content of `negative' sentiment are 23% more prone to contagion. We will publicly release our code upon acceptance.


title: K-level Reasoning for Zero-Shot Coordination in Hanabi

abstract: The standard problem setting in cooperative multi-agent settings is self-play (SP), where the goal is to train a team of agents that works well together. However, optimal SP policies commonly contain arbitrary conventions ("handshakes") and are not compatible with other, independently trained agents or humans. This latter desiderata was recently formalized by Hu et al. 2020 as the zero-shot coordination (ZSC) setting and partially addressed with their Other-Play (OP) algorithm, which showed improved ZSC and human-AI performance in the card game Hanabi. OP assumes access to the symmetries of the environment and prevents agents from breaking these in a mutually incompatible way during training. However, as the authors point out, discovering symmetries for a given environment is a computationally hard problem. Instead, we show that through a simple adaption of k-level reasoning (KLR) Costa Gomes et al. 2006, synchronously training all levels, we can obtain competitive ZSC and ad-hoc teamplay performance in Hanabi, including when paired with a human-like proxy bot. We also introduce a new method, synchronous-k-level reasoning with a best response (SyKLRBR), which further improves performance on our synchronous KLR by co-training a best response.


title: Approximation of Optimal Control Surfaces for $2\times 2$ Skew-Symmetric Evolutionary Game Dynamics

abstract: In this paper we study the problem of approximating the general solution to an optimal control problem whose dynamics arise from a $2\times 2$ skew-symmetric evolutionary game with arbitrary initial condition. Our approach uses a Fourier approximation method and generalizes prior work in the use of orthogonal function approximation for optimal control. At the same time we cast the fitting problem in the context of a non-standard feedforward neural network and derive the back-propagation operator in this context. An example of the efficacy of this approach is provided and generalizations are discussed.


title: Machine Learning to Generate Adjustable Dose Distributions in Head-and-Neck Cancer Radiation Therapy

abstract: In this work, we propose a Machine Learning model that generates an adjustable 3D dose distribution for external beam radiation therapy for head-and-neck cancer treatments. In contrast to existing Machine Learning methods that provide a single model, we create pairs of models for each organ-at-risk, namely lower-extreme and upper-extreme models. These model pairs for an organ-at-risk propose doses that give lower and higher doses to that organ-at-risk, while also encapsulating the dose trade-off to other organs-at-risk. By weighting and combining the model pairs for all organs-at-risk, we are able to dynamically create adjustable dose distributions that can be used, in real-time, to move doses between organs-at-risk, thereby customizing the dose distribution to the needs of a particular patient. We leverage a key observation that the training data set inherently contains the clinical trade-offs. We show that the adjustable distributions are able to provide reasonable clinical dose latitude in the trade-off of doses between organs-at-risk.


title: Linear arboricity of degenerate graphs

abstract: A linear forest is a union of vertex-disjoint paths, and the linear arboricity of a graph $G$, denoted by $\operatorname{la}(G)$, is the minimum number of linear forests needed to partition the edge set of $G$. Clearly, $\operatorname{la}(G) \ge \lceil\Delta(G)/2\rceil$ for a graph $G$ with maximum degree $\Delta(G)$. On the other hand, the Linear Arboricity Conjecture due to Akiyama, Exoo, and Harary from 1981 asserts that $\operatorname{la}(G) \leq \lceil(\Delta(G)+1) / 2\rceil$ for every graph $ G $. This conjecture has been verified for planar graphs and graphs whose maximum degree is at most $ 6 $, or is equal to $ 8 $ or $ 10 $. Given a positive integer $k$, a graph $G$ is $k$-degenerate if it can be reduced to a trivial graph by successive removal of vertices with degree at most $k$. We prove that for any $k$-degenerate graph $G$, $\operatorname{la}(G) = \lceil\Delta(G)/2 \rceil$ provided $\Delta(G) \ge 2k^2 -k$.


title: Beyond UVJ: Color Selection of Galaxies in the JWST Era

abstract: We present a new rest-frame color-color selection method using "synthetic $u_s-g_s$ and $g_s-i_s$'', $(ugi)_s$ colors to identify star-forming and quiescent galaxies. Our method is similar to the widely-used $U-V$ versus $V-J$ ($UVJ$) diagram. However, $UVJ$ suffers known systematics. Spectroscopic campaigns have shown that $UVJ$-selected quiescent samples at $z \gtrsim 3$ include $\sim 10-30\%$ contamination from galaxies with dust-obscured star formation and strong emission lines. Moreover, at $z>3$, $UVJ$ colors are extrapolated because the rest-frame J-band shifts beyond the coverage of the deepest bandpasses at $< 5~\mu m$ (typically $Spitzer$/IRAC 4.5 $\mu m$ or future $JWST$/NIRCam observations). We demonstrate that $(ugi)_s$ offers improvements to $UVJ$ at $z>3$, and can be applied to galaxies in the $JWST$ era. We apply $(ugi)_s$ selection to galaxies at $0.5<z<6$ from the (observed) 3D-HST and UltraVISTA catalogs, and to the (simulated) JAGUAR catalogs. We show that extrapolation can affect $(V-J)_0$ color by up to 1 magnitude, but changes $(u_s-g_s)_0$ color by $\leq$ 0.2 mag, even at $z\simeq 6$. While $(ugi)_s$-selected quiescent samples are comparable to $UVJ$ in completeness (both achieve $\sim$85-90% at $z=3-3.5$), $(ugi)_s$ reduces contamination in quiescent samples by nearly a factor of two, from $\simeq$35% to $\simeq$17% at $z=3$, and from $\simeq $60% to $\simeq $33% at $z=6$. This leads to improvements in the true-to-false-positive ratio (TP/FP), where we find TP/FP $\gtrsim$ 2.2 for $(ugi)_s$ at $z \simeq 3.5 - 6$, compared to TP/FP $<$ 1 for $UVJ$-selected samples. This indicates that contaminants will outnumber true quiescent galaxies in $UVJ$ at these redshifts, while $(ugi)_s$ will provide higher-fidelity samples.


title: An ALMA 1.3 millimeter Search for Debris Disks around Solar-type Stars in the Pleiades

abstract: Millimeter emission from debris disks around stars of different ages provides constraints on the collisional evolution of planetesimals. We present ALMA 1.3 millimeter observations of a sample of 76 Solar-type stars in the ~115 Myr old Pleiades star cluster. These ALMA observations complement previous infrared observations of this sample by providing sensitivity to emission from circumstellar dust at lower temperatures, corresponding to debris at radii comparable to the Kuiper Belt and beyond. The observations obtain a beam size of 1.5 arcsec (200 au) and a median rms noise of 54 mircoJy/beam, which corresponds to a fractional luminosity $L_{dust}/L_{star} \sim 10^{-4}$ for 40 K dust for a typical star in the sample. The ALMA images show no significant detections of the targeted stars. We interpret these limits in the context of a steady-state collisional cascade model for debris disk evolution that provides a good description of observations of the field population near the Sun but is not well-calibrated on younger populations.The ALMA non-detections of the Pleiades systems are compatible with the disk flux predictions of this model. We find no high fractional luminosity outliers from these ALMA data that could be associated with enhanced collisions resulting from activity not accounted for by steady-state evolution. However, we note that two systems (HII 1132 and HD 22680) show 24 micron excess much higher than the predictions of this model, perhaps due to unusually high dust production from dynamical events involving planets.


title: Anyonic statistics revealed by the Hong-Ou-Mandel dip for fractional excitations

abstract: The fractional quantum Hall effect (FQHE) is known to host anyons, quasiparticles whose statistics is intermediate between bosonic and fermionic. We show here that Hong-Ou-Mandel (HOM) interferences between excitations created by narrow voltage pulses on the edge states of a FQHE system at low temperature show a direct signature of anyonic statistics. The width of the HOM dip is universally fixed by the thermal time scale, independently of the intrinsic width of the excited fractional wavepackets. This universal width can be related to the anyonic braiding of the incoming excitations with thermal fluctuations created at the quantum point contact. We show that this effect could be realistically observed with periodic trains of narrow voltage pulses using current experimental techniques.


title: Deep Image Clustering with Contrastive Learning and Multi-scale Graph Convolutional Networks

abstract: Deep clustering has shown its promising capability in joint representation learning and clustering via deep neural networks. Despite the significant progress, the existing deep clustering works mostly utilize some distribution-based clustering loss, lacking the ability to unify representation learning and multi-scale structure learning. To address this, this paper presents a new deep clustering approach termed image clustering with contrastive learning and multi-scale graph convolutional networks (IcicleGCN), which bridges the gap between convolutional neural network (CNN) and graph convolutional network (GCN) as well as the gap between contrastive learning and multi-scale structure learning for the deep clustering task. Our framework consists of four main modules, namely, the CNN-based backbone, the Instance Similarity Module (ISM), the Joint Cluster Structure Learning and Instance reconstruction Module (JC-SLIM), and the Multi-scale GCN module (M-GCN). Specifically, the backbone network with two weight-sharing views is utilized to learn the representations for the two augmented samples (from each image). The learned representations are then fed to ISM and JC-SLIM for joint instance-level and cluster-level contrastive learning, respectively, during which an auto-encoder in JC-SLIM is also pretrained to serve as a bridge to the M-GCN module. Further, to enforce multi-scale neighborhood structure learning, two streams of GCNs and the auto-encoder are simultaneously trained via (i) the layer-wise interaction with representation fusion and (ii) the joint self-adaptive learning. Experiments on multiple image datasets demonstrate the superior clustering performance of IcicleGCN over the state-of-the-art. The code is available at https://github.com/xuyuankun631/IcicleGCN.


title: Causal Graphs Underlying Generative Models: Path to Learning with Limited Data

abstract: Training generative models that capture rich semantics of the data and interpreting the latent representations encoded by such models are very important problems in unsupervised learning. In this work, we provide a simple algorithm that relies on perturbation experiments on latent codes of a pre-trained generative autoencoder to uncover a causal graph that is implied by the generative model. We leverage pre-trained attribute classifiers and perform perturbation experiments to check for influence of a given latent variable on a subset of attributes. Given this, we show that one can fit an effective causal graph that models a structural equation model between latent codes taken as exogenous variables and attributes taken as observed variables. One interesting aspect is that a single latent variable controls multiple overlapping subsets of attributes unlike conventional approach that tries to impose full independence. Using a pre-trained RNN-based generative autoencoder trained on a dataset of peptide sequences, we demonstrate that the learnt causal graph from our algorithm between various attributes and latent codes can be used to predict a specific property for sequences which are unseen. We compare prediction models trained on either all available attributes or only the ones in the Markov blanket and empirically show that in both the unsupervised and supervised regimes, typically, using the predictor that relies on Markov blanket attributes generalizes better for out-of-distribution sequences.


title: Charm lifetime measurements at Belle II

abstract: Upgrades at the Belle II experiment and the SuperKEKB asymmetric-energy electron-positron collider enable precise measurements of particle decays. Even with early data, Belle II has made several world-leading measurements of particle lifetimes, which are useful as tests of effective models used in searches for physics beyond the standard model. Especially for charm hadrons, these effective models depend on careful consideration of non-perturbative effects to give an adequate description of lifetimes. Recent measurements at other experiments have stimulated particular interest in charm baryons, including several measurements that indicate the hierarchy of lifetimes is different than once believed. The measurements of the $D^+$, $D^0$, and $\Lambda_c^+$ lifetimes at Belle II are consistent with previous measurements, but with improved precision. These measurements will improve the world-average values, provide refined tests for effective models, and serve as benchmarks for future measurements.


title: Vesicle formation induced by thermal fluctuations

abstract: The process of fission and vesicle formation depends on the geometry of the membrane that will split. For instance, a flat surface finds it difficult to form vesicles because of the lack of curved regions where to start the process. Here we show that vesicle formation can be promoted by temperature, by using a membrane phase field model with Gaussian curvature. We find a phase transition between fluctuating and vesiculation phases that depends on temperature, spontaneous curvature, and the ratio between bending and Gaussian moduli. We analysed the energy dynamical behaviour of these processes and found that the main driving ingredient is the Gaussian energy term, although the curvature energy term usually helps with the process as well. We also found that the chemical potential can be used to investigate the temperature of the system. Finally we address how temperature changes the condition for spontaneous vesiculation for all geometries, making it happen in a wider range of values of the Gaussian modulus.


title: Characterizing and Optimizing End-to-End Systems for Private Inference

abstract: In two-party machine learning prediction services, the client's goal is to query a remote server's trained machine learning model to perform neural network inference in some application domain. However, sensitive information can be obtained during this process by either the client or the server, leading to potential collection, unauthorized secondary use, and inappropriate access to personal information. These security concerns have given rise to Private Inference (PI), in which both the client's personal data and the server's trained model are kept confidential. State-of-the-art PI protocols consist of a pre-processing or offline phase and an online phase that combine several cryptographic primitives: Homomorphic Encryption (HE), Secret Sharing (SS), Garbled Circuits (GC), and Oblivious Transfer (OT). Despite the need and recent performance improvements, PI remains largely arcane today and is too slow for practical use. This paper addresses PI's shortcomings with a detailed characterization of a standard high-performance protocol to build foundational knowledge and intuition in the systems community. Our characterization pinpoints all sources of inefficiency -- compute, communication, and storage. In contrast to prior work, we consider inference request arrival rates rather than studying individual inferences in isolation and we find that the pre-processing phase cannot be ignored and is often incurred online as there is insufficient downtime to hide pre-compute latency. Finally, we leverage insights from our characterization and propose three optimizations to address the storage (Client-Garbler), computation (layer-parallel HE), and communication (wireless slot allocation) overheads. Compared to the state-of-the-art PI protocol, these optimizations provide a total PI speedup of 1.8$\times$ with the ability to sustain inference requests up to a 2.24$\times$ greater rate.


title: Optically Modulated Propulsion of Electric Field Powered Photoconducting Janus Particles

abstract: Herein we demonstrate the ability to optically tune the mobility of electrically powered Janus particles (JP) that are half coated with various Zinc Oxide (ZnO) semiconducting layers, i.e. polycrystalline, amorphous and amorphous with a SiO2 passivation layer. The ZnO semiconductor photo-response enables increase in its electrical conductivity with light having wavelengths of sufficient photon energy with respect to the semiconductor bandgap. This effect, termed optically modulated electrokinetic propulsion (OMEP), can be harnessed to increase the contrast in polarizability between the dielectric and semiconducting hemispheres, which in turn, results in an increased electrokinetic mobility. The addition of optical activation to the electrical field enables an additional degree of control of JP mobility. We also demonstrate optical control of collective behavior and particle-particle interactions for dense semi-conducting Janus particle populations.


title: Rabinowitz Floer homology of negative line bundles and Floer Gysin sequence

abstract: This article is concerned with the Rabinowitz Floer homology of negative line bundles. We construct a refined version of Rabinowitz Floer homology and study its properties. In particular, we build a Gysin-type long exact sequence for this new invariant and discuss an application to the orderability problem for prequantization spaces. We also construct a short exact sequence for the ordinary Rabinowitz Floer homology and provide computational results.


title: Contrastive Adapters for Foundation Model Group Robustness

abstract: While large pretrained foundation models (FMs) have shown remarkable zero-shot classification robustness to dataset-level distribution shifts, their robustness to subpopulation or group shifts is relatively underexplored. We study this problem, and find that FMs such as CLIP may not be robust to various group shifts. Across 9 robustness benchmarks, zero-shot classification with their embeddings results in gaps of up to 80.7 percentage points (pp) between average and worst-group accuracy. Unfortunately, existing methods to improve robustness require retraining, which can be prohibitively expensive on large foundation models. We also find that efficient ways to improve model inference (e.g., via adapters, lightweight networks with FM embeddings as inputs) do not consistently improve and can sometimes hurt group robustness compared to zero-shot (e.g., increasing the accuracy gap by 50.1 pp on CelebA). We thus develop an adapter training strategy to effectively and efficiently improve FM group robustness. Our motivating observation is that while poor robustness results from groups in the same class being embedded far apart in the foundation model "embedding space," standard adapter training may not bring these points closer together. We thus propose contrastive adapting, which trains adapters with contrastive learning to bring sample embeddings close to both their ground-truth class embeddings and other sample embeddings in the same class. Across the 9 benchmarks, our approach consistently improves group robustness, raising worst-group accuracy by 8.5 to 56.0 pp over zero-shot. Our approach is also efficient, doing so without any FM finetuning and only a fixed set of frozen FM embeddings. On benchmarks such as Waterbirds and CelebA, this leads to worst-group accuracy comparable to state-of-the-art methods that retrain entire models, while only training $\leq$1% of the model parameters.


title: Passivity-based control of underactuated mechanical systems with Coulomb friction: Application to earthquake prevention

abstract: Passivity property gives a sense of energy balance. The classical definitions and theorems of passivity in dynamical systems require time invariance and locally Lipschitz functions. However, these conditions are not met in many systems. A characteristic example is nonautonomous and discontinuous systems due to presence of Coulomb friction. This paper presents an extended result for the negative feedback connection of two passive nonautonomous systems with set-valued right-hand side based on an invariance-like principle. Such extension is the base of a structural passivity-based control synthesis for underactuated mechanical systems with Coulomb friction. The first step consists in designing the control able to restore the passivity in the considered friction law, achieving stabilization of the system trajectories to a domain with zero velocities. Then, an integral action is included to improve the latter result and perform a tracking over a constant reference (regulation). At last, the control is designed considering dynamics in the actuation. These control objectives are obtained using fewer control inputs than degrees of freedom, as a result of the underactuated nature of the plant. The presented control strategy is implemented in an earthquake prevention scenario, where a mature seismogenic fault represents the considered frictional underactuated mechanical system. Simulations are performed to show how the seismic energy can be slowly dissipated by tracking a slow reference, thanks to fluid injection far from the fault, accounting also for the slow dynamics of the fluid's diffusion.


title: Holographic Meissner Effect

abstract: The holographic superconductor is the holographic dual of superconductivity, but there is no Meissner effect in the standard holographic superconductor. This is because the boundary Maxwell field is added as an external source and is not dynamical. We show the Meissner effect analytically by imposing the semiclassical Maxwell equation on the AdS boundary. Unlike in the Ginzburg-Landau (GL) theory, the extreme Type I limit cannot be reached even in the $e\to\infty$ limit where $e$ is the $U(1)$ coupling of the boundary Maxwell field. This is due to the bound current which is present even in the pure bulk Maxwell theory. In the bulk 5-dimensional case, the GL parameter and the dual GL theory are obtained analytically for the order parameter of scaling dimension 2.


title: Learning Embedded Representation of the Stock Correlation Matrix using Graph Machine Learning

abstract: Understanding non-linear relationships among financial instruments has various applications in investment processes ranging from risk management, portfolio construction and trading strategies. Here, we focus on interconnectedness among stocks based on their correlation matrix which we represent as a network with the nodes representing individual stocks and the weighted links between pairs of nodes representing the corresponding pair-wise correlation coefficients. The traditional network science techniques, which are extensively utilized in financial literature, require handcrafted features such as centrality measures to understand such correlation networks. However, manually enlisting all such handcrafted features may quickly turn out to be a daunting task. Instead, we propose a new approach for studying nuances and relationships within the correlation network in an algorithmic way using a graph machine learning algorithm called Node2Vec. In particular, the algorithm compresses the network into a lower dimensional continuous space, called an embedding, where pairs of nodes that are identified as similar by the algorithm are placed closer to each other. By using log returns of S&P 500 stock data, we show that our proposed algorithm can learn such an embedding from its correlation network. We define various domain specific quantitative (and objective) and qualitative metrics that are inspired by metrics used in the field of Natural Language Processing (NLP) to evaluate the embeddings in order to identify the optimal one. Further, we discuss various applications of the embeddings in investment management.


title: Exploring subtropical stratocumulus multiple equilibria using a mixed-layer model

abstract: Stratocumulus clouds cover about a fifth of Earths surface, and due to their albedo and low-latitude location, they have a strong effect on Earths radiation budget. Previous studies using Large Eddy Simulations have shown that multiple equilibria (both cloud-covered and cloud-free states) exist as a function of fixed-SST, with relevance to equatorward advected air masses. Multiple equilibria have also been found as a function of atmospheric CO2, with a subtropical SST nearly 10 K higher in the cloud-free state, and with suggested relevance to warm-climate dynamics. In this study, we use a mixed-layer model with an added surface energy balance and the ability to simulate both the cloud-covered and cloud-free states to study both types of multiple equilibria and the corresponding hysteresis. The model's simplicity and computational efficiency allow us to explore the mechanisms critical to the stratocumulus cloud instability and hysteresis as well as isolate key processes that allow for multiple equilibria via mechanism denial experiments not possible with a full-complexity model. For the hysteresis in fixed-SST, we find that decoupling can occur due to either enhanced entrainment warming or a reduction in cloud-top longwave cooling. The critical SST at which decoupling occurs is highly sensitive to precipitation and entrainment parameterizations. In the CO2 hysteresis, decoupling occurs in the simple model used even without the inclusion of an active SST and SST-cloud cover feedbacks, and the width of the hysteresis displays the same sensitivities as the fixed-SST case. Overall, the simple model analysis and results motivate further studies using higher complexity models.


title: Honda-Huang's work on contact convexity revisited

abstract: Following the overall strategy of the paper ``Convex hypersurfaces in contact topology" by Ko Honda and Yang Huang on contact convexity in high dimensions, we present a simplified proof of their main result.


title: Are generic dynamical properties stable under composition with rotations?

abstract: In this paper we provide a detailed topological and measure-theoretic study of Lebesgue measure-preserving circle maps that are rotated with inner and outer rotations which are independent of each other. In particular, we analyze the stability of the locally eventually onto and measure-theoretic mixing properties.


title: NASRec: Weight Sharing Neural Architecture Search for Recommender Systems

abstract: The rise of deep neural networks offers new opportunities in optimizing recommender systems. However, optimizing recommender systems using deep neural networks requires delicate architecture fabrication. We propose NASRec, a paradigm that trains a single supernet and efficiently produces abundant models/sub-architectures by weight sharing. To overcome the data multi-modality and architecture heterogeneity challenges in the recommendation domain, NASRec establishes a large supernet (i.e., search space) to search the full architectures. The supernet incorporates versatile choice of operators and dense connectivity to minimize human efforts for finding priors. The scale and heterogeneity in NASRec impose several challenges, such as training inefficiency, operator-imbalance, and degraded rank correlation. We tackle these challenges by proposing single-operator any-connection sampling, operator-balancing interaction modules, and post-training fine-tuning. Our crafted models, NASRecNet, show promising results on three Click-Through Rates (CTR) prediction benchmarks, indicating that NASRec outperforms both manually designed models and existing NAS methods with state-of-the-art performance. Our work is publicly available at https://github.com/facebookresearch/NasRec.


title: Frank Elastic Constants of Semiflexible Polymer Solutions

abstract: We derive the Frank elastic constants for nematic solutions of semiflexible polymers. We plot these results as a function of the coarse-grained Maier-Saupe quadrupole aligning strength and polymer stiffness ranging from rigid to highly flexible. The derivation uses the random phase approximation and combines the exact results for the statistics of a worm-like-chain with polymer field theory using a spherical harmonic basis. The results are evaluated using a numerical inverse Laplace transform. We present the results in terms of microscopic features such as hairpins and polymer ends so the trends can be understood independently from the derivation. Key findings are that for rigid polymers $K_{bend}>K_{splay}>K_{twist}$ while for flexible polymers $K_{splay}>K_{bend}>K_{twist}$. For rigid polymers, the Frank elastic constants grow with the polymer length. For flexible polymers the elastic constants grow with the persistence length, which becomes the characteristic length scale, with the exception of $K_{splay}$ at high alignment strengths which grows with polymer lengths due to the elimination of hairpins.


title: Current Trends in Deep Learning for Earth Observation: An Open-source Benchmark Arena for Image Classification

abstract: We present AiTLAS: Benchmark Arena -- an open-source benchmark suite for evaluating state-of-the-art deep learning approaches for image classification in Earth Observation (EO). To this end, we present a comprehensive comparative analysis of more than 500 models derived from ten different state-of-the-art architectures and compare them to a variety of multi-class and multi-label classification tasks from 22 datasets with different sizes and properties. In addition to models trained entirely on these datasets, we benchmark models trained in the context of transfer learning, leveraging pre-trained model variants, as it is typically performed in practice. All presented approaches are general and can be easily extended to many other remote sensing image classification tasks not considered in this study. To ensure reproducibility and facilitate better usability and further developments, all of the experimental resources including the trained models, model configurations, and processing details of the datasets (with their corresponding splits used for training and evaluating the models) are publicly available on the repository: https://github.com/biasvariancelabs/aitlas-arena


title: Queueing games with an endogenous number of machines

abstract: This paper studies queueing problems with an endogenous number of machines with and without an initial queue, the novelty being that coalitions not only choose how to queue, but also on how many machines. For a given problem, agents can (de)activate as many machines as they want, at a cost. After minimizing the total cost (processing costs and machine costs), we use a game theoretical approach to share to proceeds of this cooperation, and study the existence of stable allocations. First, we study queueing problems with an endogenous number of machines, and examine how to share the total cost. We provide an upper bound and a lower bound on the cost of a machine to guarantee the non-emptiness of the core (the set of stable allocations). Next, we study requeueing problems with an endogenous number of machines, where there is an existing queue. We examine how to share the cost savings compared to the initial situation, when optimally requeueing/changing the number of machines. Although, in general, stable allocation may not exist, we guarantee the existence of stable allocations when all machines are considered public goods, and we start with an initial schedule that might not have the optimal number of machines, but in which agents with large waiting costs are processed first.


title: Soft-phonon and charge-density-wave formation in nematic BaNi$_2$As$_2$

abstract: We use diffuse and inelastic x-ray scattering to study the formation of an incommensurate charge-density-wave (I-CDW) in BaNi$_2$As$_2$, a candidate system for charge-driven electronic nematicity. Intense diffuse scattering is observed around the modulation vector of the I-CDW, $Q_{I-CDW}$. It is already visible at room temperature and collapses into superstructure reflections in the long-range ordered state where a small orthorhombic distortion occurs. A clear dip in the dispersion of a low-energy transverse optical phonon mode is observed around $Q_{I-CDW}$. The phonon continuously softens upon cooling, ultimately driving the transition to the I-CDW state. The transverse character of the soft-phonon branch elucidates the complex pattern of the I-CDW satellites observed in the current and earlier studies and settles the debated unidirectional nature of the I-CDW. The phonon instability and its reciprocal space position is well captured by our ${ab}$ ${initio}$ calculations. These however indicate that neither Fermi surface nesting, nor enhanced momentum-dependent electron-phonon coupling can account for the I-CDW formation, demonstrating its unconventional nature.% and suggest alternative origins of the I-CDW.


title: Analyzing Planar Galactic Halo Distributions with Fuzzy/Cold Dark Matter Models

abstract: We perform a numerical comparison between the fuzzy dark matter model and the cold dark matter model, focusing on formation of satellite galaxy planes around massive galaxies. Such galactic dynamics with controlled initial subhalo configurations are investigated using GADGET2 for the cold dark matter and PyUltraLight for the fuzzy dark matter, respectively. We demonstrate that satellite galaxies in the fuzzy dark matter side have a tendency to form more flattened and corotating satellite systems than in the cold dark matter side mainly due to the dissipation by the gravitational cooling effect of the fuzzy dark matter. Our simulations with the fuzzy dark matter typically show the minor-to-major axis ratio $c/a$ of the satellite galaxy planes to be $0.21 \sim 0.30$; This well matches the current observed value for the Milky Way.


title: Accretion flows around exotic tidal wormholes I. Ray-tracing

abstract: This paper investigates the various spherically symmetric wormhole solutions in the presence of tidal forces and applies numerous methods, such as test particle orbital dynamics, ray-tracing and microlensing. We make the theoretical predictions on the test particle orbital motion around the tidal wormholes with the use of normalized by $\mathcal{L}^2$ effective potential. In order to obtain the ray-tracing images (of both geometrically thin and thick accretion disks, relativistic jets), we properly modify the open source $\texttt{GYOTO}$ code with python interface. We applied this techniques to probe the accretion flows nearby the Schwarzschild-like and charged Reissner-N\"ordstrom (RS) wormholes (we assumed both charged RS wormhole and special case with the vanishing electromagnetic charge, namely Damour-Solodukhin (DS) wormhole). It was shown that the photon sphere for Schwarzschild-like wormhole presents for both thin and thick accretion disks and even for the vanishing tidal forces. Moreover, it was observed that $r_{\mathrm{ph}}\to\infty$ as $\alpha\to\infty$, which constraints $\alpha$ parameter to be sufficiently small and positive in order to respect the EHT observations. On the other hand, for the case of RS wormhole, photon sphere radius shrinks as $\Lambda\to\infty$, as it was predicted by the effective potential. In addition to the accretion disks, we as well probe the relativistic jets around two wormhole solutions of our consideration. Finally, with the help of star bulb microlensing, we approximate the radius of the wormhole shadow and as we found out, for Schild WH, $R_{\mathrm{Sh}}\approx r_0$ for ZTF and grows linearly with $\alpha$. On the contrary, shadow radius for charged wormholes slowly decreases with the growing trend of DS parameter $\Lambda$.


title: A serendipity fully discrete div-div complex on polygonal meshes

abstract: In this work we address the reduction of face degrees of freedom (DOFs) for discrete elasticity complexes. Specifically, using serendipity techniques, we develop a reduced version of a recently introduced two-dimensional complex arising from traces of the three-dimensional elasticity complex. The keystone of the reduction process is a new estimate of symmetric tensor-valued polynomial fields in terms of boundary values, completed with suitable projections of internal values for higher degrees. We prove an extensive set of new results for the original complex and show that the reduced complex has the same homological and analytical properties as the original one. This paper also contains an appendix with proofs of general Poincar\'e--Korn-type inequalities for hybrid fields.


title: COOR-PLT: A hierarchical control model for coordinating adaptive platoons of connected and autonomous vehicles at signal-free intersections based on deep reinforcement learning

abstract: Platooning and coordination are two implementation strategies that are frequently proposed for traffic control of connected and autonomous vehicles (CAVs) at signal-free intersections instead of using conventional traffic signals. However, few studies have attempted to integrate both strategies to better facilitate the CAV control at signal-free intersections. To this end, this study proposes a hierarchical control model, named COOR-PLT, to coordinate adaptive CAV platoons at a signal-free intersection based on deep reinforcement learning (DRL). COOR-PLT has a two-layer framework. The first layer uses a centralized control strategy to form adaptive platoons. The optimal size of each platoon is determined by considering multiple objectives (i.e., efficiency, fairness and energy saving). The second layer employs a decentralized control strategy to coordinate multiple platoons passing through the intersection. Each platoon is labeled with coordinated status or independent status, upon which its passing priority is determined. As an efficient DRL algorithm, Deep Q-network (DQN) is adopted to determine platoon sizes and passing priorities respectively in the two layers. The model is validated and examined on the simulator Simulation of Urban Mobility (SUMO). The simulation results demonstrate that the model is able to: (1) achieve satisfactory convergence performances; (2) adaptively determine platoon size in response to varying traffic conditions; and (3) completely avoid deadlocks at the intersection. By comparison with other control methods, the model manifests its superiority of adopting adaptive platooning and DRL-based coordination strategies. Also, the model outperforms several state-of-the-art methods on reducing travel time and fuel consumption in different traffic conditions.


title: Characterization of sensitivity and responses of a 2-element prototype wavefront sensor for millimeter-wave adaptive optics attached to the Nobeyama 45 m telescope

abstract: We report the results of the performance characterization of a prototype wavefront sensor for millimetric adaptive optics (MAO) installed on the Nobeyama 45 m radio telescope. MAO is a key component to realize a future large-aperture submillimeter telescope, such as Large Submillimeter Telescope (LST) or Atacama Large Aperture Submillimeter Telescope (AtLAST). The difficulty of MAO is, however, real-time sensing of wavefront deformation with ~10 um accuracy across the aperture. Our wavefront sensor operating at 20 GHz measures the radio path length between a certain position of the primary mirror surface to the focal point where a 20 GHz coherent receiver is placed. With the 2-element prototype, we sampled two positions on the primary mirror surface (at radii of 5 m and 16 m) at a sampling rate of 10 Hz. Then an excess path length (EPL) between the two positions was obtained by differentiating the two optical paths. A power spectral density of the EPL shows three components: a low-frequency drift (1/f^n), oscillations, and a white noise. A comparison of EPL measurements under a variety of wind conditions suggests that the former two are likely induced by the wind load on the telescope structure. The power of the white noise corresponds to a 1sigma statistical error of 8 um in EPL measurements. The 8 um r.m.s. is significant with respect to the mirror surface accuracy required by the LST and AtLAST (~20-40 um r.m.s.), which demonstrates that our technique is also useful for the future large-aperture submillimeter telescopes.


title: Phenomenology of quantum eigenstates in mixed-type systems: lemon billiards with complex phase space structure

abstract: The boundary of the lemon billiards is defined by the intersection of two circles of equal unit radius with the distance 2B between their centers, as introduced by Heller and Tomsovic in Phys. Today 46 38 (1993). We study two classical and quantum lemon billiards, for the cases B = 0.1953, 0.083, which are mixed-type billiards with complex structure of phase space, without significant stickiness regions. A preliminary study of their spectra was published recently (Physics 1 1-14 (2021)). We calculate a great number $10^6$ of consecutive eigenstates and their Poincar\'e-Husimi (PH) functions, and analyze their localization properties by studying the entropy localization measure and the normalized inverse participation ratio. We also introduce an overlap index which measures the degree of the overlap of PH functions with classically regular and chaotic regions. We observe the existence of regular states associated with invariant tori and chaotic states associated with the classically chaotic regions, and also the mixed-type states. We show that in accordance with the Berry-Robnik picture and the principle of uniform semiclassical condensation of PH functions the relative fraction of mixed-type states decreases as a power law with increasing energy, thus in the strict semiclassical limit leaving only purely regular and chaotic states. Our approach offers a general very good phenomenological overview of the structural and localization properties of PH functions in quantum mixed-type Hamiltonian systems.


title: The Effect of Sideslip on Jackknife Limits During Low Speed Trailer Operation

abstract: Jackknifing refers to the serious situation where a vehicle-trailer system enters a jackknife state and the vehicle and trailer eventually collide if trailer operation is not corrected. This paper considers low speed trailer maneuvering typical of trailer backing where jackknife state limits can vary due to sideslip caused by physical interaction between the vehicle, trailer, and environment. Analysis of a kinematic model considering sideslip at the vehicle and trailer wheels indicates that vehicle-trailer systems should be divided into three categories based on the ratio of hitch length and trailer tongue length, each with distinct behaviors. The Long Trailer category may have no jackknifing state while the other two categories always have states leading to jackknifing. It is found that jackknife limits, which are the boundaries between the jackknifing state and the recoverable regions, can be divided into safe and unsafe limits, the latter of which must be avoided. Simulations and physical experiments support these results and provide insight about the implications of vehicle and trailer states with slip that lead to jackknifing. Simulations also demonstrate the benefit of considering these new slip-based jackknife limits in trailer backing control.


title: An experimental study on the settling velocity of inertial particles in different homogeneous isotropic turbulent flows

abstract: We propose an experimental study on the gravitational settling velocity of dense, sub-Kolmogorov inertial particles under different background turbulent flows. We report Phase Doppler Particle Analyzer measurements in a low-speed wind tunnel uniformly seeded with micrometer scale water droplets. Turbulence is generated with three different grids (two consisting on different active-grid protocols while the third is a regular static grid), allowing us to cover a very wide range of turbulence conditions in terms of Taylor-scale based Reynolds numbers ($Re_\lambda \in [30-520]$), Rouse numbers ($Ro \in [0-5]$) and volume fractions ($\phi_v \in[0.5\times10^{-5} - 2.0\times10^{-5}]$). We find, in agreement with previous works, that enhancement of the settling velocity occurs at low Rouse number, while hindering of the settling occurs at higher Rouse number for decreasing turbulence energy levels. The wide range of flow parameters explored allowed us to observe that enhancement decreases significantly with the Taylor Reynolds number and is significantly affected by the volume fraction $\phi_v$. We also studied the effect of large-scale forcing on settling velocity modification. The possibility of changing the inflow conditions by using different grids allowed us to test cases with fixed $Re_\lambda$ and turbulent intensity but different integral length scale. Finally, we assess the existence of secondary flows in the wind tunnel and their role on particle settling. This is achieved by characterising the settling velocity at two different positions, the centreline and close to the wall, with the same streamwise coordinate.


title: Arity Approximation of $\infty$-Operads

abstract: Let $\mathbb{E}_d$ denote the little discs operad for $1 \le d \le \infty$ and let $\mathcal{C}$ be an $\infty$-category all of whose mapping spaces are $n$-truncated. We prove that when considering $\mathbb{E}_d$-monoids in $\mathcal{C}$, all coherence diagrams of arity $>n+3$ are redundant. More generally, for an $\infty$-operad $\mathcal{O}$ we bound the arity of the relevant coherence diagrams in terms of the connectivity of certain operadic partition complexes associated to $\mathcal{O}$.


title: The propagation of air fingers into an elastic branching network

abstract: We study experimentally the propagation of an air finger through the Y-bifurcation of an elastic, liquid-filled Hele-Shaw channel, as a benchtop model of airway reopening. With channel compliance provided by an elastic upper boundary, we can impose collapsed channel configurations into which we inject air with constant volume-flux. We typically observe steady finger propagation in the main channel, which is lost ahead of the Y-bifurcation but subsequently recovered in the daughter channels. At low levels of initial collapse, steady finger shapes and bubble pressure in the daughter channels map onto those in the main channel, despite small differences in initial collapse in different parts of the Y-channel. However, at higher levels of initial collapse where the elastic sheet almost touches the bottom boundary of the channel, experimentally indistinguishable fingers in the main channel can lead to multiple states of reopening of the daughter channels. The downstream distance at which steady propagation is recovered in the daughter channels also varies considerably with injection flow rate and initial collapse because of a transition in the mechanics regulating finger propagation. We find that the characteristic time and length-scales of this recovery are largest in the regime where viscous and surface tension forces dominate at low flow rate and/or low initial collapse, and that they decrease towards a constant plateau reached in the limit where elastic and surface tension forces balance at high flow rate and/or high initial collapse. Our findings suggest that practical networks are unlikely to comprise long enough channels for steady state propagation to remain established.


title: Prediction of protein allosteric signalling pathways and functional residues through paths of optimised propensity

abstract: Allostery commonly refers to the mechanism that regulates protein activity through the binding of a molecule at a different, usually distal, site from the orthosteric site. The omnipresence of allosteric regulation in nature and its potential for drug design and screening render the study of allostery invaluable. Nevertheless, challenges remain as few computational methods are available to effectively predict allosteric sites, identify signalling pathways involved in allostery, or to aid with the design of suitable molecules targeting such sites. Recently, bond-to-bond propensity analysis has been shown successful at identifying allosteric sites for a large and diverse group of proteins from knowledge of the orthosteric sites and its ligands alone by using network analysis applied to energy-weighted atomistic protein graphs. To address the identification of signalling pathways, we propose here a method to compute and score paths of optimised propensity that link the orthosteric site with the identified allosteric sites, and identifies crucial residues that contribute to those paths. We showcase the approach with three well-studied allosteric proteins: h-Ras, caspase-1, and 3-phosphoinositide-dependent kinase-1 (PDK1). Key residues in both orthosteric and allosteric sites were identified and showed agreement with experimental results, and pivotal signalling residues along the pathway were also revealed, thus providing alternative targets for drug design. By using the computed path scores, we were also able to differentiate the activity of different allosteric modulators.


title: Lapse risk modelling in insurance: a Bayesian mixture approach

abstract: This paper focuses on modelling surrender time for policyholders in the context of life insurance. In this setup, a large lapse rate at the first months of a contract is often observed, with a decrease in this rate after some months. The modelling of the time to cancellation must account for this specific behaviour. Another stylised fact is that policies which are not cancelled in the study period are considered censored. To account for both censuring and heterogeneous lapse rates, this work assumes a Bayesian survival model with a mixture of regressions. The inference is based on data augmentation allowing for fast computations even for data sets of over a million clients. Moreover, scalable point estimation based on EM algorithm is also presented. An illustrative example emulates a typical behaviour for life insurance contracts and a simulated study investigates the properties of the proposed model. In particular, the observed censuring in the insurance context might be up to 50% of the data, which is very unusual for survival models in other fields such as epidemiology. This aspect is exploited in our simulated study.


title: Significant Digits of Primes in Subsets

abstract: Benford's Law describes the prevalence of small numbers as the leading digits of numbers in many sets of integers. We prove a variant of Benford's law for many positive-density subsets of the primes. This follows from a more general result over number fields.


title: Deformation of Nanowires and Nanotubes

abstract: In this article, we have investigated the consequences of the next to the leading order correction to the effective field theory of nanostructures. This has been done by analyzing the effects of deformed Heisenberg algebra on nanowires and nanotubes. We first deform the Schrodinger equation with cylindrical topology. Then specific solutions to the deformed Schrodinger equation with different boundary conditions are studied. These deformed solutions are used to investigate the consequences of the deformation on the energy of nanowires and nanotubes.


title: Explicit canonical heights for divisors relative to endomorphisms of $\mathbb{P}^N$

abstract: Given an endomorphism f of projective space, we exhibit explicit bounds on the difference between the naive height of a divisor and its canonical height relative to f.


title: Optimal Liouville theorem for a semilinear Ornstein-Uhlenbeck equation

abstract: The question of triviality of solutions of the semilinear Ornstein-Uhlenbeck equation, \[ \Delta w-\frac{1}{2} \langle x,\nabla w\rangle-\frac{\lambda}{p-1}w+|w|^{p-1}w=0, \] is considered. It is shown, that if $p>1$ is Sobolev subcritical or critical and $\lambda\leq 1$, then all bounded entire solutions are constant. Moreover, in the critical case, the same conclusion holds in the subclass of radial solutions provided that $n\geq 4$ and $\lambda \in \left[\frac{3 n}{2(n-1)},2\right]$.


title: Provably Adversarially Robust Nearest Prototype Classifiers

abstract: Nearest prototype classifiers (NPCs) assign to each input point the label of the nearest prototype with respect to a chosen distance metric. A direct advantage of NPCs is that the decisions are interpretable. Previous work could provide lower bounds on the minimal adversarial perturbation in the $\ell_p$-threat model when using the same $\ell_p$-distance for the NPCs. In this paper we provide a complete discussion on the complexity when using $\ell_p$-distances for decision and $\ell_q$-threat models for certification for $p,q \in \{1,2,\infty\}$. In particular we provide scalable algorithms for the \emph{exact} computation of the minimal adversarial perturbation when using $\ell_2$-distance and improved lower bounds in other cases. Using efficient improved lower bounds we train our Provably adversarially robust NPC (PNPC), for MNIST which have better $\ell_2$-robustness guarantees than neural networks. Additionally, we show up to our knowledge the first certification results w.r.t. to the LPIPS perceptual metric which has been argued to be a more realistic threat model for image classification than $\ell_p$-balls. Our PNPC has on CIFAR10 higher certified robust accuracy than the empirical robust accuracy reported in (Laidlaw et al., 2021). The code is available in our repository.


title: Sound Randomized Smoothing in Floating-Point Arithmetics

abstract: Randomized smoothing is sound when using infinite precision. However, we show that randomized smoothing is no longer sound for limited floating-point precision. We present a simple example where randomized smoothing certifies a radius of $1.26$ around a point, even though there is an adversarial example in the distance $0.8$ and extend this example further to provide false certificates for CIFAR10. We discuss the implicit assumptions of randomized smoothing and show that they do not apply to generic image classification models whose smoothed versions are commonly certified. In order to overcome this problem, we propose a sound approach to randomized smoothing when using floating-point precision with essentially equal speed and matching the certificates of the standard, unsound practice for standard classifiers tested so far. Our only assumption is that we have access to a fair coin.


title: High-resolution large-eddy simulations of simplified urban flows

abstract: High-fidelity large-eddy simulations of the flow around two rectangular obstacles are carried out at a Reynolds number of 10,000 based on the free-stream velocity and the obstacle height. The incoming flow is a developed turbulent boundary layer. Mean-velocity components, turbulence fluctuations, and the terms of the turbulent-kinetic-energy budget are analyzed for three flow regimes: skimming flow, wake interference, and isolated roughness. Three regions are identified where the flow undergoes the most significant changes: the first obstacle's wake, the region in front of the second obstacle, and that around the second obstacle. In the skimming-flow case, turbulence activity in the cavity between the obstacles is limited and mainly occurs in a small region in front of the second obstacle. In the wake-interference case, there is a strong interaction between the free-stream flow that penetrates the cavity and the wake of the first obstacle. This interaction results in more intense turbulent fluctuations between the obstacles. In the isolated-roughness case, the wake of the first obstacle is in good agreement with that of an isolated obstacle. Separation bubbles with strong turbulent fluctuations appear around the second obstacle.


title: Computing Optimal Kernels in Two Dimensions

abstract: Let $P$ be a set of $n$ points in $\Re^2$. For a parameter $\varepsilon\in (0,1)$, a subset $C\subseteq P$ is an \emph{$\varepsilon$-kernel} of $P$ if the projection of the convex hull of $C$ approximates that of $P$ within $(1-\varepsilon)$-factor in every direction. The set $C$ is a \emph{weak $\varepsilon$-kernel} of $P$ if its directional width approximates that of $P$ in every direction. Let $\mathsf{k}_{\varepsilon}(P)$ (resp.\ $\mathsf{k}^{\mathsf{w}}_{\varepsilon}(P)$) denote the minimum-size of an $\varepsilon$-kernel (resp. weak $\varepsilon$-kernel) of $P$. We present an $O(n\mathsf{k}_{\varepsilon}(P)\log n)$-time algorithm for computing an $\varepsilon$-kernel of $P$ of size $\mathsf{k}_{\varepsilon}(P)$, and an $O(n^2\log n)$-time algorithm for computing a weak $\varepsilon$-kernel of $P$ of size ${\mathsf{k}}^{\mathsf{w}}_{\varepsilon}(P)$. We also present a fast algorithm for the Hausdorff variant of this problem. In addition, we introduce the notion of \emph{$\varepsilon$-core}, a convex polygon lying inside $\mathsf{ch}(P)$, prove that it is a good approximation of the optimal $\varepsilon$-kernel, present an efficient algorithm for computing it, and use it to compute an $\varepsilon$-kernel of small size.


title: Attention, Filling in The Gaps for Generalization in Routing Problems

abstract: Machine Learning (ML) methods have become a useful tool for tackling vehicle routing problems, either in combination with popular heuristics or as standalone models. However, current methods suffer from poor generalization when tackling problems of different sizes or different distributions. As a result, ML in vehicle routing has witnessed an expansion phase with new methodologies being created for particular problem instances that become infeasible at larger problem sizes. This paper aims at encouraging the consolidation of the field through understanding and improving current existing models, namely the attention model by Kool et al. We identify two discrepancy categories for VRP generalization. The first is based on the differences that are inherent to the problems themselves, and the second relates to architectural weaknesses that limit the model's ability to generalize. Our contribution becomes threefold: We first target model discrepancies by adapting the Kool et al. method and its loss function for Sparse Dynamic Attention based on the alpha-entmax activation. We then target inherent differences through the use of a mixed instance training method that has been shown to outperform single instance training in certain scenarios. Finally, we introduce a framework for inference level data augmentation that improves performance by leveraging the model's lack of invariance to rotation and dilation changes.


title: On the distribution of Iwasawa invariants associated to multigraphs

abstract: Let $\ell$ be a prime number. The Iwasawa theory of multigraphs is the systematic study of growth patterns in the number of spanning trees in abelian $\ell$-towers of multigraphs. In this context, growth patterns are realized by certain analogues of Iwasawa invariants, which depend on the prime $\ell$ and the abelian $\ell$-tower of multigraphs. We formulate and study statistical questions about the behaviour of the Iwasawa $\mu$ and $\lambda$ invariants.


title: More on minors of Hermitian (quasi-)Laplacian matrix of the second kind for mixed graphs

abstract: A mixed graph $M_{G}$ is the graph obtained from an unoriented simple graph $G$ by giving directions to some edges of $G$, where $G$ is often called the underlying graph of $M_{G}$. In this paper, we introduce two classes of incidence matrices of the second kind of $M_{G}$, and discuss the determinants of these two matrices for rootless mixed trees and unicyclic mixed graphs. Applying these results, we characterize the explicit expressions of various minors for Hermitian (quasi-)Laplacian matrix of the second kind of $M_{G}$. Moreover, we give two sufficient conditions that the absolute values of all the cofactors of Hermitian (quasi-)Laplacian matrix of the second kind are equal to the number of spanning trees of the underlying graph $G$.


title: Warp and Weft Wiring method for rapid, modifiable, self-aligned, and bonding-free fabrication of multi electrodes microfluidic sensors

abstract: The need for rapid fabrication of microfluidic devices has become increasingly critical as microfluidics become part of biomedical sensors. Using Warp and Weft Wiring (WWW) of copper wires, this paper presents a novel low-cost method for rapid, self-aligned, bonding-free, and modifiable fabrication of multi-electrodes microfluidic sensors. All the proposed features are promising and highly recommended for the development of Point-of-Care Tests (POCTs), while most of the conventional methods have low chances of coming out of the research labs and play no role in POCTs development. To have an experimental proof of concept, the proposed chip was fabricated and then tested with two sets of experiments that showed the potential applications of water quality management, hygiene, biomedical impedance measurement, cell analysis, flow cytometry, etc.


title: On the use of graph neural networks and shape-function-based gradient computation in the deep energy method

abstract: A graph neural network (GCN) is employed in the deep energy method (DEM) model to solve the momentum balance equation in 3D for the deformation of linear elastic and hyperelastic materials due to its ability to handle irregular domains over the traditional DEM method based on a multilayer perceptron (MLP) network. Its accuracy and solution time are compared to the DEM model based on a MLP network. We demonstrate that the GCN-based model delivers similar accuracy while having a shorter run time through numerical examples. Two different spatial gradient computation techniques, one based on automatic differentiation (AD) and the other based on shape function (SF) gradients, are also accessed. We provide a simple example to demonstrate the strain localization instability associated with the AD-based gradient computation and show that the instability exists in more general cases by four numerical examples. The SF-based gradient computation is shown to be more robust and delivers an accurate solution even at severe deformations. Therefore, the combination of the GCN-based DEM model and SF-based gradient computation is potentially a promising candidate for solving problems involving severe material and geometric nonlinearities.


title: Modeling of radiative emission from shallow color centers in single crystalline diamond

abstract: Optically active defects in diamond are widely used as bright single-photon sources for quantum sensing, computing, and communication. For many applications, it is useful to place the emitter close to the diamond surface, where the radiative properties of the emitter are strongly modified by its dielectric environment. It is well-known that the radiative power from an electric dipole decreases as the emitter approaches an interface with a lower-index dielectric, leading to an increase in the radiative lifetime. For emitters in crystalline solids, modeling of this effect needs to take into account the crystal orientation and direction of the surface cut, which can greatly impact the emission characteristics. In this paper, we provide a framework for analyzing the emission rates of shallow (<100 nm) defects, in which optical transitions are derived from electric dipoles in a plane perpendicular to their spin axis. We present our calculations for the depth-dependent radiative lifetime for color centers in (100)-, (110)-, and (111)-cut diamond, which can be extended to other vacancy defects in diamond.


title: On the Selection of Tuning Parameters for Patch-Stitching Embedding Methods

abstract: While classical scaling, just like principal component analysis, is parameter-free, other methods for embedding multivariate data require the selection of one or several tuning parameters. This tuning can be difficult due to the unsupervised nature of the situation. We propose a simple, almost obvious, approach to supervise the choice of tuning parameter(s): minimize a notion of stress. We apply this approach to the selection of the patch size in a prototypical patch-stitching embedding method, both in the multidimensional scaling (aka network localization) setting and in the dimensionality reduction (aka manifold learning) setting. In our study, we uncover a new bias--variance tradeoff phenomenon.


title: Software-defined Dynamic 5G Network Slice Management for Industrial Internet of Things

abstract: This paper addresses the challenges of delivering fine-grained Quality of Service (QoS) and communication determinism over 5G wireless networks for real-time and autonomous needs of Industrial Internet of Things (IIoT) applications while effectively sharing network resources. Specifically, this work presents DANSM, a software-defined, dynamic and autonomous network slice management middleware for 5G-based IIoT use cases, such as adaptive robotic repair. Empirical studies evaluating DANSM on our testbed comprising a Free5GC-based core and UERANSIM-based simulations reveal that the software-defined DANSM solution can efficiently balance the traffic load in the data plane thereby reducing the end-to-end response time and improve the service performance by completing 34% more subtasks than a Modified Greedy Algorithm (MGA), 64% more subtasks than First Fit Descending (FFD) and 22% more subtasks than Best Fit Descending (BFD) approaches all while minimizing operational costs.