title: Braiding groups of automorphisms and almost-automorphisms of trees

abstract: We introduce "braided" versions of self-similar groups and R\"over--Nekrashevych groups, and study their finiteness properties. This generalizes work of Aroca and Cumplido, and the first author and Wu, who considered the case when the self-similar groups are what we call "self-identical". In particular we use a braided version of the Grigorchuk group to construct a new group called the braided R\"over group, which we prove is of type $F_\infty$. Our techniques involve using so called $d$-ary cloning systems to construct the groups, and analyzing certain complexes of embedded disks in a surface to understand their finiteness properties.


title: Improved Beckner's inequality for axially symmetric functions on $\mathbb{S}^4$

abstract: We show that axially symmetric solutions on $\mathbb{S}^4$ to a constant $Q$-curvature type equation (it may also be called fourth order mean field equation) must be constant, provided that the parameter $\alpha$ in front of the Paneitz operator belongs to $[\frac{473 + \sqrt{209329}}{1800}\approx0.517, 1)$. This is in contrast to the case $\alpha=1$, where a family of solutions exist, known as standard bubbles. The phenomenon resembles the Gaussian curvature equation on $ \mathbb{S}^2$. As a consequence, we prove an improved Beckner's inequality on $\mathbb{S}^4$ for axially symmetric functions with their centers of mass at the origin. Furthermore, we show uniqueness of axially symmetric solutions when $\alpha=\frac15$ by exploiting Pohozaev-type identities, and prove existence of a non-constant axially symmetric solution for $\alpha \in (\frac15, \frac12)$ via a bifurcation method.


title: Curvature-Aware Derivative-Free Optimization

abstract: The paper discusses derivative-free optimization (DFO), which involves minimizing a function without access to gradients or directional derivatives, only function evaluations. Classical DFO methods, which mimic gradient-based methods, such as Nelder-Mead and direct search have limited scalability for high-dimensional problems. Zeroth-order methods have been gaining popularity due to the demands of large-scale machine learning applications, and the paper focuses on the selection of the step size $\alpha_k$ in these methods. The proposed approach, called Curvature-Aware Random Search (CARS), uses first- and second-order finite difference approximations to compute a candidate $\alpha_{+}$. We prove that for strongly convex objective functions, CARS converges linearly provided that the search direction is drawn from a distribution satisfying very mild conditions. We also present a Cubic Regularized variant of CARS, named CARS-CR, which converges in a rate of $\mathcal{O}(k^{-1})$ without the assumption of strong convexity. Numerical experiments show that CARS and CARS-CR match or exceed the state-of-the-arts on benchmark problem sets.


title: The Tensor Brain: A Unified Theory of Perception, Memory and Semantic Decoding

abstract: We present a unified computational theory of an agent's perception and memory. In our model, perception, episodic memory, and semantic memory are realized by different operational modes of the oscillating interactions between a symbolic index layer and a subsymbolic representation layer. The two layers form a bilayer tensor network (BTN). Although memory appears to be about the past, its main purpose is to support the agent in the present and the future. Recent episodic memory provides the agent with a sense of the here and now. Remote episodic memory retrieves relevant past experiences to provide information about possible future scenarios. This aids the agent in decision-making. "Future" episodic memory, based on expected future events, guides planning and action. Semantic memory retrieves specific information, which is not delivered by current perception, and defines priors for future observations. We argue that it is important for the agent to encode individual entities, not just classes and attributes. We demonstrate that a form of self-supervised learning can acquire new concepts and refine existing ones. We test our model on a standard benchmark data set, which we expanded to contain richer representations for attributes, classes, and individuals. Our key hypothesis is that obtaining a better understanding of perception and memory is a crucial prerequisite to comprehending human-level intelligence.


title: Invertibility of Positive Toeplitz Operators and Associated Uncertainty Principle

abstract: We study invertibility and compactness of positive Toeplitz operators associated to a continuous Parseval frame on a Hilbert space. As applications, we characterize compactness of affine and Weyl-Heisenberg localization operators as well as give uncertainty principles for the associated transforms.


title: Compact Redistricting Plans Have Many Spanning Trees

abstract: In the design and analysis of political redistricting maps, it is often useful to be able to sample from the space of all partitions of the graph of census blocks into connected subgraphs of equal population. There are influential Markov chain Monte Carlo methods for doing so that are based on sampling and splitting random spanning trees. Empirical evidence suggests that the distributions such algorithms sample from place higher weight on more "compact" redistricting plans, which is a practically useful and desirable property. In this paper, we confirm these observations analytically, establishing an inverse exponential relationship between the total length of the boundaries separating districts and the probability that such a map will be sampled. This result provides theoretical underpinnings for algorithms that are already making a significant real-world impact.


title: A Symmetry-orientated Divide-and-Conquer Method for Crystal Structure Prediction

abstract: Crystal structure prediction has been a subject of topical interest, but remains a substantial challenge, especially for complex structures as it deals with the global minimization of the extremely rugged high-dimensional potential energy surface. In this manuscript, a symmetry-orientated divide-and-conquer scheme was proposed to construct a symmetry tree graph, where the entire search space is decomposed into a finite number of symmetry-dependent subspaces. An artificial intelligence-based symmetry selection strategy was subsequently devised to select the low-lying subspaces with high symmetries for global exploration and in-depth exploitation. Our approach can significantly simplify the problem of crystal structure prediction by avoiding exploration of the most complex P1 subspace on the entire search space and have the advantage of preserving the crystal symmetry during structure evolution, making it well suitable for predicting the complex crystal structures. The effectiveness of the method has been validated by successful prediction of the candidate structures of binary Lennard-Jones mixtures and high-pressure phase of ice, containing more than one hundred atoms in the simulation cell. The work, therefore, opens up an opportunity towards achieving the long-sought goal for crystal structure prediction of complex systems.


title: Bridge Data: Boosting Generalization of Robotic Skills with Cross-Domain Datasets

abstract: Robot learning holds the promise of learning policies that generalize broadly. However, such generalization requires sufficiently diverse datasets of the task of interest, which can be prohibitively expensive to collect. In other fields, such as computer vision, it is common to utilize shared, reusable datasets, such as ImageNet, to overcome this challenge, but this has proven difficult in robotics. In this paper, we ask: what would it take to enable practical data reuse in robotics for end-to-end skill learning? We hypothesize that the key is to use datasets with multiple tasks and multiple domains, such that a new user that wants to train their robot to perform a new task in a new domain can include this dataset in their training process and benefit from cross-task and cross-domain generalization. To evaluate this hypothesis, we collect a large multi-domain and multi-task dataset, with 7,200 demonstrations constituting 71 tasks across 10 environments, and empirically study how this data can improve the learning of new tasks in new environments. We find that jointly training with the proposed dataset and 50 demonstrations of a never-before-seen task in a new domain on average leads to a 2x improvement in success rate compared to using target domain data alone. We also find that data for only a few tasks in a new domain can bridge the domain gap and make it possible for a robot to perform a variety of prior tasks that were only seen in other domains. These results suggest that reusing diverse multi-task and multi-domain datasets, including our open-source dataset, may pave the way for broader robot generalization, eliminating the need to re-collect data for each new robot learning project.


title: A 4-dimensional light bulb theorem for disks

abstract: We give a 4-dimensional light bulb theorem for properly embedded disks, generalizing recent work of Gabai and Kosanovic-Teichner in certain contexts, and extending the 4-dimensional light bulb theorem for 2-spheres due to Gabai and Schneiderman-Teichner. In particular, we provide conditions under which homotopic disks properly embedded in a compact 4-manifold X with a common dual in the interior of X are smoothly isotopic rel boundary. We also provide a new geometric interpretation of the Dax invariant, to aid in its computation.


title: Unrolling SGD: Understanding Factors Influencing Machine Unlearning

abstract: Machine unlearning is the process through which a deployed machine learning model is made to forget about some of its training data points. While naively retraining the model from scratch is an option, it is almost always associated with large computational overheads for deep learning models. Thus, several approaches to approximately unlearn have been proposed along with corresponding metrics that formalize what it means for a model to forget about a data point. In this work, we first taxonomize approaches and metrics of approximate unlearning. As a result, we identify verification error, i.e., the L2 difference between the weights of an approximately unlearned and a naively retrained model, as an approximate unlearning metric that should be optimized for as it subsumes a large class of other metrics. We theoretically analyze the canonical training algorithm, stochastic gradient descent (SGD), to surface the variables which are relevant to reducing the verification error of approximate unlearning for SGD. From this analysis, we first derive an easy-to-compute proxy for verification error (termed unlearning error). The analysis also informs the design of a new training objective penalty that limits the overall change in weights during SGD and as a result facilitates approximate unlearning with lower verification error. We validate our theoretical work through an empirical evaluation on learning with CIFAR-10, CIFAR-100, and IMDB sentiment analysis.


title: Assessing Outcome-to-Outcome Interference in Sibling Fixed Effects Models

abstract: Sibling fixed effects (FE) models are useful for estimating causal treatment effects while offsetting unobserved sibling-invariant confounding. However, treatment estimates are biased if an individual's outcome affects their sibling's outcome. We propose a robustness test for assessing the presence of outcome-to-outcome interference in linear two-sibling FE models. We regress a gain-score--the difference between siblings' continuous outcomes--on both siblings' treatments and on a pre-treatment observed FE. Under certain restrictions, the observed FE's partial regression coefficient signals the presence of outcome-to-outcome interference. Monte Carlo simulations demonstrated the robustness test under several models. We found that an observed FE signaled outcome-to-outcome spillover if it was directly associated with an sibling-invariant confounder of treatments and outcomes, directly associated with a sibling's treatment, or directly and equally associated with both siblings' outcomes. However, the robustness test collapsed if the observed FE was directly but differentially associated with siblings' outcomes or if outcomes affected siblings' treatments.


title: The energy spectrum of cosmic rays beyond the turn-down around $10^{17}$ eV as measured with the surface detector of the Pierre Auger Observatory

