title: Does image resolution impact chest X-ray based fine-grained Tuberculosis-consistent lesion segmentation?

abstract: Deep learning (DL) models are state-of-the-art in segmenting anatomical and disease regions of interest (ROIs) in medical images. Particularly, a large number of DL-based techniques have been reported using chest X-rays (CXRs). However, these models are reportedly trained on reduced image resolutions for reasons related to the lack of computational resources. Literature is sparse in discussing the optimal image resolution to train these models for segmenting the Tuberculosis (TB)-consistent lesions in CXRs. In this study, we investigated the performance variations using an Inception-V3 UNet model using various image resolutions with/without lung ROI cropping and aspect ratio adjustments, and (ii) identified the optimal image resolution through extensive empirical evaluations to improve TB-consistent lesion segmentation performance. We used the Shenzhen CXR dataset for the study which includes 326 normal patients and 336 TB patients. We proposed a combinatorial approach consisting of storing model snapshots, optimizing segmentation threshold and test-time augmentation (TTA), and averaging the snapshot predictions, to further improve performance with the optimal resolution. Our experimental results demonstrate that higher image resolutions are not always necessary, however, identifying the optimal image resolution is critical to achieving superior performance.


title: Numerical investigation of progressive damage and associated seismicity on a laboratory fault

abstract: Understanding rock shear failure behavior is crucial to gain insights into slip-related geohazards such as rock avalanches, landslides, and earthquakes. However, descriptions of the progressive damage on the shear surface are still incomplete or ambiguous. In this study, we use the hybrid finite-discrete element method (FDEM) to simulate a shear experiment and obtain a detailed comprehension of shear induced progressive damage and the associated seismic activity. We built a laboratory fault model from high resolution surface scans and micro-CT imaging. Our results show that under quasi-static shear loading, the fault surface experiences local dynamic seismic activities. We found that the seismic activity is related to the stress concentration on interlocking asperities. This interlocking behavior (i) causes stress concentration at the region of contact that could reach the compressive strength, and (ii) produces tensile stress up to the tensile strength in the region adjacent to the contact area. Thus, different failure mechanisms and damage patterns including crushing and sub-vertical fracturing are observed on the rough surface. Asperity failure creates rapid local slips resulting in significant stress perturbations that alter the overall stress condition and may trigger the slip of adjacent critically stressed asperities. We found that the spatial distribution of the damaged asperities and the seismic activity is highly heterogeneous; regions with intense asperity interactions formed gouge material, while others exhibit minimal to no damage. These results emphasize the important role of surface roughness in controlling the overall shear behavior and the local dynamic seismic activities on faults.


title: Existence and Blow-up of solutions for Stochastic Modified Two-component Camassa-Holm System

abstract: In this paper, we consider the modified two-component Camassa-Holm System with multiplicative noise. For these SPDEs, we first establish the local existence and pathwise uniqueness of the pathwise solutions in Sobolev spaces $H^{s}\times H^{s}, s>\frac{3}{2}$. Then we show that strong enough noise can actually prevent blow-up with probability 1. Finally, we analyse the effects of weak noise and present conditions on the initial data that lead to the global existence and the blow-up in finite time of the solutions, and their associated probabilities are also obtained.


title: Constraints on the lunar core viscosity from tidal deformation

abstract: We use the tidal deformations of the Moon induced by the Earth and the Sun as a tool for studying the inner structure of our satellite. Based on measurements of the degree-two tidal Love numbers k2 and h2 and dissipation coefficients from the GRAIL mission, Lunar Laser Ranging and Laser Altimetry on board of the LRO spacecraft, we perform Monte Carlo samplings for 120,000 possible combinations of thicknesses and viscosities for two classes of the lunar models. The first one includes a uniform core, a low viscosity zone (LVZ) at the core-mantle boundary, a mantle and a crust. The second one has an additional inner core. All models are consistent with the lunar total mass as well as its moment of inertia. By comparing predicted and observed parameters for the tidal deformations we find that the existence of an inner core cannot be ruled out. Furthermore, by deducing temperature profiles for the LVZ and an Earth-like mantle, we obtain stringent constraints on the radius (500 +- 1) km, viscosity,21 (4.5 +- 0.8) x 10^16 Pa.s and the density (3400 +- 10) kg/m^3 of the LVZ. We also infer the first estimation for the outer core viscosity, (2.07 +- 1.03) x 10^17 Pa.s, for two different possible structures: a Moon with a 70 km thick outer core and a large inner core (290 km radius with a density of 6000 kg/m3), and a Moon with a thicker outer core (169 km thick) but a denser and smaller inner core (219 km radius for 8000 kg/m^3).


title: Deep Reinforcement Learning for Autonomous Ground Vehicle Exploration Without A-Priori Maps

abstract: Autonomous Ground Vehicles (AGVs) are essential tools for a wide range of applications stemming from their ability to operate in hazardous environments with minimal human operator input. Effective motion planning is paramount for successful operation of AGVs. Conventional motion planning algorithms are dependent on prior knowledge of environment characteristics and offer limited utility in information poor, dynamically altering environments such as areas where emergency hazards like fire and earthquake occur, and unexplored subterranean environments such as tunnels and lava tubes on Mars. We propose a Deep Reinforcement Learning (DRL) framework for intelligent AGV exploration without a-priori maps utilizing Actor-Critic DRL algorithms to learn policies in continuous and high-dimensional action spaces directly from raw sensor data. The DRL architecture comprises feedforward neural networks for the critic and actor representations in which the actor network strategizes linear and angular velocity control actions given current state inputs, that are evaluated by the critic network which learns and estimates Q-values to maximize an accumulated reward. Three off-policy DRL algorithms, DDPG, TD3 and SAC, are trained and compared in two environments of varying complexity, and further evaluated in a third with no prior training or knowledge of map characteristics. The agent is shown to learn optimal policies at the end of each training period to chart quick, collision-free exploration trajectories, and is extensible, capable of adapting to an unknown environment without changes to network architecture or hyperparameters. The best algorithm is further evaluated in a realistic 3D environment.


title: ROBUSfT: Robust Real-Time Shape-from-Template, a C++ Library

abstract: Tracking the 3D shape of a deforming object using only monocular 2D vision is a challenging problem. This is because one should (i) infer the 3D shape from a 2D image, which is a severely underconstrained problem, and (ii) implement the whole solution pipeline in real-time. The pipeline typically requires feature detection and matching, mismatch filtering, 3D shape inference and feature tracking algorithms. We propose ROBUSfT, a conventional pipeline based on a template containing the object's rest shape, texturemap and deformation law. ROBUSfT is ready-to-use, wide-baseline, capable of handling large deformations, fast up to 30 fps, free of training, and robust against partial occlusions and discontinuity in video frames. It outperforms the state-of-the-art methods in challenging datasets. ROBUSfT is implemented as a publicly available C++ library and we provide a tutorial on how to use it in https://github.com/mrshetab/ROBUSfT


title: Role of chemical potential at kinetic freeze-out using Tsallis non-extensive statistics in proton-proton collisions at the Large Hadron Collider

abstract: The charged-particle transverse momentum spectra ($p_{\rm T}$-spectra) measured by the ALICE collaboration for $pp$ collisions at $\sqrt {s} =$ 7 and 13 TeV have been studied using a thermodynamically consistent form of Tsallis non-extensive statistics. The Tsallis distribution function is fitted to the $p_{\rm T}$-spectra and the results are analyzed as a function of final state charged-particle multiplicity for various light flavor and strange particles, such as $\pi^{\pm}, K^{\pm}, p+\bar{p}, \phi, \Lambda+\bar{\Lambda}, \Xi+\bar{\Xi}, \Omega+\bar{\Omega}$. At the LHC energies, particles and antiparticles are produced in equal numbers. However, the equality of particle and antiparticle yields at the kinetic freeze-out may imply that they have the same but opposite chemical potential which is not necessarily zero. We use an alternative procedure that makes use of parameter redundancy, by introducing a finite chemical potential at the kinetic freeze-out stage. This article emphasizes the importance of the chemical potential of the system produced in $pp$ collisions at the LHC energies using the Tsallis distribution function which brings the system to a single freeze-out scenario.


title: A Theory of Theories

abstract: We take a tour through the past, present and future of Effective Field Theory, with applications ranging from LHC physics to cosmology.


title: Toeplitz operators and the full asymptotic torsion forms

abstract: This paper aims to study the asymptotic expansion of analytic torsion forms associated with a certain series of flat bundles. We prove the existence of the full expansion and give a formula for the sub-leading term, while Bismut-Ma-Zhang have studied the first-order expansion and expressed the leading term as the integral of a locally computable differential form.


title: Manifold Restricted Interventional Shapley Values

abstract: Shapley values are model-agnostic methods for explaining model predictions. Many commonly used methods of computing Shapley values, known as off-manifold methods, rely on model evaluations on out-of-distribution input samples. Consequently, explanations obtained are sensitive to model behaviour outside the data distribution, which may be irrelevant for all practical purposes. While on-manifold methods have been proposed which do not suffer from this problem, we show that such methods are overly dependent on the input data distribution, and therefore result in unintuitive and misleading explanations. To circumvent these problems, we propose ManifoldShap, which respects the model's domain of validity by restricting model evaluations to the data manifold. We show, theoretically and empirically, that ManifoldShap is robust to off-manifold perturbations of the model and leads to more accurate and intuitive explanations than existing state-of-the-art Shapley methods.


title: Modeling the vertical distribution of the Milky Way's flat subsystem objects

