title: InstaTune: Instantaneous Neural Architecture Search During Fine-Tuning

abstract: One-Shot Neural Architecture Search (NAS) algorithms often rely on training a hardware agnostic super-network for a domain specific task. Optimal sub-networks are then extracted from the trained super-network for different hardware platforms. However, training super-networks from scratch can be extremely time consuming and compute intensive especially for large models that rely on a two-stage training process of pre-training and fine-tuning. State of the art pre-trained models are available for a wide range of tasks, but their large sizes significantly limits their applicability on various hardware platforms. We propose InstaTune, a method that leverages off-the-shelf pre-trained weights for large models and generates a super-network during the fine-tuning stage. InstaTune has multiple benefits. Firstly, since the process happens during fine-tuning, it minimizes the overall time and compute resources required for NAS. Secondly, the sub-networks extracted are optimized for the target task, unlike prior work that optimizes on the pre-training objective. Finally, InstaTune is easy to "plug and play" in existing frameworks. By using multi-objective evolutionary search algorithms along with lightly trained predictors, we find Pareto-optimal sub-networks that outperform their respective baselines across different performance objectives such as accuracy and MACs. Specifically, we demonstrate that our approach performs well across both unimodal (ViT and BERT) and multi-modal (BEiT-3) transformer based architectures.


title: Integration of a new Cryogenic Liquefier into the IB-1 Cryogenic Test Facility

abstract: The increase over the last years of the testing activities related to quantum systems, SRF cavities for the PIP-II and the LCLS-II projects, as well as superconducting magnets for the HL-LHC project and Fusion research activities, has required the addition of a new Helium cryogenic plant into the existing IB-1 Industrial Cryogenic Test Facility. The new cryogenic plant is composed of a cryogenic liquefier (Cold Box) able to provide up to 340 L/h, a 4 kL Dewar and two Mycom compressors providing up to 120 g/s. AL-AT (Air Liquide Advanced Technologies) has taken part of this project by designing and manufacturing the cryogenic liquefier. This new cryogenic plant is connected through a cryogenic distribution system to a 10 kL Dewar, which is part of the existing cryogenic test facility, itself composed of another Cold Box and a Sullair compressor. The new cryogenic plant has two main operating modes: one allows to transfer liquid helium at 1.7 bar between the two Dewars, the other allows to transfer supercritical Helium at 2 bar or more between the new Cold Box and the 10 kL Dewar. The entire industrial cryogenic facility is handled by a common Inventory Control System, composed of three control valves, and 9 tanks giving a total buffer volume of more than 1000 m3. This paper presents the technical features of the new Helium cryogenic plant, as well as the main results of the liquefier commissioning phase and details of the helium transfer between the two Dewars, making the connection between the cryogenic plants at the IB-1 Industrial Cryogenic Test Facility.


title: Laplacian $\{-1,0,1\}$- and $\{-1,1\}$-diagonalizable graphs

abstract: A graph is called "Laplacian integral" if the eigenvalues of its Laplacian matrix are all integers. We investigate the subset of these graphs whose Laplacian is furthermore diagonalized by a matrix with entries coming from a fixed set, with particular emphasis on the sets $\{-1,0,1\}$ or $\{-1,1\}$. Such graphs include as special cases the recently-investigated families of "Hadamard-diagonalizable" and "weakly Hadamard-diagonalizable" graphs. As a combinatorial tool to aid in our investigation, we introduce a family of vectors that we call "balanced", which generalize totally balanced partitions, regular sequences, and complete partitions. We show that balanced vectors completely characterize which graph complements and complete multipartite graphs are $\{-1,0,1\}$-diagonalizable, and we furthermore prove results on diagonalizability of the Cartesian product, disjoint union, and join of graphs. Particular attention is paid to the $\{-1,0,1\}$- and $\{-1,1\}$-diagonalizability of the complete graphs and complete multipartite graphs. Finally, we provide a complete list of all simple, connected graphs on nine or fewer vertices that are $\{-1,0,1\}$- or $\{-1,1\}$-diagonalizable.


title: Probing Light Dark Matter with positron beams at NA64

abstract: We present the results of a missing-energy search for Light Dark Matter which has a new interaction with ordinary matter transmitted by a vector boson, called dark photon $A^\prime$. For the first time, this search is performed with a positron beam by using the significantly enhanced production of $A^\prime$ in the resonant annihilation of positrons with atomic electrons of the target nuclei, followed by the invisible decay of $A^\prime$ into dark matter. No events were found in the signal region with $(10.1 \pm 0.1)~\times~10^{9}$ positrons on target with 100 GeV energy. This allowed us to set new exclusion limits that, relative to the collected statistics, prove the power of this experimental technique. This measurement is a crucial first step toward a future exploration program with positron beams, whose estimated sensitivity is here presented.


title: Mixed Variational Flows for Discrete Variables

abstract: Variational flows allow practitioners to learn complex continuous distributions, but approximating discrete distributions remains a challenge. Current methodologies typically embed the discrete target in a continuous space - usually via continuous relaxation or dequantization - and then apply a continuous flow. These approaches involve a surrogate target that may not capture the original discrete target, might have biased or unstable gradients, and can create a difficult optimization problem. In this work, we develop a variational flow family for discrete distributions without any continuous embedding. First, we develop a measure-preserving and discrete (MAD) invertible map that leaves the discrete target invariant, and then create a mixed variational flow (MAD Mix) based on that map. Our family provides access to i.i.d. sampling and density evaluation with virtually no tuning effort. We also develop an extension to MAD Mix that handles joint discrete and continuous models. Our experiments suggest that MAD Mix produces more reliable approximations than continuous-embedding flows while being significantly faster to train.


title: Everything Perturbed All at Once: Enabling Differentiable Graph Attacks

abstract: As powerful tools for representation learning on graphs, graph neural networks (GNNs) have played an important role in applications including social networks, recommendation systems, and online web services. However, GNNs have been shown to be vulnerable to adversarial attacks, which can significantly degrade their effectiveness. Recent state-of-the-art approaches in adversarial attacks rely on gradient-based meta-learning to selectively perturb a single edge with the highest attack score until they reach the budget constraint. While effective in identifying vulnerable links, these methods are plagued by high computational costs. By leveraging continuous relaxation and parameterization of the graph structure, we propose a novel attack method called Differentiable Graph Attack (DGA) to efficiently generate effective attacks and meanwhile eliminate the need for costly retraining. Compared to the state-of-the-art, DGA achieves nearly equivalent attack performance with 6 times less training time and 11 times smaller GPU memory footprint on different benchmark datasets. Additionally, we provide extensive experimental analyses of the transferability of the DGA among different graph models, as well as its robustness against widely-used defense mechanisms.


title: Nitrogen Precooling Heat Exchanger replacement and control system upgrade in Superfluid Cryoplant at CMTF

abstract: Liquid Nitrogen precooling is used in most Cryoplants to achieve cooldown to 80 K temperature range. In one such system at Fermilab's CMTF Superfluid Cryoplant, where the Helium supply directly exchanges heat with liquid Nitrogen, freezing of Nitrogen occurred inside the heat exchanger due to heat exchanger imbalance during a Cryoplant trip. Trapped vapor pockets of N2 within the frozen heat exchanger channels were formed while warming up the heat exchanger, creating high localized pressure and subsequent damage/rupture of the heat exchanger. Replacement of the heat exchanger was done, and modifications were made in the system to rectify future occurrences. The control system was updated to bypass the heat exchanger entirely if the incoming Helium stream temperature drops below 76 K. This was done by repurposing two control valves as heat exchanger bypass valves that were previously used for a redundant 80 K adsorber in the coldbox. Additional modifications were made to further prevent return of large amount of cold Helium gas from cold end during abrupt Cryoplant shutdown. This modification has ensured high reliability of heat exchanger with prevention of freezing of Nitrogen which can damage the heat exchanger.


title: Design of Muon Campus full flow purifier for varying operational conditions and horizontal shipping

abstract: Constant ingress of impurities in Muon Campus g-2 experiment at Fermilab has resulted in reduction of efficiency of cryogenic expanders and occasional undesired downtime to flush the impurities. Due to insufficiency of current 60 g/s mobile purifier, a full flow purifier is designed to be used in Muon Campus which purifies 240 g/s of Helium throughput of 4 compressors through charcoal bed at 80 K and returns ambient Helium back to the system. The purifier is designed to be operated near liquid Nitrogen temperature during cold operations and up to 400 K during regeneration. Both warm and cold operational range of the purifier has required use of appropriate clearances in design due to expansion and contraction. The vessel of around 16 ft height which is designed to be operated vertically is to be shipped horizontally. The asymmetrical position of heavy stainless steel heat exchanger in the purifier support frame and 5g vertical load design consideration for shipping has required use of shipping supports and heat exchanger rotational stops to comply with design requirements. FEA of purifier system is performed in cold, warm and shipping cases to verify that the purifier satisfies the design requirements.


title: Streaming, Local, and Multi-Level (Hyper)Graph Decomposition

abstract: (Hyper)Graph decomposition is a family of problems that aim to break down large (hyper)graphs into smaller sub(hyper)graphs for easier analysis. The importance of this lies in its ability to enable efficient computation on large and complex (hyper)graphs, such as social networks, chemical compounds, and computer networks. This dissertation explores several types of (hyper)graph decomposition problems, including graph partitioning, hypergraph partitioning, local graph clustering, process mapping, and signed graph clustering. Our main focus is on streaming algorithms, local algorithms and multilevel algorithms. In terms of streaming algorithms, we make contributions with highly efficient and effective algorithms for (hyper)graph partitioning and process mapping. In terms of local algorithms, we propose sub-linear algorithms which are effective in detecting high-quality local communities around a given seed node in a graph based on the distribution of a given motif. In terms of multilevel algorithms, we engineer high-quality multilevel algorithms for process mapping and signed graph clustering. We provide a thorough discussion of each algorithm along with experimental results demonstrating their superiority over existing state-of-the-art techniques. The results show that the proposed algorithms achieve improved performance and better solutions in various metrics, making them highly promising for practical applications. Overall, this dissertation showcases the effectiveness of advanced combinatorial algorithmic techniques in solving challenging (hyper)graph decomposition problems.