abstract: We present a measurement of the cosmic-ray spectrum above 100\,PeV using the part of the surface detector of the Pierre Auger Observatory that has a spacing of 750~m. An inflection of the spectrum is observed, confirming the presence of the so-called \emph{second-knee} feature. The spectrum is then combined with that of the 1500\,m array to produce a single measurement of the flux, linking this spectral feature with the three additional breaks at the highest energies. The combined spectrum, with an energy scale set calorimetrically via fluorescence telescopes and using a single detector type, results in the most statistically and systematically precise measurement of spectral breaks yet obtained. These measurements are critical for furthering our understanding of the highest energy cosmic rays.


title: Column expansion identities and quadratic spanning forest identities

abstract: Column expansion identities of determinants give a source of quadratic spanning forest polynomial identities and allow us determine the dimension of the space of certain quadratic spanning forest identities, settling a conjecture of one of us with Vlasev from 2012. Furthermore, we give a combinatorial interpretation of such spanning forest identities via an edge-swapping argument previously developed by one of us in 2019. Quadratic spanning forest polynomials identities are of particular interest because they are useful for quantum field theory calculations in four dimensions.


title: Dirac Operators with Operator Data of Wigner-von Neumann Type

abstract: We consider half-line Dirac operators with operator data of Wigner-von Neumann type. If the data is a finite linear combination of Wigner-von Neumann functions, we show absence of singular continuous spectrum and provide an explicit set containing all embedded pure points depending only on the Lp decay and frequencies of the operator data. For infinite sums of Wigner-von Neumann-like terms, we bound the Hausdorff dimension of the singular part of the spectrum.


title: Mass Transfer and Stellar Evolution of the White Dwarfs in AM CVn Binaries

abstract: We calculate the stellar evolution of both white dwarfs (WDs) in AM CVn binaries with orbital periods of $P_{\mathrm{orb}} \approx 5-70$ minutes. We focus on the cases where the donor starts as a $M_{\mathrm{He}} < 0.2 \, M_{\odot}$ Helium WD and the accretor is a $M_{\mathrm{WD}} > 0.6 \, M_{\odot}$ WD. Using Modules for Experiments in Stellar Astrophysics (MESA), we simultaneously evolve both WDs assuming conservative mass transfer and angular momentum loss from gravitational radiation. This self-consistent evolution yields the important feedback of the properties of the donor on the mass transfer rate, $\dot{M}$, as well as the thermal evolution of the accreting WD. Consistent with earlier work, we find that the high $\dot{M}$'s at early times forces an adiabatic evolution of the donor for $P_{\mathrm{orb}} < 30$ minutes so that its mass-radius relation depends primarily on its initial entropy. As the donor reaches $ M_{\mathrm{He}} \approx 0.02-0.03 \, M_{\odot}$ at $P_{\mathrm{orb}} \simeq 30 $ minutes, it becomes fully convective and could lose entropy and expand much less than expected under further mass loss. However, we show that the lack of reliable opacities for the donor's surface inhibit a secure prediction for this possible cooling. Our calculations capture the core heating that occurs during the first $\approx 10^7$ years of accretion and continue the evolution into the phase of WD cooling that follows. When compared to existing data for accreting WDs, as seen by Cheng and collaborators for isolated WDs, we also find that the accreting WDs are not as cool as we would expect given the amount of time they have had to cool.


title: Millimeter Wave and Terahertz Urban Microcell Propagation Measurements and Models

abstract: Comparisons of outdoor Urban Microcell (UMi) large-scale path loss models, root mean square (RMS) delay spreads (DS), angular spreads (AS), and the number of spatial beams for extensive measurements performed at 28, 38, 73, and 142 GHz are presented in this letter. Measurement campaigns were conducted from 2011-2020 in downtown Austin, Texas, Manhattan (New York City), and Brooklyn, New York with communication ranges up to 930 m. Key similarities and differences in outdoor wireless channels are observed when comparing the channel statistics across a wide range of frequencies from millimeter-wave to sub-THz bands. Path loss exponents (PLEs) are remarkably similar over all measured frequencies, when referenced to the first meter free space path loss, and the RMS DS and AS decrease as frequency increases. The similar PLEs from millimeter-wave to THz frequencies imply that spacing between cellular base stations will not have to change as carrier frequencies increase towards THz, since wider bandwidth channels at sub-THz or THz carrier frequencies will cover similar distances because antenna gains increase quadratically with increasing frequency when the physical antenna area remain constant.


title: In-situ non-equilibrium nanomechanics in a proton-conducting ceramic at low temperatures

abstract: Nanostructured proton-conducting ceramics (PCCs) have attracted considerable interest as moderate-temperature proton conductors. Structure dynamics during proton conduction, particularly at grain boundaries, are crucial for stability and proton transport in nanostructured PCCs. A common assumption is that PCCs are structurally stable at low operating temperatures; however, material polycrystallinity, absorption, and reactive operating conditions have so far prevented verifying this assumption by nano resolved in-situ structure measurements. Here, in an archetypal PCC BaZr0.8Y0.2O3-d the premise of structural stability is demonstrated to be inaccurate at temperatures as low as 200 {\deg}C. Coherent X-ray diffraction on a nanostructured BaZr0.8Y0.2O3-d sintered pellet is adapted to image in-situ three-dimensional structural processes inside the constituent submicron grains in a humid nitrogen atmosphere at 200 {\deg}C. Direct observation reveals non-equilibrium defect generation and subsequent grain cracking on a timescale of hours, forming new, otherwise energetically unfavorable facets in BaZr0.8Y0.2O3-d. Furthermore, the structural rearrangements correlate with dynamic inhomogeneities of the lattice constant within the grains, showing potential heterogeneous H+ transport. Our results elucidate the mechanisms behind PCCs structural degradation, overturn existing assumptions about the structure dynamics in PCCs, and fill a method gap for further in-depth in-situ studies of the PCC nanostructure.


title: MPLAPACK version 2.0.1 user manual