abstract: This paper is an initial stage of consideration of the general problem of joint modeling of the vertical structure of a Galactic flat subsystem and the average surface of the disk of the Galaxy, taking into account the natural and measurement dispersions. We approximate the average surface of the Galactic disk with a polynomial model and determine its parameters by minimizing the squared deviations of objects along the normal to the model surface. The developed method allows us to simultaneously identify significant details of the Galactic warping and estimate the offset $z_\odot$ of the Sun relative to the average (non-flat) surface of the Galactic disk and the vertical scale of the object system for an arbitrary area of the disk covered by data. The method is applied to data on classical Cepheids. Significant local extremes of the average disk surface model were found: the minimum in the first Galactic quadrant and the maximum in the second. A well-known warp in the third quadrant has been confirmed. The optimal order of the model was found to be $n_\text{o}=4$. The local (near the Sun, $n_\text{o}=0$) estimate of $z_\odot = 28.1 \pm \left.6.1\right|_{\text{stat.}}\left.{}\pm1.3\right|_{\text{cal.}}$ pc is close to the non-local ($n_\text{o}=4$) $z_\odot = 27.1 \pm \left.8.8\right|_{\text{stat.}}\left.{}^{+1.3}_{-1.2}\right|_{\text{cal.}}$ pc, which suggests that the proposed method eliminates the influence of warping on the $z_\odot$ estimate. However, the non-local estimate of the vertical standard deviation of Cepheids $\sigma_{\rho} = 132.0 \pm \left.3.7\right|_{\text{stat.}}\left.{}^{+6.3}_{-5.9}\right|_{\text{cal.}}$ pc differs significantly from the local $\sigma_{\rho} = \left.76.5 \pm 4.4\right|_{\text{stat.}}\left.{}^{+3.6}_{-3.4}\right|_{\text{cal.}}$ pc, which means the need to introduce more complex models outside the Sun's vicinity.


title: Integrated Analysis of Coarse-Grained Guidance for Traffic Flow Stability

abstract: Autonomous vehicles (AVs) enable more efficient and sustainable transportation systems. Ample studies have shown that controlling a small fraction of AVs can smooth traffic flow and mitigate traffic congestion. However, deploying AVs in real-world systems poses challenges due to safety and cost concerns. A viable alternative approach that can be implemented in the near future is $\textit{coarse-grained guidance}$, where human drivers are guided by real-time instructions, updated every $\Delta$ seconds, to stabilize the traffic. While previous theoretical studies consider stability analysis for continuous AV control, this article presents the first integrated theoretical analysis that directly relates the guidance provided to the human drivers to the traffic flow stability outcome. Casting the problem into the Lyapunov stability framework, this study derives sufficient conditions for coarse-grained guidance with hold length $\Delta$ to stabilize the system, and provides extensions of the analysis to incorporate additional human driving behaviors such as magnitude error and reaction delay. Numerical simulations reveal that the theoretical analysis closely matches simulated results. The analysis further offers insights into the relationship between system parameter and stability criteria, and can be leveraged to design improved controllers with greater maximum hold length.


title: Schatten-von Neumann properties for H\"ormander classes on compact Lie groups

abstract: Let $G$ be a compact Lie group of dimension $n.$ In this work we characterise the membership of classical pseudo-differential operators on $G$ in the trace class ideal $S_{1}(L^2(G)),$ as well as in the setting of the Schatten ideals $S_{r}(L^2(G)),$ for all $r>0.$ In particular, we deduce Schatten characterisations of elliptic pseudo-differential operators of $(\rho,\delta)$-type for the large range $0\leq \delta<\rho\leq 1.$ Additional necessary and sufficient conditions are given in terms of the matrix-valued symbols of the operators, which are global functions on the phase space $G\times \widehat{G},$ with the momentum variables belonging to the unitary dual $\widehat{G}$ of $G$. In terms of the parameters $(\rho,\delta),$ on the torus $\mathbb{T}^n,$ we demonstrate the sharpness of our results showing the existence of atypical operators in the exotic class $\Psi^{-\varkappa}_{0,0}(\mathbb{T}^n),$ $\varkappa>0,$ belonging to all the Schatten ideals. Additional order criteria are given in the setting of classical pseudo-differential operators. We present also some open problems in this setting.


title: Experimental characterization of an ultra-broadband dual-mode symmetric Y-junction based on metamaterial waveguides

abstract: Silicon photonic integrated circuits routinely require 3-dB optical power dividers with minimal losses, small footprints, ultra-wide bandwidths, and relaxed manufacturing tolerances to distribute light across the chip and as a key building block to form more complex devices. Symmetric Y-junctions stand out among other power splitting devices owing to their wavelength-independent response and a straightforward design. Yet, the limited resolution of current fabrication methods results in a minimum feature size (MFS) at the tip between the two Y-junction arms that leads to significant losses for the fundamental mode. Here we propose to circumvent this limitation by leveraging subwavelength metamaterials in a new type of ultra-broadband and fabrication-tolerant Y-junction. An exhaustive experimental study over a 260 nm bandwidth (1420-1680 nm) shows excess loss below 0.3 dB for the fundamental transverse-electric mode (TE0) for a high-resolution lithographic process (MFS about 50 nm) and less than 0.5 dB for a fabrication resolution of 100 nm. Subwavelength Y-junctions with deterministically induced errors of plus-minus 10 nm further demonstrated robust fabrication tolerances. Moreover, the splitter exhibits excess loss lower than 1 dB for the first-order transverse-electric mode (TE1) within a 100 nm bandwidth (1475-1575 nm), using high-resolution lithography.


title: Numerical study of conforming space-time methods for Maxwell's equations

abstract: Time-dependent Maxwell's equations govern electromagnetics. Under certain conditions, we can rewrite these equations into a partial differential equation of second order, which in this case is the vectorial wave equation. For the vectorial wave, we investigate the numerical application and the challenges in the implementation. For this purpose, we consider a space-time variational setting, i.e. time is just another spatial dimension. More specifically, we apply integration by parts in time as well as in space, leading to a space-time variational formulation with different trial and test spaces. Conforming discretizations of tensor-product type result in a Galerkin--Petrov finite element method that requires a CFL condition for stability. For this Galerkin--Petrov variational formulation, we study the CFL condition and its sharpness. To overcome the CFL condition, we use a Hilbert-type transformation that leads to a variational formulation with equal trial and test spaces. Conforming space-time discretizations result in a new Galerkin--Bubnov finite element method that is unconditionally stable. In numerical examples, we demonstrate the effectiveness of this Galerkin--Bubnov finite element method. Furthermore, we investigate different projections of the right-hand side and their influence on the convergence rates. This paper is the first step towards a more stable computation and a better understanding of vectorial wave equations in a conforming space-time approach.


title: Local Convergence Behaviour of Generalized Gauss-Newton Multiple Shooting, Single Shooting and Differential Dynamic Programming

abstract: We revisit three classical numerical methods for solving unconstrained optimal control problems - multiple shooting, single shooting, and differential dynamic programming - and examine their local convergence behaviour. In particular, we show that all three methods converge with the same linear rate if a Gauss-Newton (GN), or more general a Generalized Gauss-Newton (GGN), Hessian approximation is used, which is the case in widely used implementations such as iLQR.


title: A Sufficient Condition for the Super-linearization of Polynomial Systems

abstract: We provide in this paper a sufficient condition for a polynomial dynamical system $\dot x(t) = f(x(t))$ to be super-linearizable, i.e., to be such that all its trajectories are linear projections of the trajectories of a linear dynamical system. The condition is expressed in terms of the hereby introduced weighted dependency graph $G$, whose nodes $v_i$ correspond to variables $x_i$ and edges $v_iv_j$ have weights $\frac{\partial f_j}{\partial x_i}$. We show that if the product of the edge weights along any cycle in $G$ is a constant, then the system is super-linearizable. The proof is constructive, and we provide an algorithm to obtain super-linearizations and illustrate it on an example.


title: Imbalanced Classification In Faulty Turbine Data: New Proximal Policy Optimization

abstract: There is growing importance to detecting faults and implementing the best methods in industrial and real-world systems. We are searching for the most trustworthy and practical data-based fault detection methods proposed by artificial intelligence applications. In this paper, we propose a framework for fault detection based on reinforcement learning and a policy known as proximal policy optimization. As a result of the lack of fault data, one of the significant problems with the traditional policy is its weakness in detecting fault classes, which was addressed by changing the cost function. Using modified Proximal Policy Optimization, we can increase performance, overcome data imbalance, and better predict future faults. When our modified policy is implemented, all evaluation metrics will increase by $3\%$ to $4\%$ as compared to the traditional policy in the first benchmark, between $20\%$ and $55\%$ in the second benchmark, and between $6\%$ and $14\%$ in the third benchmark, as well as an improvement in performance and prediction speed compared to previous methods.


title: Design, Modeling and Control of a Quadruped Robot SPIDAR: Spherically Vectorable and Distributed Rotors Assisted Air-Ground Amphibious Quadruped Robot

abstract: Multimodal locomotion capability is an emerging topic in robotics field, and various novel mobile robots have been developed to enable the maneuvering in both terrestrial and aerial domains. Among these hybrid robots, several state-of-the-art bipedal \robots enable the complex walking motion which is interlaced with flying. These robots are also desired to have the manipulation ability; however, it is difficult for the current forms to keep stability with the joint motion in midair due to the central\ized rotor arrangement. Therefore, in this work, we develop a novel air-ground amphibious quadruped robot called SPIDAR which is assisted by spherically vectorable rotors distributed in each link to enable both walking motion and transformable flight. F\irst, we present a unique mechanical design for quadruped robot that enables terrestrial and aerial locomotion. We then reveal the modeling method for this hybrid robot platform, and further develop an integrated control strategy for both walking and fl\ying with joint motion. Finally, we demonstrate the feasibility of the proposed hybrid quadruped robot by performing a seamless motion that involves static walking and subsequent flight. To the best of our knowledge, this work is the first to achieve a \quadruped robot with multimodal locomotion capability, which also shows the potential of manipulation in multiple domains.


title: Protein Co-Enrichment Analysis of Extracellular Vesicles