title: RACR-MIL: Weakly Supervised Skin Cancer Grading using Rank-Aware Contextual Reasoning on Whole Slide Images

abstract: Cutaneous squamous cell cancer (cSCC) is the second most common skin cancer in the US. It is diagnosed by manual multi-class tumor grading using a tissue whole slide image (WSI), which is subjective and suffers from inter-pathologist variability. We propose an automated weakly-supervised grading approach for cSCC WSIs that is trained using WSI-level grade and does not require fine-grained tumor annotations. The proposed model, RACR-MIL, transforms each WSI into a bag of tiled patches and leverages attention-based multiple-instance learning to assign a WSI-level grade. We propose three key innovations to address general as well as cSCC-specific challenges in tumor grading. First, we leverage spatial and semantic proximity to define a WSI graph that encodes both local and non-local dependencies between tumor regions and leverage graph attention convolution to derive contextual patch features. Second, we introduce a novel ordinal ranking constraint on the patch attention network to ensure that higher-grade tumor regions are assigned higher attention. Third, we use tumor depth as an auxiliary task to improve grade classification in a multitask learning framework. RACR-MIL achieves 2-9% improvement in grade classification over existing weakly-supervised approaches on a dataset of 718 cSCC tissue images and localizes the tumor better. The model achieves 5-20% higher accuracy in difficult-to-classify high-risk grade classes and is robust to class imbalance.


title: Intensity correlation holography for remote phase sensing and 3D imaging

abstract: Holography is an established technique for measuring the wavefront of optical signals through interferometric combination with a reference wave. Conventionally the integration time of a hologram is limited by the interferometer coherence time, thus making it challenging to prepare holograms of remote objects, especially using weak illumination. Here, we circumvent this limitation by using intensity correlation interferometry. Although the exposure time of individual holograms must be shorter than the interferometer coherence time, we show that any number of randomly phase-shifted holograms can be combined into a single intensity-correlation hologram. In a proof-of-principle experiment, we use this technique to perform phase imaging and 3D reconstruction of an object at a ~3m distance using weak illumination and without active phase stabilization.


title: Intelligent System for Assessing University Student Personality Development and Career Readiness

abstract: While academic metrics such as transcripts and GPA are commonly used to evaluate students' knowledge acquisition, there is a lack of comprehensive metrics to measure their preparedness for the challenges of post-graduation life. This research paper explores the impact of various factors on university students' readiness for change and transition, with a focus on their preparedness for careers. The methodology employed in this study involves designing a survey based on Paul J. Mayer's "The Balance Wheel" to capture students' sentiments on various life aspects, including satisfaction with the educational process and expectations of salary. The collected data from a KBTU student survey (n=47) were processed through machine learning models: Linear Regression, Support Vector Regression (SVR), Random Forest Regression. Subsequently, an intelligent system was built using these models and fuzzy sets. The system is capable of evaluating graduates' readiness for their future careers and demonstrates a high predictive power. The findings of this research have practical implications for educational institutions. Such an intelligent system can serve as a valuable tool for universities to assess and enhance students' preparedness for post-graduation challenges. By recognizing the factors contributing to students' readiness for change, universities can refine curricula and processes to better prepare students for their career journeys.


title: A hybridizable discontinuous Galerkin method for the coupled Navier-Stokes/Biot problem

abstract: In this paper we present a hybridizable discontinuous Galerkin method for the time-dependent Navier-Stokes equations coupled to the quasi-static poroelasticity equations via interface conditions. We determine a bound on the data that guarantees stability and well-posedness of the fully discrete problem and prove a priori error estimates. A numerical example confirms our analysis.


title: Flexible Handover with Real-Time Robust Dynamic Grasp Trajectory Generation

abstract: In recent years, there has been a significant effort dedicated to developing efficient, robust, and general human-to-robot handover systems. However, the area of flexible handover in the context of complex and continuous objects' motion remains relatively unexplored. In this work, we propose an approach for effective and robust flexible handover, which enables the robot to grasp moving objects with flexible motion trajectories with a high success rate. The key innovation of our approach is the generation of real-time robust grasp trajectories. We also design a future grasp prediction algorithm to enhance the system's adaptability to dynamic handover scenes. We conduct one-motion handover experiments and motion-continuous handover experiments on our novel benchmark that includes 31 diverse household objects. The system we have developed allows users to move and rotate objects in their hands within a relatively large range. The success rate of the robot grasping such moving objects is 78.15% over the entire household object benchmark.


title: Discovery of Spherules of Likely Extrasolar Composition in the Pacific Ocean Site of the CNEOS 2014-01-08 (IM1) Bolide

abstract: We have conducted an extensive towed-magnetic-sled survey during the period 14-28 June, 2023, over the seafloor centered around the calculated path of the bolide CNEOS 2014-01-08 (IM1) about 85 km north of Manus Island, Papua New Guinea. We found about 700 spherules of diameter 0.05-1.3 millimeters in our samples, of which 57 were analyzed so far. The spherules were significantly concentrated along the expected meteor path. Mass spectrometry of 47 spherules near the high-yield regions along IM1's path reveals a distinct extra-solar abundance pattern for 5 of them, while background spherules have abundances consistent with a solar system origin. The unique spherules show an excess of Be, La and U, by up to three orders of magnitude relative to the solar system standard of CI chondrites. These "BeLaU"-type spherules, never seen before, also have very low refractory siderophile elements such as Re. Volatile elements, such as Mn, Zn, Pb, are depleted as expected from evaporation losses during a meteor's airburst. In addition, the mass-dependent variations in $^{57}$Fe/$^{54}$Fe and $^{56}$Fe/$^{54}$Fe are also consistent with evaporative loss of the light isotopes during the spherules' travel in the atmosphere. The "BeLaU" abundance pattern is not found in control regions outside of IM1's path and does not match commonly manufactured alloys or natural meteorites in the solar system. This evidence points towards an association of "BeLaU"-type spherules with IM1, supporting its interstellar origin independently of the high velocity and unusual material strength implied from the CNEOS data. We suggest that the "BeLaU" abundance pattern could have originated from a highly differentiated magma ocean of a planet with an iron core outside the solar system or from more exotic sources.


title: Detection of Mild Cognitive Impairment Using Facial Features in Video Conversations

abstract: Early detection of Mild Cognitive Impairment (MCI) leads to early interventions to slow the progression from MCI into dementia. Deep Learning (DL) algorithms could help achieve early non-invasive, low-cost detection of MCI. This paper presents the detection of MCI in older adults using DL models based only on facial features extracted from video-recorded conversations at home. We used the data collected from the I-CONECT behavioral intervention study (NCT02871921), where several sessions of semi-structured interviews between socially isolated older individuals and interviewers were video recorded. We develop a framework that extracts spatial holistic facial features using a convolutional autoencoder and temporal information using transformers. Our proposed DL model was able to detect the I-CONECT study participants' cognitive conditions (MCI vs. those with normal cognition (NC)) using facial features. The segments and sequence information of the facial features improved the prediction performance compared with the non-temporal features. The detection accuracy using this combined method reached 88% whereas 84% is the accuracy without applying the segments and sequences information of the facial features within a video on a certain theme.


title: Sperner theorems for unrelated copies of some partially ordered sets in a powerset lattice and minimum generating sets of powers of distributive lattices

abstract: For a finite poset (partially ordered set) $U$ and a natural number $n$, let Sp$(U,n)$ denote the largest number of pairwise unrelated copies of $U$ in the powerset lattice (AKA subset lattice) of an $n$-element set. If $U$ is the singleton poset, then Sp$(U,n)$ was determined by E. Sperner in 1928; this result is well known in extremal combinatorics. Later, exactly or asymptotically, Sperner's theorem was extended to other posets by A. P. Dove, J. R. Griggs, G. O. H. Katona, D Nagy, J. Stahl, and W. T. Jr. Trotter. We determine Sp$(U,n)$ for all finite posets with 0 and 1, and we give reasonable estimates for the ``V-shaped'' 3-element poset and the 4-element poset with 0 and three maximal elements. For a lattice $L$, let Gmin($L$) denote the minimum size of generating sets of $L$. We prove that if $U$ is the poset of the join-irreducible elements of a finite distributive lattice $D$, then the function $k\mapsto$ Gmin($D^k)$ is the left adjoint of the function $n\mapsto$ Sp$(U,n)$. This allows us to determine Gmin($D^k)$ in many cases. E.g., for a 5-element distributive lattice $D$, Gmin($D^{2023})=18$ if $D$ is a chain and Gmin($D^{2023})=15$ otherwise. It follows that large direct powers of small distributive lattices are appropriate for our 2021 cryptographic authentication protocol.


title: Analytical One-Dimensional Model of Drop Ejection from a Micro-Size Nozzle

abstract: In this article, we construct a novel one-dimensional model of drop ejection from a micro-size nozzle due to a short pressure pulse applied to the liquid in the nozzle. The pressure pulse supplies the kinetic energy to the perturbed liquid-bulge squeezed from the nozzle, which then ballistically lengthens forming a ligament. The Plateau-Rayleigh instability forms a neck in the ligament at the nozzle, leading to detachment of the ligament from the nozzle, which then collapses in a drop. This drop formation sequence is typical for drop-on-demand printheads in which the drop is ejected from the nozzle by a short pressure pulse at the needed moment of time when it should reach the substrate. The model calculates the velocity of the droplet, length of the ligament vs. time, and the time when the ligament detaches from the nozzle as a function of the exit radius of the nozzle, the volume of the droplet, the time that the volume of the droplet is squeezed from the nozzle, viscosity, surface tension, and mass density of the liquid drop. The model also calculates a criterion for drop ejection from the nozzle.


title: Target PCA: Transfer Learning Large Dimensional Panel Data