abstract: The MPLAPACK (formerly MPACK) is a multiple-precision version of LAPACK (https://www.netlib.org/lapack/). MPLAPACK version 2.0.1 is based on LAPACK version 3.9.1 and translated from Fortran 90 to C++ using FABLE, a Fortran to C++ source-to-source conversion tool (https://github.com/cctbx/cctbx_project/tree/master/fable/). MPLAPACK version 2.0.1 provides the real and complex version of MPBLAS, and the real and complex versions of MPLAPACK support all LAPACK features: solvers for systems of simultaneous linear equations, least-squares solutions of linear systems of equations, eigenvalue problems, and singular value problems, and related matrix factorizations except for mixed-precision routines. The MPLAPACK defines an API for numerical linear algebra, similar to LAPACK. It is easy to port legacy C/C++ numerical codes using MPLAPACK. MPLAPACK supports binary64, binary128, FP80 (extended double), MPFR, GMP, and QD libraries (double-double and quad-double). Users can choose MPFR or GMP for arbitrary accurate calculations, double-double or quad-double for fast 32 or 64-decimal calculations. We can consider the binary64 version as the C++ version of LAPACK. Moreover, it comes with an OpenMP accelerated version of MPBLAS for some routines and CUDA (A100 and V100 support) for double-double versions of Rgemm and Rsyrk. The peak performances of the OpenMP version are almost proportional to the number of cores, and the performances of the CUDA version are impressive, and approximately 400-600 GFlops. MPLAPACK is available at GitHub (https://github.com/nakatamaho/mplapack/) under the 2-clause BSD license.


title: CRANE: a 10 Degree-of-Freedom, Tele-surgical System for Dexterous Manipulation within Imaging Bores

abstract: Physicians perform minimally invasive percutaneous procedures under Computed Tomography (CT) image guidance both for the diagnosis and treatment of numerous diseases. For these procedures performed within Computed Tomography Scanners, robots can enable physicians to more accurately target sub-dermal lesions while increasing safety. However, existing robots for this application have limited dexterity, workspace, or accuracy. This paper describes the design, manufacture, and performance of a highly dexterous, low-profile, 8+2 Degree-ofFreedom (DoF) robotic arm for CT guided percutaneous needle biopsy. In this article, we propose CRANE: CT Robot and Needle Emplacer. The design focuses on system dexterity with high accuracy: extending physicians' ability to manipulate and insert needles within the scanner bore while providing the high accuracy possible with a robot. We also propose and validate a system architecture and control scheme for low profile and highly accurate image-guided robotics, that meets the clinical requirements for target accuracy during an in-situ evaluation. The accuracy is additionally evaluated through a trajectory tracking evaluation resulting in <0.2mm and <0.71degree tracking error. Finally, we present a novel needle driving and grasping mechanism with controlling electronics that provides simple manufacturing, sterilization, and adaptability to accommodate different sizes and types of needles.


title: Multi-Agent Recurrent Rendezvous Using Drive-Based Motivation

abstract: Recent papers have introduced the Motivation Dynamics framework, which uses bifurcations to encode decision-making behavior in an autonomous mobile agent. In this paper, we consider the multi-agent extension of the Motivation Dynamics framework and show how the framework can be extended to encode persistent multi-agent rendezvous behaviors. We analytically characterize the bifurcation properties of the resulting system, and numerically show that it exhibits complex recurrent behavior suggestive of a strange attractor.


title: Developing the Future of Gamma-ray Astrophysics with Monolithic Silicon Pixels

abstract: This paper explores the potential of AstroPix, a project to develop Complementary Metal Oxide Semiconductor (CMOS) pixels for the next generation of space-based high-energy astrophysics experiments. Multimessenger astrophysics is a rapidly developing field whose upcoming missions need support from new detector technology such as AstroPix. ATLASPix, a monolithic silicon detector optimized for the ATLAS particle detector at CERN, is the beginning of the larger AstroPix project. Energy resolution is a driving parameter in the reconstruction of gamma-ray events, and therefore the characterization of ATLASPix energy resolution is the focus of this paper. The intrinsic energy resolution of the detector exceeded our baseline requirements of <10% at 60 keV. The digital output of ATLASPix results in energy resolutions insufficient to advance gamma-ray astronomy. However, the results from the intrinsic energy resolution indicate the digital capability of the detector can be redesigned, and the next generation of pixels for the larger AstroPix project have already been constructed. Iterations of AstroPix-type pixels are an exciting technology candidate to support new space-based missions.


title: KITTI-360: A Novel Dataset and Benchmarks for Urban Scene Understanding in 2D and 3D

abstract: For the last few decades, several major subfields of artificial intelligence including computer vision, graphics, and robotics have progressed largely independently from each other. Recently, however, the community has realized that progress towards robust intelligent systems such as self-driving cars requires a concerted effort across the different fields. This motivated us to develop KITTI-360, successor of the popular KITTI dataset. KITTI-360 is a suburban driving dataset which comprises richer input modalities, comprehensive semantic instance annotations and accurate localization to facilitate research at the intersection of vision, graphics and robotics. For efficient annotation, we created a tool to label 3D scenes with bounding primitives and developed a model that transfers this information into the 2D image domain, resulting in over 150k images and 1B 3D points with coherent semantic instance annotations across 2D and 3D. Moreover, we established benchmarks and baselines for several tasks relevant to mobile perception, encompassing problems from computer vision, graphics, and robotics on the same dataset, e.g., semantic scene understanding, novel view synthesis and semantic SLAM. KITTI-360 will enable progress at the intersection of these research areas and thus contribute towards solving one of today's grand challenges: the development of fully autonomous self-driving systems.


title: Abelian Chern-Simons Gauge Theory on The Lattice

abstract: The Abelian Chern-Simons gauge theory is constructed on the three-dimensional spacetime lattice. This proposal introduces both lattice and dual lattice, and the gauge field on the dual lattice is expressed in terms of the gauge field on the original lattice. This treatment circumvents the issue of forward/backward difference, which is the common problem that many previous proposals have, and also avoids the duplication problem, which prevents people from introducing the dual lattice. The form of the lattice action is very simple, and is symmetric with respect to the three spacetime dimensions. These features make it straightforward to calculate the expectation values of Wilson loops, and the results agree with the topological field theory in continuous spacetime. Generalizations to multiple types of lattices are also discussed.


title: Discriminative Attribution from Counterfactuals

abstract: We present a method for neural network interpretability by combining feature attribution with counterfactual explanations to generate attribution maps that highlight the most discriminative features between pairs of classes. We show that this method can be used to quantitatively evaluate the performance of feature attribution methods in an objective manner, thus preventing potential observer bias. We evaluate the proposed method on three diverse datasets, including a challenging artificial dataset and real-world biological data. We show quantitatively and qualitatively that the highlighted features are substantially more discriminative than those extracted using conventional attribution methods and argue that this type of explanation is better suited for understanding fine grained class differences as learned by a deep neural network.


title: Extended Irreducible Binary Sextic Goppa codes

abstract: Let $n (>3)$ be a prime number and $\Bbb F_{2^n}$ a finite field of $2^n$ elements. Let $L =\Bbb F_{2^n}\cup \{\infty\}$ be the support set and $g(x)$ an irreducible polynomial of degree $6$ over $\Bbb F_{2^n}$. In this paper, we obtain an upper bound on the number of extended irreducible binary Goppa codes $\Gamma(L, g)$ of degree $6$ and length $2^n+1$.


title: Targetless Extrinsic Calibration of Stereo Cameras, Thermal Cameras, and Laser Sensors in the Wild

abstract: The fusion of multi-modal sensors has become increasingly popular in autonomous driving and intelligent robots since it can provide richer information than any single sensor, enhance reliability in complex environments. Multi-sensor extrinsic calibration is one of the key factors of sensor fusion. However, such calibration is difficult due to the variety of sensor modalities and the requirement of calibration targets and human labor. In this paper, we demonstrate a new targetless cross-modal calibration framework by focusing on the extrinsic transformations among stereo cameras, thermal cameras, and laser sensors. Specifically, the calibration between stereo and laser is conducted in 3D space by minimizing the registration error, while the thermal extrinsic to the other two sensors is estimated by optimizing the alignment of the edge features. Our method requires no dedicated targets and performs the multi-sensor calibration in a single shot without human interaction. Experimental results show that the calibration framework is accurate and applicable in general scenes.


title: Safety-Critical Control Synthesis for Unknown Sampled-Data Systems via Control Barrier Functions

abstract: As the complexity of control systems increases, safety becomes an increasingly important property since safety violations can damage the plant and put the system operator in danger. When the system dynamics are unknown, safety-critical synthesis becomes more challenging. Additionally, modern systems are controlled digitally and hence behave as sampled-data systems, i.e., the system dynamics evolve continuously while the control input is applied at discrete time steps. In this paper, we study the problem of control synthesis for safety-critical sampled-data systems with unknown dynamics. We overcome the challenges introduced by sampled-data implementation and unknown dynamics by constructing a set of control barrier function (CBF)-based constraints. By satisfying the constructed CBF constraint at each sampling time, we guarantee the unknown sampled-data system is safe for all time. We formulate a non-convex program to solve for the control signal at each sampling time. We decompose the non-convex program into two convex sub-problems. We illustrate the proposed approach using a numerical case study.


title: LBV phenomenon and binarity: The environment of HR Car

abstract: Luminous blue variable stars (LBVs) are of great interest in massive-star evolution as they experience very high mass-loss episodes within short periods of time. HR Car is a famous member of this class in the Galaxy. It has a large circumstellar nebula and has also been confirmed as being in a binary system. One means of gaining information about the evolutionary status and physical nature of LBVs is studying their environments. We investigated the stellar content within ~100 pc of HR Car and also its circumstellar nebula. Very Large Telescope (VLT) Multi Unit Spectroscopic Explorer (MUSE) observations of a 2'x2' region around the star highlight the incompleteness of stellar classification for stars with magnitudes of V > 13 mag. Eight B0 to B9 stars have been identified which may lie in close spatial vicinity to HR Car. For a region with a radius of r =1.2 degree (~100 pc at a distance of 4.8 kpc) around HR Car, existing catalogs list several late O-type and early B-type stars, but only one early O-type star. Given the relatively low stellar and nebular masses in the HR Car system, no early O-type stars and only a few late O-type stars would be expected in association with HR Car. Instead, HR Car's location in a point vector diagram suggests that HR Car is not isolated, but is part of a moving group with a population of B-type stars in a spiral arm, and it has not received a strong kick from a supernova explosion of a companion star or a merger event. Potential binary evolution pathways for the HR Car system cannot be fully explored because of the unknown nature of the companion star. Furthermore, the MUSE observations reveal the presence of a fast outflow and "bullets" that have been ejected at intervals of about 400 years. These features may have been caused by recurrent mass transfer in the system.


title: Interactive Dynamic Walking: Learning Gait Switching Policies with Generalization Guarantees

abstract: In this paper, we consider the problem of adapting a dynamically walking bipedal robot to follow a leading co-worker while engaging in tasks that require physical interaction. Our approach relies on switching among a family of Dynamic Movement Primitives (DMPs) as governed by a supervisor. We train the supervisor to orchestrate the switching among the DMPs in order to adapt to the leader's intentions, which are only implicitly available in the form of interaction forces. The primary contribution of our approach is its ability to furnish certificates of generalization to novel leader intentions for the trained supervisor. This is achieved by leveraging the Probably Approximately Correct (PAC)-Bayes bounds from generalization theory. We demonstrate the efficacy of our approach by training a neural-network supervisor to adapt the gait of a dynamically walking biped to a leading collaborator whose intended trajectory is not known explicitly.


title: Diversity of growth rates maximizes phytoplankton productivity in an eddying ocean

abstract: In the subtropical gyres, phytoplankton rely on eddies for transporting nutrients from depth to the euphotic zone. But, what controls the rate of nutrient supply for new production? We show that vertical nutrient flux depends both on the vertical motion within the eddying flow and varies nonlinearly with the growth rate of the phytoplankton itself. Flux is maximized when the growth rate matches the inverse of the decorrelation timescale for vertical motion. Using a three-dimensional ocean model and a linear nutrient uptake model, we find that phytoplankton productivity is maximized for a growth rate of 1/3 day$^{-1}$, which corresponds to the timescale of submesoscale dynamics. Variability in the frequency of vertical motion across different physical features of the flow favors phytoplankton production with different growth rates. Such a growth-transport feedback can generate diversity in the phytoplankton community structure at submesoscales and higher net productivity in the presence of community diversity.


title: The Role of Lookahead and Approximate Policy Evaluation in Reinforcement Learning with Linear Value Function Approximation

abstract: Function approximation is widely used in reinforcement learning to handle the computational difficulties associated with very large state spaces. However, function approximation introduces errors which may lead to instabilities when using approximate dynamic programming techniques to obtain the optimal policy. Therefore, techniques such as lookahead for policy improvement and m-step rollout for policy evaluation are used in practice to improve the performance of approximate dynamic programming with function approximation. We quantitatively characterize, for the first time, the impact of lookahead and m-step rollout on the performance of approximate dynamic programming (DP) with function approximation: (i) without a sufficient combination of lookahead and m-step rollout, approximate DP may not converge, (ii) both lookahead and m-step rollout improve the convergence rate of approximate DP, and (iii) lookahead helps mitigate the effect of function approximation and the discount factor on the asymptotic performance of the algorithm. Our results are presented for two approximate DP methods: one which uses least-squares regression to perform function approximation and another which performs several steps of gradient descent of the least-squares objective in each iteration.


title: Evaluation of Deep Neural Network Domain Adaptation Techniques for Image Recognition

abstract: It has been well proved that deep networks are efficient at extracting features from a given (source) labeled dataset. However, it is not always the case that they can generalize well to other (target) datasets which very often have a different underlying distribution. In this report, we evaluate four different domain adaptation techniques for image classification tasks: DeepCORAL, DeepDomainConfusion, CDAN and CDAN+E. These techniques are unsupervised given that the target dataset dopes not carry any labels during training phase. We evaluate model performance on the office-31 dataset. A link to the github repository of this report can be found here: https://github.com/agrija9/Deep-Unsupervised-Domain-Adaptation.


title: An Open Problem on the Bentness of Mesnager's Functions

abstract: Let $n=2m$. In the present paper, we study the binomial Boolean functions of the form $$f_{a,b}(x) = \mathrm{Tr}_1^{n}(a x^{2^m-1 }) +\mathrm{Tr}_1^{2}(bx^{\frac{2^n-1}{3} }), $$ where $m$ is an even positive integer, $a\in \mathbb{F}_{2^n}^*$ and $b\in \mathbb{F}_4^*$. We show that $ f_{a,b}$ is a bent function if the Kloosterman sum $$K_{m}\left(a^{2^m+1}\right)=1+ \sum_{x\in \mathbb{F}_{2^m}^*} (-1)^{\mathrm{Tr}_1^{m}(a^{2^m+1} x+ \frac{1}{x})}$$ equals $4$, thus settling an open problem of Mesnager. The proof employs tools including computing Walsh coefficients of Boolean functions via multiplicative characters, divisibility properties of Gauss sums, and graph theory.


title: Hat guessing number and guaranteed subgraphs

abstract: The hat guessing number of a graph is a parameter related to the hat guessing game for graphs introduced by Winkler. In this paper, we show that graphs of sufficiently large hat guessing number must contain arbitrary trees and arbitrarily long cycles as subgraphs. More precisely, for each tree $T$, there exists a value $N = N(T)$ such that every graph that does not contain $T$ as a subgraph has hat guessing number at most $N$, and for each integer $c$, there exists a value $N' = N'(c)$ such that every graph with no cycle of length greater than $c$ has hat guessing number at most $N'$.


title: Weakly Supervised Keypoint Discovery

abstract: In this paper, we propose a method for keypoint discovery from a 2D image using image-level supervision. Recent works on unsupervised keypoint discovery reliably discover keypoints of aligned instances. However, when the target instances have high viewpoint or appearance variation, the discovered keypoints do not match the semantic correspondences over different images. Our work aims to discover keypoints even when the target instances have high viewpoint and appearance variation by using image-level supervision. Motivated by the weakly-supervised learning approach, our method exploits image-level supervision to identify discriminative parts and infer the viewpoint of the target instance. To discover diverse parts, we adopt a conditional image generation approach using a pair of images with structural deformation. Finally, we enforce a viewpoint-based equivariance constraint using the keypoints from the image-level supervision to resolve the spatial correlation problem that consistently appears in the images taken from various viewpoints. Our approach achieves state-of-the-art performance for the task of keypoint estimation on the limited supervision scenarios. Furthermore, the discovered keypoints are directly applicable to downstream tasks without requiring any keypoint labels.


title: Escape from parsimony of a double-cut-and-join genome evolution process

abstract: We analyze models of genome evolution based on both restricted and unrestricted double-cut-and-join (DCJ) operations. We compare the number of operations along the evolutionary trajectory to the DCJ distance of the genome from its ancestor at each step, and determine at what point they diverge: the process escapes from parsimony. Adapting the method developed by Berestycki and Durret (2006), we estimate the number of cycles in the breakpoint graph of a random genome at time $t$ and its ancestral genome by the number of tree components of an Erd\"os-R\'enyi random graph constructed from the model of evolution. In both models, the process on a genome of size $n$ is bound to its parsimonious estimate up to $t\approx n/2$ steps.


title: The JHU submission to VoxSRC-21: Track 3

abstract: This technical report describes Johns Hopkins University speaker recognition system submitted to Voxceleb Speaker Recognition Challenge 2021 Track 3: Self-supervised speaker verification (closed). Our overall training process is similar to the proposed one from the first place team in the last year's VoxSRC2020 challenge. The main difference is a recently proposed non-contrastive self-supervised method in computer vision (CV), distillation with no labels (DINO), is used to train our initial model, which outperformed the last year's contrastive learning based on momentum contrast (MoCo). Also, this requires only a few iterations in the iterative clustering stage, where pseudo labels for supervised embedding learning are updated based on the clusters of the embeddings generated from a model that is continually fine-tuned over iterations. In the final stage, Res2Net50 is trained on the final pseudo labels from the iterative clustering stage. This is our best submitted model to the challenge, showing 1.89, 6.50, and 6.89 in EER(%) in voxceleb1 test o, VoxSRC-21 validation, and test trials, respectively.


title: Comment on "Accumulating Evidence for the Associate Production of a Neutral Scalar with Mass around 151 GeV"

abstract: A recent paper [2109.02650] accumulates evidence for a new fundamental particle by combining several CMS and ATLAS searches for the Standard Model Higgs boson. The putative particle is a neutral scalar, $S$, with a mass of about 151 GeV. The reported significances are $5.1\sigma$ local and $4.8\sigma$ global. This nearly reaches the $5\sigma$ threshold for a discovery in high-energy physics. In this brief note we cast doubt on the strength of the evidence for a new particle. After taking into account the fact that signals were fitted to six different channels, we find that the significances are only $4.1\sigma$ local and $3.5\sigma$ global. The code and instructions for reproducing our calculations are available at https://github.com/andrewfowlie/accumulating_evidence.


title: Anyons in Quantum Hall Interferometry

abstract: The quantum Hall (QH) effect represents a unique playground where quantum coherence of electrons can be exploited for various applications, from metrology to quantum computation. In the fractional regime it also hosts anyons, emergent quasiparticles that are neither bosons nor fermions and possess fractional statistics. Their detection and manipulation represent key milestones in view of topologically protected quantum computation schemes. Exploiting the high degree of phase coherence, edge states in the QH regime have been investigated by designing and constructing electronic interferometers, able to reveal the coherence and statistical properties of the interfering constituents. Here, we review the two main geometries developed in the QH regime, the Mach-Zehnder and the Fabry-Perot interferometers. We present their basic working principles, fabrication methods, and the main results obtained both in the integer and fractional QH regime. We will also show how recent technological advances led to the direct experimental demonstration of fractional statistics for Laughlin quasiparticles in a Fabry-Perot interferometric setup.


title: Fine-grained Solar Flare Forecasting Based on the Hybrid Convolutional Neural Networks

abstract: Improving the performance of solar flare forecasting is a hot topic in solar physics research field. Deep learning has been considered a promising approach to perform solar flare forecasting in recent years. We first used the Generative Adversarial Networks (GAN) technique augmenting sample data to balance samples with different flare classes. We then proposed a hybrid convolutional neural network (CNN) model M for forecasting flare eruption in a solar cycle. Based on this model, we further investigated the effects of the rising and declining phases for flare forecasting. Two CNN models, i.e., Mrp and Mdp, were presented to forecast solar flare eruptions in the rising phase and declining phase of solar cycle 24, respectively. A series of testing results proved: 1) Sample balance is critical for the stability of the CNN model. The augmented data generated by GAN effectively improved the stability of the forecast model. 2) For C-class, M-class, and X-class flare forecasting using Solar Dynamics Observatory (SDO) line-of-sight (LOS) magnetograms, the means of true skill statistics (TSS) score of M are 0.646, 0.653 and 0.762, which improved by 20.1%, 22.3%, 38.0% compared with previous studies. 3) It is valuable to separately model the flare forecasts in the rising and declining phases of a solar cycle. Compared with model M, the means of TSS score for No-flare, C-class, M-class, X-class flare forecasting of the Mrp improved by 5.9%, 9.4%, 17.9% and 13.1%, and the Mdp improved by 1.5%, 2.6%, 11.5% and 12.2%.