abstract: Extracellular Vesicles (EVs) carry cell-derived proteins that confer functionality and selective cell uptake. However, whether proteins are packaged stochastically or co-enriched within individual EVs, and whether co-enrichment fluctuates under homeostasis and disease, has not been measured. EV abundance and protein global relative expression have been qualified by bulk analysis. Meanwhile, co-enrichment is not directly accessible via bulk measurement and has not been reported for single EV analysis. Here, we introduce the normalized index of co-enrichment (NICE) to measure protein co-enrichment. NICE was derived by (i) capturing EVs based on the expression of a membrane-bound protein, (ii) probing for the co-expression of a second protein at the population level - EV integrity underwrites the detection of single EV co-expression without the need to resolve single EVs - and (iii) normalizing measured values using two universal normalization probes. Axiomatically, NICE = 1 for stochastic inclusion or no overall co-enrichment, while for positive and negative co-enrichment NICE > 1 or < 1, respectively. We quantified the NICE of tetraspanins, growth factor receptors and integrins in EVs of eight breast cancer cell lines of varying metastatic potential and organotropism, combinatorially mapping up to 104 protein pairs. Our analysis revealed protein enrichment and co-expression patterns consistent with previous findings. For the organotropic cell lines, most protein pairs were co-enriched on EVs, with the majority of NICE values between 0.2 to 11.5, and extending from 0.037 to 80.4. Median NICE were either negative, neutral or positive depending on the cells. NICE analysis is easily multiplexed and is compatible with microarrays, bead-based and single EV assays. Additional studies are needed to deepen our understanding of the potential and significance of NICE for research and clinical uses.


title: Optimal social security timing

abstract: The optimal age that a retiree claims social security retirement benefits is in general a complicated function of many factors. However, if the beneficiary's finances and health are not the constraining factors, it is possible to formally derive mathematical models that maximize a well-defined measure of his total benefits. A model that takes into account various factors such as the increase in the benefits for delayed claims and the penalties for early retirement, the advantages of investing some of the benefits in the financial markets, and the effects of cost-of-living adjustments shows that not waiting until age 70 is almost always the better option. The optimal claiming age that maximizes the total benefits, however, depends on the expected market returns and the rate of cost-of-living adjustments, with the higher market rates in general pushing the optimal age lower. The models presented here can be easily tailored to address the particular circumstances and goals of any individual.


title: Coffee stain effect on a fibre from axisymmetric droplets

abstract: The so-called coffee stain effect has been intensively studied over the past decades, but most of the studies are focused on sessile droplets. In this paper, we analyse the origin of the difference between the deposition of suspended particles in a sessile drop and in an axisymmetric drop deposited on a fibre. First, we model the shape of a drop on a fibre and its evaporative flux with some approximations to derive analytical calculations. Then, for pinned contact lines, we solve the hydrodynamics equations in the liquid phase under the lubrication approximation to determine the flow velocity toward the contact lines. We comment these results by comparison to a sessile drop of similar evaporating conditions, and we show that the substrate curvature plays a role on the contact line depinning, the local evaporative flux, and the liquid flow field. The competition between the advection and the Brownian motion indicates that the transport of the particles toward the contact line occurs in a volume localised in the close vicinity of the contact lines for a drop on a fibre. Thus, the fibre geometry induces a weaker accumulation of particles at the contact line compared to a sessile drop, leading to the more homogeneous deposit observed experimentally.


title: Surface-Sensitive Raman Scattering by Transferable Nanoporous Plasmonic Membranes

abstract: Raman spectroscopy is a powerful technique to characterize materials. It probes non-destructively chemical composition, crystallinity, defects, strain and coupling phenomena. However, the Raman response of surfaces or thin films is often weak and obscured by dominant bulk signals. Here we overcome this limitation by placing a transferable porous gold membrane (PAuM) on top of the surface of interest. Slot-like nanopores in the membrane act as plasmonic slot antennas and enhance the Raman response of the surface or thin film underneath. Simultaneously, the PAuM suppresses the penetration of the excitation laser into the bulk, efficiently blocking the bulk Raman signal. Using graphene as a model surface, we show that these two simultaneous effects lead to an increase in the surface-to-bulk Raman signal ratio by three orders of magnitude. We find that 90% of the Raman enhancement occurs within the top 2.5 of the material, demonstrating truly surface-sensitive Raman scattering. To validate our approach, we analyze the surface of a LaNiO3 thin film. We observe a Raman mode splitting for the LaNiO3 surface-layer, which is spectroscopic evidence that the surface structure differs from the bulk. This result underpins that PAuM give direct access to Raman signatures of surfaces and their structural properties.


title: PSR~J1910$-$5959A: A rare gravitational laboratory for testing white dwarf models

abstract: PSRJ1910-5959A (J1910A) is a binary millisecond pulsar in a 0.837 day circular orbit around a helium white dwarf (HeWD) companion. This pulsar is located 6.3 arcmin away from the centre of the globular cluster NGC6752. Given the large offset, the association of the pulsar to NGC6752 has been debated. We have made use of two decades of archival Parkes 64-m "Murriyang" telescope data and recently carried out observations with the MeerKAT telescope. We obtained Pulse times of arrival using standard data reduction techniques and analysed using Bayesian pulsar timing techniques. We analysed the pulsar's total intensity and polarisation profile, to study the interstellar scattering along the line of sight, and the pulsar's geometry by applying the rotating vector model. We obtain precise measurements of several post-Keplerian parameters: the range $r=0.202(6)T_\odot$ and shape s=0.999823(4) of the Shapiro delay, from which we infer the orbital inclination to be $88.9^{+0.15}_{-0.14}\deg$ and the masses of both the pulsar and the companion to be $1.55(7)M_{\odot}$ and $0.202(6)M_{\odot}$ respectively; a secular change in the orbital period $\dot{P}_{\rm b}=-53^{+7.4}_{-6.0}\times 10^{-15}$\,s\,s$^{-1}$ that proves the association to NGGC6752 and a secular change in the projected semi-major axis of the pulsar $\dot{x}= -40.7^{+7.3}_{-8.2}\times10^{-16}$\,s\,s$^{-1}$ that is likely caused by the spin-orbit interaction from a misaligned HeWD spin, at odds with the likely isolated binary evolution of the system. We also discuss some theoretical models for the structure and evolution of WDs in NS-WD binaries by using J1910A's companion as a test bed. J1910A is a rare system for which several parameters of both the pulsar and the HeWD companion can be accurately measured. As such, it is a test bed to discriminate between alternative models for HeWD structure and cooling.


title: Reverberation of pulsar wind nebulae (II): Anatomy of the "thin-shell'' evolution

abstract: During its early evolution, a pulsar wind nebula (PWN) sweeps the inner part of the supernova ejecta and forms a thin massive shell. Later on, when the shell has been reached by the reverse shock of the supernova remnant, the evolution becomes more complex, in most cases reverting the expansion into a compression: this later phase is called "reverberation". Computations done so far to understand this phase have been mostly performed in the thin-shell approximation, where the evolution of the PWN radius is assimilated to that of the swept-up shell under the effect of both the inner pressure from the PWN, and the outer pressure from the supernova remnant. Despite the thin-shell approach seems rather justifiable, its implementations have so far been inaccurate, and its correctness, never tested. The outer pressure was naively assumed to be scaled according to the Sedov solution (or a constant fraction of it) along the entire evolution. The thin-shell assumption itself fails along the process, being the shell no longer thin in comparison with the size of the PWN. Here, through a combination of numerical models, dimensional arguments, and analytic approximations, we present a detailed analysis of the interaction of the PWN with the supernova remnant. We provide a new analytic approximation of the outer pressure, beyond the Sedov solution, and a revised "thin-shell" able to reproduce results from numerical simulations. Finally, we compute the efficiency by which the PWN is compressed during reverberation over a wide population of sources.


title: Study of neutrino mass matrices with vanishing trace and one vanishing minor

abstract: In this work we carry out a systematic texture study of the neutrino mass matrix with the ansatzes - (i) one vanishing minor and (ii) the zero sum of the mass eigenvalues with the CP phases (henceforth vanishing trace). There are six possible textures of a neutrino mass matrix with one vanishing minor. The viability of each texture is checked with $3\sigma$ values of current neutrino data by drawing scatter plots. In our analysis we are motivated to use the ratio of solar to atmospheric mass-squared differences $R_{\nu}$ for its precise measurement (and also the atmospheric mixing angle $\theta_{23}$) to constrain phenomenologically first the Dirac CP phase $\delta$ in the range of $0^{\degree}-360^{\degree}$ for a given texture with the solutions of the constraint equations. Subsequently we employ this constrained $\delta$ to determine the range of completely unknown Majorana CP Phases ($\alpha$ and $\beta$) for all the viable textures. We also check the neutrinoless double beta decay rate, $|m_{ee}|$ and the Jarlskog invariant, $J_{cp}$ for the textures. Finally the symmetry realization of all the viable textures under the flavor symmetry group $Z_5$ via seesaw mechanism is implemented along with the FN mechanism to determine mass hierarchy structure.


title: Rethinking Voxelization and Classification for 3D Object Detection

abstract: The main challenge in 3D object detection from LiDAR point clouds is achieving real-time performance without affecting the reliability of the network. In other words, the detecting network must be confident enough about its predictions. In this paper, we present a solution to improve network inference speed and precision at the same time by implementing a fast dynamic voxelizer that works on fast pillar-based models in the same way a voxelizer works on slow voxel-based models. In addition, we propose a lightweight detection sub-head model for classifying predicted objects and filter out false detected objects that significantly improves model precision in a negligible time and computing cost. The developed code is publicly available at: https://github.com/YoushaaMurhij/RVCDet.


title: Determination of the $W$-boson mass at hadron colliders

abstract: We introduce an observable relevant for the determination of the $W$-boson mass $m_W$ at hadron colliders. This observable is defined as an asymmetry around the jacobian peak of the charged-lepton transverse-momentum distribution in the charged-current Drell-Yan process. We discuss the observable's theoretical prediction, presenting results at different orders in QCD, and showing its perturbative stability. Its definition as a single scalar number and its linear sensitivity to $m_W$ allow a clean extraction of the latter and a straightforward discussion of the associated theoretical systematics: a perturbative QCD uncertainty of ${\cal O}(\pm 5)$ MeV on $m_W$ can be established by means of this observable, relying solely on charged-current Drell-Yan information. Owing to its relatively inclusive nature, the observable displays desirable properties also from the experimental viewpoint, especially for the unfolding of detector effects. We show that a measurement of this observable can lead to a competitive experimental error on $m_W$ at the LHC.


title: A Formal Disproof of the Hirsch Conjecture