abstract: This paper develops a novel method to estimate a latent factor model for a large target panel with missing observations by optimally using the information from auxiliary panel data sets. We refer to our estimator as target-PCA. Transfer learning from auxiliary panel data allows us to deal with a large fraction of missing observations and weak signals in the target panel. We show that our estimator is more efficient and can consistently estimate weak factors, which are not identifiable with conventional methods. We provide the asymptotic inferential theory for target-PCA under very general assumptions on the approximate factor model and missing patterns. In an empirical study of imputing data in a mixed-frequency macroeconomic panel, we demonstrate that target-PCA significantly outperforms all benchmark methods.


title: The EUSO-SPB2 Cherenkov Telescope -- Flight Performance and Preliminary Results

abstract: Astrophysical Very-High-Energy (VHE, >10PeV) neutrinos deliver crucial information about the sources of Ultra-High-Energy Cosmic Rays (UHECRs), the composition of UHECRs, and neutrino/particle physics at highest energies. UHE-tau neutrinos skimming the Earth's surface produce tau leptons, which can emerge from the ground, decay, and start an upward-going extensive air shower (EAS) in the Earth's atmosphere. The tau neutrino can be reconstructed by imaging the EAS. We developed an atmospheric Cherenkov Telescope flying on the Extreme Universe Space Observatory Super Pressure Balloon 2 (EUSO-SPB2) mission to test the air-shower imaging concept at highest altitudes. The EUSO-SPB2 ultra-long-duration balloon mission is a precursor of the Probe of Extreme Multi-Messenger Astrophysics (POEMMA), a candidate for an astrophysics probe-class mission. The telescope implements Schmidt optics with a 0.785 m^2 light collection area and a 512-pixel SiPM camera covering a 12.8{\deg} by 6.4{\deg} (Horizontal by Vertical) field-of-view with 0.4{\deg} resolution. The camera signals are sampled with 100MSa/s and digitized with 12-bit resolution. The objectives of the EUSO-SPB2 Cherenkov telescope include a search for UHE neutrinos below Earth's limb, UHECRs above the limb, the study of the night sky background, and studying the telescope's performance. In this presentation, I will present an overview of the Cherenkov telescope and discuss the in-flight performance of the telescope.


title: Dynamic random intersection graph: Dynamic local convergence and giant structure

abstract: Random intersection graphs containing an underlying community structure are a popular choice for modelling real-world networks. Given the group memberships, the classical random intersection graph is obtained by connecting individuals when they share at least one group. We extend this approach and make the communities dynamic by letting them alternate between an active and inactive phase. We analyse the new model, delivering results on degree distribution, local convergence, giant component, and maximum group size, paying particular attention to the dynamic description of these properties. We also describe the connection between our model and the bipartite configuration model, which is of independent interest.


title: Examining the sensitivity of FASER to Generalized Neutrino Interactions

abstract: We investigate the sensitivity of the FASER$\nu$ detector, a novel experimental setup at the LHC, to probe and constrain generalized neutrino interactions (GNI). Employing a comprehensive theoretical framework, we model the effects of generalized neutrino interactions on neutrino-nucleon deep inelastic scattering processes within the FASER$\nu$ detector. By considering all the neutrino channels produced at the LHC, we perform a statistical analysis to determine the sensitivity of FASER$\nu$ to constrain these interactions. Our results demonstrate that FASER$\nu$ can place stringent constraints on the GNI effective couplings. Additionally, we study the relation between GNI and a minimal Leptoquark model where the SM is augmented by a singlet Leptoquark with hypercharge $1/3$. We have found that the sensitivities for various combinations of the Leptoquark Yukawa couplings are approximately $\mathcal{O}(1)$, particularly when considering a Leptoquark mass in the TeV range.


title: Low-noise single-photon counting superconducting nanowire detectors at infrared wavelengths up to 29 $\mu$m

abstract: We report on the extension of the spectral sensitivity of superconducting nanowire single-photon detectors to a wavelength of 29 $\mu$m. This represents the first demonstration of a time correlated single-photon counting detector at these long infrared wavelengths. We achieve saturated internal detection efficiency from 10 to 29 $\mu$m, whilst maintaining dark count rates below 0.1 counts per second. Extension of superconducting nanowire single-photon detectors to this spectral range provides low noise and high timing resolution photon counting detection, effectively providing a new class of single-photon sensitive detector for these wavelengths. These detectors are important for applications such as exoplanet spectroscopy, infrared astrophysics, physical chemistry, remote sensing and direct dark-matter detection.


title: The science case for a far-infrared interferometer in the era of JWST and ALMA

abstract: A space-based far-infrared interferometer could work synergistically with the James Webb Space Telescope (JWST) and the Atacama Large Millimeter Array (ALMA) to revolutionize our understanding of the astrophysical processes leading to the formation of habitable planets and the co-evolution of galaxies and their central supermassive black holes. Key to these advances are measurements of water in its frozen and gaseous states, observations of astronomical objects in the spectral range where most of their light is emitted, and access to critical diagnostic spectral lines, all of which point to the need for a far-infrared observatory in space. The objects of interest - circumstellar disks and distant galaxies - typically appear in the sky at sub-arcsecond scales, which rendered all but a few of them unresolvable with the successful and now-defunct 3.5-m Herschel Space Observatory, the largest far-infrared telescope flown to date. A far-infrared interferometer with maximum baseline length in the tens of meters would match the angular resolution of JWST at 10x longer wavelengths and observe water ice and water-vapor emission, which ALMA can barely do through the Earth's atmosphere. Such a facility was conceived and studied two decades ago. Here we revisit the science case for a space-based far-infrared interferometer in the era of JWST and ALMA and summarize the measurement capabilities that will enable the interferometer to achieve a set of compelling scientific objectives. Common to all the science themes we consider is a need for sub-arcsecond image resolution.


title: The Impact of Reference-Command Preview on Human-in-the-Loop Control Behavior

abstract: This article presents results from an experiment in which 44 human subjects interact with a dynamic system to perform 40 trials of a command-following task. The reference command is unpredictable and different on each trial, but all subjects have the same sequence of reference commands for the 40 trials. The subjects are divided into 4 groups of 11 subjects. One group performs the command-following task without preview of the reference command, and the other 3 groups are given preview of the reference command for different time lengths into the future (0.5 s, 1 s, 1.5 s). A subsystem identification algorithm is used to obtain best-fit models of each subject's control behavior on each trial. The time- and frequency-domain performance, as well as the identified models of the control behavior for the 4 groups are examined to investigate the effects of reference-command preview. The results suggest that preview tends to improve performance by allowing the subjects to compensate for sensory time delay and approximate the inverse dynamics in feedforward. However, too much preview may decrease performance by degrading the ability to use the correct phase lead in feedforward.


title: Ab initio uncertainty quantification of neutrinoless double-beta decay in $^{76}$Ge

abstract: The observation of neutrinoless double-beta ($0\nu\beta\beta$) decay would offer proof of lepton number violation, demonstrating that neutrinos are Majorana particles, while also helping us understand why there is more matter than antimatter in the Universe. If the decay is driven by the exchange of the three known light neutrinos, a discovery would, in addition, link the observed decay rate to the neutrino mass scale through a theoretical quantity known as the nuclear matrix element (NME). Accurate values of the NMEs for all nuclei considered for use in $0\nu\beta\beta$ experiments are therefore crucial for designing and interpreting those experiments. Here, we report the first comprehensive ab initio uncertainty quantification of the $0\nu\beta\beta$-decay NME, in the key nucleus $^{76}$Ge. Our method employs nuclear strong and weak interactions derived within chiral effective field theory and recently developed many-body emulators. Our result, with a conservative treatment of uncertainty, is an NME of $2.60^{+1.28}_{-1.36}$, which, together with the best-existing half-life sensitivity and phase-space factor, sets an upper limit for effective neutrino mass of $187^{+205}_{-62}$ meV. The result is important for designing next-generation germanium detectors aiming to cover the entire inverted hierarchy region of neutrino masses.


title: Parameterized and Approximation Algorithms for the Maximum Bimodal Subgraph Problem

abstract: A vertex of a plane digraph is bimodal if all its incoming edges (and hence all its outgoing edges) are consecutive in the cyclic order around it. A plane digraph is bimodal if all its vertices are bimodal. Bimodality is at the heart of many types of graph layouts, such as upward drawings, level-planar drawings, and L-drawings. If the graph is not bimodal, the Maximum Bimodal Subgraph (MBS) problem asks for an embedding-preserving bimodal subgraph with the maximum number of edges. We initiate the study of the MBS problem from the parameterized complexity perspective with two main results: (i) we describe an FPT algorithm parameterized by the branchwidth (and hence by the treewidth) of the graph; (ii) we establish that MBS parameterized by the number of non-bimodal vertices admits a polynomial kernel. As the byproduct of these results, we obtain a subexponential FPT algorithm and an efficient polynomial-time approximation scheme for MBS.


title: Widely Separated MIMO Radar Using Matrix Completion

abstract: We present a low-complexity widely separated multiple-input-multiple-output (WS-MIMO) radar that samples the signals at each of its multiple receivers at reduced rates. We process the low-rate samples of all transmit-receive chains at each receiver as data matrices. We demonstrate that each of these matrices is low rank as long as the target moves slowly within a coherent processing interval. We leverage matrix completion (MC) to recover the missing samples of each receiver signal matrix at the common fusion center. Subsequently, we estimate the targets' positions and Doppler velocities via the maximum likelihood method. Our MC-WS-MIMO approach recovers missing samples and thereafter target parameters at reduced rates without discretization. Our analysis using ambiguity functions shows that antenna geometry affects the performance of MC-WS-MIMO. Numerical experiments demonstrate reasonably accurate target localization at SNR of 20 dB and sampling rate reduction to 20%.


title: Managing Software Provenance to Enhance Reproducibility in Computational Research