title: A Mattila-Sj\"{o}lin theorem for triangles

abstract: We show for a compact set $E \subset \mathbb{R}^d$, $d \geq 4$, that if the Hausdorff dimension of $E$ is larger than $\frac{2}{3}d+1$, then the set of congruence classes of triangles formed by triples of points of $E$ has nonempty interior. Here we understand the set of congruence classes of triangles formed by triples of points of $E$ as the set $$\Delta_{\text{tri}}(E) = \left \{ (t,r, \alpha) : |x-z|=t, |y-z|=r \, \text{ and }\, \alpha= \alpha(x,z,y), \ x,y,z \in E \right \},$$ where $\alpha (x,z,y)$ denotes the angle formed by $x$, $y$ and $z$ , centered at $z$. This extends the Mattila-Sj\"{o}lin theorem that establishes a non-empty interior for the distance set instead of the set of congruence classes of triangles. These theorems can be thought of as refinements and extensions of the statements in the well known Falconer distance problem.


title: SYGMA: System for Generalizable Modular Question Answering OverKnowledge Bases

abstract: Knowledge Base Question Answering (KBQA) tasks that in-volve complex reasoning are emerging as an important re-search direction. However, most KBQA systems struggle withgeneralizability, particularly on two dimensions: (a) acrossmultiple reasoning types where both datasets and systems haveprimarily focused on multi-hop reasoning, and (b) across mul-tiple knowledge bases, where KBQA approaches are specif-ically tuned to a single knowledge base. In this paper, wepresent SYGMA, a modular approach facilitating general-izability across multiple knowledge bases and multiple rea-soning types. Specifically, SYGMA contains three high levelmodules: 1) KB-agnostic question understanding module thatis common across KBs 2) Rules to support additional reason-ing types and 3) KB-specific question mapping and answeringmodule to address the KB-specific aspects of the answer ex-traction. We demonstrate effectiveness of our system by evalu-ating on datasets belonging to two distinct knowledge bases,DBpedia and Wikidata. In addition, to demonstrate extensi-bility to additional reasoning types we evaluate on multi-hopreasoning datasets and a new Temporal KBQA benchmarkdataset on Wikidata, namedTempQA-WD1, introduced in thispaper. We show that our generalizable approach has bettercompetetive performance on multiple datasets on DBpediaand Wikidata that requires both multi-hop and temporal rea-soning


title: Arranging the order of passengers on the boarding bridge to reduce the boarding time for single-aisle aircraft

abstract: Reducing the aircraft boarding time is a common problem not only for airlines, but also for passengers and airports. Group boarding is a popular boarding strategy that separates the passengers into several groups and those groups, which are then called in a certain order. Group boarding can reduce the boarding time compared with that in random order boarding; however, it is insufficient in several real scenarios because the passengers are not separated strictly into groups. In this paper, we propose a boarding strategy that arranges the order of the boarding passengers at the boarding gate. Although this approach appears more time-consuming, we show that such a rearrangement can be applied efficiently to the waiting queue of a single-aisle aircraft. We quantitatively demonstrate the boarding times for various patterns of this approach and discuss the mechanism underlying the reduction of boarding time. This strategy is a promising approach to reduce boarding times and can replace the conventional boarding strategy.