abstract: The purpose of this paper is the formal verification of a counterexample of Santos et al. to the so-called Hirsch Conjecture on the diameter of polytopes (bounded convex polyhedra). In contrast with the pen-and-paper proof, our approach is entirely computational: we implement in Coq and prove correct an algorithm that explicitly computes, within the proof assistant, vertex-edge graphs of polytopes as well as their diameter. The originality of this certificate-based algorithm is to achieve a tradeoff between simplicity and efficiency. Simplicity is crucial in obtaining the proof of correctness of the algorithm. This proof splits into the correctness of an abstract algorithm stated over proof-oriented data types and the correspondence with a low-level implementation over computation-oriented data types. A special effort has been made to reduce the algorithm to a small sequence of elementary operations (e.g., matrix multiplications, basic routines on sets and graphs), in order to make the derivation of the correctness of the low-level implementation more transparent. Efficiency allows us to scale up to polytopes with a challenging combinatorics. For instance, we formally check the two counterexamples of Matschke, Santos and Weibel to the Hirsch conjecture, respectively 20- and 23-dimensional polytopes with 36 425 and 73 224 vertices involving rational coefficients with up to 40 digits in their numerator and denominator. We also illustrate the performance of the method by computing the list of vertices or the diameter of well-known classes of polytopes, such as (polars of) cyclic polytopes involved in McMullen's Upper Bound Theorem.


title: Measuring and simulating the biophysical basis of the acoustic contrast factor of biological cells

abstract: The acoustic contrast factor (ACF) is calculated from the relative density and compressibility differences between a fluid and an object in the fluid. To name but one application, this acoustic contrast can be exploited using acoustophoretic systems to isolate cancer cells from a liquid biopsy, such as a blood sample. Knowing the ACF of a cancer cell represents a crucial step in the design of acoustophoretic systems for this purpose, potentially allowing the isolation of circulating cancer cells without labels or contact. For biological cells the static compressibility is different from the high frequency counterpart relevant for the ACF. In this study, we started by characterizing the ACF of low vs. high metastatic cell lines with known associated differences in phenotypic static E-modulus. The change in the static E-modulus, however, was not reflected in a change of the ACF, prompting a more in depth analysis of the influences on the ACF. We demonstrate that static E-modulus increased biological cells through formaldehyde fixation have an increased ACF. Conversely static E-modulus decreased biological cells treated with actin polymerization inhibitor cytochalasin D have a decreased ACF. Complementing these mechanical tests, a numerical COMSOL model was implemented and used to parametrically explore the effects of cell density, cell density ratios, dynamic compressibility and therefore the dynamic bulk modulus. Collectively the combined laboratory and numerical experiments reveal that a change in the static E-modulus alone might, but does not automatically lead to a change of the dynamic ACF for biological cells. This highlights the need for a multiparametic view of the biophysical basis of the cellular ACF, as well as the challenges in harnessing acoustophoretic systems to isolate circulating cells based on their mechanical properties alone.


title: A deep radius valley revealed by Kepler short cadence observations

abstract: The characteristics of the radius valley, i.e., an observed lack of planets between 1.5-2 Earth radii at periods shorter than about 100 days, provide insights into the formation and evolution of close-in planets. We present a novel view of the radius valley by refitting the transits of 431 planets using Kepler 1-minute short cadence observations, the vast majority of which have not been previously analysed in this way. In some cases, the updated planetary parameters differ significantly from previous studies, resulting in a deeper radius valley than previously observed. This suggests that planets are likely to have a more homogeneous core composition at formation. Furthermore, using support-vector machines, we find that the radius valley location strongly depends on orbital period and stellar mass and weakly depends on stellar age, with $\partial \log {\left(R_{p, \text{valley}} \right)}/ \partial \log{P} = -0.096_{-0.027}^{+0.023}$, $\partial \log {\left(R_{p, \text{valley}} \right)}/ \partial \log{M_{\star}} = 0.231_{-0.064}^{+0.053}$, and $\partial \log {\left(R_{p, \text{valley}} \right)}/ \partial \log{\left( \text{age} \right)} = 0.033_{-0.025}^{+0.017}$. These findings favour thermally-driven mass loss models such as photoevaporation and core-powered mass loss, with a slight preference for the latter scenario. Finally, this work highlights the value of transit observations with short photometric cadence to precisely determine planet radii, and we provide an updated list of precisely and homogeneously determined parameters for the planets in our sample.


title: On the number of Diophantine $m$-tuples in finite fields

abstract: We use a new argument to improve the error term in the asymptotic formula for the number of Diophantine $m$-tuples in finite fields, which is due to A. Dujella and M.Kazalicki (2021) and N. Mani and S. Rubinstein-Salzedo (2021).


title: The cyclotomic Grothendieck-Teichm\"{u}ller group and the motivic Galois group

abstract: We show that the level 2 case of the cyclotomic Grothendieck-Teichm\"{u}ller groups introduced by Enriquez coincides with the motivic Galois group of mixed Tate motives over $\mathbb{Z}[1/2]$.


title: Lower-temperature fabrication of airbridges by grayscale lithography to increase yield of nanowire transmons in circuit QED quantum processors

abstract: Quantum hardware based on circuit quantum electrodynamics makes extensive use of airbridges to suppress unwanted modes of wave propagation in coplanar-waveguide transmission lines. Airbridges also provide an interconnect enabling transmission lines to cross. Traditional airbridge fabrication produces a curved profile by reflowing resist at elevated temperature prior to metallization. The elevated temperature can affect the coupling energy and even yield of pre-fabricated Josephson elements of superconducting qubits, tuneable couplers and resonators. We employ grayscale lithography in place of reflow to reduce the peak airbridge processing temperature from $200$ to $150^\circ\mathrm{C}$, showing a substantial yield increase of transmon qubits with Josephson elements realized using Al-contacted InAs nanowires.


title: Multi-species generalization of the totally asymmetric simple exclusion process

abstract: Exclusion processes in one dimension first appeared in the 70s and have since dragged much attention from communities in different domains: stochastic processes, out-of-equilibriums statistical physics, and more recently integrable systems. While the state of the art for a single species totally asymmetric simple exclusion process (TASEP) can be described, from different aspects as mature, much less is known when multiple interacting species are present. Using tools from integrable systems and hydrodynamics in the first place and stochastic processes in the second place, this work attempts to study the behavior of a novel version of the model with different species of particles having hierarchical dynamics that depend on arbitrary parameters. While Burger's equation famously represents the hydrodynamic limit of TASEP with a single species, we present a counterpart coupled system of PDE representing the hydrodynamic limit for a model with two species. The solutions of these PDEs display a rich phenomenology of solutions best characterized through the underlying normal modes. This system with two species can be used as a toy model for studying driven diffusive systems with open boundaries. Using heuristics, we present results suggesting a general principle governing the boundary-induced phase diagram of systems with multiple coupled driven conserved quantities, generalizing thus the extremal current principle known for the case of a single driven quantity. Using integrability tools, we develop a formalism allowing the computation of the finite-time probability distribution of particle positions on the 1D lattice, generalizing therefore known results for TASEP and other multi-species models. We finally study the behavior and the impact of a single second-class impurity initially located at the interface separating two regions of different densities of first-class particles.


title: The Thousand-Pulsar-Array programme on MeerKAT -- VIII. The subpulse modulation of 1198 pulsars

abstract: We report on the subpulse modulation properties of 1198 pulsars using the Thousand-Pulsar-Array programme on MeerKAT. About 35% of the analysed pulsars exhibit drifting subpulses which are more pronounced towards the deathline, consistent with previous studies. We estimate that this common phenomenon is detectable in 60% of the overall pulsar population if high quality data were available for all. This large study reveals the evolution of drifting subpulses across the pulsar population in unprecedented detail. In particular, we find that the modulation period $P_3$ follows a V-shaped evolution with respect to the characteristic age $\tau_c$, such that the smallest $P_3$ values, corresponding to the Nyquist period $P_3>\sim2$, are found at $\tau_c>\sim10^{7.5}$ yr. The V-shaped evolution can be interpreted and reproduced if young pulsars possess aliased fast intrinsic $P_3$, which monotonically increase, ultimately achieving a slow unaliased $P_3$. Enhancement of irregularities in intrinsic subpulse modulation by aliasing in small $\tau_c$ pulsars would explain their observed less well defined $P_3$'s and weaker spectral features. Modelling these results as rotating subbeams, their circulation must slow down as the pulsar evolves. This is the opposite to that expected if circulation is driven by ExB drift. This can be resolved if the observed $P_3$ periodicity is due to a beat between an ExB system and the pulsar period. As a by-product, we identify the correct periods and spin-down rates for 12 pulsars, for which harmonically related values were reported in the literature.


title: Blobbed topological recursion from extended loop equations

abstract: We consider the $N\times N$ Hermitian matrix model with measure $d\mu_{E,\lambda}(M)=\frac{1}{Z} \exp(-\frac{\lambda N}{4} \mathrm{tr}(M^4)) d\mu_{E,0}(M)$, where $d\mu_{E,0}$ is the Gaussian measure with covariance $\langle M_{kl}M_{mn}\rangle=\frac{\delta_{kn}\delta_{lm}}{N(E_k+E_l)}$ for given $E_1,...,E_N>0$. It was previously understood that this setting gives rise to two ramified coverings $x,y$ of the Riemann sphere strongly tied by $y(z)=-x(-z)$ and a family $\omega^{(g)}_{n}$ of meromorphic differentials conjectured to obey blobbed topological recursion due to Borot and Shadrin. We develop a new approach to this problem via a system of six meromorphic functions which satisfy extended loop equations. Two of these functions are symmetric in the preimages of $x$ and can be determined from their consistency relations. An expansion at $\infty$ gives global linear and quadratic loop equations for the $\omega^{(g)}_{n}$. These global equations provide the $\omega^{(g)}_{n}$ not only in the vicinity of the ramification points of $x$ but also in the vicinity of all other poles located at opposite diagonals $z_i+z_j=0$ and at $z_i=0$. We deduce a recursion kernel representation valid at least for $g\leq 1$.


title: The chaotic emergence of thermalization in highly excited string decays

abstract: We analyse the most general process of a generic highly excited string that decays into a less excited, yet generic, highly excited string emitting a tachyon. We provide a simple and compact analytic description of the decay process which discriminates between and within the structure of every single microstate of the initial and final highly excited string. Taking into account the random nature of the decay process we extract the energy spectrum of highly excited strings, microstate by microstate, finding a behavior which corresponds to the greybody emission spectrum. In addition, by exploiting the analytic control of the decay process, we identify the origin of thermal effects which are triggered by the chaotic nature of the highly excited string interactions modeled by the microstates structure.