abstract: Scientific processes rely on software as an important tool for data acquisition, analysis, and discovery. Over the years sustainable software development practices have made progress in being considered as an integral component of research. However, management of computation-based scientific studies is often left to individual researchers who design their computational experiments based on personal preferences and the nature of the study. We believe that the quality, efficiency, and reproducibility of computation-based scientific research can be improved by explicitly creating an execution environment that allows researchers to provide a clear record of traceability. This is particularly relevant to complex computational studies in high-performance computing (HPC) environments. In this article, we review the documentation required to maintain a comprehensive record of HPC computational experiments for reproducibility. We also provide an overview of tools and practices that we have developed to perform such studies around Flash-X, a multi-physics scientific software.


title: A Task-Parallel Approach for Localized Topological Data Structures

abstract: Unstructured meshes are characterized by data points irregularly distributed in the Euclidian space. Due to the irregular nature of these data, computing connectivity information between the mesh elements requires much more time and memory than on uniformly distributed data. To lower storage costs, dynamic data structures have been proposed. These data structures compute connectivity information on the fly and discard them when no longer needed. However, on-the-fly computation slows down algorithms and results in a negative impact on the time performance. To address this issue, we propose a new task-parallel approach to proactively compute mesh connectivity. Unlike previous approaches implementing data-parallel models, where all threads run the same type of instructions, our task-parallel approach allows threads to run different functions. Specifically, some threads run the algorithm of choice while other threads compute connectivity information before they are actually needed. The approach was implemented in the new Accelerated Clustered TOPOlogical (ACTOPO) data structure, which can support any processing algorithm requiring mesh connectivity information. Our experiments show that ACTOPO combines the benefits of state-of-the-art memory-efficient (TTK CompactTriangulation) and time-efficient (TTK ExplicitTriangulation) topological data structures. It occupies a similar amount of memory as TTK CompactTriangulation while providing up to 5x speedup. Moreover, it achieves comparable time performance as TTK ExplicitTriangulation while using only half of the memory space.


abstract: Deep Learning is mostly responsible for the surge of interest in Artificial Intelligence in the last decade. So far, deep learning researchers have been particularly successful in the domain of image processing, where Convolutional Neural Networks are used. Although excelling at image classification, Convolutional Neural Networks are quite naive in that no inductive bias is set on the embedding space for images. Similar flaws are also exhibited by another type of Convolutional Networks - Graph Convolutional Neural Networks. However, using non-Euclidean space for embedding data might result in more robust and explainable models. One example of such a non-Euclidean space is hyperbolic space. Hyperbolic spaces are particularly useful due to their ability to fit more data in a low-dimensional space and tree-likeliness properties. These attractive properties have been previously used in multiple papers which indicated that they are beneficial for building hierarchical embeddings using shallow models and, recently, using MLPs and RNNs. However, no papers have yet suggested a general approach to using Hyperbolic Convolutional Neural Networks for structured data processing, although these are the most common examples of data used. Therefore, the goal of this work is to devise a general recipe for building Hyperbolic Convolutional Neural Networks. We hypothesize that ability of hyperbolic space to capture hierarchy in the data would lead to better performance. This ability should be particularly useful in cases where data has a tree-like structure. Since this is the case for many existing datasets \citep{wordnet, imagenet, fb15k}, we argue that such a model would be advantageous both in terms of applications and future research prospects.


title: Identifying Constitutive Parameters for Complex Hyperelastic Materials using Physics-Informed Neural Networks

abstract: Identifying constitutive parameters in engineering and biological materials, particularly those with intricate geometries and mechanical behaviors, remains a longstanding challenge. The recent advent of Physics-Informed Neural Networks (PINNs) offers promising solutions, but current frameworks are often limited to basic constitutive laws and encounter practical constraints when combined with experimental data. In this paper, we introduce a robust PINN-based framework designed to identify material parameters for soft materials, specifically those exhibiting complex constitutive behaviors, under large deformation in plane stress conditions. Distinctively, our model emphasizes training PINNs with multi-modal synthetic experimental datasets consisting of full-field deformation and loading history, ensuring algorithm robustness even with noisy data. Our results reveal that the PINNs framework can accurately identify constitutive parameters of the incompressible Arruda-Boyce model for samples with intricate geometries, maintaining an error below 5%, even with an experimental noise level of 5%. We believe our framework provides a robust modulus identification approach for complex solids, especially for those with geometrical and constitutive complexity.


title: Two-scale modelling of two-phase flows based on the Stationary Action Principle and a Geometric Method Of Moments

abstract: In this contribution, we introduce a versatile formalism to derive unified two-phase models describing both the separated and disperse regimes. It relies on the stationary action principle and interface geometric variables. The main ideas are introduced on a simplified case where all the scales and phases have the same velocity and that does not take into account large-scale capillary forces. The derivation tools yield a proper mathematical framework through hyperbolicity and signed entropy evolution. The formalism encompasses a hierarchy of small-scale reduced-order models based on a statistical description at a mesoscopic kinetic level and is naturally able to include the description of a disperse phase with polydispersity in size. This hierarchy includes both a cloud of spherical droplets and non-spherical droplets experiencing a dynamical behaviour through incompressible oscillations. The associated small-scale variables are moments of a number density function resulting from the Geometric Method Of Moments (GeoMOM). This method selects moments as small-scale geometric variables compatible with the structure and dynamics of the interface; they are defined independently of the flow topology and, therefore, this model pursues the goal of unifying the modelling of a fully-coupled two-scale flow. It is particularly showed that the resulting dynamics provides closures for the interface area density equation obtained from the averaging approach. The extension to mass transfer from one scale to the other including capillary phenomena, as well as the extension to multiple velocities are possible and proposed in complementary works.


title: Gap-Free Clustering: Sensitivity and Robustness of SDP

abstract: We study graph clustering in the Stochastic Block Model (SBM) in the presence of both large clusters and small, unrecoverable clusters. Previous convex relaxation approaches achieving exact recovery do not allow any small clusters of size $o(\sqrt{n})$, or require a size gap between the smallest recovered cluster and the largest non-recovered cluster. We provide an algorithm based on semidefinite programming (SDP) which removes these requirements and provably recovers large clusters regardless of the remaining cluster sizes. Mid-sized clusters pose unique challenges to the analysis, since their proximity to the recovery threshold makes them highly sensitive to small noise perturbations and precludes a closed-form candidate solution. We develop novel techniques, including a leave-one-out-style argument which controls the correlation between SDP solutions and noise vectors even when the removal of one row of noise can drastically change the SDP solution. We also develop improved eigenvalue perturbation bounds of potential independent interest. Our results are robust to certain semirandom settings that are challenging for alternative algorithms. Using our gap-free clustering procedure, we obtain efficient algorithms for the problem of clustering with a faulty oracle with superior query complexities, notably achieving $o(n^2)$ sample complexity even in the presence of a large number of small clusters. Our gap-free clustering procedure also leads to improved algorithms for recursive clustering.

categories: cs.LG cs.DS cs.IT math.IT math.OC stat.ML


title: Searching for TeV emission from LHAASO J0341+5258 with VERITAS and HAWC

abstract: Galactic PeVatrons are astrophysical sources accelerating particles up to a few PeV ($\sim$10$^{15}$ eV) energies. The primary signature of 100 TeV $\gamma$ rays may come from PeV protons or multi-hundred TeV (not PeV) electrons. The search for PeVatrons has been one of the key science topics for VERITAS and HAWC. In 2021, LHAASO detected 14 steady $\gamma$-ray sources with photon energies above 100 TeV, up to 1.4 PeV. This provides a clear list of PeVatron candidates for further study with VERITAS and HAWC. Most of these sources contain possible source associations, such as supernova remnants, pulsar wind nebulae, and stellar clusters. However, two sources: LHAASO J2108+5157 and LHAASO J0341+5258, do not have any such counterparts. Therefore, multiwavelength observations are required to identify the objects responsible for the UHE $\gamma$ rays, to understand the source morphology and association, and to shed light on the emission processes. Here, we will present the status of VERITAS/HAWC observations and results for the LHAASO PeVatron candidate J0341+5258, and also discuss the VERITAS PeVatron search in general.


title: Pixelated high-Q metasurfaces for in-situ biospectroscopy and AI-enabled classification of lipid membrane photoswitching dynamics

abstract: Nanophotonic devices excel at confining light into intense hot spots of the electromagnetic near fields, creating unprecedented opportunities for light-matter coupling and surface-enhanced sensing. Recently, all-dielectric metasurfaces with ultrasharp resonances enabled by photonic bound states in the continuum have unlocked new functionalities for surface-enhanced biospectroscopy by precisely targeting and reading out molecular absorption signatures of diverse molecular systems. However, BIC-driven molecular spectroscopy has so far focused on endpoint measurements in dry conditions, neglecting the crucial interaction dynamics of biological systems. Here, we combine the advantages of pixelated all-dielectric metasurfaces with deep learning-enabled feature extraction and prediction to realize an integrated optofluidic platform for time-resolved in-situ biospectroscopy. Our approach harnesses high-Q metasurfaces specifically designed for operation in a lossy aqueous environment together with advanced spectral sampling techniques to temporally resolve the dynamic behavior of photoswitchable lipid membranes. Enabled by a software convolutional neural network, we further demonstrate the real-time classification of the characteristic cis and trans membrane conformations with 98% accuracy. Our synergistic sensing platform incorporating metasurfaces, optofluidics, and deep learning opens exciting possibilities for studying multi-molecular biological systems, ranging from the behavior of transmembrane proteins to the dynamic processes associated with cellular communication.


title: AskIt: Unified Programming Interface for Programming with Large Language Models