title: Warp-Refine Propagation: Semi-Supervised Auto-labeling via Cycle-consistency

abstract: Deep learning models for semantic segmentation rely on expensive, large-scale, manually annotated datasets. Labelling is a tedious process that can take hours per image. Automatically annotating video sequences by propagating sparsely labeled frames through time is a more scalable alternative. In this work, we propose a novel label propagation method, termed Warp-Refine Propagation, that combines semantic cues with geometric cues to efficiently auto-label videos. Our method learns to refine geometrically-warped labels and infuse them with learned semantic priors in a semi-supervised setting by leveraging cycle consistency across time. We quantitatively show that our method improves label-propagation by a noteworthy margin of 13.1 mIoU on the ApolloScape dataset. Furthermore, by training with the auto-labelled frames, we achieve competitive results on three semantic-segmentation benchmarks, improving the state-of-the-art by a large margin of 1.8 and 3.61 mIoU on NYU-V2 and KITTI, while matching the current best results on Cityscapes.


title: Simulating the nematic-isotropic phase transition of liquid crystal model via generalized replica-exchange method

abstract: The nematic-isotropic (NI) phase transition of 4-cyano-4'-pentylbiphenyl (5CB) was simulated using the generalized replica-exchange method (gREM) based on molecular dynamics simulations. The effective temperature is introduced in gREM, allowing the enhanced sampling of configurations in the unstable region, which is intrinsic to the first-order phase transition. The sampling performance was analyzed with different system sizes and compared with that of the temperature replica-exchange method (tREM). It was observed that gREM is capable of sampling configurations at sufficient replica-exchange acceptance ratios even around the NI transition temperature. A bimodal distribution of the order parameter at the transition region was found, which is in agreement with the mean-field theory. In contrast, tREM is ineffective around the transition temperature owing to the potential energy gap between the nematic and isotropic phases.


title: Bound state in the continuum in slab waveguide enables low-threshold quantum-dot lasing

abstract: Colloidal quantum dots (CQDs) are a promising gain material for solution-processed, wavelength-tunable lasers, with potential application in displays, communications, and biomedical devices. In this work, we combine a CQD film with an array of titanium dioxide (TiO2) nanoantennas to achieve lasing via bound states in the continuum (BICs), which are symmetry-protected cavity modes with giant quality factors. Here, the BICs arise from slab waveguide modes in the planar film, coupled to the periodic nanoantenna array. We engineer the thickness of the CQD film and size of the nanoantennas to achieve a BIC with good spatial and spectral overlap with the CQDs, based on a 2nd-order TE-polarized waveguide mode. We obtain room-temperature lasing with a low threshold of approximately 11 kW/cm2 (peak intensity) under 5 ns-pulsed optical excitation. This work sheds light on the optical modes in solution-processed, distributed-feedback lasers, and highlights BICs as effective, versatile, surface-emitting lasing modes.


title: Rate of the enhanced dissipation for the two-jet Kolmogorov type flow on the unit sphere

abstract: We study the enhanced dissipation for the two-jet Kolmogorov type flow which is a stationary solution to the Navier-Stokes equations on the two-dimensional unit sphere given by the zonal spherical harmonic function of degree two. Based on the pseudospectral bound method developed by Ibrahim, Maekawa, and Masmoudi [15] and a modified version of the Gearhart-Pr\"uss type theorem shown by Wei [48], we derive an estimate for the resolvent of the linearized operator along the imaginary axis and show that a solution to the linearized equation rapidly decays at the rate $O(e^{-\sqrt{\nu}\,t})$ when the viscosity coefficient $\nu$ is sufficiently small as in the case of the plane Kolmogorov flow.


title: Robust Binary Linear Programming Under Implementation Uncertainty

abstract: This paper studies binary linear programming problems in the presence of uncertainties that may cause solution values to change during implementation. This type of uncertainty, termed implementation uncertainty, is modeled explicitly affecting the decision variables rather than model parameters. The binary nature of the decision variables invalidates the use of the existing models for this type of uncertainty. The robust solutions obtained are optimal for a worst-case min-max objective and allow a controlled degree of infeasibility with respect to the associated deterministic problem. Structural properties are used to reformulate the problem as a mixed-integer linear binary program. The degree of solution conservatism is controlled by combining both constraint relaxation and cardinality-constrained parameters. Solutions for optimization problems under implementation uncertainty consist of a set of robust solutions; the selection of solutions from this possibly large set is formulated as an optimization problem over the robust set. Results from an experimental study in the context of the knapsack problem suggest the methodology yields solutions that perform well in terms of objective value and feasibility. Furthermore, the selection approach can identify robust solutions that possess desirable implementation characteristics.


title: Ultra-High-Energy Cosmic Rays and Neutrinos from Relativistic Jets of Active Galactic Nuclei

abstract: In Mbarek & Caprioli (2019), we laid the groundwork for studying the espresso paradigm Caprioli (2015), a reacceleration mechanism to boost galactic CRs to UHECR levels. Our bottom-up approach uses realistic 3D MHD simulations of relativistic AGN jets and accounts for all of the crucial ingredients of a universal acceleration theory: injection, acceleration, and escape in realistic environments. Our results are consistent with the main features of UHECR spectra, i.e., power-law slopes, chemical composition, and anisotropy. In Mbarek & Caprioli (2021), we refine our model by including sub-grid particle scattering to model small-scale magnetic turbulence that cannot be resolved by MHD simulations, constraining for the first time one crucial but hard-to-model ingredient, and allowing us to establish the relative importance of espresso and stochastic shear acceleration in relativistic jets. Our framework also enables us to analyze high-energy neutrinos produced from our accelerated UHECRs considering the effects of external photon fields, and to incorporate nucleus photodisintegration. The spectra we obtain are consistent with the picture drawn by observations with Auger, Telescope Array, and IceCube observatory.


title: Not Only Domain Randomization: Universal Policy with Embedding System Identification

abstract: Domain randomization (DR) cannot provide optimal policies for adapting the learning agent to the dynamics of the environment, although it can generalize sub-optimal policies to work in a transferred domain. In this paper, we present Universal Policy with Embedding System Identification (UPESI) as an implicit system identification (SI) approach with universal policies (UPs), as a learning-based control method to execute optimal actions adaptively in environments with various dynamic properties. Previous approaches of SI for adaptive policies either conduct explicit SI, which is testified to be an ill-posed problem, or suffer from low efficiency without leveraging the simulation oracle. We propose to conduct SI in the embedding space of system dynamics by leveraging a learned forward dynamics model, and use Bayesian optimization for the SI process given transition data in a new environment. The identified embeddings are applied as additional input to the UP to enable its dynamics adaptability. Experiments demonstrate the advantageous performances of our proposed UP with embedding SI over standard DR and conventional SI approaches on both low-dimensional and high-dimensional simulation tasks.


title: Search for photoproduction of axion-like particles at GlueX

abstract: We present a search for axion-like particles, $a$, produced in photon-proton collisions at a center-of-mass energy of approximately 4 GeV, focusing on the scenario where the $a$-gluon coupling is dominant. The search uses $a\to\gamma\gamma$ and $a\to\pi^+\pi^-\pi^0$ decays, and a data sample corresponding to an integrated luminosity of 168 pb$^{-1}$ collected with the GlueX detector. The search for $a\to\gamma\gamma$ decays is performed in the mass range of $180 < m_a < 480$ MeV, while the search for $a\to\pi^+\pi^-\pi^0$ decays explores the $600 < m_a < 720$ MeV region. No evidence for a signal is found, and 90% confidence-level exclusion limits are placed on the $a$-gluon coupling strength. These constraints are the most stringent to date over much of the mass ranges considered.


title: Band Splitting and Long-lived Carrier Recombination in Ferromagnetic CrSiTe3 Nanosheets

abstract: Magnetic layered ternary chalcogenides hold great promise for future spin-optoelectronic devices in the two dimensional limit. Understanding how the properties of the materials are impacted by magnetic ordering and the spin-orbit interactions is critically needed information for the development of applications. Ultrafast transient reflectance (TR) and photocurrent (PC) spectroscopies are combined with ab initio density functional theory (DFT) calculations to investigate the band structure and photoresponse of a layered magnetic semiconductor CrSiTe3 (CST) nanosheet in the paramagnetic (PM, 300K) and ferromagnetic (FM, 10 K) phases. We observe both a decrease of the direct bandgap and emergence of a 120 meV splitting of the optical transition when the FM phase is present. DFT band structure calculations suggest that the band modifications are driven by a FM ordering-induced band splitting between the Te p and the Cr d states at the valence and conduction band edges. We find that the majority of carriers photoexcited at the direct gap recombine within picoseconds through defect-mediated recombination, but that 2-3 % of the electrons scatter into indirect conduction band valleys resulting in very long-lived electrons and holes. Those long-lived carriers contribute to the broadband PC response of CST devices that also features indirect absorption. These results provide critical insights into the dynamics and energy landscape of photoexcited electrons and holes, and how they are impacted by spin-ordering effects in layered ferromagnets.


title: DynG2G: An Efficient Stochastic Graph Embedding Method for Temporal Graphs

abstract: Dynamic graph embedding has gained great attention recently due to its capability of learning low dimensional graph representations for complex temporal graphs with high accuracy. However, recent advances mostly focus on learning node embeddings as deterministic "vectors" for static graphs yet disregarding the key graph temporal dynamics and the evolving uncertainties associated with node embedding in the latent space. In this work, we propose an efficient stochastic dynamic graph embedding method (DynG2G) that applies an inductive feed-forward encoder trained with node triplet-based contrastive loss. Every node per timestamp is encoded as a time-dependent probabilistic multivariate Gaussian distribution in the latent space, hence we can quantify the node embedding uncertainty on-the-fly. We adopted eight different benchmarks that represent diversity in size (from 96 nodes to 87,626 and from 13,398 edges to 4,870,863) and diversity in dynamics. We demonstrate via extensive experiments on these eight dynamic graph benchmarks that DynG2G achieves new state-of-the-art performance in capturing the underlying temporal node embeddings. We also demonstrate that DynG2G can predict the evolving node embedding uncertainty, which plays a crucial role in quantifying the intrinsic dimensionality of the dynamical system over time. We obtain a universal relation of the optimal embedding dimension, $L_o$, versus the effective dimensionality of uncertainty, $D_u$, and we infer that $L_o=D_u$ for all cases. This implies that the uncertainty quantification approach we employ in the DynG2G correctly captures the intrinsic dimensionality of the dynamics of such evolving graphs despite the diverse nature and composition of the graphs at each timestamp. Moreover, this $L_0 - D_u$ correlation provides a clear path to select adaptively the optimum embedding size at each timestamp by setting $L \ge D_u$.