title: A simple theory for quantum quenches in the ANNNI model

abstract: In a recent numerical study by Haldar et al. (Phys. Rev. X 11, 031062) it was shown that signatures of proximate quantum critical points can be observed at early and intermediate times after certain quantum quenches. Said work focused mainly on the case of the axial next-nearest neighbour Ising (ANNNI) model. Here we construct a simple time-dependent mean-field theory that allows us to obtain a quantitatively accurate description of these quenches at short times, which for reasons we explain remains a fair approximation at late times (with some caveats). Our approach provides a simple framework for understanding the reported numerical results as well as fundamental limitations on detecting quantum critical points through quench dynamics. We moreover explain the origin of the peculiar oscillatory behaviour seen in various observables as arising from the formation of a long-lived bound state.


title: Truncation of contact defects in reaction-diffusion systems

abstract: Contact defects are time-periodic patterns in one space dimension that resemble spatially homogeneous oscillations with an embedded defect in their core region. For theoretical and numerical purposes, it is important to understand whether these defects persist when the domain is truncated to large spatial intervals, supplemented by appropriate boundary conditions. The present work shows that truncated contact defects exist and are unique on sufficiently large spatial intervals.


title: The spectral reconstruction of inclusive rates

abstract: A recently re-discovered variant of the Backus-Gilbert algorithm for spectral reconstruction enables the controlled determination of smeared spectral densities from lattice field theory correlation functions. A particular advantage of this approach is the \emph{a priori} specification of the kernel with which the underlying spectral density is smeared, allowing for variation of its peak position, smearing width, and functional form. If the unsmeared spectral density is sufficiently smooth in the neighborhood of a particular energy, it can be obtained from an extrapolation to zero smearing-kernel width at fixed peak position. A natural application for this approach is scattering processes summed over all hadronic final states. As a proof-of-principle test, an inclusive rate is computed in the two-dimensional O(3) sigma model from a two-point correlation function of conserved currents. The results at finite and zero smearing radius are in good agreement with the known analytic form up to energies at which 40-particle states contribute, and are sensitive to the 4-particle contribution to the inclusive rate. The straight-forward adaptation to compute the $R$-ratio in lattice QCD from two-point functions of the electromagnetic current is briefly discussed.


title: Diffusion Properties of a Brownian Ratchet with Coulomb Friction

abstract: The motion of a Brownian particle in the presence of Coulomb friction and an asymmetric spatial potential was evaluated in this study. The system exhibits a ratchet effect, i.e., an average directed motion even in the absence of an external force, induced by the coupling of non-equilibrium conditions with the spatial asymmetry. Both the average motion and the fluctuations of the Brownian particle were analysed. The stationary velocity shows a non-monotonic behaviour as a function of both the temperature and the viscosity of the bath. The diffusion properties of the particle, which show several time regimes, were also investigated. To highlight the role of non-linear friction in the dynamics, a comparison is presented with a linear model of a Brownian particle driven by a constant external force, which allows for analytical treatment. In particular, the study unveils that the passage times between different temporal regimes are strongly affected by the presence of Coulomb friction.


title: Cavity-Catalyzed Hydrogen Transfer Dynamics in an Entangled Molecular Ensemble under Vibrational Strong Coupling

abstract: Microcavities have been shown to influence the reactivity of molecular ensembles by strong coupling of molecular vibrations to quantized cavity modes. In quantum mechanical treatments of such scenarios, frequently idealized models with single molecules and scaled, effective molecule-cavity interactions or alternatively ensemble models with simplified model Hamiltonians are used. In this work, we go beyond these models by applying an ensemble variant of the Pauli-Fierz Hamiltonian for vibro-polaritonic chemistry and numerically solve the underlying time-dependent Schr\"odinger equation to study the cavity-induced quantum dynamics in an ensemble of thioacetylacetone (TAA) molecules undergoing hydrogen transfer under vibrational strong coupling (VSC) conditions. Beginning with a single molecule coupled to a single cavity mode, we show that the cavity indeed enforces hydrogen transfer from an enol to an enethiol configuration with transfer rates significantly increasing with light-matter interaction strength. This positive effect of the cavity on reaction rates is different from several other systems studied so far, where a retarding effect of the cavity on rates was found. It is argued that the cavity ``catalyzes'' the reaction by transfer of virtual photons to the molecule. The same concept applies to ensembles with up to $N=20$ TAA molecules coupled to a single cavity mode, where an additional, significant, ensemble-induced collective isomerization rate enhancement is found. The latter is traced back to complex entanglement dynamics of the ensemble, which we quantify by means of von Neumann-entropies. A non-trivial dependence of the dynamics on ensemble size is found, clearly beyond scaled single-molecule models, which we interpret as transition from a multi-mode Rabi to a system-bath-type regime as $N$ increases.


title: Benchmarking Robustness in Neural Radiance Fields

abstract: Neural Radiance Field (NeRF) has demonstrated excellent quality in novel view synthesis, thanks to its ability to model 3D object geometries in a concise formulation. However, current approaches to NeRF-based models rely on clean images with accurate camera calibration, which can be difficult to obtain in the real world, where data is often subject to corruption and distortion. In this work, we provide the first comprehensive analysis of the robustness of NeRF-based novel view synthesis algorithms in the presence of different types of corruptions. We find that NeRF-based models are significantly degraded in the presence of corruption, and are more sensitive to a different set of corruptions than image recognition models. Furthermore, we analyze the robustness of the feature encoder in generalizable methods, which synthesize images using neural features extracted via convolutional neural networks or transformers, and find that it only contributes marginally to robustness. Finally, we reveal that standard data augmentation techniques, which can significantly improve the robustness of recognition models, do not help the robustness of NeRF-based models. We hope that our findings will attract more researchers to study the robustness of NeRF-based approaches and help to improve their performance in the real world.


title: Electronic Born-Oppenheimer Approximation in Nuclear-Electronic Orbital Dynamics

abstract: Within the nuclear-electronic orbital (NEO) framework, the real-time NEO time-dependent density functional theory (RT-NEO-TDDFT) approach enables the simulation of coupled electronic-nuclear dynamics. In this approach, the electrons and quantum nuclei are propagated in time on the same footing. A relatively small time step is required to propagate the much faster electronic dynamics, thereby prohibiting the simulation of long-time nuclear quantum dynamics. Herein, the electronic Born--Oppenheimer (BO) approximation within the NEO framework is presented. In this approach, the electronic density is quenched to the ground state at each time step, and the real-time nuclear quantum dynamics is propagated on an instantaneous electronic ground state defined by both the classical nuclear geometry and the nonequilibrium quantum nuclear density. Because the electronic dynamics is no longer propagated, this approximation enables the use of an order-of-magnitude larger time step, thus greatly reducing the computational cost. Moreover, invoking the electronic BO approximation also fixes the unphysical asymmetric Rabi splitting observed in previous semiclassical RT-NEO-TDDFT simulations of vibrational polaritons even for small Rabi splitting, instead yielding a stable, symmetric Rabi splitting. For intramolecular proton transfer in malonaldehyde, both RT-NEO-Ehrenfest dynamics and its BO counterpart can describe proton delocalization during the real-time nuclear quantum dynamics. Thus, the BO RT-NEO approach provides the foundation for a wide range of chemical and biological applications.


title: ALMA: Automata Learner using Modulo 2 Multiplicity Automata

abstract: We present ALMA (Automata Learner using modulo 2 Multiplicity Automata), a Java-based tool that can learn any automaton accepting regular languages of finite or infinite words with an implementable membership query function. Users can either pass as input their own membership query function, or use the predefined membership query functions for modulo 2 multiplicity automata and non-deterministic B\"uchi automata. While learning, ALMA can output the state of the observation table after every equivalence query, and upon termination, it can output the dimension, transition matrices, and final vector of the learned modulo 2 multiplicity automaton. Users can test whether a word is accepted by performing a membership query on the learned automaton. ALMA follows the polynomial-time learning algorithm of Beimel et. al. (Learning functions represented as multiplicity automata. J. ACM 47(3), 2000), which uses membership and equivalence queries and represents hypotheses using modulo 2 multiplicity automata. ALMA also implements a polynomial-time learning algorithm for strongly unambiguous B\"uchi automata by Angluin et. al. (Strongly unambiguous B\"uchi automata are polynomially predictable with membership queries. CSL 2020), and a minimization algorithm for modulo 2 multiplicity automata by Sakarovitch (Elements of Automata Theory. 2009).


title: Hybrid CGME and TCGME algorithms for large-scale general-form regularization

abstract: Two new hybrid algorithms are proposed for large-scale linear discrete ill-posed problems in general-form regularization. They are both based on Krylov subspace inner-outer iterative algorithms. At each iteration, they need to solve a linear least squares problem. It is proved that the inner linear least squares problems, which are solved by LSQR, become better conditioned as k increases, so LSQR converges faster. We also prove how to choose the stopping tolerance for LSQR in order to guarantee that the computed solutions have the same accuracy with the exact best regularized solutions. Numerical experiments are given to show the effectiveness and efficiency of our new hybrid algorithms, and comparisons are made with the existing algorithm.


title: Abelian and model structures on tame functors

abstract: In this paper, we discuss certain circumstances in which the category of tame functors inherits an abelian category structure with minimal resolutions and a model category structure with minimal cofibrant replacements. We also present a structure theorem for cofibrant objects in the category of tame functors indexed by realizations of posets of dimension $1$ with values in the category of chain complexes in an abelian category whose all objects are projectives. Moreover, we introduce a general technique to generate indecomposable objects in the abelian category of functors indexed by finite posets.


title: Null controllability of one-dimensional barotropic and non-barotropic linearized compressible Navier-Stokes system using one boundary control