abstract: Large Language Models (LLMs) exhibit a unique phenomenon known as emergent abilities, demonstrating adeptness across numerous tasks, from text summarization to code generation. While these abilities open up novel avenues in software design and crafting, their incorporation presents substantial challenges. Developers face decisions regarding the use of LLMs for directly performing tasks within applications as well as for generating and executing code to accomplish these tasks. Moreover, effective prompt design becomes a critical concern, given the necessity of extracting data from natural language outputs. To address these complexities, this paper introduces AskIt, a domain-specific language (DSL) specifically designed for LLMs. AskIt simplifies LLM integration by providing a unified interface that not only allows for direct task execution using LLMs but also supports the entire cycle of code generation and execution. This dual capability is achieved through (1) type-guided output control, (2) template-based function definitions, and (3) prompt generation for both usage modes. Our evaluations underscore AskIt's effectiveness. Across 50 tasks, AskIt generated concise prompts, achieving a 16.14 % reduction in prompt length compared to benchmarks. Additionally, by enabling a seamless transition between using LLMs directly in applications and for generating code, AskIt achieved significant efficiency improvements, as observed in our GSM8K benchmark experiments. The implementations of AskIt in TypeScript and Python are available at https://github.com/katsumiok/ts-askit and https://github.com/katsumiok/pyaskit, respectively.


title: Universal scaling in far-from-equilibrium quantum systems: an equivalent differential approach

abstract: Recent progress in out-of-equilibrium closed quantum systems has significantly advanced the understanding of mechanisms behind their evolution towards thermalization. Notably, the concept of non-thermal fixed points (NTFPs) - responsible for the emergence of spatio-temporal universal scaling in far-from-equilibrium systems - has played a crucial role in both theoretical and experimental investigations. In this work, we introduce a differential equation that has the universal scaling associated with NTFPs as a solution. The advantage of working with a differential equation, rather than only with its solution, is that we can extract several insightful properties not necessarily present in the solution alone. How the differential equation is derived allows physical interpretation of the universal exponents in terms of the time dependence of the amplitude of the distributions and their momentum scaling. Employing two limiting cases of the equation, we determined the universal exponents related to the scaling using the distributions near just two momentum values. We established a solid agreement with previous investigations by validating this approach with three distinct physical systems. This consistency highlights the universal nature of scaling due to NTFPs and emphasizes the predictive capabilities of the proposed differential equation. Moreover, under specific conditions, the equation predicts a power-law related to the ratio of the two universal exponents, leading to implications concerning particle and energy transport. This suggests that the observed power-laws in far-from-equilibrium turbulent fluids could be related to the universal scaling due to NTFPs, potentially offering new insights into the study of turbulence.


title: A General Recipe for Automated Machine Learning in Practice

abstract: Automated Machine Learning (AutoML) is an area of research that focuses on developing methods to generate machine learning models automatically. The idea of being able to build machine learning models with very little human intervention represents a great opportunity for the practice of applied machine learning. However, there is very little information on how to design an AutoML system in practice. Most of the research focuses on the problems facing optimization algorithms and leaves out the details of how that would be done in practice. In this paper, we propose a frame of reference for building general AutoML systems. Through a narrative review of the main approaches in the area, our main idea is to distill the fundamental concepts in order to support them in a single design. Finally, we discuss some open problems related to the application of AutoML for future research.


title: Universal framework for simultaneous tomography of quantum states and SPAM noise

abstract: We present a general denoising algorithm for performing simultaneous tomography of quantum states and measurement noise. This algorithm allows us to fully characterize state preparation and measurement (SPAM) errors present in any quantum system. Our method is based on the analysis of the properties of the linear operator space induced by unitary operations. Given any quantum system with a noisy measurement apparatus, our method can output the quantum state and the noise matrix of the detector up to a single gauge degree of freedom. We show that this gauge freedom is unavoidable in the general case, but this degeneracy can be generally broken using prior knowledge on the state or noise properties, thus fixing the gauge for several types of state-noise combinations with no assumptions about noise strength. Such combinations include pure quantum states with arbitrarily correlated errors, and arbitrary states with block independent errors. This framework can further use available prior information about the setting to systematically reduce the number of observations and measurements required for state and noise detection. Our method effectively generalizes existing approaches to the problem, and includes as special cases common settings considered in the literature requiring an uncorrelated or invertible noise matrix, or specific probe states.


title: On Galerkin approximations of the Navier-Stokes equations in the limit of large Grashof numbers

abstract: We examine how stationary solutions to Galerkin approximations of the Navier--Stokes equations behave in the limit as the Grashof number $G$ tends to $\infty$. An appropriate scaling is used to place the Grashof number as a new coefficient of the nonlinear term, while the body force is fixed. A new type of asymptotic expansion, as $G\to\infty$, for a family of solutions is introduced. Relations among the terms in the expansion are obtained by following a procedure that compares and totally orders positive sequences generated by the expansion. The same methodology applies to the case of perturbed body forces and similar results are obtained. We demonstrate with a class of forces and solutions that have convergent asymptotic expansions in $G$. All the results hold in both two and three dimensions, as well as for both no-slip and periodic boundary conditions.


title: Blind Cyclic Prefix-based CFO Estimation in MIMO-OFDM Systems

abstract: Low-complexity estimation and correction of carrier frequency offset (CFO) are essential in orthogonal frequency division multiplexing (OFDM). In this paper, we propose a low-overhead blind CFO estimation technique based on cyclic prefix (CP), in multi-input multi-output (MIMO)-OFDM systems. We propose to use antenna diversity for CFO estimation. Given that the RF chains for all antenna elements at a communication node share the same clock, the carrier frequency offset (CFO) between two points may be estimated by using the combination of the received signal at all antennas. We improve our method by combining the antenna diversity with time diversity by considering the CP for multiple OFDM symbols. We provide a closed-form expression for CFO estimation and present algorithms that can considerably improve the CFO estimation performance at the expense of a linear increase in computational complexity. We validate the effectiveness of our estimation scheme via extensive numerical analysis.


title: Ensuring User-side Fairness in Dynamic Recommender Systems

abstract: User-side group fairness is crucial for modern recommender systems, aiming to alleviate performance disparities among user groups defined by sensitive attributes like gender, race, or age. In the ever-evolving landscape of user-item interactions, continual adaptation to newly collected data is crucial for recommender systems to stay aligned with the latest user preferences. However, we observe that such continual adaptation often exacerbates performance disparities. This necessitates a thorough investigation into user-side fairness in dynamic recommender systems, an area that has been unexplored in the literature. This problem is challenging due to distribution shifts, frequent model updates, and non-differentiability of ranking metrics. To our knowledge, this paper presents the first principled study on ensuring user-side fairness in dynamic recommender systems. We start with theoretical analyses on fine-tuning v.s. retraining, showing that the best practice is incremental fine-tuning with restart. Guided by our theoretical analyses, we propose FAir Dynamic rEcommender (FADE), an end-to-end fine-tuning framework to dynamically ensure user-side fairness over time. To overcome the non-differentiability of recommendation metrics in the fairness loss, we further introduce Differentiable Hit (DH) as an improvement over the recent NeuralNDCG method, not only alleviating its gradient vanishing issue but also achieving higher efficiency. Besides that, we also address the instability issue of the fairness loss by leveraging the competing nature between the recommendation loss and the fairness loss. Through extensive experiments on real-world datasets, we demonstrate that FADE effectively and efficiently reduces performance disparities with little sacrifice in the overall recommendation performance.


title: BubbleDet: A Python package to compute functional determinants for bubble nucleation

abstract: We present a Python package, BubbleDet, for computing one-loop functional determinants around spherically symmetric background fields. This gives the next-to-leading order correction to both the vacuum decay rate, at zero temperature, and to the bubble nucleation rate in first-order phase transitions at finite temperature. For predictions of gravitational wave signals from cosmological phase transitions, this is expected to remove one of the leading sources of theoretical uncertainty. BubbleDet is applicable to arbitrary scalar potentials and in any dimension up to seven. It has methods for fluctuations of scalar fields, including Goldstone bosons, and for gauge fields, but is limited to cases where the determinant factorises into a product of separate determinants, one for each field degree of freedom. To our knowledge, BubbleDet is the first package dedicated to calculating functional determinants in spherically symmetric background


title: Statistical methods for resolving poor uncertainty quantification in machine learning interatomic potentials

abstract: Machine learning interatomic potentials (MLIPs) are promising surrogates for quantum mechanics evaluations in ab-initio molecular dynamics simulations due to their ability to reproduce the energy and force landscape within chemical accuracy at four orders of magnitude less cost. While developing uncertainty quantification (UQ) tools for MLIPs is critical to build production MLIP datasets using active learning, only limited progress has been made and the most robust method, ensembling, still shows low correlation between high error and high uncertainty predictions. Here we develop a rigorous method rooted in statistics for determining an error cutoff that distinguishes regions of high and low UQ performance. The statistical cutoff illuminates that a main cause of the poor UQ performance is due to the machine learning model already describing the entire dataset and not having any datapoints with error greater than the statistical error distribution. Second, we extend the statistical analysis to create an interpretable connection between the error and uncertainty distributions to predict an uncertainty cutoff separating high and low errors. We showcase the statistical cutoff in active learning benchmarks on two datasets of varying chemical complexity for three common UQ methods: ensembling, sparse Gaussian processes, and latent distance metrics and compare them to the true error and random sampling, showing that the statistical cutoff is generalizable to a variety of different UQ methods and protocols and performs similarly to using the true error. Importantly, we conclude that utilizing this uncertainty cutoff enables using significantly lower cost uncertainty quantification tools such as sparse gaussian processes and latent distances compared to ensembling approaches for generating MLIP datasets at a fraction of the cost.


title: Injective edge colorings of degenerate graphs and the oriented chromatic number