title: An Adaptive Deep Learning Framework for Day-ahead Forecasting of Photovoltaic Power Generation

abstract: Accurate forecasts of photovoltaic power generation (PVPG) are essential to optimize operations between energy supply and demand. Recently, the propagation of sensors and smart meters has produced an enormous volume of data, which supports the development of data based PVPG forecasting. Although emerging deep learning (DL) models, such as the long short-term memory (LSTM) model, based on historical data, have provided effective solutions for PVPG forecasting with great successes, these models utilize offline learning. As a result, DL models cannot take advantage of the opportunity to learn from newly-arrived data, and are unable to handle concept drift caused by installing extra PV units and unforeseen PV unit failures. Consequently, to improve day-ahead PVPG forecasting accuracy, as well as eliminate the impacts of concept drift, this paper proposes an adaptive LSTM (AD-LSTM) model, which is a DL framework that can not only acquire general knowledge from historical data, but also dynamically learn specific knowledge from newly-arrived data. A two-phase adaptive learning strategy (TP-ALS) is integrated into AD-LSTM, and a sliding window (SDWIN) algorithm is proposed, to detect concept drift in PV systems. Multiple datasets from PV systems are utilized to assess the feasibility and effectiveness of the proposed approaches. The developed AD-LSTM model demonstrates greater forecasting capability than the offline LSTM model, particularly in the presence of concept drift. Additionally, the proposed AD-LSTM model also achieves superior performance in terms of day-ahead PVPG forecasting compared to other traditional machine learning models and statistical models in the literature.

categories: cs.LG cs.AI physics.comp-ph physics.data-an


title: Higher form Yang-Mills as higher BFYM theories

abstract: The YM theory has been generalized to 2YM and 3YM theories. Similarly, we generalize the BFYM theory to "2BFYM" and "3BFYM" theories. Then, we show that these higher BFYM theories can give the formulations of the corresponding higher form YM theories. Finally, we study the gauge symmetries of these higher BFYM theories.


title: What do the two times in two-time correlation functions mean for interpreting tr-ARPES?

abstract: Time-resolved angle-resolved photoemission spectroscopy is one of the most powerful pump-probe measurements of materials driven far from equilibrium. Unlike the linear-response regime, where the frequency-dependent response function is independent of time, in a far-from-equilibrium experiment, the response function depends on two times in the time domain. In this work, we describe how one can use time-dependent frequency response functions and how they involve contributions from times that are near to each other. This implies that they should not be thought of as a frequency-dependent response at a single definite time. Instead, the Fourier uncertainty relations show that they involve contributions from ranges of times and must be interpreted in this light. We use this insight to help understand what time-resolved photoemission measurements actually measure.


title: Emergent Neural Network Mechanisms for Generalization to Objects in Novel Orientations

abstract: The capability of Deep Neural Networks (DNNs) to recognize objects in orientations outside the distribution of the training data is not well understood. We present evidence that DNNs are capable of generalizing to objects in novel orientations by disseminating orientation-invariance obtained from familiar objects seen from many viewpoints. This capability strengthens when training the DNN with an increasing number of familiar objects, but only in orientations that involve 2D rotations of familiar orientations. We show that this dissemination is achieved via neurons tuned to common features between familiar and unfamiliar objects. These results implicate brain-like neural mechanisms for generalization.


title: Runtime Safety Assurance for Learning-enabled Control of Autonomous Driving Vehicles

abstract: Providing safety guarantees for Autonomous Vehicle (AV) systems with machine-learning-based controllers remains a challenging issue. In this work, we propose Simplex-Drive, a framework that can achieve runtime safety assurance for machine-learning enabled controllers of AVs. The proposed Simplex-Drive consists of an unverified Deep Reinforcement Learning (DRL)-based advanced controller (AC) that achieves desirable performance in complex scenarios, a Velocity-Obstacle (VO) based baseline safe controller (BC) with provably safety guarantees, and a verified mode management unit that monitors the operation status and switches the control authority between AC and BC based on safety-related conditions. We provide a formal correctness proof of Simplex-Drive and conduct a lane-changing case study in dense traffic scenarios. The simulation experiment results demonstrate that Simplex-Drive can always ensure operation safety without sacrificing control performance, even if the DRL policy may lead to deviations from the safe status.


title: Scaled Affine Quantization of Ultralocal $\varphi^4_2$ a comparative Path Integral Monte Carlo study with Canonical Quantization

abstract: After the success of affine quantization in proving through Monte Carlo analysis that the covariant euclidean scalar field theory, $\varphi^r_n$, where $r$ denotes the power of the interaction term and $n = s + 1$ with $s$ the spatial dimension and $1$ adds imaginary time, such that $r \geq 2n/(n-2)$ can be acceptably quantized and the resulting theory is nontrivial, unlike what happens using canonical quantization, we show here that the same has to be expected for $r>2$ and any $n$ even for the ultralocal field theory. In particular we consider the ultralocal $\varphi^4_2$ model and study its renormalized properties for both the scaled canonical quantization version and the scaled affine quantization version through path integral Monte Carlo.


title: Lithium-ion Battery State of Health Estimation based on Cycle Synchronization using Dynamic Time Warping

abstract: The state of health (SOH) estimation plays an essential role in battery-powered applications to avoid unexpected breakdowns due to battery capacity fading. However, few studies have paid attention to the problem of uneven length of degrading cycles, simply employing manual operation or leaving to the automatic processing mechanism of advanced machine learning models, like long short-term memory (LSTM). As a result, this causes information loss and caps the full capability of the data-driven SOH estimation models. To address this challenge, this paper proposes an innovative cycle synchronization way to change the existing coordinate system using dynamic time warping, not only enabling the equal length inputs of the estimation model but also preserving all information. By exploiting the time information of the time series, the proposed method embeds the time index and the original measurements into a novel indicator to reflect the battery degradation status, which could have the same length over cycles. Adopting the LSTM as the basic estimation model, the cycle synchronization-based SOH model could significantly improve the prediction accuracy by more than 30% compared to the traditional LSTM.


title: When in Doubt: Improving Classification Performance with Alternating Normalization

abstract: We introduce Classification with Alternating Normalization (CAN), a non-parametric post-processing step for classification. CAN improves classification accuracy for challenging examples by re-adjusting their predicted class probability distribution using the predicted class distributions of high-confidence validation examples. CAN is easily applicable to any probabilistic classifier, with minimal computation overhead. We analyze the properties of CAN using simulated experiments, and empirically demonstrate its effectiveness across a diverse set of classification tasks.


title: Two-Stage Channel Estimation Approach for Cell-Free IoT With Massive Random Access

abstract: We investigate the activity detection and channel estimation issues for cell-free Internet of Things (IoT) networks with massive random access. In each time slot, only partial devices are active and communicate with neighboring access points (APs) using non-orthogonal random pilot sequences. Different from the centralized processing in cellular networks, the activity detection and channel estimation in cell-free IoT is more challenging due to the distributed and user-centric architecture. We propose a two-stage approach to detect the random activities of devices and estimate their channel states. In the first stage, the activity of each device is jointly detected by its adjacent APs based on the vector approximate message passing (Vector AMP) algorithm. In the second stage, each AP re-estimates the channel using the linear minimum mean square error (LMMSE) method based on the detected activities to improve the channel estimation accuracy. We derive closed-form expressions for the activity detection error probability and the mean-squared channel estimation errors for a typical device. Finally, we analyze the performance of the entire cell-free IoT network in terms of coverage probability. Simulation results validate the derived closed-form expressions and show that the cell-free IoT significantly outperforms the collocated massive MIMO and small-cell schemes in terms of coverage probability.


title: Justification for zeta function regularization

abstract: Using the fact that a finite sum of power series are given by the difference between two zeta functions, we justify the usage of the zeta function with a negative variable in physical problems to avoid the divergence of the infinite sum. We will show that in the case of magnetization of graphene, the zeta function with negative variable arises as a result of cut-off energy between two consecutive Landau levels. Furthermore, similar justification can be applied to the case of zero temperature Casimir force in parallel-plate geometry.


title: CateCom: a practical data-centric approach to categorization of computational models

abstract: The advent of data-driven science in the 21st century brought about the need for well-organized structured data and associated infrastructure able to facilitate the applications of Artificial Intelligence and Machine Learning. We present an effort aimed at organizing the diverse landscape of physics-based and data-driven computational models in order to facilitate the storage of associated information as structured data. We apply object-oriented design concepts and outline the foundations of an open-source collaborative framework that is: (1) capable of uniquely describing the approaches in structured data, (2) flexible enough to cover the majority of widely used models, and (3) utilizes collective intelligence through community contributions. We present example database schemas and corresponding data structures and explain how these are deployed in software at the time of this writing.


title: Scalarized neutron stars in massive scalar-tensor gravity: X-ray pulsars and tidal deformability