abstract: In this article, we study boundary null controllability properties of the linearized compressible Navier-Stokes equations in the interval $(0,2\pi)$ for both barotropic and non-barotropic fluids using only one boundary control. We consider all the possible cases of the act of control for both systems (density, velocity and temperature). These controls are acting on the boundary and are given as the difference of the values at $x=0$ and $x=2\pi$. In this setup, using a boundary control acting in density, we first prove null controllability of both the barotropic and non-barotropic systems at large time in the spaces $(\dot{L}^2(0,2\pi))^2$ and $(\dot{L}^2(0,2\pi))^3$ respectively (where the dot represents functions with mean value zero). When the control is acting in the velocity component, we prove null controllability at large time in the spaces $\dot{H}^1_{\text{per}}(0,2\pi)\times\dot{L}^2(0,2\pi)$ and $\dot{H}^1_{\text{per}}(0,2\pi)\times(\dot{L}^2(0,2\pi))^2$ respectively. Further, in both cases, we prove that these null controllability results are sharp with respect to the regularity of the initial states in velocity/ temperature case, and time in the density case. Finally, for both barotropic and non-barotropic fluids, we prove that, under some assumptions, the system cannot be approximately controllable at any time, whether there is a control acting in density, velocity or temperature.


title: Two Different Weak Modulations in ab-type RR Lyrae Variable V838 Cyg, and Potential Influence of Metal Abundance on Blazhko Modulation

abstract: Noting the weakest modulation and relatively high metal abundance of the ab-type RR Lyrae star V838 Cyg, we collected the photometric data of this star from several sky surveys to carry out an in-depth analysis. The O-C diagram shows that the pulsation period of V838 Cyg increases linearly over a long timescale. In a reanalysis of the high-precision Kepler data, we confirmed the modulation with a period of 59.45\pm0.07 days found by Benko et al., (2014), and also found an additional weak modulation with a longer period (840\pm21 days). After a series of analyses, we incline to the view that the mechanisms causing the two modulations are different: the former is more similar to the typical Blazhko effect, while the mechanism leading to the latter may be an extrinsic factor. We also collected and compared the modulation and physical parameters of other Blazhko RR Lyrae stars from several works in the literature, and find that there is a potential negative correlation between the modulation amplitude (or upper limit of amplitude) and the metal abundance. We infer that the relatively high metal abundance will promote convection in the outer stellar atmosphere, and then inhibit those factors (turbulence, shock wave, etc.) that may cause Blazhko modulation. Future observations and research work can be carried out with reference to this viewpoint. We also introduce the moire effects that appear in the Kepler long-cadence light curves and their possible interference in the previous analyses.


title: A novel way of calculating scattering integrals

abstract: The technique coined as NDIM - Negative Dimensional Integration Method by their discoverers, relies on a three-pronged basis: Gaussian integration, series expansion and analytic continuation. The technique has been successfully applied to the calculation of covariant and non covariant Feynman integrals in a generic dimensional regularization space, i.e., D-dimensional space-time for D including the negative domain values. Since the dimensionality is general, we can use specifically for one-dimensional integrals. In this work we show how this technique can be applied to tackle certain improper integrals and give an example of a particular improper integral that appears in quantum mechanical scattering process. Traditionally, improper integrals are ascribed certain values through the limiting approach or as is known, by the Cauchy principal value via residues concept technique. Here we use the NDIM approach to do the calculations and show it works fine for the improper integrals. This novel approach we believe is more straightforward and does not require to handle poles, residues, or difficult closed contours as in the traditional approach.


title: Geometry of the space of monodromy data

abstract: In a paper published by the Annales de la Facult\'e de Sciences de Toulouse, with Yousuke Ohyama, we defined and studied a space of monodromy data underlying the well known derivation of q-Painlev\'e VI equation from q-isomonodromy conditions by Jimbo and Sakai. In a recent ArXiv preprint, Nalini Joshi and Pieter Roffelsen pursued our work. However, both our article and their preprint are ambiguous on some foundational algebro-geometric matters. We proceed here to provide sound bases.


title: The dimension of harmonic measure on some AD-regular flat sets of fractional dimension

abstract: In this paper it is shown that if $E\subset\mathbb R^{n+1}$ is an $s$-AD regular compact set, with $s\in [n-\frac12,n)$, and $E$ is contained in a hyperplane or, more generally, in an $n$-dimensional $C^1$ manifold, then the Hausdorff dimension of the harmonic measure for the domain $\mathbb R^{n+1}\setminus E$ is strictly smaller than $s$, i.e., than the Hausdorff dimension of $E$.


title: Propagating photo-$z$ uncertainties: a functional derivative approach

abstract: Photometric redshifts are a key ingredient in the analysis and interpretation of large-scale structure (LSS) surveys. The accuracy and precision of these redshift estimates are directly linked to the constraining power of photometric surveys. It is hence necessary to define precision and accuracy requirements for the redshift calibration \revision{to not} infer biased results in the final analysis. For weak gravitational lensing of the LSS, the photometry culminates in the estimation of the source redshift distribution (SRD) in each of the tomographic bins used in the analysis. The focus has been on shifts of the mean of the SRDs and how well the calibration must be able to recover those. Since the estimated SRDs are usually given as a normalized histogram with corresponding errors, it would be advantageous to propagate these uncertainties accordingly to see whether the requirements of the given survey are indeed fulfilled. Here we propose the use of functional derivatives to calculate the sensitivity of the final observables, e.g. the lensing angular power spectrum, with respect to the SRD at a specific redshift. This allows the propagation of arbitrarily shaped small perturbations to the SRD, without having to run the whole analysis pipeline for each realization again. We apply our method to an EUCLID survey and demonstrate it with SRDs of the KV450 data set, recovering previous results. Lastly, we note that the moments of the SRD of order larger than two will probably not be relevant when propagating redshift uncertainties in cosmic shear analysis.


title: On the Comparison of Discounted-Sum Automata with Multiple Discount Factors

abstract: We look into the problems of comparing nondeterministic discounted-sum automata on finite and infinite words. That is, the problems of checking for automata $A$ and $B$ whether or not it holds that for all words $w$, $A(w)=B(w), A(w) \leq B(w)$, or $A(w)<B(w)$. These problems are known to be decidable when both automata have the same single integral discount factor, while decidability is open in all other settings: when the single discount factor is a non-integral rational; when each automaton can have multiple discount factors; and even when each has a single integral discount factor, but the two are different. We show that it is undecidable to compare discounted-sum automata with multiple discount factors, even if all are integrals, while it is decidable to compare them if each has a single, possibly different, integral discount factor. To this end, we also provide algorithms to check for given nondeterministic automaton $N$ and deterministic automaton $D$, each with a single, possibly different, rational discount factor, whether or not $N(w) = D(w)$, $N(w) \geq D(w)$, or $N(w) > D(w)$ for all words $w$.


title: The Low-Redshift Lyman Continuum Survey: Optically Thin and Thick Mg II Lines as Probes of Lyman Continuum Escape

abstract: The Mg II 2796, 2803 doublet has been suggested to be a useful indirect indicator for the escape of Ly-alpha and Lyman continuum (LyC) photons in local star-forming galaxies. However, studies to date have focused on small samples of galaxies with strong Mg II or strong LyC emission. Here we present the first study of Mg II probing a large dynamic range of galaxy properties, using newly obtained high signal-to-noise, moderate-resolution spectra of Mg II for a sample of 34 galaxies selected from the Low-redshift Lyman Continuum Survey. We show that the galaxies in our sample have Mg II profiles ranging from strong emission to P-Cygni profiles, and to pure absorption. We find there is a significant trend (with a possibility of spurious correlations of ~ 2%) that galaxies detected as strong LyC Emitters (LCEs) also show larger equivalent widths of Mg II emission, and non-LCEs tend to show evidence of more scattering and absorption features in Mg II We then find Mg II strongly correlates with Ly-alpha in both equivalent width and escape fraction, regardless of whether the emission or absorption dominates the Mg II profiles. Furthermore, we present that, for galaxies categorized as Mg II emitters (MgE), one can adopt the information of Mg II, metallicity, and dust to estimate the escape fraction of LyC within a factor of 3. These findings confirm that Mg II lines can be used as a tool to select galaxies as LCEs and to serve as an indirect indicator for the escape of Ly-alpha and LyC.


title: Community Detection with Known, Unknown, or Partially Known Auxiliary Latent Variables

abstract: Empirical observations suggest that in practice, community membership does not completely explain the dependency between the edges of an observation graph. The residual dependence of the graph edges are modeled in this paper, to first order, by auxiliary node latent variables that affect the statistics of the graph edges but carry no information about the communities of interest. We then study community detection in graphs obeying the stochastic block model and censored block model with auxiliary latent variables. We analyze the conditions for exact recovery when these auxiliary latent variables are unknown, representing unknown nuisance parameters or model mismatch. We also analyze exact recovery when these secondary latent variables have been either fully or partially revealed. Finally, we propose a semidefinite programming algorithm for recovering the desired labels when the secondary labels are either known or unknown. We show that exact recovery is possible by semidefinite programming down to the respective maximum likelihood exact recovery threshold.


title: An ultra-stable three-dimensional photophoretic trap in air facilitated by a single multimode fiber

abstract: Photophoretic forces - which are of thermal origin - have defined an alternative route of optical trapping of absorbing microparticles in air. Here, we show that a single multi-mode fiber facilitates significantly more robust optical traps compared to a pure Gaussian beam emanating from a single mode fiber for the trapping and manipulation of absorbing particles using photophoretic forces. We carefully study the dependency of trapping on speckle patterns generated from different modes from a multimode fiber, and experimentally observe that maximum trapping force can be obtained when the mean speckle size is comparable to the diameter of a trapped particle. We explain this observation by numerical simulations carried out to calculate the photophoretic force, and also determine stable trapping conditions from force balance equations. Interestingly, we also observe large oscillations of the trapped particle along the z-direction for multimode beams, which may be demonstrative of an effective restoring force for photophoretic trapping even in the axial direction. Our work may presage a new route for exciting applications on optical trapping and spectroscopy with photophoretic forces due to the inherent ease-of-use, portability, and flexibility of single muti-mode fiber based optical traps.


title: Finding Nontrivial Minimum Fixed Points in Discrete Dynamical Systems