abstract: Given a graph $G$, an injective edge-coloring of $G$ is a function $\psi:E(G) \rightarrow \mathbb N$ such that if $\psi(e) = \psi(e')$, then no third edge joins an endpoint of $e$ and an endpoint of $e'$. The injective chromatic index of a graph $G$, written $\chi_{inj}'(G)$, is the minimum number of colors needed for an injective edge coloring of $G$. In this paper, we investigate the injective chromatic index of certain classes of degenerate graphs. First, we show that if $G$ is a $d$-degenerate graph of maximum degree $\Delta$, then $\chi_{inj}'(G) = O(d^3 \log \Delta)$. Next, we show that if $G$ is a graph of Euler genus $g$, then $\chi_{inj}'(G) \leq (3+o(1))g$, which is tight when $G$ is a clique. Finally, we show that the oriented chromatic number of a graph is at most exponential in its injective chromatic index. Using this fact, we prove that the oriented chromatic number of a graph embedded on a surface of Euler genus $g$ has oriented chromatic number at most $O(g^{6400})$, improving the previously known upper bound of $2^{O(g^{\frac{1}{2} + \epsilon})}$ and resolving a conjecture of Aravind and Subramanian.


title: A bicomplex proportional fractional $(\vartheta,\varphi)-$weighted Cauchy-Riemann operator using Riemann-Liouville derivatives with respect to an hyperbolic-valued function

abstract: Based on the Riemann-Liouville derivatives with respect to functions taking values in the set of hyperbolic numbers, we consider a novel bicomplex proportional fractional $(\vartheta,\varphi)-$weighted Cauchy-Riemann operator, involving weights hyperbolic orthogonal bicomplex functions. This operator is defined for the first time here, and its associated fractional Borel-Pompeiu formula is proved as the main result.


title: Deep Reinforcement Learning Based Framework for Mobile Energy Disseminator Dispatching to Charge On-the-Road Electric Vehicles

abstract: The exponential growth of electric vehicles (EVs) presents novel challenges in preserving battery health and in addressing the persistent problem of vehicle range anxiety. To address these concerns, wireless charging, particularly, Mobile Energy Disseminators (MEDs) have emerged as a promising solution. The MED is mounted behind a large vehicle and charges all participating EVs within a radius upstream of it. Unfortuantely, during such V2V charging, the MED and EVs inadvertently form platoons, thereby occupying multiple lanes and impairing overall corridor travel efficiency. In addition, constrained budgets for MED deployment necessitate the development of an effective dispatching strategy to determine optimal timing and locations for introducing the MEDs into traffic. This paper proposes a deep reinforcement learning (DRL) based methodology to develop a vehicle dispatching framework. In the first component of the framework, we develop a realistic reinforcement learning environment termed "ChargingEnv" which incorporates a reliable charging simulation system that accounts for common practical issues in wireless charging deployment, specifically, the charging panel misalignment. The second component, the Proximal-Policy Optimization (PPO) agent, is trained to control MED dispatching through continuous interactions with ChargingEnv. Numerical experiments were carried out to demonstrate the demonstrate the efficacy of the proposed MED deployment decision processor. The experiment results suggest that the proposed model can significantly enhance EV travel range while efficiently deploying a optimal number of MEDs. The proposed model is found to be not only practical in its applicability but also has promises of real-world effectiveness. The proposed model can help travelers to maximize EV range and help road agencies or private-sector vendors to manage the deployment of MEDs efficiently.


title: KilonovAE: Exploring Kilonova Spectral Features with Autoencoders

abstract: Kilonovae are likely a key site of heavy r-process element production in the Universe, and their optical/infrared spectra contain insights into both the properties of the ejecta and the conditions of the r-process. However, the event GW170817/AT2017gfo is the only kilonova so far with well-observed spectra. To understand the diversity of absorption features that might be observed in future kilonovae spectra, we use the TARDIS Monte Carlo radiative transfer code to simulate a suite of optical spectra spanning a wide range of kilonova ejecta properties and r-process abundance patterns. To identify the most common and prominent absorption lines, we perform dimensionality reduction using an autoencoder, and we find spectra clusters in the latent space representation using a Bayesian Gaussian Mixture model. Our synthetic kilonovae spectra commonly display strong absorption by strontium Sr II, yttrium Y II, and zirconium Zr I - II, with strong lanthanide contributions at low electron fractions (Ye < 0.25). When a new kilonova is observed, our machine learning framework will provide context on the dominant absorption lines and key ejecta properties, helping to determine where this event falls within the larger 'zoo' of kilonovae spectra.


title: A general Heegaard Floer surgery formula

abstract: We give several new perspectives on the Heegaard Floer Dehn surgery formulas of Manolescu, Ozsv\'{a}th and Szab\'{o}. Our main result is a new exact triangle in the Fukaya category of the torus which gives a new proof of these formulas. This exact triangle is different from the one which appeared in Ozsv\'{a}th and Szab\'{o}'s original proof. This exact triangle simplifies a number of technical aspects in their proofs and also allows us to prove several new results. A first application is an extensions of the link surgery formula to arbitrary links in closed 3-manifolds, with no restrictions on the link being null-homologous. A second application is a proof that the modules for bordered manifolds with torus boundaries, defined by the author in a previous paper, are invariants. Another application is a simple proof of a version of the surgery formula which computes knot and link Floer complexes in terms of subcubes of the link surgery hypercube. As a final application, we show that the knot surgery algebra is homotopy equivalent to an endomorphism algebra of a sum of two decorated Lagrangians in the torus, mirroring a result of Auroux concerning the algebras of Lipshitz, Ozsv\'{a}th and Thurston.


title: Enabling Cooperative Hybrid Beamforming in TDD-based Distributed MIMO Systems

abstract: Distributed massive MIMO networks are envisioned to realize cooperative multi-point transmission in next-generation wireless systems. For efficient cooperative hybrid beamforming, the cluster of access points (APs) needs to obtain precise estimates of the uplink channel to perform reliable downlink precoding. However, due to the radio frequency (RF) impairments between the transceivers at the two en-points of the wireless channel, full channel reciprocity does not hold which results in performance degradation in the cooperative hybrid beamforming (CHBF) unless a suitable reciprocity calibration mechanism is in place. We propose a two-step approach to calibrate any two hybrid nodes in the distributed MIMO system. We then present and utilize the novel concept of reciprocal tandem to propose a low-complexity approach for jointly calibrating the cluster of APs and estimating the downlink channel. Finally, we validate our calibration technique's effectiveness through numerical simulation.


title: Unveiling Camouflage: A Learnable Fourier-based Augmentation for Camouflaged Object Detection and Instance Segmentation

abstract: Camouflaged object detection (COD) and camouflaged instance segmentation (CIS) aim to recognize and segment objects that are blended into their surroundings, respectively. While several deep neural network models have been proposed to tackle those tasks, augmentation methods for COD and CIS have not been thoroughly explored. Augmentation strategies can help improve the performance of models by increasing the size and diversity of the training data and exposing the model to a wider range of variations in the data. Besides, we aim to automatically learn transformations that help to reveal the underlying structure of camouflaged objects and allow the model to learn to better identify and segment camouflaged objects. To achieve this, we propose a learnable augmentation method in the frequency domain for COD and CIS via Fourier transform approach, dubbed CamoFourier. Our method leverages a conditional generative adversarial network and cross-attention mechanism to generate a reference image and an adaptive hybrid swapping with parameters to mix the low-frequency component of the reference image and the high-frequency component of the input image. This approach aims to make camouflaged objects more visible for detection and segmentation models. Without bells and whistles, our proposed augmentation method boosts the performance of camouflaged object detectors and camouflaged instance segmenters by large margins.


title: The Financial Market of Environmental Indices

abstract: This paper introduces the concept of a global financial market for environmental indices, addressing sustainability concerns and aiming to attract institutional investors. Risk mitigation measures are implemented to manage inherent risks associated with investments in this new financial market. We monetize the environmental indices using quantitative measures and construct country-specific environmental indices, enabling them to be viewed as dollar-denominated assets. Our primary goal is to encourage the active engagement of institutional investors in portfolio analysis and trading within this emerging financial market. To evaluate and manage investment risks, our approach incorporates financial econometric theory and dynamic asset pricing tools. We provide an econometric analysis that reveals the relationships between environmental and economic indicators in this market. Additionally, we derive financial put options as insurance instruments that can be employed to manage investment risks. Our factor analysis identifies key drivers in the global financial market for environmental indices. To further evaluate the market's performance, we employ pricing options, efficient frontier analysis, and regression analysis. These tools help us assess the efficiency and effectiveness of the market. Overall, our research contributes to the understanding and development of the global financial market for environmental indices.


title: Sojourns of fractional Brownian motion queues: transient asymptotics

abstract: We study the asymptotics of sojourn time of the stationary queueing process $Q(t),t\ge0$ fed by a fractional Brownian motion with Hurst parameter $H\in(0,1)$ above a high threshold $u$. For the Brownian motion case $H=1/2$, we derive the exact asymptotics of \[ P\left(\int_{T_1}^{T_2} 1(Q(t)>u+h(u))d t>x \Big{|}Q(0) >u \right) \] as $u\to\infty$, {where $T_1,T_2, x\geq 0$ and $T_2-T_1>x$}, whereas for all $H\in(0,1)$, we obtain sharp asymptotic approximations of \[ P\left( \frac 1 {v(u)} \int_{[T_2(u),T_3(u)]}1(Q(t)>u+h(u))dt>y \Bigl \lvert \frac 1 {v(u)} \int_{[0,T_1(u)]}1(Q(t)>u)dt>x\right), \quad x,y >0 \] as $u\to\infty$, for appropriately chosen $T_i$'s and $v$. Two regimes of the ratio between $u$ and $h(u)$, that lead to qualitatively different approximations, are considered.


title: Adaptive Attack Detection in Text Classification: Leveraging Space Exploration Features for Text Sentiment Classification

abstract: Adversarial example detection plays a vital role in adaptive cyber defense, especially in the face of rapidly evolving attacks. In adaptive cyber defense, the nature and characteristics of attacks continuously change, making it crucial to have robust mechanisms in place to detect and counter these threats effectively. By incorporating adversarial example detection techniques, adaptive cyber defense systems can enhance their ability to identify and mitigate attacks that attempt to exploit vulnerabilities in machine learning models or other systems. Adversarial examples are inputs that are crafted by applying intentional perturbations to natural inputs that result in incorrect classification. In this paper, we propose a novel approach that leverages the power of BERT (Bidirectional Encoder Representations from Transformers) and introduces the concept of Space Exploration Features. We utilize the feature vectors obtained from the BERT model's output to capture a new representation of feature space to improve the density estimation method.


title: Comparison between $\mu^{-}\mu^{+}$ and $e^{-}e^{+}$ colliders for charged Higgs production in 2HDM

abstract: We study the phenomenology of the charged Higgs boson at future muon colliders. We investigate both the pair production $\mu^+ \mu^- \to H^+ H^-$, the single production $\mu^+ \mu^- \to W^\pm H^\mp$, as well as the Vector Boson Fusion (VBF) $\{e^+e^-, \mu^+ \mu^-\} \to \nu \bar{\nu} H^+ H^-$. We show that the neutral Higgs exchange diagrams in the muon collider case can lead to a significant boost in the cross sections through their Yukawa couplings. Our results for the muon collider are systematically compared to the corresponding ones at $e^+e^-$ machines. It is demonstrated that the vector boson fusion (VBF) $e^+e^- \to \nu \bar{\nu} H^+ H^-$ can compete with the mentioned $2\to 2$ processes. We select benchmark points and perform signal-background analyses, taking into account detector simulations. We demonstrate the discovery region at $5\sigma$ and the excluded region at $2\sigma$ levels at a 3 TeV muon collider.


title: Edge Magnetoplasmon Dispersion and Time-Resolved Plasmon Transport in a Quantum Anomalous Hall Insulator

abstract: A quantum anomalous Hall (QAH) insulator breaks reciprocity by combining magnetic polarization and spin-orbit coupling to generate a unidirectional transmission of signals in the absence of an external magnetic field. Such behavior makes QAH materials a good platform for the innovation of circulator technologies. However, it remains elusive as to how the wavelength of the chiral edge plasmon relates to its frequency and how the plasmon wave packet is excited in the time domain in a QAH insulator. Here, we investigate the edge magnetoplasmon (EMP) resonances in Cr-(Bi,Sb)$_2$Te$_3$ by frequency and time domain measurements. From disk shaped samples with various dimensions, we obtain the dispersion relation of EMPs and extract the drift velocity of the chiral edge state. From the time-resolved transport measurements, we identify the velocity of the plasmon wave packet and observe a transition from the edge to bulk transport at an elevated temperature. We show that the frequency and time domain measurements are well modeled by loss from the microwave induced dissipative channels in the bulk area. Our results demonstrate that the EMP decay rate can be significantly reduced by applying a low microwave power and fabricating devices of larger diameter $\ge100~\mu$m. In a $R=125~\mu$m sample, a non-reciprocity of 20 dB has been realized at 1.3 GHz, shining light on using QAH insulators to develop on-chip non-reciprocal devices.


title: Convergence of non-linear diagonal frame filtering for regularizing inverse problems

abstract: Inverse problems are key issues in several scientific areas, including signal processing and medical imaging. Since inverse problems typically suffer from instability with respect to data perturbations, a variety of regularization techniques have been proposed. In particular, the use of filtered diagonal frame decompositions has proven to be effective and computationally efficient. However, existing convergence analysis applies only to linear filters and a few non-linear filters such as soft thresholding. In this paper, we analyze filtered diagonal frame decompositions with general non-linear filters. In particular, our results generalize SVD-based spectral filtering from linear to non-linear filters as a special case. As a first approach, we establish a connection between non-linear diagonal frame filtering and variational regularization, allowing us to use results from variational regularization to derive the convergence of non-linear spectral filtering. In the second approach, as our main theoretical results, we relax the assumptions involved in the variational case while still deriving convergence. Furthermore, we discuss connections between non-linear filtering and plug-and-play regularization and explore potential benefits of this relationship.


title: Bridging Distribution Learning and Image Clustering in High-dimensional Space

abstract: Distribution learning focuses on learning the probability density function from a set of data samples. In contrast, clustering aims to group similar objects together in an unsupervised manner. Usually, these two tasks are considered unrelated. However, the relationship between the two may be indirectly correlated, with Gaussian Mixture Models (GMM) acting as a bridge. In this paper, we focus on exploring the correlation between distribution learning and clustering, with the motivation to fill the gap between these two fields, utilizing an autoencoder (AE) to encode images into a high-dimensional latent space. Then, Monte-Carlo Marginalization (MCMarg) and Kullback-Leibler (KL) divergence loss are used to fit the Gaussian components of the GMM and learn the data distribution. Finally, image clustering is achieved through each Gaussian component of GMM. Yet, the "curse of dimensionality" poses severe challenges for most clustering algorithms. Compared with the classic Expectation-Maximization (EM) Algorithm, experimental results show that MCMarg and KL divergence can greatly alleviate the difficulty. Based on the experimental results, we believe distribution learning can exploit the potential of GMM in image clustering within high-dimensional space.


title: Intersectional Inquiry, on the Ground and in the Algorithm

abstract: This article makes two key contributions to methodological debates in automation research. First, we argue for and demonstrate how methods in this field must account for intersections of social difference, such as race, class, ethnicity, culture, and disability, in more nuanced ways. Second, we consider the complexities of bringing together computational and qualitative methods in an intersectional methodological approach while also arguing that in their respective subjects (machines and human subjects) and conceptual scope they enable a specific dialogue on intersectionality and automation to be articulated. We draw on field reflections from a project that combines an analysis of intersectional bias in language models with findings from a community workshop on the frustrations and aspirations produced through engagement with everyday AI-driven technologies in the context of care.


title: ACER: An AST-based Call Graph Generator Framework

abstract: We introduce ACER, an AST-based call graph generator framework. ACER leverages tree-sitter to interface with any language. We opted to focus on generators that operate on abstract syntax trees (ASTs) due to their speed and simplicitly in certain scenarios; however, a fully quantified intermediate representation usually provides far better information at the cost of requiring compilation. To evaluate our framework, we created two context-insensitive Java generators and compared them to existing open-source Java generators.


title: Multimodal Foundation Models For Echocardiogram Interpretation

abstract: Multimodal deep learning foundation models can learn the relationship between images and text. In the context of medical imaging, mapping images to language concepts reflects the clinical task of diagnostic image interpretation, however current general-purpose foundation models do not perform well in this context because their training corpus have limited medical text and images. To address this challenge and account for the range of cardiac physiology, we leverage 1,032,975 cardiac ultrasound videos and corresponding expert interpretations to develop EchoCLIP, a multimodal foundation model for echocardiography. EchoCLIP displays strong zero-shot (not explicitly trained) performance in cardiac function assessment (external validation left ventricular ejection fraction mean absolute error (MAE) of 7.1%) and identification of implanted intracardiac devices (areas under the curve (AUC) between 0.84 and 0.98 for pacemakers and artificial heart valves). We also developed a long-context variant (EchoCLIP-R) with a custom echocardiography report text tokenizer which can accurately identify unique patients across multiple videos (AUC of 0.86), identify clinical changes such as orthotopic heart transplants (AUC of 0.79) or cardiac surgery (AUC 0.77), and enable robust image-to-text search (mean cross-modal retrieval rank in the top 1% of candidate text reports). These emergent capabilities can be used for preliminary assessment and summarization of echocardiographic findings.


title: Lower Bound for Independence Covering in $C_4$-Free Graphs

abstract: An independent set in a graph $G$ is a set $S$ of pairwise non-adjacent vertices in $G$. A family $\mathcal{F}$ of independent sets in $G$ is called a $k$-independence covering family if for every independent set $I$ in $G$ of size at most $k$, there exists an $S \in \mathcal{F}$ such that $I \subseteq S$. Lokshtanov et al. [ACM Transactions on Algorithms, 2018] showed that graphs of degeneracy $d$ admit $k$-independence covering families of size $\binom{k(d+1)}{k} \cdot 2^{o(kd)} \cdot \log n$, and used this result to design efficient parameterized algorithms for a number of problems, including STABLE ODD CYCLE TRANSVERSAL and STABLE MULTICUT. In light of the results of Lokshtanov et al. it is quite natural to ask whether even more general families of graphs admit $k$-independence covering families of size $f(k)n^{O(1)}$. Graphs that exclude a complete bipartite graph $K_{d+1,d+1}$ with $d+1$ vertices on both sides as a subgraph, called $K_{d+1,d+1}$-free graphs, are a frequently considered generalization of $d$-degenerate graphs. This motivates the question whether $K_{d,d}$-free graphs admit $k$-independence covering families of size $f(k,d)n^{O(1)}$. Our main result is a resounding "no" to this question -- specifically we prove that even $K_{2,2}$-free graphs (or equivalently $C_4$-free graphs) do not admit $k$-independence covering families of size $f(k)n^{\frac{k}{4}-\epsilon}$.


title: Asymptotics for Short Maturity Asian Options in Jump-Diffusion models with Local Volatility

abstract: We present a study of the short maturity asymptotics for Asian options in a jump-diffusion model with a local volatility component, where the jumps are modeled as a compound Poisson process. The analysis for out-of-the-money Asian options is extended to models with L\'evy jumps, including the exponential L\'{e}vy model as a special case. Both fixed and floating strike Asian options are considered. Explicit results are obtained for the first-order asymptotics of the Asian options prices for a few popular models in the literature: the Merton jump-diffusion model, the double-exponential jump model, and the Variance Gamma model. We propose an analytical approximation for Asian option prices which satisfies the constraints from the short-maturity asymptotics, and test it against Monte Carlo simulations. The asymptotic results are in good agreement with numerical simulations for sufficiently small maturity.


title: MDTD: A Multi Domain Trojan Detector for Deep Neural Networks

abstract: Machine learning models that use deep neural networks (DNNs) are vulnerable to backdoor attacks. An adversary carrying out a backdoor attack embeds a predefined perturbation called a trigger into a small subset of input samples and trains the DNN such that the presence of the trigger in the input results in an adversary-desired output class. Such adversarial retraining however needs to ensure that outputs for inputs without the trigger remain unaffected and provide high classification accuracy on clean samples. In this paper, we propose MDTD, a Multi-Domain Trojan Detector for DNNs, which detects inputs containing a Trojan trigger at testing time. MDTD does not require knowledge of trigger-embedding strategy of the attacker and can be applied to a pre-trained DNN model with image, audio, or graph-based inputs. MDTD leverages an insight that input samples containing a Trojan trigger are located relatively farther away from a decision boundary than clean samples. MDTD estimates the distance to a decision boundary using adversarial learning methods and uses this distance to infer whether a test-time input sample is Trojaned or not. We evaluate MDTD against state-of-the-art Trojan detection methods across five widely used image-based datasets: CIFAR100, CIFAR10, GTSRB, SVHN, and Flowers102; four graph-based datasets: AIDS, WinMal, Toxicant, and COLLAB; and the SpeechCommand audio dataset. MDTD effectively identifies samples that contain different types of Trojan triggers. We evaluate MDTD against adaptive attacks where an adversary trains a robust DNN to increase (decrease) distance of benign (Trojan) inputs from a decision boundary.


title: Predict And Prevent DDOS Attacks Using Machine Learning and Statistical Algorithms

abstract: A malicious attempt to exhaust a victim's resources to cause it to crash or halt its services is known as a distributed denial-of-service (DDoS) attack. DDOS attacks stop authorized users from accessing specific services available on the Internet. It targets varying components of a network layer and it is better to stop into layer 4 (transport layer) of the network before approaching a higher layer. This study uses several machine learning and statistical models to detect DDoS attacks from traces of traffic flow and suggests a method to prevent DDOS attacks. For this purpose, we used logistic regression, CNN, XGBoost, naive Bayes, AdaBoostClassifier, KNN, and random forest ML algorithms. In addition, data preprocessing was performed using three methods to identify the most relevant features. This paper explores the issue of improving the DDOS attack detection accuracy using the latest dataset named CICDDoS2019, which has over 50 million records. Because we employed an extensive dataset for this investigation, our findings are trustworthy and practical. Our target class (attack class) was imbalanced. Therefore, we used two techniques to deal with imbalanced data in machine learning. The XGboost machine learning model provided the best detection accuracy of (99.9999%) after applying the SMOTE approach to the target class, outperforming recently developed DDoS detection systems. To the best of our knowledge, no other research has worked on the most recent dataset with over 50 million records, addresses the statistical technique to select the most significant feature, has this high accuracy, and suggests ways to avoid DDOS attackI.


title: Single and coupled cavity mode sensing schemes using a diagnostic field

abstract: Precise optical mode matching is of critical importance in experiments using squeezed-vacuum states. Automatic spatial-mode matching schemes have the potential to reduce losses and improve loss stability. However, in quantum-enhanced coupled-cavity experiments, such as gravitational-wave detectors, one must also ensure that the sub-cavities are also mode matched. We propose a new mode sensing scheme, which works for simple and coupled cavities. The scheme requires no moving parts, nor tuning of Gouy phases. Instead a diagnostic field tuned to the HG20/LG10 mode frequency is used. The error signals are derived to be proportional to the difference in waist position, and difference in Rayleigh ranges, between the sub-cavity eigenmodes. The two error signals are separable by 90 degrees of demodulation phase. We demonstrate reasonable error signals for a simplified Einstein Telescope optical design. This work will facilitate routine use of extremely high levels of squeezing in current and future gravitational-wave detectors.


title: Single-ancilla ground state preparation via Lindbladians

abstract: We design a quantum algorithm for ground state preparation in the early fault tolerant regime. As a Monte Carlo-style quantum algorithm, our method features a Lindbladian where the target state is stationary, and its evolution can be efficiently implemented using just one ancilla qubit. Our algorithm can prepare the ground state even when the initial state has zero overlap with the ground state, bypassing the most significant limitation of methods like quantum phase estimation. As a variant, we also propose a discrete-time algorithm, demonstrating even better efficiency and providing a near-optimal simulation cost depending on the desired evolution time and precision. Numerical simulation using Ising models and Hubbard models demonstrates the efficacy and applicability of our method.


title: Transverse Wind Velocity Recorded in Spiral-Shell Pattern

abstract: The propagation speed of a circumstellar pattern revealed in the plane of the sky is often assumed to represent the expansion speed of the wind matter ejected from a post-main-sequence star at the center. We point out that the often-adopted isotropic wind assumption and the binary hypothesis as the underlying origin for the circumstellar pattern in the shape of multilayered shells are, however, mutually incompatible. We revisit the hydrodynamic models for spiral-shell patterns induced by the orbital motion of a hypothesized binary, of which one star is losing mass at a high rate. The distributions of transverse wind velocities as a function of position angle in the plane of the sky are explored along viewing directions. The variation of the transverse wind velocity is as large as half the average wind velocity over the entire three dimensional domain in the simulated models investigated in this work. The directional dependence of the wind velocity is indicative of the overall morphology of the circumstellar material, implying that kinematic information is an important ingredient in modeling the snapshot monitoring (often in the optical and near-infrared) or the spectral imaging observations for molecular line emissions.


title: Understanding step selection analysis through numerical integration

abstract: Step selection functions (SSFs) are flexible models to jointly describe animals' movement and habitat preferences. Their popularity has grown rapidly and extensions have been developed to increase their utility, including various distributions to describe movement constraints, interactions to allow movements to depend on local environmental features, and random effects and latent states to account for within- and among-individual variability. Although the SSF is a relatively simple statistical model, its presentation has not been consistent in the literature, leading to confusion about model flexibility and interpretation. We believe that part of the confusion has arisen from the conflation of the SSF model with the methods used for parameter estimation. Notably, conditional logistic regression can be used to fit SSFs in exponential form, and this approach is often presented interchangeably with the actual model (the SSF itself). However, reliance on conditional logistic regression reduces model flexibility, and suggests a misleading interpretation of step selection analysis as being equivalent to a case-control study. In this review, we explicitly distinguish between model formulation and inference technique, presenting a coherent framework to fit SSFs based on numerical integration and maximum likelihood estimation. We provide an overview of common numerical integration techniques, and explain how they relate to step selection analyses. This framework unifies different model fitting techniques for SSFs, and opens the way for improved inference. In particular, it makes it straightforward to model movement with distributions outside the exponential family, and to apply different SSF formulations to a data set and compare them with AIC. By separating the model formulation from the inference technique, we hope to clarify many important concepts in step selection analysis.


title: The porous envelope and circumstellar wind matter of the closest carbon star, CW Leonis

abstract: Recent abrupt changes of CW Leonis may indicate that we are witnessing the moment that the central carbon star is evolving off the Asymptotic Giant Branch (AGB) and entering into the pre-planetary nebula (PPN) phase. The recent appearance of a red compact peak at the predicted stellar position is possibly an unveiling event of the star, and the radial beams emerging from the stellar position resemble the feature of the PPN Egg Nebula. The increase of light curve over two decades is also extraordinary, and it is possibly related to the phase transition. Decadal-period variations are further found in the residuals of light curves, in the relative brightness of radial beams, and in the extended halo brightness distribution. Further monitoring of the recent dramatic and decadal-scale changes of this most well-known carbon star CW Leonis at the tip of AGB is still highly essential, and will help us gain a more concrete understanding on the conditions for transition between the late stellar evolutionary phases.


title: Optimal decay for one-dimensional damped wave equations with potentials via a variant of Nash inequality

abstract: The optimality of decay properties of the one-dimensional damped wave equations with potentials belonging to a certain class is discussed. The typical ingredient is a variant of Nash inequality which involves an invariant measure for the corresponding Schr\"odinger semigroup. This enables us to find a sharp decay estimate from above. Moreover, the use of a test function method with the Nash-type inequality provides the decay estimate from below. The diffusion phenomena for the damped wave equations with potentials are also considered.


title: Scalable Composite Likelihood Estimation of Probit Models with Crossed Random Effects

abstract: Crossed random effects structures arise in many scientific contexts. They raise severe computational problems with likelihood computations scaling like $N^{3/2}$ or worse for $N$ data points. In this paper we develop a new composite likelihood approach for crossed random effects probit models. For data arranged in R rows and C columns, the likelihood function includes a very difficult R + C dimensional integral. The composite likelihood we develop uses the marginal distribution of the response along with two hierarchical models. The cost is reduced to $\mathcal{O}(N)$ and it can be computed with $R + C$ one dimensional integrals. We find that the commonly used Laplace approximation has a cost that grows superlinearly. We get consistent estimates of the probit slope and variance components from our composite likelihood algorithm. We also show how to estimate the covariance of the estimated regression coefficients. The algorithm scales readily to a data set of five million observations from Stitch Fix with $R + C > 700{,}000$.


title: Who knows what dark matter lurks in the heart of M87: The shadow knows, and so does the ringdown

abstract: We calculate the effect of dark matter on the ringdown waveform and shadow of supermassive black holes at the core of galaxies. Our main focus is on the supermassive black hole at the core of M87, which is large enough to allow for viable observational data. We compare the effects of a dark matter spike to those expected from a galactic halo of the same mass. The radial pressure is shown to be negligible for both the spike and the halo, implying that there is no difference between the isotropic case and the anisotropic case. Our calculation for the halo starts from the Hernquist density function for which the corresponding metric can be obtained analytically in closed form. The effect of the spike is orders of magnitude more significant than the halo as long as the distribution scale of the latter is within a few orders of magnitude of the value expected from observations. Our results indicate that the impact of the spike surrounding M87* on the ringdown waveform may in principle be detectable. Finally, we point out the somewhat surprising fact that existing Event Horizon Telescope observations of black hole shadows are within an order of magnitude from being able to detect, or rule out, the presence of a spike.


title: A spectrum adaptive kernel polynomial method

abstract: The kernel polynomial method (KPM) is a powerful numerical method for approximating spectral densities. Typical implementations of the KPM require an a prior estimate for an interval containing the support of the target spectral density, and while such estimates can be obtained by classical techniques, this incurs addition computational costs. We propose an spectrum adaptive KPM based on the Lanczos algorithm without reorthogonalization which allows the selection of KPM parameters to be deferred to after the expensive computation is finished. Theoretical results from numerical analysis are given to justify the suitability of the Lanczos algorithm for our approach, even in finite precision arithmetic. While conceptually simple, the paradigm of decoupling computation from approximation has a number of practical and pedagogical benefits which we highlight with numerical examples.