abstract: Neutron stars (NSs) in scalar-tensor theories of gravitation with the phenomenon of spontaneous scalarization can develop significant deviations from general relativity. Cases with a massless scalar were studied widely. Here we compare the NS scalarizations in the Damour--Esposito-Far{\`e}se theory, the Mendes-Ortiz theory, and the $\xi$-theory with a massive scalar field. Numerical solutions for slowly rotating NSs are obtained. They are used to construct the X-ray pulse profiles of a pair of extended hot spots on the surface of NSs. We also calculate the tidal deformability for NSs with spontaneous scalarization which is done for the first time with a massive scalar field. We show the universal relation between the moment of inertia and the tidal deformability. The X-ray pulse profiles, the tidal deformability, and the universal relation may help to constrain the massive scalar-tensor theories in X-ray and gravitational-wave observations of NSs, including the Neutron star Interior Composition Explorer (NICER) satellite, Square Kilometre Array (SKA) telescope, and LIGO/Virgo/KAGRA laser interferometers.


title: Scattering of low energy neutrinos and antineutrinos by atomic electrons

abstract: Studies of neutrino mixing and oscillations, solar neutrinos as background in dark matter searches involving electron detection, detection of sterile neutrino warm dark matter, and of possible electromagnetic properties of neutrinos, have generated interest in the low energy O(10 keV) scattering of electron neutrinos and antineutrinos by atomic electrons where the binding of the atomic electron cannot be ignored. Of particular interest is the ionization of atoms by neutrinos and antineutrinos. Most existing calculations are based on modifications of the free electron differential cross section, which destroy the relationship between the neutrino helicities and the orbital and spin angular momenta of the atomic electrons. The present calculations maintain the full collision dynamics by formulating the scattering in configuration space using the Bound Interaction Picture, rather than the usual formulation in the Interaction Picture in momentum space as appropriate to scattering by free electrons. Energy spectra of ionization electrons produced by scattering of neutrinos and antineutrinos with energies of 5, 10, 20, and 30 keV by hydrogen, helium and neon have been calculated using Dirac central field eigenfunctions and are presented as ratios to the spectra for scattering by free electrons. Binding effects increase strongly with atomic number, are largest for low neutrino energy and, for each neutrino energy, greatest at the high electron energy end of the spectrum. The most extreme effects of binding are for 5 keV scattering by Ne where the ratios are less than 0.1. The energy spectra are calculated for both a Coulombic final electron state and a free final electron state. The results indicate that the binding effects from the continuum state of the final electron are significant and can be comparable to those arising from the bound initial electron state.


title: Direct numerical simulations of optimal thermal convection in rotating plane layer dynamos

abstract: The heat transfer behavior of convection-driven dynamos in a rotating plane layer between two parallel plates, heated from below and cooled from the top, is investigated. At a fixed rotation rate (Ekman Number, $E=10^{-6}$) and fluid properties (thermal and magnetic Prandtl numbers, $Pr=Pr_m=1$), both dynamo convection (DC) and non-magnetic rotating convection (RC) simulations are performed to demarcate the effect of magnetic field on heat transport at different thermal forcings (Rayleigh Number, $Ra=3.83\times10^{9}-3.83\times10^{10}$). In this range, our turbulence resolving simulations demonstrate the existence of an optimum thermal forcing, at which heat transfer between the plates in DC exhibits maximum enhancement, as compared to the heat transport in the RC simulations. Unlike any global force balance reported in the literature, the present simulations reveal an increase in the Lorentz force in the \textit{thermal boundary layer}, due to stretching of magnetic field line by the vortices near the walls with no-slip boundary condition. This increase in Lorentz force mitigates turbulence suppression owing to the Coriolis force, resulting in enhanced turbulence and heat transfer


title: SiamEvent: Event-based Object Tracking via Edge-aware Similarity Learning with Siamese Networks

abstract: Event cameras are novel sensors that perceive the per-pixel intensity changes and output asynchronous event streams, showing lots of advantages over traditional cameras, such as high dynamic range (HDR) and no motion blur. It has been shown that events alone can be used for object tracking by motion compensation or prediction. However, existing methods assume that the target always moves and is the stand-alone object. Moreover, they fail to track the stopped non-independent moving objects on fixed scenes. In this paper, we propose a novel event-based object tracking framework, called SiamEvent, using Siamese networks via edge-aware similarity learning. Importantly, to find the part having the most similar edge structure of target, we propose to correlate the embedded events at two timestamps to compute the target edge similarity. The Siamese network enables tracking arbitrary target edge by finding the part with the highest similarity score. This extends the possibility of event-based object tracking applied not only for the independent stand-alone moving objects, but also for various settings of the camera and scenes. In addition, target edge initialization and edge detector are also proposed to prevent SiamEvent from the drifting problem. Lastly, we built an open dataset including various synthetic and real scenes to train and evaluate SiamEvent. Extensive experiments demonstrate that SiamEvent achieves up to 15% tracking performance enhancement than the baselines on the real-world scenes and more robust tracking performance in the challenging HDR and motion blur conditions.


title: On the Geometry of Stable Steiner Tree Instances

abstract: In this note we consider the Steiner tree problem under Bilu-Linial stability. We give strong geometric structural properties that need to be satisfied by stable instances. We then make use of, and strengthen, these geometric properties to show that $1.562$-stable instances of Euclidean Steiner trees are polynomial-time solvable. We also provide a connection between certain approximation algorithms and Bilu-Linial stability for Steiner trees.


title: The 3-adic valuations of Stirling numbers of the first kind

abstract: Let $n$ and $k$ be positive integers. The Stirling number of the first kind, denoted by $s(n,k)$, counts the number of permutations of $n$ elements with $k$ disjoint cycles. Let $p$ be a prime number and denote by $v_p(n)$ the $p$-adic valuation of $n$. In recent years, Lengyel, Komatsu, Young, Leonetti and Sanna made some progress on $v_p(s(n,k))$. Let $a\in\{1,2\}$. In this paper, by using the properties of $m$-th Stirling number of the first kind developed previously and providing a detailed 3-adic analysis, we arrive at an explicit formula on $v_3(s(a3^n, k))$ with $1\le k\le a3^n$. This gives an evidence to a conjecture of Hong and Qiu proposed in 2020. As a corollary, we show that $v_3(s(a3^n,a3^n-k))=2n-1+v_3(k)-v_3(k-1)$ if $k$ is odd and $3\le k\le a3^{n-1}+1$. This supports a conjecture of Lengyel raised in 2015.


title: Multiwavelet-based Operator Learning for Differential Equations

abstract: The solution of a partial differential equation can be obtained by computing the inverse operator map between the input and the solution space. Towards this end, we introduce a \textit{multiwavelet-based neural operator learning scheme} that compresses the associated operator's kernel using fine-grained wavelets. By explicitly embedding the inverse multiwavelet filters, we learn the projection of the kernel onto fixed multiwavelet polynomial bases. The projected kernel is trained at multiple scales derived from using repeated computation of multiwavelet transform. This allows learning the complex dependencies at various scales and results in a resolution-independent scheme. Compare to the prior works, we exploit the fundamental properties of the operator's kernel which enable numerically efficient representation. We perform experiments on the Korteweg-de Vries (KdV) equation, Burgers' equation, Darcy Flow, and Navier-Stokes equation. Compared with the existing neural operator approaches, our model shows significantly higher accuracy and achieves state-of-the-art in a range of datasets. For the time-varying equations, the proposed method exhibits a ($2X-10X$) improvement ($0.0018$ ($0.0033$) relative $L2$ error for Burgers' (KdV) equation). By learning the mappings between function spaces, the proposed method has the ability to find the solution of a high-resolution input after learning from lower-resolution data.


title: Self-Improving Voronoi Construction for a Hidden Mixture of Product Distributions

abstract: We propose a self-improving algorithm for computing Voronoi diagrams under a given convex distance function with constant description complexity. The $n$ input points are drawn from a hidden mixture of product distributions; we are only given an upper bound $m = o(\sqrt{n})$ on the number of distributions in the mixture, and the property that for each distribution, an input instance is drawn from it with a probability of $\Omega(1/n)$. For any $\varepsilon \in (0,1)$, after spending $O\bigl(mn\log^{O(1)} (mn) + m^{\varepsilon} n^{1+\varepsilon}\log(mn)\bigr)$ time in a training phase, our algorithm achieves an $O\bigl(\frac{1}{\varepsilon}n\log m + \frac{1}{\varepsilon}n2^{O(\log^* n)} + \frac{1}{\varepsilon}H\bigr)$ expected running time with probability at least $1 - O(1/n)$, where $H$ is the entropy of the distribution of the Voronoi diagram output. The expectation is taken over the input distribution and the randomized decisions of the algorithm. For the Euclidean metric, the expected running time improves to $O\bigl(\frac{1}{\varepsilon}n\log m + \frac{1}{\varepsilon}H\bigr)$.


title: Flexibility of the factorized form of the unitary coupled cluster ansatz

abstract: The factorized form of the unitary coupled cluster ansatz is a popular state preparation ansatz for electronic structure calculations of molecules on quantum computers. It often is viewed as an approximation (based on the Trotter product formula) for the conventional unitary coupled cluster operator. In this work, we show that the factorized form is quite flexible, allowing one to range from conventional configuration interaction, to conventional unitary coupled cluster, to efficient approximations that lie in between these two. The variational minimization of the energy often allows simpler factorized unitary coupled cluster approximations to achieve high accuracy, even if they do not accurately approximate the Trotter product formula. This is similar to how quantum approximate optimization algorithms can achieve high accuracy with a small number of levels.


title: Beta Decay in Medium-Mass Nuclei with the In-Medium Similarity Renormalization Group

abstract: We review the status of ab initio calculations of allowed beta decays (both Fermi and Gamow-Teller), within the framework of the valence-space in-medium similarity renormalization group approach.


title: Deep Reinforcement Learning with Adjustments

abstract: Deep reinforcement learning (RL) algorithms can learn complex policies to optimize agent operation over time. RL algorithms have shown promising results in solving complicated problems in recent years. However, their application on real-world physical systems remains limited. Despite the advancements in RL algorithms, the industries often prefer traditional control strategies. Traditional methods are simple, computationally efficient and easy to adjust. In this paper, we first propose a new Q-learning algorithm for continuous action space, which can bridge the control and RL algorithms and bring us the best of both worlds. Our method can learn complex policies to achieve long-term goals and at the same time it can be easily adjusted to address short-term requirements without retraining. Next, we present an approximation of our algorithm which can be applied to address short-term requirements of any pre-trained RL algorithm. The case studies demonstrate that both our proposed method as well as its practical approximation can achieve short-term and long-term goals without complex reward functions.


title: Spin wave propagation in uniform waveguide: effects, modulation and its application

abstract: With the advent of the post-Moore era, researches on beyond-Complementary Metal Oxide Semiconductor (CMOS) approaches have been attracted more and more attention. Magnonics, or spin wave is one of the most promising technology beyond CMOS, which magnons-quanta for spin waves-process the information analogous to electronic charges in electronics. Information transmission by spin waves, which uses the frequency, amplitude and (or) phase to encode information, has a great many of advantages such as extremely low energy loss and wide-band frequency. Moreover, using the nonlinear characteristics of spin waves for information transmission can increase the extra degree of freedom of information. This review provides a tutorial overview over the effects of spin wave propagation and recent research progress in uniform spin wave waveguide. The propagation characteristics of spin waves in uniform waveguides and some special propagation phenomena such as spin wave beam splitting and self-focusing are described by combining experimental phenomena and theoretical formulas. Furthermore, we summarize methods for modulating propagation of spin wave in uniform waveguide, and comment on the advantages and limitations of these methods. The review may promote the development of information transmission technology based on spin waves.


title: Chickens and Dukes

abstract: Following on the King Chicken Theorems originally proved by Maurer, we examine the idea of multiple flocks of chickens by bringing the chickens from tournaments to multipartite tournaments. As Kings have already been studied in multipartite settings, notably by Koh-Tan and Petrovic-Thomassen, we examine a new type of chicken more suited than Kings for these multipartite graphs: Dukes. We define an M-Duke to be a vertex from which any vertex in a different partite set is accessible by a directed path of length at most M. In analogy with Maurer's paper, we prove various structural results regarding Dukes. In particular, we prove the existence of 3-Dukes in all multipartite tournaments, and we conclude by proving that in any multipartite tournament, either there is a 1-Duke, three 2-Dukes, or four 3-Dukes.


title: Delve into the Performance Degradation of Differentiable Architecture Search

abstract: Differentiable architecture search (DARTS) is widely considered to be easy to overfit the validation set which leads to performance degradation. We first employ a series of exploratory experiments to verify that neither high-strength architecture parameters regularization nor warmup training scheme can effectively solve this problem. Based on the insights from the experiments, we conjecture that the performance of DARTS does not depend on the well-trained supernet weights and argue that the architecture parameters should be trained by the gradients which are obtained in the early stage rather than the final stage of training. This argument is then verified by exchanging the learning rate schemes of weights and parameters. Experimental results show that the simple swap of the learning rates can effectively solve the degradation and achieve competitive performance. Further empirical evidence suggests that the degradation is not a simple problem of the validation set overfitting but exhibit some links between the degradation and the operation selection bias within bilevel optimization dynamics. We demonstrate the generalization of this bias and propose to utilize this bias to achieve an operation-magnitude-based selective stop.


title: A unified differential equation solver approach for separable convex optimization: splitting, acceleration and nonergodic rate

abstract: This paper provides a self-contained ordinary differential equation solver approach for separable convex optimization problems. A novel primal-dual dynamical system with built-in time rescaling factors is introduced, and the exponential decay of a tailored Lyapunov function is established. Then several time discretizations of the continuous model are considered and analyzed via a unified discrete Lyapunov function. Moreover, two families of accelerated proximal alternating direction methods of multipliers are obtained, and nonergodic optimal mixed-type convergence rates shall be proved for the primal objective residual, the feasibility violation and the Lagrangian gap. Finally, numerical experiments are provided to validate the practical performances.


title: Finite temperature and delta-regime in the 2-flavor Schwinger model

abstract: The Schwinger model is often used as a testbed for conceptual and numerical approaches in lattice field theory. Still, some of its rich physical properties in anisotropic volumes have not yet been explored. For the multi-flavor finite temperature Schwinger model there is an approximate solution by Hosotani et al. based on bosonization. We perform lattice simulations and check the validity of this approximation in the case of two flavors. Next we exchange the r\^{o}le of the coordinates to enter the $\delta$-regime, and measure the dependence of the residual "pion" mass on the spatial size, at zero temperature. Our results show that universal features, which were derived by Leutwyler, Hasenfratz and Niedermayer referring to quasi-spontaneous symmetry breaking in $d>2$, extend even to $d=2$. This enables the computation of the Schwinger model counterpart of the pion decay constant $F_{\pi}$. It is consistent with an earlier determination by Harada et al. who considered the divergence of the axial current in a light-cone formulation, and with analytical results that we conjecture from 2d versions of the Witten--Veneziano formula and the Gell-Mann--Oakes--Renner relation, which suggest $F_{\pi} = 1/ \sqrt{2\pi}$.


title: Difference-Frequency Chirped-Pulse Dual-comb Generation in the THz Region: Temporal Magnification of the Quantum Dynamics of Water Vapor Lines by >60,000

abstract: A new difference-frequency method based on electro-optic phase modulators and two free-running lasers is reported to perform THz dual-comb spectroscopy.


title: Experimentally Validated Simulations of 50 {\mu}m X-ray PIV Tracer Particles

abstract: We evaluate Beer-Lambert (BL) ray-tracing and Monte Carlo N-Particle (MCNP) photon tracking simulations for prediction and comparison of X-ray imaging system performance. These simulation tools can aid the methodical design of laboratory-scale X-ray particle image velocimetry (XPIV) experiments and tracer particles by predicting image quality. Particle image signal-to-noise ratio (SNR) is used as the metric of system performance. Simulated and experiment data of hollow, silver-coated, glass sphere tracer particles (AGSF-33) are compared. As predicted by the simulations, the AGSF-33 particles are visible with a SNR greater than unity in 100~ms exposure time images, demonstrating their potential as X-ray PIV or particle tracking velocimetry (XPTV) tracers. The BL approach predicts the image contrast, is computationally inexpensive, and enables the exploration of a vast parameter space for system design. MCNP simulations, on the other hand, predict experiment images slightly more accurately, but are more than an order of magnitude more computationally expensive than BL simulations. For most practical XPIV system design applications, the higher computational expense of MCNP is likely not justified by the modest accuracy improvement compared to BL.


title: An Automated Approach to Causal Inference in Discrete Settings

abstract: When causal quantities cannot be point identified, researchers often pursue partial identification to quantify the range of possible values. However, the peculiarities of applied research conditions can make this analytically intractable. We present a general and automated approach to causal inference in discrete settings. We show causal questions with discrete data reduce to polynomial programming problems, and we present an algorithm to automatically bound causal effects using efficient dual relaxation and spatial branch-and-bound techniques. The user declares an estimand, states assumptions, and provides data (however incomplete or mismeasured). The algorithm then searches over admissible data-generating processes and outputs the most precise possible range consistent with available information -- i.e., sharp bounds -- including a point-identified solution if one exists. Because this search can be computationally intensive, our procedure reports and continually refines non-sharp ranges that are guaranteed to contain the truth at all times, even when the algorithm is not run to completion. Moreover, it offers an additional guarantee we refer to as $\epsilon$-sharpness, characterizing the worst-case looseness of the incomplete bounds. Analytically validated simulations show the algorithm accommodates classic obstacles, including confounding, selection, measurement error, noncompliance, and nonresponse.


title: Fast radio bursts as probes of feedback from active galactic nuclei

abstract: Fast Radio Bursts (FRBs) are a promising tool for studying the low-density universe as their dispersion measures (DM) are extremely sensitive probes of electron column density. Active Galactic Nuclei (AGN) inject energy into the intergalactic medium, affecting the DM and its scatter. To determine the effectiveness of FRBs as a probe of AGN feedback, we analysed three different AGN models from the EAGLE simulation series. We measured the mean DM-redshift relation, and the scatter around it, using $2.56 \times 10^8$ sightlines at 131 redshift ($z$) bins between $0 \leq z \leq 3$. While the DM-redshift relation itself is highly robust against different AGN feedback models, significant differences are detected in the scatter around the mean: weaker feedback leads to more scatter. We find $\sim 10^4$ localised FRBs are needed to discriminate between the scatter in standard feedback and stronger, more intermittent feedback models. The number of FRBs required is dependent on the redshift distribution of the detected population. A log-normal redshift distribution at $z=0.5$ requires approximately 50% fewer localised FRBs than a distribution centred at $z=1$. With the Square Kilometre Array expected to detect $>10^3$ FRBs per day, in the future, FRBs will be able to provide constraints on AGN feedback.


title: Two time-stepping schemes for sub-diffusion equations with singular source terms

abstract: Singular source terms in sub-diffusion equations may lead to the unboundedness of solutions, which will bring a severe reduction of convergence order of existing time-stepping schemes. In this work, we propose two efficient time-stepping schemes for solving sub-diffusion equations with a class of source terms mildly singular in time. One discretization is based on the Gr{\"u}nwald-Letnikov and backward Euler methods. First-order error estimate with respect to time is rigorously established for singular source terms and nonsmooth initial data. The other scheme derived from the second-order backward differentiation formula (BDF) is proved to possess second-order accuracy in time. Further, piecewise linear finite element and lumped mass finite element discretizations in space are applied and analyzed rigorously. Numerical investigations confirm our theoretical results.


title: Localization to delocalization transition in a double stranded helical geometry: Effects of conformation, transverse electric field and dynamics

abstract: Conformational effect on electronic localization is critically investigated for the first time considering a double-stranded helical geometry (DSHG) subjected to an electric field. In the presence of electric field the DSHG behaves like a correlated disordered system whose site potentials are modulated in a cosine form like the well known Aubry-Andre-Harper (AAH) model. The potential distribution can be modulated further by changing the orientation of the incident field. A similar kind of cosine modulation is also introduced in the inter-strand hopping integrals of the DSHG. Suitably adjusting the orientation of the electric field, we can achieve fully extended energy eigenstates or completely localized ones or a mixture of both. The effects of short-range and long-range hopping integrals along with the chirality on localization are thoroughly studied. Finally, we inspect the role of helical dynamics to make the model more realistic. The interplay between the helical geometry and electric field may open up several notable features of electronic localization and can be verified by using different chiral molecules.

categories: cond-mat.mes-hall cond-mat.dis-nn cond-mat.mtrl-sci physics.app-ph quant-ph