abstract: Networked discrete dynamical systems are often used to model the spread of contagions and decision-making by agents in coordination games. Fixed points of such dynamical systems represent configurations to which the system converges. In the dissemination of undesirable contagions (such as rumors and misinformation), convergence to fixed points with a small number of affected nodes is a desirable goal. Motivated by such considerations, we formulate a novel optimization problem of finding a nontrivial fixed point of the system with the minimum number of affected nodes. We establish that, unless P = NP, there is no polynomial time algorithm for approximating a solution to this problem to within the factor n^1-\epsilon for any constant epsilon > 0. To cope with this computational intractability, we identify several special cases for which the problem can be solved efficiently. Further, we introduce an integer linear program to address the problem for networks of reasonable sizes. For solving the problem on larger networks, we propose a general heuristic framework along with greedy selection methods. Extensive experimental results on real-world networks demonstrate the effectiveness of the proposed heuristics.


title: Complex balanced distributions for chemical reaction networks

abstract: Stationary distributions of continuous time Markov chains (CTMCs) are often a main interest, but hard to find. We consider CTMCs modelling reaction networks, and characterise complex balanced distributions (provided they exist) of reaction networks with arbitrary transition functions through conditions on the cycles of their corresponding digraph. The proof works by constructing a dynamically equivalent reaction network of disjoint cycles. We further derive a sufficient condition for the existence of a complex balanced distribution, and give precise conditions on when it is necessary. The sufficient condition holds for mass-action kinetics or if the digraph consists of only cyclic connected components. Hence, we fully characterise the existence and form of complex balanced distributions. Moreover, these results can be used to design complex balanced reaction networks with a given stationary distribution.


title: Exceptional points for associated Legendre functions of the second kind

abstract: We consider the complex $\nu$ plane structure of the associated Legendre function of the second kind $Q^{-1/2-K}_{\nu}(\cosh\rho)$. We find that for any noninteger value for $K$ $Q^{-1/2-K}_{\nu}(\cosh\rho)$ has an infinite number of poles in the complex $\nu$ plane, but for any negative integer $K$ there are no poles at all. For $K=0$ or any positive integer $K$ there is only a finite number of poles, with there only being one single pole (at $\nu=0$) when $K=0$. This pattern is characteristic of the exceptional points that appear in a wide variety of physical contexts. However, unusually for theories with exceptional points, $Q^{-1/2-K}_{\nu}(\cosh\rho)$ has an infinite number of them. Other than in the $PT$-symmetry Jordan-block case, exceptional points usually occur at complex values of parameters. While not being Jordan-block exceptional points themselves, the exceptional points associated with the $Q^{-1/2-K}_{\nu}(\cosh\rho)$ nonetheless occur at real values of $K$.


title: On the Robustness of AlphaFold: A COVID-19 Case Study

abstract: Protein folding neural networks (PFNNs) such as AlphaFold predict remarkably accurate structures of proteins compared to other approaches. However, the robustness of such networks has heretofore not been explored. This is particularly relevant given the broad social implications of such technologies and the fact that biologically small perturbations in the protein sequence do not generally lead to drastic changes in the protein structure. In this paper, we demonstrate that AlphaFold does not exhibit such robustness despite its high accuracy. This raises the challenge of detecting and quantifying the extent to which these predicted protein structures can be trusted. To measure the robustness of the predicted structures, we utilize (i) the root-mean-square deviation (RMSD) and (ii) the Global Distance Test (GDT) similarity measure between the predicted structure of the original sequence and the structure of its adversarially perturbed version. We prove that the problem of minimally perturbing protein sequences to fool protein folding neural networks is NP-complete. Based on the well-established BLOSUM62 sequence alignment scoring matrix, we generate adversarial protein sequences and show that the RMSD between the predicted protein structure and the structure of the original sequence are very large when the adversarial changes are bounded by (i) 20 units in the BLOSUM62 distance, and (ii) five residues (out of hundreds or thousands of residues) in the given protein sequence. In our experimental evaluation, we consider 111 COVID-19 proteins in the Universal Protein resource (UniProt), a central resource for protein data managed by the European Bioinformatics Institute, Swiss Institute of Bioinformatics, and the US Protein Information Resource. These result in an overall GDT similarity test score average of around 34%, demonstrating a substantial drop in the performance of AlphaFold.


title: One-particle density matrix and momentum distribution of the out-of-equilibrium 1D Tonks-Girardeau gas: Analytical results at large $N$

abstract: In one-dimensional (1D) quantum gases, the momentum distribution (MD) of the atoms is a standard experimental observable, routinely measured in various experimental setups. The MD is sensitive to correlations, and it is notoriously hard to compute theoretically for large numbers of atoms $N$, which often prevents direct comparison with experimental data. Here we report significant progress on this problem for the 1D Tonks-Girardeau (TG) gas in the asymptotic limit of large $N$, at zero temperature and driven out of equilibrium by a quench of the confining potential. We find an exact analytical formula for the one-particle density matrix $\langle \hat{\Psi}^\dagger(x) \hat{\Psi}(x') \rangle$ of the out-of-equilibrium TG gas in the $N \rightarrow \infty$ limit, valid on distances $|x-x'| $ much larger than the interparticle distance. By comparing with time-dependent Bose-Fermi mapping numerics, we demonstrate that our analytical formula can be used to compute the out-of-equilibrium MD with great accuracy for a wide range of momenta (except in the tails of the distribution at very large momenta). For a quench from a double-well potential to a single harmonic well,which mimics a `quantum Newton cradle' setup, our method predicts the periodic formation of peculiar, multiply peaked, momentum distributions.


title: Optimal randomized multilevel Monte Carlo for repeatedly nested expectations

abstract: The estimation of repeatedly nested expectations is a challenging task that arises in many real-world systems. However, existing methods generally suffer from high computational costs when the number of nestings becomes large. Fix any non-negative integer $D$ for the total number of nestings. Standard Monte Carlo methods typically cost at least $\mathcal{O}(\varepsilon^{-(2+D)})$ and sometimes $\mathcal{O}(\varepsilon^{-2(1+D)})$ to obtain an estimator up to $\varepsilon$-error. More advanced methods, such as multilevel Monte Carlo, currently only exist for $D = 1$. In this paper, we propose a novel Monte Carlo estimator called $\mathsf{READ}$, which stands for "Recursive Estimator for Arbitrary Depth.'' Our estimator has an optimal computational cost of $\mathcal{O}(\varepsilon^{-2})$ for every fixed $D$ under suitable assumptions, and a nearly optimal computational cost of $\mathcal{O}(\varepsilon^{-2(1 + \delta)})$ for any $0 < \delta < \frac12$ under much more general assumptions. Our estimator is also unbiased, which makes it easy to parallelize. The key ingredients in our construction are an observation of the problem's recursive structure and the recursive use of the randomized multilevel Monte Carlo method.


title: Probing the longitudinal spin transfer via dihadron polarization correlations in unpolarized $e^+e^-$ and $pp$ collisions

abstract: The longitudinal spin transfer represents the probability density of producing longitudinally polarized hadrons from longitudinally polarized quarks or circularly polarized gluons. It thus was usually measured in polarized reactions or high-energy collisions where weak interaction dominates. In this work, we propose the dihadron polarization correlation as a novel probe of this quantity. Such an observable does not require the fragmenting partons to be polarized and therefore can be measured in the currently available experimental facilities, such as Belle, RHIC, Tevatron, and the LHC. We make quantitative predictions for these experiments. In light of the data already harvested, the experimental investigation of this observable provides more opportunity for the quantitative study of the longitudinal spin transfer. In particular, the measurements in pp collisions can significantly constrain the fragmentation function of a circularly polarized gluon.


title: Stability of Hardy-Littlewood-Sobolev inequalities with explicit lower bounds

abstract: In this paper, we establish the stability for the Hardy-Littlewood-Sobolev (HLS) inequalities with explicit lower bounds. By establishing the relation between the stability of HLS inequalities and the stability of fractional Sobolev inequalities, we also give the stability of the fractional Sobolev inequalities with the lower bounds. This extends the stability of Sobolev inequalities with the explicit lower bounds established by Dolbeault, Esteban, Figalli, Frank and Loss in [16] to the fractional order case. Our proofs are based on the competing symmetries, the continuous Steiner symmetrization inequality for the HLS integral and the dual stability theory.


title: Investigating Conversational Search Behavior For Domain Exploration

abstract: Conversational search has evolved as a new information retrieval paradigm, marking a shift from traditional search systems towards interactive dialogues with intelligent search agents. This change especially affects exploratory information-seeking contexts, where conversational search systems can guide the discovery of unfamiliar domains. In these scenarios, users find it often difficult to express their information goals due to insufficient background knowledge. Conversational interfaces can provide assistance by eliciting information needs and narrowing down the search space. However, due to the complexity of information-seeking behavior, the design of conversational interfaces for retrieving information remains a great challenge. Although prior work has employed user studies to empirically ground the system design, most existing studies are limited to well-defined search tasks or known domains, thus being less exploratory in nature. Therefore, we conducted a laboratory study to investigate open-ended search behavior for navigation through unknown information landscapes. The study comprised of 26 participants who were restricted in their search to a text chat interface. Based on the collected dialogue transcripts, we applied statistical analyses and process mining techniques to uncover general information-seeking patterns across five different domains. We not only identify core dialogue acts and their interrelations that enable users to discover domain knowledge, but also derive design suggestions for conversational search systems.


title: (Anti-)Stokes Scattering on Kinks

abstract: At leading order, there are three inelastic scattering processes beginning with a quantum kink and a fundamental meson. Meson multiplication, in which the final state is a kink and two mesons, was treated recently. In this note we treat the other two, (anti)-Stokes scattering, in which the kink's shape mode is (de-)excited and the final state contains one meson. In the case of a general scalar kink, we find analytic formulas for the forward and backward scattering amplitudes and probabilities as functions of the momentum of the incident meson. The general results are then specialized to the kink of the $\phi^4$ double-well model.


title: Triggered Superradiance and Inversion Storage in a Hybrid Quantum System

abstract: We study the superradiant emission of an inverted spin ensemble strongly coupled to a superconducting cavity. After fast inversion, we detune the spins from the cavity and store the inversion for tens of milliseconds, during which the remaining transverse spin components disappear. Switching back on resonance enables us to study the onset of superradiance. A weak trigger pulse of a few hundred photons shifts the superradiant burst to earlier times and imprints its phase onto the emitted radiation. For long hold times, the inversion decreases below the threshold for spontaneous superradiance. There, the energy stored in the ensemble can be used to amplify microwave pulses passing through the cavity.


title: Neural Radiance Field Codebooks

abstract: Compositional representations of the world are a promising step towards enabling high-level scene understanding and efficient transfer to downstream tasks. Learning such representations for complex scenes and tasks remains an open challenge. Towards this goal, we introduce Neural Radiance Field Codebooks (NRC), a scalable method for learning object-centric representations through novel view reconstruction. NRC learns to reconstruct scenes from novel views using a dictionary of object codes which are decoded through a volumetric renderer. This enables the discovery of reoccurring visual and geometric patterns across scenes which are transferable to downstream tasks. We show that NRC representations transfer well to object navigation in THOR, outperforming 2D and 3D representation learning methods by 3.1% success rate. We demonstrate that our approach is able to perform unsupervised segmentation for more complex synthetic (THOR) and real scenes (NYU Depth) better than prior methods (29% relative improvement). Finally, we show that NRC improves on the task of depth ordering by 5.5% accuracy in THOR.


title: Effects of porous substrates on the structure of turbulent boundary layers

abstract: Three different porous substrates (with different pore sizes, s, and permeabilities, K) are used to examine their effect on the structure of boundary layer flow over them. The flow is characterised with single-point hot-wire measurements as well as planar Particle Image Velocimetry. In order to elucidate differences in shallow and deep flows past porous substrate, foams with two different thickness (h) are used (for all three substrates). A wide range of Friction Reynolds number (2000< Retau < 15000) and Permeability based Reynolds number (1<ReK< 50) are attained. For substrates with ReK=1, the flow behaviour remains similar to flow over impermeable smooth walls and as such Townsend's hypothesis remains valid. In contrast, a substantial reduction in velocity disturbances and associated length scales is achieved for permeable (ReK>1) and dense (relative to viscous scales) foam at the thick substrate limit (h/s>10), which leads to the breakdown of outer-layer similarity. As porosity is increased, a thin substrate limit is reached (h/s), and the foam becomes sparse relative to viscous scales (s+ > 100). For such foams, the flow conforms to outer-layer similarity and is more akin to flow over rough surfaces. Such substrates are unable to attenuate velocity disturbances and the dependence of substrate thickness (h/s) on spectral energy content of turbulent fluctuations ceases to exist. The present study shows that transition from thick to thin substrate flow behavior depends not only on thickness-to-pore ratio (h/s) but also on substrate density relative to viscous scales of the flow.


title: Meteorite Parent Body Aqueous Alteration Simulations of Interstellar Residue Analogs

abstract: Some families of carbonaceous chondrites are rich in prebiotic organics that may have contributed to the origin of life on Earth and elsewhere. However, the formation and chemical evolution of complex soluble organic molecules from interstellar precursors under relevant parent body conditions has not been thoroughly investigated. In this study, we approach this topic by simulating meteorite parent body aqueous alteration of interstellar residue analogs. The distributions of amines and amino acids are qualitatively and quantitatively investigated and linked to closing the gap between interstellar and meteoritic prebiotic organic abundances. We find that the abundance trend of methylamine > ethylamine> glycine > serine > alanine > \b{eta}-alanine does not change from pre- to post-aqueous alteration, suggesting that certain cloud conditions have an influential role on the distributions of interstellar-inherited meteoritic organics. However, the abundances for most of the amines and amino acids studied here varied by about 2-fold when aqueously processed for 7 days at 125 {\deg}C, and the changes in the {\alpha}- to \b{eta}-alanine ratio were consistent with those of aqueously altered carbonaceous chondrites, pointing to an influential role of meteorite parent body processing on the distributions of interstellar-inherited meteoritic organics. We detected higher abundances of {\alpha}- over \b{eta}-alanine, which is opposite to what is typically observed in aqueously altered carbonaceous chondrites; these results may be explained by at least the lack of minerals and insoluble organic matter-relevant materials in the experiments. The high abundance of volatile amines in the non-aqueously altered samples suggests that these types of interstellar volatiles can be efficiently transferred to asteroids and comets, supporting the idea of the presence of interstellar organics in solar system objects.

categories: astro-ph.EP astro-ph.GA astro-ph.SR physics.geo-ph


title: Mastering Diverse Domains through World Models

abstract: Developing a general algorithm that learns to solve tasks across a wide range of applications has been a fundamental challenge in artificial intelligence. Although current reinforcement learning algorithms can be readily applied to tasks similar to what they have been developed for, configuring them for new application domains requires significant human expertise and experimentation. We present DreamerV3, a general algorithm that outperforms specialized methods across over 150 diverse tasks, with a single configuration. Dreamer learns a model of the environment and improves its behavior by imagining future scenarios. Robustness techniques based on normalization, balancing, and transformations enable stable learning across domains. Applied out of the box, Dreamer is the first algorithm to collect diamonds in Minecraft from scratch without human data or curricula. This achievement has been posed as a significant challenge in artificial intelligence that requires exploring farsighted strategies from pixels and sparse rewards in an open world. Our work allows solving challenging control problems without extensive experimentation, making reinforcement learning broadly applicable.


title: Commensuration torques and lubricity in double moire systems

abstract: We study the commensuration torques and layer sliding energetics of alternating twist trilayer graphene (t3G) and twisted bilayer graphene on hexagonal boron nitride (t2G/BN) that have two superposed moire interfaces. Lattice relaxations for typical graphene twist angles of $\sim 1^{\circ}$ in t3G or t2G/BN are found to break the out-of-plane layer mirror symmetry, give rise to layer rotation energy local minima dips of the order of $\sim 10^{-1}$ meV/atom at double moire alignment angles, and have sliding energy landscape minima between top-bottom layers of comparable magnitude. Moire superlubricity is restored for twist angles as small as $\sim 0.03^\circ$ away from alignment resulting in suppression of sliding energies by several orders of magnitude of typically $\sim 10^{-4}$ meV/atom, hence indicating the precedence of rotation over sliding in the double moire commensuration process.


title: Simulation of ODMR Spectra from Nitrogen-Vacancy Ensembles in Diamond for Electric Field Sensing

abstract: Solid state spins in diamond, in particular negatively charged nitrogen-vacancy centers (NV), are leading contenders in the field of quantum sensing. While addressing of single NVs offers nanoscale spatial resolution, many implementations benefit from using large ensembles to increase signal magnitude and therefore sensitivity. However, sensing with ensembles brings its own challenges given the random orientation of the spin quantization axis within the diamond crystal lattice. Here, we present an open source simulation tool that models the influence of arbitrary electric and magnetic fields on the electronic and nuclear spin states of NV ensembles, and can be extended to other color centers. Specifically, the code computes the transition strengths and predicts the sensitivity under shot-noise-limited optically-detected magnetic resonance. We illustrate the use of the code in the context of electric field sensing, a promising emerging functionality of NV centers with applications in biosensing and electronics, and bring several subtle features to light that are due to the interplay between different NV orientations and the external electric and microwave fields. Moreover, we show that our code can be used to optimize sensitivity in situations where usual arguments based on neglecting terms in the full Hamiltonian would give sub-optimal results. Finally, we propose a novel sensing scheme which allows to perform full vector electrometry without the need for precise bias magnetic field alignment, thus reducing the experimental complexity and speeding up the measurement procedure.


title: Interpolating Matrix Models for WLZZ series

abstract: We suggest a two-matrix model depending on three (infinite) sets of parameters which interpolates between all the models proposed in arXiv:2206.13038, and defined there through $W$-representations. We also discuss further generalizations of these WLZZ models realized by $W$-representations associated with infinite commutative families of generators of $w_\infty$-algebra which are presumably related to more sophisticated multi-matrix models. Integrable properties of these generalizations are described by what we call the skew hypergeometric $\tau$-functions.


title: Functional observability and subspace reconstruction in nonlinear systems

abstract: Time-series analysis is fundamental for modeling and predicting dynamical behaviors from time-ordered data, with applications in many disciplines such as physics, biology, finance, and engineering. Measured time-series data, however, are often low dimensional or even univariate, thus requiring embedding methods to reconstruct the original system's state space. The observability of a system establishes fundamental conditions under which such reconstruction is possible. However, complete observability is too restrictive in applications where reconstructing the entire state space is not necessary and only a specific subspace is relevant. Here, we establish the theoretic condition to reconstruct a nonlinear functional of state variables from measurement processes, generalizing the concept of functional observability to nonlinear systems. When the functional observability condition holds, we show how to construct a map from the embedding space to the desired functional of state variables, characterizing the quality of such reconstruction. The theoretical results are then illustrated numerically using chaotic systems with contrasting observability properties. By exploring the presence of functionally unobservable regions in embedded attractors, we also apply our theory for the early warning of seizure-like events in simulated and empirical data. The studies demonstrate that the proposed functional observability condition can be assessed a priori to guide time-series analysis and experimental design for the dynamical characterization of complex systems.


title: Matching calipers and the precision of index estimation

abstract: This paper characterizes the precision of index estimation as it carries over into precision of matching. In a model assuming Gaussian covariates and making best-case assumptions about matching quality, it sharply characterizes average and worst-case discrepancies between paired differences of true versus estimated index values. In this optimistic setting, worst-case true and estimated index differences decline to zero if $p=o[n/(\log n)]$, the same restriction on model size that is needed for consistency of common index models. This remains so as the Gaussian assumption is relaxed to sub-gaussian, if in that case the characterization of paired index errors is less sharp. The formula derived under Gaussian assumptions is used as the basis for a matching caliper. Matching such that paired differences on the estimated index fall below this caliper brings the benefit that after matching, worst-case differences onan underlying index tend to 0 if $p = o\{[n/(\log n)]^{2/3}\}$. (With a linear index model, $p=o[n/(\log n)]$ suffices.) A proposed refinement of the caliper condition brings the same benefits without the sub-gaussian condition on covariates. When strong ignorability holds and the index is a well-specified propensity or prognostic score, ensuring in this way that worst-case matched discrepancies on it tend to 0 with increasing $n$ also ensures the consistency of matched estimators of the treatment effect.